Showing papers in "Journal of Geophysical Research in 2003"

Global analyses of sea surface temperature, sea ice, and night marine air temperature since the late nineteenth century

Abstract: [1] We present the Met Office Hadley Centre's sea ice and sea surface temperature (SST) data set, HadISST1, and the nighttime marine air temperature (NMAT) data set, HadMAT1. HadISST1 replaces the global sea ice and sea surface temperature (GISST) data sets and is a unique combination of monthly globally complete fields of SST and sea ice concentration on a 1° latitude-longitude grid from 1871. The companion HadMAT1 runs monthly from 1856 on a 5° latitude-longitude grid and incorporates new corrections for the effect on NMAT of increasing deck (and hence measurement) heights. HadISST1 and HadMAT1 temperatures are reconstructed using a two-stage reduced-space optimal interpolation procedure, followed by superposition of quality-improved gridded observations onto the reconstructions to restore local detail. The sea ice fields are made more homogeneous by compensating satellite microwave-based sea ice concentrations for the impact of surface melt effects on retrievals in the Arctic and for algorithm deficiencies in the Antarctic and by making the historical in situ concentrations consistent with the satellite data. SSTs near sea ice are estimated using statistical relationships between SST and sea ice concentration. HadISST1 compares well with other published analyses, capturing trends in global, hemispheric, and regional SST well, containing SST fields with more uniform variance through time and better month-to-month persistence than those in GISST. HadMAT1 is more consistent with SST and with collocated land surface air temperatures than previous NMAT data sets.

Implementation of Noah land surface model advances in the National Centers for Environmental Prediction operational mesoscale Eta model

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

Implementation of Noah land surface model advances in the National Centers for Environmental Prediction operational mesoscale Eta model : GEWEX Continental-Scale International Project, Part 3 (GCIP3)

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

An inventory of gaseous and primary aerosol emissions in Asia in the year 2000 : NASA global tropospheric experiment transport and chemical evolution over the pacific (TRACE-P): Measurements and analysis (TRACEP1)

Abstract: [i] An inventory of air pollutant emissions in Asia in the year 2000 is developed to support atmospheric modeling and analysis of observations taken during the TRACE-P experiment funded by the National Aeronautics and Space Administration (NASA) and the ACE-Asia experiment funded by the National Science Foundation (NSF) and the National Oceanic and Atmospheric Administration (NOAA). Emissions are estimated for all major anthropogenic sources, including biomass burning, in 64 regions of Asia. We estimate total Asian emissions as follows: 34.3 Tg SO 2 , 26.8 Tg NO x , 9870 Tg CO 2 , 279 Tg CO, 107 Tg CH 4 , 52.2 Tg NMVOC, 2.54 Tg black carbon (BC), 10.4 Tg organic carbon (OC), and 27.5 Tg NH 3 . In addition, NMVOC are speciated into 19 subcategories according to functional groups and reactivity. Thus we are able to identify the major source regions and types for many of the significant gaseous and particle emissions that influence pollutant concentrations in the vicinity of the TRACE-P and ACE-Asia field measurements. Emissions in China dominate the signature of pollutant concentrations in this region, so special emphasis has been placed on the development of emission estimates for China. China's emissions are determined to be as follows: 20.4 Tg SO 2 , 11.4 Tg NO x , 3820 Tg CO 2 , 116 Tg CO, 38.4 Tg CH 4 , 17.4 Tg NMVOC, 1.05 Tg BC, 3.4 Tg OC, and 13.6 Tg NH 3 . Emissions are gridded at a variety of spatial resolutions from 1° × 1° to 30 s x 30 s, using the exact locations of large point sources and surrogate GIS distributions of urban and rural population, road networks, landcover, ship lanes, etc. The gridded emission estimates have been used as inputs to atmospheric simulation models and have proven to be generally robust in comparison with field observations, though there is reason to think that emissions of CO and possibly BC may be underestimated. Monthly emission estimates for China are developed for each species to aid TRACE-P and ACE-Asia data interpretation. During the observation period of March/ April, emissions are roughly at their average values (one twelfth of annual). Uncertainties in the emission estimates, measured as 95% confidence intervals, range from a low of ±16% for SO 2 to a high of ±450% for OC.

An inventory of gaseous and primary aerosol emissions in Asia in the year 2000

Abstract: [1] An inventory of air pollutant emissions in Asia in the year 2000 is developed to support atmospheric modeling and analysis of observations taken during the TRACE-P experiment funded by the National Aeronautics and Space Administration (NASA) and the ACE-Asia experiment funded by the National Science Foundation (NSF) and the National Oceanic and Atmospheric Administration (NOAA). Emissions are estimated for all major anthropogenic sources, including biomass burning, in 64 regions of Asia. We estimate total Asian emissions as follows: 34.3 Tg SO2, 26.8 Tg NOx, 9870 Tg CO2, 279 Tg CO, 107 Tg CH4, 52.2 Tg NMVOC, 2.54 Tg black carbon (BC), 10.4 Tg organic carbon (OC), and 27.5 Tg NH3. In addition, NMVOC are speciated into 19 subcategories according to functional groups and reactivity. Thus we are able to identify the major source regions and types for many of the significant gaseous and particle emissions that influence pollutant concentrations in the vicinity of the TRACE-P and ACE-Asia field measurements. Emissions in China dominate the signature of pollutant concentrations in this region, so special emphasis has been placed on the development of emission estimates for China. China's emissions are determined to be as follows: 20.4 Tg SO2, 11.4 Tg NOx, 3820 Tg CO2, 116 Tg CO, 38.4 Tg CH4, 17.4 Tg NMVOC, 1.05 Tg BC, 3.4 Tg OC, and 13.6 Tg NH3. Emissions are gridded at a variety of spatial resolutions from 1° × 1° to 30 s × 30 s, using the exact locations of large point sources and surrogate GIS distributions of urban and rural population, road networks, landcover, ship lanes, etc. The gridded emission estimates have been used as inputs to atmospheric simulation models and have proven to be generally robust in comparison with field observations, though there is reason to think that emissions of CO and possibly BC may be underestimated. Monthly emission estimates for China are developed for each species to aid TRACE-P and ACE-Asia data interpretation. During the observation period of March/April, emissions are roughly at their average values (one twelfth of annual). Uncertainties in the emission estimates, measured as 95% confidence intervals, range from a low of ±16% for SO2 to a high of ±450% for OC.

