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)

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

Evaluation of the Analysis Influence on Transport in Reanalysis Regional Water Cycles

Regional water cycles of reanalyses do not follow theoretical assumptions applicable to pure simulated budgets. The data analysis changes the wind, temperature and moisture, perturbing the theoretical balance. Of course, the analysis is correcting the model forecast error, so that the state fields should be more aligned with observations. Recently, it has been reported that the moisture convergence over continental regions, even those with significant quantities of radiosonde profiles present, can produce long term values not consistent with theoretical bounds. Specifically, long averages over continents produce some regions of moisture divergence. This implies that the observational analysis leads to a source of water in the region. One such region is the Unite States Great Plains, which many radiosonde and lidar wind observations are assimilated. We will utilize a new ancillary data set from the MERRA reanalysis called the Gridded Innovations and Observations (GIO) which provides the assimilated observations on MERRA's native grid allowing more thorough consideration of their impact on regional and global climatology. Included with the GIO data are the observation minus forecast (OmF) and observation minus analysis (OmA). Using OmF and OmA, we can identify the bias of the analysis against each observing system and gain a better understanding of the observations that are controlling the regional analysis. In this study we will focus on the wind and moisture assimilation

Regional Climate and Variability of the Summertime Continental United States in Reanalyses

Understanding climate variability at regional scales is an important for research and societal needs. Atmospheric retrospective-analyses (or reanalyses) integrate multitudes of observing systems with numerical models to produce continuous data that include variables not easily observed, if at all. The breadth of variables as well as observational influence included in reanalyses make them ideal for investigating climate variability. In this paper, we assess NASA s Modern Era Retrospective-analysis for Research and Application (MERRA) regional variability in North America, specifically the United States, in conjunction with current satellite data reanalyses. Emphasis is placed on summertime precipitation because 1) it is a difficult parameter to capture in the most difficult season, 2) significant observational resources exist to benchmark comparisons, and 3) accurate assessment of precipitation variability is crucial to a multitude of sectors and applications. Likewise, we have also begun to evaluate surface air temperature. While precipitation biases are identified, year to year variability of the precipitation variations, in many cases, are quite reasonable. However, some spurious long term trends and sudden shifts in the time series are also identified. In surface air temperature, analysis of station observations provides ERA Interim a clear overall advantage. However, in a number of regions, all the reanalyses are quite comparable in variability and trend. In other regions, significant precipitation biases may occur, which has implications for the ancillary process data in a reanalysis, such as surface fluxes. We also characterize the reanalyses ability to capture variability related to ENSO. In general, the summertime variations of precipitation in the reanalyses are more highly correlated (positively) to ENSO (using ENSO34) than are the observations. The Northwestern US shows the largest positive correlations to ENSO34, and reanalyses agree with that, and correlate highly with the gauge observations there. Additional evaluation of the data assimilation using the MERRA Gridded Innovations and Observations (GIO) data, which consists of the assimilated observations as well as forecast an analysis error fields. This work represents an initial evaluation of using the MERRA reanalysis to study regional climate variations of the United States, and identifying its applied limitation

The Impact of the Evolving Satellite Data Record on Reanalysis Water and Energy Fluxes During the Past 30 Years

Retrospective analyses (reanalyses) use a fixed assimilation model to take diverse observations and synthesize consistent, time-dependent fields of state variables and fluxes (e.g. temperature, moisture, momentum, turbulent and radiative fluxes). Because they offer data sets of these quantities at regular space / time intervals, atmospheric reanalyses have become a mainstay of the climate community for diagnostic purposes and for driving offline ocean and land models. Of course, one weakness of these data sets is the susceptibility of the flux products to uncertainties because of shortcomings in parameterized model physics. Another issue, perhaps less appreciated, is the fact that the discreet changes in the evolving observational system, particularly from satellite sensors, may also introduce artifacts in the time series of quantities. In this paper we examine the ability of the NASA MERRA (Modern Era Retrospective Analysis for Research and Applications) and other recent reanalyses to determine variability in the climate system over the satellite record (~ the last 30 years). In particular we highlight the effect on reanalyses of discontinuities at the junctures of the onset of passive microwave imaging (Special Sensor Microwave Imager) in late 1987 as well as improved sounding and imaging with the Advanced Microwave Sounding Unit, AMSU-A, in 1998. We examine these data sets from two perspectives. The first is the ability to capture modes of variability that have coherent spatial structure (e.g. ENSO events and near-decadal coupling to SST changes) and how these modes are contained within trends in near global averages of key quantities. Secondly, we consider diagnostics that measure the consistency in energetic scaling in the hydrologic cycle, particularly the fractional changes in column-integrated water vapor versus precipitation as they are coupled to radiative flux constraints. These results will be discussed in the context of implications for science objectives and priorities of the NASA Energy and Water Cycle Study, NEWS

The Effect of Satellite Observing System Changes on MERRA Water and Energy Fluxes

