SEARCH · Search NASA

Results for “Earth System Models (ESMs)”

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.

SourceResearch connectionFrom yearThrough year

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

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS2023

Modeling Dust Mineralogical Composition: Sensitivity to Soil Mineralogy Atlases and Their Expected Climate Impacts

Soil dust aerosols are a key component of the climate system, as they interact with short- and long-wave radiation, alter cloud formation processes, affect atmospheric chemistry and play a role in biogeochemical cycles by providing nutrient inputs such as iron and phosphorus. The influence of dust on these processes depends on its physicochemical properties, which, far from being homogeneous, are shaped by its regionally varying mineral composition. The relative amount of minerals in dust depends on the source region and shows a large geographical variability. However, many state-of-the-art Earth system models (ESMs), upon which climate analyses and projections rely, still consider dust mineralogy to be invariant. The explicit representation of minerals in ESMs is more hindered by our limited knowledge of the global soil composition along with the resulting size-resolved airborne mineralogy than by computational constraints. In this work we introduce an explicit mineralogy representation within the state-of-the-art Multiscale Online Nonhydrostatic AtmospheRe CHemistry (MONARCH) model. We review and compare two existing soil mineralogy datasets, which remain a source of uncertainty for dust mineralogy modeling and provide an evaluation of multiannual simulations against available mineralogy observations. Soil mineralogy datasets are based on measurements performed after wet sieving, which breaks the aggregates found in the parent soil. Our model predicts the emitted particle size distribution (PSD) in terms of its constituent minerals based on brittle fragmentation theory (BFT), which reconstructs the emitted mineral aggregates destroyed by wet sieving. Our simulations broadly reproduce the most abundant mineral fractions independently of the soil composition data used. Feldspars and calcite are highly sensitive to the soil mineralogy map, mainly due to the different assumptions made in each soil dataset to extrapolate a handful of soil measurements to arid and semi-arid regions worldwide. For the least abundant or more difficult-to-determine minerals, such as iron oxides, uncertainties in soil mineralogy yield differences in annual mean aerosol mass fractions of up to ∼ 100 %. Although BFT restores coarse aggregates including phyllosilicates that usually break during soil analysis, we still identify an overestimation of coarse quartz mass fractions (above 2 µm in diameter). In a dedicated experiment, we estimate the fraction of dust with undetermined composition as given by a soil map, which makes up ∼ 10 % of the emitted dust mass at the global scale and can be regionally larger. Changes in the underlying soil mineralogy impact our estimates of climate-relevant variables, particularly affecting the regional variability of the single-scattering albedo at solar wavelengths or the total iron deposited over oceans. All in all, this assessment represents a baseline for future model experiments including new mineralogical maps constrained by high-quality spaceborne hyperspectral measurements, such as those arising from the NASA Earth Surface Mineral Dust Source Investigation (EMIT) mission.

Soil dust↗

NASA NTRS2012

Global Modeling, Field Campaigns, Upscaling and Ray Desjardins

