Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “snow contribution”

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 163 records · Page 9

Surface Deformation and Gravity Changes from Surface and Internal Loads

Air and space borne remote sensing have made it possible to monitor the mass and energy transport at various scales within the cryosphere-hydrosphere-atmosphere system. The recent surface mass balance (the rate of net gain of snow and ice at a geographic point) map for the Antarctic ice sheet is constructed by interpolating sparse in situ observations (about 1,800 points) with empirically calibrated satellite data of passive back emission of microwaves. The digital elevation model obtained from satellite radar altimetry is used to improve the delineation of the ice flow drainage basins. As important as these results are, the uncertainty remains up to about 2 mm/yr of eustatic sea level change with the net imbalance. In other words, we are still unable to determine even the sign of the contribution of the Antarctic ice sheet to contemporary sea level change. The problem is more likely with the discharge rather than accumulation.

Hager, Bradford H.↗

Application of Satellite Gravimetry for Water Resource Vulnerability Assessment

The force of Earth's gravity field varies in proportion to the amount of mass near the surface. Spatial and temporal variations in the gravity field can be measured via their effects on the orbits of satellites. The Gravity Recovery and Climate Experiment (GRACE) is the first satellite mission dedicated to monitoring temporal variations in the gravity field. The monthly gravity anomaly maps that have been delivered by GRACE since 2002 are being used to infer changes in terrestrial water storage (the sum of groundwater, soil moisture, surface waters, and snow and ice), which are the primary source of gravity variability on monthly to decadal timescales after atmospheric and oceanic circulation effects have been removed. Other remote sensing techniques are unable to detect water below the first few centimeters of the land surface. Conventional ground based techniques can be used to monitor terrestrial water storage, but groundwater, soil moisture, and snow observation networks are sparse in most of the world, and the countries that do collect such data rarely are willing to share them. Thus GRACE is unique in its ability to provide global data on variations in the availability of fresh water, which is both vital to life on land and vulnerable to climate variability and mismanagement. This chapter describes the unique and challenging aspects of GRACE terrestrial water storage data, examples of how the data have been used for research and applications related to fresh water vulnerability and change, and prospects for continued contributions of satellite gravimetry to water resources science and policy.

Rodell, Matthew↗

SMAP Data Assimilation at the GMAO

The NASA Soil Moisture Active Passive (SMAP) mission has been providing L-band (1.4 GHz) passive microwave brightness temperature (Tb) observations since April 2015. These observations are sensitive to surface(0-5 cm) soil moisture. Several of the key applications targeted by SMAP, however, require knowledge of deeper-layer, root zone (0-100 cm) soil moisture, which is not directly measured by SMAP. The NASA Global Modeling and Assimilation Office (GMAO) contributes to SMAP by providing Level 4 data, including the Level 4 Surface and Root Zone Soil Moisture(L4_SM) product, which is based on the assimilation of SMAP Tb observations in the ensemble-based NASA GEOS-5 land surface data assimilation system. The L4_SM product offers global data every three hours at 9 km resolution, thereby interpolating and extrapolating the coarser- scale (40 km) SMAP observations in time and in space (both horizontally and vertically). Since October 31, 2015, beta-version L4_SM data have been available to the public from the National Snow and Ice Data Center for the period March 31, 2015, to near present, with a mean latency of approx. 2.5 days.

Land Surface Model↗

National Climate Assessment - Land Data Assimilation System (NCA-LDAS) Data and Services at NASA GES DISC

