NOAA National Centers for Environmental Prediction

Abstract ERA‐40 is a re‐analysis of meteorological observations from September 1957 to August 2002 produced by the European Centre for Medium‐Range Weather Forecasts (ECMWF) in collaboration with many institutions. The observing system changed considerably over this re‐analysis period, with assimilable data provided by a succession of satellite‐borne instruments from the 1970s onwards, supplemented by increasing numbers of observations from aircraft, ocean‐buoys and other surface platforms, but with a declining number of radiosonde ascents since the late 1980s. The observations used in ERA‐40 were accumulated from many sources. The first part of this paper describes the data acquisition and the principal changes in data type and coverage over the period. It also describes the data assimilation system used for ERA‐40. This benefited from many of the changes introduced into operational forecasting since the mid‐1990s, when the systems used for the 15‐year ECMWF re‐analysis (ERA‐15) and the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) re‐analysis were implemented. Several of the improvements are discussed. General aspects of the production of the analyses are also summarized. A number of results indicative of the overall performance of the data assimilation system, and implicitly of the observing system, are presented and discussed. The comparison of background (short‐range) forecasts and analyses with observations, the consistency of the global mass budget, the magnitude of differences between analysis and background fields and the accuracy of medium‐range forecasts run from the ERA‐40 analyses are illustrated. Several results demonstrate the marked improvement that was made to the observing system for the southern hemisphere in the 1970s, particularly towards the end of the decade. In contrast, the synoptic quality of the analysis for the northern hemisphere is sufficient to provide forecasts that remain skilful well into the medium range for all years. Two particular problems are also examined: excessive precipitation over tropical oceans and a too strong Brewer‐Dobson circulation, both of which are pronounced in later years. Several other aspects of the quality of the re‐analyses revealed by monitoring and validation studies are summarized. Expectations that the ‘second‐generation’ ERA‐40 re‐analysis would provide products that are better than those from the firstgeneration ERA‐15 and NCEP/NCAR re‐analyses are found to have been met in most cases. © Royal Meteorological Society, 2005. The contributions of N. A. Rayner and R. W. Saunders are Crown copyright.

This powerful new land surface modeling system integrates data from advanced observing systems to support improved forecast model initialization and hydrometeorological investigations.

The Global Precipitation Climatology Project (GPCP) Version-2 Monthly Precipitation Analysis is described. This globally complete, monthly analysis of surface precipitation at 2.5 latitude 2.5 longitude resolution is available from January 1979 to the present. It is a merged analysis that incorporates precipitation estimates from low-orbit satellite microwave data, geosynchronous-orbit satellite infrared data, and surface rain gauge observations. The merging approach utilizes the higher accuracy of the low-orbit microwave observations to calibrate, or adjust, the more frequent geosynchronous infrared observations. The dataset is extended back into the premicrowave era (before mid-1987) by using infrared-only observations calibrated to the microwave-based analysis of the later years. The combined satellite-based product is adjusted by the rain gauge analysis. The dataset archive also contains the individual input fields, a combined satellite estimate, and error estimates for each field. This monthly analysis is the foundation for the GPCP suite of products, including those at finer temporal resolution. The 23-yr GPCP climatology is characterized, along with time and space variations of precipitation.

