Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Earth System Model”

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 1,171 records · Page 65

A Lagrangian perspective on tropical anvil cloud lifecycle in present and future climate

The evolution of tropical anvil clouds from their origin in deep convective cores to their slow decay determines the climatic effects of clouds in tropical convective regions. Despite the relevance of anvil clouds for climate and responses of clouds to global warming, processes dominating their evolution are not well understood. Currently available observational data reveal instantaneous snapshots of anvil cloud properties, but cannot provide a process-based perspective on anvil evolution. We therefore conduct simulations with the high resolution version of the Exascale Earth System Model in which we track mesoscale convective systems over the Tropical Western Pacific and compute trajectories that follow air parcels detrained from peaks of convective activity. With this approach we gain new insight into the anvil cloud evolution both in present day and future climate. Comparison with geostationary satellite data shows that the model is able to simulate maritime mesoscale convective systems reasonably well. Trajectory results indicate that anvil cloud lifetime is about 15 hours with no significant change in a warmer climate. The anvil ice mixing ratio is larger in a warmer climate due to a larger source of ice by detrainment and larger depositional growth leading to a more negative net cloud radiative effect along detrained trajectories. However, the increases in sources are counteracted by increases in sinks of ice, particularly snow formation and sedimentation. Furthermore, we find that the mean anvil cloud feedback along trajectories is positive and consistent with results from more traditional cloud feedback calculation methods.

climate, weather, mesocale cellular convection↗

Full Implementation of Matrix Approach to Biogeochemistry Module of CLM5

Abstract Earth system models (ESMs) have been rapidly developed in recent decades to advance our understanding of climate change‐carbon cycle feedback. However, those models are massive in coding, require expensive computational resources, and have difficulty in diagnosing their performance. It is highly desirable to develop ESMs with modularity and effective diagnostics. Toward these goals, we implemented a matrix approach to the Community Land Model version 5 (CLM5) to represent carbon and nitrogen cycles. Specifically, we reorganized 18 balance equations each for carbon and nitrogen cycles among the 18 vegetation pools in the original CLM5 into two matrix equations. Similarly, 140 balance equations each for carbon and nitrogen cycles among the 140 soil pools were reorganized into two additional matrix equations. The vegetation carbon and nitrogen matrix equations are connected to soil matrix equations via litterfall. The matrix equations fully reproduce simulations of carbon and nitrogen dynamics by the original model. The computational cost for forwarding simulation of the CLM5 matrix model was 26% more expensive than the original model, largely due to calculation of additional diagnostic variables, but the spin‐up computational cost was significantly saved. We showed a case study on modeled soil carbon storage under two forcing data sets to illustrate the diagnostic capability that the matrix approach uniquely offers to understand simulation results of global carbon and nitrogen dynamics. The successful implementation of the matrix approach to CLM5, one of the most complex land models, demonstrates that most, if not all, the biogeochemical models can be reorganized into the matrix form to gain high modularity, effective diagnostics, and accelerated spin‐up.

58 GEOSCIENCES↗

Assessing CLUBB PDF Closure Assumptions for a Continental Shallow-to-Deep Convective Transition Case Over Multiple Spatial Scales

Assumed-PDF (probability density function) higher-order turbulence closures (APHOCs) are now widely used for parameterizing boundary layer turbulence and shallow convection in Earth system models (ESMs). A better understanding of the resolution-dependent behavior of APHOCs is essential for improving the performance of next-generation ESMs with intended horizontal resolutions finer than 10 km. In this study, we evaluate the PDF family of Analytic double-Gaussian 1 implemented in Cloud Layers Unified By Binormals (CLUBB) over a range of spatial scales (Dx) from 2 to 100 km. A 120-km-wide large eddy simulation (LES) for a continental convection case during 2016 Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems (HI-SCALE) field campaign serves as benchmark to evaluate the PDF closure using an off-line approach. We find during the shallow convection period, the CLUBB PDF closure tends to produce positive biases of cloud properties and liquid water flux near cloud base for all scales of analysis. It produces negative biases for these variables near cloud top that are more severe for Dx larger than 25 km. Results show that replacing the CLUBB-parameterized moisture and temperature skewnesses with LES-derived ones can fix most of the biases if clipping of input moments is allowed to prevent the occurrence of unrealizable solutions. Overall, the performance of the PDF closure is better for smaller Dx = 2–5 km than for larger Dx = 50–100 km; for a given grid spacing, it is better when the convective clouds become deeper in the late afternoon. Likely causes for the resolution dependence and implications for improving the PDF closure are discussed.

