Probabilistic Modeling of Climate Change Impacts on Renewable Energy and Storage Requirements for NM's Energy Transition Act.
Abstract not provided.
SEARCH · Engineering Papers
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.
Abstract not provided.
Abstract not provided.
We use a nascent data-driven causal discovery method to find and compare causal relationships in observed data and climate model output. We consider ten different features in the Arctic climate collected from public databases on observational and Energy Exascale Earth System Model (E3SM) data. In identifying and analyzing the resulting causal networks, we make meaningful comparisons between observed and climate model interdependencies. This work demonstrates our ability to apply the PCMCI causal discovery algorithm to Arctic climate data, that there are noticeable similarities between observed and simulated Arctic climate dynamics, and that further work is needed to identify specific areas for improvement to better align models with natural observations.
Abstract Model evidence for the “pattern effect” assumes that global climate models (GCMs) faithfully simulate how clouds respond to varying sea surface temperature (SST) patterns and associated meteorological perturbations. We exploit time-invariant satellite-based estimates of the sensitivity of marine low clouds to meteorological perturbations to estimate how these clouds responded to time-varying SST patterns and meteorology between 1870 and 2014. GCMs and reanalyses provide estimates of the historical meteorological changes. Observations suggest that increasing estimated inversion strength (EIS) between 1980 and 2014 produced a negative low cloud feedback, opposite to the positive feedback expected from increasing CO 2 . This indicates that the processes responsible for marine cloud changes from 1980 to the near present are distinct from those associated with an increase in CO 2 . We also observationally constrain the difference between the historical near-global marine low cloud feedback, , and that arising from increasing CO 2 , . We find that this cloud feedback pattern effect depends strongly on time period and reanalysis dataset, and that varying changes in EIS and SST with warming explain much of its variability. Between 1980 and 2014, we estimate that (90% confidence) assuming meteorological changes from the Multiple Reanalysis Ensemble, implying a total pattern effect (that arising from all climate feedbacks) of 1.86 ± 0.45 W m −2 K −1 . This observational evidence corroborates previous quantitative estimates of the pattern effect, which heretofore relied largely upon GCM-based cloud changes. However, disparate historical meteorological changes across individual reanalyses contribute to considerable uncertainty in its magnitude.
This study explores the added value (AV) of a regional Earth system model (ESM) compared to an atmosphere-only regional climate model (RCM) in simulating West African monsoon (WAM) rainfall. The primary goals are to foster discussions on the suitability of coupled RCMs for WAM projections and deepen our understanding of ocean–atmosphere coupling’s influence on the WAM system. The study employs results from dynamical downscaling of the ERA-Interim reanalysis and Max Plank Institute ESM, low resolution (MPI-ESM-LR), by two RCMs, atmosphere only (REMO) and REMO coupled with Max Planck Institute Ocean Model (MPIOM) (ROM), at ~25-km horizontal resolution. Results show that in regions distant from coupling domain boundaries such as West Africa (WA), constraint conditions from ERA-Interim are more beneficial than coupling effects. REMO, reliant on oceanic sea surface temperatures (SSTs) from observations and influenced by ERA-Interim, is biased under coupling conditions, although coupling offers potential advantages in representing heat and mass fluxes. Contrastingly, as intended, coupling improves SSTs and monsoon fluxes’ relationships under ESM-forced conditions. In this latter case, the coupling features a dipole-like spatial structure of AV, improving precipitation over the Guinea Coast but degrading precipitation over half of the Sahel. Our extensive examination of physical processes and mechanisms underpinning the WAM system supports the plausibility of AV. Additionally, we found that the monsoonal dynamics over the ocean respond to convective activity, with the Sahara–Sahel surface temperature gradient serving as the maintenance mechanism. While further efforts are needed to enhance the coupled RCM, we advocate for its use in the context of WAM rainfall forecasts and projections.
Abstract Large ensembles of model simulations are frequently used to reduce the impact of internal variability when evaluating climate models and assessing climate change induced trends. However, the optimal number of ensemble members required to distinguish model biases and climate change signals from internal variability varies across models and metrics. Here we analyze the mean, variance and skewness of precipitation and sea surface temperature in the eastern equatorial Pacific region often used to describe the El Niño–Southern Oscillation (ENSO), obtained from large ensembles of Coupled model intercomparison project phase 6 climate simulations. Leveraging established statistical theory, we develop and assess equations to estimate, a priori, the ensemble size or simulation length required to limit sampling‐based uncertainties in ENSO statistics to within a desired tolerance. Our results confirm that the uncertainty of these statistics decreases with the square root of the time series length and/or ensemble size. Moreover, we demonstrate that uncertainties of these statistics are generally comparable when computed using either pre‐industrial control or historical runs. This suggests that pre‐industrial runs can sometimes be used to estimate the expected uncertainty of statistics computed from an existing historical member or ensemble, and the number of simulation years (run duration and/or ensemble size) required to adequately characterize the statistic. This advance allows us to use existing simulations (e.g., control runs that are performed during model development) to design ensembles that can sufficiently limit diagnostic uncertainties arising from simulated internal variability. These results may well be applicable to variables and regions beyond ENSO.
The land ice contribution to global mean sea level rise has not yet been predicted using ice sheet and glacier models for the latest set of socio-economic scenarios, nor using coordinated exploration of uncertainties arising from the various computer models involved. Two recent international projects generated a large suite of projections using multiple models, but primarily used previous-generation scenarios and climate models, and could not fully explore known uncertainties. Here we estimate probability distributions for these projections under the new scenarios using statistical emulation of the ice sheet and glacier models. In this work, we find that limiting global warming to 1.5 degrees Celsius would halve the land ice contribution to twenty-first-century sea level rise, relative to current emissions pledges. The median decreases from 25 to 13 centimetres sea level equivalent (SLE) by 2100, with glaciers responsible for half the sea level contribution. The projected Antarctic contribution does not show a clear response to the emissions scenario, owing to uncertainties in the competing processes of increasing ice loss and snowfall accumulation in a warming climate. However, under risk-averse (pessimistic) assumptions, Antarctic ice loss could be five times higher, increasing the median land ice contribution to 42 centimetres SLE under current policies and pledges, with the 95th percentile projection exceeding half a metre even under 1.5 degrees Celsius warming. This would severely limit the possibility of mitigating future coastal flooding. Given this large range (between 13 centimetres SLE using the main projections under 1.5 degrees Celsius warming and 42 centimetres SLE using risk-averse projections under current pledges), adaptation planning for twenty-first-century sea level rise must account for a factor-of-three uncertainty in the land ice contribution until climate policies and the Antarctic response are further constrained.
Anthropogenic climate change is unfolding rapidly, yet its regional manifestation can be obscured by internal variability. A primary goal of climate science is to identify the externally forced climate response from among the noise of internal variability. Separating the forced response from internal variability can be addressed in climate models by using a large ensemble to average over different possible realizations of internal variability. However, with only one realization of the real world, it is a major challenge to isolate the forced response directly in observations. In the Forced Component Estimation Statistical Method Intercomparison Project (ForceSMIP), contributors used existing and newly developed statistical and machine learning methods to estimate the forced response over 1950–2022 within individual realizations of the climate system. Participants used neural networks, linear inverse models, fingerprinting methods, and low-frequency component analysis, among other approaches. These methods were trained using large ensembles from multiple climate models and then applied to observations. Here, we evaluate method performance within large ensembles and investigate the estimates of the forced response in observations. Our results show that many different types of methods are skillful for estimating the forced response in climate models, though the relative skill of individual methods varies depending on the variable and evaluation metric. Methods with comparable skill in models can give a wide range of estimates of the forced response pattern in observations, illustrating the epistemic uncertainty in forced response estimates. ForceSMIP gives new insights into the forced response in observations, its uncertainty, and methods for its estimation.
Abstract Since the late 1970s, observations have shown a widening of the tropical Hadley cell (HC) circulation. State‐of‐the‐art climate models reproduce the general trend along with a projected continuous expansion. Discrepancies in expansion rates of observation‐ and model‐based studies have been attributed to differences in applied methods, natural variability and model shortcomings. Furthermore, the driving influence of tropical or extratropical processes on these changes is not well understood. All of this highlights the dynamical mechanisms and the region of origin controlling the tropical width are still insufficiently understood. Here we examine the influence of systematic model biases of the atmosphere‐only Unified Model (UM) onto the simulation of the Southern Hemisphere (SH) tropical edge. We utilize nudged experiments with prescribed sea surface temperatures, where potential temperature and horizontal winds are relaxed back to ERA‐Interim reanalysis for a 20‐year period in selected regions. Correcting model biases in the tropics and extratropics separately allows us to dissect the dominant remote impacts of present model errors onto the SH tropical edge simulation. The experiments are applied to established tropical width metrics ranging from near‐surface to upper‐level metrics capturing the poleward flank of the HC. We find both regions work remotely to reduce errors in the UM fields and location of the tropical edge. Surprisingly, correcting the extratropical biases, south of 45°S, more consistently improves the tropical width across the metrics and seasons than nudging the tropics (10°N–10°S). These findings demonstrate the substantial role of extratropical influences in locating the SH tropical edge.
Since the early 2010s, anthropogenic aerosols have started decreasing in East Asia (EA) while have continued to increase in South Asia (SA). Yet the climate impacts of this Asian aerosol dipole (AAD) pattern remain largely unknown. Using a state-of-the-art climate model, we demonstrate that the climate response is distinctly different between the SA aerosol increases and EA aerosol decreases. The SA aerosol increases lead to ~2.7 times stronger land summer precipitation change within the forced regions than the EA aerosol decreases. Contrastingly, the SA aerosol increases, within the tropical monsoon regime, produce weak and tropically confined responses, while the EA aerosol decreases yield a pronounced northern hemisphere warming aided by extratropical mean westerly and positive air-sea feedbacks over the western North Pacific. By scaling the observed instantaneous shortwave radiative forcing, we reveal that the recent AAD induces a pronounced northern hemisphere extratropical (beyond 30°N) warming (0.024 ± 0.010 °C decade -1 ), particularly over Europe (0.049 ± 0.009 °C decade -1 ). These findings highlight the importance of the pattern effect of forcings in driving global climate and have important implications for decadal prediction.
The Observations for Model Intercomparison Project (Obs4MIPs) was initiated in 2010 to facilitate the use of observations in climate model evaluation and research, with a particular target being the Coupled Model Intercomparison Project (CMIP), a major initiative of the World Climate Research Programme (WCRP). To this end, Obs4MIPs (1) targets observed variables that can be compared to CMIP model variables; (2) utilizes dataset formatting specifications and metadata requirements closely aligned with CMIP model output; (3) provides brief technical documentation for each dataset, designed for nonexperts and tailored towards relevance for model evaluation, including information on uncertainty, dataset merits, and limitations; and (4) disseminates the data through the Earth System Grid Federation (ESGF) platforms, making the observations searchable and accessible via the same portals as the model output. Taken together, these characteristics of the organization and structure of obs4MIPs should entice a more diverse community of researchers to engage in the comparison of model output with observations and to contribute to a more comprehensive evaluation of the climate models. At present, the number of obs4MIPs datasets has grown to about 80; many are undergoing updates, with another 20 or so in preparation, and more than 100 are proposed and under consideration. A partial list of current global satellite-based datasets includes humidity and temperature profiles; a wide range of cloud and aerosol observations; ocean surface wind, temperature, height, and sea ice fraction; surface and top-of-atmosphere longwave and shortwave radiation; and ozone (O 3 ), methane (CH 4 ), and carbon dioxide (CO 2 ) products. A partial list of proposed products expected to be useful in analyzing CMIP6 results includes the following: alternative products for the above quantities, additional products for ocean surface flux and chlorophyll products, a number of vegetation products (e.g., FAPAR, LAI, burned area fraction), ice sheet mass and height, carbon monoxide (CO), and nitrogen dioxide (NO 2 ). While most existing obs4MIPs datasets consist of monthly-mean gridded data over the global domain, products with higher time resolution (e.g., daily) and/or regional products are now receiving more attention. Along with an increasing number of datasets, obs4MIPs has implemented a number of capability upgrades including (1) an updated obs4MIPs data specifications document that provides additional search facets and generally improves congruence with CMIP6 specifications for model datasets, (2) a set of six easily understood indicators that help guide users as to a dataset's maturity and suitability for application, and (3) an option to supply supplemental information about a dataset beyond what can be found in the standard metadata. With the maturation of the obs4MIPs framework, the dataset inclusion process, and the dataset formatting guidelines and resources, the scope of the observations being considered is expected to grow to include gridded in situ datasets as well as datasets with a regional focus, and the ultimate intent is to judiciously expand this scope to any observation dataset that has applicability for evaluation of the types of Earth system models used in CMIP.
It remains a challenge for climate models to correctly capture the relationship between precipitation and ENSO. This studyexamines the linkage between the simulated precipitation climatology and ENSO-related precipitation anomaly during borealwinter based on the multi-model ensemble from the Atmospheric Model Intercomparison Project Phase 5 (AMIP5) and perturbedparameter ensemble (PPE) with the Beijing Climate Center (BCC) atmospheric model. The AMIP5 models have large biases insimulating the tropical precipitation anomaly during El Niño, such as the shifts of the Inter-Tropical Convergence Zone (ITCZ)and South Pacific Convergence Zone (SPCZ). The inter-model differences show that the precipitation change in response to seasurface temperature (SST) change increases with enhanced precipitation climatology. The ENSO-related precipitation anomalycan also be related to the spatial distribution of the mean-state precipitation. The simulated ITCZ/SPCZ displacements aresignificantly correlated with the spatial precipitation-SST relationship in the mean state. Models with stronger mean-stateprecipitation-SST relationship also produce stronger SPCZ/ITCZ displacements. The mean-state tropical precipitation also hasstrong impacts on the ENSO-related precipitation anomalies over East Asia. In the BCC PPE, the connection between the mean-state precipitation and ENSO-related precipitation anomaly is overall consistent with that in AMIP5. Parameters associated withthe low-cloud and deep convection processes are the most influential ones for the precipitation simulations in BCC. Comparedwith the version in AMIP5, the new BCC model can better simulate the precipitation climatology and the relationship betweenENSO and precipitation over southern China. These results have important implications for model development and model-basedclimate predictions.
Abstract Hyperparameter optimization (HPO) is an important step in machine learning (ML) model development, but common practices are archaic—primarily relying on manual or grid searches. This is partly because adopting advanced HPO algorithms introduces added complexity to the workflow, leading to longer computation times. This poses a notable challenge to ML applications, as suboptimal hyperparameter selections curtail the potential of ML model performance, ultimately obstructing the full exploitation of ML techniques. In this article, we present a two-step HPO method as a strategic solution to curbing computational demands and wait times, gleaned from practical experiences in applied ML parameterization work. The initial phase involves a preliminary evaluation of hyperparameters on a small subset of the training dataset, followed by a reevaluation of the top-performing candidate models postretraining with the entire training dataset. This two-step HPO method is universally applicable across HPO search algorithms, and we argue it has attractive efficiency gains. As a case study, we present our recent application of the two-step HPO method to the development of neural network emulators for aerosol activation. Although our primary use case is a data-rich limit with many millions of samples, we also find that using up to 0.0025% of the data—a few thousand samples—in the initial step is sufficient to find optimal hyperparameter configurations from much more extensive sampling, achieving up to 135× speedup. The benefits of this method materialize through an assessment of hyperparameters and model performance, revealing the minimal model complexity required to achieve the best performance. The assortment of top-performing models harvested from the HPO process allows us to choose a high-performing model with a low inference cost for efficient use in global climate models (GCMs).
A nocturnal Amazonian low-level jet (ALLJ) was recently diagnosed using reanalysis data. This work assesses the ability of CESM1.2.2 to reproduce the jet and explores the mechanisms by which the ALLJ influences convection in the Amazon. The coupled CESM simulates the nocturnal ALLJ realistically, while CAM5 does not. A low-level cold air temperature bias in the eastern Amazon exists in CAM5, and thus the ALLJ is weaker than observed. However, a cold SST bias over the equatorial North Atlantic in the coupled model offsets the cold air temperature bias, producing a realistic ALLJ. Climate models significantly underestimate March–May (MAM) precipitation over the eastern Amazon. Here, we ran two sensitivity experiments using the coupled CESM by adding bottom-heavy diabatic heating at noon and midnight for 2.5 h along the coastal Amazon during MAM to mimic the occurrence of shallow precipitating convection. When heating is added during the early afternoon, coastal convection deepens and the ALLJ transports moisture inland from the ocean, preconditioning the environment for deep convective development during the ensuing hours. The increased convection over the eastern Amazon also moderately alleviates the equatorial Atlantic westerly wind bias, leading to deepening of the east Atlantic thermocline in the following months and partially improving the simulated June–August (JJA) Atlantic cold tongue in the coupled model. When heating is added at night, coastal convection does not strengthen as much and the ALLJ transports less moisture. Improvements in the simulated Atlantic winds and SST are negligible. Therefore, diurnal circulations matter to the organization of convection and rain across the Amazon, with impacts over the equatorial Atlantic.
Abstract This work documents version two of the Department of Energy's Energy Exascale Earth System Model (E3SM). E3SMv2 is a significant evolution from its predecessor E3SMv1, resulting in a model that is nearly twice as fast and with a simulated climate that is improved in many metrics. We describe the physical climate model in its lower horizontal resolution configuration consisting of 110 km atmosphere, 165 km land, 0.5° river routing model, and an ocean and sea ice with mesh spacing varying between 60 km in the mid‐latitudes and 30 km at the equator and poles. The model performance is evaluated with Coupled Model Intercomparison Project Phase 6 Diagnosis, Evaluation, and Characterization of Klima simulations augmented with historical simulations as well as simulations to evaluate impacts of different forcing agents. The simulated climate has many realistic features of the climate system, with notable improvements in clouds and precipitation compared to E3SMv1. E3SMv1 suffered from an excessively high equilibrium climate sensitivity (ECS) of 5.3 K. In E3SMv2, ECS is reduced to 4.0 K which is now within the plausible range based on a recent World Climate Research Program assessment. However, a number of important biases remain including a weak Atlantic Meridional Overturning Circulation, deficiencies in the characteristics and spectral distribution of tropical atmospheric variability, and a significant underestimation of the observed warming in the second half of the historical period. An analysis of single‐forcing simulations indicates that correcting the historical temperature bias would require a substantial reduction in the magnitude of the aerosol‐related forcing.
Abstract In numerical weather prediction and climate models, boundary‐layer clouds are controlled by a wide range of subgrid‐scale processes. However, understanding the nature of these processes and their role in the evolution of the cloud size distribution as a whole has been elusive. To address this issue, we adopt a novel empirical framework from the field of population dynamics to model the evolution of cloud size statistics by using the shallow cumulus properties obtained from a large‐eddy simulation (LES). Our approach involves representing the cloud size distribution and the total cloud area using a revised Lotka‐Volterra model and ridge linear model, respectively. The physical interpretation of the total cloud area and coefficients obtained from the optimization of the models reveals three stages probably interpreted by dominant processes: the formation of new clouds, the growth of single clouds, and a steady state with organized transitions involving the growth and decay of multiple clouds. Furthermore, we showcase the potential of this framework to serve as a component of scale‐aware parameterizations of shallow‐convective clouds in atmospheric models.
This project brought DOE’s climate modeling effort with the Energy Exascale Earth System Model into the Polar Amplification Model Intercomparison Project (PAMIP), which is part of the sixth and latest Coupled Model Intercomparison Project, CMIP6. PAMIP examines the causes and consequences of polar amplification, when external forcing results in a larger temperature increase in high latitudes than the global average, in a coordinated set of model experiments with a common modeling protocol. Our teams from UC Irvine and the University of Toronto have designed, carried out, analyzed, and disseminated PAMIP output from the Energy Exascale Earth System Model (E3SM) and the Community Earth System Model (CESM). PAMIP’s ongoing stream of significant new results have advanced progress in the community’s understanding and led to new outstanding research questions that have motivated further work. PAMIP has led to improved consensus on the atmospheric response to sea ice loss. The important finding is that for a similar sea ice anomaly forcing, the simulated atmospheric response in the troposphere is remarkably consistent among the 16 models’ runs analyzed. The zonal-mean tropospheric response consists of a very robust equatorward shift of the westerly flow in mid-latitudes. However, while the multi-model mean response is robust, it has a weak amplitude relative to internal variability. We identified a weakness in the models (including E3SM) in terms of their response to sea-ice forcing that is related to eddy forcing (or nonlinear dynamical effects) at mid- to high latitudes. In fact, E3SM is an outlier in terms of the models participating, and in that sense, it turned out to be a vital participant model. We found that reductions in energy transport due sea-ice loss and involving dry air only are compensated by increases in moist energy transport from warmer sea surface temperature in midlatitudes. This leads to a large spread in energy transport into the Arctic and is a potential source of spread in Arctic amplification. We identified an important role that climate modes, including tropical modes of variability (El Nino and the Southern Oscillation (ENSO); the Quasi-biennial Oscillation (QBO)) play in the response to sea-ice anomalies, including in ocean coupled experiments. Similarly, we identified and quantified the contribution of sea-ice thickness to the atmospheric response compared to the response to sea-ice extent only. We found that it is important to run large ensembles and even with an ensemble size of 100 simulations the response is largely influenced by internal variability. We demonstrated convincingly that Ural blocking, not sea-ice loss, provides the weakening of the stratospheric polar vortex in fall/early winter and a negative phase of the North Atlantic Oscillation that can last for up to two months. However, sea-ice anomalies can influence the background flow so that the response to Ural blocking is more persistent under low sea-ice conditions in the Barents/Kara Sea than high sea ice Atmospheric model hierarchies that progressively add individual processes have a long history in providing dynamical insight for modeling the atmosphere. Similarly, coupled model hierarchies that progressively add individual ocean processes can provide insights into the workings of the coupled climate system, however such hierarchies have not been available except for a non-dynamical slab ocean model. Because of the missing processes, surface flux corrections must be added to produce a target climate. In this project, we managed to overcome this problem and develop a globally coupled ocean model hierarchy in CESM that can turn on and off the processes of mixed-layer entrainment and Ekman flow. We used the hierarchy to study the impact of Arctic sea-ice loss on the climate system. We find that the effect of mixed-layer entrainment on ocean heat uptake influences the atmospheric circulation by shifting the latitudinal positions of the mid-latitude westerly jet and the Intertropical Convergence Zone (ITCZ). In quadrupled CO 2 experiments, we studied how air-sea coupling affects the response of tropical rainfall under global warming. In order to identify the importance of individual ocean processes, we used the hierarchy of ocean models to separate the effects of seasonal mixed-layer entrainment, wind-driven Ekman flows, and frictional flows. We showed that including Ekman and frictional flows allows our simulation to produce the Pacific Ocean's enhanced equatorial warming pattern and equatorward ITCZ contraction noted in previous climate simulations. We also showed that the frictional flow, which has yet to receive much attention, is as important as the Ekman flow in generating equatorial heat convergence.
Climate models project higher temperatures and a substantial net water deficit during the summer months over Germany until the end of the 21st century. The agricultural sector is particularly vulnerable to changing weather and climate conditions. In this study, we examine the role of irrigation for potato yields as a climate change adaptation measure in Northeast Lower Saxony (Germany). The region represents the largest irrigated area in Germany and is one of the main growing areas of potatoes. We develop a statistical multivariate regression model to investigate whether regional climate and irrigation data at the monthly and county level are suitable to describe potato yield variability from 1978 to 2018. Irrigation is estimated based on the climatic water balance and crop water demand. These estimates are calibrated against irrigation data from local agencies to account for realistic monthly water withdrawals. We analyze future yield changes based on climate model projections for the 21st century and different irrigation scenarios. We find that the potato yield will decrease significantly by 18% on average until 2050 with climate projection uncertainty of ±6.4% under the high emission scenario RCP 8.5 and without irrigation. Irrigation at current levels could compensate for about 54% of these yield losses. Doubling the amount of irrigation minimizes yield losses by as much as 80%. Under the low (middle) emission scenario RCP 2.6 (RCP 4.5), potato yield are projected to decrease moderately by 4.5% ± 3.4% (5.4% ± 2.8%) with irrigation at current levels. We show that the yield losses mainly originate from higher temperatures during the summer. Our analysis indicates that the capacity of irrigation to minimize yield losses is limited. Strengthening climate change mitigation efforts may be more effective in preventing yield losses. However, the effect of elevated CO 2 might buffer the high yield losses, especially in RCP 8.5. The studies supports farmers and decision makers in assessing future risks in potato yields. It can guide policymakers in allocating economic resources to cope with climate change.