Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “global models”

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 523 records · Page 29

Contextual Active Online Model Selection with Expert Advice

How can we collect the most useful labels to learn a model selection policy, when presented with arbitrary heterogeneous data streams? In this paper, we formulate this task as a contextual active model selection problem, where at each round the learner receives an unlabeled data point along with a context. The goal is to output the best model for any given context without obtaining an excessive amount of labels. In particular, we focus on the task of selecting pre-trained classifiers, and propose a contextual active model selection algorithm (CAMS), which relies on a novel uncertainty sampling query criterion defined on a given policy class for adaptive model selection. In comparison to prior art, our algorithm does not assume a globally optimal model. We provide rigorous theoretical analysis for the regret and query complexity under both adversarial and stochastic settings. Our experiments on several benchmark classification datasets demonstrate the algorithm’s effectiveness in terms of both regret and query complexity. Notably, to achieve the same accuracy, CAMS incurs less than 10% of the label cost when compared to the best online model selection baselines on CIFAR10.

Liu, Xuefeng↗

Mesoscale Convective Systems in a Superparameterized E3SM Simulation at High Resolution

Accurately representing mesoscale convective systems (MCSs) is crucial to simulating the energy and water cycles in global climate models. Using a novel MCS identification and tracking algorithm applied to observations and model simulations, we evaluate how well the Energy Exascale Earth System Model simulates MCSs over the central US. Simulations performed by E3SM at 25 km grid spacing with (SP-E3SM) and without (E3SM) superparameterization are compared and evaluated against observations using multiple metrics highlighting important MCS characteristics. Compared to E3SM, SP-E3SM better simulates MCS number and MCS precipitation amount, diurnal cycle, propagation, and the probability distribution of precipitation rate in both spring and summer. The improvement from SP is partly contributed by improvement in simulating the large-scale environments, featuring enhanced atmospheric low-level moisture and larger moisture transport to the central US relative to E3SM. However, SP-E3SM still underestimates MCS precipitation amount, particularly in summer. This underestimation is closely related to the negative bias in MCS precipitation intensity, although the drier environments simulated during summer also contributes to underestimation of MCS frequency. Without SP, the larger bias in MCS precipitation amount is closely related to the negative bias in MCS frequency, while the substantial dry bias in the large-scale environments also contributes to underestimation of MCS intensity. Our results suggest that SP improves MCS simulation by improving modeling of the large-scale environments and convection initiation, which are both major limiting factors in E3SM even at 25 km grid spacing where deep convection is represented by a cumulus parameterization.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of CMIP6 GCMs Over the CONUS for Downscaling Studies

Despite the necessity of Global Climate Models (GCMs) sub-selection in downscaling studies, an objective approach for their selection is currently lacking. Building on the previously established concepts in GCMs evaluation frameworks, we develop a weighted averaging technique to remove the redundancy in the evaluation criteria and rank 37 GCMs from the sixth phase of the Coupled Models Intercomparison Project over the contiguous United States. GCMs are rated based on their average performance across 66 evaluation measures in the historical period (1981–2014) after each metric is weighted between zero and one, depending on its uniqueness. The robustness of the outcome is tested by repeating the process with the empirical orthogonal function analysis in which each GCM is ranked based on its sum of distances from the reference in the principal component space. The two methodologies work in contrasting ways to remove the metrics redundancy but eventually develop similar GCMs rankings. A disparity in GCMs' behavior related to their sensitivity to the size of the evaluation suite is observed, highlighting the need for comprehensive multi-variable GCMs evaluation at varying timescales for determining their skillfulness over a region. The sub-selection goal is to use a representative set of skillful models over the region of interest without substantial overlap in their future climate responses and modeling errors in representing historical climate. Additional analyses of GCMs' independence and spread in their future projections provide the necessary information to objectively select GCMs while keeping all aspects of necessity in view.

54 ENVIRONMENTAL SCIENCES↗

Future projections of storm surge in Hurricane Katrina and sensitivity to meteorological forcing resolution

