Search NASA⌕ Search

SEARCH · Search NASA

Results for “model reanalysis”

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 397 records · Page 22

Satellite-based Assessment of Global Warm Cloud Properties Associated with Aerosols, Atmospheric Stability, and Diurnal Cycle

This study examines variability in marine low cloud properties derived from semi-global observations by the Tropical Rainfall Measuring Mission (TRMM) satellite, as linked to the aerosol index (AI) and lower-tropospheric stability (LTS). AI is derived from the Moderate Resolution Imaging Spectroradiometer (Terra MODIS) sensor and the Goddard Chemistry Aerosol Radiation and Transportation (GOCART) model, and is used to represent column-integrated aerosol concentrations. LTS is derived from the NCEP/NCAR reanalysis, and represents the background thermodynamic environment in which the clouds form. Global statistics reveal that cloud droplet size tends to be smallest in polluted (high-AI) and strong inversion (high-LTS) environments. Statistical quantification shows that cloud droplet size is better correlated with AI than it is with LTS. Simultaneously, the cloud liquid water path (CLWP) tends to decrease as AI increases. This correlation does not support the hypothesis or assumption that constant or increased CLWP is associated with high aerosol concentrations. Global variability in corrected cloud albedo (CCA), the product of cloud optical depth and cloud fraction, is very well explained by LTS, while both AI and LTS are needed to explain local variability in CCA. Most of the local correlations between AI and cloud properties are similar to the results from the global statistics, while weak anomalous aerosol-cloud correlations appear locally in the regions where simultaneous high (low) AI and low (high) LTS compensate each other. Daytime diurnal cycles explain additional variability in cloud properties. CCA has the largest diurnal cycle in high-LTS regions. Cloud droplet size and CLWP have weak diurnal cycles that differ between clean and polluted environments. The combined results suggest that investigations of marine low cloud radiative forcing and its relationship to hypothesized aerosol indirect effects must consider the combined effects of aerosols, thermodynamics, and the diurnal cycle.

Matsui, Toshihisa↗

Principle Component Analysis of the Evolution of the Saharan Air Layer and Dust Transport: Comparisons between a Model Simulation and MODIS Retrievals

The onset and evolution of Saharan Air Layer (SAL) episodes during June-September 2002 are diagnosed by applying principal component analysis to the NCEP reanalysis temperature anomalies at 850 hPa, where the largest SAL-induced temperature anomalies are located. The first principal component (PC) represents the onset of SAL episodes, which are associated with large warm anomalies located at the west coast of Africa. The second PC represents two opposite phases of the evolution of the SAL. The positive phase of the second PC corresponds to the southwestward extension of the warm anomalies into the tropical-subtropical North Atlantic Ocean, and the negative phase corresponds to the northwestward extension into the subtropical to mid-latitude North Atlantic Ocean and the southwest Europe. A dust transport model (CARMA) and the MODIS retrievals are used to study the associated effects on dust distribution and deposition. The positive (negative) phase of the second PC corresponds to a strengthening (weakening) of the offshore flows in the lower troposphere around 10deg - 20degN, causing more (less) dust being transported along the tropical to subtropical North Atlantic Ocean. The variation of the offshore flow indicates that the subseasonal variation of African Easterly Jet is associated with the evolution of the SAL. Significant correlation is found between the second PC time series and the daily West African monsoon index, implying a dynamical linkage between West African monsoon and the evolution of the SAL and Saharan dust transport.

Wong, S.↗

Non-seasonal Gravitational Effects of Interannual Meteorological Oscillations in Atmosphere and Ocean