Global frequency and distribution of lightning as observed from space by the Optical Transient Detector

Abstract: of uncertainty for the OTD global totals represents primarily the uncertainty (and variability) in the flash detection efficiency of the instrument The OTD measurements have been used to construct lightning climatology maps that demonstrate the geographical and seasonal distribution of lightning activity for the globe An analysis of this annual lightning distribution confirms that lightning occurs mainly over land areas, with an average land/ocean ratio of 10:1 The Congo basin, which stands out year-round, shows a peak mean annual flash density of 80 fl km 2 yr 1 in Rwanda, and includes an area of over 3 million km 2 exhibiting flash densities greater than 30 fl km 2 yr 1 (the flash density of central Florida) Lightning is predominant in the northern Atlantic and western Pacific Ocean basins year-round where instability is produced from cold air passing over warm ocean water Lightning is less frequent in the eastern tropical Pacific and Indian Ocean basins where the air mass is warmer A dominant Northern Hemisphere summer peak occurs in the annual cycle, and evidence is found for a tropically driven semiannual cycle INDEX TERMS: 3304 Meteorology and Atmospheric Dynamics: Atmospheric electricity; 3309 Meteorology and Atmospheric Dynamics: Climatology (1620); 3324 Meteorology and Atmospheric Dynamics: Lightning; 3394 Meteorology and Atmospheric Dynamics: Instruments and techniques;

Mineral Dust Entrainment and Deposition (DEAD) model: Description and 1990s dust climatology

Abstract: 17 ± 2 Tg; and optical depth at 0.63 mm, 0.030 ± 0.004. This emission, burden, and optical depth are significantly lower than some recent estimates. The model underestimates transport and deposition of East Asian and Australian dust to some regions of the Pacific Ocean. An underestimate of long-range transport of particles larger than 3 mm contributes to this bias. Our experiments support the hypothesis that dust emission ‘‘hot spots’’ exist in regions where alluvial sediments have accumulated and may be disturbed. INDEX TERMS: 0305 Atmospheric Composition and Structure: Aerosols and particles (0345, 4801); 0322 Atmospheric Composition and Structure: Constituent sources and sinks; 4801 Oceanography: Biological and Chemical: Aerosols (0305); 5415 Planetology: Solid Surface Planets: Erosion and weathering; KEYWORDS: mineral dust aerosol, aerosol climatology, mineral deposition, aerosol scavenging, saltation sandblasting, ecosystem fertilization

A global simulation of tropospheric ozone and related tracers: Description and evaluation of MOZART, version 2

Abstract: [1] We have developed a global three-dimensional chemical transport model called Model of Ozone and Related Chemical Tracers (MOZART), version 2. This model, which will be made available to the community, is built on the framework of the National Center for Atmospheric Research (NCAR) Model of Atmospheric Transport and Chemistry (MATCH) and can easily be driven with various meteorological inputs and model resolutions. In this work, we describe the standard configuration of the model, in which the model is driven by meteorological inputs every 3 hours from the middle atmosphere version of the NCAR Community Climate Model (MACCM3) and uses a 20-min time step and a horizontal resolution of 2.8° latitude × 2.8° longitude with 34 vertical levels extending up to approximately 40 km. The model includes a detailed chemistry scheme for tropospheric ozone, nitrogen oxides, and hydrocarbon chemistry, with 63 chemical species. Tracer advection is performed using a flux-form semi-Lagrangian scheme with a pressure fixer. Subgrid-scale convective and boundary layer parameterizations are included in the model. Surface emissions include sources from fossil fuel combustion, biofuel and biomass burning, biogenic and soil emissions, and oceanic emissions. Parameterizations of dry and wet deposition are included. Stratospheric concentrations of several long-lived species (including ozone) are constrained by relaxation toward climatological values. The distribution of tropospheric ozone is well simulated in the model, including seasonality and horizontal and vertical gradients. However, the model tends to overestimate ozone near the tropopause at high northern latitudes. Concentrations of nitrogen oxides (NOx) and nitric acid (HNO3) agree well with observed values, but peroxyacetylnitrate (PAN) is overestimated by the model in the upper troposphere at several locations. Carbon monoxide (CO) is simulated well at most locations, but the seasonal cycle is underestimated at some sites in the Northern Hemisphere. We find that in situ photochemical production and loss dominate the tropospheric ozone budget, over input from the stratosphere and dry deposition. Approximately 75% of the tropospheric production and loss of ozone occurs within the tropics, with large net production in the tropical upper troposphere. Tropospheric production and loss of ozone are three to four times greater in the northern extratropics than the southern extratropics. The global sources of CO consist of photochemical production (55%) and direct emissions (45%). The tropics dominate the chemistry of CO, accounting for about 75% of the tropospheric production and loss. The global budgets of tropospheric ozone and CO are generally consistent with the range found in recent studies. The lifetime of methane (9.5 years) and methylchloroform (5.7 years) versus oxidation by tropospheric hydroxyl radical (OH), two useful measures of the global abundance of OH, agree well with recent estimates. Concentrations of nonmethane hydrocarbons and oxygenated intermediates (carbonyls and peroxides) generally agree well with observations.

Variations in the light absorption coefficients of phytoplankton, nonalgal particles, and dissolved organic matter in coastal waters around Europe