In this study, we investigated whether and how the storm surge induced by Hurricane Katrina could change if it occurs in a future warmer climate, and the sensitivity of the changes to atmospheric forcing resolution. Climate model simulations of Hurricane Katrina at 27 km, 4.5 km, and 3 km resolutions were used to drive storm surge simulations in historical and future climates using the ADvanced CIRCulation (ADCIRC) model. We found that peak surge height increased significantly in the future with all forcing resolutions. However, the future projection is 22% greater in the 3 km forcing, typical of regional climate models, compared to the 27 km forcing, typical of state-of-the-art global climate models. Additionally, the spatial extent of the future change is highly sensitive to forcing resolution, extending most broadly under the 27 km forcing. Furthermore, we found that storm surge duration decreases in the future with all forcing resolutions due to increasing TC translation speed and decreasing ocean lifetime. However, the future change in the surge duration is sensitive to the forcing resolution, decreasing by 31% in the 27 km forcing and 6% in the 3 km forcing.

54 ENVIRONMENTAL SCIENCES↗

A global DSCOVR/EPIC-based hourly/daily shortwave radiation/PAR dataset

Accurate characterization of the global downward shortwave (SW) and photosynthetically active radiation (PAR) is fundamental for Earth system modeling and global change research. Combined with a machine-learning method, we used the Earth Polychromatic Imaging Camera (EPIC) data onboard the Deep Space Climate Observatory (DSCOVR) to generate the global SW/PAR dataset at both hourly and daily temporal scales and 0.1degree×0.1degree spatial resolution. Covering a period from June, 2015 to June (about 4 years), 2019, the dataset consists of total SW/PAR and their direct and diffuse components. Our dataset will be valuable for understanding the diurnal cycles and seasonal patterns of ecosystem processes, as well as the solar energy utilizations

14 SOLAR ENERGY↗

Resolution dependence of the turbulent atmospheric boundary layer in global storm-resolving climate simulations

The current generation of state-of-the-art global climate models are being run at increasingly higher horizontal resolutions, with the goal of resolving organised deep convection explicitly. How a kilometre-scale resolution impacts the representation of the atmospheric boundary layer is, however, not well known. Using statistical analysis on global fields as well as high-frequency data at selected locations, produced with the Integrated Forecasting System (IFS) model for the Next Generation Earth-system Models (nextGEMS) project, we investigate the horizontal resolution dependence of some boundary-layer processes. We find that a change in resolution from 9 to 2.8 km causes no substantial changes to boundary-layer properties and processes at most of the locations studied, although some global changes are detected that indicate circulation changes. Small changes to the boundary-layer depth and structure are found in the Tropics. The short simulation length and lack of data for optimal boundary-layer analysis limits the conclusions, especially in relation to the connection between the boundary layer and the atmosphere general circulation.

54 ENVIRONMENTAL SCIENCES↗

Hydrological projections over the Upper Indus Basin at 1.5 °C and 2.0 °C temperature increase

We analyse an ensemble of statistically downscaled Global Climate Models (GCMs) to investigate future water availability in the Upper Indus Basin (UIB) of Pakistan for the time horizons when the global and/or regional warming levels cross Paris Agreement (PA) targets. The GCMs data is obtained from the 5th Phase of Coupled Model Inter-Comparison Project under two Representative Concentration Pathways (RCP4.5 and RCP8.5). Based on the five best performing GCMs, we note that global 1.5 °C and 2.0 °C warming thresholds are projected in 2026 and 2047 under RCP4.5 and 2022 and 3036 under RCP8.5 respectively while these thresholds are reached much earlier over Pakistan i.e. 2016 and 2030 under RCP4.5 and 2012 and 2025 under RCP8.5 respectively. Interestingly, the GCMs with the earliest emergence at the global scale are not necessarily the ones with the earliest emergence over Pakistan, highlighting spatial non-linearity in GCMs response. The emergence of 2.0 °C warming at global scale across 5 GCMs ranges from 2031 (CCSM4) to 2049 (NorESM) under RCP8.5. Precipitation generally exhibits a progressive increasing trend with stronger changes at higher warming or radiative forcing levels. Hydrological simulations representing the historical, 1.5 °C and 2.0 °C global and region warming time horizons indicate a robust but seasonally varying increase in the inflows. The highest inflows in the baseline and future are witnessed in July. However, the highest future increase in inflows is projected in October under RCP4.5 (37.99% and 65.11% at 1.5 °C and 2.0 °C) and in April under RCP8.5 (37% and 62.05% at 1.5 °C and 2.0 °C). These hydrological changes are driven by increases in the snow and glacial melt contribution, which are more pronounced at 2.0 °C warming level. These findings should help for effective water management in Pakistan over the coming decades.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of CMIP6 models in simulating the statistics of extreme precipitation over Eastern Africa