Via mass transport, the interannual meteorological oscillations (ENSO, NAO, etc.) produce slight variations in the global gravitational field, which can be observed by satellite laser ranging (SLR) and anticipated space missions such as GRACE. The meteorological oscillation modes are identified using empirical orthogonal function and principal component decomposition using de-seasoned surface fields produced from the 40-years of NCEP reanalysis. This decompositional analysis requires an accounting for issues associated with area-weighting and non-zero mean. The atmospheric contributions to gravitational variation, mode by mode and in combinations, are then computed via regional integration. To address the corresponding oceanic contributions, we examine the output of the numerical POCM ocean circulation model (1992-1998) in two ways: ocean bottom pressure field and ocean surface topography minus steric effect. The resulting lowest-degree Stokes coefficients for both NCEP reanalysis and POCM are then compared with space geodetic SLR observations, to identify the importance of each meteorological oscillations in gravitational variation signals.

Chao, Benjamin F.↗

Comparison of TOPEX/Poseidon Sea Level and Linear Model Results forced by Various Wind Products for the Tropical Pacific

The goal of this paper is to compare TOPEX/Posaidon (T/P) sea level with sea level results from linear ocean model experiments forced by several different wind products for the tropical Pacific. During the period of this study (October 1992 - October 1995), available wind products include satellite winds from the ERS-1 scatterometer product of [HALP 97] and the passive microwave analysis of SSMI winds produced using the variational analysis method (VAM) of [ATLA 91]. In addition, atmospheric GCM winds from the NCEP reanalysis [KALN 96], ECMWF analysis [ECMW94], and the Goddard EOS-1 (GEOS-1) reanalysis experiment [SCHU 93] are available for comparison. The observed ship wind analysis of FSU [STRI 92] is also included in this study. The linear model of [CANE 84] is used as a transfer function to test the quality of each of these wind products for the tropical Pacific. The various wind products are judged by comparing the wind-forced model sea level results against the T/P sea level anomalies. Correlation and RMS difference maps show how well each wind product does in reproducing the T/P sea level signal. These results are summarized in a table showing area average correlations and RMS differences. The large-scale low-frequency temporal signal is reproduced by all of the wind products, However, significant differences exist in both amplitude and phase on regional scales. In general, the model results forced by satellite winds do a better job reproducing the T/P signal (i.e. have a higher average correlation and lower RMS difference) than the results forced by atmospheric model winds.

Hackert, Eric C.↗

Diagnostic Studies with GLA Fields

Assessments of the NASA Goddard Earth Observing System-1 Data Assimilation System(GEOS-1 DAS), regarding heating rates, energetics, and angular momentum quantities were made. These diagnostics can be viewed as measures of climate variability. Comparisons with the NOAA/NCEP reanalysis system of momentum and energetics diagnostics are included. Water vapor and angular momentum are diagnosed in many models, including those of NASA, as part of the Atmospheric Model Intercomparison Project. 'Me GEOS-I and NOAA/NCEP global atmospheric angular momentum values are coherent on time scales down to about three days. Furthermore, they agree with the series of Earth angular momentum, as measured by tiny fluctuations in the rotation rate of the Earth, as variations in the length of day. The torques that effect such changes in atmospheric and Earth momentum are dominated by the influence of particular mountain systems, including the Rockies, Himalayas, and Andes, upon mountain torques on time scales shorter than about two weeks. Other project areas included collaboration with Goddard Space Flight Center to examine the impact of mountainous areas and the treatments of parameterizations on diagnoses of the atmosphere. Relevant preprints are included herein.

Salstein, David A.↗

Diagnostic Studies With GLA Fields

Assessments of the NASA Goddard Earth Observing System-1 Data Assimilation System (GEOS-1 DAS) regarding heating rates, energetics and angular momentum quantities were made. These diagnostics can be viewed as measures of climate variability. Comparisons with the NOAA/NCEP reanalysis system of momentum and energetics diagnostics are included. Water vapor and angular momentum are diagnosed in many models, including those of NASA, as part of the Atmospheric Model Intercomparison Project. Relevant preprints are included herein.

Salstein, David A.↗

Evolution of Tropical and Extratropical Precipitation Anomalies During the 1997 to 1999 ENSO Cycle

