Go directly to Data Access Information
Relevant references concerning the DAO assimilation system (GEOS-1 DAS)
DATA ASSIMILATION OFFICE code 910.3 NASA/Goddard Space Flight Center Greenbelt MD 20771
DATA ASSIMILATION OFFICE code 910.3 NASA/Goddard Space Flight Center Greenbelt MD 20771
The DAO has produced this benchmark TOGA COARE dataset using version 1 of the GEOS-1 assimilation system. The objectives in producing these data are threefold. First, it is believed that the absence of spin-up in the hydrological cycle and the use of a fixed assimilation system will make these data extremely useful for a wide range of process studies. Secondly, these data will serve as the benchmark against which future assimilations, which test the impact of new data types, will be compared. Thirdly, making the data available to the larger scientific community will attract valuable feedback on the quality and limitations of the assimilated data. This feedback will help guide future development of the system.
These parameters consist of various surface boundary conditions, prognostic and diagnostic surface quantities and vertically integrated/averaged fields. The fields are made available every six hours, though it should be noted the boundary condition fields are interpolated linearly in time from monthly mean or climatological (surface wetness) fields.
These parameters are "instantaneous snapshots" reported at 6 hour intervals. Each product is reported as a full three-dimensional field at 18 pressure levels (1000, 950, 900, 850, 800, 700, 600, 500, 400, 300, 250, 200, 150, 100, 70, 50, 30, 20mb).
These parameters are 3-hour averages, reported at the endpoint of the 3-hour interval.
These parameters are 3-hour averages, reported at the endpoint of the 3-hour interval.
These parameters are 3-hour averages, reported at the endpoint of the 3-hour interval.
These parameters are 3-hour averages, reported at the endpoint of the 3-hour interval.
These parameters are 6-hour averages, reported at the midpoint of the 6-hour interval. Each product is reported as a full three-dimensional field at 18 pressure levels (1000, 950, 900, 850, 800, 700, 600, 500, 400, 300, 250, 200, 150, 100, 70, 50, 30, 20mb).
Each parameter is described further in the Parameter Table, and in more detail in the document: geos1.0_gcm.doc.ps. We request that if a user discovers additional problems that the user will alert the DAO by e-mail at:
toga@dao.gsfc.nasa.govRESOLUTION
The tropospheric version of the OI analysis scheme being used for the control assimilation has been carried out at a horizontal resolution of 2 degree latitude by 2.5 degree longitude at 14 upper-air pressure levels (20, 30, 50, 70, 100, 150, 200, 250, 300, 400, 500, 700, 850, 1000 mb) and at sea level. The analysis increments are computed every 6 hours using observations from a +/- 3 hour data window centered on the analysis times (00, 06, 12, and 18 UTC). The innovation vector (observation minus background forecast) used as input to the OI is computed using a single forecast valid at the analysis time.
THE ASSIMILATION
The OI scheme is multivariate in geopotential height and winds and employs a damped cosine function for the horizontal correlation of model prediction error. The height-wind cross-correlation model is geostrophic and scaled to zero at the equator. The multivariate surface analysis scheme over the oceans adopts an Ekman balance for the pressure-wind analysis. The moisture analysis for mixing ratio employs only rawinsonde moisture data. Although, it should be noted, that all TOGA COARE upper air moisture data was withheld from the assimilation. All grid point analyses are done using up to 75 nearby observations from within a radius of 1600 km.The assimilation system does not include an initialization scheme and relies on the damping properties of a Matsuno time differencing scheme to control initial imbalances generated by the insertion of observations. However, the initial imbalances and spinup have been greatly reduced over earlier versions by the introduction of an incremental analysis update (IAU) procedure [Bloom et al., 1991]. As shown in the figure of the IAU procedure,
![]()
standard OI analysis increments are computed at the analysis times (0, 6, 12, 18Z). The increments are then inserted gradually into the AGCM by rerunning the forecast and adding a fraction of the increment at each model time step. Over the 6 hour period centered at the analysis time the full effect of the increment is realized. The assimilation thus effectively consists of a continuous AGCM forecast with additional heat, momentum, moisture and mass source terms updated every 6 hours from observations. It is the output of this IAU procedure that makes up the dataset provided here. An important difference between the IAU scheme and the usual Newtonian nudging procedure is that the IAU forcing terms are held constant over the insertion period, while in Newtonian nudging they are proportional to the difference between a target analysis and the instantaneous current model state.
BUDGET CALCULATIONS
The implication of the IAU procedure for performing budget (e.g., moisture or heat) calculations with the assimilated data is that, in order to balance the budget, one must include the analysis increments and the filling of negative specific humidities as additional forcing terms. If the assimilating model had no bias the mean analysis increments would, of course, be zero and the increments would have no contribution to the mean budget equations. The current model (GEOS-1), does have a bias, and while the initial timeseries dataset provided by the DAAC does not include the analysis increments, they (the increments) can be obtained as the residual of the other terms in the budget equations. The filling of negative specific humidities (QFILL) are provided in the timeseries.
Note that all diagnostic parameters are calculated from AGCM's physical parameterizations in a manner consistent with the prognostic fields (See Optimal Interpolation).
The input observational database is one that has been accumulated at the GLA. These data have been obtained directly from NMC, and do not include data which came in after the cut-off time for the operational NMC system. In addition to this source, some TOVS temperature soundings have come directly from the National Oceanic and Atmospheric Administration (NOAA) National Environmental Satellite, Data and Information Service (NESDIS), and gaps have been filled with data from the National Climatic Data Center (NCDC) in Ashville, NC. For the global sea level pressure and near surface wind analysis over the oceans, data from surface land synoptic reports (sea level pressure only), ships and buoys are used. The upper-air analyses of height, wind and moisture incorporate the data from rawinsondes, dropwindsondes, rocketsondes, aircraft winds, cloud tracked winds, and thicknesses from the historical TOVS soundings produced by NOAA NESDIS. The satellite heights are computed using a reference level which depends on the analyzed sea level pressure. The only situation where data are used as a proxy for un-collected observations are 1000 mb height observations which are generated above pressure reports from ships. These serve to further couple the surface and upper-air analyses.For a description of how these data are used in the model, see the description of optimal interpolation.
NAME MBytes PARAMETERS SFCPROG 58 PHIS,ALBD,GWET,PSMPTOP,GTMP,SLP,LWI,UBAR,VBAR UWND 21 UWND VWND 21 VWND HGHT 21 HGHT TMPU 21 TMPU SPHU 21 SPHU QQ 21 QQ HGHTE 17 HGHTE DIAG1 21 PSMPTOP,PREACC,PRECON,EVAP,VINTUQ,VINTVQ,VINTUT,VINTVT,QINT DIAG2 24 PSMPTOP,UFLUX,VFLUX,HFLUX,CT,CU,WINDS,USTAR,Z0,PBL DIAG3 26 PSMPTOP,RADLWG,RADSWG,OLR,OLRCLR,LWGCLR,RADSWT,OSR,OSRCLR,SWGCLR,CLDFRC DIAG4 31 PSMPTOP,TG,TS,QS,DPDT,U2M,V2M,T2M,Q2M,U10M,V10M,T10M,Q10M TURBU 21 TURBU TURBV 21 TURBV TURBT 21 TURBT TURBQ 21 TURBQ MOISTT 21 MOISTT MOISTQ 21 MOISTQ RADLW 21 RADLW RADSW 21 RADSW OMEGA 21 OMEGA QFILL 21 QFILL
For more detailed information see geos1.0_gcm.doc.ps.PROGNOSTIC VARIABLES
Surface Prognostic
NAME UNIT DESCRIPTION PHIS (m/s)^2 surface geopotential height (z * gravity) ALBD 0-1 proportion of solar radiation reflected by the surface (0-1) GWET 0-1 ground wetness from off-line bucket model [Schemm et al, 1992] PSMPTOP mb surface pressure minus top of the atmosphere pressure (PTOP= 10mb) GTMP K surface ground temperature SLP mb sea level pressure LWI flag water=1, land=2, permanent ice=3, sea ice=4 flags UBAR m/s vertically integrated (barotropic) zonal wind VBAR m/s vertically integrated (barotropic) meridional windUpper Air Prognostic
NAME UNIT DESCRIPTION UWND m/s zonal wind speed VWND m/s meridional wind speed HGHT m geopotential height TMPU K temperature profiles SPHU g/kg specific humidity profiles QQ (m/s)^2 turbulent kinetic energy HGHTE m std dev of the height errorDIAGNOSTIC VARIABLES
Surface Diag 1 Precip/Evap
NAME UNIT DESCRIPTION PSMPTOP mb surface pressure minus top of the atmosphere pressure (PTOP= 10mb) PREACC mm/dy total precipitation PRECON mm/dy convective precipitation EVAP mm/dy surface evaporation VINTUQ m/s g/kg vertically averaged zonal wind * specific humidity (U*Q) VINTVQ m/s g/kg vertically averaged meridional wind * specific humidity (V*Q) VINTUT m/s K vertically averaged zonal wind * temperature (U*T) VINTVT m/s K vertically averaged meridional wind * temperature (V*T) QINT g/cm^2 vertically integrated ("total")precipitable waterSurface Diag 2 Momentum/Heat Flux
NAME UNIT DESCRIPTION PSMPTOP mb surface pressure minus top of the atmosphere pressure (PTOP= 10mb) UFLUX N/m^2 zonal momentum surface stress VFLUX N/m^2 meridional momentum surface stress HFLUX W/m^2 surface flux of sensible heat CT - surface drag coefficient for temperature and specific humidity CU - surface drag coefficient for winds WINDS m/s surface wind speed USTAR m/s friction velocity Z0 m surface roughness PBL mb planetary boundary layer depthSurface Diag 3 Radiation
NAME UNIT DESCRIPTION PSMPTOP mb surface pressure minus top of the atmosphere pressure (PTOP= 10mb) RADLWG W/m^2 net upward LW radiation at the surface RADSWG W/m^2 net downward SW radiation at the surface OLR W/m^2 outgoing longwave radiation OLRCLR W/m^2 outgoing longwave radiation clear sky LWGCLR W/m^2 surface longwave radiation clear sky RADSWT W/m^2 incident SW radiation at top of the atmosphere OSR W/m^2 outgoing shortwave radiation OSRCLR W/m^2 outgoing shortwave radiation clear sky SWGCLR W/m^2 surface SW radiation clear sky. CLDFRC 0-1 2-dimensional total cloud fraction (0-1)Surface Diag 4 Near Surface
NAME UNIT DESCRIPTION PSMPTOP mb surface pressure minus top of the atmosphere pressure (PTOP= 10mb) TG K ground temperature TS K surface air temperature QS g/kg saturation specific humidity at the surface DPDT mb/dy surface pressure tendency U2M m/s zonal wind at 2m V2M m/s meridional wind at 2m T2M K temperature at 2m Q2M kg/kg specific humidity at 2m U10M m/s zonal wind at 10m V10M m/s meridional wind at 10m T10M K temperature at 10m Q10M kg/kg specific humidity at 10mUpper Air Diagnostic
NAME UNIT DESCRIPTION TURBU m/s/dy zonal momentum changes due to turbulence TURBV m/s/dy meridional momentum changes due to turbulence TURBT K/dy temperature changes due to turbulence TURBQ g/kg/dy moisture changes due to turbulence MOISTT K/dy temperature changes due to moist processes MOISTQ g/kg/dy moisture changes due to moist processes RADLW K/dy temperature changes due to LW radiation RADSW K/dy temperature changes due to SW radiation OMEGA mb/dy vertical velocity QFILL g/kg/dy filling of negative specific humidities