Abstract: [1] We measured the absorption properties of phytoplankton, nonalgal particles (NAP), and colored dissolved organic matter (CDOM) at about 350 stations in various coastal waters around Europe including the English Channel, Adriatic Sea, Baltic Sea, Mediterranean Sea, and North Sea. For comparison, we also collected data in the open ocean waters of North Atlantic. The exponential slope of the CDOM absorption spectrum varied within a narrow range around 0.0176 nm−1 (SD = 0.0020 nm−1). When data from all the regions were considered altogether, the relationship between phytoplankton absorption and chlorophyll concentration was generally similar to the one previously established for open oceanic waters. Our coastal data, however, show that significant departures from the general trend may occur due to peculiar pigment composition and cell size. In some coastal areas, high phaeopigment concentrations gave rise to especially high blue-to-red ratio of phytoplankton absorption. The NAP absorption covaried with the particle dry weight. Most absorption spectra of these particles were well described by an exponential function with a slope averaging 0.0123 nm−1 (SD = 0.0013 nm−1). In some highly turbid waters, the spectra exhibited a signature possibly associated with iron oxides. In the Baltic Sea, NAP absorption systematically showed lower values at wavelengths shorter than 440 nm than predicted from the fitted exponential function. Overall, the variability in the absorption properties of European coastal waters showed some consistent patterns despite the high diversity of the examined waters. Distinct features were identified in the phytoplankton and NAP components. An absorption budget is presented and parameterizations are proposed.

An overview of ACE-Asia: Strategies for quantifying the relationships between Asian aerosols and their climatic impacts

Abstract: [1] The International Global Atmospheric Chemistry Program (IGAC) has conducted a series of Aerosol Characterization Experiments (ACE) that integrate in situ measurements, satellite observations, and models to reduce the uncertainty in calculations of the climate forcing due to aerosol particles. ACE-Asia, the fourth in this series of experiments, consisted of two focused components: (1) An intensive field study that sought to quantify the spatial and vertical distribution of aerosol concentrations and properties, the processes controlling their formation, evolution, and fate, and the column-integrated radiative effect of the aerosol (late March through May 2001). (2) A longer-term network of ground stations that used in situ and column-integrated measurements to quantify the chemical, physical, and optical properties of aerosols in the ACE-Asia study area and to assess their spatial and temporal (seasonal and interannual) variability (2000–2003). The approach of the ACE-Asia science team was to make simultaneous measurements of aerosol chemical, physical, and optical properties and their radiative impacts in a variety of air masses, often coordinated with satellite overpasses. Three aircraft, two research ships, a network of lidars, and many surface sites gathered data on Asian aerosols. Chemical transport models (CTMs) were integrated into the program from the start, being used in a forecast mode during the intensive observation period to identify promising areas for airborne and ship observations and then later as tools for integrating observations. The testing and improvement of a wide range of aerosol models (including microphysical, radiative transfer, CTM, and global climate models) was one important way in which we assessed our understanding of the properties and controlling processes of Asian aerosols. We describe here the scientific goals and objectives of the ACE-Asia experiment, its observational strategies, the types of observations made by the mobile platforms and stationary sites, the models that will integrate our understanding of the climatic effect of aerosol particles, and the types of data that have been generated. Eight scientific questions focus the discussion. The intensive observations took place during a season of unusually heavy dust, so we have a large suite of observations of dust and its interaction with air pollutants. Further information about ACE-Asia can be found on the project Web site at http://saga.pmel.noaa.gov/aceasia/.

Ambient aerosol sampling using the Aerodyne Aerosol Mass Spectrometer

Abstract: The Aerodyne Aerosol Mass Spectrometer (AMS) has been designed to measure size-resolved mass distributions and total mass loadings of volatile and semivolatile chemical species in/on submicron particles. This paper describes the application of this instrument to ambient aerosol sampling. The AMS uses an aerodynamic lens to focus the particles into a narrow beam, a roughened cartridge heater to vaporize them under high vacuum, and a quadrupole mass spectrometer to analyze the vaporized molecules. Particle size is measured via particle time-of-flight. The AMS is operated in two modes: (1) a continuous mass spectrum mode without size information; and (2) a size distribution measurement mode for selected m/z settings of the quadrupole. Single particles can also be detected and sized if they have enough mass of a chemical component. The AMS was deployed at a ground sampling site near downtown Atlanta during August 1999, as part of the Environmental Protection Agency/Southern Oxidant Study Particulate Matter “Supersite” experiment, and at a suburban location in the Boston area during September 1999. The major observed components of the aerosol at both sites were sulfate and organics with a minor fraction of nitrate, consistent with prior studies and colocated instruments. Different aerosol chemical components often had different size distributions and time evolutions. More than half of the sulfate mass was contained in 2% of the ambient particles in one of the sampling periods. Trends in mass concentrations of sulfate and nitrate measured with the AMS in Atlanta compare well with those measured with ion chromatography-based instruments. A marked diurnal cycle was observed for aerosol nitrate in Atlanta. A simple model fit is used to illustrate the integration of data from several chemical components measured by the AMS together with data from other particle instruments into a coherent representation of the ambient aerosol.

Subduction factory 1. Theoretical mineralogy, densities, seismic wave speeds, and H 2 O contents

Abstract: [1] We present a new compilation of physical properties of minerals relevant to subduction zones and new phase diagrams for mid-ocean ridge basalt, lherzolite, depleted lherzolite, harzburgite, and serpentinite. We use these data to calculate H2O content, density and seismic wave speeds of subduction zone rocks. These calculations provide a new basis for evaluating the subduction factory, including (1) the presence of hydrous phases and the distribution of H2O within a subduction zone; (2) the densification of the subducting slab and resultant effects on measured gravity and slab shape; and (3) the variations in seismic wave speeds resulting from thermal and metamorphic processes at depth. In considering specific examples, we find that for ocean basins worldwide the lower oceanic crust is partially hydrated (<1.3 wt % H2O), and the uppermost mantle ranges from unhydrated to � 20% serpentinized (� 2.4 wt % H2O). Anhydrous eclogite cannot be distinguished from harzburgite on the basis of wave speeds, but its � 6% greater density may render it detectable through gravity measurements. Subducted hydrous crust in cold slabs can persist to several gigapascals at seismic velocities that are several percent slower than the surrounding mantle. Seismic velocities and VP/VS ratios indicate that mantle wedges locally reach 60–80% hydration. INDEX TERMS: 3040 Marine Geology and Geophysics: Plate tectonics (8150, 8155, 8157, 8158); 3660 Mineralogy and Petrology: Metamorphic petrology; 3919 Mineral Physics: Equations of state; 5199 Physical Properties of Rocks: General or miscellaneous; 8123 Tectonophysics: Dynamics, seismotectonics; KEYWORDS: subduction, seismic velocities, mineral physics, H2O

