Explore the vast range of titles available.

Recent research highlights the role of land surface processes in heat waves, droughts, and other extreme events. Here we use an earth system model (ESM) from the Geophysical Fluid Dynamics Laboratory (GFDL) to investigate the regional impacts of historical anthropogenic land useland cover change (LULCC) on combined extremes of temperature and humidity. A bivariate assessment allows us to consider aridity and moist enthalpy extremes, quantities central to human experience of near-surface climate conditions. We show that according to this model, conversion of forests to cropland has contributed to much of the upper central US and central Europe experiencing extreme hot, dry summers every 2-3 years instead of every 10 years. In the tropics, historical patterns of wood harvesting, shifting cultivation and regrowth of secondary vegetation have enhanced near surface moist enthalpy, leading to extensive increases in the occurrence of humid conditions throughout the tropics year round. These critical land use processes and practices are not included in many current generation land models, yet these results identify them as critical factors in the energy and water cycles of the midlatitudes and tropics.

The authors assess the uptake, transport, and storage of oceanic anthropogenic carbon and heat over the period 1861–2005 in a new set of coupled carbon–climate Earth system models conducted for the fifth phase of the Coupled Model Intercomparison Project (CMIP5), with a particular focus on the Southern Ocean. Simulations show that the Southern Ocean south of 30°S, occupying 30% of global surface ocean area, accounts for 43% ± 3% (42 ± 5 Pg C) of anthropogenic CO₂ and 75% ± 22% (23 ± 9 × 1022J) of heat uptake by the ocean over the historical period. Northward transport out of the Southern Ocean is vigorous, reducing the storage to 33 ± 6 Pg anthropogenic carbon and 12 ± 7 × 1022J heat in the region. The CMIP5 models, as a class, tend to underestimate the observation-based global anthropogenic carbon storage but simulate trends in global ocean heat storage over the last 50 years within uncertainties of observation-based estimates. CMIP5 models suggest global and Southern Ocean CO₂ uptake have been largely unaffected by recent climate variability and change. Anthropogenic carbon and heat storage show a common broad-scale pattern of change, but ocean heat storage is more structured than ocean carbon storage. The results highlight the significance of the Southern Ocean for the global climate and as the region where models differ the most in representation of anthropogenic CO₂ and, in particular, heat uptake.

Using simulations performed with 18 coupled atmosphere–ocean global climate models from phase 5 of the Coupled Model Intercomparison Project (CMIP5), projections of the Northern Hemisphere snowfall under the representative concentration pathway (RCP4.5) scenario are analyzed for the period 2006–2100. These models perform well in simulating twentieth-century snowfall, although there is a positive bias in many regions. Annual snowfall is projected to decrease across much of the Northern Hemisphere during the twenty-first century, with increases projected at higher latitudes. On a seasonal basis, the transition zone between negative and positive snowfall trends corresponds approximately to the −10°C isotherm of the late twentieth-century mean surface air temperature, such that positive trends prevail in winter over large regions of Eurasia and North America. Redistributions of snowfall throughout the entire snow season are projected to occur—even in locations where there is little change in annual snowfall. Changes in the fraction of precipitation falling as snow contribute to decreases in snowfall across most Northern Hemisphere regions, while changes in total precipitation typically contribute to increases in snowfall. A signal-to-noise analysis reveals that the projected changes in snowfall, based on the RCP4.5 scenario, are likely to become apparent during the twenty-first century for most locations in the Northern Hemisphere. The snowfall signal emerges more slowly than the temperature signal, suggesting that changes in snowfall are not likely to be early indicators of regional climate change.

This contribution describes the ocean biogeochemical component of the Geophysical Fluid Dynamics Laboratory's Earth System Model 4.1 (GFDL‐ESM4.1), assesses GFDL‐ESM4.1's capacity to capture observed ocean biogeochemical patterns, and documents its response to increasing atmospheric CO2. Notable differences relative to the previous generation of GFDL ESM's include enhanced resolution of plankton food web dynamics, refined particle remineralization, and a larger number of exchanges of nutrients across Earth system components. During model spin‐up, the carbon drift rapidly fell below the 10 Pg C per century equilibration criterion established by the Coupled Climate‐Carbon Cycle Model Intercomparison Project (C4MIP). Simulations robustly captured large‐scale observed nutrient distributions, plankton dynamics, and characteristics of the biological pump. The model overexpressed phosphate limitation and open ocean hypoxia in some areas but still yielded realistic surface and deep carbon system properties, including cumulative carbon uptake since preindustrial times and over the last decades that is consistent with observation‐based estimates. The model's response to the direct and radiative effects of a 200% atmospheric CO2 increase from preindustrial conditions (i.e., years 101–120 of a 1% CO2 yr−1 simulation) included (a) a weakened, shoaling organic carbon pump leading to a 38% reduction in the sinking flux at 2,000 m; (b) a two‐thirds reduction in the calcium carbonate pump that nonetheless generated only weak calcite compensation on century time‐scales; and, in contrast to previous GFDL ESMs, (c) a moderate reduction in global net primary production that was amplified at higher trophic levels. We conclude with a discussion of model limitations and priority developments. Plain Language Summary This paper describes and evaluates the ocean biogeochemical component of the Geophysical Fluid Dynamics Laboratory's Earth System Model 4.1 (GFDL‐ESM4.1). GFDL‐ESM4.1 was developed to study the past, present, and future evolution of the Earth system under scenarios for natural and anthropogenic drivers of Earth system change, including greenhouse gases and aerosols. The response of the ocean's vast carbon and heat reservoirs to accumulating greenhouse gases greatly reduces their atmospheric and terrestrial impacts, but also puts ocean environments and the marine resources they support at risk. Relative to previous models, GFDL‐ESM4.1 improves the representation of (a) ocean food webs connecting plankton and fish; (b) biological processes influencing the sequestration of carbon in the deep ocean; and (c) land‐atmosphere‐ocean nutrient exchanges. While simulations have biases, they capture many critical aspects of the global ocean carbon cycle and ocean ecosystem, including the observed uptake of anthropogenic carbon over the last ~150 yr. Projections suggest that continued CO2 increases could significantly decrease ocean productivity and the ocean's capacity to sequester atmospheric carbon. Key Points Enhanced plankton food webs, remineralization, and Earth system linkages yield skillful global carbon cycle and ecosystem simulations Ocean productivity estimates improved, but phosphate limitation and hypoxia overestimated in some areas High CO2 and associated warming substantially reduce the biological pump and ocean productivity across trophic levels