54 ENVIRONMENTAL SCIENCES↗

A Baseline for Global Weather and Climate Simulations at 1 km Resolution

In an attempt to advance the understanding of the Earth's weather and climate by representing deep convection explicitly, we present a global, four-month simulation (November 2018 to February 2019) with ECMWF's hydrostatic Integrated Forecasting System (IFS) at an average grid spacing of 1.4 km. The impact of explicitly simulating deep convection on the atmospheric circulation and its variability is assessed by comparing the 1.4 km simulation to the equivalent well-tested and calibrated global simulations at 9 km grid spacing with and without parametrized deep convection. The explicit simulation of deep convection at 1.4 km results in a realistic large-scale circulation, better representation of convective storm activity, and stronger convective gravity wave activity when compared to the 9 km simulation with parametrized deep convection. Comparison of the 1.4 km simulation to the 9 km simulation without parametrized deep convection shows that switching off deep convection parametrization at a too coarse resolution (i.e., 9 km) generates too strong convective gravity waves. Based on the limited statistics available, improvements to the Madden-Julian Oscillation or tropical precipitation are not observed at 1.4 km, suggesting that other Earth system model components and/or their interaction are important for an accurate representation of these processes and may well need adjusting at deep convection resolving resolutions. Overall, the good agreement of the 1.4 km simulation with the 9 km simulation with parametrized deep convection is remarkable, despite one of the most fundamental parametrizations being turned off at 1.4 km resolution and despite no adjustments being made to the remaining parametrizations.

54 ENVIRONMENTAL SCIENCES↗

Satellite Observations of the Tropical Terrestrial Carbon Balance and Interactions With the Water Cycle During the 21st Century

A constellation of satellites are now in orbit providing information about terrestrial carbon and water storage and fluxes. These combined observations show that the tropical biosphere has changed significantly in the last two decades from the combined effects of climate variability and land use. Large areas of forest have been cleared in both wet and dry forests, increasing the source of carbon to the atmosphere. Concomitantly, tropical fire emissions have declined, at least until 2016, from changes in land-use practices and rainfall, increasing the net carbon sink. Measurements of carbon stocks and fluxes from disturbance and recovery and of vegetation photosynthesis show significant regional variability of net biosphere exchange (NBE) and gross primary productivity (GPP) across the tropics and are tied to seasonal and interannual changes in water fluxes and storage. Comparison of satellite based estimates of evapotranspiration (ET), photosynthesis, and the deuterium content of water vapor with patterns of total water storage and rainfall demonstrate the presence of vegetation-atmosphere interactions and feedback mechanisms across tropical forests. However, these observations of stocks, fluxes and inferred interactions between them do not point unambiguously to either positive or negative feedbacks in carbon and water exchanges. These ambiguities highlight the need for assimilation of these new measurements with Earth System models for a consistent assessment of process interactions, along with focused field campaigns that integrate ground, aircraft and satellite measurements, to quantify the controlling carbon and water processes and their feedback mechanisms.

58 GEOSCIENCES↗

The Response of Precipitation Extremes to the Twentieth- and Twenty-First-Century Global Temperature Change in a Comprehensive Suite of CESM1 Large Ensemble Simulation: Revisiting the Role of Forcing Agents Vs. the Role of Forcing Magnitudes