We report the Eastern Africa region experiences frequent extreme precipitation events that can cause destruction of property and environment, and loss of lives. Thus, there is a need to understand how these events may change in the future and how well the global climate models that are used to make projections can simulate precipitation extremes in this region before they can be used in downscaling or flood and drought impact assessment studies. In this work, we evaluated the ability of sixteen Coupled Model Intercomparison Project Phase 6 (CMIP6) models to simulate present-day precipitation extremes over the Eastern Africa region during the two rainy seasons (March–May and September–November). We used nine extreme precipitation indices (including seven (one) indices of wet (dry) extremes) defined by the Expert Team on Climate Change Detection and Indices. The CMIP6 models were evaluated against two gridded observation datasets: Global Precipitation Climatology Project One-Degree Daily Dataset and Tropical Rainfall Measuring Mission Multi-satellite Precipitation Analysis 3B42. Three model performance metrics (percentage bias, normalized root-mean-square error, and pattern correlation coefficient) were employed to further assess the strengths and weakness of the models. Our results show that the multi-model ensemble mean generally provides a better representation of observed precipitation and related extremes compared to individual models when considering all metrics and seasons. Several consistent biases are evident across CMIP6 models, which tend to overestimate the total-wet day precipitation and consecutive wet days, and underestimate very wet days and maximum 5-day precipitation in both seasons. Furthermore, no single model consistently performs best, model performance varies with the season and index under consideration and is generally independent of horizontal resolution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Assessment of CMIP6 Performance and Projected Temperature and Precipitation Changes Over South America

We evaluate the performance of a large ensemble of Global Climate Models (GCMs) from the Coupled Model Intercomparison Project Phase 6 (CMIP6) over South America for a recent past reference period and examine their projections of twenty-first century precipitation and temperature changes. The future changes are computed for two time slices (2040–2059 and 2080–2099) relative to the reference period (1995–2014) under four Shared Socioeconomic Pathways (SSPs, SSP1–2.6, SSP2–4.5, SSP3–7.0 and SSP5–8.5). The CMIP6 GCMs successfully capture the main climate characteristics across South America. However, they exhibit varying skill in the spatiotemporal distribution of precipitation and temperature at the sub-regional scale, particularly over high latitudes and altitudes. Future precipitation exhibits a decrease over the east of the northern Andes in tropical South America and the southern Andes in Chile and Amazonia, and an increase over southeastern South America and the northern Andes—a result generally consistent with earlier CMIP (3 and 5) projections. However, most of these changes remain within the range of variability of the reference period. In contrast, temperature increases are robust in terms of magnitude even under the SSP1–2.6. Future changes mostly progress monotonically from the weakest to the strongest forcing scenario, and from the mid-century to late-century projection period. There is an increase in the seasonality of the intra-annual precipitation distribution, as the wetter part of the year contributes relatively more to the annual total. Furthermore, an increasingly heavy-tailed precipitation distribution and a rightward shifted temperature distribution provide strong indications of a more intense hydrological cycle as greenhouse gas emissions increase. The relative distance of an individual GCM from the ensemble mean does not substantially vary across different scenarios. We found no clear systematic linkage between model spread about the mean in the reference period and the magnitude of simulated sub-regional climate change in the future period. Overall, these results could be useful for regional climate change impact assessments across South America.

54 ENVIRONMENTAL SCIENCES↗

A nonhydrostatic formulation for MPAS-Ocean