Subduction factory 2. Are intermediate-depth earthquakes in subducting slabs linked to metamorphic dehydration reactions?

Abstract: [1] New thermal-petrologic models of subduction zones are used to test the hypothesis that intermediate-depth intraslab earthquakes are linked to metamorphic dehydration reactions in the subducting oceanic crust and mantle. We show that there is a correlation between the patterns of intermediate-depth seismicity and the locations of predicted hydrous minerals: Earthquakes occur in subducting slabs where dehydration is expected, and they are absent from parts of slabs predicted to be anhydrous. We propose that a subductingoceanicplatecanconsistoffourpetrologicallyandseismicallydistinctlayers:(1) hydrated, fine-grained basaltic upper crust dehydrating under equilibrium conditions and producing earthquakes facilitated by dehydration embrittlement; (2) coarse-grained, locally hydrated gabbroic lower crust that produces some earthquakes during dehydration but transformschieflyaseismicallytoeclogiteatdepthsbeyondequilibrium;(3)locallyhydrated uppermost mantle dehydrating under equilibrium conditions and producing earthquakes; and (4) anhydrous mantle lithosphere transforming sluggishly and aseismically to denser minerals. Fluid generated through dehydration reactions can move via at least three distinct flowpaths:percolationthroughlocal,transient,reaction-generatedhigh-permeabilityzones; flow through mode I cracks produced by the local stress state; and postseismic flow through fault zones. INDEX TERMS: 7218 Seismology: Lithosphere and upper mantle; 7230 Seismology: Seismicity and seismotectonics; 8123 Tectonophysics: Dynamics, seismotectonics; 8135 Tectonophysics: Evolution of the Earth: Hydrothermalsystems (8424); 3660 Mineralogyand Petrology: Metamorphicpetrology;

Models‐3 Community Multiscale Air Quality (CMAQ) model aerosol component 1. Model description

Abstract: [1] An initial evaluation of the Models-3 Community Multiscale Air Quality (CMAQ) model aerosol component reveals CMAQ's varying ability to simulate observed visibility indices and aerosol species concentrations. The visibility evaluation, using National Weather Service observations from 139 airports for 11–15 July 1995, shows that CMAQ reasonably captured the general spatial and temporal patterns of visibility degradation, including major gradients, maxima and minima. However, CMAQ's two visibility prediction methods, Mie theory approximation and mass reconstruction, both underpredict visibility degradation (i.e., overpredict visibility). The mean bias, normalized mean bias (NMB), mean error and normalized mean error (NME) for the Mie calculations are −5.9 dv, −21.7%, 7.0 dv and 25.4%, respectively. For the reconstruction simulations, these statistics are −9.8 dv, −35.5%, 10.0 dv and 36.2%, respectively. Most simulated values (∼90% Mie and ∼85% reconstruction) fall within a factor of two of the observations, although r2 = 0.25 (Mie) and r2 = 0.24 (reconstruction). The speciated aerosol evaluation uses observations of sulfate, nitrate, PM2.5, PM10 and organic carbon obtained from 18 stations of the Interagency Monitoring of Protected Visual Environments (IMPROVE) network in June 1995. This evaluation reveals that, with the exception of sulfate (mean bias: 0.15 μg/m3, NMB: 3.1%), the model consistently underpredicts aerosol concentrations of nitrate (−0.10 μg/m3, −33.1%), PM2.5 (−3.9 μg/m3, −30.1%), PM10 (−5.66 μg/m3, −29.2%) and organic carbon (−0.78 μg/m3, −33.7%). Sulfate was simulated best by the model (r2 = 0.63, mean error = 1.75 μg/m3, NME = 36.2%), followed by PM2.5 (0.55, 5.00 μg/m3, 38.5%), organic carbon (0.25, 0.94 μg/m3, 40.6%), PM10 (0.13, 9.85 μg/m3, 50.8%) and nitrate (0.01, 0.33 μg/m3, 104.3%). Except for nitrate, 75–80% of simulated concentrations fall within a factor of two of the IMPROVE observations.

A generalized single‐channel method for retrieving land surface temperature from remote sensing data

Abstract: [1] Many papers have developed algorithms to retrieve land surface temperature from at-sensor and land surface emissivity data. These algorithms have been specified for different thermal sensors on board satellites, i.e., the algorithm used for one thermal sensor (or a combination of thermal sensors) cannot be used for other thermal sensor. The main goal of this paper is to propose a generalized single-channel algorithm that only uses the total atmospheric water vapour content and the channel effective wavelength (assuming that emissivity is known), and can be applied to thermal sensors characterized with a FWHM (Full-Width Half-Maximum) of around 1 μm actually operative on board satellites. The main advantage of this algorithm compared with the other single-channel methods is that in-situ radiosoundings or effective mean atmospheric temperature values are not needed, whereas the main advantage of this algorithm compared with split-window and dual-angle methods is that it can be applied to different thermal sensors using the same equation and coefficients. The validation for different test sites shows root mean square deviations lower than 2 K for AVHRR channel 4 (λ ≈ 10.8 μm) and ATSR-2 channel 2 (λ ≈ 11 μm), and lower than 1.5 K for Landsat Thematic Mapper (TM) band 6 (λ ≈ 11.5 μm).

P wave tomography of the mantle under the Alpine-Mediterranean area

Abstract: [1] We study the upper mantle P wave velocity structure below the Euro-Mediterranean area, down to 1000 km depth, by seismic travel time tomography. We invert summary residuals constructed with both regional and teleseismic first arrival data reported by the International Seismological Centre (ISC) (1964–1995), introducing some alternative strategies in the travel time tomographic approach and a new scheme to correct teleseismic data for global mantle structure. Our high-resolution model PM0.5 is parameterized with three-dimensional (3-D) linear splines on a grid of nodes with 0.5° spacing in both horizontal directions and 50 km vertical spacing. We obtain about 26% root-mean-square (RMS) reduction of residuals by inversion in addition to roughly 31% reduction after summary rays formation and selection. Sensitivity analyses are performed through several test inversions to explore the resolution characteristics of the model at different spatial scales. The distribution of large-scale fast anomalies suggests that two different stages of a convection process presently coexist in very close regions. The mantle dynamics of western central Europe is dominated by blockage of subducted slabs at the 660 km discontinuity and ponding of seismically fast material in the transition zone. Contrarily, in the eastern Mediterranean, fast velocity material sinks into the lower mantle, suggesting that the flow of the cold downwelling here is not blocked by the 660 km discontinuity. On a smaller scale, the existence of tears in the subducted slab (lithospheric detachment) all along both margins of the Adriatic plate, as proposed by some authors, is not supported by our tomographic images.

