Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Ensemble models”

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 343 records · Page 19

Impacts of 1.5 °C and 2.0 °C global warming above pre-industrial on potential winter wheat production of China

Keeping global temperatures below 2.0 °C above pre-industrial condition and pursuing efforts toward the more ambitious 1.5 °C goal in the late 21st century was the main target from the Paris Agreement in 2015. Here we assessed the likely challenges for the China’s winter wheat production under 1.5 °C and 2.0 °C increase of global temperature, with four wheat crop models (CERES-Wheat, Nwheat, WheatGrow, and APSIM-Wheat) and the latest climate projections from the Half a degree Additional warming, Prognosis and Projected Impacts project (HAPPI). Instead of using average “winter type” wheat cultivar, and same management and soil inputs for whole region, location-specific winter wheat cultivars with local agronomic information were calibrated for each of the representative wheat growing area of China, allowing a better spatial agronomic representation of the whole wheat planting area. The mean growing season temperature (GST) during the winter wheat vegetative stage was projected to increase by 0.6 to 1.4 °C for the 1.5 °C scenario, and 0.9 to 1.8 °C for the 2.0 °C scenario, while during the reproductive stage was decreased between 0 and 0.9 °C for the 1.5 °C scenario and -0.3 and 1.1 °C for the 2.0 °C scenario. Growing season duration (GSD) for the whole period was shortened by 6 to 15 days for the 1.5 °C scenario and 8 to 18 days for the 2.0 °C scenario, as a result of higher GST under global warming. Increase in GST and decrease in GSD was more obvious in the Southwest Subregion (SWS) than subregions in the north. The shortening GSD for the whole wheat growth period was mostly from the shortening vegetative period, as no appreciable difference in number of days from anthesis to maturity was found for the whole regions. Although there is variability among models, the indication is that wheat yields were projected to increase in the North Subregion (NS), the Huang-Huai Subregion (HHS), and the Middle-lower Researches of Yangzi River Subregion (MYS), but to decrease in the SWS under two warming scenarios. The effects of elevated CO2 concentration were mostly beneficial and tended to offset the negative impacts of increasing temperature at both global warming scenarios, with a rate of 7-14% yield increase per 100-ppm, except for locations with GST of baseline higher than 11 °C. Aggregating to regional wheat production, the total winter wheat production of China was projected to increase by 2.8% (1.6% to 3.0%, 25th percentile to 75th percentile) and 8.3% (7.0% to 9.6%, 25th percentile to 75th percentile) under 1.5 °C and 2.0 °C scenarios, and most of increase was observed in the north subregions due to the largest wheat planting area. Our results will lay the foundation for developing adaptation strategies to future climate change to ensure China and global wheat supply and food security.

Winter wheat↗

Quantifying the role of ocean coupling in Arctic amplification and sea-ice loss over the 21st century

The enhanced warming of the Arctic, relative to other parts of the Earth, a phenomenon known as Arctic amplification, is one of the most striking features of climate change, and has important climatic impacts for the entire Northern Hemisphere. Several mechanisms are believed to be responsible for Arctic amplification; however, a quantitative understanding of their relative importance is still missing. Here, using ensembles of model integrations, we quantify the contribution of ocean coupling, both its thermodynamic and dynamic components, to Arctic amplification over the 20th and 21st centuries. We show that ocean coupling accounts for ~80% of the amplification by 2100. In particular, we show that thermodynamic coupling is responsible for future amplification and sea-ice loss as it overcomes the effect of dynamic coupling which reduces the amplification and sea-ice loss by ~35%. Our results demonstrate the utility of targeted numerical experiments to quantify the role of specific mechanisms in Arctic amplification, for better constraining climate projections.

Ocean heat transport↗

Dissemination of Global Flood Information Through Nasa Disasters Program Global Initiative for Flood Forecasting and Alerting (GIFFT)

