Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Modeling and simulation”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS2021

Land-Atmosphere Coupling at the U.S. Southern Great Plains: A Comparison on Local Convective Regimes Between ARM Observations, Reanalysis, and Climate Model Simulations

Using the 9-yr warm-season observations at the Atmospheric Radiation Measurement Southern Great Plains site, we assess the land-atmosphere (L-A) coupling in North American Regional Reanalysis (NARR) and two climate models: hindcasts with the Community Atmosphere Model version 5.1 by Cloud-Associated Parameterizations Testbed (CAM5-CAPT) and nudged runs with the Energy Exascale Earth System Model Atmosphere Model version 1 Regionally Refined Model (EAMv1-RRM). We focus on three local convective regimes and diagnose model behaviors using the Local Coupling metrics (Santanello et al. 2018). NARR agrees well with observations except a slightly warmer and drier surface with higher downwelling shortwave radiation and lower evaporative fraction. On clear-sky days, it shows warmer and drier early-morning conditions in both models with significant underestimates in surface evaporation by EAMv1-RRM. On the majority of the ARM-observed shallow cumulus days, there is no or little low-level clouds in either model. When captured in models, the simulated shallow cumulus shows much less cloud fraction and lower cloud bases than observed. On the days with late-afternoon deep convection, models tend to present a stable early-morning lower atmosphere more frequently than the observations, suggesting that the deep convection is triggered more often by elevated instabilities. Generally, CAM5-CAPT can reproduce the local L-A coupling processes to some extent due to the constrained early-morning conditions and large-scale winds. EAMv1-RRM exhibits large precipitation deficits and warm and dry biases towards mid-to-late summers, which may be an amplification through a positive L-A feedback among initial atmosphere and land states, convection triggering and large-scale circulations.

land-atmosphere↗
DOE OSTI2021

Land–Atmosphere Coupling at the U.S. Southern Great Plains: A Comparison on Local Convective Regimes between ARM Observations, Reanalysis, and Climate Model Simulations

Using the 9-yr warm-season observations at the Atmospheric Radiation Measurement Southern Great Plains site, in this work we assess the land–atmosphere (LA) coupling in the North American Regional Reanalysis (NARR) and two climate models: hindcasts with the Community Atmosphere Model version 5.1 by Cloud-Associated Parameterizations Testbed (CAM5-CAPT) and nudged runs with the Energy Exascale Earth System Model Atmosphere Model version 1 Regionally Refined Model (EAMv1-RRM). We focus on three local convective regimes and diagnose model behaviors using the local coupling metrics. NARR agrees well with observations except a slightly warmer and drier surface with higher downwelling shortwave radiation and lower evaporative fraction. On clear-sky days, it shows warmer and drier early-morning conditions in both models with significant underestimates in surface evaporation by EAMv1-RRM. On the majority of the ARM-observed shallow cumulus days, there is no or little low-level clouds in either model. When captured in models, the simulated shallow cumulus shows much less cloud fraction and lower cloud bases than observed. On the days with late-afternoon deep convection, models tend to present a stable early-morning lower atmosphere more frequently than the observations, suggesting that the deep convection is triggered more often by elevated instabilities. Generally, CAM5-CAPT can reproduce the local LA coupling processes to some extent due to the constrained early-morning conditions and large-scale winds. EAMv1-RRM exhibits large precipitation deficits and warm and dry biases toward mid-to-late summers, which may be an amplification through a positive LA feedback among initial atmosphere and land states, convection triggering and large-scale circulations.

58 GEOSCIENCES↗
NASA NTRS2015

Benchmark Data Set for Wheat Growth Models: Field Experiments and AgMIP Multi-Model Simulations.

The data set includes a current representative management treatment from detailed, quality-tested sentinel field experiments with wheat from four contrasting environments including Australia, The Netherlands, India and Argentina. Measurements include local daily climate data (solar radiation, maximum and minimum temperature, precipitation, surface wind, dew point temperature, relative humidity, and vapor pressure), soil characteristics, frequent growth, nitrogen in crop and soil, crop and soil water and yield components. Simulations include results from 27 wheat models and a sensitivity analysis with 26 models and 30 years (1981-2010) for each location, for elevated atmospheric CO2 and temperature changes, a heat stress sensitivity analysis at anthesis, and a sensitivity analysis with soil and crop management variations and a Global Climate Model end-century scenario.