A method for computing horizontal pressure-gradient force in an oceanic model with a nonaligned vertical coordinate

Abstract: [1] Discretization of the pressure-gradient force is a long-standing problem in terrain-following (or σ) coordinate oceanic modeling. When the isosurfaces of the vertical coordinate are not aligned with either geopotential surfaces or isopycnals, the horizontal pressure gradient consists of two large terms that tend to cancel; the associated pressure-gradient error stems from interference of the discretization errors of these terms. The situation is further complicated by the nonorthogonality of the coordinate system and by the common practice of using highly nonuniform stretching for the vertical grids, which, unless special precautions are taken, causes both a loss of discretization accuracy overall and an increase in interference of the component errors. In the present study, we design a pressure-gradient algorithm that achieves more accurate hydrostatic balance between the two components and does not lose as much accuracy with nonuniform vertical grids at relatively coarse resolution. This algorithm is based on the reconstruction of the density field and the physical z coordinate as continuous functions of transformed coordinates with subsequent analytical integration to compute the pressure-gradient force. This approach allows not only a formally higher order of accuracy, but it also retains and expands several important symmetries of the original second-order scheme to high orders [Mellor et al., 1994; Song, 1998], which is used as a prototype. It also has built-in monotonicity constraining algorithm that prevents appearance of spurious oscillations of polynomial interpolant and, consequently, insures numerical stability and robustness of the model under the conditions of nonsmooth density field and coarse grid resolution. We further incorporate an alternative method of dealing with compressibility of seawater, which escapes pressure-gradient errors associated with interference of the nonlinear nature of equation of state and difficulties to achieve accurate polynomial fits of resultant in situ density profiles. In doing so, we generalized the monotonicity constraint to guarantee nonnegative physical stratification of the reconstructed density profile in the case of compressible equation of state. To verify the new method, we perform traditional idealized (Seamount) and realistic test problems.

Interannual and seasonal variability of biomass burning emissions constrained by satellite observations

Abstract: We present a methodology for estimating the seasonal and interannual variation of biomass burning designed for use in global chemical transport models. The average seasonal variation is estimated from 4 years of fire-count data from the Along Track Scanning Radiometer (ATSR) and 1-2 years of similar data from the Advanced Very High Resolution Radiometer (AVHRR) World Fire Atlases. We use the Total Ozone Mapping Spectrometer (TOMS) Aerosol Index (AI) data product as a surrogate to estimate interannual variability in biomass burning for six regions: Southeast Asia, Indonesia and Malaysia, Brazil, Central America and Mexico, Canada and Alaska, and Asiatic Russia. The AI data set is available from 1979 to the present with an interruption in satellite observations from mid-1993 to mid-1996; this data gap is filled where possible with estimates of area burned from the literature for different regions. Between August 1996 and July 2000, the ATSR fire-counts are used to provide specific locations of emissions and a record of interannual variability throughout the world. We use our methodology to estimate mean seasonal and interannual variations for emissions of carbon monoxide from biomass burning, and we find that no trend is apparent in these emissions over the last two decades, but that there is significant interannual variability.

Relativistic electron pitch-angle scattering by electromagnetic ion cyclotron waves during geomagnetic storms

Abstract: During magnetic storms, relativistic electrons execute nearly circular orbits about the Earth and traverse a spatially confined zone within the duskside plasmapause where electromagnetic ion cyclotron (EMIC) waves are preferentially excited. We examine the mechanism of electron pitch-angle diffusion by gyroresonant interaction with EMIC waves as a cause of relativistic electron precipitation loss from the outer radiation belt. Detailed calculations are carried out of electron cyclotron resonant pitch-angle diffusion coefficients Dααfor EMIC waves in a multi-ion (H+, He+, O+) plasma. A simple functional form for Dαα is used, based on quasi-linear theory that is valid for parallel-propagating, small-amplitude electromagnetic waves of general spectral density. For typical observed EMIC wave amplitudes (l-10nT), the rates of resonant pitch-angle diffusion are close to the limit of "strong" diffusion, leading to intense electron precipitation. In order for gyroresonance to take place, electrons must possess a minimum kinetic energy Emin which depends on the value of the ratio (electron plasma frequency/ electron gyrofrequency); Emin also depends on the properties of the EMIC wave spectrum and the ion composition. Geophysically interesting scattering, with Emin comparable to 1 MeV, can only occur in regions where (electron plasma frequency/electron gyrofrequency) ≥ 10, which typically occurs within the duskside plasmapause. Under such conditions, electrons with energy ≥ 1 MeV can be removed from the outer radiation belt by EMIC wave scattering during a magnetic storm over a time-scale of several hours to a day.

Imaging spectroscopy: Earth and planetary remote sensing with the USGS Tetracorder and expert systems