The response of precipitation extremes (PEs) to global warming is found to be nonlinear in Community Earth System Model version 1 (CESM1) and other global climate models (Pendergrass et al., 2019), which led to the concern that it is not accurate to approximate the response of PE to a single forcing as the difference between simulations with all forcing agents and those that exclude one specific forcing. This calls into question previous model-based results that the sensitivity of PE with warming due to aerosol forcing is significantly larger than that due to greenhouse gases (GHGs). We reevaluate the PE sensitivity to GHGs and aerosols using available CESM1 ensemble simulations. We show that although the PE response to warming is nonlinear in CESM1, especially for the high warming projected in the twenty-first-century, PE sensitivity to aerosols is still significantly stronger than that due to GHGs when the comparison is made within similar warming regimes to avoid the bias induced by the nonlinear behavior. But the difference is smaller than previously estimated. We also conclude that the additivity assumption is largely valid to isolate the PE response due to aerosol forcing from the paired simulations including the “all forcing” experiment when the warming regime is small (e.g., 1°C–2°C in the next few decades when aerosol forcing is projected to decline and becomes a major source of uncertainty for model projection).

58 GEOSCIENCES↗

Characterizing Drought Behavior in the Colorado River Basin Using Unsupervised Machine Learning

Drought is a pressing issue for the Colorado River Basin (CRB) due to the social and economic value of water resources in the region and the significant uncertainty of future drought under climate change. Here, we use climate simulations from various Earth System Models (ESMs) to force the Variable Infiltration Capacity hydrologic model and project multiple drought indicators for the sub-watersheds within the CRB. We apply an unsupervised machine learning (ML) based on Non-Negative Matrix Factorization using K-means clustering (NMFk) to synthesize the simulated historical, future, and change in drought indicators. The unsupervised ML approach can identify sub-watersheds where key changes to drought indicator behavior occur, including shifts in snowpack, snowmelt timing, precipitation, and evapotranspiration. While changes in future precipitation vary across ESMs, the results indicate that the Upper CRB will experience increasing evaporative demand and surface-water scarcity, with some locations experiencing a shift from a radiation-limited to a water-limited evaporation regime in the summer. Large shifts in peak runoff are observed in snowmelt-dominant sub-watersheds, with complete disappearance of the snowmelt signal for some sub-watersheds. The work demonstrates the utility of the NMFk algorithm to efficiently identify behavioral changes of drought indicators across space and time and to quickly analyze and interpret hydro climate model results.

54 ENVIRONMENTAL SCIENCES↗

Tree Canopies Reflect Mycorrhizal Composition

Mycorrhizae alter global patterns of CO 2 fertilization, carbon storage, and elemental cycling, yet knowledge of their global distributions is currently limited by the availability of forest inventory data. Here, we show that maps of tree-mycorrhizal associations (hereafter “mycorrhizal maps”) can be improved by the novel technology of imaging spectroscopy because mycorrhizal signatures propagate up from plant roots to impact forest canopy chemistry. We analyzed measurements from 143 airborne imaging spectroscopy surveys over 112,975 individual trees collected across 13 years. Results show remarkable accuracy in capturing ground truth observations of mycorrhizal associations from canopy signals across disparate landscapes (R 2 = 0.92, p < 0.01). Additionally, upcoming imaging spectroscopy satellite missions can reveal new insights into landscape-scale variations in water, nitrogen, phosphorus, carotenoid/anthocyanin, and cellulose/lignin composition. Applied globally, this approach could improve the spatial precision of mycorrhizal distributions by a factor of roughly 104 and facilitate the incorporation of dynamic shifts in forest composition into Earth system models.

58 GEOSCIENCES↗

Contribution of the Two Types of Ekman Pumping Induced Eddy Heat Flux to the Total Vertical Eddy Heat Flux

