Search NASA⌕ Search

SEARCH · Search NASA

Results for “Coupled models of the climate system”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

Impact of Tropical Cyclone Wind Forcing on the Global Climate in a Fully Coupled Climate Model

Tropical cyclones (TCs) alter upper-ocean temperature and influence ocean heat content via enhanced turbulent mixing. A better understanding of the role of TCs within the climate system requires a fully coupled modeling framework, where TC-induced ocean responses feed back to the atmosphere and subsequently to the climate mean state and variability. Here, we investigate the impacts of TC wind forcing on the global ocean and the associated feedbacks within the climate system using the fully coupled Community Earth System Model version 1.3 (CESM1.3). Using the low-resolution version of CESM1.3 (1° atmosphere and ocean grid spacing) with no intrinsic TCs, we conduct a suite of sensitivity experiments by inserting TC winds extracted from a high-resolution (0.25° atmosphere grid spacing) TC-permitting simulation into the low-resolution model. Results from the low-resolution TC experiment are compared to a low-resolution control simulation to diagnose TCs’ impact. We found that the added TC winds can increase ocean heat content by affecting ocean vertical mixing, air–sea enthalpy fluxes, and cloud amount. Here, the added TCs can influence mean SST, precipitation, ocean subsurface temperature, and ocean mixed layer depth. We found a strengthening of the wind-driven subtropical cells and a weakening of the Atlantic meridional overturning circulation due to the changes of surface buoyancy fluxes. TCs in the model cause anomalous equatorward ocean heat convergence in the deep tropics and an increase of poleward ocean heat transport out of the subtropics. Our modeling results provide new insights into the multiscale interactions between TCs and the coupled climate system.

54 ENVIRONMENTAL SCIENCES↗

On the Relationship Between the Brewer-Dobson Circulation and the Southern Annular Mode During Austral Summer in Couple Chemistry-Climate Model Simulations

The Brewer-Dobson circulation (BDC) is the mean transport circulation in the stratosphere. It consists of an upwelling branch in the tropics, poleward flows from the tropics to the extratropics, and downward flows in the extratropics. The BDC plays a crucial role in the distribution of important stratospheric trace gases, such as ozone and water vapor. Therefore changes in the strength of the BDC under global warming could have significant impact on stratospheric ozone depletion and recovery. For example, all climate models that are used by the World Meteorological Organization to predict ozone evolution in the 21 st century project a strengthening of the BDC that leads to ozone superrecovery in the mid-latitudes. On the other hand, ozone changes could also affect the strength of the BDC. This work investigates an outstanding question: whether and how changes in the Brewer-Dobson circulation are connected to climate change in the troposphere, in particular, the annular modes. The annular modes are the leading variability in the extratropical troposphere, which describes a seesaw pattern of circulation fluctuations between the polar and middle latitudes. Using simulations from the Goddard Earth Observing System Coupled Chemistry Climate Model (GEOS CCM), we found the strengthening of the BDC in the summer Southern Hemisphere is strongly correlated with a shift of the Southern Hemisphere Annular Mode (SAM) toward its positive phase for the last 4 decades of the 20th century. This relationship is only present in model runs that simulate the stratospheric ozone depletion. Therefore it is concluded that the BDC-SAM relationship is driven by Antarctic ozone depletion. The ozone hole significantly cools the Antarctic stratosphere in late spring/early summer, which leads to a delayed breakdown of the polar vortex: strong circumpolar eastward flows that usually shift to westward winds in late spring. The prolonged persistence of stratospheric eastward flow enhances upward propagation of tropospheric waves into the stratosphere and strengthens the BDC. The increased wave flux in the stratosphere in turn drives a SAM trend toward its positive phase. Our results also show that the BDC-SAM relationship is robust on the interannual timescale

Li, Feng↗

Ensemble Spread Behavior in Coupled Climate Models: Insights From the Energy Exascale Earth System Model Version 1 Large Ensemble

