Search NASA⌕ Search

SEARCH · Search NASA

Results for “circulation 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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Are 100 Ensemble Members Enough to Capture the Remote Atmospheric Response to +2°C Arctic Sea Ice Loss?

This study presents results from the Polar Amplification Multimodel Intercomparison Project (PAMIP) single-year time-slice experiments that aim to isolate the atmospheric response to Arctic sea ice loss at global warming levels of +2°C. Using two general circulation models (GCMs), the ensemble size is increased up to 300 ensemble members, beyond the recommended 100 members. After partitioning the response in groups of 100 ensemble members, the reproducibility of the results is evaluated, with a focus on the response of the midlatitude jet streams in the North Atlantic and North Pacific. Both atmosphere-only and coupled ocean–atmosphere PAMIP experiments are analyzed. Substantial differences in the midlatitude response are found among the different experiment subsets, suggesting that 100-member ensembles are still significantly influenced by internal variability, which can mislead conclusions. Despite an overall stronger response, the coupled ocean–atmosphere runs exhibit greater spread due to additional ENSO-related internal variability when the ocean is interactive. The lack of consistency in the response is true for anomalies that are statistically significant according to Student’s t and false discovery rate tests. This is problematic for the multimodel assessment of the response, as some of the spread may be attributed to different model sensitivities whereas it is due to internal variability. Here, we propose a method to overcome this consistency issue that allows for more robust conclusions when only 100 ensemble members are used.

58 GEOSCIENCES↗
DOE OSTI2021

Using TRMM Latent Heat as a Source to Estimate Convection Induced Gravity Wave Momentum Flux in the Lower Stratosphere

Due to a lack of observations, it is a big challenge to quantify the momentum flux transported from the troposphere into the lower stratosphere by gravity waves generated by convection. This limits our understanding of the dynamics of general circulation in the stratosphere and the realism of simulations of the Quasi Biennial Oscillation. In the past decade, some general circulation models have linked the momentum flux to the latent heating from subgrid scale convective precipitation with some success. However, there is still a large uncertainty in the sources of subgrid scale convection in these models. This study applies the links between the momentum flux and convective latent heating directly to the convective precipitation derived from 16 years of TRMM precipitation radar observations in the tropics and subtropics. The total and directional momentum flux at 100 hPa are derived at individual convective pixels by using the TRMM Spectrum Latent Heating product and large-scale wind profiles from the ERA-Interim dataset. For the first time, we are able to estimate the geographical distribution of momentum flux at 100 hPa in tropics and subtropics from the observed convective sources. The diurnal, seasonal, and interannual variations of the derived momentum flux are presented and discussed. These results could provide a reference to validate the dynamic coupling between the tropical troposphere and the stratosphere in the general circulation models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

The Role of Radiative Interactions in Tropical Cyclone Development under Realistic Boundary Conditions

Abstract The impact of radiative interactions on tropical cyclone (TC) climatology is investigated using a global, TC-permitting general circulation model (GCM) with realistic boundary conditions. In this model, synoptic-scale radiative interactions are suppressed by overwriting the model-generated atmospheric radiative cooling rates with their monthly varying climatological values. When radiative interactions are suppressed, the global TC frequency is significantly reduced, indicating that radiative interactions are a critical component of TC development even in the presence of spatially varying boundary conditions. The reduced TC activity is primarily due to a decrease in the frequency of pre-TC synoptic disturbances (“seeds”), whereas the likelihood that the seeds undergo cyclogenesis is less affected. When radiative interactions are suppressed, TC genesis shifts toward coastal regions, whereas TC lysis locations stay almost unchanged; together the distance between genesis and lysis is shortened, reducing TC duration. In a warmer climate, the magnitude of TC reduction from suppressing radiative interactions is diminished due to the larger contribution from latent heat release with increased sea surface temperatures. These results highlight the importance of radiative interactions in modulating the frequency and duration of TCs.

Meteorology & Atmospheric Sciences↗
DOE OSTI2023

Adaptation to climate change through strategic integration of long fallow into cropping system in a dryland Mediterranean-type environment