The 1997-1999 ENSO period was very powerful, but also well observed. Multiple satellite rainfall estimates combined with gauge observations allow for a quantitative analysis of precipitation anomalies in the tropics and elsewhere accompanying the 1997-99 ENSO cycle. An examination of the evolution of the El Nino and accompanying precipitation anomalies revealed that a dry Maritime Continent preceded the formation of positive SST anomalies in the eastern Pacific Ocean. 30-60 day oscillations in the winter of 1996/97 may have contributed to this lag relationship. Furthermore, westerly wind burst events may have maintained the drought over the Maritime Continent. The warming of the equatorial Pacific was then followed by an increase in convection. A rapid transition from El Nino to La Nina occurred in May 1998, but as early as October-November 1997 precipitation indices captured substantial changes in Pacific rainfall anomalies. The global precipitation patterns for this event were in good agreement with the strong consistent ENSO-related precipitation signals identified in earlier studies. Differences included a shift in precipitation anomalies over Africa during the 1997-98 El Nino and unusually wet conditions over northeast Australia during the later stages of the El Nino. Also, the typically wet region in the north tropical Pacific was mostly dry during the 1998-99 La Nina. Reanalysis precipitation was compared to observations during this time period and substantial differences were noted. In particular, the model had a bias towards positive precipitation anomalies and the magnitudes of the anomalies in the equatorial Pacific were small compared to the observations. Also, the evolution of the precipitation field, including the drying of the Maritime Continent and eastward progression of rainfall in the equatorial Pacific was less pronounced for the model compared to the observations.

Curtis, Scott↗

Predictability of the 1997 and 1998 South Asian Summer Monsoons on the Intraseasonal Time Scale Based on 10 AMIP2 Model Runs

Predictability of the 1997 and 1998 South Asian summer monsoons is examined using National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) reanalyses, and 100 two-year simulations with ten different Atmospheric General Circulation Models (AGCMs) with prescribed sea surface temperature (SST). We focus on the intraseasonal variations of the south Asian summer monsoon associated with the Madden-Julian Oscillation (MJO). The NCEP/NCAR reanalysis shows a clear coupling between SST anomalies and upper level velocity potential anomalies associated with the MJO. We analyze several MJO events that developed during the 1997 and 1998 focusing of the coupling with the SST. The same analysis is carried out for the model simulations. Remarkably, the ensemble mean of the two-year AGCM simulations show a signature of the observed MJO events. The ensemble mean simulated MJO events are approximately in phase with the observed events, although they are weaker, the period of oscillation is somewhat longer, and their onset is delayed by about ten days compared with the observations. Details of the analysis and comparisons among the ten AMIP2 (Atmospheric Model Intercomparison Project) models will be presented in the conference.

Wu, Man Li C.↗

The Global Character of the Flux of Downward Longwave Radiation

Four different types of estimates of the surface downwelling longwave radiative flux (DLR) are reviewed. One group of estimates synthesizes global cloud, aerosol, and other information in a radiation model that is used to calculate fluxes. Because these synthesis fluxes have been assessed against observations, the global-mean values of these fluxes are deemed to be the most credible of the four different categories reviewed. The global, annual mean DLR lies between approximately 344 and 350 W/sq m with an error of approximately +/-10 W/sq m that arises mostly from the uncertainty in atmospheric state that governs the estimation of the clear-sky emission. The authors conclude that the DLR derived from global climate models are biased low by approximately 10 W/sq m and even larger differences are found with respect to reanalysis climate data. The DLR inferred from a surface energy balance closure is also substantially smaller that the range found from synthesis products suggesting that current depictions of surface energy balance also require revision. The effect of clouds on the DLR, largely facilitated by the new cloud base information from the CloudSat radar, is estimated to lie in the range from 24 to 34 W/sq m for the global cloud radiative effect (all-sky minus clear-sky DLR). This effect is strongly modulated by the underlying water vapor that gives rise to a maximum sensitivity of the DLR to cloud occurring in the colder drier regions of the planet. The bottom of atmosphere (BOA) cloud effect directly contrast the effect of clouds on the top of atmosphere (TOA) fluxes that is maximum in regions of deepest and coldest clouds in the moist tropics.