The Global Initiative for Flood Forecasting and Alerting (GIFFT) is an initiative to disseminate flood products from different sensors to global stakeholders via the Pacific Disaster Center’s (PDC) DisasterAWARE®, the NASA Disasters Mapping Portal and potentially other mechanisms. GIFFT will also generate an integrated product(s) using ensembled hydrologic models and Earth observation derived flood outputs to stakeholders globally. This initiative combines the results of several Disasters Program funded projects that focus on flood forecasting, flood severity estimation, post-event flood mapping, flooddepth estimation and resulting flood severity using remote sensing-based flood information derived from SAR and optical imagery. PDC’s DisasterAWARE platform is a global multi-hazard alerting platform that currently incorporatesModel of Models (MoM) outputs as flood “incidents,” visually depicting potential floods in the context ofpopulation and infrastructure that may become affected. Automated procedures are being developed tocategorize MoM outputs as DisasterAWARE “hazards,” allowing for their dissemination to users alongwith other flood products that assess potential impacts. PDC’s stakeholders include UN, DoD, FEMA andtheir equivalents around the world. Additionally, the DisasterALERT mobile platform has 2 million usersworldwide and is a fully functional and sustained platform used to disseminate alerts and analyticalproducts for 18 different hazard types. GIFFT has incorporated HydroSAR products that will bedistributed via DisasterAWARE.

remote sensing↗

IAC Abstract for Orbital Debris Phase II Report

Orbital debris poses a growing risk for spacecraft operations. This risk has most often been measured by proxies, like the number of debris or expected collisions, but few have attempted quantify the consequences of orbital debris or the costs and benefits of different methods to address the problem. This has made it difficult to assess the trade-offs and relative value among proposed actions for the mitigation, tracking, and remediation of debris. In March of 2023, the NASA Office of Technology, Policy, and Strategy (OTPS) took a step toward quantifying these relative values when we released Cost and Benefit Analysis of Orbital Debris Remediation. Since that time, OTPS has continued to build on this work. We have increased the fidelity of our risk model. For example, we added an ensemble of risk models that use varying assumptions of debris populations and lethality, included debris as small as 1mm, and incorporated atmospheric drag to naturally deorbit debris. We have also expanded our cost-benefit calculations to include action that mitigate and track debris, both large and small. The analysis has several important findings regarding the conditions under which debris remediation is more cost-effective than mitigation and tracking, the value of increased shielding of spacecraft and faster timelines for post-mission disposal, and how changes in the portfolio have interdependent effects. This presentation summarizes the body of work, with a particular focus on the risks and benefits to the global spacecraft population.

orbital debris↗

Towards an improved understanding of the Antarctic coastal zone and its contribution to future global sea level

Understanding the coastal zone of the Antarctic Ice Sheet, where it interacts with the Southern Ocean and warmer air masses, is crucial for predicting Antarctica's influence on the global climate. This region has multiple tipping mechanisms that could trigger large, rapid, and potentially irreversible changes in the coming centuries. The Antarctic Ice Sheet remains the largest source of uncertainty in future sea-level projections. Insufficient knowledge of bed topography beneath the ice shelves and the coastal ice sheet is not yet well documented, but is a major source of this uncertainty. This review assesses current knowledge of the coastal zone and highlights methods to investigate it, including aerogeophysical surveys, ground- and ship-based measurements, satellite observations, and computer modeling. An ensemble analysis of published bed topography datasets identifies significant data gaps and their regional distribution, framed in the context of current ice-sheet behavior and potential instability. We propose scientific priorities and guidelines for future aerogeophysical surveys, advocating for a comprehensive, coordinated international effort to build a next-generation dataset of Antarctic bed properties. Such an initiative would significantly advance understanding of the role of coastal processes in ice-sheet dynamics, reducing uncertainties in sea-level rise projections and enhancing predictions of future ocean and climate changes.

Kenichi Matsuoka↗

Towards an Improved Understanding of the Antarctic Coastal Zone and Its Contribution to Future Global Sea Level