In the early 1980's, it became apparent that land surface radiation and energy budgets were unrealistically represented in Global Circulation models (GCM's), Shortly thereafter, it became clear that the land carbon budget was also poorly represented in Earth System Models (ESM's), A number of scientific communities, including GCM/ESM modelers, micrometeorologists, satellite data specialists and plant physiologists, came together to design field experiments that could be used to develop and validate the contemporary prototype land surface models. These experiments were designed to measure land surface fluxes of radiation, heat, water vapor and CO2 using a network of flux towers and other plot-scale techniques, coincident with satellite measurements of related state variables, The interdisciplinary teams involved in these experiments quickly became aware of the scale gap between plot-scale measurements (approx 10 - 100m), satellite measurements (100m - 10 km), and GCM grid areas (l0 - 200km). At the time, there was no established flux measurement capability to bridge these scale gaps. Then, a Canadian science learn led by Ray Desjardins started to actively participate in the design and execution of the experiments, with airborne eddy correlation providing the radically innovative bridge across the scale gaps, In a succession of brilliantly executed field campaigns followed up by convincing scientific analyses, they demonstrated that airborne eddy correlation allied with satellite data was the most powerful upscaling tool available to the community, The rest is history: the realism and credibility of weather and climate models has been enormously improved enormously over the last 25 years with immense benefits to the public and policymakers.

Sellers, P. J.↗

NASA NTRS2014

HEMCO v1.0: A Versatile, ESMF-Compliant Component for Calculating Emissions in Atmospheric Models

We describe the Harvard-NASA Emission Component version 1.0 (HEMCO), a stand-alone software component for computing emissions in global atmospheric models. HEMCO determines emissions from different sources, regions, and species on a user-defined grid and can combine, overlay, and update a set of data inventories and scale factors, as specified by the user through the HEMCO configuration file. New emission inventories at any spatial and temporal resolution are readily added to HEMCO and can be accessed by the user without any preprocessing of the data files or modification of the source code. Emissions that depend on dynamic source types and local environmental variables such as wind speed or surface temperature are calculated in separate HEMCO extensions. HEMCO is fully compliant with the Earth System Modeling Framework (ESMF) environment. It is highly portable and can be deployed in a new model environment with only few adjustments at the top-level interface. So far, we have implemented HEMCO in the NASA Goddard Earth Observing System (GEOS-5) Earth system model (ESM) and in the GEOS-Chem chemical transport model (CTM). By providing a widely applicable framework for specifying constituent emissions, HEMCO is designed to ease sensitivity studies and model comparisons, as well as inverse modeling in which emissions are adjusted iteratively. The HEMCO code, extensions, and the full set of emissions data files used in GEOS-Chem are available at http: //wiki.geos-chem.org/HEMCO.

GEOS-5↗

NASA NTRS2015

Inconsistent Strategies to Spin up Models in CMIP5: Implications for Ocean Biogeochemical Model Performance Assessment

During the fifth phase of the Coupled Model Intercomparison Project (CMIP5) substantial efforts were made to systematically assess the skill of Earth system models. One goal was to check how realistically representative marine biogeochemical tracer distributions could be reproduced by models. In routine assessments model historical hindcasts were compared with available modern biogeochemical observations. However, these assessments considered neither how close modeled biogeochemical reservoirs were to equilibrium nor the sensitivity of model performance to initial conditions or to the spin-up protocols. Here, we explore how the large diversity in spin-up protocols used for marine biogeochemistry in CMIP5 Earth system models (ESMs) contributes to model-to-model differences in the simulated fields. We take advantage of a 500-year spin-up simulation of IPSL-CM5A-LR to quantify the influence of the spin-up protocol on model ability to reproduce relevant data fields. Amplification of biases in selected biogeochemical fields (O2, NO3, Alk-DIC) is assessed as a function of spin-up duration. We demonstrate that a relationship between spin-up duration and assessment metrics emerges from our model results and holds when confronted with a larger ensemble of CMIP5 models. This shows that drift has implications for performance assessment in addition to possibly aliasing estimates of climate change impact. Our study suggests that differences in spin-up protocols could explain a substantial part of model disparities, constituting a source of model-to- model uncertainty. This requires more attention in future model intercomparison exercises in order to provide quantitatively more correct ESM results on marine biogeochemistry and carbon cycle feedbacks.

histories↗

NASA NTRS2017

The Impact of Anthropogenic Land Use and Land Cover Change on Regional Climate Extremes

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.

Near-surface↗

NASA NTRS2020

The interactive global fire module pyrE (v1.0)

Fires affect the composition of the atmosphere and Earth’s radiation balance by emitting a suite of reactive gases and particles. An interactive fire module in an Earth system model (ESM) allows us to study the natural and anthropogenic drivers, feedbacks, and interactions of open fires. To do so, we have developed pyrE, the NASA GISS (Goddard Institute for Space Studies) interactive fire emissions *module. The pyrE module is driven by environmental variables like flammability and cloud-to-ground lightning, calculated by the GISS ModelE ESM, and parameterized by anthropogenic impacts based on population density data. Fire emissions are generated from the flaming phase in pyrE (active fires). Using pyrE, we examine fire occurrence, regional fire suppression, burned area, fire emissions, and how it all affects atmospheric composition. To do so, we evaluate pyrE by comparing it to satellite-based datasets of fire count, burned area, fire emissions, and aerosol optical depth (AOD). We demonstrate pyrE’s ability to simulate the daily and seasonal cycles of open fires and resulting emissions. Our results indicate that interactive fire emissions are biased low by 32 %–42 %, depending on emitted species, compared to the GFED4s (Global Fire Emissions Database) inventory. The bias in emissions drives underestimation in column densities, which is diluted by natural and anthropogenic emissions sources and production and loss mechanisms. Regionally, the resulting AOD of a simulation with interactive fire emissions is underestimated mostly over Indonesia compared to a simulation with GFED4s emissions and to MODIS AOD. In other parts of the world pyrE’s performance in terms of AOD is marginal to a simulation with prescribed fire emissions.

Fires↗

NASA NTRS2020

Fast Responses on Pre-Industrial Climate from Present-Day Aerosols in A CMIP6 Multi-Model Study

In this work, we use Coupled Model Intercomparison Project Phase 6 (CMIP6) simulations from 10 Earth system models (ESMs) and general circulation models (GCMs) to study the fast climate responses on pre-industrial climate, due to present-day aerosols. All models carried out two sets of simulations: a control experiment with all forcings set to the year 1850 and a perturbation experiment with all forcings identical to the control, except for aerosols with precursor emissions set to the year 2014. In response to the pattern of all aerosols effective radiative forcing (ERF), the fast temperature responses are characterized by cooling over the continental areas, especially in the Northern Hemisphere, with the largest cooling over East Asia and India, sulfate being the dominant aerosol surface temperature driver for present-day emissions. In the Arctic there is a warming signal for winter in the ensemble mean of fast temperature responses, but the model-to-model variability is large, and it is presumably linked to aerosol-induced circulation changes. The largest fast precipitation responses are seen in the tropical belt regions, generally characterized by a reduction over continental regions and presumably a southward shift of the tropical rain belt. This is a characteristic and robust feature among most models in this study, associated with weakening of the monsoon systems around the globe (Asia, Africa and America) in response to hemispherically asymmetric cooling from a Northern Hemisphere aerosol perturbation, forcing possibly the Intertropical Convergence Zone (ITCZ) and tropical precipitation to shift away from the cooled hemisphere despite that aerosols' effects on temperature and precipitation are only partly realized in these simulations as the sea surface temperatures are kept fixed. An interesting feature in aerosol-induced circulation changes is a characteristic dipole pattern with intensification of the Icelandic Low and an anticyclonic anomaly over southeastern Europe, inducing warm air advection towards the northern polar latitudes in winter.

Pre-industrial climate↗

NASA NTRSDate not supplied

Regional variability of dust single scattering albedo due to mineral composition

Nearly all Earth System Models (ESMs) assume globally homogeneous dust aerosols, neglecting regional variations of the imaginary refractive index (IRI) due to varying mineral composition. This has led to a range of single scattering albedo (SSA) and direct radiative forcing (DRF) estimates, as models assign global properties using dust measurements from different regions. We use model and observational data to assess to what extent regionally varying mineral composition affects visible-band SSA and short-wave DRF of dust. We run global simulations with NASA GISS ModelE2.1, using optical properties for minerals based on laboratory-derived empirical relationships between dust IRI at visible wavelengths and the mineral content of iron oxides. When allowing mineral variations instead of homogeneous dust, we find regional differences in dust SSA up to ~0.06 and consequent variations in DRF up to ~6 and ~5 W/m2, at surface and top-of-atmosphere respectively. We compare model SSA with AERONET inversion data (Version 3.0, Level 2), filtered by aerosol size and optical properties to identify pure dust scenes. To investigate possible contamination by biomass burning aerosols, we use the dataset of Schuster et al. (2016, doi:10.5194/acp-16-1565-2016), who separately calculated the contribution of iron oxides and carbonaceous species to AERONET extinction and absorption optical depths. Our results show that: (1) the range of observed SSA (even without carbonaceous species) is larger than that resulting from homogeneous dust; (2) residual amounts of fine mode black and brown carbon may affect AERONET SSA even in seasons with limited biomass burning; (3) model SSA from our mineral scheme exceeds the AERONET range, possibly due also to an uncertain treatment of goethite. In summary, homogeneous dust cannot explain the variability of observed SSA, so regionally varying optical properties based on mineral content are necessary. Higher accuracy in soil mineralogy maps is needed to better reproduce SSA variations at regional scales. Also, distinct soil maps for hematite and goethite are required, given the high sensitivity of model SSA to their extreme optical properties.

V. Obiso↗

NASA NTRSDate not supplied

NASA Earth: Synthetic Spectranomics - Deep Learning of Surface 3-D Geometry, Chemistry, and Hyperspectra to Inform Next-generation Land Models

Machine and deep learning (ML/DL) have transformed our approach to Earth observation and system modeling (EOSM), unifying both in view of ML/DL models as a form of data assimilation (DA). Trained on diverse Earth observation records, detailed physical models, or hybrids of both in physics-informed machine learning, ML/DL may improve upon existing Earth system model (ESM) formulations while creating entirely new classes of models. One important application of deep learning is observation synthesis, allowing ESM developers to prepare for the increased spatial, temporal, and spectral/polar resolution of proposed future observing systems. This may involve the retrospective application of learning algorithms to existing observational records, with or without physical radiative transfer models, to co-inform mission planning and ESM development while providing a degree of data continuity for new missions.

Adam Erickson↗

NASA NTRS2021

AeroCom Phase III Multi-Model Evaluation of the Aerosol Life Cycle and Optical Properties Using Ground and Space-Based Remote Sensing as Well as Surface In Situ Observations

Within the framework of the AeroCom (Aerosol Comparisons between Observations and Models) initiative, the state-of-the-art modelling of aerosol optical properties is assessed from 14 global models participating in the phase III control experiment (AP3). The models are similar to CMIP6/AerChemMIP Earth System Models (ESMs) and provide a robust multi-model ensemble. Inter-model spread of aerosol species lifetimes and emissions appears to be similar to that of mass extinction coefficients (MECs), suggesting that aerosol optical depth (AOD) uncertainties are associated with a broad spectrum of parameterised aerosol processes. Total AOD is approximately the same as in AeroCom phase I (AP1) simulations. However, we find a 50 % decrease in the optical depth (OD) of black carbon (BC), attributable to a combination of decreased emissions and lifetimes. Relative contributions from sea salt (SS) and dust (DU) have shifted from being approximately equal in AP1 to SS contributing about 2∕3 of the natural AOD in AP3. This shift is linked with a decrease in DU mass burden, a lower DU MEC, and a slight decrease in DU lifetime, suggesting coarser DU particle sizes in AP3 compared to AP1. Relative to observations, the AP3 ensemble median and most of the participating models underestimate all aerosol optical properties investigated, that is, total AOD as well as fine and coarse AOD (AODf, AODc), Ångström exponent (AE), dry surface scattering (SCdry), and absorption (ACdry) coefficients. Compared to AERONET, the models underestimate total AOD by ca. 21 % ± 20 % (as inferred from the ensemble median and interquartile range). Against satellite data, the ensemble AOD biases range from −37 % (MODIS-Terra) to −16 % (MERGED-FMI, a multi-satellite AOD product), which we explain by differences between individual satellites and AERONET measurements themselves. Correlation coefficients (R) between model and observation AOD records are generally high (R>0.75), suggesting that the models are capable of capturing spatio-temporal variations in AOD. We find a much larger underestimate in coarse AODc (∼ −45 % ± 25 %) than in fine AODf (∼ −15 % ± 25 %) with slightly increased inter-model spread compared to total AOD. These results indicate problems in the modelling of DU and SS. The AODc bias is likely due to missing DU over continental land masses (particularly over the United States, SE Asia, and S. America), while marine AERONET sites and the AATSR SU satellite data suggest more moderate oceanic biases in AODc. Column AEs are underestimated by about 10 % ± 16 %. For situations in which measurements show AE > 2, models underestimate AERONET AE by ca. 35 %. In contrast, all models (but one) exhibit large overestimates in AE when coarse aerosol dominates (bias ca. +140 % if observed AE < 0.5). Simulated AE does not span the observed AE variability. These results indicate that models overestimate particle size (or underestimate the fine-mode fraction) for fine-dominated aerosol and underestimate size (or overestimate the fine-mode fraction) for coarse-dominated aerosol. This must have implications for lifetime, water uptake, scattering enhancement, and the aerosol radiative effect, which we can not quantify at this moment. Comparison against Global Atmosphere Watch (GAW) in situ data results in mean bias and inter-model variations of −35 % ± 25 % and −20 % ± 18 % for SCdry and ACdry, respectively. The larger underestimate of SCdry than ACdry suggests the models will simulate an aerosol single scattering albedo that is too low. The larger underestimate of SCdry than ambient air AOD is consistent with recent findings that models overestimate scattering enhancement due to hygroscopic growth. The broadly consistent negative bias in AOD and surface scattering suggests an underestimate of aerosol radiative effects in current global aerosol models. Considerable inter-model diversity in the simulated optical properties is often found in regions that are, unfortunately, not or only sparsely covered by ground-based observations. This includes, for instance, the Sahara, Amazonia, central Australia, and the South Pacific. This highlights the need for a better site coverage in the observations, which would enable us to better assess the models, but also the performance of satellite products in these regions. Using fine-mode AOD as a proxy for present-day aerosol forcing estimates, our results suggest that models underestimate aerosol forcing by ca. −15 %, however, with a considerably large interquartile range, suggesting a spread between −35 % and +10 %.

space-based remote sensing↗

NASA NTRS2020

Bias in CMIP6 Models as Compared to Observed Regional Dimming and Brightening

Anthropogenic aerosol emissions have increased considerably over the last century, but climate effects and quantification of the emissions are highly uncertain as one goes back in time. This uncertainty is partly due to a lack of observations in the pre-satellite era, making the observations we do have before 1990 additionally valuable. Aerosols suspended in the atmosphere scatter and absorb incoming solar radiation and thereby alter the Earth's surface energy balance. Previous studies show that Earth system models (ESMs) do not adequately represent surface energy fluxes over the historical era. We investigated global and regional aerosol effects over the time period 1961–2014 by looking at surface downwelling shortwave radiation (SDSR). We used observations from ground stations as well as multiple experiments from eight ESMs participating in the Coupled Model Intercomparison Project Version 6 (CMIP6). Our results show that this subset of models reproduces the observed transient SDSR well in Europe but poorly in China. We suggest that this may be attributed to missing emissions of sulfur dioxide in China, sulfur dioxide being a precursor to sulfate, which is a highly reflective aerosol and responsible for more reflective clouds. The emissions of sulfur dioxide used in the models do not show a temporal pattern that could explain observed SDSR evolution over China. The results from various aerosol emission perturbation experiments from DAMIP, RFMIP and AerChemMIP show that only simulations containing anthropogenic aerosol emissions show dimming, even if the dimming is underestimated. Simulated clear-sky and all-sky SDSR do not differ greatly, suggesting that cloud cover changes are not a dominant cause of the biased SDSR evolution in the simulations. Therefore we suggest that the discrepancy between modeled and observed SDSR evolution is partly caused by erroneous aerosol and aerosol precursor emission inventories. This is an important finding as it may help interpret whether ESMs reproduce the historical climate evolution for the right or wrong reason.

surface downwelling shortwave radiation (SDSR)↗

NASA NTRS2020

Harmonization of global land use change and management for the period 850–2100 (LUH2) for CMIP6

Human land use activities have resulted in large changes to the biogeochemical and biophysical properties of the Earth's surface, with consequences for climate and other ecosystem services. In the future, land use activities are likely to expand and/or intensify further to meet growing demands for food, fiber, and energy. As part of the World Climate Research Program Coupled Model Intercomparison Project (CMIP6), the international community has developed the next generation of advanced Earth system models (ESMs) to estimate the combined effects of human activities (e.g., land use and fossil fuel emissions) on the carbon–climate system. A new set of historical data based on the History of the Global Environment database (HYDE), and multiple alternative scenarios of the future (2015–2100) from Integrated Assessment Model (IAM) teams, is required as input for these models. With most ESM simulations for CMIP6 now completed, it is important to document the land use patterns used by those simulations. Here we present results from the Land-Use Harmonization 2 (LUH2) project, which smoothly connects updated historical reconstructions of land use with eight new future projections in the format required for ESMs. The harmonization strategy estimates the fractional land use patterns, underlying land use transitions, key agricultural management information, and resulting secondary lands annually, while minimizing the differences between the end of the historical reconstruction and IAM initial conditions and preserving changes depicted by the IAMs in the future. The new approach builds on a similar effort from CMIP5 and is now provided at higher resolution (0.25°×0.25°) over a longer time domain (850–2100, with extensions to 2300) with more detail (including multiple crop and pasture types and associated management practices) using more input datasets (including Landsat remote sensing data) and updated algorithms (wood harvest and shifting cultivation); it is assessed via a new diagnostic package. The new LUH2 products contain > 50 times the information content of the datasets used in CMIP5 and are designed to enable new and improved estimates of the combined effects of land use on the global carbon–climate system.

land use land cover change↗

NASA NTRSDate not supplied

Variability in Modelled Airborne Dust Mineralogy Derived from Global Soil Composition Uncertainties

Dust aerosols consist on a variety of minerals with different physic and chemical properties. As such, they interact with short and long wave radiation, potentially form clouds, act as nutrients modulating biogeochemical cycles, or influence atmospheric chemistry, differently. Most current state-of-the-art Earth System Models (ESMs) neglect the complexity in dust composition, mainly due to computational constraints, but also to the existing uncertainties in the size resolved composition of parent soils, the resulting distribution of minerals in airborne dust, and the scarcity of observations to constrain them. Within this work, we assess the variability of global dust composition due to uncertainties in the characterization of the parent soil mineralogy. To that end, we consider two available global soil mineralogy atlases, developed by Claquin et al. (1999) –C1999- and Journet et al. (2014) –J2014-, which represent respectively 8 and 12 relevant minerals for climate (namely: illite, smectite/montmorillonite, kaolinite, calcite, gypsum, hematite, quartz, and feldspars in C1999, and those plus chlorite, vermiculite, goethite, and mica in J2014). Thanks to a recently developed feature of the MONARCH atmospheric-chemistry model, we are able to explicitly resolve the minerals’ atmospheric cycle. Therefore, we define two global experiments to assess changes on airborne dust composition attributed to the soil mineralogy assumptions and provide a measure of their variability. We also perform a preliminary evaluation of the global mineralogy results against available observations of mineral fractions in surface dust concentration. Our results will inform the climate modelling community about the potential variability in dust composition, an aspect that will gain relevance as ESMs continue growing in complexity and new processes to better characterize aerosols’ forcing or biogeochemical cycles are added. Further observational constraints, such as those that will derive from the EMIT NASA mission on soil composition or the FRAGMENT experimental campaigns on airborne dust characterization, will be key in the near future to improve our understanding of the impact of dust mineralogy on fundamental climate features.

Dust aerosols↗

NASA NTRS2021

Investigations on the Anthropogenic Reversal of the Natural Ozone Gradient Between Northern and Southern Midlatitudes

Our quantitative understanding of natural tropospheric ozone concentrations is limited by the paucity of reliable measurements before the 1980s. We utilize the existing measurements to compare the long-term ozone changes that occurred within the marine boundary layer at northern and southern midlatitudes. Since 1950 ozone concentrations have increased by a factor of 2.1 ± 0.2 in the Northern Hemisphere (NH) and are presently larger than in the Southern Hemisphere (SH), where only a much smaller increase has occurred. These changes are attributed to increased ozone production driven by anthropogenic emissions of photochemical ozone precursors that increased with industrial development. The greater ozone concentrations and increases in the NH are consistent with the predominant location of anthropogenic emission sources in that hemisphere. The available measurements indicate that this interhemispheric gradient was much smaller and was likely reversed in the pre-industrial troposphere with higher concentrations in the SH. Six Earth system model (ESM) simulations indicate similar total NH increases (1.9 with a standard deviation of 0.3), but they occurred more slowly over a longer time period, and the ESMs do not find higher pre-industrial ozone in the SH. Several uncertainties in the ESMs may cause these model–measurement disagreements: the assumed natural nitrogen oxide emissions may be too large, the relatively greater fraction of ozone injected by stratosphere–troposphere exchange to the NH may be overestimated, ozone surface deposition to ocean and land surfaces may not be accurately simulated, and model treatment of emissions of biogenic hydrocarbons and their photochemistry may not be adequate.

natural tropospheric ozone concentrations↗

NASA NTRS2021

Reduced Complexity Model Intercomparison Project Phase 2: Synthesizing Earth System Knowledge for Probabilistic Climate Projections

Over the last decades, climate science has evolved rapidly across multiple expert domains. Our best tools to capture state-of-the-art knowledge in an internally self-consistent modelling framework are the increasingly complex fully coupled Earth System Models (ESMs). However, computational limitations and the structural rigidity of ESMs mean that the full range of uncertainties across multiple domains are difficult to capture with ESMs alone. The tools of choice are instead more computationally efficient reduced complexity models (RCMs), which are structurally flexible and can span the response dynamics across a range of domain-specific models and ESM experiments. Here we present Phase 2 of the Reduced Complexity Model Intercomparison Project (RCMIP Phase 2), the first comprehensive intercomparison of RCMs that are probabilistically calibrated with key benchmark ranges from specialized research communities. Unsurprisingly, but crucially, we find that models which have been constrained to reflect the key benchmarks better reflect the key benchmarks. Under the low-emissions SSP1-1.9 scenario, across the RCMs, median peak warming projections range from 1.3 to 1.7°C (relative to 1850-1900, using an observationally-based historical warming estimate of 0.8°C between 1850-1900 and 1995-2014). Further developing methodologies to constrain these projection uncertainties seems paramount given the international community's goal to contain warming to below 1.5°C above pre-industrial in the long-term. Our findings suggest that users of RCMs should carefully evaluate their RCM, specifically its skill against key benchmarks and consider the need to include projections benchmarks either from ESM results or other assessments to reduce divergence in future projections.

Climate↗

NASA NTRS2022

The Simulation and Subseasonal Forecasting of Hydrological Variables: Insights from a Simple Water Balance Model

Past work has shown that a land surface model’s (LSM’s) implicit (not explicitly coded) relationships between soil moisture and both evapotranspiration (ET) and runoff largely determine the LSM’s hydrological behavior. Here we estimate the relationships that appear to be operating in the real world and compare them to those of the LSM component of a state-of-the-art Earth system model (ESM). The two sets of relationships are determined by calibrating them within a simple water balance model (WBM): once using stream gauge observations from small, unregulated rivers over the eastern half of the U.S., and once using the runoffs generated by the LSM as part of a state-of-the-art atmospheric reanalysis. Hydrological simulations and subseasonal hydrological forecasts performed with the two calibrated versions of the WBM provide two key results. First, the version calibrated to the LSM-generated runoffs does successfully reproduce, to first order, the hydrological behavior of the full LSM within its ESM environment. Second, of the two WBM versions, the one calibrated to the observations reproduces more accurately a broad collection of fully independent streamflow observations as well as a similarly broad collection of in-situ soil moisture measurements. Taken together, the two results suggest that the observations-calibrated ET and runoff efficiency functions do successfully represent, at least to some degree, soil moisture controls over hydrological variability in Nature and can serve as potentially useful targets for further LSM development.

Water Balance Model↗

NASA NTRS2023

How Does the Use of Different Soil Mineralogical Atlases Impact Soluble Iron Deposition Estimates?

Atmospheric deposition is the primary input of soluble iron (Fe) in remote regions of the open ocean. Fe constitutes a fundamental micro-nutrient for marine biota, affecting biogeochemical cycles and ultimately carbon uptake. Fe is emitted to the atmosphere primarily from dust sources. Fe solubility at emission is origin-dependent, being mostly insoluble when associated with dust minerals. Earth System Models (ESMs) commonly assume that dust aerosols have a globally uniform composition, neglecting known regional variations in soil mineralogy. This work assesses the implications of soil composition uncertainties on the bio-available Fe delivery to the ocean by using a state-of-the-art ESM with a detailed atmospheric Fe cycle (EC-Earth3-Iron). We run two 1-year simulations considering two different Soil Mineralogy Atlases (SMAs). Our results reveal a non-negligible impact of soil mineralogy uncertainties on soluble Fe estimates. By only changing the source of information for the soil composition and updating the Fe content in minerals, we find differences of 31.5% in mineral Fe emissions and 26.7% in soluble Fe deposition budgets.

soluble iron↗

NASA NTRS2023

Approaching a Thermal Tipping Point in the Eurasian Boreal Forest at its Southern Margin

Climate change is increasing the intensity and frequency of extreme heat events. Ecological responses to extreme heat will depend on vegetation physiology and thermal tolerance. Here we report that Larix sibirica, a foundation species across boreal Eurasia, is vulnerable to extreme heat at its southern range margin due to its low thermal tolerance (T crit of photosynthesis: ~ 37–48 °C). Projections from CMIP6 Earth System Models (ESMs) suggest that leaf temperatures might exceed the 25 th percentile of Larix sibirica’s T crit by two to three days per year within the next two to three decades (by 2050) under high emission scenarios (SSP3-7.0 and SSP5-8.5). This degree of warming will threaten the biome’s continued ability to assimilate and sequester carbon. This work highlights that under high emission trajectories we may approach an abrupt ecological tipping point in southern boreal Eurasian forests substantially sooner than ESM estimates that do not consider plant thermal tolerance traits.

Boreal ecology↗