Abstract: [1] Imaging spectroscopy is a tool that can be used to spectrally identify and spatially map materials based on their specific chemical bonds. Spectroscopic analysis requires significantly more sophistication than has been employed in conventional broadband remote sensing analysis. We describe a new system that is effective at material identification and mapping: a set of algorithms within an expert system decision-making framework that we call Tetracorder. The expertise in the system has been derived from scientific knowledge of spectral identification. The expert system rules are implemented in a decision tree where multiple algorithms are applied to spectral analysis, additional expert rules and algorithms can be applied based on initial results, and more decisions are made until spectral analysis is complete. Because certain spectral features are indicative of specific chemical bonds in materials, the system can accurately identify and map those materials. In this paper we describe the framework of the decision making process used for spectral identification, describe specific spectral feature analysis algorithms, and give examples of what analyses and types of maps are possible with imaging spectroscopy data. We also present the expert system rules that describe which diagnostic spectral features are used in the decision making process for a set of spectra of minerals and other common materials. We demonstrate the applications of Tetracorder to identify and map surface minerals, to detect sources of acid rock drainage, and to map vegetation species, ice, melting snow, water, and water pollution, all with one set of expert system rules. Mineral mapping can aid in geologic mapping and fault detection and can provide a better understanding of weathering, mineralization, hydrothermal alteration, and other geologic processes. Environmental site assessment, such as mapping source areas of acid mine drainage, has resulted in the acceleration of site cleanup, saving millions of dollars and years in cleanup time. Imaging spectroscopy data and Tetracorder analysis can be used to study both terrestrial and planetary science problems. Imaging spectroscopy can be used to probe planetary systems, including their atmospheres, oceans, and land surfaces.

Intercomparison of remote sounding instruments

Abstract: [1] When intercomparing measurements made by remote sounders, it is necessary to make due allowance for the differing characteristics of the observing systems, particularly their averaging kernels and error covariances. We develop the methods required to do this, applicable to any kind of retrieval method, not only to optimal estimators. We show how profiles and derived quantities such as the total column of a constituent may be properly compared, yielding different averaging kernels. We find that the effect of different averaging kernels can be reduced if the retrieval or the derived quantity of one instrument is simulated using the retrieval of the other. We also show how combinations of measured signals can be found, which can be compared directly. To illustrate these methods, we apply them to two real instruments, calculating the expected amplitudes and variabilities of the diagnostics for a comparison of CO measurements made by a ground-based Fourier Transform spectrometer (FTIR) and the “measurement of pollution in the troposphere” instrument (MOPITT), which is mounted on the EOS Terra platform. The main conclusions for this case are the following: (1) Direct comparison of retrieved profiles is not satisfactory, because the expected standard deviation of the difference is around half of the expected natural variability of the true atmospheric profiles. (2) Comparison of the MOPITT profile retrieval with a simulation using FTIR is much more useful, though still not ideal, with expected standard deviation of differences of around 20% of the expected natural variability. (3) Direct comparison of total columns gives an expected standard deviation of about 9%, while comparison of MOPITT with a simulation derived from FTIR improved this to 8%. (4) There is only one combination of measured signals that can be usefully compared. The difference is expected to have a standard deviation of about 5.5% of the expected natural variability, which is mostly due to noise.

Emission from international sea transportation and environmental impact

Abstract: [1] Emission generated by the international merchant fleet has been suggested to represent a significant contribution to the global anthropogenic emissions. To analyze the impacts of these emissions, we present detailed model studies of the changes in atmospheric composition of pollutants and greenhouse compounds due to emissions from cargo and passenger ships in international trade. Global emission inventories of NOx, SO2, CO, CO2, and volatile organic compounds (VOC) are developed by a bottom-up approach combining ship-type specific engine emission modeling, oil cargo VOC vapor modeling, alternative global distribution methods, and ship operation data. Calculated bunker fuel consumption is found in agreement with international sales statistics. The Automated Mutual-assistance Vessel Rescue system (AMVER) data set is found to best reflect the distributions of cargo ships in international trade. A method based on the relative reporting frequency weighted by the ship size for each vessel type is recommended. We have exploited this modeled ship emissions inventory to estimate perturbations of the global distribution of ozone, methane, sulfate, and nitrogen compounds using a global 3-D chemical transport model with interactive ozone and sulfate chemistry. Ozone perturbations are highly nonlinear, being most efficient in regions of low background pollution. Different data sets (e.g., AMVER, The Comprehensive Ocean-Atmosphere Data Set (COADS)) lead to highly different regional perturbations. A maximum ozone perturbation of approximately 12 ppbv is obtained in the North Atlantic and in the North Pacific during summer months. Global average sulfate loading increases with 2.9%, while the increase is significantly larger over parts of western Europe (up to 8%). In contrast to the AMVER data, the COADS data give particularly large enhancements over the North Atlantic. Ship emissions reduce methane lifetime by approximately 5%. CO2 and O3 give positive radiative forcing (RF), and CH4 and sulfate give negative forcing. The total RF is small (0.01–0.02 W/m2) and connected with large uncertainties. Increase in acidification is 3–10% in certain coastal areas. The approach presented here is clearly useful for characterizing the present impact of ship emission and will be valuable for assessing the potential effect of various emission-control options.

Spectral discrimination of coarse and fine mode optical depth

Abstract: [1] The recognition that the aerosol particle size distribution (PSD) is effectively bimodal permits the extraction of the fine and coarse mode optical depths (tf and tc) from the spectral shape of the total aerosol optical depth (ta = tf + tc). This purely optical technique avoids intermediate computations of the PSD and yields a direct optical output that is commensurate in complexity with the spectral information content of ta. The separation into tf and tc is a robust process and yields aerosol optical statistics, which are more intrinsic than those, obtained from a generic analysis of ta. Partial (optical) validation is provided by (1) demonstrating the physical coherence of the simple model employed, (2) demonstrating that tc variation is coherent with photographic evidence of thin cloud events and that tf variation is coherent with photographic evidence of clear sky and haze events, and (3) showing that the retrieved values of tf and tc are wellcorrelated, if weakly biased, relative to formal inversions of combined solar extinction and sky radiance data. The spectral inversion technique permitted a closer scrutiny of a standard (temporally based) cloud-screening algorithm. Perturbations of monthly or longer-term statistics associated with passive or active shortcomings of operational cloud screening were inferred to be small to occasionally moderate over a sampling of cases. Diurnal illustrations were given where it was clear that such shortcomings can have a significant impact on the interpretation of specific events; (1) commission errors in tf due to the exclusion of excessively high-frequency fine mode events and (2) omission errors in tc due to the inclusion of insufficiently high-frequency thin homogeneous cloud events. INDEX TERMS: 0305 Atmospheric Composition and Structure: Aerosols and particles (0345, 4801); 0360 Atmospheric Composition and Structure: Transmission and scattering of radiation; 3360 Meteorology and Atmospheric Dynamics: Remote sensing; 4801 Oceanography: Biological and Chemical: Aerosols (0305); KEYWORDS: Sun photometry, aerosol, optical depth, fine mode, coarse mode, cloud screening