Based on eddy-rich (0.1° horizontal resolution for ocean) Community Earth System Model simulations, we compute the Ekman pumping induced vertical eddy heat flux using two formulations: the classical and the Stern-Ekman pumping. The two mechanisms yield similar patterns in the eddy-induced vertical heat flux, showing strong upward heat flux in ocean frontal regions such as the western boundary current regions. Using the Gulf Stream Extension region as an example, we estimate (1) the Stern-Ekman pumping induced eddy heat flux is 23% and 12% of the total vertical eddy heat flux at 50 m depth during summer and winter, respectively; (2) the classical Ekman pumping accounts for 82% and 88% of the Stern-Ekman pumping induced eddy heat flux at 50 m during summer and winter, respectively. The second finding indicates eddy current feedback on wind stress is the primary cause for the upward eddy heat flux generated by Ekman pumping.

58 GEOSCIENCES↗

Coupled Climate Responses to Recent Australian Wildfire and COVID-19 Emissions Anomalies Estimated in CESM2

In this work, multiple 50-member ensemble simulations with the Community Earth System Model version 2 are performed to estimate the coupled climate responses to the 2019–2020 Australian wildfires and COVID-19 pandemic policies. The climate response to the pandemic is found to be weak generally, with global-mean net top-of-atmosphere radiative anomalies of +0.23 ± 0.14 W m –2 driving a gradual global warming of 0.05 ± 0.04 K by the end of 2022. While regional anomalies are detectable in aerosol burdens and clear-sky radiation, few significant anomalies exist in other fields due to internal variability. In contrast, the simulated response to Australian wildfires is a strong and rapid cooling, peaking globally at –0.95 ± 0.15 W m –2 in late 2019 with a global cooling of 0.06 ± 0.04 K by mid-2020. Transport of fire aerosols throughout the Southern Hemisphere increases albedo and drives a strong interhemispheric radiative contrast, with simulated responses that are consistent generally with those to a Southern Hemisphere volcanic eruption.

54 ENVIRONMENTAL SCIENCES↗

Dissecting Anvil Cloud Response to Sea Surface Warming

Abstract We derive an anvil cloud diagnostic from the continuity equation of cloud ice and apply it to the output of convection‐permitting Energy Exascale Earth System Model (E3SM) simulations run in radiative‐convective equilibrium mode. This diagnostic shows that anvil cloud fraction can be reliably diagnosed as a product of cloud detrainment and lifetime. Detrainment is found to be approximated well by a product of clear sky convergence and cloud ice mixing ratio, while cloud lifetime is dominated by sedimentation. Taken together, this diagnostic expresses anvil cloud fraction as a function of five physically measurable quantities. Of these, clear‐sky convergence changes drive the anvil cloud reduction with warming while an increase in cloud ice mixing ratio buffers the decrease. Accordingly, this study provides a theoretical foundation upon which the Stability‐Iris hypothesis can be tested.

54 ENVIRONMENTAL SCIENCES↗

Linear Response Function Reveals the Most Effective Remote Forcing in Causing September Arctic Sea Ice Melting in CESM

We apply the linear response function method to investigate the most excitable mode of the September Arctic sea ice and its corresponding remote oceanic forcing in climate models. The method is useful in identifying the fundamental, dynamical and causal relationship between the remote forcing and sea ice response. We find that the most excitable mode of the September Arctic sea ice preferentially takes place over the Pacic side of the Arctic and its remote forcing corresponds to a dipole pattern of precipitation anomaly in the tropics with an increase of precipitation over the western and central tropical and subtropical Pacic ocean while a decrease over the Maritime Continent. The tropical precipitation anomaly likely drives a Rossby wave train propagating toward higher latitudes and leads to a ridge anomaly over the Pacic side of the Arctic, resulting in poleward atmospheric heat transport, enhanced downward longwave radiation and thus melting of the sea ice. In addition, a good agreement is found with the leading tropical Arctic teleconnection mode detected in a pre-industrial control simulation, which supports the usefulness and robustness of the linear response function method and the importance of tropical forcing in shaping the sea ice variability.

Wu, Yutian↗

