Showing papers in "Journal of Geophysical Research in 1997"

Radiative transfer for inhomogeneous atmospheres: RRTM, a validated correlated-k model for the longwave

Abstract: A rapid and accurate radiative transfer model (RRTM) for climate applications has been developed and the results extensively evaluated. The current version of RRTM calculates fluxes and cooling rates for the longwave spectral region (10–3000 cm−1) for an arbitrary clear atmosphere. The molecular species treated in the model are water vapor, carbon dioxide, ozone, methane, nitrous oxide, and the common halocarbons. The radiative transfer in RRTM is performed using the correlated-k method: the k distributions are attained directly from the LBLRTM line-by-line model, which connects the absorption coefficients used by RRTM to high-resolution radiance validations done with observations. Refined methods have been developed for treating bands containing gases with overlapping absorption, for the determination of values of the Planck function appropriate for use in the correlated-k approach, and for the inclusion of minor absorbing species in a band. The flux and cooling rate results of RRTM are linked to measurement through the use of LBLRTM, which has been substantially validated with observations. Validations of RRTM using LBLRTM have been performed for the midlatitude summer, tropical, midlatitude winter, subarctic winter, and four atmospheres from the Spectral Radiance Experiment campaign. On the basis of these validations the longwave accuracy of RRTM for any atmosphere is as follows: 0.6 W m−2 (relative to LBLRTM) for net flux in each band at all altitudes, with a total (10–3000 cm−1) error of less than 1.0 W m−2 at any altitude; 0.07 K d−1 for total cooling rate error in the troposphere and lower stratosphere, and 0.75 K d−1 in the upper stratosphere and above. Other comparisons have been performed on RRTM using LBLRTM to gauge its sensitivity to changes in the abundance of specific species, including the halocarbons and carbon dioxide. The radiative forcing due to doubling the concentration of carbon dioxide is attained with an accuracy of 0.24 W m−2, an error of less than 5%. The speed of execution of RRTM compares favorably with that of other rapid radiation models, indicating that the model is suitable for use in general circulation models.

Precise point positioning for the efficient and robust analysis of GPS data from large networks