The crop-growing region of Western Australia characterized by a Mediterranean-type climate is projected to become warmer and drier. Appropriate selection of crop sequences will be of importance to cope with these climatic changes for this largest grain-producing region of Australia. Through linking a widely used crop model (APSIM), 26 General Circulation Models (GCMs) with one Shared Socioeconomic Pathway (SSP585) and economic analysis, we explored how the climate change would affect dryland wheat cropping and whether/how long fallow (the practice of leaving a field out of production for an entire growing season) could be integrated into wheat cropping system in Western Australia. The potential adaptation of long fallow into wheat system was assessed with four fixed rotations (fallow-wheat, fallow-wheat-wheat, fallow-wheat-wheat-wheat, and fallow-wheat-wheat-wheat-wheat) and four flexible sowing rule-based rotations (the land was fallowed if sowing rule was not met), compared with continuous wheat. The simulation results at four representing locations show that climate change would have negative impacts on both yield and economic return of continuous wheat cropping in Western Australia. Wheat after fallow out-yielded and out-profited wheat after wheat under future climate. But integrating fallow into wheat cropping systems with the above fixed rotations would lead to yield and economic loss. By contrast, cropping systems in which fallowing took place when sowing condition could not be met at a certain time would achieve comparable yield and economic return to continuous wheat, with wheat yield being only 5% less than continuous wheat and the gross margin being $12ha -1 more than continuous wheat averaged across locations. We highlight strategic integration of long fallow into cropping system in a dryland Mediterranean-type environment would have a great potential to cope with future climate change. These findings can be extended into other Mediterranean-type cropping regions in Australia and beyond.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

On the Dependence of Simulated Convection on Domain Size in CRMs

Abstract We present a heuristic model to explain the suppression of deep convection in convection‐resolving models (CRMs) with a small number of grid columns, such as those used in super‐parameterized or multi‐scale modeling framework (MMF) general circulation models (GCM) of the atmosphere. Domains with few grid columns require greater instability to sustain convection because they force a large convective fraction, driving strong compensating subsidence warming. Updraft dilution, which is stronger for reduced horizontal grid spacing, enhances this effect. Thus, suppression of deep convection in CRMs with few grid columns can be reduced by increasing grid spacing. Radiative‐convective equilibrium simulations using standalone CRM simulations with the System for Atmospheric Modeling (SAM) and using GCM‐coupled CRM simulations with the Energy Exascale Earth System Model (E3SM)‐MMF confirm the heuristic model results.

CRM↗
DOE OSTI2021

Retrospective Analysis and Bayesian Model Averaging of CMIP6 Precipitation in the Nile River Basin

The Nile River basin is one of the global hotspots vulnerable to climate change impacts because of a fast-growing population and geopolitical tensions. Previous studies demonstrated that general circulation models (GCMs) frequently show disagreement in the sign of change in annual precipitation projections. Here, we first evaluate the performance of 20 GCMs from phase six of the Coupled Model Intercomparison Project (CMIP6) benchmarked against a high-spatial-resolution precipitation dataset dating back to 1983 from Precipitation Estimation from Remotely Sensed Information Using Artificial Neural Networks–Climate Data Record (PERSIANN-CDR). Next, a Bayesian model averaging (BMA) approach is adopted to derive probability distributions of precipitation projections in the Nile basin. Retrospective analysis reveals that most GCMs exhibit considerable (up to 64% of mean annual precipitation) and spatially heterogenous bias in simulating annual precipitation. Moreover, it is shown that all GCMs underestimate interannual variability; thus, the ensemble range is underdispersive and is a poor indicator of uncertainty. The projected changes from the BMA model show that the value and sign of change vary considerably across the Nile basin. Specifically, it is found that projected changes in the two headwaters basins, namely, the Blue Nile and Upper White Nile, are 0.03% and -1.65%, respectively; both are statistically insignificant at α = 0.05. The uncertainty range estimated from the BMA model shows that the probability of a precipitation decrease is much higher in the Upper White Nile basin whereas projected change in the Blue Nile is highly uncertain both in magnitude and sign of change.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Convection-Permitting Simulations With the E3SM Global Atmosphere Model

