A study of the total atmospheric sulfur dioxide load using ground-based measurements and the satellite derived Sulfur Dioxide Index

We present characteristics of the sulfur dioxide (SO₂) loading over Thessaloniki, Greece, and seven other selected sites around the world using SO₂ total column measurements from Brewer spectrophotometers together with satellite estimates of the Version 8 TOMS Sulfur Dioxide Index (SOI) over the same locations, retrieved from Nimbus 7 TOMS (1979–1993), Earth Probe TOMS (1996–2003) and OMI/Aura (2004–2006). Traditionally, the SOI has been used to quantify the SO₂ quantities emitted during great volcanic eruptions. Here, we investigate whether the SOI can give an indication of the total SO₂ load for areas and periods away from eruptive volcanic activity by studying its relative changes as a correlative measure to the SO₂ total column. We examined time series from Thessaloniki and another seven urban and non-urban stations, five in the European Union (Arosa, De Bilt, Hohenpeissenberg, Madrid, Rome) and two in India (Kodaikanal, New Delhi). Based on the Brewer data, Thessaloniki shows high SO₂ total columns for a European Union city but values are still low if compared to highly affected regions like those in India. For the time period 1983–2006 the SO₂ levels above Thessaloniki have generally decreased with a rate of 0.028 Dobson Units (DU) per annum, presumably due to the European Union's strict sulfur control policies. The seasonal variability of the SO₂ total column exhibits a double peak structure with two maxima, one during winter and the second during summer. The winter peak can be attributed to central heating while the summer peak is due to synoptic transport from sources west of the city and sources in the north of Greece. A moderate correlation was found between the seasonal levels of Brewer total SO₂ and SOI for Thessaloniki, Greece (R = 0.710–0.763) and Madrid, Spain (R = 0.691) which shows that under specific conditions the SOI might act as an indicator of the SO₂ total load.     

Adjoint inverse modeling of NOx emissions over eastern China using satellite observations of NO2 vertical column densities

A four-dimensional variational data assimilation system for optimization of NO_x emissions (RC4-NO_x) was developed. A parameterized NO_x chemistry scheme was introduced into the RC4-NO_x system, and key parameters such as chemical production and loss terms of NO_x were calculated in advance using the Community Multiscale Air Quality (CMAQ) modeling system. RC4-NO_x was applied to optimize NO_x emissions over eastern China (EC) in July 1996, 1999, and 2002 using Global Ozone Monitoring Experiment (GOME) satellite observations of NO₂ vertical column densities (VCDs) and a priori emissions from the Regional Emission Inventory in Asia (REAS). After assimilation, RC4-NO_x generally reproduced the spatial distribution, regional averaged values, and time evolution of GOME NO₂ VCDs. Over EC, a priori emissions were reduced by 20% in 1996 and by 8% in 1999, whereas a posteriori emissions were almost the same as a priori emissions in 2002. A priori emissions in the Beijing region were reduced by optimization over the whole simulation period. A posteriori emissions over the Yangtze Delta were larger than a priori emissions in 2002, although they were smaller in both 1996 and 1999. As in other areas, a priori emissions over the North China Plain were reduced in 1996; but those over the eastern part of the plain were increased in 1999, and the area of increased emissions moved slightly westward in 2002. In each region, the growth rates of a posteriori emissions during both 1996–1999 and 1999–2002 became generally larger than those of a priori emissions, and the trends of a posteriori emissions became similar to those of GOME NO₂ VCDs. Our inverse modeling analysis indicates that the rate of increase of NO_x emissions over EC from 1996 to 2002 was much larger for a posteriori emissions (49%) than for a priori emissions (19%).     

Spatiotemporal modelling of ozone distribution in the State of California

This paper is concerned with the spatiotemporal mapping of monthly 8-h average ozone (O₃) concentrations over California during a 15-years period. The basic methodology of our analysis is based on the spatiotemporal random field (S/TRF) theory. We use a S/TRF decomposition model with a dominant seasonal O₃ component that may change significantly from site to site. O₃ seasonal patterns are estimated and separated from stochastic fluctuations. By means of Bayesian Maximum Entropy (BME) analysis, physically meaningful and sufficiently detailed space–time maps of the seasonal O₃ patterns are generated across space and time. During the summer and winter months the seasonal O₃ concentration maps exhibit clear and progressively changing geographical patterns over time, suggesting the existence of relationships in accordance with the typical physiographic and climatologic features of California. BME mapping accuracy can be superior to that of other techniques commonly used by EPA; its framework can rigorously assimilate useful data sources that were previously unaccounted for; the generated maps offer valuable assessments of the spatiotemporal O₃ patterns that can be helpful in the identification of physical mechanisms and their interrelations, the design of human exposure and population health models, and in risk assessment. As they focus on the seasonal patterns, the maps are not contingent on short-time and locally prevalent weather conditions, which are of no interest in a global and non-forecasting framework. Moreover, the maps offer valuable insight about the space–time O₃ concentration patterns and are, thus, helpful for disentangling the influence of explanatory factors or even for identifying some influential ones that could have been otherwise overlooked.     

