Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reanalysis data”

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 811 records · Page 45

Estimating UV-B, UV-Erithemic, and UV-A Irradiances From Global Horizontal Irradiance and MERRA-2 Ozone Column Information

The ground ultraviolet (UV) solar radiation is relevant due to its impacts on plastics degradation (mainly UVA) and on human health (UVB and erithemic UV (UVE)). UV ground measurements are not as ubiquitous as the relatively common global horizontal irradiance (GHI) measurements. Three simple models that estimate the UVA, UVB, and UVE components of solar irradiance from GHI and ozone column information are locally adjusted and validated. Five one-minute datasets from three sites in southeastern South America and two in the United States are used for simultaneous solar irradiance and UV data. All sites correspond to temperate mid-latitude regions. Simultaneous atmospheric total ozone column information is obtained from the reanalysis modern-era retrospective analysis for research and applications (MERRA-2) database for each site. Aside from locally adjusted models, average models with a single set of coefficients are also evaluated. For instance, the best average model is able to estimate UVE with a typical uncertainty below 12% and mean biases between +-3%, relative to the average of the measurements. Similar results are reported for the UVB and UVA components. These results, which can be useful in regions with similar climate and geography, provide a simple way to estimate UV irradiance under all-sky conditions with known uncertainty. This is an alternative to global satellite-based UV estimates, which can have high uncertainties at specific locations. Because MERRA-2 information has a global coverage, when coupled with good satellite-based estimates for GHI, UV irradiances can be estimated by this method over a large territory.

environmental UV radiation↗

Characterizing the Influence of the General Circulation on Marine Boundary Layer Clouds

The seasonal and intraseasonal variability of boundary layer cloud in the subtropical eastern oceans are studied using combined data from the International Satellite Cloud Climatology Project (ISCCP) and the European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis. Spectral analysis reveals that most of the time variability of cloud properties occurs on seasonal to annual time scales. The variance decreases one to two orders of magnitude for each decade of time scale decrease, indicating that daily to monthly time scales have smaller, but non-negligible variability. The length of these dominant time scales suggests that the majority of the variability is influenced by the general circulation and its interaction with boundary layer turbulence, rather than a product of boundary layer turbulence alone. Previous datasets have lacked the necessary resolution in either time or in space to properly characterize variability on synoptic scales; this is remedied by using global satellite-retrieved cloud properties. We characterize the intraseasonal subtropical cloud variability in both hemispheres and in different seasons. In addition to cloud fraction, we examine variability of cloud optical thickness - cloud top pressure frequency distributions. Despite the large concentration of research on the variability of Northern Hemisphere (NH) regions during summer, it is noted that the largest amplitude intraseasonal variability in the NH regions occurs during local winter. The effect of intraseasonal variability on the calculation and interpretation of seasonal results is investigated. Decreases in seasonally averaged cloud cover, optical thickness and cloud top pressure from the May-through-September season to the November-through-March season are most apparent in the NH regions. Further analysis indicates that these changes are due to an increase in frequency, but a decrease in the persistence of synoptic events. In addition, changes in cloud top pressure and optical thickness characteristics from the summer to winter seasons indicate that the NH subtropics undergo a change in dynamic regime with season. This change appears in the cloud fields as a shift from the more commonly seen lower-altitude, thicker optical thickness clouds to higher-altitude, thinner clouds. The latter cloud-type is associated with the lower sea level pressure, upward vertical velocity phase of the synoptic wave. Intraseasonal changes in cloud properties in the Southern Hemisphere and NH summer are much smaller in amplitude. Although they also appear to be linked to changes in the large-scale dynamics, similarly to NH winter variations, the relationships are more ambiguous due to the small amplitudes and longer time scales. We attempt to interpret some of these relationships using the results of the Betts and Ridgway (1989) box model. However, these results cannot consistently explain the patterns when results from all regions are considered, implying that this model may not adequately explain all the processes involved in the variability.

Rozendaal, Margaret A.↗

Regional Model Nesting Within GFS Daily Forecasts Over West Africa