The National Climate Assessment-Land Data Assimilation System (NCA-LDAS) is an Integrated Terrestrial Water Analysis, and is one of NASAs contributions to the NCA of the United States. The NCA-LDAS has undergone extensive development, including multi-variate assimilation of remotely-sensed water states and anomalies as well as evaluation and verification studies, led by the Goddard Space Flight Centers Hydrological Sciences Laboratory (HSL). The resulting NCA-LDAS data have recently been released to the general public and include those from the Noah land-surface model (LSM) version 3.3 (Noah-3.3) and the Catchment LSM version Fortuna-2.5 (CLSM-F2.5). Standard LSM output variables including soil moistures temperatures, surface fluxes, snow cover depth, groundwater, and runoff are provided, as well as streamflow using a river routing system. The NCA-LDAS data are archived at and distributed by the NASA Goddard Earth Sciences Data and Information Services Center (GES DISC). The data can be accessed via HTTP, OPeNDAP, Mirador search and download, and NASA Earth data Search. To further facilitate access and use, the NCA-LDAS data are integrated into the NASA Giovanni, for quick visualization and analysis, and into the Data Rods system, for retrieval of time series of long time periods. The temporal and spatial resolutions of the NCA-LDAS data are, respectively, daily-averages and 0.125x0.125 degree, covering North America (25N 53N; 125W 67W) and the period January 1979 to December 2015. The data files are in self-describing, machine-independent, CF-compliant netCDF-4 format.

data assimiliation↗

A search for interstellar scandium

The resonance line of Sc II at 3642.785 A has been sought in the interstellar spectrum of Zeta Oph, resulting in an upper limit N(Sc II) not greater than 8.6 billion/sq cm and a depletion of this element at least by a factor of 200, if Sc III contributes a negligible column density. The present upper limit may be somewhat below the level expected under the grain-formation model for interstellar depletions, and is consistent with the abundance expected if the depletions are due to gas-grain collisions. The new limit on this element helps to support the hypothesis that species within a specific range of ionization potential are preferentially depleted.

Snow, T. P., Jr.↗

Comparison of Satellite-Derived and In-Situ Observations of Ice and Snow Surface Temperatures over Greenland

The most practical way to get a spatially broad and continuous measurements of the surface temperature in the data-sparse cryosphere is by satellite remote sensing. The uncertainties in satellite-derived LSTs must be understood to develop internally-consistent decade-scale land-surface temperature (LST) records needed for climate studies. In this work we assess satellite-derived "clear-sky" LST products from the Moderate Resolution Imaging Spectroradiometer (MODIS) and the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER), and LSTs derived from the Enhanced Thematic Mapper Plus (ETM+) over snow and ice on Greenland. When possible, we compare satellite-derived LSTs with in-situ air-temperature observations from Greenland Climate Network (GC-Net) automatic-weather stations (AWS). We find that MODIS, ASTER and ETM+ provide reliable and consistent LSTs under clear-sky conditions and relatively-flat terrain over snow and ice targets over a range of temperatures from -40 to 0 C. The satellite-derived LSTs agree within a relative RMS uncertainty of approx.0.5 C. The good agreement among the LSTs derived from the various satellite instruments is especially notable since different spectral channels and different retrieval algorithms are used to calculate LST from the raw satellite data. The AWS record in-situ data at a "point" while the satellite instruments record data over an area varying in size from: 57 X 57 m (ETM+), 90 X 90 m (ASTER), or to 1 X 1 km (MODIS). Surface topography and other factors contribute to variability of LST within a pixel, thus the AWS measurements may not be representative of the LST of the pixel. Without more information on the local spatial patterns of LST, the AWS LST cannot be considered valid ground truth for the satellite measurements, with RMS uncertainty approx.2 C. Despite the relatively large AWS-derived uncertainty, we find LST data are characterized by high accuracy but have uncertain absolute precision.

Hall, Dorothy K.↗

Reanalyses and Essential Climate Variables

Reanalyses are a potentially powerful climate data collection driven by observations but also subjected to model bias. Additionally, reanalyses can produce and use essential climate variables in a consistent method. For example, snow cover and soil moisture (among other variables) will eventually be assimilated into the reanalyses, but also provide crucial validation data. Sea surface temperature can be prescribed or assimilated in a coupled reanalysis. The strength of reanalysis lies in the ancillary data that is produced from the modeling components but not routinely observed thereby providing more complete Earth system information. The weakness in this concept is that the model derived data can be affected by model bias and may also change relative to the available observing system. Here, we will review the status of existing reanalyses and the ECVs being considered for the workshop. Purpose of Michael Bosilovich's contribution to the workshop: Michael Bosilovich will represent US reanalysis community in this international discussion of Essential Climate Variables (ECVs) and the relative nature of reanalyses to ECVs.