climate change impact↗
NASA NTRS2000

Coupled Model Simulation of Snowfall Events Over the Black Hills

Although many long-term simulations of snow accumulation and oblation have been made using stand-alone land surface models and surface models coupled with GCMs, less research has focused on short-term event simulations. Actually, accurate event simulations of snow-related processes are the basis for successful long-term simulation. Three advantages of event simulations of snowfall and snow melting are availability of: (1) intensive observation data from field experiments for validation; (2) more physically-realistic precipitation schemes for use in atmospheric models to simulate snowfall; and (3) a more detailed analysis of the snow melting processes. In addition to the complexities of snow related processes themselves, terrain-induced effects on snowfall/snow melting make simulations of snow events more difficult. Climatological observations indicate that terrain features such as the Black Hills of South Dakota and Wyoming can exert important effects on snow accumulation and snow oblation processes. One of the primary effects is that the orography causes forced uplift of airflow and causes atmospheric waves to form both upwind and downwind of it. Airflow often splits around the obstacle, converging on the lee side. This convergence may lead to precipitation enhancement. It also provides an elevated heat and moisture source that enhances atmospheric instability. During the period of April 5-May 5, 1999, the Upper Missouri River Basin Pilot Project (UMRBPP) made intensive observations on precipitation events occurring in the Black Hills. Two moderate snowfall events were captured during the period. The resulting high temporal and spatial resolution data provides opportunities to investigate terrain effects on snowfall amount, distribution, and melting. Successful simulation of snowfall amount, distribution, and evolution using atmospheric models is important to subsequent modeling of snow melting using snow sub-models in land surface schemes. In this paper, a coupled model system, consisting of an atmosphere model (ARPS) and a land-surface model (revised NCAR LSM), is used to simulate one of these UMRBPP cases.

Wang, Jianzhong↗
NASA NTRS1985

Transport processes in the stratosphere: Model simulations and comparisons with satellite observations

A three dimensional atmospheric model was used to study transport processes and to simulate the distribution of chemically active species in the stratosphere. The present results are part of a long term simulation of the seasonally varying structure of stratospheric trace constituents. The occurrence of a midwinter stratospheric warming and the associated transport of O3 and HNO3 during the simulation are described. Comparison of the simulated distributions of O3 and HNO3 are made with data from the Limb Infrared Monitor of the Stratosphere (LIMS) experiment. In addition, distributions of Ertel's potential vorticity on isentropic surfaces (IPV) are evaluated as a diagnostic for interpreting transport processes. Comparisons are made with IPV distributions inferred from LIMS temperature data.

Grose, W. L.↗
NASA NTRS2007

Statistical Analyses of Satellite Cloud Object Data from CERES: Comparison with Cloud-Resolving Model Simulations of Tropical Convective Clouds - Part III

The present study evaluates the ability of a cloud-resolving model (CRM) to simulate the physical properties of tropical deep convective cloud objects identified from a Clouds and the Earth s Radiant Energy System (CERES) data product. The emphasis of this study is the comparisons among the small-, medium- and large-size categories of cloud objects observed during March 1998 and between the large-size categories of cloud objects observed during March 1998 (strong El Ni o) and March 2000 (weak La Ni a). Results from the CRM simulations are analyzed in a way that is consistent with the CERES retrieval algorithm and they are averaged to match the scale of the CERES satellite footprints. Cloud physical properties are analyzed in terms of their summary histograms for each category. It is found that there is a general agreement in the overall shapes of all cloud physical properties between the simulated and observed distributions. Each cloud physical property produced by the CRM also exhibits different degrees of disagreement with observations over different ranges of the property. The simulated cloud tops are generally too high and cloud top temperatures are too low except for the large-size category of March 1998. The probability densities of the simulated top-of-the-atmosphere (TOA) albedos for all four categories are underestimated for high albedos, while those of cloud optical depth are overestimated at its lowest bin. These disagreements are mainly related to uncertainties in the cloud microphysics parameterization and inputs such as cloud ice effective size to the radiation calculation. Summary histograms of cloud optical depth and TOA albedo from the CRM simulations of the large-size category of cloud objects do not differ significantly between the March 1998 and 2000 periods, consistent with the CERES observations. However, the CRM is unable to reproduce the significant differences in the observed cloud top height while it overestimates the differences in the observed outgoing longwave radiation and cloud top temperature between the two periods. Comparisons between the CRM results and the observations for most parameters in March 1998 consistently show that both the simulations and observations have larger differences between the large- and small-size categories than between the large- and medium-size, or between the medium- and small-size categories. However, the simulated cloud properties do not change as much with size as observed. These disagreements are likely related to the spatial averaging of the forcing data and the mismatch in time and in space between the numerical weather prediction model from which the forcing data are produced and the CERES observed cloud systems.