The study uses the RM3, the regional climate model at the Center for Climate Systems Research of Columbia University and the NASA/Goddard Institute for Space Studies (CCSR/GISS). The paper evaluates 30 48-hour RM3 weather forecasts over West Africa during September 2006 made on a 0.5 grid nested within 1 Global Forecast System (GFS) global forecasts. September 2006 was the Special Observing Period #3 of the African Monsoon Multidisciplinary Analysis (AMMA). Archived GFS initial conditions and lateral boundary conditions for the simulations from the US National Weather Service, National Oceanographic and Atmospheric Administration were interpolated four times daily. Results for precipitation forecasts are validated against Tropical Rainfall Measurement Mission (TRMM) satellite estimates and data from the Famine Early Warning System (FEWS), which includes rain gauge measurements, and forecasts of circulation are compared to reanalysis 2. Performance statistics for the precipitation forecasts include bias, root-mean-square errors and spatial correlation coefficients. The nested regional model forecasts are compared to GFS forecasts to gauge whether nesting provides additional realistic information. They are also compared to RM3 simulations driven by reanalysis 2, representing high potential skill forecasts, to gauge the sensitivity of results to lateral boundary conditions. Nested RM3/GFS forecasts generate excessive moisture advection toward West Africa, which in turn causes prodigious amounts of model precipitation. This problem is corrected by empirical adjustments in the preparation of lateral boundary conditions and initial conditions. The resulting modified simulations improve on the GFS precipitation forecasts, achieving time-space correlations with TRMM of 0.77 on the first day and 0.63 on the second day. One realtime RM3/GFS precipitation forecast made at and posted by the African Centre of Meteorological Application for Development (ACMAD) in Niamey, Niger is shown.

Druyan, Leonard M.↗

MERRA-2: File Specification

The second Modern-Era Retrospective analysis for Research and Applications (MERRA-2) is a NASA atmospheric reanalysis that begins in 1980. It replaces the original MERRA reanalysis (Rienecker et al., 2011) using an upgraded version of the Goddard Earth Observing System Model, Version 5 (GEOS-5) data assimilation system. The file collections for MERRA-2 are described in detail in this document, including some important changes from those of the MERRA dataset (Lucchesi, 2012).

MERRA-2↗

Exploring Climatology and Long-Term Variations of Aerosols from NASA Reanalysis MERRA-2 with Giovanni

Dust plays important roles in energy cycle and climate variations. The dust deposition is the major source of iron in the open ocean, which is an essential micronutrient for phytoplankton growth and therefore may influence the ocean uptake of atmospheric CO2. Mineral dust can also act as fertilizer for forests over long time periods. Over 35 years of simulated global aerosol products from NASA atmospheric reanalysis, second Modern-Era Retrospective analysis for Research and Applications (MERRA-2) are available from NASA Goddard Earth Science Data and Information Services Center (GES DISC). The MERRA-2 covers the period 1980-present, continuing as an ongoing climate analysis. Aerosol assimilation is included throughout the period, using MODIS, MISR, AERONET, and AVHRR (in the pre-EOS period). The aerosols are assimilated by using MERRA-2 aerosol model, which interact directly with the radiation parameterization, and radiatively coupled with atmospheric model dynamics in the Goddard Earth Observing System Model, Version 5 (GEOS-5). Dust deposition data along with other major aerosol compositions (e.g. black carbon, sea salt, and sulfate, etc.) are simulated as dry and wet deposition, respectively. The hourly and monthly data are available at spatial resolution of 0.5ox0.625o (latitude x longitude). Quick data exploration of climatology and interannual variations of MERRA-2 aerosol can be done through the online visualization and analysis tool, Giovanni. This presentation, using dust deposition as an example, demonstrates a number of MERRA-2 data services at GES DISC. Global distributions of dust depositions, and their seasonal and inter-annual variations are investigated from MERRA-2 monthly aerosol products.

data services↗

File Specification for MERRA-2 Stratospheric Composition Reanalysis of Aura MLS (M2-SCREAM)