Bosilovich, Michael↗

Copernicus studies of interstellar material in the Perseus II complex. III - The line of sight to Zeta Persei

Ultraviolet spectrophotometric data obtained with Copernicus are used to analyze the distribution, composition, density, temperature, and kinematics of the interstellar material along the line of sight to Zeta Persei. The far-UV extinction curve for the star is evaluated along with the kinematics of the interstellar gas, observations of atomic and molecular hydrogen, curves of growth for neutral and ionized species, atomic abundances and depletions, ionization equilibria, and observations of CO and OH lines. The results show that there are apparently three clouds along the line of sight to Zeta Persei: a main cloud at approximately +13 km/s which contains most of the material and forms all the neutral and molecular lines as well as most of the ionic lines, a second component at +22 km/s which must contribute to the strong UV lines of most ions, and a third component at roughly +2 km/s which gives rise to a strong Si III line at 1206 A. It is also found that the UV extinction curve has a somewhat steep far-UV rise, indicating the presence of a substantial number of small grains, and that about 30% of the hydrogen nuclei over the entire line of sight are in molecular form.

Snow, T. P., Jr.↗

Diagnosing Changes In Earth's Top-Of-Atmosphere Radiation Budget Since 2000

Earth’s radiation budget (ERB) is a key indicator of the climate system as it describes the delicate balance between how much of the sun’s energy reaching Earth is absorbed and how much thermal infrared radiation is emitted back to space. It is now well established that there is a positive net imbalance in ERB of 0.5-1 Wm-2. This imbalance, referred to as Earth’s Energy Imbalance (EEI), is shaped by natural and anthropogenic climate forcings and the climate system’s response to the forcing.When averaged over a long period, a positive EEI leads to an increase in surface temperature, sea levelrise, ocean heat content, melting of snow and sea-ice, and an acceleration of the hydrological cycle.EEI is also influenced by internal climate variability at interannual (e.g., El Niño/Southern Oscillation(ENSO)) and longer timescales (e.g., Pacific Decadal Oscillation (PDO)). Here we use CERES data toshow that since 2000, EEI has been increasing at a rate of 0.04+-0.03 Wm-2 per yr, corresponding to 0.8 Wm-1 over the 20-year CERES record. To understand the underlying causes for this increase, we analyze the CERES data in the context of internal climate variability indices and regional pattern changes. We also perform a partial radiative perturbation analysis that decomposes the TOA SW and LW radiation variations in terms of contributions by clouds, aerosols, surface albedo, temperature,water vapor and other gases (e.g., ozone, carbon dioxide, methane, etc.).

Norman G Loeb↗

An Update on Oceanic Precipitation Rate and its Zonal Distribution in Light of Advanced Observations from Space

This study contributes to the estimation of the global mean and zonal distribution of oceanic precipitation rate using complementary information from advanced precipitation measuring sensors and provides an independent reference to assess current precipitation products. Precipitation estimates from the Tropical Rainfall Measuring Mission (TRMM) precipitation radar (PR) and CloudSat cloud profiling radar (CPR) were merged, as the two complementary sensors yield an unprecedented range of sensitivity to quantify rainfall from drizzle through the most intense rates. At higher latitudes, where TRMM PR does not exist, precipitation estimates from Aqua's Advanced Microwave Scanning Radiometer for Earth Observing System (AMSR-E) complemented CloudSat CPR to capture intense precipitation rates. The high sensitivity of CPR allows estimation of snow rate, an important type of precipitation at high latitudes, not directly observed in current merged precipitation products. Using the merged precipitation estimate from the CloudSat, TRMM, and Aqua platforms (this estimate is abbreviated to MCTA), the authors' estimate for 3-yr (2007-09) nearglobal (80degS-80degN) oceanic mean precipitation rate is approx. 2.94mm/day. This new estimate of mean global ocean precipitation is about 9% higher than that of the corresponding Climate Prediction Center (CPC) Merged Analysis of Precipitation (CMAP) value (2.68mm/day) and about 4% higher than that of the Global Precipitation Climatology Project (GPCP; 2.82mm/day). Furthermore, MCTA suggests distinct differences in the zonal distribution of precipitation rate from that depicted in GPCPand CMAP, especially in the Southern Hemisphere.