Understanding the coastal zone of the Antarctic Ice Sheet (AIS), where it interacts with the Southern Ocean and warmer air masses, is crucial for predicting Antarctica's influence on the global climate and sea level. This region has multiple tipping mechanisms that could trigger large, rapid, and potentially irreversible changes in the AIS, the Southern Ocean and their global connections in the coming centuries. The AIS remains the largest source of uncertainty in future sea-level projections. Bed topography beneath the ice shelves and the coastal ice sheet is not yet well documented, and is a major source of this uncertainty. This review assesses current knowledge of the coastal zone and highlights methods to investigate it, including aerogeophysical surveys, ground- and ship-based measurements, satellite observations, and computer modeling. An ensemble analysis of published bed topography data sets identifies significant data gaps and their regional distribution, framed in the context of current ice-sheet behavior and potential instability. We propose scientific priorities and guidelines for future aerogeophysical surveys, advocating for a comprehensive, coordinated international effort to build a next-generation data set of Antarctic bed properties. Such an initiative would significantly advance understanding of the role of coastal processes in ice-sheet dynamics, reducing uncertainties in sea-level rise projections and improving predictions of future ocean and climate changes.

Kenichi Matsuoka↗

Evaluation of ACCMIP Outgoing Longwave Radiation from Tropospheric Ozone Using TES Satellite Observations.

We use simultaneous observations of tropospheric ozone and outgoing longwave radiation (OLR) sensitivity to tropospheric ozone from the Tropospheric Emission Spectrometer (TES) to evaluate model tropospheric ozone and its effect on OLR simulated by a suite of chemistry-climate models that participated in the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP). The ensemble mean of ACCMIP models show a persistent but modest tropospheric ozone low bias (5-20 ppb) in the Southern Hemisphere (SH) and modest high bias (5-10 ppb) in the Northern Hemisphere (NH) relative to TES ozone for 2005-2010. These ozone biases have a significant impact on the OLR. Using TES instantaneous radiative kernels (IRK), we show that the ACCMIP ensemble mean tropospheric ozone low bias leads up to 120mW/ sq. m OLR high bias locally but zonally compensating errors reduce the global OLR high bias to 39+/- 41mW/ sq. m relative to TES data. We show that there is a correlation (Sq. R = 0.59) between the magnitude of the ACCMIP OLR bias and the deviation of the ACCMIP preindustrial to present day (1750-2010) ozone radiative forcing (RF) from the ensemble ozone RF mean. However, this correlation is driven primarily by models whose absolute OLR bias from tropospheric ozone exceeds 100mW/ sq. m. Removing these models leads to a mean ozone radiative forcing of 394+/- 42mW/ sq. m. The mean is about the same and the standard deviation is about 30% lower than an ensemble ozone RF of 384 +/- 60mW/ sq. m derived from 14 of the 16 ACCMIP models reported in a companion ACCMIP study. These results point towards a profitable direction of combining satellite observations and chemistry-climate model simulations to reduce uncertainty in ozone radiative forcing.

troposphere↗

Meridional Propagation of the MJO/ISO and Asian Monsoon Variability

In this study we examine the links between tropical heating, the Madden Julian Oscillation (MJO)/Intraseasonal Oscillation (ISO), and the Asian monsoon. We are particularly interested in isolating the nature of the poleward propagation of the ISO/MJO in the monsoon region. We examine both observations and idealized "MJO heating" experiments employing the NASA Seasonal-Interannual Prediction Project (NSIPP) atmospheric general circulation model (AGCM). In the idealized 10-member ensemble simulations, the model is forced by climatological SST and an idealized eastward propagating heating profile that is meant to mimic the canonical heating associated with the MJO in the Indian Ocean and western Pacific. In order to understand the impact of SST on the off equatorial convection (or Rossby-wave response), a second set of 10-member ensemble simulations is carried out with the climatological SSTs shifted in time by 6-months. The observational analysis highlights the strong link between the Indian summer monsoon and the tropical ISO/MJO activity and heating. This includes the well-known meridional propagation that affects the summer monsoons of both hemispheres. The AGCM experiments with the idealized eastward propagating MJO-like heating reproduce the observed meridional propagation including the observed seasonal differences. The impact of the SSTs are to enhance the magnitude of the propagation into the summer hemispheres. The results suggest that the winter/summer differences associated with the MJO/ISO are auxiliary features that depend on the MJO's environment (basic state and boundary conditions) and are not the result of fundamental differences in the MJO itself.