This document provides a brief description of the output collections from the MERRA-2 Stratospheric Composition Reanalysis of the Aura Microwave Limb Sounder (M2-SCREAM) produced at NASA’s Global Modeling and Assimilation Office. These data are generated by assimilating retrievals from the Microwave Limb Sounder (MLS) and the Ozone Monitoring Instrument (OMI) into the Global Earth Observing System (GEOS) Constituent Data Assimilation System (CoDAS) driven by meteorological fields from MERRA-2. M2-SCREAM assimilates hydrochloric acid (HCl), nitric acid (HNO3), stratospheric water vapor (H2O), nitrous oxide (N2O) and ozone with a system equipped with a version of the GEOS general circulation model and a stratospheric chemistry model, StratChem. Assimilated fields are provided globally at 0.5° by 0.625° resolution at three-hourly frequencies from 2004/09/01 to 2021. Assimilation uncertainties for each of the assimilated constituents are calculated from the CoDAS statistical output (Wargan et al., 2022) and provided as global full-resolution three-dimensional monthly files.

Krzysztof Wargan↗

Impact of Climatic Variability on Atmospheric Mass Distribution and GRACE-Derived Gravity Fields

During the period we calculated the atmospheric data sets related to its mass and angular momentum distribution. For mass, we determined the various harmonics from the NCEP-NCAR reanalysis, especially the low-order harmonics that are useful in studying the gravitation distribution as will be determined from the GRACE mission. Atmospheric mass is also related to the atmospheric loading on the solid Earth; we cooperated with scientists who needed the atmospheric mass information for understanding its contributions to the overall loading, necessary for vertical and horizontal coordinate estimation. We calculated atmospheric angular momentum from the NCEP-NCAR reanalyses and 4 operational meteorological centers, based on the motion (wind) terms and the mass (surface pressure) terms. These are associated with motions of the planet, including its axial component causing changes in the length of day, more related to the winds, and the equatorial component related to motions of the pole, more related to the mass. Tasks related to the ocean mass and angular momentum were added to the project as well. For these we have noted the ocean impact on motions of the pole as well as the torque mechanisms that relate the transfer of angular momentum between oceans and solid earth. The activities of the project may be summarized in the following first manuscript written in December 2002, for a symposium that Dr. Salstein attended on Geodynamics. We have continued to assess ocean angular momentum (OAM) quantities derived from bottom pressure and velocity fields estimated with our finite-difference barotropic (single layer) model. Three years of output (1993-95) from a run without any data constraints was compared to output from a corresponding run that was constrained by altimeter data using a Kalman filter and smoother scheme. Respective OAM time series were combined with corresponding atmospheric series and compared to observed polar motion. The constrained OAM series provided slightly better variance reduction than the unconstrained series. Analysis provided a check on the estimation scheme and pointed to further work to improve the determination of OAM using this method. A significant effort was also devoted to quantifying effects of uncertainties in high frequency winds on the mean and seasonal momentum exchange between atmosphere and oceans.

Salstein, David A.↗

Assessing the large-scale drivers of precipitation in the northeastern United States via linear orthogonal decomposition

This study examines the linear orthogonal modes associated with monthly precipitation in the northeastern United States, from CESM1 LENS (35 ensemble members, 1979–2005) and two reanalysis datasets (ERA5, 1950–2018 and NOAA-CIRES-DOE 20CRv3, 1950–2015). Calendar months are aggregated together, and any linear trends in data are removed. Using region-averaged precipitation anomaly time series and monthly anomalies for several global 2D atmospheric fields, a linear orthogonal decomposition method is implemented to iteratively extract time series (based on field and geographic location) of absolute maximum correlation. Linear modes associated with this method are then projected onto the full set of 2D fields to provide physical insight into the mechanisms involved in generating precipitation. In this region, the first mode is associated with vapor transport from the Atlantic seaboard, the second mode is characterized by westward vapor transport associated with extratropical cyclones, and the third mode captures vapor transport from the Gulf of Mexico during the fall and winter. However, the third mode is less robust in the spring and summer. Results are generally consistent across the datasets, and applying multiple linear regression with the linear modes to predict the precipitation anomalies produces R-squared values of around 0.54–0.65 for CESM1 LENS, and around 0.58–0.88 for reanalysis, with the lowest values generally in the spring and late summer. The influence of low-frequency climate variability on the modes is considered for CESM1 LENS, and the modes in late winter can be predicted with some success via a combination of several, prominent large-scale teleconnection patterns.

54 ENVIRONMENTAL SCIENCES↗

High-Resolution Fire Weather Index Data for the Conterminous US (1980–2099), Version 1