Atmospheric chemistry of perfluorobutenes (CF3CFCFCF3 and CF3CF2CFCF2): Kinetics and mechanisms of reactions with OH radicals and chlorine atoms,IR spectra,global warming potentials,and oxidation to perfluorocarboxylic acids

Relative rate techniques were used to determine k(Cl + CF₃CFCFCF₃) = (7.27 ± 0.88) × 10^?12, k(Cl + CF₃CF₂CFCF₂) = (1.79 ± 0.41) × 10^?11, k(OH + CF₃CFCFCF₃) = (4.82 ± 1.15) × 10^?13, and k(OH + CF₃CF₂CFCF₂) = (1.94 ± 0.27) × 10^?12 cm³ molecule^?1 s^?1 in 700 Torr of air or N₂ diluent at 296 K. The chlorine atom- and OH radical-initiated oxidation of CF₃CFCFCF₃ in 700 Torr of air gives CF₃C(O)F in molar yields of 196 ± 11 and 218 ± 20%, respectively. Chlorine atom-initiated oxidation of CF₃CF₂CFCF₂ gives molar yields of 97 ± 9% CF₃CF₂C(O)F and 97 ± 9% COF₂. OH radical-initiated oxidation of CF₃CF₂CFCF₂ gives molar yields of 110 ± 15% CF₃CF₂C(O)F and 99 ± 8% COF₂. The atmospheric fate of CF₃CF₂C(O)F and CF₃C(O)F is hydrolysis to give CF₃CF₂C(O)OH and CF₃C(O)OH. The atmospheric lifetimes of CF₃CFCFCF₃ and CF₃CF₂CFCF₂ are determined by reaction with OH radicals and are approximately 24 and 6 days, respectively. The contribution of CF₃CFCFCF₃ and CF₃CF₂CFCF₂ to radiative forcing of climate change will be negligible.     

Ammonia emissions from a representative in-use fleet of light and medium-duty vehicles in the California South Coast Air Basin

We used Fourier Transform Infrared Spectroscopy (FTIR) to measure tailpipe ammonia emissions from a representative fleet of 41 light and medium-duty vehicles recruited in the California South Coast Air Basin. A total of 121 chassis dynamometer emissions tests were conducted on these vehicles and the test results were examined to determine the effects of several key variables on ammonia emissions. Variables included vehicle type, driving cycle, emissions technology, ammonia precursor emissions (i.e. CO and NOx) and odometer readings/model year as a proxy for catalyst age. The mean ammonia emissions factor was 46 mg km^?1 (σ = 48 mg km^?1) for the vehicle fleet. Average emission factors for specific vehicle groups are also reported in this study. Results of this study suggest vehicles with the highest ammonia emission rates possess the following characteristics: medium-duty vehicles, older emissions technologies, mid-range odometer readings, and higher CO emissions. In addition, vehicles subjected to aggressive driving conditions are likely to be higher ammonia emitters. Since the vehicles we studied were representative of recent model year vehicles and technologies in urban airsheds, the results of our study will be useful for developing ammonia emissions inventories in Los Angeles and other urban areas where California-certified vehicles are driven. However, efforts should also be made to continue emissions testing on in-use vehicles to ensure greater confidence in the ammonia emission factors reported here.     

Effects of post-sampling conditions on ambient carbon aerosol filter measurements

Ambient carbonaceous material collected on quartz filters is prone to measurement artifacts due to material gained or lost during post-sampling field latency, shipping, and storage. In seventeen sampling events over a one year period, ambient PM_2.5 aerosols were collected on quartz filters (without denuders) and subjected to various filter treatments to assess the potential for and extent of artifacts. The filter treatments simulated post-sampling environments that filters may be exposed to and included: storage at 40 °C for up to 96 h, storage at ?16 °C for 48 h, and storage at room temperature (～21 °C) for 48 h. Carbon mass on the filters was measured using a thermal-optical method. The total carbon (TC), total organic carbon (TOC) and total elemental carbon (TEC) as well as carbon thermal fraction masses were obtained. Statistical analyses were performed to identify significant differences in carbon fraction concentrations between filters analyzed immediately after sampling and after being subjected to treatment.TOC and TC concentrations decreased by on average 15 ± 5% and 10 ± 4%, respectively, for filters maintained at 40 °C for 96 h but did not change for filters stored at room temperature or frozen for 48 h. TEC did not change for any of the filter treatments. The mass concentration for the organic carbon thermal fraction that evolves at the lowest temperature step (OC1) decreased with increasing storage time at 40 °C with average losses of 70 ± 7% after 96 h. Therefore, OC1 is not a stable measurement due to post-sampling conditions that may be encountered. This work demonstrates that TOC and TC can have substantial measurement artifacts on filters subjected to field latency and other non-temperature controlled post-sampling handling, compared to the carbon loadings on the filter at the end of the sampling period.     

The role of satellite and decentralized strategies in water resources management

