Precession modulation of the South Pacific westerly wind belt over the past million years.

The southern westerly wind belt (SWW) interacts with the Antarctic Circumpolar Current and strongly impacts the Southern Ocean carbon budget, and Antarctic ice-sheet dynamics across glacial-interglacial cycles. We investigated precipitation-driven sediment input changes to the Southeast Pacific off the southern margin of the Atacama Desert over the past one million years, revealing strong precession (19/23-ka) cycles. Our simulations with 2 ocean-atmosphere general circulation models suggest that observed cyclic rainfall changes are linked to meridional shifts in water vapor transport from the tropical Pacific toward the southern Atacama Desert. These changes reflect a precessional modulation of the split in the austral winter South Pacific jet stream. For precession maxima, we infer significantly enhanced rainfall in the southern Atacama Desert due to a stronger South Pacific split jet with enhanced subtropical/subpolar jets, and a weaker midlatitude jet. Conversely, we derive dry conditions in northern Chile related to reduced subtropical/subpolar jets and an enhanced midlatitude jet for precession minima. The presence of precessional cycles in the Pacific SWW, and lack thereof in other basins, indicate that orbital-scale changes of the SWW were not zonally homogeneous across the Southern Hemisphere, in contrast to the hemispherewide shifts of the SWW suggested for glacial terminations. The strengthening of the jet is unique to the South Pacific realm and might have affected winter-controlled changes in the mixed layer depth, the formation of intermediate water, and the buildup of sea-ice around Antarctica, with implications for the global overturning circulation and the oceanic storage of atmospheric CO2.

The Southern Hemisphere westerly wind belt (SWW) plays an important role in mid- and high-latitude atmosphere–ocean dynamic. As such, variations in the high-altitude jet stream of the SWW are critical for communicating climate signals between low and high latitudes (1). Together with Southern Ocean buoyancy (2), the SWW represents the major forcing of the Antarctic Circumpolar Current (3). These atmosphere–ocean processes control the upwelling of carbon-rich deep-water masses in the Southern Ocean, thereby affecting atmospheric CO2 changes over glacial cycles (4, 5) and the Holocene (6, 7). In addition, the strength of the midlatitude and subpolar westerlies affects warm Circumpolar Deep-Water incursions onto the West Antarctic continental shelf (8), impacting the stability of the West Antarctic Ice Sheet via submarine basal melting (9). In contrast to the Atlantic and Indian Ocean sectors, during austral winter, the Pacific SWW is characterized by a split of the high-altitude jet stream into strong subtropical and subpolar jets, and a weaker midlatitude jet (10–12). Variations in these jets affect the strength and latitudinal position of the South Pacific surface westerlies, altering the forcing on the wind-driven ocean circulation in the South Pacific.

Beyond the short instrumental record (13), proxy data-based reconstructions of the SWW traditionally focus on the Holocene (6, 7) and the Last Glacial Maximum (LGM) (14, 15). These studies indicate an overall northward shift and/or intensification of the SWW during the LGM and a similar trend from the Early to Late Holocene. At millennial timescales across the last glacial termination and during Marine Isotope Stage 3, it has been argued that the SWW plays a role in the interhemispheric transfer of climate signals between the North Atlantic region and the Southern Hemisphere mid- and high latitudes (4, 16–18). Most previous studies assume zonally symmetric shifts or intensity changes of the westerlies across the Southern Hemisphere at millennial and orbital timescales, but modern variability of the SWW at interannual to decadal timescales also exhibits pronounced zonally asymmetric changes (19, 20), such as modulation of the South Pacific Split Jet (11, 20) or SWW changes related to the El Niño-Southern Oscillation (21).

Evidence for changes in the SWW on orbital timescales and across multiple glacial–interglacial cycles is very limited (22–25). Here, we present an ∼950-ka-long marine sediment record, recovered from the Chilean continental slope at 27.5°S off the southern Atacama Desert, documenting austral winter rain changes at the northernmost SWW in the Southeast Pacific. Our study provides evidence that the austral winter SWW varied zonally asymmetrically on precessional timescales with important implications for seasonal Southern Ocean atmosphere–ocean changes.

Winter Rain Changes in Northern Chile and the Strength of the Pacific Subropical Jet (STJ)

