Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ocean model comparisons”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

The Gulf Stream front - A comparison between Seasat altimeter observations and theory

A quantitative comparison is made between a model which gives a universal shape and width of the sea surface height rise anomaly across the Gulf Stream and six Gulf Stream altimeter profiles representing a range of conditions. The closeness of the fit obtained reinforces the validity of both the model and the role of satellite altimetry in ocean dynamics. Further comparisons of the data with the theories of Stommel (1966) and Charney (1955) show that the former is incomplete, accepting only one altimeter-derived parameter.

Kao, T. W.↗
NASA NTRS2023

The Sensitivity of the Equatorial Pacific ODZ to Particulate Organic Matter Remineralization in A Climate Model Under Pre-Industrial Conditions

Marine oxygen plays a fundamental role in regulating the transfer of organic carbon and nutrients to their dissolved inorganic forms, serving as the terminal electron acceptor for heterotrophic respiration. Oxygen can become limiting to these processes in coastal and open-ocean oxygen deficient zones (ODZs). The maintenance of ODZs depends on the balance between physical processes such as ventilation and biogeochemical processes such as remineralization. These two processes act in opposing directions on ODZs, with ventilation responsible for transporting oxygen from the surface, where it is near saturation with the partial pressure of O2 in the atmosphere, to deeper oxygen-depleted layers, and remineralization acting to consume oxygen through biogeochemical processes such as microbial degradation throughout the water column. Remineralization is represented in all CMIP6 models, but the actual parameterizations, as well as their magnitude, are widely varying. In this study, we examine the sensitivity to remineralization of the equatorial Pacific ODZ O2 ol/kg) in a model; the NASA GISS Ocean Biogeochemical Model (NOBM-G) embedded in the NASA-GISS coupled atmosphere-ocean model. We find that increasing the remineralization rate shoals the ODZ onset depth (i.e. the regionally averaged upper bound of the ODZ), and decreases ODZ volume and regional-averaged thickness. Changes in biological consumption and vertical convergence of oxygen are identified as the primary contributors to changes in ODZ onset depth. ODZ thickness is mainly influenced by the shoaling of the bottom boundary of the ODZ, which decreases with maximum remineralization rate due to decreasing oxygen consumption and increasing horizontal oxygen convergence in deep waters. On the other hand, vertically-averaged ODZ area has a more complex, non-monotonic relationship with maximum remineralization rate. While these findings suggest an important role for remineralization in determining the shape of the ODZ in our model, the relative importance of remineralization vs. other physical parameterizations remains to be established. While our results reflect the structure of our ocean biogeochemical model, the relationship of remineralization to ODZ characteristics in other models should be examined to better inform model inter-comparisons.

Oxygen-deficient zone↗
NASA NTRS1981

The effect of zonal gradients of sea surface temperature on the Indian Ocean winter monsoon

Several global climate simulations by the 7-layer, 8 x 10 GISS climate model were designed to test the contributions of various surface boundary conditions to the global climate. The model was run with the sun fixed at a perpetual January. In a comparison of run #5, in which realistic January surface boundary conditions were used, with run #4, which was the same except that a zonally symmetric climatological January sea surface temperature (SST) field was used, one of the results was that run #5 provided a better simulation of the Indian Ocean monsoon. A further comparison of the wind fields over the Indian Ocean that were generated by these two model runs is presented.

Cohen, C.↗
NASA NTRS1983

A baroclinic quasigeostrophic open ocean model

A baroclinic quasigeostrophic open ocean model is presented, calibrated by a series of test problems, and demonstrated to be feasible and efficient for application to realistic mid-oceanic mesoscale eddy flow regimes. Two methods of treating the depth dependence of the flow, a finite difference method and a collocation method, are tested and intercompared. Sample Rossby wave calculations with and without advection are performed with constant stratification and two levels of nonlinearity, one weaker than and one typical of real ocean flows. Using exact analytical solutions for comparison, the accuracy and efficiency of the model is tabulated as a function of the computational parameters and stability limits set; typically, errors were controlled between 1 percent and 10 percent RMS after two wave periods. Further Rossby wave tests with realistic stratification and wave parameters chosen to mimic real ocean conditions were performed to determine computational parameters for use with real and simulated data. Finally, a prototype calculation with quasiturbulent simulated data was performed successfully, which demonstrates the practicality of the model for scientific use.