Luo, Yali↗
NASA NTRS1991

Vertical eddy heat fluxes from model simulations

Vertical eddy fluxes of heat are calculated from simulations with a variety of climate models, ranging from three-dimensional GCMs to a one-dimensional radiative-convective model. The models' total eddy flux in the lower troposphere is found to agree well with Hantel's analysis from observations, but in the mid and upper troposphere the models' values are systematically 30 percent to 50 percent smaller than Hantel's. The models nevertheless give very good results for the global temperature profile, and the reason for the discrepancy is unclear. The model results show that the manner in which the vertical eddy flux is carried is very sensitive to the parameterization of moist convection. When a moist adiabatic adjustment scheme with a critical value for the relative humidity of 100 percent is used, the vertical transports by large-scale eddies and small-scale convection on a global basis are equal: but when a penetrative convection scheme is used, the large-scale flux on a global basis is only about one-fifth to one-fourth the small-scale flux. Comparison of the model results with observations indicates that the results with the latter scheme are more realistic. However, even in this case, in mid and high latitudes the large and small-scale vertical eddy fluxes of heat are comparable in magnitude above the planetary boundary layer.

Stone, Peter H.↗
DOE OSTI2018

AM Modeling & Simulation

A predictive simulation capability for AM is important to improve quality, decrease processing time and minimize costs. Has significant impact as LANL continues to be a leader in modeling and simulation for AM. This improved simulation capability will allow important AM specific physics questions to be answered for the first time. This approach, in combination with LANL's HPC resources will enable tractable, predictive simulation of AM and other solidification processes.

42 ENGINEERING↗
DOE OSTI2016

Bayesian optimization of the Community Land Model simulated biosphere-atmosphere exchange using CO 2 observations from a dense tower network and aircraft campaigns over Oregon

Here, the vast forests and natural areas of the Pacific Northwest comprise one of the most productive ecosystems in the northern hemisphere. The heterogeneous landscape of Oregon poses a particular challenge to ecosystem models. We present a framework using a scaling factor Bayesian inversion to improve the modeled atmosphere-biosphere exchange of carbon dioxide. Observations from 5 CO/CO 2 towers, eddy covariance towers, and airborne campaigns were used to constrain the Community Land Model CLM4.5 simulated terrestrial CO 2 exchange at a high spatial and temporal resolution (1/24°, 3-hourly). To balance aggregation errors and the degrees of freedom in the inverse modeling system, we applied an unsupervised clustering approach for the spatial structuring of our model domain. Data from flight campaigns were used to quantify the uncertainty introduced by the Lagrangian particle dispersion model that was applied for the inversions. The average annual statewide net ecosystem productivity (NEP) was increased by 32% to 29.7 TgC per year by assimilating the tropospheric mixing ratio data. The associated uncertainty was decreased by 28.4% to 29%, on average over the entire Oregon model domain with the lowest uncertainties of 11% in western Oregon. The largest differences between posterior and prior CO 2 fluxes were found for the Coast Range ecoregion of Oregon that also exhibits the highest availability of atmospheric observations and associated footprints. In this area, covered by highly productive Douglas-fir forest, the differences between the prior and posterior estimate of NEP averaged 3.84 TgC per year during the study period from 2012 through 2014.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Evaluation of top-down and bottom-up global terrestrial respiration estimates and their mismatch with model simulations