Sedimentological data are based on sediment core GeoB3375-1 (26), combined with a composite record from the Center for Marine Environmental Sciences (MARUM) deep-sea drill-rig (MeBo) site GeoB15016 (27) (). The sites are located ∼25–35 km WSW off the mouth of the Copiapó River. Sediments at both sites are mainly composed of siliciclastic material and biogenic carbonate (CaCO3 contents range from ∼5 to ∼40%). Lithologically, the sediments are composed of homogeneous calcareous microfossil-bearing clayey silt with subtle color changes (). We use high-resolution records of logarithm-normalized changes in iron to calcium ratios [log (Fe/Ca) following (28)] to document variations in the amount of terrigenous sediments and grain-size analyses of the siliciclastic sediment fraction to infer changes in the mode of sediment input (Fig. 1 and ). Consistent with earlier observations [e.g., the ∼120-ka record from core GeoB3375-1 (26, 29)], iron-rich sediments [higher log (Fe/Ca)] are finer grained and more poorly sorted, with lower carbonate percent throughout the ∼950-ka record. Conversely, the siliciclastic fraction of more carbonate-rich sediments [lower log (Fe/Ca)] is coarser grained and better sorted.

Fig. 1.

Sedimentological data, stratigraphy, and spectral pattern. (A) Benthic foraminifera oxygen isotope data (blue: core GeoB3375-1, red: GeoB15016) versus LR04 stack (59). (B) log (Fe/Ca) record documenting changes in fluvial sediment input (). (C) Precession index. (D) Zoom into last glacial/interglacial cycle with radiocarbon dates and tie points between log (Fe/Ca) and the precession index. Black curve shows intermediate depth Pacific stack (60). (E) Cross-spectral analyses (Blackman–Tukey method; 80% confidence interval) precession versus log (Fe/Ca).

Presently, the Chilean Norte Chico at the southern border region of the Atacama Desert is characterized by prevailing arid conditions and dry valleys with only sporadic water runoff and eolian sediment input to the ocean (30, 31) with mean grain-size and -sorting characteristics similar to those found in our sediment records. Sediments derived from volcanic source rocks in the Andes are mostly trapped in the alluvial basins east of the Coastal Range and do not reach the ocean (). During more humid periods in the past, the dry valleys were activated and the supply of fluvially derived sediments to our coring sites was enhanced (26, 29). Therefore, the well-defined variations in sediment composition off northern Chile indicate changes in the relative contributions of fluvial [higher log (Fe/Ca), finer mean grain size, poorer sorted] and eolian [lower log (Fe/Ca), coarser mean grain size, better sorted] sediment supply. The only modern perennial river in the area is the Copiapó River (), characterized by a primary discharge peak in summer and a secondary peak in winter (32, 33). The main peak is due to snow melt accumulated in the previous winter, while the winter peak is directly recording winter precipitation and/or lagged groundwater discharge originating from the summer snow melt (32, 33).

The log (Fe/Ca) record is characterized by pronounced orbital-scale variations in all Milankovitch orbital bands (Fig. 1). The most striking feature is the strong spectral power in the precessional band (∼19–23-ka cycles), consistent with earlier findings on core GeoB3375-1 (26, 29, 34). Our record extends this cyclicity back to ∼950 ka BP, showing high cross-spectral coherency at both the 19- and 23-ka precessional bands (Fig. 1). Average sedimentation rates of ∼6 cm/ka and reliable age control, primarily based on benthic foraminiferal oxygen isotope (δ18O) stratigraphy (27) (Fig. 1 and ), enable us to resolve these precessional cycles. The phasing between log (Fe/Ca) and precession is constrained by radiocarbon dating for the past ∼45 ka, indicating that higher log (Fe/Ca) and thus, enhanced rainfall occurs during precession maxima (Fig. 1). Based on this radiocarbon-based phasing, we revisited the age model of GeoB15016 and iteratively performed a fine-tuning of the log (Fe/Ca) record to precession back to 950 ka (Fig. 1 and ). The resulting age model is consistent with the phasing of the 41-ky obliquity component in our record ().

We interpret the changes in terrigenous input off northern Chile as variations in river runoff caused by winter rain variability, which is related to the strength and location of the Pacific STJ. Today, precipitation along the Pacific margin of South America is nearly entirely controlled by synoptic disturbances embedded within the SWW and zonal wind strength changes are well correlated to rainfall (35). During austral summer, the STJ is weak and the Pacific westerlies are dominated by a midlatitude eddy-driven jet. Consistent with this westerly configuration, subtropical dry climates extend south to ∼40°S and high rainfall occurs only over the southernmost South America. During austral winter, the STJ strengthens over Australia and southwestern Pacific and extends eastward. As a result, precipitation extends northward to the southern margin of the Atacama Desert, directly onshore from our core locations. In particular, extreme winter rainfall events over subtropical Chile are caused by warm winter storms from atmospheric river events (35). There is no evidence for significant river discharge originating from summer rainfall at the southern boundary of the Atacama Desert (32).