Miller, R. N.↗
NASA NTRS2007

Evaluating the Capacity of Global CO2 Flux and Atmospheric Transport Models to Incorporate New Satellite Observations

As we enter the new era of satellite remote sensing for CO2 and other carbon cyclerelated quantities, advanced modeling and analysis capabilities are required to fully capitalize on the new observations. Model estimates of CO2 surface flux and atmospheric transport are required for initial constraints on inverse analyses, to connect atmospheric observations to the location of surface sources and sinks, and ultimately for future projections of carbon-climate interactions. For application to current, planned, and future remotely sensed CO2 data, it is desirable that these models are accurate and unbiased at time scales from less than daily to multi-annual and at spatial scales from several kilometers or finer to global. Here we focus on simulated CO2 fluxes from terrestrial vegetation and atmospheric transport mutually constrained by analyzed meteorological fields from the Goddard Modeling and Assimilation Office for the period 1998 through 2006. Use of assimilated meteorological data enables direct model comparison to observations across a wide range of scales of variability. The biospheric fluxes are produced by the CASA model at lxi degrees on a monthly mean basis, modulated hourly with analyzed temperature and sunlight. Both physiological and biomass burning fluxes are derived using satellite observations of vegetation, burned area (as in GFED-2), and analyzed meteorology. For the purposes of comparison to CO2 data, fossil fuel and ocean fluxes are also included in the transport simulations. In this presentation we evaluate the model's ability to simulate CO2 flux and mixing ratio variability in comparison to in situ observations at sites in Northern mid latitudes and the continental tropics. The influence of key process representations is inferred. We find that the model can resolve much of the hourly to synoptic variability in the observations, although there are limits imposed by vertical resolution of boundary layer processes. The seasonal cycle and its interannual variations generally respond adequately, but discrepancies in the tropics suggest the need for a refinement of the soil moisture dependence of the respiration flux in CASA. Examples and inferences for interpretation of satellite data will be discussed. In general, the fidelity of these simulations leads us to anticipate incorporation of real-time, highly resolved remote sensing and other observations into quantitative analyses that will reduce uncertainty in the terrestrial CO2 sink and revolutionize our understanding of the key processes controlling atmospheric CO2 and its evolution with time.

Kawa, S. R.↗
NASA NTRS2010

Comparing Global Atmospheric CO2 Flux and Transport Models with Remote Sensing (and Other) Observations

We report recent progress derived from comparison of global CO2 flux and transport models with new remote sensing and other sources of CO2 data including those from satellite. The overall objective of this activity is to improve the process models that represent our understanding of the workings of the atmospheric carbon cycle. Model estimates of CO2 surface flux and atmospheric transport processes are required for initial constraints on inverse analyses, to connect atmospheric observations to the location of surface sources and sinks, to provide the basic framework for carbon data assimilation, and ultimately for future projections of carbon-climate interactions. Models can also be used to test consistency within and between CO2 data sets under varying geophysical states. Here we focus on simulated CO2 fluxes from terrestrial vegetation and atmospheric transport mutually constrained by analyzed meteorological fields from the Goddard Modeling and Assimilation Office for the period 2000 through 2009. Use of assimilated meteorological data enables direct model comparison to observations across a wide range of scales of variability. The biospheric fluxes are produced by the CASA model at 1x1 degrees on a monthly mean basis, modulated hourly with analyzed temperature and sunlight. Both physiological and biomass burning fluxes are derived using satellite observations of vegetation, burned area (as in GFED-3), and analyzed meteorology. For the purposes of comparison to CO2 data, fossil fuel and ocean fluxes are also included in the transport simulations. In this presentation we evaluate the model's ability to simulate CO2 flux and mixing ratio variability in comparison to remote sensing observations from TCCON, GOSAT, and AIRS as well as relevant in situ observations. Examples of the influence of key process representations are shown from both forward and inverse model comparisons. We find that the model can resolve much of the synoptic, seasonal, and interannual variability in the observations, although reasons for persistent discrepancies in northern hemisphere vegetation uptake are examined. At this time, we do not find any serious shortcomings in the model transport representation, but this is still the subject of close scrutiny. In general, the fidelity of these simulations leads us to anticipate incorporation of real-time, highly resolved remote sensing and other observations into quantitative analyses that will reduce uncertainty in CO2 fluxes and revolutionize our understanding of the key processes controlling atmospheric CO2 and its evolution with time.