AbstractAssessing uncertainty in future climate projections requires understanding both internal climate variability and external forcing. For this reason, single‐model initial condition large ensembles (SMILEs) run with Earth System Models (ESMs) have recently become popular. Here we present a new 20‐member SMILE with the Energy Exascale Earth System Model version 1 (E3SMv1‐LE), which uses a “macro” initialization strategy choosing coupled atmosphere/ocean states based on inter‐basin contrasts in ocean heat content (OHC). The E3SMv1‐LE simulates tropical climate variability well, albeit with a muted warming trend over the twentieth century due to overly strong aerosol forcing. The E3SMv1‐LE's initial climate spread is comparable to other (larger) SMILEs, suggesting that maximizing inter‐basin ocean heat contrasts may be an efficient method of generating ensemble spread. We also compare different ensemble spread across multiple SMILEs, using surface air temperature and OHC. The Community Earth system Model version 1, the only ensemble which utilizes a “micro” initialization approach perturbing only atmospheric initial conditions, yields lower spread in the first ∼30 years. The E3SMv1‐LE exhibits a relatively large spread, with some evidence for anthropogenic forcing influencing spread in the late twentieth century. However, systematic effects of differing “macro” initialization strategies are difficult to detect, possibly resulting from differing model physics or responses to external forcing. Notably, the method of standardizing results affects ensemble spread: control simulations for most models have either large background trends or multi‐centennial variability in OHC. This spurious disequlibrium behavior is a substantial roadblock to understanding both internal climate variability and its response to forcing.

Stevenson, Samantha↗

Near-term tropical cyclone risk and coupled Earth system model biases

Most current climate models predict that the equatorial Pacific will evolve under greenhouse gas–induced warming to a more El Niño-like state over the next several decades, with a reduced zonal sea surface temperature gradient and weakened atmospheric Walker circulation. Yet, observations over the last 50 y show the opposite trend, toward a more La Niña-like state. Recent research provides evidence that the discrepancy cannot be dismissed as due to internal variability but rather that the models are incorrectly simulating the equatorial Pacific response to greenhouse gas warming. This implies that projections of regional tropical cyclone activity may be incorrect as well, perhaps even in the direction of change, in ways that can be understood by analogy to historical El Niño and La Niña events: North Pacific tropical cyclone projections will be too active, North Atlantic ones not active enough, for example. Other perils, including severe convective storms and droughts, will also be projected erroneously. While it can be argued that these errors are transient, such that the models’ responses to greenhouse gases may be correct in equilibrium, the transient response is relevant for climate adaptation in the next several decades. Given the urgency of understanding regional patterns of climate risk in the near term, it would be desirable to develop projections that represent a broader range of possible future tropical Pacific warming scenarios—including some in which recent historical trends continue—even if such projections cannot currently be produced using existing coupled earth system models.

Science & Technology - Other Topics↗

WCRP Task Team for the Intercomparison of ReAnalyses (TIRA): Motivation and Progress

Reanalyses have proven to be an important resource for weather and climate related research, as well as societal applications at large. New atmospheric reanalyses in various forms are produced every few years. In addition, land and ocean communities are producing disciplinary uncoupled reanalyses, and regional reanalyses are emerging in addition to global reanalyses. Systems of data information are being developed to facilitate the intercomparison of reanalyses among themselves and with observations and climate change experiments (e.g. CREATE-IP and WRIT). Current research and development in reanalysis is directed at (1) extending the length of reanalyzed period and (2) use of coupled Earth system models for climate reanalysis. While the World Climate Research Programme (WCRP) involvement in the reanalyses communities through its Data Advisory Council (WDAC) has been substantial, for example in organizing international conferences on reanalyses. However, these efforts have previously not included an analysis of the differences among reanalyses and their inherent uncertainties which are important questions for both users and developers of reanalyses. Therefore, the WCRP Task Team for the Intercomparison of ReAnalyses (TIRA) task team has been formed as collaborative effort to address reanalyses role in WCRP. This constitutes a logical progression that fills the needs of the community and contributes to the WCRP mission. The primary charge to TIRA is to develop a reanalysis intercomparison project plan that will attain the following objectives.1)To foster understanding and estimation of uncertainties in reanalysis data by intercomparison and other means)To communicate new developments and best practices among the reanalyses producing centers3)To enhance the understanding of data and assimilation issues and their impact on uncertainties, leading to improved reanalyses for climate assessment4)To communicate the strengths and weaknesses of reanalyses, their fitness for purpose, and best practices in the use of reanalysis datasets by the scientific communityThis contribution outlines the need for a task team on reanalyses, their intercomparison, the objectives of the team and progress thus far. Results for some pilot intercomparisons will also be discussed.

Bosilovich, Michael↗