This paper describes the first implementation of the Δx = 3.25 km version of the Energy Exascale Earth System Model (E3SM) global atmosphere model and its behavior in a 40-day prescribed-sea-surface-temperature simulation (January 20 through February 28, 2020). This simulation was performed as part of the DYnamics of the Atmospheric general circulation Modeled On Non-hydrostatic Domains (DYAMOND) Phase 2 model intercomparison. Effective resolution is found to be ~ 6x the horizontal dynamics grid resolution despite using a coarser grid for physical parameterizations. Despite this new model being in an immature and untuned state, moving to 3.25 km grid spacing solves several long-standing problems with the E3SM model. In particular, Amazon precipitation is much more realistic, the frequency of light and heavy precipitation is improved, agreement between the simulated and observed diurnal cycle of tropical precipitation is excellent, and the vertical structure of tropical convection and coastal stratocumulus look good. In addition, the new model is able to capture the frequency and structure of important weather events (e.g., tropical cyclones, extratropical cyclones including atmospheric rivers, and cold air outbreaks). Interestingly, this model does not get rid of the erroneous southern branch of the intertropical convergence zone nor the tendency for strongest convection to occur over the Maritime Continent rather than the West Pacific, both of which are classic climate model biases. Several other problems with the simulation are identified, underscoring the fact that this model is a work in progress.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A review of future weather data for assessing climate change impacts on buildings and energy systems

The effectiveness of climate change impact assessments and the development of adaptation strategies depend on the availability of high-quality future weather data. However, significant gaps exist between the needs of the energy research community and the focus of the climate modeling community, primarily due to a historical lack of communication and collaboration between the two groups. Here, to address this issue, this work provides a comprehensive overview of the critical aspects involved in creating future weather data for building and energy system modeling, including emissions scenarios, general circulation models, downscaling methods, categories of future weather data, and uncertainties in climate simulations. Moreover, it critically evaluates the applicability and suitability of various types of future weather data in five key application scenarios: energy use analysis, resilience analysis, HVAC design, utility-scale analysis, and renewable energy analysis. Finally, this work presents recommendations for high-level actions and research directions to foster collaboration between the energy research and climate modeling communities and to promote the integration of future weather data into energy codes and the design practices of buildings and energy systems.

Climate change↗
DOE OSTI2020

Clouds and Convective Self-Aggregation in a Multimodel Ensemble of Radiative-Convective Equilibrium Simulations

The Radiative-Convective Equilibrium Model Intercomparison Project (RCEMIP) is an intercomparison of multiple types of numerical models configured in radiative-convective equilibrium (RCE). RCE is an idealization of the tropical atmosphere that has long been used to study basic questions in climate science. Here, we employ RCE to investigate the role that clouds and convective activity play in determining cloud feedbacks, climate sensitivity, the state of convective aggregation, and the equilibrium climate. RCEMIP is unique among intercomparisons in its inclusion of a wide range of model types, including atmospheric general circulation models (GCMs), single column models (SCMs), cloud-resolving models (CRMs), large eddy simulations (LES), and global cloud-resolving models (GCRMs). The first results are presented from the RCEMIP ensemble of more than 30 models. While there are large differences across the RCEMIP ensemble in the representation of mean profiles of temperature, humidity, and cloudiness, in a majority of models anvil clouds rise, warm, and decrease in area coverage in response to an increase in sea surface temperature (SST). Nearly all models exhibit self-aggregation in large domains and agree that self-aggregation acts to dry and warm the troposphere, reduce high cloudiness, and increase cooling to space. The degree of self-aggregation exhibits no clear tendency with warming. There is a wide range of climate sensitivities, but models with parameterized convection tend to have lower climate sensitivities than models with explicit convection. In models with parameterized convection, aggregated simulations have lower climate sensitivities than unaggregated simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Observations of Gravity Waves in the Middle Atmosphere of Mars

Gravity waves are ubiquitous throughout the atmosphere of Mars. Their propagation and dissipation influence the circulation and thermal structure of the middle and upper atmosphere. Yet there have been few studies of gravity wave characteristics in the middle atmosphere, a region that is critical for their propagation from generation in the lower atmosphere to dissipation and associated exchanges of momentum and energy in the upper atmosphere. There have been few studies because few atmospheric profiles span the middle atmosphere with the kilometer-scale or smaller vertical resolution that is Required to characterize gravity waves. Here we report the characterization of gravity waves in the middle atmosphere using 12 high-resolution atmospheric density profiles. Four of these were acquired from a ground-based stellar occultation from 1976 that yielded immersion and emersion profiles on opposite sides of the planet, and the remaining eight were measured during atmospheric entry by landers and rovers. Predominant wavelengths were 3–14 km, and amplitudes were generally 0.8%–2.5%. Where static stability is large and positive, gravity waves grow efficiently. In other instances, static stability is not large and positive over a wide altitude range, and gravity wave amplitudes do not behave as neatly. These observations of gravity waves in the middle atmosphere of Mars can be used to test gravity wave parameterizations in large-scale general circulation models and to investigate predictions for how gravity wave propagation and dissipation influence the circulation and thermal structure of the middle and upper atmosphere.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Overestimated Eddy Kinetic Energy in the Eddy‐Rich Regions Simulated by Eddy‐Resolving Global Ocean–Sea Ice Models