profiling radar↗

Quantifying the Early Snowmelt Event of 2015 in the Cascade Mountains, USA by Developing and Validating MODIS-Based Snowmelt Timing Maps

Spring snowmelt serves as the major hydrological contribution to many watersheds of the US West. Since the 1970s the conterminous western USA has seen an earlier arrival of spring snowmelt. The extremely low snowpack and early melt of 2015 in the Cascade Mountains may be a harbinger of winters to come, underscoring the interest in advancements in spring snowmelt monitoring. Target-of-opportunity and point measurements of snowmelt using meteorological stations or stream gauges are common sources of these data, however, there have been few attempts to identify snowmelt timing using remote sensing. In this study, we describe the creation of snowmelt timing maps (STMs) which identify the day of year that each pixel of a remotely sensed image transitions from "snow-covered" to "no snow" during the spring melt season, controlling for cloud coverage and ephemeral spring snow storms. Derived from the 500 meter MODerate-resolution Imaging Spectroradiometer (MODIS) standard snow map, MOD10A2, this new dataset provides annual maps of snowmelt timing, with corresponding maps of cloud interference and interannual variability in snow coverage from 2001-2015. We first show that the STMs agree strongly with in-situ snow telemetry (SNOTEL) meteorological station measurements in terms of snowmelt timing. We then use the STMs to investigate the early snowmelt event of 2015 in the Cascade Mountains, USA, highlighting the protected areas of Mt. Rainier, Crater Lake, and Lassen Volcanic National Parks. In 2015 the Cascade Mountains experienced snowmelt 41 days earlier than the 2001-2015 average, with 25 percent of its land area melting more than 65 days earlier than average. The upper elevations of the Cascade Mountains experienced the greatest snowmelt anomaly. Our results are relevant to land managers and biologists as they plan adaptation strategies for mitigating the effects of climate change throughout temperate mountains.

snowmelt timing↗

Global Precipitation Measurement (GPM) Mission Applications: Activities, Challenges, and Vision

Global Precipitation Measurement (GPM) is an international satellite mission to provide nextgeneration observations of rain and snow worldwide every three hours. NASA and the Japan Aerospace Exploration Agency (JAXA) will launch a "Core" satellite carrying advanced instruments that will set a new standard for precipitation measurements from space. The data they provide will be used to unify precipitation measurements made by an international network of partner satellites to quantify when, where, and how much it rains or snows around the world. The GPM mission will help advance our understanding of Earth's water and energy cycles, improve the forecasting of extreme events that cause natural disasters, and extend current capabilities of using satellite precipitation information to directly benefit society. Building upon the successful legacy of the Tropical Rainfall Measuring Mission (TRMM), GPM's next-generation global precipitation data will lead to scientific advances and societal benefits within a range of hydrologic fields including natural hazards, ecology, public health and water resources. This talk will highlight some examples from TRMM's IS-year history within these applications areas as well as discuss some existing challenges and present a look forward for GPM's contribution to applications in hydrology.

Kirschbaum, Dalia↗

Climate Variability in the Amundsen and Bellingshausen Seas

Satellite data reveal a 20% decline in sea ice extent in the Amundsen and Bellingshausen Seas in the two decades following 1973. This change is negatively correlated with surface air temperatures on the west side of the Antarctic Peninsula, which have increased -0.5 C / decade- since the mid-1940s, The recession was strongest during summer, when monthly average minima in 1991-92 removed much of the incipient multiyear ice over the continental shelf. This would have lowered the regional-mean ice thickness, impacting snow ice formation, brine production, and vertical heat flux. The northern ice edge contracted by 1 deg of latitude in all seasons from 1973-79 to 1987-93, returning toward mean conditions in 1993-95. The decline included multiyear cycles of several years in length, superimposed on high interannual variability. A review of atmospheric forcing shows winds consistent with mean and extreme ice extents, and suggests links to larger-scale circulation changes in the South Pacific. Historical ocean measurements are sparse in this sector, but mixed-layer depths and upper pycnoclines beneath the sea ice resemble those in the Weddell Sea. Weaker surface currents or changes in the upwelling of Circumpolar Deep Water on the continental shelf could have contributed to the anomaly persistence.