Stephens, Graeme L.↗

Reconstruction of Historical Surface Mass Balance, 1984–2017 from GreenTrACS Multi-offset Ground-penetrating Radar

We present continuous estimates of snow and firn density, layer depth and accumulation from a multi-channel, multi-offset, ground-penetrating radar traverse. Our method uses the electromagnetic velocity, estimated from waveform travel-times measured at common-midpoints between sources and receivers. Previously, common-midpoint radar experiments on ice sheets have been limited to point observations. We completed radar velocity analysis in the upper ~2 m to estimate the surface and average snow density of the Greenland Ice Sheet. We parameterized the Herron and Langway (1980) firn density and age model using the radar-derived snow density, radar-derived surface mass balance (2015–2017) and reanalysis-derived temperature data. We applied structure-oriented filtering to the radar image along constant age horizons and increased the depth at which horizons could be reliably interpreted. We reconstructed the historical instantaneous surface mass balance, which we averaged into annual and multidecadal products along a 78 km traverse for the period 1984–2017. We found good agreement between our physically constrained parameterization and a firn core collected from the dry snow accumulation zone, and gained insights into the spatial correlation of surface snow density.

glacier geophysics↗

Use of Forecast Atmosphere for Earth Entry, Descent, and Landing Modeling

Flight mechanics simulations are used to characterize the performance of Earth entry, descent, and landing vehicles. Atmospheric prediction models are often a key component of these simulations. Global atmospheric models are used for both engineering design due to their ability to define atmospheric properties over a large swath of locations and times as well as used for the availability of atmospheric uncertainties in the models that can be used in statistical performance analysis. However for day-of-flight operations, the use of forecast atmospheres based on more current measurements are beneficial to the accuracy of the performance prediction, including landing locations. Reanalysis of these forecasts can also be used for post-flight analysis, including trajectory reconstruction. This paper describes how forecast atmospheres can be beneficial for Earth entry, descent, and landing analysis and ways these models can be implemented in simulations.

Soumyo Dutta↗

Basin-Scale River Runoff Estimation from GRACE Gravity Satellites, Climate Models and In Situ Observations: a Case Study in the Amazon Basin

River runoff is estimated as a water budget residual using Gravity Recovery and Climate Experiment (GRACE) terrestrial water storage time series, ERA5 reanalysis data, and precipitation observations for January 2003 through December 2015 for the Obidos upstream drainage basin and for the entire Amazon basin. Estimated runoff based on the water budget agrees remarkably well with in situgauge observations at Obidos, especially at seasonal time scales, with nearly perfect phase agreementbut slightly larger seasonal amplitude. The discrepancy in the seasonal amplitude maybe attributed to underestimation of river gauge runoff during the wet season when water overflows the riverbanks. The ERA5 model appears to overestimate long-term mean evapotranspiration in the Amazon by ~2 cm/monthb ased on comparisons with precipitation and runoff observations. Using precipitation data based on satellites and gauge observations relative to gauge observations alone improved agreement between water budget runoff estimates and in situ runoff observations. Seasonal variations in ERA5 simulated runoff are about twice as large as those from in situobservations and show a large phase lag as well. Water budget based runoff for the entire Amazon (~ 7,200 km3averaged 29over the 13 yr period) is significantly larger than observed runoff (~5,700 km3) at notably larger than previous estimates for the entire Amazon. These differences may be partly related to submarine runoff from the Amazon basin that cannot be captured by surface gauges.

Jianli Chen↗

Using Two-Stream Theory to Capture Fluctuations of Satellite-Perceived TOA SW Radiances Reflected from Clouds Over Ocean