The NCEP Climate Forecast System Reanalysis (CFSR) was completed for the 31-yr period from 1979 to 2009, in January 2010. The CFSR was designed and executed as a global, high-resolution coupled atmosphere–ocean–land surface–sea ice system to provide the best estimate of the state of these coupled domains over this period. The current CFSR will be extended as an operational, real-time product into the future. New features of the CFSR include 1) coupling of the atmosphere and ocean during the generation of the 6-h guess field, 2) an interactive sea ice model, and 3) assimilation of satellite radiances by the Gridpoint Statistical Interpolation (GSI) scheme over the entire period. The CFSR global atmosphere resolution is ~38 km (T382) with 64 levels extending from the surface to 0.26 hPa. The global ocean's latitudinal spacing is 0.25° at the equator, extending to a global 0.5° beyond the tropics, with 40 levels to a depth of 4737 m. The global land surface model has four soil levels and the global sea ice model has three layers. The CFSR atmospheric model has observed variations in carbon dioxide (CO2) over the 1979–2009 period, together with changes in aerosols and other trace gases and solar variations. Most available in situ and satellite observations were included in the CFSR. Satellite observations were used in radiance form, rather than retrieved values, and were bias corrected with “spin up” runs at full resolution, taking into account variable CO2 concentrations. This procedure enabled the smooth transitions of the climate record resulting from evolutionary changes in the satellite observing system. CFSR atmospheric, oceanic, and land surface output products are available at an hourly time resolution and a horizontal resolution of 0.5° latitude × 0.5° longitude. The CFSR data will be distributed by the National Climatic Data Center (NCDC) and NCAR. This reanalysis will serve many purposes, including providing the basis for most of the NCEP Climate Prediction Center's operational climate products by defining the mean states of the atmosphere, ocean, land surface, and sea ice over the next 30-yr climate normal (1981–2010); providing initial conditions for historical forecasts that are required to calibrate operational NCEP climate forecasts (from week 2 to 9 months); and providing estimates and diagnoses of the Earth's climate state over the satellite data period for community climate research. Preliminary analysis of the CFSR output indicates a product that is far superior in most respects to the reanalysis of the mid-1990s. The previous NCEP–NCAR reanalyses have been among the most used NCEP products in history; there is every reason to believe the CFSR will supersede these older products both in scope and quality, because it is higher in time and space resolution, covers the atmosphere, ocean, sea ice, and land, and was executed in a coupled mode with a more modern data assimilation system and forecast model.

Abstract The Modern-Era Retrospective Analysis for Research and Applications (MERRA) was undertaken by NASA’s Global Modeling and Assimilation Office with two primary objectives: to place observations from NASA’s Earth Observing System satellites into a climate context and to improve upon the hydrologic cycle represented in earlier generations of reanalyses. Focusing on the satellite era, from 1979 to the present, MERRA has achieved its goals with significant improvements in precipitation and water vapor climatology. Here, a brief overview of the system and some aspects of its performance, including quality assessment diagnostics from innovation and residual statistics, is given. By comparing MERRA with other updated reanalyses [the interim version of the next ECMWF Re-Analysis (ERA-Interim) and the Climate Forecast System Reanalysis (CFSR)], advances made in this new generation of reanalyses, as well as remaining deficiencies, are identified. Although there is little difference between the new reanalyses in many aspects of climate variability, substantial differences remain in poorly constrained quantities such as precipitation and surface fluxes. These differences, due to variations both in the models and in the analysis techniques, are an important measure of the uncertainty in reanalysis products. It is also found that all reanalyses are still quite sensitive to observing system changes. Dealing with this sensitivity remains the most pressing challenge for the next generation of reanalyses. Production has now caught up to the current period and MERRA is being continued as a near-real-time climate analysis. The output is available online through the NASA Goddard Earth Sciences Data and Information Services Center (GES DISC).

A weekly 18 spatial resolution optimum interpolation (OI) sea surface temperature (SST) analysis has been produced at the National Oceanic and Atmospheric Administration (NOAA) using both in situ and satellite data from November 1981 to the present. The weekly product has been available since 1993 and is widely used for weather and climate monitoring and forecasting. Errors in the satellite bias correction and the sea ice to SST conversion algorithm are discussed, and then an improved version of the OI analysis is developed. The changes result in a modest reduction in the satellite bias that leaves small global residual biases of roughly 20.038C. The major improvement in the analysis occurs at high latitudes due to the new sea ice algorithm where local differences between the old and new analysis can exceed 18C. Comparisons with other SST products are needed to determine the consistency of the OI. These comparisons show that the differences among products occur on large time- and space scales with monthly rms differences exceeding 0.58C in some regions. These regions are primarily the mid- and high-latitude Southern Oceans and the Arctic where data are sparse, as well as high-gradient areas such as the Gulf Stream and Kuroshio where the gradients cannot be properly resolved on a 18 grid. In addition, globally averaged differences of roughly 0.058C occur among the products on decadal scales. These differences primarily arise from the same regions where the rms differences are large. However, smaller unexplained differences also occur in other regions of the midlatitude Northern Hemisphere where in situ data should be adequate. 1.