Improving tide-estuary representation in MPAS-Ocean (Final Report)

The overarching goal of this research is to produce a skill-assessed tide-estuary representation, including BGC (biogeochemistry) components of the ecosystem, into the MPAS-O global climate modeling system. Specifically, in the original proposal we proposed one-way nesting, into global MPAS-O, of a high-resolution model for TES (tidal and estuary systems) including the coupling of the BGC components. We have completed the proposed tasks during the past 4 years. Specifically, after extensive consultation with the MPAS-O group during the first year, we have decided to go beyond the original proposal and to two-way nest SCHISM into MPAS-O. As a result, SCHISM-CoSiNE (with CoSiNE being the BGC model) is now fully integrated into MPAS-O as a subroutine (more precisely, an analysis member) and can be invoked as MPAS-O is running with two-way exchange of information between the TES and global models (with proper spatial interpolation done on the fly). As originally proposed, we have also developed the new finite-volume (FV) model MPAS-OI (aka ELCIRC-sub) for TES. We have conducted several tests using the coupled MPAS-SCHISM- CoSiNE model and also MPAS-OI for Columbia River, San Francisco Bay and Chesapeake Bay. Four journal publications have been produced so far as the result of the project and 1 is under review (cf. Publication section). With these achievements, we will be ready to conduct fully coupled simulations of TES and global ocean anywhere, as soon as some missing capabilities in MPAS-O (e.g. tides, improved ability to handle the boundary conditions, and streamlining of the global BGC model) are worked out. The significance of the outcome of this project cannot be overstated. It marks the first time, as far as we know, that a global climate model is tightly coupled to a bona fide nearshore and terrestrial model (SCHISM system); in fact SCHISM is the first bona fide compound flooding model that seamlessly covers the coastal ocean, estuaries, rivers and watersheds (Ye et al. 2020; Zhang et al. 2020). The coupled model will reduce the errors and uncertainties of the MPAS-O global climate model system. It also complements and extends on-going E3SM project, as the latter seeks in the long-term to include components required to simulate impacts of sea-level change and storm surge on coastal regions, including wave models and focusing resolution in coastal and storm-track regions, as an important component of the Cryosphere-ocean system. Note that a spectral wave model (the unstructured-grid version of Wave Watch III, aka Wind Wave Model) is already incorporated inside the SCHISM system and is therefore directly accessible in the coupled system.

58 GEOSCIENCES↗

WCRP Task Team for the Intercomparison of Reanalyses (TIRA): Motivation and Progress

Reanalyses have proven to be an important resource for weather and climate related research, as well as societal applications at large. Several centers have emerged to produce new atmospheric reanalyses in various forms every few years. In addition, land and ocean communities are producing disciplinary uncoupled reanalyses. Current research and development in reanalysis is directed at (1) extending the length of reanalyzed period and (2) use of coupled Earth system models for climate reanalysis. While WCRPs involvement in the reanalyses communities through its Data Advisory Council (WDAC) has been substantial, for example in organizing international conferences on reanalyses, a central team of reanalyses expertise is not in place in the WCRP structure. The differences among reanalyses and their inherent uncertainties are some of the most important questions for both users and developers of reanalyses. Therefore, a collaborative effort to systematically assess and intercompare reanalyses would be a logical progression that fills the needs of the community and contributes to the WCRP mission. The primary charge to the TIRA is to develop a reanalysis intercomparison project plan that will attain the following objectives.1)To foster understanding and estimation of uncertainties in reanalysis data by intercomparison and other means 2)To communicate new developments and best practices among the reanalyses producing centers 3)To enhance the understanding of data and assimilation issues and their impact on uncertainties, leading to improved reanalyses for climate assessment 4)To communicate the strengths and weaknesses of reanalyses, their fitness for purpose, and best practices in the use of reanalysis datasets by the scientific community. This presentation outlines the need for a task team on reanalyses, their intercomparison, the objectives of the team and progress thus far.

WDA↗

Effects of Greenhouse Gas Increase and Stratospheric Ozone Depletion on Stratospheric Mean Age of Air in 1960-2010