This dataset presents a suite of high-resolution fire weather index datasets calculated from observation (gridMet, Livneh, Daymet V4), reanalysis (AgERA5), downscaled hydro-climate projections over the conterminous United States (CONUS) based on multiple selected Global Climate Models (GCMs) from the Coupled Model Intercomparison Project Phase 6 (CMIP6). Aside from the daily FWI datasets, we also include a set of FWI extreme indicators at annual, seasonal, and monthly scales, including 1) fwixx: maximum FWI; 2) fwisa: mean FWI. For annual fwisa, it refers to the season with the maximum seasonal average; 3) fwils: Length of fire season over a specified period, where fire season is defined as the days exceeding the median value of the normalized FWI during the reference period (1980-1984); 4) fwixd: Number of extreme fire weather days over a specified period, where extreme day is defined as the day with FWI > the 95th percentile of the FWI during the reference period (1980-1984). All FWI datasets cover 1980-2020 baseline and the model simulated products including the downscaled products additionally include 2021-2099 near-future periods under the high-end (SSP585) emission scenario.

54 ENVIRONMENTAL SCIENCES↗

Dead Fuel Moisture Content Reanalysis Dataset for California (2000–2020)

This study presents a novel reanalysis dataset of dead fuel moisture content (DFMC) across California from 2000 to 2020 at a 2 km resolution. Utilizing a data assimilation system that integrates a simplified time-lag fuel moisture model with 10-h fuel moisture observations from remote automated weather stations (RAWS) allowed predictions of 10-h fuel moisture content by our method with a mean absolute error of 0.03 g/g compared to the widely used Nelson model, with a mean absolute error prediction of 0.05 g/g. For context, the values of DFMC in California are commonly between 0.05 g/g and 0.30 g/g. The presented product provides gridded hourly moisture estimates for 1-h, 10-h, 100-h, and 1000-h fuels, essential for analyzing historical fire activity and understanding climatological trends. The methodology presented here demonstrates significant advancements in the accuracy and robustness of fuel moisture estimates, which are critical for fire forecasting and management.

Farguell, Angel (ORCID:000000032395220X)↗

Using automated machine learning for the upscaling of gross primary productivity

Estimating gross primary productivity (GPP) over space and time is fundamental for understanding the response of the terrestrial biosphere to climate change. Eddy covariance flux towers provide in situ estimates of GPP at the ecosystem scale, but their sparse geographical distribution limits larger-scale inference. Machine learning (ML) techniques have been used to address this problem by extrapolating local GPP measurements over space using satellite remote sensing data. However, the accuracy of the regression model can be affected by uncertainties introduced by model selection, parameterization, and choice of explanatory features, among others. Recent advances in automated ML (AutoML) provide a novel automated way to select and synthesize different ML models. In this work, we explore the potential of AutoML by training three major AutoML frameworks on eddy covariance measurements of GPP at 243 globally distributed sites. We compared their ability to predict GPP and its spatial and temporal variability based on different sets of remote sensing explanatory variables. Explanatory variables from only Moderate Resolution Imaging Spectroradiometer (MODIS) surface reflectance data and photosynthetically active radiation explained over 70 % of the monthly variability in GPP, while satellite-derived proxies for canopy structure, photosynthetic activity, environmental stressors, and meteorological variables from reanalysis (ERA5-Land) further improved the frameworks' predictive ability. We found that the AutoML framework Auto-sklearn consistently outperformed other AutoML frameworks as well as a classical random forest regressor in predicting GPP but with small performance differences, reaching an r 2 of up to 0.75. We deployed the best-performing framework to generate global wall-to-wall maps highlighting GPP patterns in good agreement with satellite-derived reference data. This research benchmarks the application of AutoML in GPP estimation and assesses its potential and limitations in quantifying global photosynthetic activity.

54 ENVIRONMENTAL SCIENCES↗

Insights From Dayflow: A Historical Streamflow Reanalysis Dataset for the Conterminous United States

