Search for: All records

Abstract Historical observations of Earth’s climate underpin our knowledge and predictions of climate variability and change. However, the observations are incomplete and uncertain, and existing datasets based on these observations typically do not assimilate observations simultaneously across different components of the climate system, yielding inconsistencies that limit understanding of coupled climate dynamics. Here, we use coupled data assimilation, which synthesizes observational and dynamical constraints across all climate fields simultaneously, to reconstruct globally resolved sea surface temperature (SST), near-surface air temperature (T), sea level pressure (SLP), and sea ice concentration (SIC), over 1850–2023. We use a Kalman filter and forecasts from an efficient emulator, the linear inverse model (LIM), to assimilate observations of SST, landT, marine SLP, and satellite-era SIC. We account for model error by training LIMs on eight CMIP6 models, and we use the LIMs to generate eight independent reanalyses with 200 ensemble members, yielding 1600 total members. Key findings in the tropics include post-1980 trends in the Walker circulation that are consistent with past variability, whereas the tropical SST contrast (the difference between warmer and colder SSTs) shows a distinct strengthening since 1975. El Niño–Southern Oscillation (ENSO) amplitude exhibits substantial low-frequency variability and a local maximum in variance over 1875–1910. In polar regions, we find a muted cooling trend in the Southern Ocean post-1980 and substantial uncertainty. Changes in Antarctic sea ice are relatively small between 1850 and 2000, while Arctic sea ice declines by 0.5 ± 0.1 (1σ) million km²during the 1920s. Significance StatementThe key advance in our reconstruction is that the ocean, atmosphere, and sea ice are dynamically consistent with each other and with observations across all components, thus forming a true climate reanalysis. Existing climate datasets are typically derived separately for each component (e.g., atmosphere, ocean, and sea ice), leading to spurious trends and inconsistencies in coupled climate variability. We use coupled data assimilation to unify observations and coupled dynamics across components. We combine forecasts from climate models with observations from ocean vessels and weather stations to produce monthly state estimates spanning 1850–2023 and a novel quantification of globally resolved uncertainty. This reconstruction provides insights into historical variability and trends while motivating future efforts to reduce uncertainties in the climate record. 

Abstract In the equatorial and subtropical east Pacific Ocean, strong ocean‐atmosphere coupling results in large‐amplitude interannual variability. Recent literature debates whether climate models reproduce observed short and long‐term surface temperature trends in this region. We reconcile the debate by reevaluating a large range of trends in initial condition ensembles of 15 climate models. We confirm that models fail to reproduce long‐term trends, but also find that many models do not reproduce the observed decadal‐scale swings in the East to West gradient of the equatorial Pacific. Models with high climate sensitivity are less likely to reproduce observed decadal‐scale swings than models with a modest climate sensitivity, possibly due to an incorrect balance of cloud feedbacks driven by changing inversion strength versus surface warming. Our findings suggest that two not well understood problems of the current generation of climate models are connected and we highlight the need to increase understanding of decadal‐scale variability. 

Abstract Coupled global climate models (GCMs) generally fail to reproduce the observed sea‐surface temperature (SST) trend pattern since the 1980s. The model‐observation discrepancies may arise in part from the lack of realistic Antarctic ice‐sheet meltwater input in GCMs. Here we employ two sets of CESM1‐CAM5 simulations forced by anomalous Antarctic meltwater fluxes over 1980–2013 and through the 21st century. Both show a reduced global warming rate and an SST trend pattern that better resembles observations. The meltwater drives surface cooling in the Southern Ocean and the tropical southeast Pacific, in turn increasing low‐cloud cover and driving radiative feedbacks to become more stabilizing (corresponding to a lower effective climate sensitivity). These feedback changes can contribute as substantially as ocean heat uptake efficiency changes in reducing the global warming rate. Accurately projecting historical and future warming thus requires improved representation of Antarctic meltwater and its impacts. 

Abstract Observed surface temperature trends over recent decades are characterized by (a) intensified warming in the Indo‐Pacific Warm Pool and slight cooling in the eastern equatorial Pacific, consistent with Walker circulation strengthening, and (b) Southern Ocean cooling. In contrast, state‐of‐the‐art coupled climate models generally project enhanced warming in the eastern equatorial Pacific, Walker circulation weakening, and Southern Ocean warming. Here we investigate the ability of 16 climate model large ensembles to reproduce observed sea‐surface temperature and sea‐level pressure trends over 1979–2020 through a combination of externally forced climate change and internal variability. We find large‐scale differences between observed and modeled trends that are very unlikely (<5% probability) to occur due to internal variability as represented in models. Disparate trends in the ratio of Indo‐Pacific Warm Pool to tropical‐mean warming, which shows little multi‐decadal variability in models, hint that model biases in the response to historical forcing constitute part of the discrepancy. 