Wu, Man Li C.↗

Mid-Pliocene El Niño/Southern Oscillation Suppressed By Pacific Intertropical Convergence Zone Shift

The El Niño/Southern Oscillation (ENSO), the dominant driver of year-to-year climate variability in the equatorial Pacific Ocean, impacts climate pattern across the globe. However, the response of the ENSO system to past and potential future temperature increases is not fully understood. Here we investigate ENSO variability in the warmer climate of the mid-Pliocene (~3.0–3.3 Ma), when surface temperatures were ~2–3 °C above modern values, in a large ensemble of climate models—the Pliocene Model Intercomparison Project. We show that the ensemble consistently suggests a weakening of ENSO variability, with a mean reduction of 25% (±16%). We further show that shifts in the equatorial Pacific mean state cannot fully explain these changes. Instead, ENSO was suppressed by a series of off-equatorial processes triggered by a northward displacement of the Pacific intertropical convergence zone: weakened convective feedback and intensified Southern Hemisphere circulation, which inhibit various processes that initiate ENSO. The connection between the climatological intertropical convergence zone position and ENSO we find in the past is expected to operate in our warming world with important ramifications for ENSO variability.

Atmospheric dynamics↗

Evaluation of Clouds, Radiation, and Precipitation in CMIP6 Models Using Global Weather States Derived from ISCCP-H Cloud Property Data

A clustering methodology is applied to cloud optical depth cloud top pressure (TAU-PC) histograms from the new, 1-degree resolution, ISCCP-H dataset, to derive an updated global Weather State (WS) dataset. Then, PC-TAU histograms from current-climate CMIP6 model simulations are assigned to the ISCCP-H WSs along with their concurrent radiation and precipitation properties, to evaluate model cloud, radiation, and precipitation properties in the context of the Weather States. The new ISCCP-H analysis produces WSs that are very similar to those previously found in the lower resolution ISCCP-D dataset. The main difference lies in the splitting of the ISCCP-D thin stratocumulus WS between the ISCCP-H shallow cumulus and stratocumulus WSs, which results in the reduction by one of the total WS number. The evaluation of the CMIP6 models against the ISCCP-H Weather States, shows that, in the ensemble mean, the models are producing an adequate representation of the frequency and geographical distribution of the WSs, with measurable improvements compared to the WSs derived for the CMIP5 ensemble. However, the frequency of shallow cumulus clouds continues to be underestimated, and, in some WSs the good agreement of the ensemble mean with observations comes from averaging models that significantly overpredict and underpredict the ISCCP-H WS frequency. In addition, significant biases exist in the internal cloud properties of the model WSs, such as the model underestimation of cloud fraction in middle-top clouds and secondarily in midlatitude storm and stratocumulus clouds, that result in an underestimation of cloud SW cooling in those regimes.

Clouds↗

On the Role of SST Forcing in the 2011 and 2012 Extreme U.S. Heat and Drought: A Study in Contrasts

This study compares the extreme heat and drought that developed over the United States in 2011 and 2012 with a focus on the role of SST forcing. Experiments with the NASA GEOS-5 atmospheric general circulation model show that the winter/spring response over the U.S. to the Pacific SST is remarkably similar for the two years despite substantial differences in the tropical Pacific SST. As such, the pronounced winter and early spring temperature differences between the two years (warmth confined to the south in 2011 and covering much of the continent in 2012) primarily reflect differences in the contributions from the Atlantic and Indian Oceans, with both acting to cool the east and upper mid-west during 2011, while during 2012 the Indian Ocean reinforced the Pacific-driven continental-wide warming and the Atlantic played a less important role. During late spring and summer of 2011 the tropical Pacific SST force a continued warming and drying over the southern U.S., though considerably weaker than observed. Nevertheless, the observed anomalies fall within the models intra-ensemble spread. In contrast, the rapid development of intense heat and drying over the central U.S. during June and July of 2012 falls outside the models intra-ensemble spread. The response to the SST (a northward expansion of a modest summer warming linked to the Atlantic) gives little indication that 2012 would produce record-breaking precipitation deficits and heat in the central Great Plains. A diagnosis of the 2012 observed circulation anomalies shows that the most extreme heat and drought was tied to the development of a stationary Rossby wave and an associated anomalous upper tropospheric high maintained by weather transients.