Abstract The performance of eddy‐resolving global ocean–sea ice models in simulating mesoscale eddies is evaluated using six eddy‐resolving experiments forced by different atmospheric reanalysis products. Interestingly, eddy‐resolving ocean general circulation models (OGCMs) tend to simulate more (less) energetic eddy‐rich (eddy‐poor) regions with a smaller (larger) spatial extent than satellite observation, which finally shows that larger (smaller) mesoscale energy intensity (EI) is simulated in the eddy‐rich (eddy‐poor) regions. Quantitatively, there is an approximately 27%–60% overestimation of EI in the eddy‐rich regions, which are mainly located in the Kuroshio–Oyashio Extension, the Gulf Stream, and the Antarctic Circumpolar Currents regions, although the global mean EI is underestimated by 25%–45%. Apparently, the eddy kinetic energy in the eddy‐poor region is underestimated. Further analyses based on coherent mesoscale eddy properties show that the overestimation in the eddy‐rich regions is mainly attributed to mesoscale eddies’ intensity and is more prominent when mesoscale eddies are in their growth stage.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Nudging allows direct evaluation of coupled climate models with in situ observations: a case study from the MOSAiC expedition

Comparing the output of general circulation models to observations is essential for assessing and improving the quality of models. While numerical weather prediction models are routinely assessed against a large array of observations, comparing climate models and observations usually requires long time series to build robust statistics. Here, we show that by nudging the large-scale atmospheric circulation in coupled climate models, model output can be compared to local observations for individual days. We illustrate this for three climate models during a period in April 2020 when a warm air intrusion reached the MOSAiC (Multidisciplinary drifting Observatory for the Study of Arctic Climate) expedition in the central Arctic. Radiosondes, cloud remote sensing and surface flux observations from the MOSAiC expedition serve as reference observations. The climate models AWI-CM1/ECHAM and AWI-CM3/IFS miss the diurnal cycle of surface temperature in spring, likely because both models assume the snowpack on ice to have a uniform temperature. CAM6, a model that uses three layers to represent snow temperature, represents the diurnal cycle more realistically. During a cold and dry period with pervasive thin mixed-phase clouds, AWI-CM1/ECHAM only produces partial cloud cover and overestimates downwelling shortwave radiation at the surface. AWI-CM3/IFS produces a closed cloud cover but misses cloud liquid water. Our results show that nudging the large-scale circulation to the observed state allows a meaningful comparison of climate model output even to short-term observational campaigns. We suggest that nudging can simplify and accelerate the pathway from observations to climate model improvements and substantially extends the range of observations suitable for model evaluation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Horizontal Resolution Sensitivity of the Simple Convection-Permitting E3SM Atmosphere Model in a Doubly-Periodic Configuration

We develop a doubly periodic version of the Simple Convection-Permitting E3SM Atmosphere Model (SCREAM) to provide an efficient configuration for this global convection permitting model (GCPM), akin to a single column model often found in conventional general circulation models. The design details are explained, in addition to the extensive case library associated with the doubly periodic SCREAM (DP-SCREAM) configuration. We demonstrate that doubly periodic cloud resolving models are useful tools to explore the horizontal resolution sensitivity of GCPMs, in addition to replicating biases seen in the global models. Using DP-SCREAM, we show that SCREAM exhibits behaviors of a scale aware model as it is able to naturally partition between sub-grid scale (SGS) and resolved vertical transport across the gray zone of turbulence. We show that SCREAM is reasonably scale insensitive when run at resolutions from 1 to 5 km, but can exhibit sensitivity, particularly for the shallow convective regime, when run at resolutions approaching that of large eddy simulations. We conclude that SGS parameterization improvements are likely needed to reduce this scale sensitivity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Influence of sea-ice anomalies on Antarctic precipitation using source attribution in the Community Earth System Model