There are no continuous continental rainfall records from subtropical Chile south of the Atacama Desert at orbital timescales beyond the last glacial cycle. Discontinuous regional terrestrial records (36, 37) show several more humid intervals over the past ∼200 ka. The best documented is the LGM precession maximum at ∼21 ka, with glaciological evidence for increased precipitation (38) in line with our data, whereas older intervals of paleosoil development only partially correspond to higher log (Fe/Ca) and thus enhanced continental rainfall in our marine record.

In order to mechanistically explain the remarkably strong precessional cycles in our paleoprecipitation records, we evaluated simulations of 2 coupled ocean–atmosphere general circulation models, EC-Earth 2.2 (39) and the Geophysical Fluid Dynamics Laboratory (GFDL) Climate Model 2.1 (40), both orbitally forced by changes to the phase of precession. For each model, simulations were done with the orbital conditions set to either precession maximum or precession minimum, with other boundary conditions being the same (). Both model simulations show enhanced austral winter precipitation over subtropical Chile during precession maxima (Fig. 2 and ). The enhanced rainfall is tied to an increase in lower tropospheric water vapor transport from the central tropical South Pacific to subtropical Chile (see vectors in Fig. 2). We attribute the altered moisture transport to a stronger STJ over the western sector of the South Pacific, which in turn allows for increased synoptic disturbances further east (see shading in Fig. 2 and ). The disturbances “tap” into the warm, moisture-laden air of the tropical Pacific, allowing for the greater export of moisture to subtropical Chile. The subtropical anticyclone over the southeastern Pacific is also weakened and shifted equatorward during precession maxima (). This is consistent with observations that link the occurrence of austral winter precipitation over subtropical arid Chile to a weaker and equatorward-shifted subtropical anticyclone (41). Overall, the simulated changes to the large-scale circulation over the South Pacific all point to a wetter subtropical Chile during precession maxima as compared to precession minima, in agreement with our paleodata.

Fig. 2.

Simulated June–August climate changes resulting from opposite phases of precession, maximum minus minimum. Shown are the results for EC-Earth; the equivalent figures for GFDL CM2.1 are shown in . June–August averages are shown as those are the months of maximum precipitation over subtropical Chile in the simulation. (A) Precipitation (shaded; units are mm/d) and 850-mb moisture transport (wind speed and direction qv; reference vector is 0.05 m/s). For clarity, vectors with zonal components less than +0.005 m/s are not shown. The circle marks the location of the sediment cores. The figure shows increased moisture transport from the central tropical Pacific to subtropical Chile, leading to increases in rainfall there. (B) Upper tropospheric winds (200 mb; contour interval 2 m/s) and 2–8-d filtered eddy kinetic energy (shaded; units are m2/s2). The latter field is a measure of synoptic eddy activity or “storminess;” we highlight the increased storminess over the southeastern tropical Pacific, coincident with the increased moisture transport and rainfall to subtropical Chile, (C) SST (shaded; units are K). In B, the approximate locations of the climatological subtropical, midlatitude, and subpolar jets are labeled; and in C, the locations of the climatological cold tongue and warm pool regions.

Both models show a cooling of the eastern equatorial Pacific cold tongue during precession maxima (Fig. 2 and ), tied to a phase change in the cold tongue annual sea-surface temperature (SST) cycle (42). The key feature relevant to rainfall over subtropical Chile is the zonal SST gradient across the equatorial Pacific. The climatological wintertime Southern Hemisphere storm track has a pronounced zonally asymmetric configuration: as viewed from the South Pole, it resembles a structure that spirals clockwise and poleward starting north of New Zealand, across the South Pacific and Patagonia, peaking over the Indian Ocean, then decays along the Antarctic edge and terminating at the Ross Sea (42). However, Inatsu and Hoskins (43) showed that when the zonal asymmetry in tropical SST is removed, the Southern Hemisphere storm track is pushed poleward and eastward, and the spiral structure is destroyed. The spiral structure in the storm track reflects the presence of the split jet and is conducive to increased synoptic disturbances over the southeastern tropical Pacific (Fig. 2 and ). This link suggests that increased rainfall over subtropical Chile occurs when the austral winter equatorial zonal SST gradient is large as seen in both model simulations (Fig. 2 and ). Proxy data-based reconstructions of tropical Pacific SST confirm variability at precessional timescales both in the western Pacific warm pool (44) and in the eastern tropical Pacific (44) (). However, interpretations of austral winter SST gradient changes are complicated because the warm pool and cold tongue SST records are based on different proxies (Mg/Ca- versus alkenone-based) and reflect opposite seasons (45).

Orbital-Scale Variability of the SWW