Gridded fields (analyses) of global monthly precipitation have been constructed on a 2.5° latitude–longitude grid for the 17-yr period from 1979 to 1995 by merging several kinds of information sources with different characteristics, including gauge observations, estimates inferred from a variety of satellite observations, and the NCEP–NCAR reanalysis. This new dataset, which the authors have named the CPC Merged Analysis of Precipitation (CMAP), contains precipitation distributions with full global coverage and improved quality compared to the individual data sources. Examinations showed no discontinuity during the 17-yr period, despite the different data sources used for the different subperiods. Comparisons of the CMAP with the merged analysis of Huffman et al. revealed remarkable agreements over the global land areas and over tropical and subtropical oceanic areas, with differences observed over extratropical oceanic areas. The 17-yr CMAP dataset is used to investigate the annual and interannual variability in large-scale precipitation. The mean distribution and the annual cycle in the 17-yr dataset exhibit reasonable agreement with existing long-term means except over the eastern tropical Pacific. The interannual variability associated with the El Niño-Southern Oscillation phenomenon resembles that found in previous studies, but with substantial additional details, particularly over the oceans. With complete global coverage, extended period and improved quality, the 17-yr dataset of the CMAP provides very useful information for climate analysis, numerical model validation, hydrological research, and many other applications. Further work is under way to improve the quality, extend the temporal coverage, and to refine the resolution of the merged analysis.

The National Centers for Environmental Prediction (NCEP) and National Center for Atmospheric Research (NCAR) have cooperated in a project (denoted "reanalysis") to produce a retroactive record of more than 50 years of global analyses of atmospheric fields in support of the needs of the research and climate monitoring communities. This effort involved the recovery of land surface, ship, rawinsonde, pibal, aircraft, satellite, and other data. These

In 1997, during the late stages of production of NCEP–NCAR Global Reanalysis (GR), exploration of a regional reanalysis project was suggested by the GR project's Advisory Committee, “particularly if the RDAS [Regional Data Assimilation System] is significantly better than the global reanalysis at capturing the regional hydrological cycle, the diurnal cycle and other important features of weather and climate variability.” Following a 6-yr development and production effort, NCEP's North American Regional Reanalysis (NARR) project was completed in 2004, and data are now available to the scientific community. Along with the use of the NCEP Eta model and its Data Assimilation System (at 32-km–45-layer resolution with 3-hourly output), the hallmarks of the NARR are the incorporation of hourly assimilation of precipitation, which leverages a comprehensive precipitation analysis effort, the use of a recent version of the Noah land surface model, and the use of numerous other datasets that are additional or improved compared to the GR. Following the practice applied to NCEP's GR, the 25-yr NARR retrospective production period (1979–2003) is augmented by the construction and daily execution of a system for near-real-time continuation of the NARR, known as the Regional Climate Data Assimilation System (R-CDAS). Highlights of the NARR results are presented: precipitation over the continental United States (CONUS), which is seen to be very near the ingested analyzed precipitation; fits of tropospheric temperatures and winds to rawinsonde observations; and fits of 2-m temperatures and 10-m winds to surface station observations. The aforementioned fits are compared to those of the NCEP–Department of Energy (DOE) Global Reanalysis (GR2). Not only have the expectations cited above been fully met, but very substantial improvements in the accuracy of temperatures and winds compared to that of GR2 are achieved throughout the troposphere. Finally, the numerous datasets produced are outlined and information is provided on the data archiving and present data availability.

Abstract The second version of the NCEP Climate Forecast System (CFSv2) was made operational at NCEP in March 2011. This version has upgrades to nearly all aspects of the data assimilation and forecast model components of the system. A coupled reanalysis was made over a 32-yr period (1979–2010), which provided the initial conditions to carry out a comprehensive reforecast over 29 years (1982–2010). This was done to obtain consistent and stable calibrations, as well as skill estimates for the operational subseasonal and seasonal predictions at NCEP with CFSv2. The operational implementation of the full system ensures a continuity of the climate record and provides a valuable up-to-date dataset to study many aspects of predictability on the seasonal and subseasonal scales. Evaluation of the reforecasts show that the CFSv2 increases the length of skillful MJO forecasts from 6 to 17 days (dramatically improving subseasonal forecasts), nearly doubles the skill of seasonal forecasts of 2-m temperatures over the United States, and significantly improves global SST forecasts over its predecessor. The CFSv2 not only provides greatly improved guidance at these time scales but also creates many more products for subseasonal and seasonal forecasting with an extensive set of retrospective forecasts for users to calibrate their forecast products. These retrospective and real-time operational forecasts will be used by a wide community of users in their decision making processes in areas such as water management for rivers and agriculture, transportation, energy use by utilities, wind and other sustainable energy, and seasonal prediction of the hurricane season.