We document the configuration and emergent simulation features from the Geophysical Fluid Dynamics Laboratory (GFDL) OM4.0 ocean/sea ice model. OM4 serves as the ocean/sea ice component for the GFDL climate and Earth system models. It is also used for climate science research and is contributing to the Coupled Model Intercomparison Project version 6 Ocean Model Intercomparison Project. The ocean component of OM4 uses version 6 of the Modular Ocean Model and the sea ice component uses version 2 of the Sea Ice Simulator, which have identical horizontal grid layouts (Arakawa C‐grid). We follow the Coordinated Ocean‐sea ice Reference Experiments protocol to assess simulation quality across a broad suite of climate‐relevant features. We present results from two versions differing by horizontal grid spacing and physical parameterizations: OM4p5 has nominal 0.5° spacing and includes mesoscale eddy parameterizations and OM4p25 has nominal 0.25° spacing with no mesoscale eddy parameterization. Modular Ocean Model version 6 makes use of a vertical Lagrangian‐remap algorithm that enables general vertical coordinates. We show that use of a hybrid depth‐isopycnal coordinate reduces the middepth ocean warming drift commonly found in pure z* vertical coordinate ocean models. To test the need for the mesoscale eddy parameterization used in OM4p5, we examine the results from a simulation that removes the eddy parameterization. The water mass structure and model drift are physically degraded relative to OM4p5, thus supporting the key role for a mesoscale closure at this resolution. Key Points Documentation is provided for a new generation of NOAA‐GFDL CMIP6/OMIP ocean ice climate models Dynamical core and physical parameterizations are described and key features of interannual CORE simulations are assessed Using hybrid vertical coordinates reduces spurious ocean heat drift

The authors describe carbon system formulation and simulation characteristics of two new global coupled carbon–climate Earth System Models (ESM), ESM2M and ESM2G. These models demonstrate good climate fidelity as described in part I of this study while incorporating explicit and consistent carbon dynamics. The two models differ almost exclusively in the physical ocean component; ESM2M uses the Modular Ocean Model version 4.1 with vertical pressure layers, whereas ESM2G uses generalized ocean layer dynamics with a bulk mixed layer and interior isopycnal layers. On land, both ESMs include a revised land model to simulate competitive vegetation distributions and functioning, including carbon cycling among vegetation, soil, and atmosphere. In the ocean, both models include new biogeochemical algorithms including phytoplankton functional group dynamics with flexible stoichiometry. Preindustrial simulations are spun up to give stable, realistic carbon cycle means and variability. Significant differences in simulation characteristics of these two models are described. Because of differences in oceanic ventilation rates, ESM2M has a stronger biological carbon pump but weaker northward implied atmospheric CO₂ transport than ESM2G. The major advantages of ESM2G over ESM2M are improved representation of surface chlorophyll in the Atlantic and Indian Oceans and thermocline nutrients and oxygen in the North Pacific. Improved tree mortality parameters in ESM2G produced more realistic carbon accumulation in vegetation pools. The major advantages of ESM2M over ESM2G are reduced nutrient and oxygen biases in the southern and tropical oceans.