Updated emissions from ocean shipping

Abstract: [1] Marine vessel inventories demonstrate that ship emissions cannot be neglected in assessing environmental impacts of air pollution, although significant uncertainty in these inventories remains. We address this uncertainty by employing a bottom-up estimate of fuel consumption and vessel activity for internationally registered fleets, including cargo vessels, other commercial vessels, and military vessels. We identify model bias in previous work, which assumed internationally registered ships primarily consume international marine fuels. Updated results suggest fuel consumption is ∼289 million metric tons per year, more than twice the quantity reported as international fuel. According to our analysis, fuel used by internationally registered fleets is apparently allocated to both international and domestic fuel statistics; this implies either that ships operate along domestic routes much of the time or that marine fuel sales to these ships may be misassigned. If the former is true, then allocation of emissions to international shipping routes may underestimate near-coastal emissions from ships. Our updated inventories increases previous ship emissions inventories for all pollutants; for example, global NOx emissions (∼6.87 Tg N) are more than doubled. This work also produces detailed sensitivity analyses of inputs to these estimates, identifying uncertainty in vessel duty-cycle as critical to overall emissions estimates. We discuss implications for assessing ship emissions impacts.

Real‐time and retrospective forcing in the North American Land Data Assimilation System (NLDAS) project

Abstract: [1] The accuracy of forcing data greatly impacts the ability of land surface models (LSMs) to produce realistic simulations of land surface processes. With this in mind, the multi-institutional North American Land Data Assimilation System (NLDAS) project has produced retrospective (1996–2002) and real-time (1999–present) data sets to support its LSM modeling activities. Featuring 0.125° spatial resolution, hourly temporal resolution, nine primary forcing fields, and six secondary validation/model development fields, each data set is based on a backbone of Eta Data Assimilation System/Eta data and is supplemented with observation-based precipitation and radiation data. Hourly observation-based precipitation data are derived from a combination of daily National Center for Environmental Prediction Climate Prediction Center (CPC) gauge-based precipitation analyses and hourly National Weather Service Doppler radar-based (WSR-88D) precipitation analyses, wherein the hourly radar-based analyses are used to temporally disaggregate the daily CPC analyses. NLDAS observation-based shortwave values are derived from Geostationary Operational Environmental Satellite radiation data processed at the University of Maryland and at the National Environmental Satellite Data and Information Service. Extensive quality control and validation efforts have been conducted on the NLDAS forcing data sets, and favorable comparisons have taken place with Oklahoma Mesonet, Atmospheric Radiation Measurement Program/cloud and radiation test bed, and Surface Radiation observation data. The real-time forcing data set is constantly evolving to make use of the latest advances in forcing-related data sets, and all of the real-time and retrospective data are available online at http://ldas.gsfc.nasa.gov for visualization and downloading in both full and subset forms.

Quantitative sampling using an Aerodyne aerosol mass spectrometer 1. Techniques of data interpretation and error analysis

Abstract: Received 22 March 2002; revised 2 July 2002; accepted 5 August 2002; published 4 February 2003. [1] The aerosol mass spectrometer (AMS), manufactured by Aerodyne Research, Inc., has been shown to be capable of delivering quantitative information on the chemical composition and size of volatile and semivolatile fine airborne particulate matter with high time resolution. Analytical and software tools for interpreting the data from this instrument and generating meaningful, quantitative results have been developed and are presented here with a brief description of the instrument. These include the conversion of detected ion rates from the quadrupole mass spectrometer during the mass spectrum (MS) mode of operation to atmospheric mass concentrations of chemical species (in m gm � 3 ) by applying calibration data. It is also necessary to correct for variations in the electron multiplier performance, and a method involving the measurement of the instrument’s response to gas phase signals is also presented. The techniques for applying particle velocity calibration data and transforming signals from time of flight (TOF) mode to chemical mass distributions in terms of aerodynamic diameter (dM/dlog(Da) distributions) are also presented. It is also possible to quantify the uncertainties in both MS and TOF data by evaluating the ion counting statistics and variability of the background signal, respectively. This paper is accompanied by part 2 of this series, in which these methods are used to process and analyze AMS results on ambient aerosol from two U.K. cities at different times of the year. INDEX TERMS: 0305 Atmospheric Composition and Structure: Aerosols and particles (0345, 4801); 0394 Atmospheric Composition and Structure: Instruments and techniques; 0399 Atmospheric Composition and Structure: General or miscellaneous; KEYWORDS: aerosols, chemical composition, mass spectrometry, analysis techniques

Interplanetary coronal mass ejections in the near‐Earth solar wind during 1996–2002

Abstract: We summarize the occurrence of interplanetary coronal mass injections (ICMEs) in the near-Earth solar wind during 1996-2002, corresponding to the increasing and maximum phases of solar cycle 23. In particular, we give a detailed list of such events. This list, based on in-situ observations, is not confined to subsets of ICMEs, such as magnetic clouds or those preceded by halo CMEs observed by the SOHO/LASCO coronagraph, and provides an overview of 214 ICMEs in the near-Earth solar wind during this period. The ICME rate increases by about an order of magnitude from solar minimum to solar maximum (when the rate is approximately 3 ICMEs/solar rotation period). The rate also shows a temporary reduction during 1999, and another brief, deeper reduction in late 2000-early 2001, which only approximately track variations in the solar 10 cm flux. In addition, there are occasional periods of several rotations duration when the ICME rate is enhanced in association with high solar activity levels. We find an indication of a periodic variation in the ICME rate, with a prominent period of approximately 165 days similar to that previously reported in various solar phenomena. It is found that the fraction of ICMEs that are magnetic clouds has a solar cycle variation, the fraction being larger near solar minimum. For the subset of events that we could associate with a CME at the Sun, the transit speeds from the Sun to the Earth were highest after solar maximum.