Terrestrial respiration is one of the most poorly understood processes in the global carbon cycle, making respiration predictions uncertain. However, expanding observations and machine learning approaches have led to a proliferation of estimates. We compiled total ecosystem and heterotrophic respiration estimates derived from top-down atmospheric inversions and bottom-up upscaling of ecosystem observations and compared them with dynamic vegetation models (DGVM) simulations over the 1980-2020 period. Our analysis revealed a convergence in mean annual global total ecosystem respiration estimates between top-down 97.1 (± SD 6.8) PgC yr-1 and bottom-up 98.5 (+/-13.4) PgC yr-1, which were both significantly lower than the ensemble mean from DGVMs estimates 133.7 (±4.7) PgC yr-1. We also found similar temporal trends between top-down estimates with a mean of 0.075 (±0.05) PgC yr-2, and bottom-up estimates of 0.05 (±0.05) PgC yr-2, which were 5 to 7 times smaller than the ensemble mean trend of 0.34 PgC yr-2 simulated by DGVMs. Global heterotrophic respiration showed much less agreement, ranging from top-down estimates of 42.7 (±4.0) PgC yr-1 to bottom-up estimates of 51.5 (±4.0) PgC yr-1 and a significantly larger ensemble model mean estimate of 60.8 PgC yr-1 (±1.9). The temporal trends in observation-based bottom-up estimates of heterotrophic respiration of 0.03 PgC yr-2 were five times lower than the model ensemble mean 0.15 PgC yr-2. Large regional disagreements in heterotrophic respiration estimates and simulations were evident in tropical and boreal latitudes. Therefore, improved regional and heterotrophic respiration estimates are necessary to reduce uncertainties regarding the future vulnerability of soil carbon.

Ballantyne, Ashley↗
NASA NTRSDate not supplied

Development of Modeling & Simulation Capability to Analyze Supply Chain

Supply chain has brought up many challenges to the current aerospace industry (e.g. limited sources, lack of inventory, and cyber security). The purpose of this research is to develop a web-based modeling and simulation platform for supply chains relating to Advanced Air Mobility and vertical take-off and landing vehicles (VTOLs).

Cybersecurity↗
NASA NTRS2002

Cloud-Resolving Model Simulations of LBA Convective Systems: Easterly and Westerly Regimes

The 3D Goddard Cumulus Ensemble (GCE) model was used to simulate convection that occurred during the TRMM LBA field experiment in Brazil. Convection in this region can be categorized into two different regimes. Low-level easterly flow results in moderate to high CAPE and a drier environment. Convection is more intense like that seen over continents. Low-level westerly flow results in low CAPE and a moist environment. Convection is weaker and more widespread characteristic of oceanic or monsoon-like systems. The GCE model has been used to study both regimes in order to provide cloud data sets that are representative of both environments in support of TRMM rainfall and heating algorithm development. Two different case are presented: Jan 26,1999, an easterly regime case, and Feb 23,1999, a westerly regime case. The Jan 26 case is an organized squall line and is initialized with a standard cold pool. The sensitivity to mid-level sounding moisture and wind shear will also be shown. The Feb 23 case is less-organized with only transient lines and is initialized with either warm bubbles or prescribed surface fluxes. Heating profiles, rainfall statistics and storm characteristics are compared and validated for the two cases against observations collected during the experiment.

Lang, Stephen E.↗
NASA NTRS2008

Forecasting Lightning Threat using Cloud-Resolving Model Simulations

Two new approaches are proposed and developed for making time and space dependent, quantitative short-term forecasts of lightning threat, and a blend of these approaches is devised that capitalizes on the strengths of each. The new methods are distinctive in that they are based entirely on the ice-phase hydrometeor fields generated by regional cloud-resolving numerical simulations, such as those produced by the WRF model. These methods are justified by established observational evidence linking aspects of the precipitating ice hydrometeor fields to total flash rates. The methods are straightforward and easy to implement, and offer an effective near-term alternative to the incorporation of complex and costly cloud electrification schemes into numerical models. One method is based on upward fluxes of precipitating ice hydrometeors in the mixed phase region at the-15 C level, while the second method is based on the vertically integrated amounts of ice hydrometeors in each model grid column. Each method can be calibrated by comparing domain-wide statistics of the peak values of simulated flash rate proxy fields against domain-wide peak total lightning flash rate density data from observations. Tests show that the first method is able to capture much of the temporal variability of the lightning threat, while the second method does a better job of depicting the areal coverage of the threat. Our blended solution is designed to retain most of the temporal sensitivity of the first method, while adding the improved spatial coverage of the second. Exploratory tests for selected North Alabama cases show that, because WRF can distinguish the general character of most convective events, our methods show promise as a means of generating quantitatively realistic fields of lightning threat. However, because the models tend to have more difficulty in predicting the instantaneous placement of storms, forecasts of the detailed location of the lightning threat based on single simulations can be in error. Although these model shortcomings presently limit the precision of lightning threat forecasts from individual runs of current generation models,the techniques proposed herein should continue to be applicable as newer and more accurate physically-based model versions, physical parameterizations, initialization techniques and ensembles of forecasts become available.