The relative impacts of greenhouse gas (GHG) increase and stratospheric ozone depletion on stratospheric mean age of air in the 1960-2010 period are quantified using the Goddard Earth Observing System Chemistry-Climate Model. The experiment compares controlled simulations using a coupled atmosphere-ocean version of the Goddard Earth Observing System Chemistry-Climate Model, in which either GHGs or ozone depleting substances, or both factors evolve over time. The model results show that GHGs and ozone-depleting substances have about equal contributions to the simulated mean age decrease, but GHG increases account for about two thirds of the enhanced strength of the lower stratospheric residual circulation. It is also found that both the acceleration of the diabatic circulation and the decrease of the mean age difference between downwelling and upwelling regions are mainly caused by GHG forcing. The results show that ozone depletion causes an increase in the mean age of air in the Antarctic summer lower stratosphere through two processes: (1) a seasonal delay in the Antarctic polar vortex breakup that inhibits young midlatitude air from mixing with the older air inside the vortex, and (2) enhanced Antarctic downwelling that brings older air from middle and upper stratosphere into the lower stratosphere.

stratospheric ozone depletion↗

Arctic Warming in Response to Regional Aerosol Emissions Reductions

This study examines the Arctic surface air temperature response to regional aerosol emissions reductions using three fully coupled chemistry–climate models: National Center for Atmospheric Research-Community Earth System Model version 1, Geophysical Fluid Dynamics Laboratory-Coupled Climate Model version 3 (GFDL-CM3) and Goddard Institute for Space Studies-ModelE version 2. Each of these models was used to perform a series of aerosol perturbation experiments, in which emissions of different aerosol types (sulfate, black carbon (BC), and organic carbon) in different northern mid-latitude source regions, and of biomass burning aerosol over South America and Africa, were substantially reduced or eliminated. We find that the Arctic warms in nearly every experiment, the only exceptions being the U.S. and Europe BC experiments in GFDL-CM3 in which there is a weak and insignificant cooling. The Arctic warming is generally larger than the global mean warming (i.e. Arctic amplification occurs), particularly during non-summer months. The models agree that changes in the poleward atmospheric moisture transport are the most important factor explaining the spread in Arctic warming across experiments: the largest warming tends to coincide with the largest increases in moisture transport into the Arctic. In contrast, there is an inconsistent relationship (correlation) across experiments between the local radiative forcing over the Arctic and the simulated Arctic warming, with this relationship being positive in one model (GFDL-CM3) and negative in the other two. Our results thus highlight the prominent role of poleward energy transport in driving Arctic warming and amplification, and suggest that the relative importance of poleward energy transport and local forcing/feedbacks is likely to be model dependent.

Arctic↗

Predictability of weather and climate in a coupled ocean-atmosphere model: A dynamical systems approach

A dynamical systems approach is used to quantify the instantaneous and time-averaged predictability of a low-order moist general circulation model. Specifically, the effects on predictability of incorporating an active ocean circulation, implementing annual solar forcing, and asynchronously coupling the ocean and atmosphere are evaluated. The predictability and structure of the model attractors is compared using the Lyapunov exponents, the local divergence rates, and the correlation, fractal, and Lyapunov dimensions. The Lyapunov exponents measure the average rate of growth of small perturbations on an attractor, while the local divergence rates quantify phase-spatial variations of predictability. These local rates are exploited to efficiently identify and distinguish subtle differences in predictability among attractors. In addition, the predictability of monthly averaged and yearly averaged states is investigated by using attractor reconstruction techniques.

Nese, Jon M.↗

Earth System Reanalysis in Support of Climate Model Improvements

Recent climate model developments, established through increased model resolution, have led to substantial improvements in model simulations of the time-evolving, coupled Earth system and its subcomponents. However, regardless of resolution, climate models will always produce climate features and variability that differ from the real world and will be prone to biases. This is due to many remaining uncertainties, such as in parametric and structural model uncertainty, in the initial conditions prescribed, and in the prescribed (scenario) forcing which varies on decadal to centennial timescales. Further model improvements are expected to arise specifically from improved representation of physical processes realized through model-data fusion. This will create an unprecedented opportunity to better exploit a large array of Earth observations, from in situ measurements to weather radars and satellite observations, as the resolved scales of the models approach those of the observations. For this, climate DA will be the central tool to bring models and observations into consistency, by improving initial conditions, inferring uncertain model parameters and structure, and quantifying uncertainty. Generally, there will be advantages and complementarities of adjoint-based smoother approaches, ensemble-based filter approaches, or new ML-inspired approaches. Yet, the ever-increasing model resolution will present growing challenges arising from computational cost, calling for new ways of performing data assimilation and model optimization. Using the complementarity in a hybrid approach, blending tools and concepts from variational, ensemble and ML methods might be what is required in the future. In this context ML could be important to handle non-linear responses, and to better approximate non-Gaussian distributions.