The physical climate formulation and simulation characteristics of two new global coupled carbon–climate Earth System Models, ESM2M and ESM2G, are described. These models demonstrate similar climate fidelity as the Geophysical Fluid Dynamics Laboratory’s previous Climate Model version 2.1 (CM2.1) while incorporating explicit and consistent carbon dynamics. The two models differ exclusively in the physical ocean component; ESM2M uses Modular Ocean Model version 4p1 with vertical pressure layers while ESM2G uses Generalized Ocean Layer Dynamics with a bulk mixed layer and interior isopycnal layers. Differences in the ocean mean state include the thermocline depth being relatively deep in ESM2M and relatively shallow in ESM2G compared to observations. The crucial role of ocean dynamics on climate variability is highlighted in El Niño–Southern Oscillation being overly strong in ESM2M and overly weak in ESM2G relative to observations. Thus, while ESM2G might better represent climate changes relating to total heat content variability given its lack of long-term drift, gyre circulation, and ventilation in the North Pacific, tropical Atlantic, and Indian Oceans, and depth structure in the overturning and abyssal flows, ESM2M might better represent climate changes relating to surface circulation given its superior surface temperature, salinity, and height patterns, tropical Pacific circulation and variability, and Southern Ocean dynamics. The overall assessment is that neither model is fundamentally superior to the other, and that both models achieve sufficient fidelity to allow meaningful climate and earth system modeling applications. This affords the ability to assess the role of ocean configuration on earth system interactions in the context of two state-of-the-art coupled carbon–climate models.

Previous studies have demonstrated the importance of enhanced vegetation growth under future elevated atmospheric CO ₂ for 21st century climate warming. Surprisingly no study has completed an analogous assessment for the historical period, during which emissions of greenhouse gases increased rapidly and land-use changes (LUC) dramatically altered terrestrial carbon sources and sinks. Using the Geophysical Fluid Dynamics Laboratory comprehensive Earth System Model ESM2G and a reconstruction of the LUC, we estimate that enhanced vegetation growth has lowered the historical atmospheric CO ₂ concentration by 85 ppm, avoiding an additional 0.31 ± 0.06 °C warming. We demonstrate that without enhanced vegetation growth the total residual terrestrial carbon flux (i.e., the net land flux minus LUC flux) would be a source of 65–82 Gt of carbon (GtC) to atmosphere instead of the historical residual carbon sink of 186–192 GtC, a carbon saving of 251–274 GtC.

The diurnal cycle of precipitation and precipitation variances at different time scales are analyzed in this study based on multiple high-resolution 3-h precipitation datasets. The results are used to evaluate nine CMIP6 models and a series of GFDL-AM4.0 model simulations, with the goal of examining the impact of SST diurnal cycle, varying horizontal resolutions, and different microphysics schemes on these two precipitation features. It is found that although diurnal amplitudes are reasonably simulated, models generally generate too early diurnal peaks over land, with a diurnal phase peaking around noon instead of the observed late afternoon (or early evening) peak. As for precipitation variances, irregular subdaily fluctuations dominate the total variance, followed by variance of daily mean precipitation and variance associated with the mean diurnal cycle. While the spatial and zonal distributions of precipitation variances are generally captured by the models, significant biases are present in tropical regions, where large mean precipitation biases are observed. The comparisons based on AM4.0 model simulations demonstrate that the inclusion of ocean coupling, adoption of a new microphysics scheme, and increasing of horizontal resolution have limited impacts on these two simulated features, emphasizing the need for future investigation into these model deficiencies at the process level. Conducting routine examinations of these metrics would be a crucial first step toward better simulation of precipitation intermittence in future model development. Last, distinct differences in these two features are found among observational datasets, highlighting the urgent need for a detailed evaluation of precipitation observations, especially at subdaily time scales, as model evaluation heavily relies on high-quality observations.

Results from the fully and biogeochemically coupled simulations in which CO2 increases at a rate of 1 % yr−1 (1pctCO2) from its preindustrial value are analyzed to quantify the magnitude of carbon–concentration and carbon–climate feedback parameters which measure the response of ocean and terrestrial carbon pools to changes in atmospheric CO2 concentration and the resulting change in global climate, respectively. The results are based on 11 comprehensive Earth system models from the most recent (sixth) Coupled Model Intercomparison Project (CMIP6) and compared with eight models from the fifth CMIP (CMIP5). The strength of the carbon–concentration feedback is of comparable magnitudes over land (mean ± standard deviation = 0.97 ± 0.40 PgC ppm−1) and ocean (0.79 ± 0.07 PgC ppm−1), while the carbon–climate feedback over land (−45.1 ± 50.6 PgC ∘C−1) is about 3 times larger than over ocean (−17.2 ± 5.0 PgC ∘C−1). The strength of both feedbacks is an order of magnitude more uncertain over land than over ocean as has been seen in existing studies. These values and their spread from 11 CMIP6 models have not changed significantly compared to CMIP5 models. The absolute values of feedback parameters are lower for land with models that include a representation of nitrogen cycle. The transient climate response to cumulative emissions (TCRE) from the 11 CMIP6 models considered here is 1.77 ± 0.37 ∘C EgC−1 and is similar to that found in CMIP5 models (1.63 ± 0.48 ∘C EgC−1) but with somewhat reduced model spread. The expressions for feedback parameters based on the fully and biogeochemically coupled configurations of the 1pctCO2 simulation are simplified when the small temperature change in the biogeochemically coupled simulation is ignored. Decomposition of the terms of these simplified expressions for the feedback parameters is used to gain insight into the reasons for differing responses among ocean and land carbon cycle models.