We present the impact tests that preceded the most recent operational upgrades to the land surface model used in the National Centers for Environmental Prediction (NCEP) mesoscale Eta model, whose operational domain includes North America. These improvements consist of changes to the “Noah” land surface model (LSM) physics, most notable in the area of cold season processes. Results indicate improved performance in forecasting low‐level temperature and humidity, with improvements to (or without affecting) the overall performance of the Eta model quantitative precipitation scores and upper air verification statistics. Remaining issues that directly affect the Noah LSM performance in the Eta model include physical parameterizations of radiation and clouds, which affect the amount of available energy at the surface, and stable boundary layer and surface layer processes, which affect surface turbulent heat fluxes and ultimately the surface energy budget.

A new technique is presented in which half-hourly global precipitation estimates derived from passive microwave satellite scans are propagated by motion vectors derived from geostationary satellite infrared data. The Climate Prediction Center morphing method (CMORPH) uses motion vectors derived from half-hourly interval geostationary satellite IR imagery to propagate the relatively high quality precipitation estimates derived from passive microwave data. In addition, the shape and intensity of the precipitation features are modified (morphed) during the time between microwave sensor scans by performing a time-weighted linear interpolation. This process yields spatially and temporally complete microwave-derived precipitation analyses, independent of the infrared temperature field. CMORPH showed substantial improvements over both simple averaging of the microwave estimates and over techniques that blend microwave and infrared information but that derive estimates of precipitation from infrared data when passive microwave information is unavailable. In particular, CMORPH outperforms these blended techniques in terms of daily spatial correlation with a validating rain gauge analysis over Australia by an average of 0.14, 0.27, 0.26, 0.22, and 0.20 for April, May, June–August, September, and October 2003, respectively. CMORPH also yields higher equitable threat scores over Australia for the same periods by an average of 0.11, 0.14, 0.13, 0.14, and 0.13. Over the United States for June–August, September, and October 2003, spatial correlation was higher for CMORPH relative to the average of the same techniques by an average of 0.10, 0.13, and 0.13, respectively, and equitable threat scores were higher by an average of 0.06, 0.09, and 0.10, respectively.

This first paper of the two-part series describes the objectives of the community efforts in improving the Noah land surface model (LSM), documents, through mathematical formulations, the augmented conceptual realism in biophysical and hydrological processes, and introduces a framework for multiple options to parameterize selected processes (Noah-MP). The Noah-MP's performance is evaluated at various local sites using high temporal frequency data sets, and results show the advantages of using multiple optional schemes to interpret the differences in modeling simulations. The second paper focuses on ensemble evaluations with long-term regional (basin) and global scale data sets. The enhanced conceptual realism includes (1) the vegetation canopy energy balance, (2) the layered snowpack, (3) frozen soil and infiltration, (4) soil moisture-groundwater interaction and related runoff production, and (5) vegetation phenology. Sample local-scale validations are conducted over the First International Satellite Land Surface Climatology Project (ISLSCP) Field Experiment (FIFE) site, the W3 catchment of Sleepers River, Vermont, and a French snow observation site. Noah-MP shows apparent improvements in reproducing surface fluxes, skin temperature over dry periods, snow water equivalent (SWE), snow depth, and runoff over Noah LSM version 3.0. Noah-MP improves the SWE simulations due to more accurate simulations of the diurnal variations of the snow skin temperature, which is critical for computing available energy for melting. Noah-MP also improves the simulation of runoff peaks and timing by introducing a more permeable frozen soil and more accurate simulation of snowmelt. We also demonstrate that Noah-MP is an effective research tool by which modeling results for a given process can be interpreted through multiple optional parameterization schemes in the same model framework. Copyright © 2011 by the American Geophysical Union.