The Model for Prediction Across Scales-Ocean (MPAS-Ocean) is an open-source, global ocean model and is one component of a family of climate models within the MPAS framework, including atmosphere, sea-ice, and land-ice models. Here, in this work, a new formulation for the ocean model is presented that solves the nonhydrostatic, incompressible Boussinesq equations on an unstructured, staggered, z-level grid. The introduction of this nonhydrostatic capability is necessary for the resolution of internal wave dynamics and large eddy simulations. Compared to the standard, hydrostatic formulation, a nonhydrostatic pressure solver and a vertical momentum equation are added, where the PETSc (Portable Extensible Toolkit for Scientific Computation) library is used for the inversion of a large sparse system for the nonhydrostatic pressure. Numerical results on a stratified seiche, internal solitary wave, overflow and lock-exchange test cases are presented, and the parallel efficiency of the code is evaluated using up to 1024 processors.

3D Poisson equation↗

A nonhydrostatic formulation for MPAS-Ocean

The Model for Prediction Across Scales-Ocean (MPAS-Ocean) is an open-source, global ocean model and is one component of a family of climate models within the MPAS framework, including atmosphere, sea-ice, and land-ice models. Here, in this work, a new formulation for the ocean model is presented that solves the nonhydrostatic, incompressible Boussinesq equations on an unstructured, staggered, z-level grid. The introduction of this nonhydrostatic capability is necessary for the resolution of internal wave dynamics and large eddy simulations. Compared to the standard, hydrostatic formulation, a nonhydrostatic pressure solver and a vertical momentum equation are added, where the PETSc (Portable Extensible Toolkit for Scientific Computation) library is used for the inversion of a large sparse system for the nonhydrostatic pressure. Numerical results on a stratified seiche, internal solitary wave, overflow and lock-exchange test cases are presented, and the parallel efficiency of the code is evaluated using up to 1024 processors.

3D Poisson equation↗

Winter Warming in North America Induced by Urbanization in China

Rapid urbanization in China in the past several decades has been shown to affect surface air temperature (SAT) locally, but its global impact is unclear. Using a global climate model, we have investigated the potential impact of urbanization in China on SAT at the global scale. The expansion of urban areas in China from 1985 to 2017 leads to winter surface air warming by as much as 0.8°C in North America (NA), largely caused by warm temperature advection due to global-scale atmospheric circulation changes. The East Asian jet stream moves northward in eastern and northeastern China where the expansion of urban areas has increased SAT and sensible heat flux, forcing change of the atmospheric stationary waves and so the structure of the jet streams.

54 ENVIRONMENTAL SCIENCES↗

DESI DR2 Galaxy luminosity functions

We present galaxy luminosity functions (LFs) for the Dark Energy Spectroscopic Instrument (DESI) DR2 Bright Galaxy Survey (BGS) in the g, r, z, and $w1$ bands over $0.002\lt z\lt 0.6$. Our analysis uses updated k-corrections and evolutionary corrections, including new polynomial kcorrection fits derived from BGS Year 1 data that supersede earlier GAMA-based prescriptions. Exploiting the statistical power of DESI, we measure LFs to very faint magnitudes, reaching $^{0.1}M_r-5\log h\sim -10$. Independent measurements from the North and South survey regions show excellent agreement around the LF knee, but the very small statistical uncertainties reveal that simple analytic forms fail to capture the full LF shape. The bright end departs from a pure exponential decline, while the faint end exhibits complex, non-powerlaw behaviour, including a pronounced upturn at $^{0.1}M_r-5\log h\gtrsim -15$, which is stronger for red galaxies than for blue. We show that our LFs are largely complete for galaxies with surface brightness $\mu _{50}\lt 25$, and that an apparent steepening fainter than $-13$ is driven primarily by local overdensity and fragmentation of large galaxies. A systematic North–South offset at the brightest magnitudes is traced to red galaxies and may reflect shallower North photometry underestimating extended earlytype profiles, although this remains inconclusive. We therefore also provide LFs based on model Petrosian magnitudes. Redshift splitting reveals small but significant residuals, indicating limitations of a simple global evolutionary model. Using the redshift limits of J. Loveday et al. (2012), we find excellent agreement with GAMA, with substantially reduced statistical errors. These measurements provide a precise reference for studies of environmental and population-dependent LFs and for testing galaxy formation models.

79 ASTRONOMY AND ASTROPHYSICS↗

Historical evaluation and future projections of compound heatwave and drought extremes over the conterminous United States in CMIP6 *