Vegetation Response to Rising CO 2 Amplifies Contrasts in Water Resources Between Global Wet and Dry Land Areas

Rising atmospheric CO 2 impacts on vegetation physiological processes can alter land feedbacks on precipitation and water resources, but understanding of regional differences in these changes is uncertain. We investigate the impact of rising CO 2 on land water resources for different wetness levels using four Earth system models. We find an overall tendency of runoff to increase across all wetness levels. However, runoff increases in wet regions are much larger than those in dry regions, especially in wet seasons. This substantial amplification of contrasts between wet and dry regions increases at 3% per 100 ppm increase in CO 2 relative to the historical period, reaching 18% for a quadrupling of CO 2 , quantified by a new wetting contrast index (WCI). Physiological effects suppress evapotranspiration more in wet than dry regions, which has a larger contribution than radiative forcing to the amplification of runoff contrast, reshaping the spatial distribution of future land water resources.

58 GEOSCIENCES↗

The Oceanic Barrier Layer in the Eastern Indian Ocean as a Predictor for Rainfall Over Indonesia and Australia

Barrier layers in the tropics trap heat in a shallow and stable near-surface layer and limit entrainment of cooler water from below. Both processes act to increase sea surface temperature and enhance atmospheric convection. The high resolution fully coupled pre-industrial Energy Exascale Earth System Model version 0 (E3SMv0) is used to investigate the relationship between barrier layers in the eastern Indian Ocean during the wet season with local atmospheric convection and remote rainfall. A partial least squares regression reveals a significant relationship between Australasian rainfall and the barrier layer thickness (BLT) west of Sumatra, occurring one month earlier. Furthermore, the largest positive regression coefficients are over northern Australia. The region west of Sumatra is strategically located where the East-Asian monsoon moisture flows toward northern Australia. Thickening of the west Sumatra BLT intensifies evaporation and local convection and amplifies the moisture transported to Australia acting to increase the terrestrial rainfall.

58 GEOSCIENCES↗

Does Regional Hydroclimate Change Scale Linearly With Global Warming?

Many aspects of climate change scale linearly with global warming. However, nonlinear changes are possible, especially in the context of hydroclimate, and under emissions scenarios with stabilized global temperature, as aspired to by current climate targets. In CMIP5 and 6, a progressively larger land area shows nonlinear changes as a function of global warming when considering precipitation, evaporation, and soil moisture, with the latter showing nonlinearity over ~50% of global land. Here, using ensemble simulations with the Community Earth System Model 1, in which individual forcing factors are held constant, we illustrate how nonadditive responses to anthropogenic greenhouse gases and industrial and fire-related aerosols can yield complex soil moisture changes in certain regions. This complexity contributes to uncertainty in regional soil moisture projections and suggests that the timing of, as well as model response uncertainty to, future aerosol reductions will have significant impacts on regional hydroclimate change as global temperatures stabilize.

54 ENVIRONMENTAL SCIENCES↗

The Role of Coupled Feedbacks in the Decadal Variability of the Southern Hemisphere Eddy-Driven Jet

The Southern Hemisphere summertime eddy-driven jet and storm tracks have shifted poleward over the recent few decades. In previous studies, explanations have mainly stressed the influence of external forcing in driving this trend. Here we examine the role of internal tropical SST variability in controlling the austral summer jet’s poleward migration, with a focus on interdecadal time scales. The role of external forcing and internal variability are isolated by using a hierarchy of Community Earth System Model version 1 (CESM1) simulations, including the pre-industrial control, large ensemble, and pacemaker runs. Model simulations suggest that in the early twenty-first century, both external forcing and internal tropical Pacific SST variability are important in driving a positive southern annular mode (SAM) phase and a poleward migration of the eddy-driven jet. Additionally, Tropical Pacific SST variability, associated with the negative phase of the interdecadal Pacific oscillation (IPO), acts to shift the jet poleward over the southern Indian and southwestern Pacific Oceans and intensify the jet in the southeastern Pacific basin, while external forcing drives a significant poleward jet shift in the South Atlantic basin. In response to both external forcing and decadal Pacific SST variability, the transient eddy momentum flux convergence belt in the middle latitudes experiences a poleward migration due to the enhanced meridional temperature gradient, leading to a zonally symmetric southward migration of the eddy-driven jet. This mechanism distinguishes the influence of the IPO on the midlatitude circulation from the dynamical impact of ENSO, with the latter mainly promoting the subtropical wave-breaking critical latitude poleward and pushing the midlatitude jet to higher latitudes.