SST Forcing↗

Validation of Model Forecasts of the Ambient Solar Wind

Independent and automated validation is a vital step in the progression of models from the research community into operational forecasting use. In this paper we describe a program in development at the CCMC to provide just such a comprehensive validation for models of the ambient solar wind in the inner heliosphere. We have built upon previous efforts published in the community, sharpened their definitions, and completed a baseline study. We also provide first results from this program of the comparative performance of the MHD models available at the CCMC against that of the Wang-Sheeley-Arge (WSA) model. An important goal of this effort is to provide a consistent validation to all available models. Clearly exposing the relative strengths and weaknesses of the different models will enable forecasters to craft more reliable ensemble forecasting strategies. Models of the ambient solar wind are developing rapidly as a result of improvements in data supply, numerical techniques, and computing resources. It is anticipated that in the next five to ten years, the MHD based models will supplant semi-empirical potential based models such as the WSA model, as the best available forecast models. We anticipate that this validation effort will track this evolution and so assist policy makers in gauging the value of past and future investment in modeling support.

Macneice, P. J.↗

Investigating Sources of Variability and Error in Simulations of Carbon Dioxide in an Urban Region

Greenhouse gas (GHG) emissions estimation methods that use atmospheric trace gas observations, including inverse modeling techniques, perform better when carbon dioxide (CO2) fluxes are more accurately transported and dispersed in the atmosphere by a numerical model. In urban areas, transport and dispersion is particularly difficult to simulate using current mesoscale meteorological models due, in part, to added complexity from surface heterogeneity and fine spatial/temporal scales. It is generally assumed that the errors in GHG estimation methods in urban areas are dominated by errors in transport and dispersion. Other significant errors include, but are not limited to, those from assumed emissions magnitude and spatial distribution. To assess the predictability of simulated trace gas mole fractions in urban observing systems using a numerical weather prediction model, we employ an Eulerian model that combines traditional meteorological variables with multiple passive tracers of atmospheric CO2 from anthropogenic inventories and a biospheric model. The predictability of the Eulerian model is assessed by comparing simulated atmospheric CO2 mole fractions to observations from four in situ tower sites (three urban and one rural) in the Washington DC/Baltimore, MD area for February 2016. Four different gridded fossil fuel emissions inventories along with a biospheric flux model are used to create an ensemble of simulated atmospheric CO2 observations within the model. These ensembles help to evaluate whether the modeled observations are impacted more by the underlying emissions or transport. The spread of modeled observations using the four emission fields indicates the model's ability to distinguish between the different inventories under various meteorological conditions. Overall, the Eulerian model performs well; simulated and observed average CO2 mole fractions agree within 1% when averaged at the three urban sites across the month. However, there can be differences greater than 10% at any given hour, which are attributed to complex meteorological conditions rather than differences in the inventories themselves. On average, the mean absolute error of the simulated compared to actual observations is generally twice as large as the standard deviation of the modeled mole fractions across the four emission inventories. This result supports the assumption, in urban domains, that the predicted mole fraction error relative to observations is dominated by errors in model meteorology rather than errors in the underlying fluxes in winter months. As such, minimizing errors associated with atmospheric transport and dispersion may help improve the performance of GHG estimation models more so than improving flux priors in the winter months. We also find that the errors associated with atmospheric transport in urban domains are not restricted to certain times of day. This suggests that atmospheric inversions should use CO2 observations that have been filtered using meteorological observations rather than assuming that meteorological modeling is most accurate at certain times of day (such as using only mid-afternoon observations).

CO2↗

Reduced El NiñO Variability in the Mid-Pliocene According to the Pliomip2 Ensemble