Kawa, S. R.↗
NASA NTRS1993

Comparison of remote sensing measurements with a two-scale polarimetric emission and scattering model for sea surfaces

Recently, it has been observed that the brightness temperatures of sea surfaces correlate with the azimuth angle of the ocean wind vector, indicating that ocean wind direction can not only be retrieved from the microwave backscatter, but also from the brightness temperature measurements. In this paper, comparison of the theoretical result calculated from a two-scale emission and scattering model with the Seasat and SSM/I model functions and aircraft measurements is presented and potential applications of microwave polarimetry are discussed. In our two-scale model, the modified reflectivities of large scale surfaces are calculated by extending the small perturbation method to the second order for small scale perturbation with anisotropic directional spectrum. It was found that the modified reflectivities derived from the second-order scattered field agree excellently well with the results obtained from a Monte Carlo simulation technique which numerically calculates the polarimetric reflectivities of one-dimensional random rough surfaces with a power-law spectrum. Without the second-order correction, the modified reflectivities of the rough surfaces are significantly over-estimated and sign errors are observed in the third and fourth Stokes parameters for thermal emissions. The surface spectrum parameters and two-scale cutoff are selected so that the calculated scattering and emission signatures agree with the reported model functions for Seasat and SSM/I. Subsequently, the polarimetric signatures of sea surfaces are illustrated to indicate the possibility of reducing the number of azimuthal looks required for spaceborne sensors in the remote sensing of ocean wind by using the polarimetric information. Furthermore, it is found that contrary to the dependence of backscattering coefficients on incidence angles, polarimetric brightness temperatures display a stronger wind direction dependence in the near nadir-looking direction than away from nadir. Finally, we discuss the potential of a single-look polarimetric radiometer for wind retrieval.

Yueh, S. H.↗
NASA NTRS1992

Comparison of the observed and calculated clear sky greenhouse effect - Implications for climate studies

The clear sky greenhouse effect is defined in terms of the outgoing longwave clear sky flux at the top of the atmosphere. Recently, interest in the magnitude of the clear sky greenhouse effect has increased due to the archiving of the clear sky flux quantity through the Earth Radiation Budget Experiment (ERBE). The present study investigates to what degree of accuracy this flux can be analyzed by using independent atmospheric and surface data in conjunction with a detailed longwave radiation model. The conclusion from this comparison is that for most regions over oceans the analyzed fluxes agree to within the accuracy of the ERBE-retrieved fluxes (+/- 5 W/sq m). However, in regions where deep convective activity occurs, the ERBE fluxes are significantly higher (10-15 W/sq m) than the calculated fluxes. This bias can arise from either cloud contamination problems or variability in water vapor amount. It is argued that the use of analyzed fluxes may provide a more consistent clear sky flux data set for general circulation modeling validation. Climate implications from the analyzed fluxes are explored. Finally, results for obtaining longwave surface fluxes over the oceans are presented.

Kiehl, J. T.↗
NASA NTRS1996

The Nature of Global Large-scale Sea Level Variability in Relation to Atmospheric Forcing: A Modeling Study