McCaul, Eugene W., Jr.↗
DOE OSTI2015

Macroscopic impacts of cloud and precipitation processes on maritime shallow convection as simulated by a large eddy simulation model with bin microphysics

This paper discusses impacts of cloud and precipitation processes on macrophysical properties of shallow convective clouds as simulated by a large eddy model applying warm-rain bin microphysics. Simulations with and without collision–coalescence are considered with cloud condensation nuclei (CCN) concentrations of 30, 60, 120, and 240 mg -1 . Simulations with collision–coalescence include either the standard gravitational collision kernel or a novel kernel that includes enhancements due to the small-scale cloud turbulence. Simulations with droplet collisions were discussed in Wyszogrodzki et al. (2013) focusing on the impact of the turbulent collision kernel. The current paper expands that analysis and puts model results in the context of previous studies. Despite a significant increase of the drizzle/rain with the decrease of CCN concentration, enhanced by the effects of the small-scale turbulence, impacts on the macroscopic cloud field characteristics are relatively minor. Model results show a systematic shift in the cloud-top height distributions, with an increasing contribution of deeper clouds for stronger precipitating cases. We show that this is consistent with the explanation suggested in Wyszogrodzki et al. (2013); namely, the increase of drizzle/rain leads to a more efficient condensate offloading in the upper parts of the cloud field. A second effect involves suppression of the cloud droplet evaporation near cloud edges in low-CCN simulations, as documented in previous studies (e.g., Xue and Feingold, 2006). We pose the question whether the effects of cloud turbulence on drizzle/rain formation in shallow cumuli can be corroborated by remote sensing observations, for instance, from space. Although a clear signal is extracted from model results, we argue that the answer is negative due to uncertainties caused by the temporal variability of the shallow convective cloud field, sampling and spatial resolution of the satellite data, and overall accuracy of remote sensing retrievals.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2001

Two-Dimensional Model Simulations of Interannual Variability in the Tropical Stratosphere

Meteorological data from the United Kingdom Meteorological Office (UKMO) and constituent data from the Upper Atmospheric Research Satellite (UARS) are used to construct yearly zonal mean dynamical fields for the 1990s for use in the GSFC 2-D chemistry and transport model. This allows for interannual dynamical variability to be included in the model constituent simulations. In this study, we focus on the tropical stratosphere. We find that the phase of quasi-biennial oscillation (QBO) signals in equatorial CH4, and profile and total column 03 data is resolved quite well using this empirically- based 2-D model transport framework. However. the QBO amplitudes in the model constituents are systematically underestimated relative to the observations at most levels. This deficiency is probably due in part to the limited vertical resolutions of the 2-D model and the UKMO and UARS input data sets. We find that using different heating rate calculations in the model affects the interannual and QBO amplitudes in the constituent fields, but has little impact on the phase. Sensitivity tests reveal that the QBO in transport dominates the ozone interannual variability in the lower stratosphere. with the effect of the temperature QBO being dominant in the tipper stratosphere via the strong temperature dependence of the ozone loss reaction rates. We also find that the QBO in odd nitrogen radicals, which is caused by the QBO modulated transport of NOy, plays a significant but not dominant role in determining the ozone QBO variability in the middle stratosphere. The model mean age of air is in good overall agreement with that determined from tropical lower,middle stratospheric OMS balloon observations of SF6 and CO2. The interannual variability of tile equatorial mean age in the model increases with altitude and maximizes near 40 km, with a range, of 4-5 years over the 1993-2000 time period.

Fleming, Eric L.↗
DOE OSTI2025

Advancing the Understanding of Cloud Microphysical Processes and Aerosol Indirect Effects in High-Latitude Mixed-Phase Clouds by Linking ARM Measurements with Climate Model Simulations (Final Report)