Jacobs, Stanley S.↗

An Evaluation of A Physics-Based Firn Model and A Semi-Empirical Firn Model Across the Greenland Ice Sheet (1980–2020)

The Greenland Ice Sheet's (GrIS) firn layer buffers the ice sheet's contribution to sea level rise by storing meltwater in its pore space. However, available pore space and meltwater retention capability is lost due to ablation of the firn layer and refreezing of meltwater as near-surface ice slabs in the firn. Understanding how firn properties respond to climate is important for constraining the GrIS's future contribution to sea level rise in a warming climate. Observations of firn density provide detailed information about firn properties, but they are spatially and temporally limited. Here we use two firn models, the physics-based SNOWPACK model and the Community Firn Model configured with a semi-empirical densification equation (CFM-GSFC), to quantify firn properties across the GrIS from 1980 through 2020. We use an identical forcing (Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2) atmospheric reanalysis) for SNOWPACK and the CFM-GSFC in order to isolate firn model differences. To evaluate the models, we compare simulated firn properties, including firn air content (FAC), to measurements from the Surface Mass Balance and Snow on Sea Ice Working Group (SUMup) dataset of snow and firn density. Both models perform well (mean absolute percentage errors of 14 % in SNOWPACK and 16 % in the CFM-GSFC), though their performance is hindered by the spatial resolution of the atmospheric forcing. In the ice-sheet-wide simulations, the 1980–1995 average spatially integrated FAC (i.e., air volume in the firn) for the upper 100 m is 34 645 km 3 from SNOWPACK and 28 581 km 3 from the CFM-GSFC. The discrepancy in the magnitude of the modeled FAC stems from differences in densification with depth and variations in the sensitivity of the models to atmospheric forcing. In more recent years (2005–2020), both models simulate substantial depletion of pore space. During this period, the spatially integrated FAC across the entire GrIS decreases by 3.2 % (−66.6 km 3 yr −1 ) in SNOWPACK and 1.5 % (−17.4 km 3 yr −1 ) in the CFM-GSFC. These differing magnitudes demonstrate how model differences propagate throughout the FAC record. Over the full modeled record (1980–2020), SNOWPACK simulates a loss of pore space equivalent to 3 mm of sea level rise buffering, while the CFM-GSFC simulates a loss of 1 mm. The greatest depletion in FAC is along the margins and especially along the western margin where observations and models show the formation of near-surface, low-permeability ice slabs that may inhibit meltwater storage.

Megan Thompson-Munson↗

Impacts of aerosol size-dependent below-cloud scavenging on tropospheric aerosol in the NASA GEOS model

Cloud scavenging is the major sink for a suite of tropospheric aerosols, and thus largely affects the global distribution and lifetime of aerosols as well as the formation of clouds. Simplified model parameterizations of cloud scavenging of aerosols substantially contribute to large uncertainties in the simulated aerosol burdens, spatial distributions, lifetimes, as well as aerosol radiative forcing. Here we evaluate the wet scavenging scheme in the current NASA GEOS model, develop more physically-based parameterizations of below-cloud scavenging, and assess the impacts of the new scavenging parameterizations on the GEOS simulations of aerosols. The current below-cloud scavenging (BCS; washout) scheme uses a first-order removal with a constant scavenging coefficient. We replaced the current scheme with a size-dependent parameterization of Croft et al.(Atmos. Chem. Phys., 2009; hereafter referred to as C09). The C09 BCS parameterizations calculate aerosol below-cloud scavenging coefficients based on a size-dependent collision efficiency between aerosol particles and raindrops (or snow crystals). A look-up table of compiled collision efficiencies provides values for sixty aerosol sizes and nine rain flux rates. Snow-aerosol collision efficiencies are a function of the sixty aerosol sizes only. We assume a log-normal size distribution for aerosols (bulk) in GOCART. The collision efficiencies are then calculated online using a bilinear interpolation from the look-up table for the given aerosol sizes and precipitation flux rates in the model.