The relation between large-scale sea level variability and ocean circulation is studied using a numerical model. A global primitive equaiton model of the ocean is forced by daily winds and climatological heat fluxes corresponding to the period from January 1992 to February 1996. The physical nature of the temporal variability from periods of days to a year, are examined based on spectral analyses of model results and comparisons with satellite altimetry and tide gauge measurements.

ocean circulation circulation models ocean models ↗
NASA NTRS1993

The earth's radiation budget and its relation to atmospheric hydrology. III - Comparison of observations over the oceans with a GCM

The ability of the Colorado State University general circulation model (GCM) to simulate interactions between the hydrological cycle and the radiative processes on earth was examined by comparing various sensitivity relationships established by the model with those observed on earth, and the observed and calculated seasonal cycles of the greenhouse effect and cloud radiative forcing. Results showed that, although the GCM model used was able to simulate well some aspects of the observed sensitivities, there were many serious quantitative differences, including problems in the simulation of the column vapor in the tropics and an excessively strong clear-sky greenhouse effect in the mid-latitudes. These differences led to an underestimation by the model of the sensitivity of the clear-sky greenhouse to changes in sea surface temperature.

Stephens, Graeme L.↗
NASA NTRS1994

Temperature-moisture biases in ECMWF analyses based on clear sky longwave simulations constrained by SSMI and MSU measurements and comparisons to ERBE estimates

Clear sky longwave radiation fluxes for the summer of 1988 and winter of 1989 have been simulated with a radiative transfer model that includes detailed treatment of atmospheric gas absorption. The input data to the model are humidity and temperature profiles from European Center for Medium-Range Weather Forecasts (ECMWF) analyses and surface temperature measurements from the International Satellite Cloud Climatology Project. To reduce the inherent uncertainties in humidity profiles, the ECMWF moisture fields have been adjusted, based on the climatological relationship between the moisture profile and total precipitable water (PW), so that the ECMWF PW is equal to that derived from Special Sensor Microwave Imager (SSMI) data. SSMI and ECMWF PW patterns are generally similar, but significant differences in magnitude are found in the low latitudes, particularly over the subtropical oceans in the Southern Hemisphere. At the same time the tropospheric brightness temperature (T(sub B)) at 53.74 GHz has been computed using a microwave transfer model with ECMWF temperature and humidity fields as inputs. The comparison of computed T(sub B) with the Microwave Sounding Unit (MSU) channel 2 T(sub B) revealed differences larger than 1.5 K over most of the oceans, suggesting that the ECMWF model atmosphere has a cold bias. Maximum biases greater than 3 K are found in low latitudes. After removing the bias from the ECMWF temperature field, a new set of temperature profiles yielding brightness temperature equal to MSU channel 2 T(sub B) has been obtained. Simulation results clearly demonstrate that inclusion of satellite estimates of PW and T(sub B) enhances the accuracy of the clear sky LW flux simulation, substantially reducing the differences from Earth Radiation Budget Experiment (ERBE) estimates. Global averages indicate that the model-derived top-of-atmosphere (TOA) clear sky fluxes from the adjusted moisture and temperature profiles are in very good agreement with satellite-measured ERBE values, differing by only 1.4 W/sq m. By contrast, the calculation using only ECMWF analyses without any adjustments gives a significant underestimation of TOA clear sky fluxes up to 10.7 W/sq m compared to the satellite measurements, suggesting that the constraint method can significantly improve the accuracy in the radiation budget simulation. The results also indicate that the ECMWF-based simulation errors are, at least for the two seasons studied, mainly attributed to the error in the temperature field rather than to errors in the moisture field.

Sohn, Byung-Ju↗
NASA NTRSDate not supplied

Early Comparison of Ocean Analyses Using JEDI/SOCA 3 DVar vs. GEOS-S2S- 3 ODAS LETKF

