SRB_REL3.1_LONGWAVE_3HRLY - GEWEX Longwave 3-Hourly Data Set README File 1.0 Introduction This README file provides information on the SRB_REL3.1_LONGWAVE_3HRLY data set. The data set contains 3-hourly global fields of six longwave (LW) surface and Top of Atmosphere (TOA) radiative parameters, in addition to a day/night flag, derived with the Longwave algorithm of the NASA World Climate Research Programme/Global Energy and Water-Cycle Experiment (WCRP/GEWEX) Surface Radiation Budget (SRB) Project. If users have questions, please contact the Langley Atmospheric Science Data Center (ASDC) User and Data Services Office at: Atmospheric Science Data Center User and Data Services Office Mail Stop 157D NASA Langley Research Center Hampton, Virginia 23681-2199 U.S.A. E-mail: support-asdc@earthdata.nasa.gov Phone: (757)864-8656 FAX: (757)864-8807 URL: http://eosweb.larc.nasa.gov This readme includes the following sections: 1.0 Introduction 2.0 Data Set Description 2.1 Data Quality 2.2 Input Information 2.3 Grid Description 2.4 Points of Contact 3.0 Format and Packaging 4.0 Science Parameters Information 5.0 Sample Read Software Description 6.0 Implementing the Sample Read Software 7.0 Sample Output 8.0 Additional Derivable Parameters References 2.0 Data Set Description There are a total of seven parameters in these files as follows: 1. Day/Night flag (daynite; 1=Day, 0=Night) 2. TOA Upward Clear-Sky Flux/Clear-sky Outgoing Longwave Radiation (OLR) (clr_toa_up) 3. Surface Clear-sky Upward Longwave Flux (clr_sfc_up) 4. Surface Clear-sky Downward Longwave Flux (clr_sfc_down) 5. TOA All-Sky Upward Longwave Flux/OLR (toa_up) 6. Surface All-Sky Upward Longwave Flux (sfc_up) 7. Surface All-Sky Downward Longwave Flux (sfc_down) These parameters were derived originally on a 3-hourly temporal resolution (i.e., a global instantaneous gridded field every 3 hours), at UT hours 00, 03, 06, 09, 12, 15, 18, and 21 for every day of the month. The current version of the data set is Release 3.1. The GEWEX LW algorithm uses the Fu et al. (1997) thermal infrared radiative transfer code requiring atmospheric profile information, cloud, and surface properties. The sources for these inputs are briefly described below. A detailed description of the algorithm is currently being prepared for publication. Please contact the Dr. Paul W. Stackhouse Jr. at the address below for further details. Version History: Release 2.1: 12 year data set (July 1983-October 1995), on nested grid (described in Section 2.3), using GEOS-1 meteorological data. Release 2.5: 22 year data set (July 1983-June 2005). Using GEOS-4 meteorological inputs for the data set in place of GEOS-1. Release 3.0: 24.5 year data set (July 1983-December 2007). This version includes improved cloud properties in areas in missing and sun glint regions where ISCCP cloud retrievals aren't performed. Additionally, the IR radiative parameterization of ice clouds has been updated (Fu et al. 1998). The water vapor continuum has been updated (Kratz and Rose, 1999). An error in the ozone profile assignment is corrected. The surface vegetation type maps have been updated. This affects the surface emissivity values (Rutan et al. 2009). CO2 concentration value now varies month to month, based on monthly trend values from ftp://ftp.cmdl.noaa.gov/ccg/co2/trends/co2_mm_gl.txt Release 3.1: Corrections to nonphysical fluxes have been made to Rel. 3.0. Negative TOA fluxes in the 3-hourly data files were found to occur about 7 grid box times per month (out of 44016 grid boxes x 248 hours per month), with an additional 5 to 10 values per month identified as being unphysically low. The problem was found to be an numerical instability occurring due to an optimization switch in the Fortran compiler. The downwelling fluxes were also affected. So, the nonphysical values were replaced with a recomputation of those grid boxes using the same code but built without an optimization. Additionally, 3-hourly values of NaN's were detected and were traced to en error in the temperature profiles. The frequency of this occurrence was far more rare and was found to be mostly clustered for three months in 1987 and July 1990. After correcting the temperature profiles, the grid boxes were recomputed and replaced in the 3-hourly files using the same "non optimized" version of the code. The daily, monthly and 3-hourly monthly files were reprocessed using the improved 3-hourly files. Differences in the monthly averages proved to be small and mostly < 2 W m^-2 on a grid box level and < 0.01 W m^-2 on a global mean. Differences on the grid box level in the 3-hourly monthly and daily averages were mostly < 10 W m^-2. 2.1 Data Quality An assessment of the quality of these monthly average fluxes was accomplished by comparisons with corresponding ground-measured fluxes over a period from January 1992 to December 2007 from a number of sites of the Baseline Surface Radiation Network (BSRN). From the aggregate data set for all sites and years, mean bias was determined to be about 0.78 W/m**2 (0.2%, model fluxes higher), and the root mean square difference is 30.1 W/m**2 (9.8%). Uncertainties associated with operational BSRN measurements during this period are believed to be about +/- 5 W/m**2 (1.5%, Ellsworth Dutton, NOAA, BSRN Manager). Thus, mean bias for the present results is within the uncertainty of BSRN measurements. Errors for individual 3-hourly values may be different from the above values because those are subject to bias and random errors due to local meteorological conditions. 2.1.2. Indian Ocean Gap Artifact There is a visible and common artifact in much of the data set period, due to a lack of coverage from geostationary satellites over an area centered on 70 degrees east longitude. This situation , commonly called the Indian Ocean gap, occurs for all of the July 1983 - June 1998 time period, except for April 1988 - March 1989, when data from the INSAT satellite is available to cover the gap. In July of 1998, Meteosat-5 was moved over the gap area, eliminating the gap. When the Indian Ocean gap occurs, the gap area is covered by polar orbiting satellites, which can result in only one or two daytime overpasses per day. Geosynchronous temporal sampling during the daytime is 3-5 times per day depending upon the latitude (between 55 degrees North and South) and the time or year. In addition, the limbs of the geostationary satellites which bound the gap may suffer from spuriously high cloud amounts, due to large view angles. This results in an abrupt drop-off of cloud fraction in the gap as compared to the gap boundary. Downward longwave radiation is lower in the gap, creating an appearance of a flux discontinuity. For 3-hourly averaged fluxes, a discontinuity of magnitude less than 20 W/m**2 for TOA and 5 W/m**2 surface fluxes may appear in the Indian Ocean gap region. 2.2 Input Information Inputs to the algorithm were obtained from the following sources: Cloud parameters were derived from the International Satellite Cloud Climatology Project (Rossow and Schiffer, 1999) DX data product. The cloud pixels were separated into categories of high, middle and low where middle and low clouds could be composed of ice or water, while high clouds were composed of ice only. Cloud fractions and cloud optical depths were determined within these categories. Cloud particle sizes were assumed and cloud physical thicknesses were also assigned based upon information from literature. Random overlap is used between the high, middle and low layers to better approximate undercast conditions. Temperature and moisture profiles were obtained from the 4-D data assimilation Goddard EOS Data Assimilation System, level-4 (GEOS-4) obtained from the Global Modeling and Assimilation Office (GMAO) at NASA Goddard Space Flight Center (GSFC) (Bloom et al., 2005) Column ozone values for most of the duration of this dataset (July 1983 to December 2004) were obtained primarily from the Total Ozone Mapping Spectrometer (TOMS) archive. For the early period (July 1983-November 1994), TOMS data came from NIMBUS-7 and Meteor-3 satellites. There was an interruption of about 20 months (December 1994-July 1996) after which TOMS data from EP-TOMS became available in August 1996 and continued until December 2004. All gaps in TOMS data, including those over the polar night areas every year, were filled with column ozone values from TIROS Operational Vertical Sounder (TOVS) data. Column ozone data continued to be available beyond December 2004 from OMI instrument aboard Aura satellite but TOVS data, which is essential for filling the gaps in OMI data, developed some unexplained gaps of its own and became unusable. Beginning in January 2005, GEWEX/SRB started using a daily analysis ozone product from NOAA Climate Predictions Center (CPC), known as the Stratospheric Monitoring-group Ozone Blended Analysis (SMOBA). Surface emissivities were taken from a map developed at NASA LaRC (Wilber et al. 1999). 2.3 Grid Description The fluxes are generated on a nested grid, which contains 44016 cells. The grid has a resolution of 1 degree latitude globally, and longitudinal resolution ranging from 1 degree in the tropics and subtropics to 120 degrees at the poles. The first cell is Latitude 89-90 degrees South, Longitude 0-120 degrees East. The cells start at the Greenwich meridian and proceed east around the globe, then shift one degree to the north. The number of cells per latitude band starting at the South Pole are: 3, 45, 45, 45, 45, 45, 45, 45, 45, 45, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 180, 180, 180, 180, 180, 180, 180, 180, 180, 180, 180, 180, 180, 180, 180, 180, 180, 180, 180, 180, 180, 180, 180, 180, 180, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 360, 180, 180, 180, 180, 180, 180, 180, 180, 180, 180, 180, 180, 180, 180, 180, 180, 180, 180, 180, 180, 180, 180, 180, 180, 180, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 45, 45, 45, 45, 45, 45, 45, 45, 45, 3 The read software described below contains a subroutine to regrid the fluxes to 1 degree latitude by 1 degree longitude equal-angle grid using replication. 2.4 Points of Contact Scientific contact: Dr. Paul W. Stackhouse Jr. Mail Stop 420 21 Langley Boulevard NASA Langley Research Center Hampton, VA 23681-2199 U.S.A. E-mail: Paul.W.Stackhouse@nasa.gov Production Contact: Atmospheric Science Data Center User and Data Services Office Mail Stop 157D NASA Langley Research Center Hampton, VA 23681-2199 U.S.A. E-mail: support-asdc@earthdata.nasa.gov 3.0 Format and Packaging Each data file contains an entire month of 3-hourly global fields of the parameters described in Section 4.0 on an approximately 1 deg x 1 deg equal-area grid described in Section 2.3. The files contain binary data and are named according to the following convention: srb_rel3.1_longwave_3hrly_yyyymm.binary, where srb Project name, Surface Radiation Budget rel3.1 Release number for these data (Release 3.1) longwave Name of the algorithm, GEWEX Longwave 3hrly Time resolution of the data file yyyy 4-digit year mm 2-digit month binary file format 4.0 Science Parameters Information The files contain global fields of 3-hourly values of the day/night flag and the six radiative parameters on the nested grid. Each file has 7 records, containing one global field for every time period in each record. The parameters are: Name: Day/Night flag (Day=1, Night=0) Units: none Type: real Range: 0.0 or 1.0 Fill Values: n/a Scale Factor: None Name: Top-of-Atmosphere Clear-sky Upward LW Flux Units: Watts per square meter Type: Real Range: 50 to 600 Fill Values: -999.0 Scale Factor: None Name: Surface Clear-sky Upward LW Flux Units: Watts per square meter Type: Real Range: 50 to 800 Fill Values: -999.0 Scale Factor: None Name: Surface Clear-sky Downward LW Flux Units: Watts per square meter Type: Real Range: 50 to 600 Fill Values: -999.0 Scale Factor: None Name: Top-of-Atmosphere All-sky Upward LW Flux Units: Watts per square meter Type: Real Range: 50 to 600 Fill Values: -999.0 Scale Factor: None Name: Surface All-sky Upward LW Flux Units: Watts per square meter Type: Real Range: 50 to 800 Fill Values: -999.0 Scale Factor: None Name: Surface All-sky Downward LW Flux Units: Watts per square meter Type: Real Range: 50 to 600 Fill Values: -999.0 Scale Factor: None 5.0 Sample Read Software Description Sample read software written in Fortran-90, read_longwave_3hrly.f90, was developed for reading these data. The software constitutes the name of the input data file, accesses and reads it, using the information provided in the namelist file (longwave_3hrly.nml). The input files are read as direct-access binary on the nested (44016 box) grid. The software reads one or more of the 7 parameter fields, regrids them to an equal-angle 1 deg x 1 deg grid, and writes the output as ascii or binary format. The choice of file format (ascii or binary) and of the location of the output files is also made through the namelist file. A sample namelist file that would be used to read the July 1992 data file and write all parameters to an ascii format output file is presented below: &time_vars yr=1992 mon=7 ascii=.true. binary=.false. path_in='**** input file path here****' path_out='**** output file path here****' little_endian=.false. daynite=.true. clr_toa_up=.true. clr_sfc_up=.true. clr_sfc_down=.true. toa_up=.true. sfc_up=.true. sfc_down=.true. / There is a choice to convert the input fields from big endian to little endian byte order with the logical variable "little_endian" in the namelist. This applies to operating systems where byte order is stored opposite that of the machines used to create the data set, such as Linux and Macs with the Intel chip. If possible, a better choice for doing the conversion in these cases would be to use a compiler option. If using a compiler option, do not set little_endian to true. Both input and output fields have the same orientation: they start at the Greenwich meridian-south pole and go east and north from there. A limitation of this code is that it provides a complete global field of the specified parameters in the above orientation. The user should be easily able to extract values for any box or lat-lon region from these fields. 6.0 Implementing the Sample Read Software The sample read software can be compiled with any Fortran 90 or 95 compiler. To compile: % f90 -o run_longwave_3hrly read_longwave_3hrly.f90 The providers used a gfortran compiler but any F90/F95 compiler should work. Edit the namelist file to select month and year to be processed, choose the parameters to be read and the format of the output file. Run the software: % run_longwave_3hrly 6.1 Read Software Incompatibilities With some Fortran compilers the RECL keyword in the OPEN statement assumes record lengths are specified in 4 byte increments. If that is the case, then the following statement in the read program: open(10,file=file_in,status='old',access='direct', & form='unformatted', recl=44016*4) should be modified to: open(10,file=file_in,status='old',access='direct', & form='unformatted', recl=44016) The same should be done with the output binary file: open (15,file=outfile, form='unformatted', access='direct', & recl=nlon*nlat*4, status='replace') should be modified to: open (15,file=outfile, form='unformatted', access='direct', & recl=nlon*nlat, status='replace') 7.0 Sample Output The seven tables of numbers below show the values of the parameters contained in these files for latitude bands 45-51 (starting at the south pole) and longitude boxes 100-104 (starting at the Greenwich meridian) for hour 06 of day 14 of the month. Values for only a small lat-lon box for a single time are printed to the screen. When the software is run, the following information appears on the screen: ***************************************************************** * * * * * Data Set srb_rel3.1_longwave_3hrly Read Software * * * * Version: 1.0 * * * * Contact: Atmospheric Science Data Center * * User and Data Services Office * * Mail Stop 157D * * NASA Langley Research Center * * Hampton, Virginia 23681-2199 * * U.S.A. * * * * E-mail: support-asdc@earthdata.nasa.gov * * Phone: (757)864-8656 * * FAX: (757)864-8807 * * * ***************************************************************** srb_rel3.1_longwave_3hrly_199207.binary input file is opened Variable daynite_ Hour = 06 Day = 14 lon # = 100 101 102 103 104 lat band # 45 1.000 1.000 1.000 1.000 1.000 lat band # 46 1.000 1.000 1.000 1.000 1.000 lat band # 47 1.000 1.000 1.000 1.000 1.000 lat band # 48 1.000 1.000 1.000 1.000 1.000 lat band # 49 1.000 1.000 1.000 1.000 1.000 lat band # 50 1.000 1.000 1.000 1.000 1.000 lat band # 51 1.000 1.000 1.000 1.000 1.000 file daynite_3hrly_199207.ascii has been written Variable clr_toa_up_ Hour = 06 Day = 14 lon # = 100 101 102 103 104 lat band # 45 240.281 242.069 242.069 245.189 245.189 lat band # 46 244.500 246.679 249.253 250.589 251.494 lat band # 47 250.905 253.226 254.133 254.682 254.964 lat band # 48 256.507 256.579 256.846 256.986 257.350 lat band # 49 258.414 258.592 258.802 259.310 259.993 lat band # 50 260.622 260.591 261.185 262.363 263.396 lat band # 51 262.874 263.752 264.714 266.368 267.508 file clr_toa_up_3hrly_199207.ascii has been written Variable clr_sfc_up_ Hour = 06 Day = 14 lon # = 100 101 102 103 104 lat band # 45 353.652 354.267 354.267 354.576 354.576 lat band # 46 358.097 358.122 358.270 358.314 357.913 lat band # 47 362.008 361.735 361.592 361.303 360.583 lat band # 48 364.498 364.115 363.757 363.186 362.334 lat band # 49 366.010 365.756 365.313 364.633 363.901 lat band # 50 367.753 367.797 367.506 366.954 366.517 lat band # 51 369.956 370.400 370.440 370.254 370.241 file clr_sfc_up_3hrly_199207.ascii has been written Variable clr_sfc_down_ Hour = 06 Day = 14 lon # = 100 101 102 103 104 lat band # 45 276.191 280.274 280.274 279.927 279.927 lat band # 46 283.525 282.858 281.696 280.624 279.804 lat band # 47 284.484 283.048 281.900 281.084 280.297 lat band # 48 284.935 283.926 283.898 282.981 281.975 lat band # 49 287.017 286.431 285.606 284.561 282.776 lat band # 50 288.540 286.697 285.764 284.581 283.399 lat band # 51 287.916 286.486 285.103 284.102 283.406 file clr_sfc_down_3hrly_199207.ascii has been written Variable toa_up_ Hour = 06 Day = 14 lon # = 100 101 102 103 104 lat band # 45 159.160 144.281 144.281 166.914 166.914 lat band # 46 158.646 154.003 177.925 183.555 197.353 lat band # 47 234.818 238.107 219.529 192.144 173.403 lat band # 48 240.041 207.717 168.884 171.972 169.141 lat band # 49 194.644 167.088 181.239 191.057 193.637 lat band # 50 180.981 185.219 209.293 185.553 213.443 lat band # 51 244.872 223.532 220.893 222.663 245.977 file toa_up_3hrly_199207.ascii has been written Variable sfc_up_ Hour = 06 Day = 14 lon # = 100 101 102 103 104 lat band # 45 354.182 354.662 354.662 355.124 355.124 lat band # 46 358.574 358.544 359.067 358.883 358.626 lat band # 47 362.967 362.641 362.538 362.030 361.034 lat band # 48 365.455 365.000 364.258 363.563 362.720 lat band # 49 366.885 366.162 365.843 365.165 364.480 lat band # 50 368.342 368.355 368.041 367.527 367.106 lat band # 51 370.905 371.342 371.282 371.257 371.240 file sfc_up_3hrly_199207.ascii has been written Variable sfc_down_ Hour = 06 Day = 14 lon # = 100 101 102 103 104 lat band # 45 312.925 307.631 307.631 317.876 317.876 lat band # 46 316.690 312.018 337.011 319.919 329.020 lat band # 47 351.081 345.611 347.409 331.261 311.234 lat band # 48 351.214 345.114 318.466 308.846 308.414 lat band # 49 347.628 314.396 322.268 321.387 322.615 lat band # 50 329.416 325.296 322.719 323.970 323.688 lat band # 51 353.538 351.463 342.864 352.995 351.608 file sfc_down_3hrly_199207.ascii has been written 8.0 Additional Derivable Parameters The net LW flux at the top-of-atmosphere (TOA) is simply the TOA upward LW flux. The net LW flux at the surface can be defined as: Net LW Flux = Downward LW Flux - Upward LW Flux and is, therefore, generally a negative number. Net fluxes can be computed for the clear-sky and all-sky conditions. The estimates of clear-sky and all-sky fluxes also allow the estimation of the contribution by clouds to the all-sky fluxes. This is commonly referred to as the cloud radiative forcing (CRF) and is computed according to: CRF = Flux (all-sky) - Flux (clear-sky) Thus, the cloud radiative forcing on the downward longwave flux is generally positive because clouds act to increase the emission to the surface. In this way, the effect of the cloud emission on the fluxes can be estimated for each flux component. Lastly, providing TOA and surface fluxes allows one to derive the net radiative flux of the atmosphere. This is given by the relation Net Atmos. Flux = Net TOA Flux - Net Surface Flux For the LW, this flux is negative meaning that the atmosphere is cooling over the LW wavelengths. References: Bloom, Stephen, A. daSilva. D. Dee, M. Bosilovich, J-D. Chern, S. Pawson, S. Schubert, M. Sienkiewicz, I. Stajner, W-W. Tan, and M-L Wu, 2005: Documentation and Validation of the Goddard Earth Observing System (GEOS) Data Assimilation System, Version 4, NASA Technical Report,Report Number: NASA/TM-2005104606/ VOL26/VER4, Rept- 2005-01264-0/VOL26/VER4 Fu, Qiang, K. N. Liou, M. C. Cribb, T. P. Charlock, and A. Grossman, 1997: Multiple Scattering Parameterization in Thermal Infrared Radiative Transfer. J. Atmos. Sci. , Vol. 54, 2799-2812, doi: 10.1175/1520-0469(1997)054<2799:MSPITI>2.0.CO;2 Fu, Qiang, P. Yang, and W. B. Sun, 1998: An Accurate Parameterization of the Infrared Radiative Properties of Cirrus Clouds for Climate Models. J. Climate, Vol. 11, 2223-2237, doi: 10.1175/1520-0442(1998)011<2223:AAPOTI>2.0.CO;2 Kratz, David P. and Rose, Fred G., 1999: Accounting for Molecular Absorption Within the Spectral Range of the CERES Window Channel. J. Quant. Spectrosc. Radiat. Transfer, Vol. 61, 83-95. Rossow, William B. and R. A. Schiffer, 1999: Advances in Understanding Clouds from ISCCP. BAMS, Vol. 80, 2261-2287, doi: 10.1175/1520-0477(1999)080<2261:AIUCFI>2.0.CO;2. Rutan, D., F. Rose, M. Roman, N. Manalo-Smith, C. Schaaf, and T. Charlock (2009), Development and assessment of broadband surface albedo from Clouds and the Earth's Radiant Energy System Clouds and Radiation Swath data product, J. Geophys. Res., 114, D08125, doi:10.1029/2008JD010669. Wilber, Anne C., Kratz, D. P., Gupta, S. K., 1999: Surface Emissivity Maps for Use in Satellite Retrievals of Longwave Radiation, NASA Technical Report, Report Number: L-17861, NAS 1.60:209362, NASA/TP-1999-209362.