Abstract Independently, both droughts and heatwaves can induce severe impacts on human and natural systems. However, when these two climate extremes occur concurrently in a given region, their compound impacts are often more pronounced. With the improvement in both the spatiotemporal resolution and representation of complex climate processes in the global climate models (GCMs), they are increasingly used to study future changes in these extremes and associated regional impacts. However, GCM selection for such impact assessments is generally based on historical performance and/or future mean changes, without considering individual or compound extremes. In contrast, this study evaluates historical performance and projected changes in heatwaves, droughts, and compound heatwave-droughts using an ensemble of GCMs from the latest Phase 6 of Coupled Models Intercomparison Project at a regional scale across the conterminous United States. Additionally, we explore the inter-model differences in the projected changes that are associated with various characteristics of extremes and the choice of drought indices. Our analysis reveals considerable variation among the GCMs, as well as substantial differences in the projected changes based on the choice of drought indices and region of interest. For example, the projected increases in both the frequency and intensity of drought and associated compound extreme days, based on the standardized precipitation evapotranspiration index far exceed those derived from the standard precipitation index. Further, the largest changes in the frequency of compound extremes are projected over the Southwest, South Central, and parts of the Southeast while the smallest changes are projected over the Northeast. Overall, this study provides important insights for the interpretation and selection of GCMs for future assessment studies that are crucial for the development of regional adaptation strategies in the face of climate change.

54 ENVIRONMENTAL SCIENCES↗

Quantifying the Known Unknown: Including Marine Sources of Greenhouse Gases in Climate Modeling

Researchers have recently estimated that Arctic submarine permafrost currently traps 60 billion tons of methane and contains 560 billion tons of organic carbon in seafloor sediments and soil, a giant pool of carbon with potentially large feedbacks on the climate system. Unlike terrestrial permafrost, the submarine permafrost system has remained a “known unknown” because of the difficulty in acquiring samples and measurements. Consequently, this potentially large carbon stock never yet considered in global climate models or policy discussions, represents a real wildcard in our understanding of Earth’s climate. This report summarizes our group’s effort at developing a numerical modeling framework designed to produce a first-of-its-kind estimate of Arctic methane gas releases from the marine sediments to the water column, and potentially to the atmosphere, where positive climate feedback may occur. Newly developed modeling capability supported by the Laboratory Directed Research and Development (LDRD) program at Sandia National Laboratories now gives us the ability to probabilistically map gas distribution and quantity in the seabed by using a hybrid approach of geospatial machine learning, and predictive numerical thermodynamic ensemble modeling. The novelty in this approach is its ability to produce maps of useful data in regions that are only sparsely sampled, a common challenge in the Arctic, and a major obstacle to progress in the past. By applying this model to the circum-Arctic continental shelves and integrating the flux of free gas from in situ methanogenesis and dissociating gas hydrates from the sediment column under climate forcing, we can provide the most reliable estimate of a spatially and temporally varying source term for greenhouse gas flux that can be used by global oceanographic circulation and Earth system models (such as DOE’s E3SM). The result will allow us to finally tackle the wildcard of the submarine permafrost carbon system, and better inform us about the severity of future national security threats that sustained climate change poses.

54 ENVIRONMENTAL SCIENCES↗

Burning conditions and transportation pathways determine biomass-burning aerosol properties in the Ascension Island marine boundary layer