At NASA’s Global Modeling and Assimilation Office (GMAO), data assimilation (DA) for the next-generation Goddard Earth Observing System Subseasonal-to-Seasonal (GEOS-S 2 S) coupled-model forecast system will transition to the Joint Effort for DA Integration (JEDI) system, which includes the marine DA component SOCA (Sea ice, Ocean, and Coupled Assimilation). It is envisioned that incorporating SOCA into GMAO Earth System Modeling will allow a more systematic approach to assimilating new data types (e.g., SWOT KaRIn), increasing resolution (e.g., 1/12°), and facilitating weakly (and eventually strongly) coupled DA (air-sea-ice, etc.). To prepare for this transition, testing is underway to compare SOCA results at high ocean resolution against the current GEOS-S 2 S Version 3 Ocean DA System (ODAS) 1 using similar initial conditions and observations.

David Russell↗
NASA NTRS1986

Toward the application of the Kalman filter to regional open ocean modeling

A partial differential equation model is defined for ocean meteorological prediction and synoptic analysis. The Kalman filter used for data assimilation is described and applied to the one-dimensional linear barotropic quasi-geostrophic model with periodic and open boundary conditions. The model accounts for eddy scale dynamics in the ocean. The assumptions made in the forecast model are discussed, along with comparisons of the error variances expected with the filter and from an objective analysis method. The effectiveness of the Kalman filter is demonstrated and subsequent efforts to extend the filter to two dimensions are indicated.

Miller, R. N.↗
NASA NTRS2017

Use of A-Train Aerosol Observations to Constrain Direct Aerosol Radiative Effects (DARE) Comparisons with Aerocom Models and Uncertainty Assessments

We describe a technique for combining multiple A-Train aerosol data sets, namely MODIS spectral AOD (aerosol optical depth), OMI AAOD (absorption aerosol optical depth) and CALIOP aerosol backscatter retrievals (hereafter referred to as MOC retrievals) to estimate full spectral sets of aerosol radiative properties, and ultimately to calculate the 3-D distribution of direct aerosol radiative effects (DARE). We present MOC results using almost two years of data collected in 2007 and 2008, and show comparisons of the aerosol radiative property estimates to collocated AERONET retrievals. Use of the MODIS Collection 6 AOD data derived with the dark target and deep blue algorithms has extended the coverage of the MOC retrievals towards higher latitudes. The MOC aerosol retrievals agree better with AERONET in terms of the single scattering albedo (ssa) at 441 nm than ssa calculated from OMI and MODIS data alone, indicating that CALIOP aerosol backscatter data contains information on aerosol absorption. We compare the spatio-temporal distribution of the MOC retrievals and MOC-based calculations of seasonal clear-sky DARE to values derived from four models that participated in the Phase II AeroCom model intercomparison initiative. Overall, the MOC-based calculations of clear-sky DARE at TOA over land are smaller (less negative) than previous model or observational estimates due to the inclusion of more absorbing aerosol retrievals over brighter surfaces, not previously available for observationally-based estimates of DARE. MOC-based DARE estimates at the surface over land and total (land and ocean) DARE estimates at TOA are in between previous model and observational results. Comparisons of seasonal aerosol property to AeroCom Phase II results show generally good agreement best agreement with forcing results at TOA is found with GMI-MerraV3. We discuss sampling issues that affect the comparisons and the major challenges in extending our clear-sky DARE results to all-sky conditions. We present estimates of clear-sky and all-sky DARE and show uncertainties that stem from the assumptions in the spatial extrapolation and accuracy of aerosol and cloud properties, in the diurnal evolution of these properties, and in the radiative transfer calculations.

aerosol↗
NASA NTRS1992

Precipitable water and surface humidity over global oceans from special sensor microwave imager and European Center for Medium Range Weather Forecasts