On a global scale, we do not find conclusive proxy evidence for orbital-scale SWW variations in other Southern Hemisphere sectors. Sediment records from the northern margin of the SWW, linked to changes in the SWW south of Africa, suggest northward-extended westerlies during glacials (beyond precessional timescales), which increased winter rain and fluvial supply in SW Africa (25). An 800-ka SST record from the Agulhas Current off southeast Africa (Fig. 3) indicates that these changes were accompanied by northward shifts of the subtropical front (22), which would be consistent with a strengthening of the STJ in the Western Indian Ocean–Atlantic sector. Spectral analysis of the SST record does not reveal significant variability on precessional timescales and is instead dominated by 100-ka glacial-interglacial cycles (). Some evidence exists for processional changes of winter rain in adjacent South Africa based on speleothem and pollen records, but the precession cycles are only identified back to ∼300 ka and are not persistent through time, due to overlapping summer and winter rainfall changes at these locations (46, 47). Taken together, these proxy results are consistent with the reconstructed and modeled strength of the STJ, which is enhanced during precession maxima over the South Pacific, but weaker over the Atlantic and Western Indian Ocean sectors (). This essentially increases the zonal asymmetry in the subtropical westerlies, which are generally stronger in the South Pacific sector and weaker in the other basins.

Fig. 3.

Comparison to other SWW-related records. (A) Precession index. (B) log (Fe/Ca) record as a proxy for changes in the strength of the Pacific STJ. (C) ssNa flux record of the EDC ice core (53, 61). (D) Southeast Pacific SST changes (PS75/34) indicating shifts in the subantarctic front and the midlatitude SWW (23). (E) SST south of SE Africa (MD96/2077) documenting variations in the subtropical front in the western Indian Ocean/Atlantic (22). (F) Temperature record of the EDC ice core (62). Gray bars across A–C mark precession maxima (reduced austral winter and higher austral summer insolation); Gray bars across D–F mark glacial marine isotope stages numbered below; MBT = Mid-Brunhes Event; MPT = Middle Pleistocene Transition.

Within the context of modern climatology, we infer that this unique Pacific pattern is related to a long-term (i.e., precessional) modulation of the South Pacific split jet, not present in the other SWW sectors. A number of South Pacific paleoclimate records over the last glacial termination have been related to a split-jet response to millennial-scale North Atlantic cooling intervals (11). This interhemispheric teleconnection is thought to operate through a southward shift of the Intertropical Convergence Zone and a weakening of the Asian monsoon induced by North Atlantic cooling and is similar to the previously suggested interhemispheric processes across glacial terminations (4, 16, 48). The occurrence of precessional cycles in the Pacific STJ supports the concept of a strong role of orbital precession in modulating interhemispheric temperature gradients and latitudinal shifts of the SWW (48).

Traditionally, precessional forcing has been primarily related to summer insolation changes, particularly as it affects changes in monsoonal rainfall in both hemispheres (49). During precession maxima, higher austral summer insolation enhanced the South American summer monsoon (50), whereas our proxy-reconstructed and modeled enhancement of the STJ during precession maxima occurred in austral winter, the opposite season. At the same precession phase, boreal summer insolation and Northern Hemisphere summer monsoons were reduced (49). This pattern either documents a dynamical link of the Pacific STJ to changes in the monsoonal systems, potentially involving the Northern Hemisphere westerly jet (51), or an independent response of both to precessional forcing.

Our 2 model simulations suggest a close link of a strengthened STJ to enhanced austral winter zonal SST gradients in the equatorial Pacific (Fig. 2 and ). This link indicates a strong tropical forcing on the extratropical westerlies, similar to what is seen in the response to El Niño-Southern Oscillation variation (1). The strong meridional teleconnections in the South Pacific region appear to be reduced during austral summer. An ∼700-ka alkenone SST record from the subantarctic SE Pacific has been interpreted in terms of northward-shifting Southern Ocean fronts and the SWW during glacials and more southward locations during interglacials (23) (PS75/34; Fig. 3). Alkenones are strongly biased toward austral summer at high latitudes (23, 52) and the SST record does not reveal substantial precessional cyclicity. Instead, the SST record is dominated by ∼41- and ∼100-ka cycles, similar to the subtropical front record from the western Indian Ocean/Atlantic sector ().

Connection to High-Latitude Atmosphere–Ocean Processes