Abstract: Networks of dozens to hundreds of permanently operating precision Global Positioning System (GPS) receivers are emerging at spatial scales that range from 10(exp 0) to 10(exp 3) km. To keep the computational burden associated with the analysis of such data economically feasible, one approach is to first determine precise GPS satellite positions and clock corrections from a globally distributed network of GPS receivers. Their, data from the local network are analyzed by estimating receiver- specific parameters with receiver-specific data satellite parameters are held fixed at their values determined in the global solution. This "precise point positioning" allows analysis of data from hundreds to thousands of sites every (lay with 40-Mflop computers, with results comparable in quality to the simultaneous analysis of all data. The reference frames for the global and network solutions can be free of distortion imposed by erroneous fiducial constraints on any sites.

A finite-volume, incompressible Navier Stokes model for studies of the ocean on parallel computers

Abstract: The numerical implementation of an ocean model based on the incompressible Navier Stokes equations which is designed for studies of the ocean circulation on horizontal scales less than the depth of the ocean right up to global scale is described. A "pressure correction" method is used which is solved as a Poisson equation for the pressure field with Neumann boundary conditions in a geometry as complicated as that of the ocean basins. A major objective of the study is to make this inversion, and hence nonhydrostatic ocean modeling, efficient on parallel computers. The pressure field is separated into surface, hydrostatic, and nonhydrostatic components. First, as in hydrostatic models, a two-dimensional problem is inverted for the surface pressure which is then made use of in the three-dimensional inversion for the nonhydrostatic pressure. Preconditioned conjugate-gradient iteration is used to invert symmetric elliptic operators in both two and three dimensions. Physically motivated preconditioners are designed which are efficient at reducing computation and minimizing communication between processors. Our method exploits the fact that as the horizontal scale of the motion becomes very much larger than the vertical scale, the motion becomes more and more hydrostatic and the three- dimensional Poisson operator becomes increasingly anisotropic and dominated by the vertical axis. Accordingly, a preconditioner is used which, in the hydrostatic limit, is an exact integral of the Poisson operator and so leads to a single algorithm that seamlessly moves from nonhydrostatic to hydrostatic limits. Thus in the hydrostatic limit the model is "fast," competitive with the fastest ocean climate models in use today based on the hydrostatic primitive equations. But as the resolution is increased, the model dynamics asymptote smoothly to the Navier Stokes equations and so can be used to address small- scale processes. A "finite-volume" approach is employed to discretize the model in space in which property fluxes are defined normal to faces that delineate the volumes. The method makes possible a novel treatment of the boundary in which cells abutting the bottom or coast may take on irregular shapes and be "shaved" to fit the boundary. The algorithm can conveniently exploit massively parallel computers and suggests a domain decomposition which allocates vertical columns of ocean to each processing unit. The resulting model, which can handle arbitrarily complex geometry, is efficient and scalable and has been mapped on to massively parallel multiprocessors such as the Connection Machine (CM5) using data-parallel FORTRAN and the Massachusetts Institute of Technology data-flow machine MONSOON using the implicitly parallel language Id. Details of the numerical implementation of a model which has been designed for the study of dynamical processes in the ocean from the convective, through the geostrophic eddy, up to global scale are set out. The "kernel" algorithm solves the incompressible Navier Stokes equations on the sphere, in a geometry as complicated as that of the ocean basins with ir- regular coastlines and islands. (Here we use the term "Navier Stokes" to signify that the full nonhydrostatic equations are being employed; it does not imply a particular constitutive relation. The relevant equations for modeling the full complex- ity of the ocean include, as here, active tracers such as tem- perature and salt.) It builds on ideas developed in the compu- tational fluid community. The numerical challenge is to ensure that the evolving velocity field remains nondivergent. Most

Radiative forcing and climate response

Abstract: We examine the sensitivity of a climate model to a wide range of radiative forcings, including changes of solar irradiance, atmospheric CO2, O3, CFCs, clouds, aerosols, surface albedo, and a “ghost” forcing introduced at arbitrary heights, latitudes, longitudes, seasons, and times of day. We show that, in general, the climate response, specifically the global mean temperature change, is sensitive to the altitude, latitude, and nature of the forcing; that is, the response to a given forcing can vary by 50% or more depending upon characteristics of the forcing other than its magnitude measured in watts per square meter. The consistency of the response among different forcings is higher, within 20% or better, for most of the globally distributed forcings suspected of influencing global mean temperature in the past century, but exceptions occur for certain changes of ozone or absorbing aerosols, for which the climate response is less well behaved. In all cases the physical basis for the variations of the response can be understood. The principal mechanisms involve alterations of lapse rate and decrease (increase) of large-scale cloud cover in layers that are preferentially heated (cooled). Although the magnitude of these effects must be model-dependent, the existence and sense of the mechanisms appear to be reasonable. Overall, we reaffirm the value of the radiative forcing concept for predicting climate response and for comparative studies of different forcings; indeed, the present results can help improve the accuracy of such analyses and define error estimates. Our results also emphasize the need for measurements having the specificity and precision needed to define poorly known forcings such as absorbing aerosols and ozone change. Available data on aerosol single scatter albedo imply that anthropogenic aerosols cause less cooling than has commonly been assumed. However, negative forcing due to the net ozone change since 1979 appears to have counterbalanced 30–50% of the positive forcing due to the increase of well-mixed greenhouse gases in the same period. As the net ozone change includes halogen-driven ozone depletion with negative radiative forcing and a tropospheric ozone increase with positive radiative forcing, it is possible that the halogen-driven ozone depletion has counterbalanced more than half of the radiative forcing due to well-mixed greenhouse gases since 1979.

Marine gravity anomaly from Geosat and ERS 1 satellite altimetry

Abstract: Closely spaced satellite altimeter profiles collected during the Geosat Geodetic Mission (-6 km) and the ERS 1 Geodetic Phase (8 km) are easily converted to grids of vertical gravity gradient and gravity anomaly. The long-wavelength radial orbit error is suppressed below the noise level of the altimeter by taking the along-track derivative of each profile. Ascending and descending slope profiles are then interpolated onto separate uniform grids. These four grids are combined to form comparable grids of east and north vertical deflection using an iteration scheme that interpolates data gaps with minimum curvature. The vertical gravity gradient is calculated directly from the derivatives of the vertical deflection grids, while Fourier analysis is required to construct gravity anomalies from the two vertical deflection grids. These techniques are applied to a combination of high-density data from the dense mapping phases of Geosat and ERS 1 along with lower-density but higher-accuracy profiles from their repeat orbit phases. A comparison with shipboard gravity data shows the accuracy of the satellite- derived gravity anomaly is about 4-7 mGal for random ship tracks. The accuracy improves to 3 mGal when the ship track follows a Geosat Exact Repeat Mission track line. These data provide the first view of the ocean floor structures in many remote areas of the Earth. Some applications include inertial navigation, prediction of seafloor depth, planning shipboard surveys, plate tectonics, isostasy of volcanoes and spreading ridges, and petroleum exploration.

Operational remote sensing of tropospheric aerosol over land from EOS moderate resolution imaging spectroradiometer

Abstract: Daily distribution of the aerosol optical thickness and columnar mass concentration will be derived over the continents, from the EOS moderate resolution imaging spectroradiometer (MODIS) using dark land targets. Dark land covers are mainly vegetated areas and dark soils observed in the red and blue channels; therefore the method will be limited to the moist parts of the continents (excluding water and ice cover). After the launch of MODIS the distribution of elevated aerosol concentrations, for example, biomass burning in the tropics or urban industrial aerosol in the midlatitudes, will be continuously monitored. The algorithm takes advantage of the MODIS wide spectral range and high spatial resolution and the strong spectral dependence of the aerosol opacity for most aerosol types that result in low optical thickness in the mid-IR (2.1 and 3.8 pm). The main steps of the algorithm are (1) identification of dark pixels in the mid-IR; (2) estimation of their reflectance at 0.47 and 0.66 pm; and (3) derivation of the optical thickness and mass concentration of the accumulation mode from the detected radiance. To differentiate between dust and aerosol dominated by accumulation mode particles, for example, smoke or sulfates, ratios of the aerosol path radiance at 0.47 and 0.66 pm are used. New dynamic aerosol models for biomass burning aerosol, dust and aerosol from industrial/urban origin, are used to determine the aerosol optical properties used in the algorithm. The error in the retrieved aerosol optical thicknesses, r,, is expected to be AT, = 0.05 5 0.27,. Daily values are stored on a resolution of 10 X 10 pixels (1 km nadir resolution). Weighted and gridded 8-day and monthly composites of the optical thickness, the aerosol mass concentration and spectral radiative forcing are generated for selected scattering angles to increase the accuracy. The daily aerosol information over land and oceans (Tunr& et al., this issue), combined with continuous aerosol remote sensing from the ground, will be used to study aerosol climatology, to monitor the sources and sinks of specific aerosol types, and to study the interaction of aerosol with water vapor and clouds and their radiative forcing of climate. The aerosol information will also be used for atmospheric corrections of remotely sensed surface reflectance. In this paper, examples of applications and validations are provided.

Hydrostatic, quasi‐hydrostatic, and nonhydrostatic ocean modeling

Abstract: Ocean models based on consistent hydrostatic, quasi-hydrostatic, and nonhydrostatic equation sets are formulated and discussed. The quasi-hydrostatic and nonhydrostatic sets are more accurate than the widely used hydrostatic primitive equations. Quasi-hydrostatic models relax the precise balance between gravity and pressure gradient forces by including in a consistent manner cosine-of-latitude Coriolis terms which are neglected in primitive equation models. Nonhydrostatic models employ the full incompressible Navier Stokes equations; they are required in the study of small-scale phenomena in the ocean which are not in hydrostatic balance. We outline a solution strategy for the Navier Stokes model on the sphere that performs efficiently across the whole range of scales in the ocean, from the convective scale to the global scale, and so leads to a model of great versatility. In the hydrostatic limit the Navier Stokes model involves no more computational effort than those models which assume strict hydrostatic balance on all scales. The strategy is illustrated in simulations of laboratory experiments in rotating convection on scales of a few centimeters, simulations of convective and baroclinic instability of the mixed layer on the 1- to 10-km scale, and simulations of the global circulation of the ocean.

Observing Earth's atmosphere with radio occultation measurements using the Global Positioning System

Abstract: The implementation of the Global Positioning System (GPS) network of satellites and the development of small, high-performance instrumentation to receive GPS signals have created an opportunity for active remote sounding of the Earth's atmosphere by radio occultation at comparatively low cost. A prototype demonstration of this capability has now been provided by the GPS/MET investigation. Despite using relatively immature technology, GPS/MET has been extremely successful [Ware et al., 1996; Kursinski et al., 1996], although there is still room for improvement. The aim of this paper is to develop a theoretical estimate of the spatial coverage, resolution, and accuracy that can be expected for atmospheric profiles derived from GPS occultations. We consider observational geometry, attenuation, and diffraction in defining the vertical range of the observations and their resolution. We present the first systematic, extensive error analysis of the spacecraft radio occultation technique using a combination of analytical and simulation methods to establish a baseline accuracy for retrieved profiles of refractivity, geopotential, and temperature. Typically, the vertical resolution of the observations ranges from 0.5 km in the lower troposphere to 1.4 km in the middle atmosphere. Results indicate that useful profiles of refractivity can be derived from ∼60 km altitude to the surface with the exception of regions less than 250 m in vertical extent associated with high vertical humidity gradients. Above the 250 K altitude level in the troposphere, where the effects of water are negligible, sub-Kelvin temperature accuracy is predicted up to ∼40 km depending on the phase of the solar cycle. Geopotential heights of constant pressure levels are expected to be accurate to ∼10 m or better between 10 and 20 km altitudes. Below the 250 K level, the ambiguity between water and dry atmosphere refractivity becomes significant, and temperature accuracy is degraded. Deep in the warm troposphere the contribution of water to refractivity becomes sufficiently large for the accurate retrieval of water vapor given independent temperatures from weather analyses [Kursinski et al., 1995]. The radio occultation technique possesses a unique combination of global coverage, high precision, high vertical resolution, insensitivity to atmospheric particulates, and long-term stability. We show here how these properties are well suited for several applications including numerical weather prediction and long-term monitoring of the Earth's climate.

Element transport from slab to volcanic front at the Mariana arc

Abstract: We present a comprehensive geochemical data set for the most recent volcanics from the Mariana Islands, which provides new constraints on the timing and nature of fluxes from the subducting slab. The lavas display many features typical of island arc volcanics, with all samples showing large negative niobium anomalies and enrichments in alkaline earth elements and lead (e.g., high Ba/La and Pb/Ce). Importantly, many of these key ratios correlate with a large range in 238U excesses, (238U/230Th) = 0.97–1.56. Geochemical features show island to island variations; lavas from Guguan have the largest 238U-excesses, Pb/Ce and Ba/La ratios, while Agrigan lavas have small 238U excesses, the least radiogenic 143Nd/144Nd, and the largest negative cerium and niobium anomalies. These highly systematic variations enable two discrete slab additions to the subarc mantle to be identified. The geochemical features of the Agrigan lavas are most consistent with a dominant subducted sediment contribution. The added sedimentary component is not identical to bulk subducted sediment and notably shows a marked enrichment of Th relative to Nb. This is most readily explained by melt fractionation of the sediment with residual rutile and transfer of sedimentary material as a melt phase. For most of the highly incompatible elements, the sedimentary contribution dominates the total elemental budgets of the lavas. The characteristics best exemplified by the Guguan lavas are attributed to a slab-derived aqueous fluid phase, and Pb and Sr isotope compositions point toward the subducted, altered oceanic crust as a source of this fluid. Variable addition of the sedimentary component, but near-constant aqueous fluid flux along arc strike, can create the compositional trends observed in the Mariana lavas. High field strength element ratios (Ta/Nb and Zr/Nb) of the sediment poor Guguan lavas are higher than those of most mid-oceanic ridge basalts and suggest a highly depleted subarc mantle prior to any slab additions. The 238U-230Th systematics indicate >350 kyr between sediment and mantle melting but <30 kyr between slab dehydration and eruption of the lavas. This necessitates rapid magma migration rates and suggests that the aqueous fluid itself may trigger major mantle melting.

Global distribution of UV-absorbing aerosols from Nimbus 7/TOMS data

Abstract: Global distributions of UV-absorbing aerosols are obtained using measured differences between the 340 and the 380 nm radiances from the Nimbus 7 Total Ozone Mapping Spectrometer (TOMS) for the years 1979–1993. Time series are shown for major sources of biomass burning and desert dust giving the frequency of occurrence and areal coverage over land and oceans. Minor sources of UV-absorbing aerosols in the atmosphere are also discussed (volcanic ash and oil fires). Relative values of year-to-year variability of UV-absorbing aerosol amounts are shown for major aerosol source regions: (1) central South America (Brazil) near 10°S latitude; (2) Africa near 0°–20°S and 0° to 10°N latitude; (3) Saharan Desert and sub-Saharan region (Sahel), Arabian Peninsula, and the northern border region of India; (4) agricultural burning in Indonesia, Eastern China, and Indochina, and near the mouth of the Amazon River; and (5) coal burning and dust in northeastern China. The first three of these regions dominate the injection of UV-absorbing aerosols into the atmosphere each year and cover areas far outside of their source regions from advection of UV-absorbing particulates by atmospheric wind systems. During the peak months, smoke and dust from these sources are transported at altitudes above 1 km with an optical depth of at least 0.1 and can cover about 10% of the Earth's surface. Boundary layer absorbing aerosols are not readily seen by TOMS because the small amount of underlying Rayleigh scattering leads to a small signal. Significant portions of the observed dust originate from agricultural regions frequently within arid areas, such as in the Sahel region of Africa, especially from the dry lake-bed near Lake Chad (13.5°N, 14°E), and intermittently dry drainage areas and streams. In addition to drought cycle effects, this suggests there may be an anthropogenic component to the amount of dust injected into the atmosphere each year. Detection of absorbing aerosols and calculation of optical depths are affected by the presence of large-scale and subpixel clouds in the TOMS field of view.

A unified directional spectrum for long and short wind-driven waves

Abstract: Review of several recent ocean surface wave models finds that while comprehensive in many regards, these spectral models do not satisfy certain additional, but fundamental, criteria. We propose that these criteria include the ability to properly describe diverse fetch conditions and to provide agreement with in situ observations of Cox and Munk [1954] and Jahne and Riemer [1990] and Hara et al. [1994] data in the high-wavenumber regime. Moreover, we find numerous analytically undesirable aspects such as discontinuities across wavenumber limits, nonphysical tuning or adjustment parameters, and noncentrosymmetric directional spreading functions. This paper describes a two-dimensional wavenumber spectrum valid over all wavenumbers and analytically amenable to usage in electromagnetic models. The two regime model is formulated based on the Joint North Sea Wave Project (JONSWAP) in the long-wave regime and on the work of Phillips [1985] and Kitaigorodskii [1973] at the high wavenumbers. The omnidirectional and wind-dependent spectrum is constructed to agree with past and recent observations including the criteria mentioned above. The key feature of this model is the similarity of description for the high- and low-wavenumber regimes; both forms are posed to stress that the air-sea interaction process of friction between wind and waves (i.e., generalized wave age, u/c) is occurring at all wavelengths simultaneously. This wave age parameterization is the unifying feature of the spectrum. The spectrum's directional spreading function is symmetric about the wind direction and has both wavenumber and wind speed dependence. A ratio method is described that enables comparison of this spreading function with previous noncentrosymmetric forms. Radar data are purposefully excluded from this spectral development. Finally, a test of the spectrum is made by deriving roughness length using the boundary layer model of Kitaigorodskii. Our inference of drag coefficient versus wind speed and wave age shows encouraging agreement with Humidity Exchange Over the Sea (HEXOS) campaign results.

A new mechanism for regional atmospheric chemistry modeling

Abstract: A new gas-phase chemical mechanism for the modeling of regional atmospheric chemistry, the “Regional Atmospheric Chemistry Mechanism” (RACM) is presented. The mechanism is intended to be valid for remote to polluted conditions and from the Earth's surface through the upper troposphere. The RACM mechanism is based upon the earlier Regional Acid Deposition Model, version 2 (RADM2) mechanism [Stockwell et al., 1990] and the more detailed Euro-RADM mechanism [Stockwell and Kley, 1994]. The RACM mechanism includes rate constants and product yields from the most recent laboratory measurements, and it has been tested against environmental chamber data. A new condensed reaction mechanism is included for biogenic compounds: isoprene, α-pinene, and d-limonene. The branching ratios for alkane decay were reevaluated, and in the revised mechanism the aldehyde to ketone ratios were significantly reduced. The relatively large amounts of nitrates resulting from the reactions of unbranched alkenes with NO3 are now included, and the production of HO from the ozonolysis of alkenes has a much greater yield. The aromatic chemistry has been revised through the use of new laboratory data. The yield of cresol production from aromatics was reduced, while the reactions of HO, NO3, and O3 with unsaturated dicarbonyl species and unsaturated peroxynitrate are now included in the RACM mechanism. The peroxyacetyl nitrate chemistry and the organic peroxy radical-peroxy radical reactions were revised, and organic peroxy radical +NO3 reactions were added.

Remote sensing of aerosol properties over oceans using the MODIS/EOS spectral radiances

Abstract: Spectral radiances measured at the top of the atmosphere in a wide spectral range (0.55–2.13 μm) are used to monitor the aerosol optical thickness and the aerosol size distribution (integrated on the vertical column) of the ambient (undisturbed) aerosol over the oceans. Even for the moderate resolution imaging spectrometer (MODIS) wide spectral range, only three parameters that describe the aerosol loading and size distribution can be retrieved. These three parameters are not always unique. For instance, the spectral radiance of an aerosol with a bilognormal size distribution can be simulated very well with a single lognormal aerosol with an appropriate mean radius and width of distribution. Preassumptions on the general structure of the size distribution are therefore required in the inversion of MODIS data. The retrieval of the aerosol properties is performed using lookup table computations. The volume size distribution in the lookup table is described with two lognormal modes: a single mode to describe the accumulation mode particles (radius 1.0 μm). Note that two accumulation modes may be present, one dominated by gas phase processes and a second dominated by cloud phase processes. The coarse mode can also be split into several partially overlapping modes describing maritime salt particles and dust. The aerosol parameters we expect to retrieve are η, the fractional contribution of the accumulation mode to scattering; τ, the spectral optical thickness; and rm, the mean particle size of the dominant mode. Additional radiative quantities such as asymmetry parameter and effective radius are derived subsequently. The impact of the surface conditions, wind speed and chlorophyll content on the retrieval is estimated, the impact of potential sources of error like the calibration of the instrument is also tested. The algorithm has been applied successfully to actual data sets provided by the Thematic Mapper on Landsat 5 and by the MODIS airborne simulator on the ER-2 and tested against ground and airborne measurements. A first estimate of the general accuracy is Δτ = ±0.05±0.05τ (at 550 nm), Δrm = 0.3rm, Δη = ±0.25.

Leaf area index of boreal forests: theory, techniques, and measurements

Abstract: Leaf area index (LAI) is a key structural characteristic of forest ecosystems because of the role of green leaves in controlling many biological and physical processes in plant canopies. Accurate LAI estimates are required in studies of ecophysiology, atmosphere-ecosystem interactions, and global change. The objective of this paper is to evaluate LAI values obtained by several research teams using different methods for a broad spectrum of boreal forest types in support of the international Boreal Ecosystem-Atmosphere Study (BOREAS). These methods include destructive sampling and optical instruments: the tracing radiation and architecture of canopies (TRAC), the LAI-2000 plant canopy analyzer, hemispherical photography, and the Sunfleck Ceptometer. The latter three calculate LAI from measured radiation transmittance (gap fraction) using inversion models that assume a random spatial distribution of leaves. It is shown that these instruments underestimate LAI of boreal forest stands where the foliage is clumped. The TRAC quantifies the clumping effect by measuring the canopy gap size distribution. For deciduous stands the clumping index measured from TRAC includes the clumping effect at all scales, but for conifer stands it only resolves the clumping effect at scales larger than the shoot (the basic collection of needles). To determine foliage clumping within conifer shoots, a video camera and rotational light table system was used. The major difficulties in determining the surface area of small conifer needles have been largely overcome by the use of an accurate volume displacement method. Hemispherical photography has the advantage that it also provides a permanent image record of the canopies. Typically, LAI falls in the range from 1 to 4 for jack pine and aspen forests and from 1 to 6 for black spruce. Our comparative studies provide the most comprehensive set of LAI estimates available for boreal forests and demonstrate that optical techniques, combined with limited direct foliage sampling, can be used to obtain quick and accurate LAI measurements.

Atmospheric effects in interferometric synthetic aperture radar surface deformation and topographic maps

Abstract: Interferogram images derived from repeat-pass spaceborne synthetic aperture radar systems exhibit artifacts due to the time and space variations of atmospheric water vapor Other tropospheric variations, such as pressure and temperature, also induce distortions, but the effects are smaller in magnitude and more evenly distributed throughout the interferogram than the wet troposphere term Spatial and temporal changes of 20% in relative humidity lead to 10 cm errors in deformation products, and perhaps 100 m of error in derived topographic maps for those pass pairs with unfavorable baseline geometries In wet regions such as Hawaii, these are by far the dominant errors in the Spaceborne Imaging Radar-C and X Band Synthetic Aperature Radar (SIR-C/X-SAR) interferometric products The unknown time delay from tropospheric distortion is independent of frequency, and thus multiwavelength measurements, such as those commonly used to correct radar altimeter and Global Positioning System (GPS) ionospheric biases, cannot be used to rectify the error In the topographic case, the errors may be mitigated by choosing interferometric pairs with relatively long baselines, as the error amplitude is inversely proportional to the perpendicular component of the interferometer baseline For the SIR-C/X-SAR Hawaii data we found that the best (longest) baseline pair produced a map supporting 100 m contouring, whereas the poorest baseline choice yielded an extremely noisy topographic map even at this coarse contour interval In the case of deformation map errors the result is either independent of baseline parameters or else very nearly so Here the only solution is averaging of independent interferograms, so in order to create accurate deformation products in wet regions many multiple passes may be required Rules for designing optimal data acquisition and processing sequences for interferometric analyses in nondesert parts of the world are (1) to use the longest radar wavelengths possible, within ionospheric scintillation and Faraday rotation limits, (2) for topography, maximize interferometer baseline within decorrelation limits* and (3) for surface deformation, use multiple observations and average the derived products Following the above recipe yields accuracies of 10 m for digital elevation models and 1 cm for deformation maps even in very wet regions, such as Hawaii

Sensitivity of a global sea ice model to the treatment of ice thermodynamics and dynamics

Abstract: The sensitivity of a global thermodynamic-dynamic sea ice model coupled to a one-dimensional upper ocean model to degradations of the model physics is investigated. The thermodynamic component of the sea ice model takes into consideration the presence of snow on top of sea ice, the storage of sensible and latent heat inside the snow-ice system, the influence of the subgrid-scale snow and ice thickness distributions on sea ice thermodynamics, the transformation of snow into snow ice when snow depth increases to the point where the snow-ice interface sinks below the waterline, and the existence of leads and polynyas (areas of open water) within the ice cover. Ice dynamics is treated basically as by Hibler [1979]. Regarding the upper ocean model, it is made up of an integral mixed layer model and of a diffusive model of the pycnocline. Advection of heat and salt by oceanic currents is implicitly accounted for by restoring the temperatures and salinities of the water column to annual mean data. It is very important to note that a single set of parameter values is employed to simultaneously simulate the Arctic and Antarctic ice regimes. A control run carried out with the model demonstrates that it does reasonably well in simulating the seasonal waxing and waning of both ice packs. The sensitivity study focuses on physical processes pertaining to (1) the vertical growth and decay of sea ice (thermal inertia of snow and ice, heat conduction, and snow cover), (2) the lateral growth and decay of sea ice (leads and polynyas), and (3) the sea ice dynamics (ice motion and shear strength). A total of nine sensitivity experiments have been performed. Each experiment consisted of removing a particular parameterization from the control run computer code. It appears that the thermal inertia of the snow-ice system is negligible in the Antarctic but not in the Arctic, where the total heat content of sea ice is chiefly dictated by internal storage of latent heat in brine pockets, sensible heat storage being of very minor consequence. It is also found that the inclusion of a prognostic snow layer and of a scheme of snow ice formation is important for sea ice modeling in the southern hemisphere. Furthermore, our results suggest that the thermodynamic effect of the subgrid-scale snow and ice thickness distributions, the existence of open water areas within the ice cover, and the ice motion play a crucial role in determining the seasonal behavior of both ice packs. The ice shear strength seems to be of lesser importance, although it has a nonnegligible effect in both hemispheres. We can therefore conclude that all these processes should be represented in global climate models.

Oxygen 18/16 variability in Greenland snow and ice with 10 -3- to 105-year time resolution

Abstract: The 3-km-long Greenland Ice Sheet Project 2 (GISP2) ice core presents a 100,000 +- year detailed oxygen isotope profile covering almost a full glacial-interglacial cycle. Measuranents of isotopic fluctuations in snow, frost, and atmospheric water vapor samples collected during summer field seasons (up to 20%0) are compatible with the large and abrupt 80/160 changes observed in accumulated tim. Snow pit 1580 profiles from the GISP2 summit area, however, show rapid smoothing of the 180/160 signal near the surface. Beyond about 2-m depth the smoothedi5180 signal is fairly well preserved and can be interpreted in terms of average local weather conditions and climate. The longer climate fluctuations also have regional and often global significance. In the older part of the record, corresponding to marine isotope stages (MIS) 5a to 5d, the effect of orbital climate forcing via the 19- and 23-kyr precession cycles and the 41-kyr obliquity cycle is obvious. From the end ofMIS 5a, at about 75,000 years B.P., till the end of the glacial at the Younger Dryas-Preboreal transition, at 11,650 years B.P., the O180/160 record shows frequent, rapid switches between intermediate interstadial and low stadial values. Fourier spectra of the oscillations that are superimposed on the orbitally induced changes contain a strong periodicity at 1.5 kyr, a broad peak at 4.0 kyr, and additional shorter periods. Detailed comparison of the GISP2 180/160 record with the Vostok, Antarctica, 15D record; Pacific Ocean foraminiferal 180/160; Grande Pile, France, tree pollen; and insolation indicates that a counterpart to many of the rapid 180/160 fluctuations of GISP2 can be found in the other records, and that the GISP2 isotopic changes clearly are the local expression of climate changes of worldwide extent. Correlation of events on the independent GISP2 and SPECMAP time scales for the interval 10,000-50,000 years B.P. shows excellent chronometric agreement, except possibly for the event labeled 3.1. The glacial to interglacial transition evidently started simultaneously in the Arctic and the Antarctic, but its development and its expression in Greenland isotopes was later suppressed by the influence of meltwater, especially from the Barents Sea ice sheet, on deep water formation and ocean circulation. Meltwaters from different ice sheets bordering the North Atlantic also influenced ocean circulation during the Bolling-Allerod interstadial complex and the Younger Dryas and led to a distinct development of European climate and Greenland 180/160 values. The Holocene interval with long-term stable mean isotopic values contains several fluauations with periods from years to millennia. Dominant is a 6.3-year oscillation with amplitude up to 3 to 4%0. Periodicities of 11 and 210 years, also found in the solar-modulated records of the cosmogenic isotopes 1oBe and 14C, suggest solar processes as the cause of these cycles. Depression of180/160 values (cooling) by volcanic eruptions is observed in stacked GISP21580 records, but the effect is small and not likely to trigger major climate changes.

Major features and forcing of high-latitude northern hemisphere atmospheric circulation using a 110,000-year-long glaciochemical series

Abstract: The Greenland Ice Sheet Project 2 glaciochemical series (sodium, potassium, ammonium, calcium, magnesium, sulfate, nitrate, and chloride) provides a unique view of the chemistry of the atmosphere and the history of atmospheric circulation over both the high latitudes and mid-low latitudes of the northern hemisphere. Interpretation of this record reveals a diverse array of environmental signatures that include the documentation of anthropogenically derived pollutants, volcanic and biomass burning events, storminess over marine surfaces, continental aridity and biogenic source strength plus information related to the controls on both high- and low-frequency climate events of the last 110,000 years. Climate forcings investigated include changes in insolation of the order of the major orbital cycles that control the long-term behavior of atmospheric circulation patterns through changes in ice volume (sea level), events such as the Heinrich events (massive discharges of icebergs first identified in the marine record) that are found to operate on a 6100-year cycle due largely to the lagged response of ice sheets to changes in insolation and consequent glacier dynamics, and rapid climate change events (massive reorganizations of atmospheric circulation) that are demonstrated to operate on 1450-year cycles. Changes in insolation and associated positive feedbacks related to ice sheets may assist in explaining favorable time periods and controls on the amplitude of massive rapid climate change events. Explanation for the exact timing and global synchroneity of these events is, however, more complicated. Preliminary evidence points to possible solar variability-climate associations for these events and perhaps others that are embedded in our ice-core-derived atmospheric circulation records.

The transition from brittle faulting to cataclastic flow in porous sandstones: Mechanical deformation

Abstract: Triaxial compression experiments were conducted to investigate the inelastic and failure behavior of six sandstones with porosities ranging from 15% to 35%. A broad range of effective pressures was used so that the transition in failure mode from brittle faulting to cataclastic flow could be observed. In the brittle faulting regime, shear-induced dilation initiates in the prepeak stage at a stress level C' which increases with effective mean stress. Under elevated effective pressures, a sample fails by cataclastic flow. Strain hardening and shear-enhanced compaction initiates at a stress level C* which decreases with increasing effective mean stress. The critical stresses C' and C* were marked by surges in acoustic emission. In the stress space, C* maps out an approximately elliptical yield envelope, in accordance with the critical state and cap models. Using plasticity theory, the flow rule associated with this yield envelope was used to predict porosity changes which are comparable to experimental data. In the brittle faulting regime the associated flow rule predicts dilatancy to increase with decreasing effective pressure in qualitative agreement with the experimental observations. The data were also compared with prediction of a nonassociative model on the onset of shear localization. Experimental data suggest that a quantitative measure of brittleness is provided by the grain crushing pressure (which decreases with increasing porosity and grain size). Geologic data on tectonic faulting in siliciclastic formations (of different porosity and grain size) are consistent with the laboratory observations.

Tropical Atlantic sea surface temperature variability and its relation to El Niño‐Southern Oscillation

Abstract: Past analyses of tropical Atlantic sea surface temperature variability have suggested a dipole behavior between the northern and southern tropics, across the Intertropical Convergence Zone (ITCZ). By analyzing an improved 43-year (1950-1992) record of SST (Smith et al., 1996) and other data derived from the Comprehensive Ocean-Atmosphere Data Set (COADS), it is shown that the regions north and south of the ITCZ are statistically independent of each other at the seasonal to interannual timescales dominating the data, confirming the conclusions of Houghton and Tourre (1992). Some dipole behavior does develop weakly during the boreal spring season, when there is a tendency for SST anomaly west of Angola to be opposite of that in the tropical North Atlantic. It is further shown that tropical Atlantic SST variability is corre- lated with Pacific E1 Nifio-Southern Oscillation (ENSO) variability in several regions. The ma- jor region affected is the North Atlantic area of NE trades west of 40oW along 10oN - 20oN and extending into the Caribbean. There, about 50-80% of the anomalous SST variability is associ- ated with the Pacific ENSO, with Atlantic warmings occurring 4-5 months after the mature phases of Pacific warm events. An analysis of local surface flux fields derived from COADS data shows that the ENSO-related Atlantic warmings occur as a result of reductions in the surface NE trade wind speeds, which in turn reduce latent and sensible heat losses over the region in ques- tion, as well as cooling due to entrainment. This ENSO connection is best developed during the boreal spring following the most frequent season of maximum ENSO anomalies in the Pacific. A region of secondary covariability with ENSO occurs along the northern edge of the mean ITCZ position and appears to be associated with northward migrations of the ITCZ when the North At- lantic warmings occur. Although easterly winds are intensified in the western equatorial Atlantic in response to Pacific warm events, they do not produce strong local changes in SST. Contrary to expectations from studies based on equatorial dynamics, these teleconnected wind anomalies do not give rise to significant correlations of SST in the Gulf of Guinea with the Pacific ENSO. As the teleconnection sequence matures, strong SE trades at low southern latitudes follow the de- velopment of the North Atlantic SST anomaly and precede by several months the appearance of weak negative SST anomalies off Angola and stronger positive anomalies extending eastward from southern Brazil along 15o-30oS.

Characterization of tropospheric aerosols over the oceans with the NOAA advanced very high resolution radiometer optical thickness operational product

Abstract: The National Oceanic and Atmospheric Administration (NOAA) advanced very high resolution radiometer (AVHRR) is an instrument on a polar orbiting satellite that provides information on global aerosol distributions. The remote sensing algorithm is based on measurements of backscattered solar radiation which yield a measure of the radiatively equivalent aerosol optical thickness τ A sat (EAOT) over the oceans. Seasonally composited EAOT data for the period July 1989 to June 1991 reveal many spatially coherent plume-like patterns that can usually be interpreted in terms of known (or reasonably hypothesized) sources in association with climatological wind fields. The largest and most persistent areas of high EAOT values are associated with wind-blown dust and biomass burning sources; especially prominent are sources in Africa, the middle East, and the Asian subcontinent. Prominent plumes over the midlatitude North Atlantic are attributed to pollution emissions from North America and Europe. Large plumes attributed to pollution aerosols and dust from sources in Asia are clearly visible over the western and central North Pacific. On a global scale the annually averaged northern hemisphere EAOT values are about 1.7 times greater than those in the southern hemisphere. Considering each hemisphere separately, EAOT values in summer are about twice those in winter. Within the midlatitude band 30°-60° (i.e., where anthropogenic emissions are greatest) the summer/winter ratio is about 3. The temporal variability of monthly mean EAOT in specific ocean regions often shows characteristic seasonal patterns that are usually consistent with aerosol measurements made in the marine boundary layer. Nonetheless, there are features in the EAOT distributions that can not be readily interpreted at this time. The AVHRR EAOT distributions demonstrate that satellite products can serve as a useful tool for the planning and implementation of focused aerosol research programs and that they will be especially important in studies of climate-related processes.

A well‐calibrated ocean algorithm for special sensor microwave / imager

Abstract: I describe an algorithm for retrieving geophysical parameters over the ocean from special sensor microwave/imager (SSM/I) observations. This algorithm is based on a model for the brightness temperature T(sub B) of the ocean and intervening atmosphere. The retrieved parameters are the near-surface wind speed W, the columnar water vapor V, the columnar cloud liquid water L, and the line-of-sight wind W(sub LS). I restrict my analysis to ocean scenes free of rain, and when the algorithm detects rain, the retrievals are discarded. The model and algorithm are precisely calibrated using a very large in situ database containing 37,650 SSM/I overpasses of buoys and 35,108 overpasses of radiosonde sites. A detailed error analysis indicates that the T(sub B) model rms accuracy is between 0.5 and 1 K and that the rms retrieval accuracies for wind, vapor, and cloud are 0.9 m/s, 1.2 mm, and 0.025 mm, respectively. The error in specifying the cloud temperature will introduce an additional 10% error in the cloud water retrieval. The spatial resolution for these accuracies is 50 km. The systematic errors in the retrievals are smaller than the rms errors, being about 0.3 m/s, 0.6 mm, and 0.005 mm for W, V, and L, respectively. The one exception is the systematic error in wind speed of -1.0 m/s that occurs for observations within +/-20 deg of upwind. The inclusion of the line-of-sight wind W(sub LS) in the retrieval significantly reduces the error in wind speed due to wind direction variations. The wind error for upwind observations is reduced from -3.0 to -1.0 m/s. Finally, I find a small signal in the 19-GHz, horizontal polarization (h(sub pol) T(sub B) residual DeltaT(sub BH) that is related to the effective air pressure of the water vapor profile. This information may be of some use in specifying the vertical distribution of water vapor.

Digital isochrons of the world's ocean floor

Abstract: We have created a digital age grid of the ocean floor with a grid node interval of 6 arc min using a self-consistent set of global isochrons and associated plate reconstruction poles. The age at each grid node was determined by linear interpolation between adjacent isochrons in the direction of spreading. Ages for ocean floor between the oldest identified magnetic anomalies and continental crust were interpolated by estimating the ages of passive continental margin segments from geological data and published plate models. We have constructed an age grid with error estimates for each grid cell as a function of (1) the error of ocean floor ages identified from magnetic anomalies along ship tracks and the age of the corresponding grid cells in our age grid, (2) the distance of a given grid cell to the nearest magnetic anomaly identification, and (3) the gradient of the age grid: i.e., larger errors are associated with high age gradients at fracture zones or other age discontinuities. Future applications of this digital grid include studies of the thermal and elastic structure of the lithosphere, the heat loss of the Earth, ridge-push forces through time, asymmetry of spreading, and providing constraints for seismic tomography and mantle convection models.

A new functional form to study the solar wind control of the magnetopause size and shape

Abstract: In this study a new functional form, r = r 0 [2/(1 + cos θ)] α , is used to fit. the size and shape of the magnetopause using crossings from ISEE 1 and 2, Active Magnetospheric Particle Tracer Explorers/Ion Release Module (AMPTE/IRM), and IMP 8 satellites. This functional form has two parameters, τ 0 and α, representing the standoff distance and the level of tail flaring. The value r is the radial distance at an angle (θ) between the Sun-Earth line and the direction of τ. It is found that r 0 varies with the interplanetary magnetic field (IMF) B z component and has a break in the slope at B z = 0 nT. The best-fit value of τ 0 decreases with increasing southward IMF B z . For northward IMF B z , the best-fit value of τ 0 increases slightly with increasing B z . The best-fit value of α increases monotonically with decreasing IMF B z . The dynamic pressure (D p ) also changes τ 0 and α. The parameters D p and τ 0 are related by a power law of -1/(6.6±0.8). The best-fit value of α is slightly larger for larger dynamic pressure, which implies that D p also has a role in flux transfer from the dayside to the nightside, but the size of this effect is small. An explicit function for the size and shape of the magnetopause, in terms of D p and B z , is obtained by using multiple parameter fitting in a form that is useful for operational space applications such as predicting when satellites at geosynchronous orbit will be found in the magnetosheath.

Atmospheric correction of ocean color imagery in the Earth Observing System era

Abstract: Sensors that can be used for the observation of ocean color in NASA's Earth Observing System era (SeaWiFS, MODIS, and MISR) have been designed with 2–4 times the radiometric sensitivity of the proof-of-concept ocean color instrument CZCS (coastal zone color scanner). To realize an improvement in the retrieval of biologically important ocean parameters, e.g., the concentration of the photosynthetic pigment chlorophyll a, from this increased sensitivity, significantly better atmospheric correction than was applied to CZCS is required. Atmospheric correction improvement necessitates the inclusion of the effects of multiple scattering, which are strongly dependent on the aerosol size distribution, concentration, and absorption properties. We review the basic concepts of atmospheric correction over the oceans and provide the details of the algorithms currently being developed for SeaWiFS, MODIS, and MISR. An alternate correction algorithm that could be of significant value in the coastal zone is described for MISR. Related issues such as the influence of aerosol vertical structure in the troposphere, polarization of the light field, sea surface roughness, and oceanic whitecaps on the sea surface are evaluated and plans for their inclusion in the algorithm are described. Unresolved issues, such as the presence of stratospheric aerosol, the appropriateness of the aerosol models used in the assessment of multiple scattering, and the identification of, and difficulties associated with the correction for, the presence of absorbing aerosols, e.g., urban pollution or mineral dust, are identified, and suggestions are provided for their resolution.

Modeling phase functions for dustlike tropospheric aerosols using a shape mixture of randomly oriented polydisperse spheroids

Abstract: Laboratory and in situ measurements show that scattering properties of natural nonspherical particles can be significantly different from those of volume-or surface-equivalent spheres, thus suggesting that Mie theory may not be suitable for interpreting satellite reflectance measurements for dustlike tropospheric aerosols. In this paper we use the rigorous T-matrix method to extensively compute light scattering by shape distributions of polydisperse, randomly oriented spheroids with refractive indices and size distributions representative of naturally occurring dust aerosols. Our calculations show that even after size and orientation averaging, a single spheroidal shape always produces a unique, shape-specific phase function distinctly different from those produced by other spheroidal shapes. However, phase functions averaged over a wide aspect-ratio distribution of prolate and oblate spheroids are smooth, featureless, and nearly flat at side-scattering angles and closely resemble those measured for natural soil and dust particles. Thus, although natural dust particles are, of course, not perfect spheroids, they are always mixtures of highly variable shapes, and their phase function can be adequately modeled using a wide aspect-ratio distribution of prolate and oblate spheroidal grains. Our comparisons of nonspherical versus projected-area-equivalent spherical particles show that spherical-nonspherical differences in the scattering phase function can be large and therefore can cause significant errors in the retrieved aerosol optical thickness if Mie theory is used to analyze reflectance measurements of nonspherical aerosols. On the other hand, the differences in the total optical cross sections, single-scattering albedo, asymmetry parameter of the phase function, and backscattered fraction are much smaller and in most cases do not exceed 10%. This may suggest that for a given aerosol optical thickness the influence of particle shape on the aerosol radiative forcing is negligibly small. Spherical-nonspherical differences in the extinction-to-backscatter ratio are very large and should be explicitly taken into account in inverting lidar measurements of dustlike aerosols.

Global Positioning System measurements of present-day crustal movements in the Arabia-Africa-Eurasia plate collision zone

Abstract: We present and interpret Global Positioning System (GPS) measurements of crustal motions for the period 1988–1994 at 54 sites extending east-west from the Caucasus mountains of southern Russia, Georgia, and Armenia to the Aegean coast of Turkey and north-south from the southern edge of the Eurasian plate (Pontus block) to the northern edge of the Arabian platform. Viewed from a Eurasia-fixed reference frame, sites on the northern Arabian platform move N38±13°W at 20±3 mm/yr, roughly consistent with the velocity implied by NUVEL 1A circuit closure (N23±7°W at 24±2 mm/yr). The motion of Arabia appears to be transferred directly to the region of Turkey north of the suture. However, eastern Turkey is characterized by distributed deformation while central/western Turkey is characterized by coherent plate motion involving westward displacement and counterclockwise rotation of the Anatolian plate. Internal deformation within the central part of the Anatolian plate is less than 2 mm/yr. The Anatolian plate is decoupled from Eurasia along the right-lateral, strike-slip North Anatolian fault (NAF). This different response in eastern and western Turkey to the collision of Arabia may result from the different boundary conditions, the Hellenic arc forming a “free” boundary to the west and the Asian continent and oceanic lithosphere of the Black and Caspian Seas forming a resistant boundary to the north and east. We derive a best fitting Euler vector for Anatolia-Eurasia motion of 29.2±0.8°N, 32.9±0.4°, 1.3±0.1°/m.y. The mapped surface trace of the NAF corresponds well to a small circle about this pole. The new Euler vector implies an upper bound for NAF slip rate of 30±2 mm/yr (i.e., assuming all relative motion is accommodated along the NAF). Using the NUVEL 1A Euler vector for Arabia-Eurasia and the GPS Euler vector for Anatolia-Eurasia, we determine an Arabia-Anatolia Euler vector of 31±2°N, 45±2°E, 0.9±0.1 °/m.y. and an upper bound on the East Anatolian fault slip rate of 15±3 mm/yr. The Aegean Trough region of western Turkey deviates significantly from coherent plate rotation. In addition to rotating with Anatolia, this region shows roughly N-S extension at a rate of 14±5 mm/yr. Taken together with satellite laser ranging results along the Hellenic arc, the contemporary pattern of deformation indicates increasing motions toward the arc, suggesting that the westward displacement and counterclockwise rotation of Anatolia is driven both by “pushing” from the Arabian plate and by “pulling” or basal drag associated with the foundering African plate along the Hellenic subduction zone.

Contribution of different aerosol species to the global aerosol extinction optical thickness: Estimates from model results

Abstract: We combine global distributions of aerosol loading resulting from transport models for soil dust, sulfate, sea salt, and carbonaceous aerosol. From the aerosol distributions we estimate optical thicknesses and compare them with Sun photometer measurements and satellite retrievals, thereby revealing problems with both model results and comparisons with such measurements. Globally, sulfate, dust, and carbonaceous particles appear to contribute equally to the total aerosol optical thickness. Owing to the different optical properties of different aerosol types, aerosol composition should be taken into consideration for estimating the aerosol climate effect as well as for aerosol retrievals from satellite measurements.

Atmospheric correction of visible to middle-infrared EOS-MODIS data over land surfaces: Background, operational algorithm and validation

Abstract: The NASA moderate resolution imaging spectroradiometer (MODIS) instrument will provide a global and improved source of information for the study of land surfaces with a spatial resolution of up to 250 m. Prior to the derivation of various biophysical parameters based on surface reflectances, the top of the atmosphere signals need to be radiometrically calibrated and corrected for atmospheric effects. The present paper describes in detail the state of the art techniques that will be used for atmospheric correction of MODIS bands 1 through 7, centered at 648, 858, 470, 555, 1240, 1640, and 2130 nm, respectively. Previous operational correction schemes have assumed a standard atmosphere with zero or constant aerosol loading and a uniform, Lambertian surface. The MODIS operational atmospheric correction algorithm, reported here, uses aerosol and water vapor information derived from the MODIS data, corrects for adjacency effects and takes into account the directional properties of the observed surface. This paper also describes the operational implementation of these techniques and its optimization. The techniques are applied to remote sensing data from the Landsat Thematic Mapper (TM), the NOAA advanced very high resolution radiometer (AVHRR), and the MODIS airborne simulator (MAS) and validated against ground-based measurements from the Aerosol Robotic Network (AERONET).

Validity of the temperature reconstruction from water isotopes in ice cores

Abstract: Well-documented present-day distributions of stable water isotopes (HDO and H218O) show the existence, in middle and high latitudes, of a linear relationship between the mean annual isotope content of precipitation (δD and δ18O) and the mean annual temperature at the precipitation site. Paleoclimatologists have used this relationship, which is particularly well obeyed over Greenland and Antarctica, to infer paleotemperatures from ice core data. There is, however, growing evidence that spatial and temporal isotope/surface temperature slopes differ, thus complicating the use of stable water isotopes as paleothermometers. In this paper we review empirical estimates of temporal slopes in polar regions and relevant information that can be inferred from isotope models: simple, Rayleigh-type distillation models and (particularly over Greenland) general circulation models (GCMs) fitted with isotope tracer diagnostics. Empirical estimates of temporal slopes appear consistently lower than present-day spatial slopes and are dependent on the timescale considered. This difference is most probably due to changes in the evaporative origins of moisture, changes in the seasonality of the precipitation, changes in the strength of the inversion layer, or some combination of these changes. Isotope models have not yet been used to evaluate the relative influences of these different factors. The apparent disagreement in the temporal and spatial slopes clearly makes calibrating the isotope paleothermometer difficult. Nevertheless, the use of a (calibrated) isotope paleothermometer appears justified; empirical estimates and most (though not all) GCM results support the practice of interpreting ice core isotope records in terms of local temperature changes.