Abstract Reconstructed historical streamflow time series can supplement limited streamflow gauge observations. However, there are common challenges of typical modeling approaches: process‐based hydrologic models can be data/computation‐intensive, and statistics‐based models can be region/stream‐specific. Here we present a nationally scalable modeling framework integrating the simulated runoff from the Variable Infiltration Capacity (VIC) model with the Routing Application for Parallel computatIon of Discharge (RAPID) routing model leveraging high‐performance computing. We demonstrate an efficient method of assimilating streamflow at US Geological Survey (USGS) streamflow monitoring sites using a simple hierarchical approach in the VIC‐RAPID framework. The result is a reconstructed 36‐year (1980–2015) daily and monthly streamflow dataset (Dayflow) at ∼2.7 million NHDPlusV2 stream reaches in the conterminous US (CONUS). We perform a comprehensive evaluation at 7,526 USGS sites and characterize their error statistics. The results demonstrate that 49% of the USGS sites demonstrate Kling–Gupta Efficiency (KGE) > 0.5 and 58% of the sites show percentage bias within ±20% for the daily naturalized streamflow. Streamflow data assimilation across CONUS shows an overall improvement over naturalized streamflow, notably in the western semiarid‐to‐arid regions. Comparison to other national and global streamflow reanalysis datasets such as the National Water Model and Global Reach‐scale A priori Discharge Estimates for SWOT demonstrates improved KGE, reduced bias, and directions for Dayflow improvements. Investigations of error statistics with key hydrologic, hydroclimatic, and geomorphologic basin characteristics reveal region‐specific patterns which may help improve future framework applications. Overall, Dayflow may enable a better understanding of hydrologic conditions in a changing environment, especially in locations currently not represented by streamflow monitoring networks.

54 ENVIRONMENTAL SCIENCES↗

Unified 0.25-degree gridded infrastructure-critical extreme weather for the United States from 1979 to 2100

Extreme weather events can severely disrupt critical infrastructure, triggering cascading effects on power, transportation, and essential services. However, standard weather and climate datasets often lack specialized variables necessary for hazard assessments. We present a unified dataset of infrastructure-critical weather and climate variables across the United States at 0.25° resolution, covering daily or sub-daily intervals from 1979 to 2100. The dataset includes temperature, dew point, wind gusts, precipitation partitioned by rain, snow, and freezing rain or ice pellets, lightning, and wildfire metrics. Historical conditions (1979-2023) are synthesized from observations and reanalysis products, while future projections are derived from 14 CMIP6 global climate models (historical, SSP245, and SSP585 experiments). Physically based and data-driven methods are used to estimate variables not directly provided by existing models. By integrating these variables into a single unified dataset, we enable consistent, high-resolution assessments of weather-related infrastructure risks across past and future periods, supporting wide-ranging applications in energy, transportation, water resources, emergency management, and beyond.

Climate and Earth system modelling↗

Understanding the compound flood risk along the coast of the contiguous United States

Abstract. Compound flooding is a type of flood event caused by multiple flood drivers. The associated risk has usually been assessed using statistics-based analyses or hydrodynamics-based numerical models. This study proposes a compound flood (CF) risk assessment (CFRA) framework for coastal regions in the contiguous United States (CONUS). In this framework, a large-scale river model is coupled with a global ocean reanalysis dataset to (a) evaluate the CF exposure related to the coastal backwater effects on river basins, and (b) generate spatially distributed data for analyzing the CF hazard using a bivariate statistical model of river discharge and storm surge. The two kinds of risk are also combined to achieve a holistic understanding of the continental-scale CF risk. The estimated CF risk shows remarkable inter- and intra-basin variabilities along the CONUS coast with more variabilities in the CF hazard over the US west and Gulf coastal basins. Different risk assessment methods present significantly different patterns in a few key regions such as the San Francisco Bay area, the lower Mississippi River, and Puget Sound. Our results highlight the need to weigh different CF risk measures and avoid using single statistics-based or hydrodynamics-based CFRAs. Uncertainty sources in these CFRAs include the use of gauge observations, which cannot account for the flow physics or resolve the spatial variability of risks, and underestimations of the flood extremes and the dependence of CF drivers in large-scale models, highlighting the importance of understanding the CF risks for developing a more robust CFRA.

54 ENVIRONMENTAL SCIENCES↗

Dispersion of the Volcanic Sulfate Cloud from the Mount Pinatubo Eruption