The mid-Pliocene warm period (3.264–3.025 Ma) is the most recent geological period during which atmospheric CO 2 levels were similar to recent historical values (∼400 ppm). Several proxy reconstructions for the mid-Pliocene show highly reduced zonal sea surface temperature (SST) gradients in the tropical Pacific Ocean, indicating an El Niño-like mean state. However, past modelling studies do not show these highly reduced gradients. Efforts to understand mid-Pliocene climate dynamics have led to the Pliocene Model Intercomparison Project (PlioMIP). Results from the first phase (PlioMIP1) showed clear El Niño variability (albeit significantly reduced) and did not show the greatly reduced time-mean zonal SST gradient suggested by some of the proxies. In this work, we study El Niño–Southern Oscillation (ENSO) variability in the PlioMIP2 ensemble, which consists of additional global coupled climate models and updated boundary conditions compared to PlioMIP1. We quantify ENSO amplitude, period, spatial structure and “flavour”, as well as the tropical Pacific annual mean state in mid-Pliocene and pre-industrial simulations. Results show a reduced ENSO amplitude in the model-ensemble mean (−24 %) with respect to the pre-industrial, with 15 out of 17 individual models showing such a reduction. Furthermore, the spectral power of this variability considerably decreases in the 3–4-year band. The spatial structure of the dominant empirical orthogonal function shows no particular change in the patterns of tropical Pacific variability in the model-ensemble mean, compared to the pre-industrial. Although the time-mean zonal SST gradient in the equatorial Pacific decreases for 14 out of 17 models (0.2 ∘C reduction in the ensemble mean), there does not seem to be a correlation with the decrease in ENSO amplitude. The models showing the most “El Niño-like” mean state changes show a similar ENSO amplitude to that in the pre-industrial reference, while models showing more “La Niña-like” mean state changes generally show a large reduction in ENSO variability. The PlioMIP2 results show a reasonable agreement with both time-mean proxies indicating a reduced zonal SST gradient and reconstructions indicating a reduced, or similar, ENSO variability.

mid-Pliocene warm period↗

Inconsistent Strategies to Spin up Models in CMIP5: Implications for Ocean Biogeochemical Model Performance Assessment

During the fifth phase of the Coupled Model Intercomparison Project (CMIP5) substantial efforts were made to systematically assess the skill of Earth system models. One goal was to check how realistically representative marine biogeochemical tracer distributions could be reproduced by models. In routine assessments model historical hindcasts were compared with available modern biogeochemical observations. However, these assessments considered neither how close modeled biogeochemical reservoirs were to equilibrium nor the sensitivity of model performance to initial conditions or to the spin-up protocols. Here, we explore how the large diversity in spin-up protocols used for marine biogeochemistry in CMIP5 Earth system models (ESMs) contributes to model-to-model differences in the simulated fields. We take advantage of a 500-year spin-up simulation of IPSL-CM5A-LR to quantify the influence of the spin-up protocol on model ability to reproduce relevant data fields. Amplification of biases in selected biogeochemical fields (O2, NO3, Alk-DIC) is assessed as a function of spin-up duration. We demonstrate that a relationship between spin-up duration and assessment metrics emerges from our model results and holds when confronted with a larger ensemble of CMIP5 models. This shows that drift has implications for performance assessment in addition to possibly aliasing estimates of climate change impact. Our study suggests that differences in spin-up protocols could explain a substantial part of model disparities, constituting a source of model-to- model uncertainty. This requires more attention in future model intercomparison exercises in order to provide quantitatively more correct ESM results on marine biogeochemistry and carbon cycle feedbacks.

histories↗

An Assessment of Multimodel Simulations for the Variability of Western North Pacific Tropical Cyclones and Its Association with ENSO