Previous estimates of land-atmosphere interaction (the impact of soil moisture on precipitation) have been limited by a lack of observational data and by the model dependence of computational estimates. To counter the second limitation, a dozen climate-modeling groups have recently performed the same highly controlled numerical experiment as part of a coordinated comparison project. This allows a multimodel estimation of the regions on Earth where precipitation is affected by soil moisture anomalies during Northern Hemisphere summer. Potential benefits of this estimation may include improved seasonal rainfall forecasts.

In this paper, the incorporation of a simple atmospheric boundary layer diffusion scheme into the NCEP Medium-Range Forecast Model is described. A boundary layer diffusion package based on the Troen and Mahrt nonlocal diffusion concept has been tested for possible operational implementation. The results from this approach are compared with those from the local diffusion approach, which is the current operational scheme, and verified against FIFE observations during 9–10 August 1987. The comparisons between local and nonlocal approaches are extended to the forecast for a heavy rain case of 15–17 May 1995. The sensitivity of both the boundary layer development and the precipitation forecast to the tuning parameters in the nonlocal diffusion scheme is also investigated. Special attention is given to the interaction of boundary layer processes with precipitation physics. Some results of parallel runs during August 1995 are also presented.

The analysis of univariate or multivariate time series provides crucial information to describe, understand, and predict climatic variability. The discovery and implementation of a number of novel methods for extracting useful information from time series has recently revitalized this classical field of study. Considerable progress has also been made in interpreting the information so obtained in terms of dynamical systems theory. In this review we describe the connections between time series analysis and nonlinear dynamics, discuss signal‐to‐noise enhancement, and present some of the novel methods for spectral analysis. The various steps, as well as the advantages and disadvantages of these methods, are illustrated by their application to an important climatic time series, the Southern Oscillation Index. This index captures major features of interannual climate variability and is used extensively in its prediction. Regional and global sea surface temperature data sets are used to illustrate multivariate spectral methods. Open questions and further prospects conclude the review.

The primary focus of this review is tropical‐extratropical interactions and especially the issues involved in determining the response of the extratropical atmosphere to tropical forcing associated with sea surface temperature (SST) anomalies. The review encompasses observations, empirical studies, theory and modeling of the extratropical teleconnections with a focus on developments over the Tropical Oceans‐Global Atmosphere (TOGA) decade and the current state of understanding. In the tropical atmosphere, anomalous SSTs force anomalies in convection and large‐scale overturning with subsidence in the descending branch of the local Hadley circulation. The resulting strong upper tropospheric divergence in the tropics and convergence in the subtropics act as a Rossby wave source. The climatological stationary planetary waves and associated jet streams, especially in the northern hemisphere, can make the total Rossby wave sources somewhat insensitive to the position of the tropical heating that induces them and thus can create preferred teleconnection response patterns, such as the Pacific‐North American (PNA) pattern. However, a number of factors influence the dispersion and propagation of Rossby waves through the atmosphere, including zonal asymmetries in the climatological state, transients, and baroclinic and nonlinear effects. Internal midlatitude sources can amplify perturbations. Observations, modeling, and theory have clearly shown how storm tracks change in response to changes in quasi‐stationary waves and how these changes generally feedback to maintain or strengthen the dominant perturbations through vorticity and momentum transports. The response of the extratropical atmosphere naturally induces changes in the underlying surface, so that there are changes in extratropical SSTs and changes in land surface hydrology and moisture availability that can feedback and influence the total response. Land surface processes are believed to be especially important in spring and summer. Anomalous SSTs and tropical forcing have tended to be strongest in the northern winter, and teleconnections in the southern hemisphere are weaker and more variable and thus more inclined to be masked by natural variability. Occasional strong forcing in seasons other than winter can produce strong and identifiable signals in the northern hemisphere and, because the noise of natural variability is less, the signal‐to‐noise ratio can be large. The relative importance of tropical versus extratropical SST forcings has been established through numerical experiments with atmospheric general circulation models (AGCMs). Predictability of anomalous circulation and associated surface temperature and precipitation in the extratropics is somewhat limited by the difficulty of finding a modest signal embedded in the high level of noise from natural variability in the extratropics, and the complexity and variety of the possible feedbacks. Accordingly, ensembles of AGCM runs and time averaging are needed to identify signals and make predictions. Strong anomalous tropical forcing provides opportunities for skillful forecasts, and the accuracy and usefulness of forecasts is expected to improve as the ability to forecast the anomalous SSTs improves, as models improve, and as the information available from the mean and the spread of ensemble forecasts is better utilized.