aerosols↗

Atmospheric River Precipitation Contributed to Rapid Increases in Surface Height of the West Antarctic Ice Sheet in 2019

Estimating the relative contributions of the atmospheric and dynamic components of ice-sheet mass balance is critical for improving projections of future sea level rise. Existing estimates of changes in Antarctic ice-sheet height, which can be used to infer changes in mass, are only accurate at multiyear time scales. However, NASA's Ice, Cloud, and land Elevation Satellite-2 (ICESat-2) laser altimetry mission now allows us to accurately measure changes in ice-sheet height at subannual time scales. Here, we use ICESat-2 data to estimate height changes over Antarctica between April 2019 and June 2020. These data show widespread increases in surface height over West Antarctica during the 2019 austral winter. Using climate reanalysis data, we show that 41% of increases in height during winter were from snow accumulation via extreme precipitation events—63% of these events were associated with landfalling atmospheric rivers (ARs) which occurred only 5.1% of the time.

Susheel Adusumilli↗

Advances in Land Data Assimilation at the NASA Goddard Space Flight Center

Research in land surface data assimilation has grown rapidly over the last decade. In this presentation we provide a brief overview of key research contributions by the NASA Goddard Space Flight Center (GSFC). The GSFC contributions to land assimilation primarily include the continued development and application of the Land Information System (US) and the ensemble Kalman filter (EnKF). In particular, we have developed a method to generate perturbation fields that are correlated in space, time, and across variables and that permit the flexible modeling of errors in land surface models and observations, along with an adaptive filtering approach that estimates observation and model error input parameters. A percentile-based scaling method that addresses soil moisture biases in model and observational estimates opened the path to the successful application of land data assimilation to satellite retrievals of surface soil moisture. Assimilation of AMSR-E surface soil moisture retrievals into the NASA Catchment model provided superior surface and root zone assimilation products (when validated against in situ measurements and compared to the model estimates or satellite observations alone). The multi-model capabilities of US were used to investigate the role of subsurface physics in the assimilation of surface soil moisture observations. Results indicate that the potential of surface soil moisture assimilation to improve root zone information is higher when the surface to root zone coupling is stronger. Building on this experience, GSFC leads the development of the Level 4 Surface and Root-Zone Soil Moisture (L4_SM) product for the planned NASA Soil-Moisture-Active-Passive (SMAP) mission. A key milestone was the design and execution of an Observing System Simulation Experiment that quantified the contribution of soil moisture retrievals to land data assimilation products as a function of retrieval and land model skill and yielded an estimate of the error budget for the SMAP L4_SM product. Terrestrial water storage observations from GRACE satellite system were also successfully assimilated into the NASA Catchment model and provided improved estimates of groundwater variability when compared to the model estimates alone. Moreover, satellite-based land surface temperature (LST) observations from the ISCCP archive were assimilated using a bias estimation module that was specifically designed for LST assimilation. As with soil moisture, LST assimilation provides modest yet statistically significant improvements when compared to the model or satellite observations alone. To achieve the improvement, however, the LST assimilation algorithm must be adapted to the specific formulation of LST in the land model. An improved method for the assimilation of snow cover observations was also developed. Finally, the coupling of LIS to the mesoscale Weather Research and Forecasting (WRF) model enabled investigations into how the sensitivity of land-atmosphere interactions to the specific choice of planetary boundary layer scheme and land surface model varies across surface moisture regimes, and how it can be quantified and evaluated against observations. The on-going development and integration of land assimilation modules into the Land Information System will enable the use of GSFC software with a variety of land models and make it accessible to the research community.

Reichle, Rolf↗

Critical Load Exceedances for North America and Europe using an Ensemble of Models and an Investigation of Causes for Environmental Impact Estimate Variability: An AQMEII4 Study