54 ENVIRONMENTAL SCIENCES↗

Quantifying Carbon Cycle Extremes and Attributing Their Causes Under Climate and Land Use and Land Cover Change From 1850 to 2300

The increasing atmospheric carbon dioxide (CO 2 ) mole fraction affects global climate through radiative (trapping longwave radiation) and physiological effects (reduction of plant transpiration). We use the simulations of the Community Earth System Model (CESM1-BGC) forced with Representative Concentration Pathway 8.5 to investigate climate-vegetation feedbacks from 1850 to the year 2300. Human-induced land use and land cover change (LULCC), through biogeochemical and biogeophysical processes, alter the climate and modify photosynthetic activity. The changing characteristics of extreme anomalies in photosynthesis, referred to as carbon cycle extremes, increase the uncertainty of terrestrial ecosystems to act as a net carbon sink. However, the role of LULCC in altering carbon cycle extremes under business-as-usual (continuously rising) CO 2 emissions is unknown. Here we show that LULCC magnifies the intensity, frequency, and extent of carbon cycle extremes, resulting in a net reduction in expected photosynthetic activity in the future. We found that large temporally contiguous negative carbon cycle extremes are due to a persistent decrease in soil moisture, which is triggered by declines in precipitation. With LULCC and global warming, vegetation exhibits increased vulnerability to hot and dry environmental conditions, increasing the frequency of fire events and resulting in considerable losses in photosynthetic activity. While most regions show strengthening of negative carbon cycle extremes, a few locations show a weakening effect driven by declining vegetation cover or benign climate conditions for photosynthesis. Increasing hot, dry, and fire-driven carbon cycle extremes are essential for improving carbon cycle modeling and estimation of ecosystem responses to LULCC and rising CO 2 mole fractions. Moreover, large aberrations in vegetation productivity represent potential and growing threats to human lives, wildlife, and food security.

54 ENVIRONMENTAL SCIENCES↗

Maximum Supersaturation in the Marine Boundary Layer Clouds Over the North Atlantic

Abstract The maximum supersaturation ( S x ) in clouds is a key parameter affecting the cloud's microphysical and radiative properties. We investigate the S x of the marine boundary layer clouds by combining airborne and surface observations in the Eastern North Atlantic. The cloud droplet number concentration ( N c ) in the least diluted cloud cores agrees well with the number concentration of particles larger than the Hoppel Minimum (HM) ( N >HM ) below clouds, indicating that the HM represents the average size threshold above which particles are activated to form cloud droplets. The S x values derived from surface observations vary from 0.10% to 0.50% from June 2017 to June 2018, with a clear seasonal variation exhibiting higher values during winter. Most of the S x variance (∼60%) can be explained by the cloud condensation nuclei (CCN) concentration and updraft velocity ( w ), with the CCN concentration playing a more important role than w in explaining the variation of S x . The influence of CCN concentration on S x leads to a buffered response of N c to aerosol perturbations. The response of N c to low aerosol concentration during winter is further buffered by the high w . The global Community Earth System Model (CESM) simulated S x values in the Azores have a positive bias compared to measured S x , likely due to overestimated w and underestimated CCN concentration. The CESM simulated S x exhibits higher values further north over the North Atlantic Ocean, which is attributed to stronger w . The suppression of S x by aerosol is also evident in regions with high CCN concentrations.

54 ENVIRONMENTAL SCIENCES↗