Shortwave (SW) fluxes estimated from broadband radiometry rely on empirically gathered and hemispherically resolved fields of outgoing top-of-atmosphere (TOA) radiances. This study aims to provide more accurate and precise fields of TOA SW radiances reflected from clouds over ocean by introducing a novel semiphysical model predicting radiances per narrow sun-observer geometry. This model was statistically trained using CERES-measured radiances paired with MODIS-retrieved cloud parameters as well as reanalysis-based geophysical parameters. By using radiative transfer approximations as a framework to ingest the above parameters, the new approach incorporates cloud-top effective radius and above-cloud water vapor in addition to traditionally used cloud optical depth, cloud fraction, cloud phase, and surface wind speed. A two-stream cloud albedo – serving to statistically incorporate cloud optical thickness and cloud-top effective radius – and Cox–Munk ocean reflectance were used to describe an albedo over each CERES footprint. Effective-radius-dependent asymmetry parameters were obtained empirically and separately for each viewing-illumination geometry. A simple equation of radiative transfer, with this albedo and attenuating above-cloud water vapor as inputs, was used in its log-linear form to allow for statistical optimization. We identified the two-stream functional form that minimized radiance residuals calculated against CERES observations and outperformed the state-of-the-art approach for most observer geometries outside the sun-glint and solar zenith angles between 20 and 70∘, reducing the median SD of radiance residuals per solar geometry by up to 13.2 % for liquid clouds, 1.9 % for ice clouds, and 35.8 % for footprints containing both cloud phases. Geometries affected by sun glint (constituting between 10 % and 1 % of the discretized upward hemisphere for solar zenith angles of 20 and 70∘, respectively), however, often showed weaker performance when handled with the new approach and had increased residuals by as much as 60 % compared to the state-of-the-art approach. Overall, uncertainties were reduced for liquid-phase and mixed-phase footprints by 5.76 % and 10.81 %, respectively, while uncertainties for ice-phase footprints increased by 0.34 %. Tested for a variety of scenes, we further demonstrated the plausibility of scene-wise predicted radiance fields. This new approach may prove useful when employed in angular distribution models and may result in improved flux estimates, in particular dealing with clouds characterized by small or large droplet/crystal sizes.

Two- stream theory↗

Improving the Timeliness of Winter Wheat Production Forecast in the United States of America, Ukraine and China Using MODIS Data and NCAR Growing Degree Day Information

Wheat is themost important cereal crop traded on international markets and winter wheat constitutes approximately 80% of global wheat production. Thus, accurate and timely production forecasts are critical for making informed agricultural policies and investments, as well as increasing market efficiency and stability. Becker-Reshef et al. (2010) developed an empirical generalized model for forecasting winter wheat production. Their approach combined BRDF-corrected daily surface reflectance from Moderate resolution Imaging Spectroradiometer (MODIS) Climate Modeling Grid (CMG) with detailed official crop statistics and crop typemasks. It is based on the relationship between the Normalized Difference Vegetation Index (NDVI) at the peak of the growing season, percent wheat within the CMG pixel (area within the CMG pixel occupied by wheat), and the final yields. This method predicts the yield approximately one month to six weeks prior to harvest. In this study, we include Growing Degree Day (GDD) information extracted from NCEP/NCAR reanalysis data in order to improve the winter wheat production forecast by increasing the timeliness of the forecasts while conserving the accuracy of the original model. We apply this modified model to three major wheat-producing countries: the Unites States (US), Ukraine and China from 2001 to 2012. We show that a reliable forecast can be made between one month to a month and a half prior to the peak NDVI (meaning two months to two and a half months prior to harvest), while conserving an accuracy of 10% in the production forecast.

Production↗

Equatorial Annual Oscillation with QBO-driven 5-year Modulation in NCEP Data