54 ENVIRONMENTAL SCIENCES↗

On a simplified solution of climate-carbon dynamics in idealized flat10MIP simulations

Abstract. Idealized experiments with coupled climate-carbon Earth system models (ESMs) provide a basis for understanding the response of the carbon cycle to external forcing and for quantifying climate-carbon feedbacks. Here, we analyze globally-averaged results from idealized esm-flat10 experiments and show that most models exhibit a quasi-linear relationship between cumulative carbon uptake on land and in the ocean during a period of constant fossil fuel emissions of 10 Pg C yr−1. We hypothesize that this relationship does not depend on emission pathways. Further, as a simplification, we quantify the relationship between cumulative ocean carbon uptake and changes in ocean heat content using a linear approximation. In this way, changes in oceanic heat content and atmospheric CO2 concentration become interdependent variables, reducing the coupled temperature-CO2 system to just one differential equation. The equation can be solved analytically or numerically for the atmospheric CO2 concentration as a function of fossil fuel emissions. This approach leads to a simplified description of global carbon and climate dynamics, which could be used for applications beyond existing analytical frameworks.

Brovkin, Victor↗

Reducing uncertainty of polar to midlatitude linkages using DOE’s E3SM in a coordinated model-experiment setting

This project brought DOE’s climate modeling effort with the Energy Exascale Earth System Model into the Polar Amplification Model Intercomparison Project (PAMIP), which is part of the sixth and latest Coupled Model Intercomparison Project, CMIP6. PAMIP examines the causes and consequences of polar amplification, when external forcing results in a larger temperature increase in high latitudes than the global average, in a coordinated set of model experiments with a common modeling protocol. Our teams from UC Irvine and the University of Toronto have designed, carried out, analyzed, and disseminated PAMIP output from the Energy Exascale Earth System Model (E3SM) and the Community Earth System Model (CESM). PAMIP’s ongoing stream of significant new results have advanced progress in the community’s understanding and led to new outstanding research questions that have motivated further work. PAMIP has led to improved consensus on the atmospheric response to sea ice loss. The important finding is that for a similar sea ice anomaly forcing, the simulated atmospheric response in the troposphere is remarkably consistent among the 16 models’ runs analyzed. The zonal-mean tropospheric response consists of a very robust equatorward shift of the westerly flow in mid-latitudes. However, while the multi-model mean response is robust, it has a weak amplitude relative to internal variability. We identified a weakness in the models (including E3SM) in terms of their response to sea-ice forcing that is related to eddy forcing (or nonlinear dynamical effects) at mid- to high latitudes. In fact, E3SM is an outlier in terms of the models participating, and in that sense, it turned out to be a vital participant model. We found that reductions in energy transport due sea-ice loss and involving dry air only are compensated by increases in moist energy transport from warmer sea surface temperature in midlatitudes. This leads to a large spread in energy transport into the Arctic and is a potential source of spread in Arctic amplification. We identified an important role that climate modes, including tropical modes of variability (El Nino and the Southern Oscillation (ENSO); the Quasi-biennial Oscillation (QBO)) play in the response to sea-ice anomalies, including in ocean coupled experiments. Similarly, we identified and quantified the contribution of sea-ice thickness to the atmospheric response compared to the response to sea-ice extent only. We found that it is important to run large ensembles and even with an ensemble size of 100 simulations the response is largely influenced by internal variability. We demonstrated convincingly that Ural blocking, not sea-ice loss, provides the weakening of the stratospheric polar vortex in fall/early winter and a negative phase of the North Atlantic Oscillation that can last for up to two months. However, sea-ice anomalies can influence the background flow so that the response to Ural blocking is more persistent under low sea-ice conditions in the Barents/Kara Sea than high sea ice Atmospheric model hierarchies that progressively add individual processes have a long history in providing dynamical insight for modeling the atmosphere. Similarly, coupled model hierarchies that progressively add individual ocean processes can provide insights into the workings of the coupled climate system, however such hierarchies have not been available except for a non-dynamical slab ocean model. Because of the missing processes, surface flux corrections must be added to produce a target climate. In this project, we managed to overcome this problem and develop a globally coupled ocean model hierarchy in CESM that can turn on and off the processes of mixed-layer entrainment and Ekman flow. We used the hierarchy to study the impact of Arctic sea-ice loss on the climate system. We find that the effect of mixed-layer entrainment on ocean heat uptake influences the atmospheric circulation by shifting the latitudinal positions of the mid-latitude westerly jet and the Intertropical Convergence Zone (ITCZ). In quadrupled CO 2 experiments, we studied how air-sea coupling affects the response of tropical rainfall under global warming. In order to identify the importance of individual ocean processes, we used the hierarchy of ocean models to separate the effects of seasonal mixed-layer entrainment, wind-driven Ekman flows, and frictional flows. We showed that including Ekman and frictional flows allows our simulation to produce the Pacific Ocean's enhanced equatorial warming pattern and equatorward ITCZ contraction noted in previous climate simulations. We also showed that the frictional flow, which has yet to receive much attention, is as important as the Ekman flow in generating equatorial heat convergence.