Exceedances of critical loads for deposition of sulphur (S) and nitrogen (N) to different ecosystems were estimated using European and North American ensembles of air quality models, under Phase 4 of the Air Quality Model Evaluation International Initiative (AQMEII4), to identify where risk of ecosystem harm is expected to occur based on model deposition estimates. The ensembles were driven by common emissions and lateral boundary condition inputs. Model output was regridded to common North American and Europe 0.125° resolution domains, which were then used to calculate critical load exceedances. New, targeted deposition diagnostics implemented in AQMEII4 allowed an unprecedented level of post-simulation analysis to be carried out and facilitated the identification of specific causes of model-to-model variability in critical load exceedance estimates. New datasets for North American critical loads for acidity for forest soil water and aquatic ecosystems were combined with the ensemble deposition predictions to show a substantial decrease in the area and number of locations in exceedance between 2010 and 2016 (forest soils: 13.2 % to 6.1 %; aquatic ecosystems: 21.2 % to 11.4 %). All models agreed in the direction of the ensemble exceedance change between 2010 and 2016. The North American ensemble also predicted a decrease in both severity and total area in exceedance between the years 2010 and 2016 for eutrophication-impacted ecosystems in the USA (sensitive epiphytic lichen: 81.5 % to 75.8 %). The exceedances for herbaceous community richness also decreased between 2010 and 2016, from 13.9 % to 3.9 %. The uncertainty associated with the North American eutrophication results is high; there were sharp differences between the models in both predictions of total N deposition and the change in N deposition, and hence in the predicted eutrophication exceedances between the two years. The European ensemble was used to predict relatively static exceedances of critical loads with respect to acidification (4.48 % to 4.32 % from 2009 to 2010) while eutrophication exceedance increased slightly (60.2 % to 62.2 %). While most models showed the same changes in critical load exceedances as the ensemble between the two years, the spatial extent and magnitude of exceedances varied significantly between the models. The reasons for this variation were examined in detail by first ranking the relative contribution of different sources of sulphur and nitrogen deposition in terms of deposited mass and model-to-model variability in that deposited mass, followed by their analysis using AQMEII4 diagnostics, along with evaluation of the most recent literature. All models in both the North American and European ensembles had net annual negative biases with respect to observed wet deposition of sulphate, nitrate and ammonium. Diagnostics and recent literature suggest that this bias may stem from insufficient cloud scavenging of aerosols and gases, and may be improved through the incorporation of multiphase hydrometeor scavenging within the modelling frameworks. The inability of North American models to predict the timing of the seasonal peak in wet ammonium ion deposition (observed maximum was in April, while all models predicted a June maximum) may also relate to the need for multiphase hydrometeor scavenging (absence of snow scavenging in all models employed here). High variability in the relative importance of particulate sulphate, nitrate and ammonium deposition fluxes between models was linked to the use of updated particle dry deposition parameterizations in some models. However, recent literature and further development of some of the models within the ensemble suggests these particulate biases may also be ameliorated via the incorporation of multiphase hydrometeor scavenging. Annual sulphur and nitrogen deposition prediction variability was linked to SO 2 and HNO 3 dry deposition parameterizations, and diagnostic analysis showed that the cuticle and soil deposition pathways dominate the deposition mass flux of these species. Further work improving parameterizations for these deposition pathways should reduce variability in model acidifying gas deposition estimates. The absence of base cation chemistry in some models was shown to be a major factor in positive biases in fine mode particulate ammonium and particle nitrate concentrations. Models employing ammonia bidirectional fluxes had both the largest and the smallest magnitude biases, depending on the model and bidirectional flux algorithm employed. A careful analysis of bidirectional flux models suggests that those with poor NH 3 performance may underestimate the extent of NH 3 emissions fluxes from forested areas. Based on these results, an increased process-research focus is therefore recommended for the following model processes and on observations which may assist in model evaluation and improvement: multiphase hydrometeor scavenging combined with updated particle dry deposition, cuticle and soil deposition pathway algorithms for acidifying gases, base cation chemistry and emissions, and NH 3 bidirectional fluxes. Comparisons with satellite observations suggest that oceanic NH 3 emissions sources should be included in regional chemical transport models. The choice of land use database employed within any given model was shown to significantly influence deposition totals in several instances, and employing a common land use database across chemical transport models and critical load calculations is recommended for future work.

critical loads↗