Orbital-scale variability in many paleoclimate records from high- and midlatitudes is primarily dominated by ∼41-ka (obliquity) and ∼100-ka (eccentricity) cycles. In contrast, our austral winter STJ records show pronounced precessional cycles and only reduced spectral power at the obliquity and eccentricity orbital band (), the latter most likely reflecting more global climate changes in response to glacial/interglacial cycles. Consistent with modern interannual variations (35) at precessional timescales, our model results indicate out-of-phase variations of the STJ/subpolar jet (stronger during precession maxima) compared to the intercalated midlatitude branch (weaker during precession maxima) of the upper-tropospheric Pacific SWW during austral winter (Fig. 2). Austral winter changes of atmosphere–ocean processes in the Southern Ocean are not well constrained as paleoproxy data are rare. One exception is changes in sea-salt sodium (ssNa) flux recorded in the ∼800-ka EPICA Dome C (EDC) ice-core record (53) (Fig. 3), which have been interpreted in terms of large-scale winter sea-ice extent around Antarctica. Spectral analyses document strong power in the precessional 19/23-ka bands compared to the dominating 100-ka cyclicity in the ssNa record (). These spectral characteristics are different from the local Antarctic temperature record based on the same ice core, which primarily shows ∼41- and ∼100-ka cycles (Fig. 3 and ). The occurrence of strong precessional cycles in both Antarctic ssNa and our Pacific SWW records suggests a mechanistic link. Overall, the ssNa flux is higher during precession maxima (including the LGM) coinciding with our stronger simulated STJ/subpolar jet and weaker midlatitude jet. However, ssNa fluxes at orbital timescales in the EDC ice core may not only reflect large-scale changes in the extent of sea ice around Antarctica (54). Alternatively, the supply of sea spray might be controlled by changes in midlatitude and/or subpolar cyclogenesis (55). Cyclonic activity along the sea-ice edge might act as an effective source for sea spray being uplifted to sufficient altitudes for transport to the high-altitude EDC site (56).

Conclusions

Our combined proxy and climate model results document substantial orbital-scale changes in the South Pacific jet stream configuration over the past 1 Ma, primarily paced by precession. The amplitudes of ∼41- and ∼100-ka cycles, often dominating midlatitude and high-latitude orbital variability, are strongly reduced in our records, suggesting prevailing tropical forcing of the South Pacific westerlies during austral winter. For precession maxima (including the LGM at ∼21 ka), we plausibly infer a stronger South Pacific split jet with enhanced STJ and subpolar jets and weaker midlatitude jet during austral winter (). Conversely, we derive reduced STJ and subpolar jets, and enhanced midlatitude jet, for precession minima. The presence of precessional cycles in the Pacific SWW, and lack thereof in other basins, indicate that orbital-scale changes of the SWW were not zonally homogeneous across the Southern Hemisphere, in contrast to the hemispherewide shifts of the SWW suggested for glacial terminations and on millennial timescales (4, 16, 48).

The unique response of the Pacific sector of the SWW to precessional forcing in the tropics, shown by model and proxy results, strongly supports a tropical Pacific driver at orbital timescales in the South Pacific, at least during austral winter. The atmosphere–ocean pattern during precession maxima including the LGM are similar to modern La Niña-like states. During austral summer, the available proxy data support zonally more symmetric shifts of the SWW paced by ∼41- and ∼100-ka cycles. These patterns resemble modern changes in the Southern Annular Mode (57). Our reconstructed austral winter precessional cycles appear to persist across Late and Middle Pleistocene ∼100-ka-paced glacial–interglacial cycles and major climate reorganizations such as the Mid-Brunhes Event and the Middle Pleistocene Transition. Austral winter SWW changes along the Southern Ocean ACC might have affected winter-controlled changes in the mixed layer depth, the formation of intermediate water, and the built-up of sea ice around Antarctica (20) during past climates. Our results support recent model simulations, suggesting that precessional forcing of wind stress across the Southern Ocean may play an important role in atmosphere–ocean changes that impact atmospheric CO2 concentrations (58). The strong meridional teleconnections in the Pacific SWW might therefore have affected global overturning circulation and the storage of atmospheric CO2 in the deep Pacific Ocean. Furthermore, they might have influenced the exchange of heat, salt, and nutrients between ocean basins, as well as between high and low latitudes.

Methods

We updated the published age models of the sediment records (). High-resolution major element records were obtained with an X-ray fluorescence core scanner (XRF) and quantification using inductively coupled plasma optical emission spectrometry. Grain-size distributions were obtained on a Coulter laser-granulometer (LS-200). We used Blackman–Tukey spectral analyses to derive the relative power of the major cycles in the Milankovitch band. We performed model simulations with EC-Earth 2.2 (39) fully coupled ocean–atmosphere general circulation model as well as with the GFDL Climate Model version 2.1 (GFDL CM2.1) (40).

Data Availability Statement.

All data discussed in the paper have been made available to readers.