Abstract. We conduct sensitivity experiments using a general circulation model that has an explicit water source tagging capability forced by prescribed composites of pre-industrial sea-ice concentrations (SICs) and corresponding sea surface temperatures (SSTs) to understand the impact of sea-ice anomalies on regional evaporation, moisture transport and source–receptor relationships for Antarctic precipitation in the absence of anthropogenic forcing. Surface sensible heat fluxes, evaporation and column-integrated water vapor are larger over Southern Ocean (SO) areas with lower SICs. Changes in Antarctic precipitation and its source attribution with SICs have a strong spatial variability. Among the tagged source regions, the Southern Ocean (south of 50° S) contributes the most (40 %) to the Antarctic total precipitation, followed by more northerly ocean basins, most notably the South Pacific Ocean (27%), southern Indian Ocean (16 %) and South Atlantic Ocean (11 %). Comparing two experiments prescribed with high and low pre-industrial SICs, respectively, the annual mean Antarctic precipitation is about 150 Gt yr-1 (or 6 %) more in the lower SIC case than in the higher SIC case. This difference is larger than the model-simulated interannual variability in Antarctic precipitation (99 Gt yr-1). The contrast in contribution from the Southern Ocean, 102 Gt yr-1, is even more significant compared to the interannual variability of 35 Gt yr-1 in Antarctic precipitation that originates from the Southern Ocean. The horizontal transport pathways from individual vapor source regions to Antarctica are largely determined by large-scale atmospheric circulation patterns. Vapor from lower-latitude source regions takes elevated pathways to Antarctica. In contrast, vapor from the Southern Ocean moves southward within the lower troposphere to the Antarctic continent along moist isentropes that are largely shaped by local ambient conditions and coastal topography. This study also highlights the importance of atmospheric dynamics in affecting the thermodynamic impact of sea-ice anomalies associated with natural variability on Antarctic precipitation. Our analyses of the seasonal contrast in changes of basin-scale evaporation, moisture flux and precipitation suggest that the impact of SIC anomalies on regional Antarctic precipitation depends on dynamic changes that arise from SIC–SST perturbations along with internal variability. The latter appears to have a more significant effect on the moisture transport in austral winter than in summer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A new spatially distributed added value index for regional climate models: the EURO-CORDEX and the CORDEX-CORE highest resolution ensembles

The added value of using regional climate models (RCMs) to downscale data from general circulation models (GCMs) has often been questioned and researched. Although several studies have used different methods to identify (and in some cases quantify) the added value, there is still a need to find a general metric that quantifies the added value of any variable. This paper builds on past studies to propose a new metric of added value in the simulation of present-day climate which measures the difference in the probability density functions (PDFs) at each grid-cell between a model and an observation source, and then compares the results of the RCM and GCM in order to spatially compute the added value index. The same method is also adapted to quantify the climate change downscaling signal in a way that is consistent with the present-day metric. These new metrics are tested on the daily precipitation output from the EURO-CORDEX and CORDEX-CORE projection ensembles and reveal an overall positive added value of RCMs, especially at the tail-end of the distribution. Higher added value is obtained in areas of complex topography and coast-lines, as well as in tropical regions. Areas with large added value in present-day climate are consistent with areas of significant climate change downscaling signal in the RCP 8.5 far future simulations, and when the analysis is repeated at a low-resolution. The use of different resolution observations shows that the added value tends to decrease when models are compared to low-resolution observation datasets.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Ensemble modeling of global lake evaporation under climate change

Global projections of lake evaporation are typically based on simulations using single mechanistic models. However, because of its complex interactions with various lake physical properties, environmental and anthropogenic drivers, lake evaporation is highly variable and sensitive to the choice of model used. In this study, we present a multi-model analysis to investigate differences across global simulations of lake evaporation during the warm-season using three different lake models driven by outputs from four general circulation models (GCM) (i.e. 12 model combinations in total) for historic and future scenarios. Our results suggest substantial differences among lake-climate model simulations of lake evaporation. These differences varied throughout the 20th and 21st century, with model driver data explaining 74% of the variance in future projections of warm- season lake evaporation. Our projections indicate that, by the end of the 21st century (2070–2099), global annual lake evaporation rates will increase by 10–27% under Representative Concentration Pathways (RCPs) 2.6–8.5. Finally, we highlight the importance of using a multi-model approach for the prediction of future global lake evaporation responses to climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Improved Boundary Conditions for Coupled Geospace Models: An Application in Modeling Spacecraft Surface Charging Environment