An analysis is presented of the zonal wind and temperature variations supplied by the National Center for Environmental Prediction (NCEP), which have been assimilated in the Reanalysis and the Climate Prediction Center (CCP) data sets. The derived zonal-mean variations are employed. Stimulated by modeling studies, the data are separated into the hemispherically symmetric and anti-symmetric components, and spectral analysis is applied to study the annual 12-month oscillation and Quasi-biennial Oscillation (QBO). For data samples that cover as much as 40 years, the results reveal a pronounced 5-year modulation of the symmetric AO in the lower stratosphere, which is confined to equatorial latitudes. This modulation is also inferred for the temperature variations but extends to high latitudes, qualitatively consistent with published model results. A comparison between different data samples indicates that the signature of the 5-year oscillation is larger when the QBO of 30 months is more pronounced. Thus there is circumstantial evidence that this periodicity of the QBO is involved in generating the oscillation. The spectral analysis shows that there is a weak anti-symmetric 5-year oscillation in the zonal winds, which could interact with the large antisymmetric A0 to produce the modulation of the symmetric AO as was shown in earlier modeling studies. According to these studies, the 30-month QBO tends to be synchronized by the equatorial Semi-annual Oscillation (SAO), and this would explain why the inferred 5-year modulation is observed to persist and is phase locked over several cycles.

Mayr, H. G.↗

Toward The Development of Hailstorm Climatologies Derived From Reanalyses and Infared/Passive Microwave Satellite Imagers

Geostationary satellite imagers, such as those of the Geostationary Operational Environmental Satellite (GOES) and Meteosat series, provide both historical and near-real-time observations of cloud top patterns that are commonly associated with severe convection. Environmental conditions favorable for severe weather are thought to be represented well by reanalyses. Predicting exactly where convection and costly storm hazards like hail will occur using models or satellite imagery alone, however, is extremely challenging. The multivariate combination of satellite-observed cloud patterns with reanalysis environmental parameters, linked to United States Next Generation Weather Radar- (NEXRAD-) estimated Maximum Expected Size of Hail (MESH) using a deep neural network (DNN), enables estimation of potentially severe hail likelihood for any observed storm cell. These estimates are specifically designed to make hail likelihood distinctions based on satellite-indicated points of deep convection within environments favorable for storm development. We seek an approach that can be used to estimate climatological hailstorm frequency and risk throughout the historical satellite data record. This presentation demonstrates that statistical distributions of convective parameters from satellite and reanalysis show separation between non-severe/severe hailstorm classes for predictors including overshooting cloud top temperature and area characteristics, convective available potential energy, vertical wind shear, 500 hPa temperature, mid-level lapse rate, precipitable water, and convective inhibition. These complex, multivariate predictor relationships are exploited within a DNN to produce a hail likelihood metric with a critical success index of 0.504 and Heidke skill score of 0.403, which is exceptional among recent analogous hail studies. Furthermore, applications of the DNN to select case studies demonstrate good qualitative agreement between hail likelihood and MESH. These hail classifications are aggregated across an 11-year GOES-12/13 image database to derive a hail frequency and severity climatology, which denotes the Central Plains, the Midwest, and northwestern Mexico as being the most hail-prone regions within the domain studied. Opportunities for training and applying DNN-based hailstorm predictions to recently developed GOES-8/10/12/13/16 and Meteosat Second Generation convective storm detection and characterization climatologies over South America and South Africa, respectively, will also be presented.

Kristopher Bedka↗

The NASA Orbital Debris Engineering Model 3.1: Development, Verification, and Validation