We simulate the transport of the volcanic cloud from the 1991 eruption of Mount Pinatubo with the GEOS-5 general circulation model. Our simulations are in good agreement with observational data. We tested the importance of initial condition corresponding to the specific meteorological situation at the time of the eruption by employing reanalysis from MERRA. We found no significant difference in the transport of the cloud. We show how the inclusion of the interaction between volcanic sulfate aerosol and radiation is essential for a reliable simulation of the transport of the volcanic cloud. The absorption of long wave radiation by the volcanic sulfate induces a rising of the volcanic cloud up to the middle stratosphere, combined with divergent motion from the latitude of the eruption to the tropics. Our simulations indicate that the cloud diffuses to the northern hemisphere through a lower stratospheric pathway, and to mid- and high latitudes of the southern hemisphere through a middle stratospheric pathway, centered at about 30 hPa. The direction of the middle stratospheric pathway depends on the season. We did not detect any significant change of the mixing between tropics and mid- and high latitudes in the southern hemisphere.

Aquila, Valentina↗

Covariability in the Monthly Mean Convective and Radiative Diurnal Cycles in the Amazon

The diurnal cycle of convective clouds greatly influences the radiative energy balance in convectively active regions of Earth, through both direct presence, and the production of anvil and stratiform clouds. Previous studies show that the frequency and properties of convective clouds can vary on monthly timescales as a result of variability in the monthly mean atmospheric state. Furthermore, the radiative budget in convectively active regions also varies by up to 7 Wm-2 in convectively active regions. These facts suggest that convective clouds connect atmospheric state variability and radiation variability beyond clear sky effects alone. Previous research has identified monthly covariability between the diurnal cycle of CERES-observed top-of-atmosphere radiative fluxes and multiple atmospheric state variables from reanalysis over the Amazon region. ASVs that enhance (reduce) deep convection, such as CAPE (LTS), tend to shift the daily OLR and cloud albedo maxima earlier (later) in the day by 2-3 hr. We first test the analysis method using multiple reanalysis products for both the dry and wet seasons to further investigate the robustness of the preliminary results. We then use CloudSat data as an independent cloud observing system to further evaluate the relationships of cloud properties to variability in radiation and atmospheric states. While CERES can decompose OLR variability into clear sky and cloud effects, it cannot determine what variability in cloud properties lead to variability in the radiative cloud effects. Cloud frequency, cloud top height, and cloud microphysics all contribute to the cloud radiative effect, all of which are observable by CloudSat. In addition, CloudSat can also observe the presence and variability of deep convective cores responsible for the production of anvil clouds. We use these capabilities to determine the covariability of convective cloud properties and the radiative diurnal cycle.

Dodson, Jason B.↗

Understanding Arctic Surface Temperature Differences in Reanalyses

Reanalyses in the Arctic are widely used for model evaluation and for understanding contemporary climate change. Nevertheless, differences among reanalyses in fundamental meteorological variables including surface air temperature are large. For example, the 1980-2009 mean surface air temperature for the north polar cap (70ÂdegN-90ÂdegN) among global reanalyses span a range of 2.4 K, which approximates the average warming trend from these reanalyses over the 30-year period of 2.1 K. Understanding these differences requires evaluation over the three principal surface domains of the Arctic: glaciated land, the unglaciated terrestrial surface, and sea ice/ocean. An examination is conducted of contemporary global reanalyses of the ECMWF Interim project, NASA MERRA, MERRA-2, JRA-55, and NOAA CFSR using available in situ data and assessments of the surface energy budget. Overly-simplistic representations of the Greenland Ice Sheet surface are found to be associated with local warm air temperature biases in winter. A review of progress made in the development of the MERRA-2 land-ice representation is presented. Large uncertainty is also found in temperatures over the Arctic tundra and boreal forest zone. But a key focus of temperature differences for northern high latitudes is the Arctic Ocean. Near-surface air temperature differences over the Arctic Ocean are found to be related to discrepancies in sea ice and sea surface temperature boundary data, which are severely compromised in current reanalyses. Issues with the modeled representation of sea ice cover are an additional factor in reanalysis temperature trends. Differences in the representation of the surface energy budget among the various reanalyses are also reviewed.

Surface Temperatur↗

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↗