Abstract Spacecraft surface charging in the inner magnetosphere often occurs in the pre‐midnight to the dawn sector when electron fluxes of tens of keV increase. Inner magnetosphere ring current models can be used to simulate/predict the spacecraft surface charging environment, with their outer boundary conditions specified either based on observations or provided by other models, such as MHD models. In the latter approach, using MHD quantities, the flux spectrum at the outer boundary is commonly assumed to follow a Kappa or Maxwellian distribution function. Such a method however often departs greatly from the realistic spectrum at E < tens of keV, a crucial energy range in the surface charging anomaly. In order to achieve a better representation of the surface charging environment, we propose to combine the MHD‐parameterized flux spectrum with an empirical electron flux model of E < 40 keV to set the electron flux boundary condition. Results indicate that as opposed to the case where the MHD‐parameterized flux distribution is solely used at the model boundary, simulations with the new boundary condition yields a more intense surface charging environment. The integrated electron flux between 10 < E < 50 keV, a measure of the severity of the surface charging environment, is significantly enhanced by 1‐2 orders of magnitude, leading to a much better agreement with Van Allen Probes measurements. This study hence demonstrates a reasonable solution to the setting of outer boundary conditions for inner magnetosphere models and is recommended for coupled geospace circulation models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Climate projections over the Great Lakes Region: using two-way coupling of a regional climate model with a 3-D lake model

Abstract. Warming trends in the Laurentian Great Lakes and surrounding areas have been observed in recent decades, and concerns continue to rise about the pace and pattern of future climate change over the world's largest freshwater system. To date, most regional climate models used for Great Lakes projections either neglected the lake-atmosphere interactions or are only coupled with a 1-D column lake model to represent the lake hydrodynamics. This study presents a Great Lakes climate change projection that has employed the two-way coupling of a regional climate model with a 3-D lake model (GLARM) to resolve 3-D hydrodynamics essential for large lakes. Using the three carefully selected Coupled Model Intercomparison Project Phase 5 (CMIP5) general circulation models (GCMs), we show that the GLARM ensemble average substantially reduces surface air temperature and precipitation biases of the driving GCM ensemble average in present-day climate simulations. The improvements are not only displayed from an atmospheric perspective but are also evident in the accurate simulations of lake temperature and ice coverage. We further present the GLARM projected climate change for the mid-21st century (2030–2049) and the late 21st century (2080–2099) in the Representative Concentration Pathway (RCP) 4.5 and RCP 8.5 scenarios. Under RCP 8.5, the Great Lakes basin is projected to warm by 1.3–2.1 ∘C by the mid-21st century and 4.1–5.0 ∘C by the end of the century relative to the early century (2000–2019). Moderate mitigation (RCP 4.5) reduces the mid-century warming to 0.8–1.8 ∘C and late-century warming to 1.8–2.7 ∘C. Annual precipitation in GLARM is projected to increase for the entire basin, varying from 0 % to 13 % during the mid-century and from 9 % to 32 % during the late century in different scenarios and simulations. The most significant increases are projected in spring and fall when current precipitation is highest and a minimal increase in winter when it is lowest. Lake surface temperatures (LSTs) are also projected to increase across the five lakes in all of the simulations, but with strong seasonal and spatial variability. The most significant LST increases occur in Lakes Superior and Ontario. The strongest warming is projected in spring that persists into the summer, resulting from earlier and more intense stratification in the future. In addition, diminishing winter stratification in the future suggests the transition from dimictic lakes to monomictic lakes by the end of the century. In contrast, a relatively smaller increase in LSTs during fall and winter is projected with heat transfer to the deep water due to the strong mixing and energy required for ice melting. Correspondingly, the highest monthly mean ice cover is projected to reduce to 3 %–15 % and 10 %–40 % across the lakes by the end of the century in RCP 8.5 and RCP 4.5, respectively. In the coastal regions, ice duration is projected to decrease by up to 60 d.

Xue, Pengfei (ORCID:000000025702421X)↗