The key objectives of this project were to advance our understanding of cloud microphysical characteristics and aerosol indirect effects on mixed-phase clouds in high latitudes. To improve the representation of ice and mixed-phase clouds in Earth System Models (ESMs), we propose an integrated observation and modeling study of cloud macro- and microphysical properties, including spatial heterogeneities, mass partitioning between ice crystals and supercooled liquid water, effects of ice nucleating particles (INPs), and efficiency of secondary ice production (SIP), etc. Specifically, we took four main approaches in this project: (1) examining macro- and microphysical properties of ice and mixed-phase clouds based on in-situ and ground-based observations from multiple field campaigns funded by the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) program, including the Mixed-Phase Arctic Cloud Experiment (M-PACE), Indirect and Semi-Direct Aerosol Campaign (ISDAC), Ice Nucleating Particle Sources at Oliktok Point (INPOP), ARM West Antarctic Radiation Experiment (AWARE), Measurements of Aerosols, Radiation, and Clouds over the Southern Ocean (MARCUS), and Macquarie Island Cloud and Radiation Experiment (MICRE); (2) evaluating the DOE Energy Exascale Earth System Model (E3SM) simulations based on observations, particularly for ice and mixed-phase cloud microphysical properties; (3) examining the impacts of INPs on ice and mixed-phase clouds. Specifically, a series of comparisons were conducted using observations over the Arctic, Southern Ocean, and Antarctica, including comparisons between the lower and higher southern latitudes as well as comparisons between the northern and southern hemispheres. In addition, aerosol indirect effects from distinct sources of dust particles were examined; and (4) investigating the impacts of SIP. Ultimately, these results helped to improve cloud microphysics and aerosol-cloud interaction parameterizations in the E3SM model. Overall, the project provided improved understanding regarding various factors, including thermodynamic, dynamic, and aerosol conditions, on the micro- and macrophysical properties of ice and mixed-phase clouds in the high latitudes. Resulting analysis helped to provide an improved physical basis for refining the current cloud microphysics parameterizations related to ice and mixed-phase clouds in E3SM.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2011

Developing Flexible Discrete Event Simulation Models in an Uncertain Policy Environment

On February 1st, 2010 U.S. President Barack Obama submitted to Congress his proposed budget request for Fiscal Year 2011. This budget included significant changes to the National Aeronautics and Space Administration (NASA), including the proposed cancellation of the Constellation Program. This change proved to be controversial and Congressional approval of the program's official cancellation would take many months to complete. During this same period an end-to-end discrete event simulation (DES) model of Constellation operations was being built through the joint efforts of Productivity Apex Inc. (PAl) and Science Applications International Corporation (SAIC) teams under the guidance of NASA. The uncertainty in regards to the Constellation program presented a major challenge to the DES team, as to: continue the development of this program-of-record simulation, while at the same time remain prepared for possible changes to the program. This required the team to rethink how it would develop it's model and make it flexible enough to support possible future vehicles while at the same time be specific enough to support the program-of-record. This challenge was compounded by the fact that this model was being developed through the traditional DES process-orientation which lacked the flexibility of object-oriented approaches. The team met this challenge through significant pre-planning that led to the "modularization" of the model's structure by identifying what was generic, finding natural logic break points, and the standardization of interlogic numbering system. The outcome of this work resulted in a model that not only was ready to be easily modified to support any future rocket programs, but also a model that was extremely structured and organized in a way that facilitated rapid verification. This paper discusses in detail the process the team followed to build this model and the many advantages this method provides builders of traditional process-oriented discrete event simulations.

Miranda, David J.↗
DOE OSTI2016

Using GRACE and climate model simulations to predict mass loss of Alaskan glaciers through 2100

Glaciers in Alaska are currently losing mass at a rate of ~–50 Gt a –1 , one of the largest ice loss rates of any regional collection of mountain glaciers on Earth. Existing projections of Alaska's future sea-level contributions tend to be divergent and are not tied directly to regional observations. Here we develop a simple, regional observation-based projection of Alaska's future sea-level contribution. We compute a time series of recent Alaska glacier mass variability using monthly GRACE gravity fields from August 2002 through December 2014. We also construct a three-parameter model of Alaska glacier mass variability based on monthly ERA-Interim snowfall and temperature fields. When these three model parameters are fitted to the GRACE time series, the model explains 94% of the variance of the GRACE data. Using these parameter values, we then apply the model to simulated fields of monthly temperature and snowfall from the Community Earth System Model, to obtain predictions of mass variations through 2100. Here, we conclude that mass loss rates may increase between –80 and –110 Gt a –1 by 2100, with a total sea-level rise contribution of 19 ± 4 mm during the 21st century.

54 ENVIRONMENTAL SCIENCES↗