The observed rate of global warming since the 1970s has been proposed as a strong constraint on equilibrium climate sensitivity (ECS) and transient climate response (TCR)—key metrics of the global climate response to greenhouse-gas forcing. Using CMIP5/6 models, we show that the inter-model relationship between warming and these climate sensitivity metrics (the basis for the constraint) arises from a similarity in transient and equilibrium warming patterns within the models, producing an effective climate sensitivity (EffCS) governing recent warming that is comparable to the value of ECS governing long-term warming under CO ${}_2$forcing. However, CMIP5/6 historical simulations do not reproduce observed warming patterns. When driven by observed patterns, even high ECS models produce low EffCS values consistent with the observed global warming rate. The inability of CMIP5/6 models to reproduce observed warming patterns thus results in a bias in the modeled relationship between recent global warming and climate sensitivity. Correcting for this bias means that observed warming is consistent with wide ranges of ECS and TCR extending to higher values than previously recognized. These findings are corroborated by energy balance model simulations and coupled model (CESM1-CAM5) simulations that better replicate observed patterns via tropospheric wind nudging or Antarctic meltwater fluxes. Because CMIP5/6 models fail to simulate observed warming patterns, proposed warming-based constraints on ECS, TCR, and projected global warming are biased low. The results reinforce recent findings that the unique pattern of observed warming has slowed global-mean warming over recent decades and that how the pattern will evolve in the future represents a major source of uncertainty in climate projections. 

Abstract. Antarctic sea ice has exhibited significant variability over the satellite record, including a period of prolonged and gradual expansion, as well as a period of sudden decline. A number of mechanisms have been proposed to explain this variability, but how each mechanism manifests spatially and temporally remains poorly understood. Here, we use a statistical method called low-frequency component analysis to analyze the spatiotemporal structure of observed Antarctic sea ice concentration variability. The identified patterns reveal distinct modes of low-frequency sea ice variability. The leading mode, which accounts for the large-scale, gradual expansion of sea ice, is associated with the Interdecadal Pacific Oscillation and resembles the observed sea surface temperature trend pattern that climate models have trouble reproducing. The second mode is associated with the central Pacific El Niño–Southern Oscillation (ENSO) and the Southern Annular Mode and accounts for most of the sea ice variability in the Ross Sea. The third mode is associated with the eastern Pacific ENSO and Amundsen Sea Low and accounts for most of the pan-Antarctic sea ice variability and almost all of the sea ice variability in the Weddell Sea. The third mode is also related to periods of abrupt Antarctic sea ice decline that are associated with a weakening of the circumpolar westerlies, which favors surface warming through a shoaling of the ocean mixed layer and decreased northward Ekman heat transport. Broadly, these results suggest that climate model biases in long-term Antarctic sea ice and large-scale sea surface temperature trends are related to each other and that eastern Pacific ENSO variability is a key ingredient for abrupt Antarctic sea ice changes. 

Abstract. Antarctic sea ice gradually increased from the late 1970s until 2016, when it experienced an abrupt decline. A number of mechanisms have been proposed for both the gradual increase and abrupt decline of Antarctic sea ice, but how each mechanism manifests spatially and temporally remains poorly understood. Here, we use a statistical method called low-frequency component analysis to analyze the spatial-temporal structure of observed Antarctic sea-ice concentration variability. The identified patterns reveal distinct modes of low-frequency sea ice variability. The leading mode, which accounts for the large-scale, gradual expansion of sea ice, is associated with the Interdecadal Pacific Oscillation and resembles the observed sea-surface temperature trend pattern that climate models have trouble reproducing. The second mode is associated with the central Pacific El Niño–Southern Oscillation (ENSO) and the Southern Annular Mode, and accounts for most of the sea ice variability in the Ross Sea. The third mode is associated with the eastern Pacific ENSO and Amundsen Sea Low, and accounts for most of the pan-Antarctic sea-ice variability and almost all of the sea ice variability in the Weddell Sea. This mode is associated with periods of abrupt Antarctic sea-ice decline and is related to a weakening of the circumpolar westerlies, which favors surface warming through a shoaling of the ocean mixed layer and decreased northward Ekman heat convergence. Broadly, these results suggest that climate model biases in long-term Antarctic sea-ice and global sea-surface temperature trends are related to each other and that eastern Pacific ENSO variability causes abrupt sea ice changes. 

Abstract. The Arctic sea ice cover is strongly influenced by internal variability on decadal time scales, affecting both short-term trends and the timing of the first ice-free summer. Several mechanisms of variability have been proposed, but how these mechanisms manifest both spatially and temporally remains unclear. The relative contribution of internal variability to observed Arctic sea ice changes also remains poorly quantified. Here, we use a novel technique called low-frequency component analysis to identify the dominant patterns of winter and summer decadal Arctic sea-ice variability in the satellite record. The identified patterns account for most of the observed regional sea ice variability and trends, and thus help to disentangle the role of forced and internal sea ice changes over the satellite record. In particular, we identify a mode of decadal ocean-atmosphere-sea ice variability, characterized by an anomalous atmospheric circulation over the central Arctic, that accounts for approximately 30 % of the accelerated decline in pan-Arctic summer sea-ice area between 2000 and 2012. For winter sea ice, we find that internal variability has dominated decadal trends in the Bering Sea, but has contributed less to trends in the Barents and Kara Seas. These results, which detail the first purely observation-based estimate of the contribution of internal variability to Arctic sea ice trends, suggest a lower estimate of the contribution from internal variability than most model-based assessments.