54 ENVIRONMENTAL SCIENCES↗

Characteristics of Future Warmer Base States in CESM2

Abstract Simulations of 21st century climate with Community Earth System Model version 2 (CESM2) using the standard atmosphere (CAM6), denoted CESM2(CAM6), and the latest generation of the Whole Atmosphere Community Climate Model (WACCM6), denoted CESM2(WACCM6), are presented, and a survey of general results is described. The equilibrium climate sensitivity (ECS) of CESM2(CAM6) is 5.3°C, and CESM2(WACCM6) is 4.8°C, while the transient climate response (TCR) is 2.1°C in CESM2(CAM6) and 2.0°C in CESM2(WACCM6). Thus, these two CESM2 model versions have higher values of ECS than the previous generation of model, the CESM (CAM5) (hereafter CESM1), that had an ECS of 4.1°C, though the CESM2 versions have lower values of TCR compared to the CESM1 with a somewhat higher value of 2.3°C. All model versions produce credible simulations of the time evolution of historical global surface temperature. The higher ECS values for the CESM2 versions are reflected in higher values of global surface temperature increase by 2,100 in CESM2(CAM6) and CESM2(WACCM6) compared to CESM1 between comparable emission scenarios for the high forcing scenario. Future warming among CESM2 model versions and scenarios diverges around 2050. The larger values of TCR and ECS in CESM2(CAM6) compared to CESM1 are manifested by greater warming in the tropics. Associated with a higher climate sensitivity, for CESM2(CAM6) the first instance of an ice‐free Arctic in September occurs for all scenarios and ensemble members in the 2030–2050 time frame, but about a decade later in CESM2(WACCM6), occurring around 2040–2060.

54 ENVIRONMENTAL SCIENCES↗

Prediction of the Arctic Oscillation in Boreal Winter by Dynamical Seasonal Forecasting Systems

This study assesses the prediction skill of the boreal winter Arctic Oscillation (AO) in the state-of-the-art dynamical ensemble prediction systems (EPSs): the UKMO GloSea4, the NCEP CFSv2, and the NASA GEOS-5. Long-term reforecasts made with the EPSs are used to evaluate representations of the AO, and to examine skill scores for the deterministic and probabilistic forecast of the AO index. The reforecasts reproduce the observed changes in the large-scale patterns of the Northern Hemispheric surface temperature, upper-level wind, and precipitation according to the AO phase. Results demonstrate that all EPSs have better prediction skill than the persistence prediction for lead times up to 3-month, suggesting a great potential for skillful prediction of the AO and the associated climate anomalies in seasonal time scale. It is also found that the deterministic and probabilistic forecast skill of the AO in the recent period (1997-2010) is higher than that in the earlier period (1983-1996).

Coupled Models↗

Coupled Data Assimilation for Integrated Earth System Analysis and Prediction: Goals, Challenges, and Recommendations