The Global Precipitation Climatology Project (GPCP) has released the GPCP Version 1 Combined Precipitation Data Set, a global, monthly precipitation dataset covering the period July 1987 through December 1995. The primary product in the dataset is a merged analysis incorporating precipitation estimates from low-orbit-satellite microwave data, geosynchronous-orbit-satellite infrared data, and rain gauge observations. The dataset also contains the individual input fields, a combination of the microwave and infrared satellite estimates, and error estimates for each field. The data are provided on 2.5 x2.5 latitude-longitude global grids. Preliminary analyses show general agreement with prior studies of global precipitation and extends prior studies of El Nino-Southern Oscillation precipitation patterns. At the regional scale there are systematic differences with standard climatologies.

An algorithm is developed to construct global gridded fields of monthly precipitation by merging estimates from five sources of information with different characteristics, including gauge-based monthly analyses from the Global Precipitation Climatology Centre, three types of satellite estimates [the infrared-based GOES Precipitation Index, the microwave (MW) scattering-based Grody, and the MW emission-based Chang estimates], and predictions produced by the operational forecast model of the European Centre for Medium-Range Weather Forecasts. A two-step strategy is used to: 1) reduce the random error found in the individual sources and 2) reduce the bias of the combined analysis. First, the three satellite-based estimates and the model predictions are combined linearly based on a maximum likelihood estimate, in which the weighting coefficients are inversely proportional to the squares of the individual random errors determined by comparison with gauge observations and subjective assumptions. This combined analysis is then blended with an analysis based on gauge observations using a method that presumes that the bias of the gauge-based field is small where sufficient gauges are available and that the gradient of the precipitation field is best represented by the combination of satellite estimates and model predictions elsewhere. The algorithm is applied to produce monthly precipitation analyses for an 18-month period from July 1987 to December 1988. Results showed substantial improvements of the merged analysis relative to the individual sources in describing the global precipitation field. The large-scale spatial patterns, both in the Tropics and the extratropics, are well represented with reasonable amplitudes. Both the random error and the bias have been reduced compared to the individual data sources, and the merged analysis appears to be of reasonable quality everywhere. However, the actual quality of the merged analysis depends strongly on our uncertain and incomplete knowledge of the error structures of the individual data sources.

Results are presented from the second phase of the multiinstitution North American Land Data Assimilation System (NLDAS‐2) research partnership. In NLDAS, the Noah, Variable Infiltration Capacity, Sacramento Soil Moisture Accounting, and Mosaic land surface models (LSMs) are executed over the conterminous U.S. (CONUS) in realtime and retrospective modes. These runs support the drought analysis, monitoring and forecasting activities of the National Integrated Drought Information System, as well as efforts to monitor large‐scale floods. NLDAS‐2 builds upon the framework of the first phase of NLDAS (NLDAS‐1) by increasing the accuracy and consistency of the surface forcing data, upgrading the land surface model code and parameters, and extending the study from a 3‐year (1997–1999) to a 30‐year (1979–2008) time window. As the first of two parts, this paper details the configuration of NLDAS‐2, describes the upgrades to the forcing, parameters, and code of the four LSMs, and explores overall model‐to‐model comparisons of land surface water and energy flux and state variables over the CONUS. Focusing on model output rather than on observations, this study seeks to highlight the similarities and differences between models, and to assess changes in output from that seen in NLDAS‐1. The second part of the two‐part article focuses on the validation of model‐simulated streamflow and evaporation against observations. The results depict a higher level of agreement among the four models over much of the CONUS than was found in the first phase of NLDAS. This is due, in part, to recent improvements in the parameters, code, and forcing of the NLDAS‐2 LSMs that were initiated following NLDAS‐1. However, large inter‐model differences still exist in the northeast, Lake Superior, and western mountainous regions of the CONUS, which are associated with cold season processes. In addition, variations in the representation of sub‐surface hydrology in the four LSMs lead to large differences in modeled evaporation and subsurface runoff. These issues are important targets for future research by the land surface modeling community. Finally, improvement from NLDAS‐1 to NLDAS‐2 is summarized by comparing the streamflow measured from U.S. Geological Survey stream gauges with that simulated by four NLDAS models over 961 small basins.