Because reanalysis data sets offer state variables and fluxes at regular space / time intervals, atmospheric reanalyses have become a mainstay of the climate community for diagnostic purposes and for driving offline ocean and land models. Although one weakness of these data sets is the susceptibility of the flux products to uncertainties because of shortcomings in parameterized model physics, another issue, perhaps less appreciated, is the fact that continual but discreet changes in the evolving observational system, particularly from satellite sensors, may also introduce artifacts in the time series of quantities. In this paper we examine the ability of the NASA MERRA (Modern Era Retrospective Analysis for Research and Applications) and other recent reanalyses to determine variability in the climate system over the satellite record (approximately the last 30 years). In particular we highlight the effect on the reanalysis of discontinuities at the junctures of the onset of passive microwave imaging (Special Sensor Microwave Imager) in late 1987 as well as improved sounding and imaging with the Advanced Microwave Sounding Unit, AMSU-A, in 1998. We first examine MERRA fluxes from the perspective of how physical modes of variability (e.g. ENSO events, Pacific Decadal Variability) are contaminated by artificial step-like trends induced by the onset of new moisture data these two satellite observing systems. Secondly, we show how Redundancy Analysis, a statistical regression methodology, is effective in relating these artifact signals in the moisture and temperature analysis increments to their presence in the physical flux terms (e.g. precipitation, radiation). This procedure is shown to be effective greatly reducing the artificial trends in the flux quantities

Comparing Evaporative Sources of Terrestrial Precipitation and Their Extremes in MERRA Using Relative Entropy

A quasi-isentropic back trajectory scheme is applied to output from the Modern Era Retrospective-analysis for Research and Applications and a land-only replay with corrected precipitation to estimate surface evaporative sources of moisture supplying precipitation over every ice-free land location for the period 1979-2005. The evaporative source patterns for any location and time period are effectively two dimensional probability distributions. As such, the evaporative sources for extreme situations like droughts or wet intervals can be compared to the corresponding climatological distributions using the method of relative entropy. Significant differences are found to be common and widespread for droughts, but not wet periods, when monthly data are examined. At pentad temporal resolution, which is more able to isolate floods and situations of atmospheric rivers, values of relative entropy over North America are typically 50-400 larger than at monthly time scales. Significant differences suggest that moisture transport may be the key to precipitation extremes. Where evaporative sources do not change significantly, it implies other local causes may underlie the extreme events

Development of Gridded Innovations and Observations Supplement to MERRA-2

Atmospheric reanalysis have become an important source of data for weather and climate research, owing to the continuity of the data, but especially because of the multitude of observational data included (radiosondes, commercial aircraft, retrieved data products and radiances). However, the presence of assimilated observations can vary based on numerous factors, and so it is difficult or impossible for a researcher to say with any degree of certainty how many and what type of observations contributed to the reanalysis data they are using at any give point in time or space. For example, quality control, transmission interruptions, and station outages can occasionally affect data availability. While orbital paths can be known, drift in certain instruments and the large number of available instruments makes it challenging to know which satellite is observing any region at any point in the diurnal cycle. Furthermore, there is information from the statistics generated by the data assimilation that can help understand the model and the quality of the reanalysis. Typically, the assimilated observations and their innovations are in observation-space data formats and have not been made easily available to reanalysis users.A test data set has been developed to make the MERRA-2 assimilated observations available for rapid and general use, by simplifying the data format. The observations are binned to a grid similar as MERRA-2 and saved as netCDF. This data collection includes the mean and number of observations in the bin as well as its variance. The data will also include the innovations from the data assimilation, the forecast departure and the analysis increment, as well as bias correction (for satellite radiances). We refer to this proof-of-concept data as the MERRA-2 Gridded Innovations and Observations (GIO). In this paper, we present the data format and its strengths and limitations with some initial testing and validation of the methodology

Uncertainty in Tropical Ocean Latent Heat Flux Variability During the Last 25 Years

When averaged over the tropical oceans (30deg N/S), latent heat flux anomalies derived from passive microwave satellite measurements as well as reanalyses and climate models driven with specified seal-surface temperatures show considerable disagreement in their decadal trends. These estimates range from virtually no trend to values over 8.4 W/sq m decade. Satellite estimates also tend to have a larger interannual signal related to El Nino/Southern Oscillation (ENSO) events than do reanalyses or model simulations. An analysis of wind speed and humidity going into bulk aerodynamic calculations used to derive these fluxes reveals several error sources. Among these are apparent remaining intercalibration issues affecting passive microwave satellite 10 m wind speeds and systematic biases in retrieval of near-surface humidity. Likewise, reanalyses suffer from discontinuities in availability of assimilated data that affect near surface meteorological variables. The results strongly suggest that current latent heat flux trends are overestimated

Where Does the Irrigation Water Go? An Estimate of the Contribution of Irrigation to Precipitation Using MERRA

Irrigation is an important human activity that may impact local and regional climate, but current climate model simulations and data assimilation systems generally do not explicitly include it. The European Centre for Medium-Range Weather Forecasts (ECMWF) Interim Re-Analysis (ERA-Interim) shows more irrigation signal in surface evapotranspiration (ET) than the Modern-Era Retrospective Analysis for Research and Applications (MERRA) because ERA-Interim adjusts soil moisture according to the observed surface temperature and humidity while MERRA has no explicit consideration of irrigation at the surface. But, when compared with the results from a hydrological model with detailed considerations of agriculture, the ET from both reanalyses show large deficiencies in capturing the impact of irrigation. Here, a back-trajectory method is used to estimate the contribution of irrigation to precipitation over local and surrounding regions, using MERRA with observation-based corrections and added irrigation-caused ET increase from the hydrological model. Results show substantial contributions of irrigation to precipitation over heavily irrigated regions in Asia, but the precipitation increase is much less than the ET increase over most areas, indicating that irrigation could lead to water deficits over these regions. For the same increase in ET, precipitation increases are larger over wetter areas where convection is more easily triggered, but the percentage increase in precipitation is similar for different areas. There are substantial regional differences in the patterns of irrigation impact, but, for all the studied regions, the highest percentage contribution to precipitation is over local land