The purpose of this report is to identify fundamental issues for coupled data assimilation (CDA), such as gaps in science and limitations in forecasting systems, in order to provide guidance to the World Meteorological Organization (WMO) on how to facilitate more rapid progress internationally. Coupled Earth system modeling provides the opportunity to extend skillful atmospheric forecasts beyond the traditional two-week barrier by extracting skill from low-frequency state components such as the land, ocean, and sea ice. More generally, coupled models are needed to support seamless prediction systems that span timescales from weather, subseasonal to seasonal (S2S), multiyear, and decadal. Therefore, initialization methods are needed for coupled Earth system models, either applied to each individual component (called Weakly Coupled Data Assimilation - WCDA) or applied the coupled Earth system model as a whole (called Strongly Coupled Data Assimilation - SCDA). Using CDA, in which model forecasts and potentially the state estimation are performed jointly, each model domain benefits from observations in other domains either directly using error covariance information known at the time of the analysis (SCDA), or indirectly through flux interactions at the model boundaries (WCDA). Because the non-atmospheric domains are generally under-observed compared to the atmosphere, CDA provides a significant advantage over single-domain analyses. Next, we provide a synopsis of goals, challenges, and recommendations to advance CDA: Goals: (a) Extend predictive skill beyond the current capability of NWP (e.g. as demonstrated by improving forecast skill scores), (b) produce physically consistent initial conditions for coupled numerical prediction systems and reanalyses (including consistent fluxes at the domain interfaces), (c) make best use of existing observations by allowing observations from each domain to influence and improve the full earth system analysis, (d) develop a robust observation-based identification and understanding of mechanisms that determine the variability of weather and climate, (e) identify critical weaknesses in coupled models and the earth observing system, (f) generate full-field estimates of unobserved or sparsely observed variables, (g) improve the estimation of the external forcings causing changes to climate, (h) transition successes from idealized CDA experiments to real-world applications. Challenges: (a) Modeling at the interfaces between interacting components of coupled Earth system models may be inadequate for estimating uncertainty or error covariances between domains, (b) current data assimilation methods may be insufficient to simultaneously analyze domains containing multiple spatiotemporal scales of interest, (c) there is no standardization of observation data or their delivery systems across domains, (d) the size and complexity of many large-scale coupled Earth system models makes it is difficult to accurately represent uncertainty due to model parameters and coupling parameters, (e) model errors lead to local biases that can transfer between the different Earth system components and lead to coupled model biases and long-term model drift, (e) information propagation across model components with different spatiotemporal scales is extremely complicated, and must be improved in current coupled modeling frameworks, (h) there is insufficient knowledge on how to represent evolving errors in non-atmospheric model components (e.g. as sea ice, land and ocean) on the timescales of NWP.

CDA↗

A11I-2104 Evaluation of the NASA GEOS Chemistry-Climate Model Coupled Atmosphere-Ocean Configuration for its Suitability to Simulations of Atmospheric Composition in a Changing Climate

The NASA Goddard Earth Observing System Chemistry-Climate Model (GEOSCCM) is a full tropospheric-stratospheric chemistry enabled configuration of the NASA GEOS Earth system model. Among the objectives for research performed with GEOSCCM is to understand the impacts of climate change on recovery of the stratospheric ozone layer over the twenty-first century. Simulations performed in the past with GEOSCCM did not have an interactive ocean and were driven with projected sea surface temperature and sea ice boundary conditions provided by external models. This method limits the application of GEOSCCM for climate study because it doesn’t allow the feedback of composition changes to ocean. For a full representation of the chemistry-climate feedbacks in the Earth System, we have developed a configuration of GEOSCCM that makes use of the MOM5 ocean general circulation model coupled to the GEOS atmospheric general circulation model. This configuration leverages development of the GEOS Sub-seasonal-to-Seasonal (S2S) prediction system, and includes interactive, radiatively coupled aerosols and two-moment, aerosol-aware cloud microphysics. W e have assessed the baseline model climate sensitivity with an Monday, 11 December 2023 08:30 - 12:50 Poster Hall A-C - South (Exhibition Level, South, MC) experiment performed under pre-industrial (e.g., year 1850) greenhouse gas conditions and a second experiment under 4xCO2 conditions, similar to the typical CMIP protocol. The GEOS model is demonstrated to have a climate sensitivity of a 2.6 K increase in global mean surface temperature for an equivalent doubling of CO2 from pre-industrial conditions, in line with current CMIP models. This configuration of the GEOS model is suitable for application to multidecadal to century-long simulations of climate system response to changing greenhouse gas levels. We report here our evaluation of the climate diagnostics of this configuration and discuss future directions for work with this model configuration, including an ongoing twenty-first century projection experiment based on the Chemistry-Climate Model Intercomparison project protocol.

models↗