Existing and projected water shortages and related factors have helped focus attention on the need for water reuse. With recent technological advances in wastewater treatment, it is now possible to produce reclaimed water of any quality. Thus, the use of reclaimed water will depend on the reuse opportunities and the cost of the required infrastructure. Historically, centralized wastewater treatment facilities have served the needs of organized societies since the mid 1800s. However, as there are limited options for expansion of most existing centralized facilities, the use of satellite and decentralized wastewater management systems offers significant advantages including being close both to the source of wastewater generation and to potential water reuse applications. The comparative advantages of satellite and decentralized wastewater management systems for a number of water reuse applications are presented and discussed in this paper. Selected case studies are presented to demonstrate the utility of satellite and decentralized wastewater management. Specific issues associated with the application of such systems in existing and in new developments are examined and discussed.     

Mass conservative,positive definite integrator for atmospheric chemical dynamics

Air quality models compute the transformation of species in the atmosphere undergoing chemical and physical changes. The numerical algorithms used to predict these transformations should obey mass conservation and positive definiteness properties. Among all physical phenomena, the chemical kinetics solver provides the greatest challenge to attain these two properties. In general, most chemical kinetics solvers are mass conservative but not positive definite. In this article, a new numerical algorithm for the computation of chemical kinetics is presented. The integrator is called Split Single Reaction Integrator (SSRI). It is both mass conservative and positive definite. It solves each chemical reaction exactly and uses operator splitting techniques (symmetric split) to combine them into the entire system.The method can be used within a host integrator to fix the negative concentrations while preserving the mass, or it can be used as a standalone integrator that guarantees positive definiteness and mass conservation. Numerical results show that the new integrator, used as a standalone integrator, is second order accurate and stable under large fixed time steps when other conventional integrators are unstable.     

Modeling Postfire Response and Recovery using the Hydrologic Engineering Center Hydrologic Modeling System (HEC‐HMS)1

Abstract: This paper investigates application of the Army Corps of Engineers’ Hydrologic Engineering Center Hydrologic Modeling System (HEC‐HMS) to a burned watershed in San Bernardino County, California. We evaluate the HEC‐HMS’ ability to simulate discharge in prefire and postfire conditions in a semi arid watershed and the necessary parameterizations for modeling hydrologic response during the immediate, and subsequent recovery, period after a wildfire. The model is applied to City Creek watershed, which was 90% burned during the Old Fire of October 2003. An optimal spatial resolution for the HEC‐HMS model was chosen based on an initial sensitivity analysis of subbasin configurations and related model performance. Five prefire storms were calibrated for the selected model resolution, defining a set of parameters that reasonably simulate prefire conditions. Six postfire storms, two from each of the following rainy (winter) seasons were then selected to simulate postfire response and evaluate relative changes in parameter values and model behavior. There were clear trends in the postfire parameters [initial abstractions (I_a), curve number (CN), and lag time] that reveal significant (and expected) changes in watershed behavior. CN returns to prefire (baseline) values by the end of Year 2, while I_a approaches baseline by the end of the third rainy season. However, lag time remains significantly lower than prefire values throughout the three‐year study period. Our results indicate that recovery of soil conditions and related runoff response is not entirely evidenced by the end of the study period (three rainy seasons postfire). Understanding the evolution of the land surface and related hydrologic properties during the highly dynamic postfire period, and accounting for these changes in model parameterizations, will allow for more accurate and reliable discharge simulations in both the immediate, and subsequent, rainy seasons following fire.     

Mathematical analysis of the whole core injection method accuracy for measuring phenanthrene biodegradation rates in undisturbed marine sediments

Rates of 14C-phenanthrene mineralization in contaminated, undisturbed marine sediments were measured using the whole core injection method to assess microbial natural attenuation activity as a function of sediment depth. Submerged sediments were sampled from Eagle Harbor, a marine superfund site in Puget Sound. Experiments show significant biodegradation activities (0.0012-0.0036 day(-1)) in the sediment horizons from 0 to 10 cm. The purpose and scope of this paper is to evaluate the range of experimental conditions giving valid results; a mathematical simulation described competing contaminant 14C-phenanthrene diffusion and simultaneous biodegradation (Monod kinetics), both retarded by sorption. The effect of aging was examined with two sorption models in presumed pseudo-homogenous sediments having effective properties. The simulation predictions provide quantitative guidelines for the successful use of the whole core injection method. (1) The effective Monod constant KS' in sediment is increased by a large partition coefficient KP between sediment and water and makes the apparent 14C-phenanthrene biodegradation approach first-order kinetics. (2) When KS'>1 mg(-1) l(-1), the measured 14C-phenanthrene biodegradation extent is biased by inadequately distributed injected tracer only when less than 7% of the sediment horizon is initially probed and mixed with injected tracer. (3) A short incubation time (<20 days) is necessary when a mobile indicator, e.g., gaseous 14CO2, is used. For longer incubation times, predictions show that a 14CO2 indicator diffuses to adjacent horizons, thus smearing the depth profile of biodegradation. (4) This method employing a radiolabeled tracer provides accurate biodegradation rates for freshly contaminated sediments, and represents an upper limit to the natural phenanthrene biodegradation extents if the contaminant is aged over 50 days.