Analytical solutions are derived for steady state groundwater flow in a heterogeneous, anisotropic, semiconfined aquifer. The aquifer consists of a number of horizontal layers, while each layer consists of a number of homogeneous cells with different hydraulic conductivity tensors. An exact solution for an arbitrary number of layers and cells is derived using a multilayer approach. Hypothetical models are used to study how flow is affected by heterogeneity in the anisotropy. When the major principal direction of the transmissivity tensor varies between layers, streamlines have the shape of spirals. Bundles of spiraling streamlines rotating in the same direction are called "groundwater whirls." Analytical and numerical models of hypothetical cases produce the same flow patterns of groundwater whirls. Projected streamlines may be represented by stream function contours, which allows for easy interpretation of the main characteristics of complex whirl patterns. Typical whirl patterns may be identified in relatively simple models. Clockwise and counterclockwise whirl axes occur at the interfaces of adjacent layers. One or more clockwise whirls may occur within one counterclockwise whirl, and vice versa. The practical consequence of groundwater whirls is that the exchange of water between aquifer layers is intensified. This may have a significant impact on contaminant spreading.