An assessment of simulations of the interannual variability of tropical cyclones (TCs) over the western North Pacific (WNP) and its association with El Niño-Southern Oscillation (ENSO), as well as a subsequent diagnosis for possible causes of model biases generated from simulated large-scale climate conditions, are documented in the paper. The model experiments are carried out by the Hurricane Work Group under the U.S. Climate Variability and Predictability Research Program (CLIVAR) using five global climate models (GCMs) with a total of 16 ensemble members forced by the observed sea surface temperature and spanning the 28-yr period from 1982 to 2009. The results show GISS and GFDL model ensemble means best simulate the interannual variability of TCs, and the multimodel ensemble mean (MME) follows. Also, the MME has the closest climate mean annual number of WNP TCs and the smallest root-mean-square error to the observation. Most GCMs can simulate the interannual variability of WNP TCs well, with stronger TC activities during two types of El Niño-namely, eastern Pacific (EP) and central Pacific (CP) El Niño-and weaker activity during La Niña. However, none of the models capture the differences in TC activity between EP and CP El Niño as are shown in observations. The inability of models to distinguish the differences in TC activities between the two types of El Niño events may be due to the bias of the models in response to the shift of tropical heating associated with CP El Niño.

climate models↗

Hybrid Data Assimilation without Ensemble Filtering

The Global Modeling and Assimilation Office is preparing to upgrade its three-dimensional variational system to a hybrid approach in which the ensemble is generated using a square-root ensemble Kalman filter (EnKF) and the variational problem is solved using the Grid-point Statistical Interpolation system. As in most EnKF applications, we found it necessary to employ a combination of multiplicative and additive inflations, to compensate for sampling and modeling errors, respectively and, to maintain the small-member ensemble solution close to the variational solution; we also found it necessary to re-center the members of the ensemble about the variational analysis. During tuning of the filter we have found re-centering and additive inflation to play a considerably larger role than expected, particularly in a dual-resolution context when the variational analysis is ran at larger resolution than the ensemble. This led us to consider a hybrid strategy in which the members of the ensemble are generated by simply converting the variational analysis to the resolution of the ensemble and applying additive inflation, thus bypassing the EnKF. Comparisons of this, so-called, filter-free hybrid procedure with an EnKF-based hybrid procedure and a control non-hybrid, traditional, scheme show both hybrid strategies to provide equally significant improvement over the control; more interestingly, the filter-free procedure was found to give qualitatively similar results to the EnKF-based procedure.

Kalman Filter↗

Multi-source Estimates of Land / Ocean Moisture Transport Variability over the Satellite Era

It is widely appreciated that atmospheric transport of water from the world’s oceans is a process key to planetary energy balance as well as Earth's habitability. What is not yet clear is the extent of variability in moisture transports, the relative importance of interdecadal variability versus climate change signals, and importantly, our ability to quantify these changes. This work assesses variations in moisture transport variability during the satellite era (~1980 to present) by comparing several different estimates. (i) The most direct estimate is the vertically integrated flux convergence of moisture from reanalyses which use observed wind and moisture information. (ii) One alternative estimate comes from P-ET over land taken from global hydrologic models constrained with precipitation and near-surface meteorology. Here we use an ensemble of six models. An adjunct to this method is to employ satellite derive ET (e.g., GLEAM or DOLCE). (iii) Complementary to this is E-P over the global oceans derived from satellite estimates of P such as TRMM, GPM and GPCP and SeaFlux V3 or J-OFURO3 estimates of E, all relying heavily upon microwave measurements. Transport between land and oceans must essentially balance at monthly scales, i.e., vanish globally. (iv) a fourth perspective comes from estimate of terrestrial RO + storage rate, delta S. G-RUN Ensemble which uses observed streamflow and P measurements to calibrate a statistical model provides the former while GRACE, GRACE-FO provide total water storage anomalies used to calculate storage rate changes. GRACE REC uses GRACE data to train a precipitation-driven statistical model to extend storage estimates before the GRACE era. (All of these alternatives to reanalysis estimates also consider the small atmospheric column water vapor contribution.) We examine the transport changes from these three different methodologies, their relative accuracies and discuss the origin of their differences. Regional trends in moisture flux divergence and their role in multi-decadal trends are considered. Interannual variability arising in connection with ENSO variability is a dominant signal, driven largely by P changes. Trends since 1980 include reductions in moisture delivery to the western U.S., eastern Brazil, and central Africa with recovery of moisture convergence to the Sahel and parts of eastern North America.

Franklin Robertson↗