Global fields of precipitable water W from the special sensor microwave imager were compared with those from the European Center for Medium Range Weather Forecasts (ECMWF) model. They agree over most ocean areas; both data sets capture the two annual cycles examined and the interannual anomalies during an ENSO episode. They show significant differences in the dry air masses over the eastern tropical-subtropical oceans, particularly in the Southern Hemisphere. In these regions, comparisons with radiosonde data indicate that overestimation by the ECMWF model accounts for a large part of the differences. As a check on the W differences, surface-level specific humidity Q derived from W, using a statistical relation, was compared with Q from the ECMWF model. The differences in Q were found to be consistent with the differences in W, indirectly validating the Q-W relation. In both W and Q, SSMI was able to discern clearly the equatorial extension of the tongues of dry air in the eastern tropical ocean, while both ECMWF and climatological fields have reduced spatial gradients and weaker intensity.

Liu, W. T.↗
NASA NTRS2020

Comparison of Surface Fluxes Derived from CYGNSS and Simulated by WRF Model: An MJO Case Study

This study focuses on ocean surface fluxes, mainly the latent heat flux, and their impact on MJO propagation and associated precipitation structures over the Indian Ocean and Maritime Continent. The Coupled-Ocean-Atmosphere-Wave-Sediment Transport (COAWST) model is used to simulate two MJO events during the 2017-2018 season: the December 10 - January 20, 2017 case, which maintained its strong precipitation signal over the Maritime Continent, and the March 1 - 20, 2018 case, which was weaker and did not propagate through the Maritime Continent. Both simulated MJO events show positive biases in surface rainfall compared with GPM IMERG data. During the MJO suppressed phase, the simulations rain more often than the observations. During the active phase, the westward propagating precipitation structures are more organized and much stronger compared with the observations, sometimes forming westward propagating cyclones that weakened the eastward precipitation signals. Two aspects of the surface flux interactions are investigated: the impact of the domain mean surface fluxes, and the impact of storm scale circulations and their interactions with local surface fluxes. Both aspects affect water vapor budget, atmosphere instability and mean flow, through which convection initiation, organization, and propagation are influenced. Model sensitivity tests with different radiation, microphysics, PBL schemes and nudging schemes indicate that in the control simulations, higher SST and surface fluxes, especially during the suppressed period, are the main reason of rainfall overestimation compared with IMERG data. The strong westward propagating signals are caused by both increased atmosphere instability and reduced mean wind shear. Unfortunately, the small differences in mean SST and surface fluxes between different model sensitivity tests are all within the satellite observation error margin, and cannot be directly corroborated by observations. One of the advantages of CYGNSS satellites is that they observe ocean surface wind and heat fluxes underneath strong rainfall events such as the convective systems associated with MJO active phases. Currently we are comparing CYGNSS level 2 surface fluxes retrievals and the model simulations in order to better understand the second aspect of the MJO and surface fluxes interactions, and how this affects MJO strengths and propagations. The interactive atmosphere-ocean-wave model also provides cases that directly comparing satellite observables (the bistatic radar cross section) and the model simulations (through CYGNSS satellite simulator). These discrepancies are more prominent in coupled ocean simulations, mainly due to higher SST and enhanced surface fluxes.

Li, Xiaowen↗
NASA NTRS1978

Analysis and geological interpretation of gravity data from GEOS-3 altimeter

A number of detailed gravimetric geoids of portions of the world's oceans from marine gravity measurements were constructed. The geoids were constructed by computing 1 x 1 deg or 10 x 10 deg averages of free-air anomaly data and subtracting these values from currently used satellite derived Earth models. The resulting difference gravity anomalies are then integrated over a sphere using a simplified form of Stoke's equation to obtain a difference geoid. This difference geoid is added to the satellite derived model to obtain a 1 x 1 deg or 10 x 10 deg total gravimetric geoid. The geoid undulations are studied by comparison of the altimeter measurements with the morphology of the ocean floor. Utilizing a combination of altimetry data, gravity and seismic reflection data, geophysical models of the earth can be constructed.

Talwani, M.↗
NASA NTRS2024

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↗