The NASA Orbital Debris Program Office has developed the Orbital Debris Engineering Model (ORDEM) primarily as a tool for spacecraft designers and other users to understand the long-term risk of collisions with orbital debris. The newest version, ORDEM 3.1, incorporates the latest and highest fidelity datasets available to build and validate representative orbital debris populations encompassing low Earth orbit (LEO) to geosynchronous orbit (GEO) altitudes for the years 2016-2050. ORDEM 3.1 models fluxes for object sizes > 10 μm within or transiting LEO and > 10 cm in GEO. The deterministic portion of the populations in ORDEM 3.1 is based on the U.S. Space Surveillance Network (SSN) catalog, which provides coverage down to approximately 10 cm in LEO and 1 m in GEO. Observational datasets from radar, in situ, and optical sources provide a foundation from which the model populations are statistically extrapolated to smaller sizes and orbit regions that are not well-covered by the SSN catalog, yet may pose the greatest threat to operational spacecraft. Objects in LEO ranging from approximately 5 mm to 10 cm are modeled using observational data from ground-based radar, namely the Haystack Ultrawideband Satellite Imaging Radar (HUSIR – formerly known as Haystack). The LEO population smaller than approximately 3 mm in size is characterized based on a reanalysis of in situ data from impacts to the windows and radiators of the U.S. Space Transportation System orbiter vehicle, i.e., the Space Shuttle. Data from impacts on the Hubble Space Telescope are also used to validate the sub-millimeter model populations in LEO. Debris in GEO with sizes ranging from 10 cm to 1 m is modeled using optical measurement data from the Michigan Orbital DEbris Survey Telescope (MODEST). Specific, major debris-producing events, including the Fengyun-1C, Iridium 33, and Cosmos 2251 debris clouds, and unique populations, such as sodium-potassium droplets, have been re-examined and are modeled and added to the ORDEM environment separately. The debris environment greater than 1 mm is forecast using NASA’s LEO-to- GEO ENvironment Debris model (LEGEND). Future explosions of intact objects and collisions involving objects greater than 10 cm are assessed statistically, and the NASA Standard Satellite Breakup Model is used to generate fragments from these events. Fragments smaller than 10 cm are further differentiated based on material density categories, i.e., high-, medium-, and low-density, to better characterize the potential debris risk posed to spacecraft. The future projection of the sub-millimeter environment is computed using a special small-particle degradation model where small particles are created from intact spacecraft and rocket bodies. This work discusses the development, features, and capabilities of the ORDEM 3.1 model; the ne new data analyses used to build the model populations; and sample verification and validation results.

Matney, M.↗

How Accurate is Land/Ocean Moisture Transport Variability in Reanalyses?

Quantifying the global hydrological cycle and its variability across various time scales remains a challenge to the climate community. Direct measurements of evaporation (E), evapotranspiration (ET), and precipitation (P) are not feasible on a global scale, nor is the transport of water vapor over the global oceans and sparsely populated land areas. Expanding satellite data streams have enabled development of various water (and energy) flux products, complementing reanalyses and facilitating observationally constrained modeling. But the evolution of the global observing system has produced additional complications--improvements in satellite sensor resolution and accuracy have resulted in "epochs" of observational quasi-uniformity that can adversely affect reanalysis trends. In this work we focus on vertically integrated moisture flux convergence (VMFC) variations within the period 1979 - present integrated over global land. We show that VMFC in recent reanalyses (e.g. ERA-I, NASA MERRA, NOAA CFSR and JRA55) suffers from observing system changes, though differently in each product. Land Surface Models (LSMs) forced with observations-based precipitation, radiation and near-surface meteorology share closely the interannual P-ET variations of the reanalyses associated with ENSO events. (VMFC over land and P-ET estimates are equivalent quantities since atmospheric storage changes are small on these scales.) But the long-term LSM trend over the period since 1979 is approximately one-fourth that of the reanalyses. Additional reduced observation reanalyses assimilating only surface pressure and /or specifying seasurface temperature also have a much smaller trend in P-ET like the LSMs. We explore the regional manifestation of the reanalysis P-ET / VMFC problems, particularly over land. Both principal component analysis and a simple time series changepoint analysis highlight problems associated with data poor regions such as Equatorial Africa and, for one reanalysis, the Equatorial Andes region. Onset of the availability of passive microwave Special Sensor Microwave Imager (SSMI) moisture data in July 1987 and the transition from the Microwave Sounder Unit (MSU) to an advanced version (AMSU) have significant impacts on VMFC variability. Simple accounting for these errors of leading importance results in modified reanalysis VMFC estimates that agree much better with the LSM results. Regional details of the modified reanalysis VMFC and LSM P-ET are related to changes in Pacific Decadal Variability as manifest in SST changes after the late 1990s.

Robertson, F. R.↗