It was created by the fortran read program called: single_level_read.f.
This version of the AGCM is run without a land surface model. For the assimilation described here, soil moisture is computed off-line based on a simple bucket model using monthly mean observed surface air temperature and precipitation [Schemm et al., 1992]. The climatology calculated using these monthly means was used in this assimilation run. The snow line and surface albedo are prescribed and vary with the season. The sea surface temperature is updated according to the observed monthly mean values provided by the Climate Analysis Center at NMC and the Center for Ocean, Land and Atmosphere (COLA) at the University of Maryland. Long-term plans call for the incorporation of both land-surface and ocean models.
These files should be retrieved via anonymous FTP.
toga_doc.txt This file contains information about this dataset.
volume_4.ps.Z This compressed postscript file contains information about the assimilation system (GEOS-1 DAS).
geos1.0_gcm.doc.ps This file contains information about the gcm and how the diagnostics are calculated.
BAMS Article (available online) decribing the DAO 5-year assimilation.
Arakawa, A., and M.J. Suarez, 1983: Vertical differencing of the primitive equations in sigma coordinates. Mon. Wea. Rev., 111, 34-45.
Bates, J. R., S. Moorthi and R. W. Higgins, 1993: A global multilevel atmospheric model using a vector semi-Lagrangian finite-difference scheme. Part I: Adiabatic formulation. Mon. Wea. Rev., 121, 244-263.
Bengtsson, L. and J. Shukla, 1988: Integration of space and in situ observations to study global climate change. Bull. Amer. Meteor. Soc., 69, 1130-1143.
Bloom, S. C., L. L. Takacs and E. Brin, 1991: A scheme to incorporate analysis increments gradually in the GLA assimilation system. Ninth Conference on Numerical Weather Prediction, Denver, CO.
Burridge, D. M. and J. Haseler, 1977: A model for medium range weather forecasting-adiabatic formulation. Tech. Rep. No. 4, European Centre for Medium Range Weather Forecasts, Bracknell, Berks, U. K., 46pp.
Chou, M.D., 1984: Broadband water vapor transmission functions for atmospheric IR flux computations. J. Atmos. Sci. , 41, 1775-1778.
Chou, M.D.,and L. Peng, 1983: A parameterization of the absorption in the 15m CO2 spectral region with application to climate sensitivity studies. J. Atmos. Sci. , 40, 2183-2192.
Davies, R., 1982: Documentation of the solar radiation parameterization in the GLAS climate model. NASA Tech. Memo. 83961, 57pp., Goddard Space Flight Center, Greenbelt, MD 20771.
Harshvardhan, R. Davies, D.A. Randall, and T.G. Corsetti, 1987: A fast radiation parameterization for atmospheric circulation models. J. Geophys. Res., 92, 1009-1016.
Helfand, H.M. and J.C. Labraga, 1988: Design of a non-singular level 2.5 second-order closure model For the prediction of atmospheric turbulence. J. Atmos. Sci., 45, 113-132.
Helfand, H. M. and S. D. Schubert, 1993: Contribution of the Great Plains low-level jet to the simulated continental moisture budget of the United States. J. Climate, (submitted).
Joseph, J.H., W.J. Wiscombe and J.E. Weinman, 1976: The delta-Eddington approximation for radiative flux transfer. J. Atmos Sci., 33,2452-2459.
Kalnay, E. and R. Jenne, 1991: Summary of the NMC/NCAR reanalysis workshop of April 1991. Bull. Amer. Meteor. Soc., 72, 1897-1904.
King, M.D. and Harshvardhan, 1986: Comparitive accuracy of selected multiple scattering approximations. J. Atmos. Sci., 43, 784-801.
Lacis, A.A. and J.E. Hansen, 1974: A parameterization for the absorption of solar radiation in the Earth's atmosphere. J. Atmos. Sci. , 31, 118-133.
Lanczos, C., 1966: Discourse on Fourier Series. Hafner Publishing.
Moorthi, S. and M. J. Suarez, 1992: Relaxed Arakawa-Schubert: A parameterization of moist convection for general circulation models. Mon. Wea. Rev., 120,978-1002.
National Research Council, 1991: Four-dimensional model assimilation of data: a strategy for earth system sciencs. Report from the Panel on Model-Assimilated Data Sets for Atmospheric and Oceanic Research, National Academy Press, Washington D.C., 1991.
Rodgers, C.D., 1968: Some extensions and applications of the new random model for molecular band transmission. Quart. J. R. Met. Soc., 94, 99-102.
Rood, R. B., J. E. Nielsen, R. S. Stolarski, A. R. Douglass, J. A. Kaye and D. J. Allen, 1992: Episodic total ozone minima and associated effects on heterogeneous chemistry and lower stratospheric transport. J. Geophys. Res., 97, 7979-7996.
Rosenfield, J.E., M.R. Schoeberl, and M.A. Geller, 1987: A computation of the stratospheric diabatic circulation using an accurate radiative transfer model. J. Atmos. Sci. , 44, 859-876.
Sadourney, R., 1975: The dynamics of finite difference models of the shallow water equations. J. Atmos. Sci. , 32, 680-689.
Schemm, J.-K., S. Schubert, J. Terry and S. Bloom, 1992: Estimates of monthly mean soil moisture for 1979-89, NASA Tech. Memo. 104571, pp 252, Oct. 1992.
Sud, Y. and A. Molod, 1988: The roles of dry convection, cloud-radiation feedback processes and the influence of recent improvements in the parameterization of convection in the GLA AGCM. Mon. Wea. Rev., 116., 2366-2387.
Susskind, J., J. Rosenfield, D. Reuter and M. T. Chahine, 1984: Remote sensing of weather and climate parameters from HIRS2/MSU on TIROS-N, J. Geophys. Res., 89, 4677- 4697.
PARAMETER (IM=89,JNP=27)
REAL FIELD(IM,JNP)
DO 1 ITIMES=1,NTIMES
DO 2 IXX=1,NXX
READ(8) FIELD
2 CONTINUE
1 CONTINUE
For each DATA FILE there is an associated TABLE FILE which provides which values to use for these parameters. The TABLE FILE has the same name as the DATA FILE except the part of the name ".prs." is replaced by ".tabl.".
In the TABLE FILE:
A program written by the DAAC to read the upper air data products, based on this template is available: UPPER AIR SAMPLE READ PROGRAM.
Assign statements (for Cray) 1) cray data (iau and restart files) assign -a $data1 fort.21 2) ieee data (all other files) assign -a $data1 -N ieee -F f77 fort.21
Goddard DAAC code 610.2 NASA/Goddard Space Flight Center GREENBELT MD 20771
FAX:301-614-5268
AGCM Atmospheric General Circulation Model COARE Coupled Ocean-Atmosphere Response Experiment COLA Center for Ocean, Land and Atmosphere at the University of Maryland DAAC Distributed Active Archive Center DAO Data Assimilation Office ECMWF European Center for Medium Range Forecasts EOS Earth Observing System GEOS-1 Goddard Earth Observing System version 1 model GEWEX Global Energy and Water Cycle Experiment GLA Goddard Lab for Atmospheres GSFC Goddard Space Flight Center HIRS2 High-Resolution Infrared Sounder 2 IAU Incremental Analysis Update IMS Information Management System KF Kalman Filtering LW Longwave (same as IR) NASA National Aeronautics and Space Administration NCAR National Center for Atmospheric Research NCDC National Climatic Data Center NESDIS National Environmental Satellite Data and Information Service NH Northern Hemisphere NMC National Meteorological Center NOAA National Oceanic and Atmospheric Administration NRC National Research Council OI Optimal Interpolation PBL Planetary Boundary Layer RAS Relaxed Arakawa-Schubert RHS Right Hand Side (e.g. of an equation) SPADE Stratospheric Photochemistry, Aerosols, and Dynamics Expedition STRATAN Stratospheric version of GEOS-1 SH Southern Hemisphere SW Shortwave (same as VISIBLE) TIROS Television Infrared Observing Satellite TOA Top Of the Atmosphere TOGA Tropical Ocean Global Atmosphere TOVS TIROS Operational Vertical Sounder UARS Upper Atmosphere Research Satellite