CERES Aqua Edition1B SSF Data Quality Summary Surface Fluxes - Accuracy and Validation

One of the principal objectives for the CERES data products is to provide improved estimates of surface fluxes (net and downward) for shortwave (SW) and longwave (LW) radiation. To achieve this objective, considerable effort has been focused upon obtaining consistent fluxes at the surface, within the atmosphere, and at the top of the atmosphere, all of which are produced as part of the CERES CRS data product using the SSF as input data. Validated CRS surface fluxes, however, are just now becoming available. Thus, a second effort was initiated which uses much simpler algorithms either:

• to directly tie surface fluxes to broadband CERES TOA fluxes such as in Li et al. (1993) and Darnell et al. (1992) for SW fluxes, and Inamdar and Ramanathan (1997) for clear-sky LW surface fluxes.

• or to use simple radiative parameterizations (Gupta 1989 and Gupta, Darnell, and Wilber 1992) to estimate surface fluxes, especially for the case of surface downward LW fluxes which are effectively decoupled from the TOA fluxes for cloudy sky conditions.

Consequently, these simpler SSF surface flux parameterizations are more comparable to results used in past analyses of surface radiation data sets based on ERBE or geostationary data. In general, however, they are not expected to be as precise as the CERES CRS surface fluxes, though they do represent an independent method to get to the more difficult surface flux estimates.

The CERES SSF data product provides 4 surface flux algorithm results:

1. Shortwave Flux Model A, Daytime only, Clear-sky only
   • Net surface fluxes use Li et al. (1993).
   • Downward surface fluxes use Li et al. (1993) for net and Li and Garand (1994) for surface albedo.
   
   
2. Shortwave Flux Model B, Daytime only, Clear and All-sky
   • Net and downward surface fluxes use the Langley Parameterized Shortwave Algorithm (LPSA) (Darnell et al. 1992; Gupta et al. 1999).
   
   
3. Longwave Flux Model A, Daytime and Nighttime, Clear-sky only
   • Net and downward surface fluxes uses Inamdar and Ramanathan (1997).
   
   
4. Longwave Flux Model B, Daytime and Nighttime, Clear and All-sky
   • Net and downward surface fluxes use the Langley Parameterized Longwave Algorithm (LPLA) (Gupta 1989 and Gupta, Darnell, and Wilber 1992).

For Aqua surface fluxes, clear-sky conditions are defined for CERES footprints with an imager determined cloud cover percentage less than 0.1%. Thus, to be consistent with the angular distribution models, our validation effort has also taken clear-sky to be defined as a CERES footprint with an imager determined cloud cover percentage less than 0.1%. The SSF surface fluxes are being validated using both theoretical analyses and simultaneous matching of satellite data to a range of surface sites. Preliminary results are discussed in the sections which follow.

The CERES SSF surface flux estimates are derived using the Aqua data starting with July 2000 and running through March 2005. The coincident surface fluxes are nominally gathered from the Atmospheric Radiation Measurement (ARM) networks which include the Southern Great Plains (SGP), Tropical Western Pacific (TWP) and North Slope Alaska (NSA) sites, the Climate Modeling and Diagnostic Laboratory (CMDL) network, the Baseline Surface Radiation Network (BSRN) and the Surface Radiation Budget Network (SURFRAD). Unless otherwise noted, surface site fluxes are 1 minute averages and are compared to the CERES footprint which includes the surface site.

The validation results reported in this data quality statement compare Aqua Edition 1B.

Clear-sky Shortwave Downward Flux Validation: Model A and B

For the shortwave, two models have been used to produce the surface fluxes. Both of these shortwave models are part of our validation effort; however, Model A currently produces fluxes only for clear-sky conditions while Model B produces fluxes for both clear and all-sky conditions. When the column ozone exceeds 500 DU, Model B net and downward SW surface flux values are not computed. Instead they are set to the CERES fill value.

Validation studies of the TRMM Edition 2B surface fluxes demonstrated that shortwave Model A overestimated surface insolation at the ARM Central Facility by approximately 30 Wm^-2. Considering that such biases were not observed for pristine high-latitude surface sites, it was hypothesized that the effects of aerosols could be the cause. Thus, an aerosol correction factor based on the Masuda et al. (1995) method and using the GFDL climatological aerosols (Haywood et al., 1999) was incorporated into shortwave Model A. The use of the Masuda et al. (1995) method with the GFDL climatological aerosols was shown earlier to produce a significant improvement to shortwave Model A.

Unlike previous versions of the SSF Data Quality Summary, this version groups together surface sites with similar characteristics: Continental, Desert, Coastal, Island and Polar, rather than grouping together surface sites from a single source. This will alow for a better understanding of which surface and climatological types are the most problematic.

The following table for the clear-sky cases compares shortwave Models A and B to the surface measured fluxes. Biases are defined to be CERES derived surface fluxes minus surface measured fluxes.

Results are also presented for the all-sky Model B case. To reduce the considerable variance introduced by broken cloud fields, the surface data is averaged over the 60 minutes centered on the time of the satellite overpass. Note, the variance introduced by broken cloud fields is far greater than that introduced by the temporal averaging.

Clear-sky Longwave Downward Flux Validation: Model A

Longwave Model A uses CERES-derived window and non-window TOA fluxes as well as the meteorological profiles to obtain surface fluxes for clear sky conditions. Biases are defined to be CERES derived surface fluxes minus surface measured fluxes.

Theoretical studies and validation studies employing data from Central Equatorial Pacific Experiment (CEPEX), reported by Inamdar and Ramanathan (1997), are consistent with our results. The parameterization over land surfaces was initially developed using a limited set of emissivity data available from IRIS measurements aboard NIMBUS 4 (Prabhakara and Dalu 1976). The current version of longwave Model A, however, was developed using the global emissivity maps developed by Wilber et al. (1999) and thus can be applied to the extra-tropics as well as to the tropics. Other possible sources of errors include:

1. Specification of the true radiating temperature (especially land surfaces);

2. Errors in scene identification;

3. Emissions from aerosols in the boundary layer. For instance, Inamdar and Ramanathan (1997) noted that sensitivity studies had revealed that thick haze in the boundary layer (visibilities less than 15 km) could increase the downward emissions by about 3 - 5 W m^-2.

All-sky Longwave Downward Flux Validation: Model B

Longwave Model B uses the meteorological profiles and CERES MODIS-derived cloud properties, but not the CERES-derived TOA fluxes, to obtain surface fluxes for clear and all-sky conditions. Biases are defined to be CERES derived surface fluxes minus surface measured fluxes.