Transport and Chemical Evolution over the Pacific (TRACE-P) aircraft mission: Design, execution, and first results

Abstract: The NASA Transport and Chemical Evolution over the Pacific (TRACE-P) aircraft mission was conducted in February-April 2001 over the NW Pacific (1) to characterize the Asian chemical outflow and relate it quantitatively to its sources and (2) to determine its chemical evolution. It used two aircraft, a DC-8 and a P-3B, operating out of Hong Kong and Yokota Air Force Base (near Tokyo), with secondary sites in Hawaii, Wake Island, Guam, Okinawa, and Midway. The aircraft carried instrumentation for measurements of long-lived greenhouse gases, ozone and its precursors, aerosols and their precursors, related species, and chemical tracers. Five chemical transport models (CTMs) were used for chemical forecasting. Customized bottom-up emission inventories for East Asia were generated prior to the mission to support chemical forecasting and to serve as a priori for evaluation with the aircraft data. Validation flights were conducted for the Measurements Of Pollution In The Troposphere (MOPITT) satellite instrument and revealed little bias (6 plus or minus 2%) in the MOPITT measurements of CO columns. A major event of transpacific Asian pollution was characterized through combined analysis of TRACE-P and MOPITT data. The TRACE-P observations showed that cold fronts sweeping across East Asia and the associated warm conveyor belts (WCBs) are the dominant pathway for Asian outflow to the Pacific in spring. The WCBs lift both anthropogenic and biomass burning (SE Asia) effluents to the free troposphere, resulting in complex chemical signatures. The TRACE-P data are in general consistent with a priori emission inventories, lending confidence in our ability to quantify Asian emissions from socioeconomic data and emission factors. However, the residential combustion source in rural China was found to be much larger than the a priori, and there were also unexplained chemical enhancements (HCN, CH3Cl, OCS, alkylnitrates) in Chinese urban plumes. The Asian source of CCl4 was found to be much higher than government estimates. Measurements of HCN and CH3CN indicated a dominant biomass burning source and ocean sink for both gases. Large fractions of sulfate and nitrate were found to be present in dust aerosols. Photochemical activity in the Asian outflow was strongly reduced by aerosol attenuation of UV radiation, with major implications for the concentrations of HOx, radicals. New particle formation, apparently from ternary nucleation involving NH3, was observed in Chinese urban plumes.

Global monitoring of air pollution over land from the Earth Observing System‐Terra Moderate Resolution Imaging Spectroradiometer (MODIS)

Abstract: [1] Moderate Resolution Imaging Spectroradiometer (MODIS) measurements (7 channels: 0.47-2.1 μm, 250-500 m resolutions) provide us with new insights into the characteristics of global aerosols. MODIS retrieves not only aerosol loading but also the fraction of fine mode particle. In this paper we demonstrate MODIS capability for use in monitoring global, regional, and local air pollution. Three case studies in northern Italy, Los Angeles, and Beijing showed the conclusive results of applying MODIS-derived aerosol optical depths (T a ) to regional and local air pollution in terms of accuracy (ΔT a = ±0.05 ± 0.2τ a ) and spatial sensitivity of the retrievals. Under stagnant condition, accumulated aerosol abundance can reach T a > 1 (at 0.55 μm) before being removed by wind or precipitation. The correlation found between Aerosol Robotic Network (AERONET) daily averaged T a and 24-hour PM 10 (particulate matter with diameter <10 μm) concentration (μg/m 3 ) in northern Italy is encouraging with correlation coefficient ∼0.82. The derivation of PM concentration from satellite measurements may be possible once we know the detailed aerosol vertical distribution. To compare aerosol loading in different regions of the globe, we choose the two most populated regions (eastern China and India) and the two most industrialized regions (the eastern United States/Canada and western Europe). The time series of MODIS monthly mean T a from July 2000 to May 2001 depicts a strong seasonal variation with maxima in the spring/ summer and minima in the winter. The clear separation between (1) the eastern United States/Canada and western Europe and (2) eastern China and India shows that the T a values in (2) are 50% to 2-3 times higher compared to those in (1). The enhancements of aerosol loading were due to smoke as originated from Montana/Idaho forest fires transported to the eastern United States in late August 2000 and dust outbreaks from Taklimakan and Gobi Deserts to eastern China as well as smoke from Southeast Asia to southern China in February-April 2001.

Global inventory of nitrogen oxide emissions constrained by space‐based observations of NO2 columns

Abstract: [1] We use tropospheric NO2 columns from the Global Ozone Monitoring Experiment (GOME) satellite instrument to derive top-down constraints on emissions of nitrogen oxides (NOx ≡ NO + NO2), and combine these with a priori information from a bottom-up emission inventory (with error weighting) to achieve an optimized a posteriori estimate of the global distribution of surface NOx emissions. Our GOME NO2 retrieval improves on previous work by accounting for scattering and absorption of radiation by aerosols; the effect on the air mass factor (AMF) ranges from +10 to −40% depending on the region. Our AMF also includes local information on relative vertical profiles (shape factors) of NO2 from a global 3-D chemical transport model (GEOS-CHEM); assumption of a globally uniform shape factor, as in most previous retrievals, would introduce regional biases of up to 40% over industrial regions and a factor of 2 over remote regions. We derive a top-down NOx emission inventory from the GOME data by using the local GEOS-CHEM relationship between NO2 columns and NOx emissions. The resulting NOx emissions for industrial regions are aseasonal, despite large seasonal variation in NO2 columns, providing confidence in the method. Top-down errors in monthly NOx emissions are comparable with bottom-up errors over source regions. Annual global a posteriori errors are half of a priori errors. Our global a posteriori estimate for annual land surface NOx emissions (37.7 Tg N yr−1) agrees closely with the GEIA-based a priori (36.4) and with the EDGAR 3.0 bottom-up inventory (36.6), but there are significant regional differences. A posteriori NOx emissions are higher by 50–100% in the Po Valley, Tehran, and Riyadh urban areas, and by 25–35% in Japan and South Africa. Biomass burning emissions from India, central Africa, and Brazil are lower by up to 50%; soil NOx emissions are appreciably higher in the western United States, the Sahel, and southern Europe.