African biomass-burning aerosol (BBA) in the southeast Atlantic Ocean (SEA) marine boundary layer (MBL) is an important contributor to Earth’s radiation budget yet its representation remains poorly constrained in regional and global climate models. Data from the Layered Atlantic Smoke Interactions with Clouds (LASIC) field campaign on Ascension Island (-7.95° N, -14.36° E) detail how fire source regions (burning conditions and fuel type), transport pathways, and longer-term chemical processing affect the chemical, microphysical, and optical properties of the BBA in the remote MBL between June and September of 2017. Ten individual plume events characterize the seasonal evolution of BBA characteristics. Inefficient burning conditions, determined by the mass ratio of refractory black carbon to above-background carbon monoxide (rBC:ΔCO), enhance organic- and sulfate-rich aerosol concentrations in June–July. In contrast, the heart of the burning season exhibited higher rBC:ΔCO values indicative of efficient burning conditions, correlating with more rBC-enriched BBA. Toward the end of the burning season, a mix of burning conditions results in increased variation of the BBA properties. The BBA transit to Ascension Island was predominantly through slow-moving pathways in the MBL and lower free troposphere (FT), facilitating prolonged chemical transformations through heterogeneous and aqueous phase processes. Heterogeneous oxidation can persist for up to 10 days, resulting in a considerable decrease in organic aerosol (OA) mass. OA to rBC mass ratios (OA:rBC) in the MBL between 2 and 5 contrast to higher values of 5 to 15 observed in the nearby FT. Conversely, early-season aqueous-phase processes primarily contributed to aerosol oxidation and some aerosol production, but not appreciable aerosol removal. These two chemical processes yield more light-absorbing BBA in the MBL than in the FT and explain the notably low scattering albedo at 530 nm (SSA 530 ) values (< 0.80) at Ascension Island. This study establishes a robust correlation between SSA 530 and OA:rBC across both MBL and FT, underscoring the dependency of optical properties on chemical composition. These findings highlight how the interplay between chemical composition and atmospheric processing can be improved in global and regional climate models. Questions remain on the mixing of aerosols with different pathway histories, and on what accounts for the doubling of the mass absorption coefficient in the boundary layer.

54 ENVIRONMENTAL SCIENCES↗

A novel network-based approach to determining measurement representation error for model evaluation of aerosol microphysical properties

Atmospheric aerosol size and abundance influence radiative effects and climate change. To date, efforts to constrain global climate models’ radiative forcing with in situ aerosol observations have been hamstrung by uncertainty. One source of error, the regional “representation error,” arises when accurate but sparse single-point measurements of atmospheric aerosol distributions are compared with a model value, assuming that the single-point measurement is representative of the model domain. The Portable Optical Particle Spectrometer network in the Southern Great Plains (POPSnet-SGP) campaign has demonstrated that a network of nearly autonomous aerosol instruments operating at ambient temperature and relative humidity (with low measurement error) may be used to quantify measurement representation error and investigate the factors introducing heterogeneity in aerosol distributions across a rural, continental background region. Measurements were made using Portable Optical Particle Spectrometer (POPS) instruments at several sites for five months across the Department of Energy’s Aerosol Radiation Measurement Southern Great Plains (ARM-SGP) User Facility in the central USA. Measurement representation error decreased with longer averaging periods (20-40 % between 1 sec and 1 day), varied between sites by 10 – 20 % for aerosol concentration 140 – 2500 nm in diameter (N_140), and was higher for aerosols > 400 nm in diameter (N_400). Our measurements also show the influence of local meteorology on aerosol surface area (A_140) and size distributions: A_140 is positively correlated with wind speed and relative humidity, negatively correlated with precipitation, and lower given westerly winds. Based on this study, we conclude that the POPSnet approach provides considerably more insight into the spatial variability in the aerosol population that can be used to constrain climate models than would be available from similar networks of PM 2.5 monitors.

54 ENVIRONMENTAL SCIENCES↗

MERIT-Plus Dataset: Delineation of endorheic basins in 5 and 15 min upscaled river networks

Endorheic drainage basins, those inland basins not connected directly to ocean, are essential for hydrological modeling of global and regional water balances, land surface water storage, gravity anomalies, sea level rise, etc. Within many hydrological model frameworks, river basins are defined by digital river networks through their flow direction and connectivity datasets. Here we present an improvement to gridded flow direction data and its derivatives produced from upscaled global 5 and 15 arc minute MERIT networks. We explicitly label endorheic and exorheic drainage basins and alter the delineation of endorheic basins by merging small inland watersheds to the adjacent host basins. The resulting datasets have a significantly reduced number of endorheic basins while preserving the total land portion of those basins since most of the merged catchments were inside other larger endorheic areas. We developed and present here the endorheic basin delineation method. This method performs an analysis of the contributing river and basin geometry relative to the location of the flow end point (i.e. potential endorheic lake), proximity of the latter to the drainage basin boundary and the elevation difference between the basin's lowest point and potential spillover location at the basin boundary. The new digital river network was validated using the University of New Hampshire Water Balance Model by comparing the water balance of endorheic inland depressions with modeled accumulation of water in their inland lakes based on the observed historical climate drivers used by WBM.

Earth Systems↗