Lichens used as monitors of atmospheric pollution around Agadir (Southwestern Morocco) - a case study predating lead-free gasoline

More than 30 epiphytic lichens, collected in Agadir (Morroco) and along a 150-km transect from the Atlantic Ocean eastward, were analyzed for their metal content and lead isotopic composition. This dataset was used to evaluate atmospheric metal contamination and the impact of the city on the surrounding area. The concentrations of Cu, Pb, and Zn (average ± 1 SD) were 20.9 ± 15.2 μg g−1, 13.8 ± 9.0 μg g−1, and 56.6 ± 26.6 μg g−1, respectively, with the highest values observed in lichens collected within the urban area. The 206Pb/207Pb and 208Pb/207Pb ratios in the lichens varied from 1.146 to 1.186 and from 2.423 to 2.460, respectively. Alkyllead-gasoline sold in Morocco by the major petrol companies gave isotopic ratios of 206Pb/207Pb = 1.076–1.081 and 208Pb/207Pb = 2.348–2.360. These new, homogeneous values for gasoline-derived lead improve and update the scarce isotopic database of potential lead sources in Morocco, and may be of great value to future environmental surveys on the presence of lead in natural reservoirs, where it persists over time (e.g., soils and sediments). The interest of normalizing metal concentrations in lichens to concentrations of a lithogenic element is demonstrated by the consistency of the results thus obtained with lead isotopic ratios. Leaded gasoline contributed less than 50% of the total amount of lead accumulated in lichens, even in areas subject to high vehicular traffic. This strongly suggests that the recent banishment of leaded gasoline in Morocco will not trigger a drastic improvement in air quality, at least in Agadir.

Post-frying oil drainage from potato chips and french fries: a comparative study of atmospheric and vacuum drainage

This paper represents a study of the transient changes occurring in temperature, and moisture and oil contents during the so called “post-frying drainage”—which is the duration for which a product is held in the head space of the fryer after it is removed from the oil. Since most of the oil adhering to the product penetrates into the structure during this period, this paper examines the effects of applying vacuum during drainage (1.33 kPa) to maintain the product temperature consistently above the water saturation temperature corresponding to the prevailing pressure (11 °C), which potentially eliminates water condensation and prevents the occluded surface oil from penetrating into the product structure. Draining under vacuum significantly lowers the oil content of potato chips by 38% compared to atmospheric drainage. This phenomenon can be further confirmed by confocal laser scanning microscopy (CLSM) images, which show that the boundary between the core and the crust regions is clearly visible in the case of vacuum drainage, whereas in the case of atmospheric drainage, the oil is distributed throughout the structure. Unfortunately, the same approach did not reduce the oil content of French fries—the oil content of vacuum-drained product was found similar to the product obtained by draining under atmospheric pressure. This is because the reduction in oil content only occurs when there is net moisture evaporation from the product and the evaporation rate is sufficient to force out the oil from the product; this was clearly not the case with French fries. The CLSM images show that the oil distribution in the products drained under atmospheric pressure and vacuum was similar.

The effects of a changing pollution climate on throughfall deposition and cycling in a forested area in southern England

This study compares two sets of measurements of the composition of bulk precipitation and throughfall at a site in southern England with a 20-year gap between them. During this time, SO2 emissions from the UK fell by 82%, NOx emissions by 35% and NH3 emissions by 7%. These reductions were partly reflected in bulk precipitation, with deposition reductions of 56% in SO4,38% in NO3, 32% in NH4, and 73% in H+. In throughfall under Scots pine, the effects were more dramatic, with an 89% reduction in SO4 deposition and a 98% reduction in H+ deposition. The mean pH under these trees increased from 2.85 to 4.30. Nitrate and ammonium deposition in throughfall increased slightly, however. In the earlier period, the Scots pines were unable to neutralise the high flux of acidity associated with sulphur deposition, even though this was not a highly polluted part of the UK, and deciduous trees (oak and birch) were only able to neutralise it in summer when the leaves were present. In the later period, the sulphur flux had reduced to the point where the acidity could be neutralised by all species — the neutralisation mechanism is thus likely to be largely leaching of base cations and buffering substances from the foliage. The high fluxes are partly due to the fact that these are 60–80 year old trees growing in an open forest structure. The increase in NO3 and NH4 in throughfall in spite of decreased deposition seems likely due to a decrease in foliar uptake, perhaps due to the increasing nitrogen saturation of the catchment soils. These changes may increase the rate of soil microbial activity as nitrogen increases and acidity declines, with consequent effects on water quality of the catchment drainage stream.

Advances in sequential data assimilation and numerical weather forecasting: an Ensemble Transform Kalman-Bucy Filter, a study on clustering in deterministic Ensemble Square Root Filters, and a test of a new time stepping scheme in an atmospheric model

This dissertation deals with aspects of sequential data assimilation (in particular ensemble Kalman filtering) and numerical weather forecasting. In the first part, the recently formulated Ensemble Kalman-Bucy (EnKBF) filter is revisited. It is shown that the previously used numerical integration scheme fails when the magnitude of the background error covariance grows beyond that of the observational error covariance in the forecast window. Therefore, we present a suitable integration scheme that handles the stiffening of the differential equations involved and doesn’t represent further computational expense. Moreover, a transform-based alternative to the EnKBF is developed: under this scheme, the operations are performed in the ensemble space instead of in the state space. Advantages of this formulation are explained. For the first time, the EnKBF is implemented in an atmospheric model. The second part of this work deals with ensemble clustering, a phenomenon that arises when performing data assimilation using of deterministic ensemble square root filters in highly nonlinear forecast models. Namely, an M-member ensemble detaches into an outlier and a cluster of M-1 members. Previous works may suggest that this issue represents a failure of EnSRFs; this work dispels that notion. It is shown that ensemble clustering can be reverted also due to nonlinear processes, in particular the alternation between nonlinear expansion and compression of the ensemble for different regions of the attractor. Some EnSRFs that use random rotations have been developed to overcome this issue; these formulations are analyzed and their advantages and disadvantages with respect to common EnSRFs are discussed. The third and last part contains the implementation of the Robert-Asselin-Williams (RAW) filter in an atmospheric model. The RAW filter is an improvement to the widely popular Robert-Asselin filter that successfully suppresses spurious computational waves while avoiding any distortion in the mean value of the function. Using statistical significance tests both at the local and field level, it is shown that the climatology of the SPEEDY model is not modified by the changed time stepping scheme; hence, no retuning of the parameterizations is required. It is found the accuracy of the medium-term forecasts is increased by using the RAW filter.

Tropical cyclone climatology in a 10-km global atmospheric GCM: toward weather-resolving climate modeling

Northern Hemisphere tropical cyclone (TC) activity is investigated in multiyear global climate simulations with theECMWFIntegrated Forecast System (IFS) at 10-km resolution forced by the observed records of sea surface temperature and sea ice. The results are compared to analogous simulationswith the 16-, 39-, and 125-km versions of the model as well as observations. In the North Atlantic, mean TC frequency in the 10-km model is comparable to the observed frequency, whereas it is too low in the other versions. While spatial distributions of the genesis and track densities improve systematically with increasing resolution, the 10-km model displays qualitatively more realistic simulation of the track density in the western subtropical North Atlantic. In the North Pacific, the TC count tends to be too high in thewest and too low in the east for all resolutions. These model errors appear to be associated with the errors in the large-scale environmental conditions that are fairly similar in this region for all model versions. The largest benefits of the 10-km simulation are the dramatically more accurate representation of the TC intensity distribution and the structure of the most intense storms. The model can generate a supertyphoon with a maximum surface wind speed of 68.4 m s21. The life cycle of an intense TC comprises intensity fluctuations that occur in apparent connection with the variations of the eyewall/rainband structure. These findings suggest that a hydrostatic model with cumulus parameterization and of high enough resolution could be efficiently used to simulate the TC intensity response (and the associated structural changes) to future climate change.

Climate change and coupling of macronutrient cycles along the atmospheric, terrestrial, freshwater and estuarine continuum

This paper provides an introduction to the Special Issue on “Climate Change and Coupling of Macronutrient Cycles along the Atmospheric, Terrestrial, Freshwater and Estuarine Continuum”, dedicated to Colin Neal on his retirement. It is not intended to be a review of this vast subject, but an attempt to synthesize some of the major findings from the 22 contributions to the Special Issue in the context of what is already known. The major research challenges involved in understanding coupled macronutrient cycles in these environmental media are highlighted, and the difficulties of making credible predictions of the effects of climate change are discussed. Of particular concern is the possibility of interactions which will enhance greenhouse gas concentrations and provide positive feedback to global warming.

Physicochemical changes occuring in oil when atmospheric frying is combined with post-frying vacuum application

The aim of this study was to investigate the effect of atmospheric frying followed by drainage under vacuum on the stability of oil, compared to similar frying with drainage at atmospheric pressure. Changes in the oil were assessed by the free fatty acid (FFA) content, p-anisidine value (AnV), colour, viscosity, fatty acid profile and concentration of tocols. The rate of FFA formation in the case of vacuum drainage was found to be about half that of atmospheric drainage. Oil deterioration by oxidation and polymerisation was also reduced by the use of vacuum drainage. The AnV of the oil after vacuum drainage was lower by about 12%, the total colour difference was improved by 14% and viscosity was slightly reduced after 5 days of frying, compared to the values for oil that had been drained at atmospheric pressure. There was a reduction in the loss of polyunsaturated fatty acids in the case of vacuum drainage after 5 days of frying but differences in retention of tocols were only evident in the first two days of frying.

On the momentum fluxes associated with mountain waves in directionally sheared flows

The direct impact of mountain waves on the atmospheric circulation is due to the deposition of wave momentum at critical levels, or levels where the waves break. The first process is treated analytically in this study within the framework of linear theory. The variation of the momentum flux with height is investigated for relatively large shears, extending the authors’ previous calculations of the surface gravity wave drag to the whole atmosphere. A Wentzel–Kramers–Brillouin (WKB) approximation is used to treat inviscid, steady, nonrotating, hydrostatic flow with directional shear over a circular mesoscale mountain, for generic wind profiles. This approximation must be extended to third order to obtain momentum flux expressions that are accurate to second order. Since the momentum flux only varies because of wave filtering by critical levels, the application of contour integration techniques enables it to be expressed in terms of simple 1D integrals. On the other hand, the momentum flux divergence (which corresponds to the force on the atmosphere that must be represented in gravity wave drag parameterizations) is given in closed analytical form. The momentum flux expressions are tested for idealized wind profiles, where they become a function of the Richardson number (Ri). These expressions tend, for high Ri, to results by previous authors, where wind profile effects on the surface drag were neglected and critical levels acted as perfect absorbers. The linear results are compared with linear and nonlinear numerical simulations, showing a considerable improvement upon corresponding results derived for higher Ri.

Assessing wind profile effects on the global atmospheric torque

The impact of a new approach to the evaluation of surface gravity wave drag (GWD) is assessed. This approach uses linear theory, but incorporates the effects of wind profile shear and curvature, by means of a second-order WKB approximation. While the theory predicts the possibility of either drag enhancement or reduction, depending on the wind profile, results obtained with the ERA-40 reanalysis data clearly indicate the predominance of local drag enhancement. However, the global impact of shear on the atmospheric axial GWD torque comes mostly from regions with predominantly easterly flow, contributing to a slight reduction of the bias found in different studies of the global angular momentum budget. The relative correction due to shear on linear GWD is found not to depend too strongly on the levels chosen for the computation of the low-level wind derivatives.

Effect of mountain breeze in ventilating the city under calm and neutral atmospheric environment: Interaction of airflow structures and ventilation efficiency

Many modern cities locate in the mountainous areas, like Hong Kong, Phoenix City and Los Angles. It is confirmed in the literature that the mountain wind system developed by differential heating or cooling can be very beneficial in ventilating the city nearby and alleviating the UHI effect. However, the direct interaction of mountain wind with the natural-convection circulation due to heated urban surfaces has not been studied, to our best knowledge. This kind of unique interaction of two kinds of airflow structures under calm and neutral atmospheric environment is investigated in this paper by CFD approach. A physical model comprising a simple mountain and three long building blocks (forming two street canyons) is firstly developed. Different airflow structures are identified within the conditions of different mountain-building height ratios (R=Hm/Hb) by varying building height but fixing mountain height. It is found that the higher ventilation rate in the street canyons is expected in the cases of smaller mountain-building ratios, indicating the stronger natural convection due to increasing heated building surfaces. However, there is the highest air change rate (ACH) in the lowest-building-height case and most of the air is advective into the street canyon through the top open area, highlighting the important role played by the mountain wind. In terms of the ventilation efficiency, it is shown that the smallest R case enjoys the best air change efficiency followed by the highest R case, while the worst ventilative street canyons occur at the middle R case. In the end, a gap across the streets is introduced in the modeling. The existence of the gap can greatly channel the mountain wind and distribute the air into streets nearby. Thus the ACH can be doubled and air quality can be significantly improved.

Arctic Oscillation and Antarctic Oscillation in internal atmospheric variability

n this study, we investigated the features of Arctic Oscillation (AO) and Antarctic Oscillation (AAO), that is, the annular modes in the extratropics, in the internal atmospheric variability attained through an ensemble of integrations by an atmospheric general circulation model (AGCM) forced with the global observed SSTs. We focused on the interannual variability of AO/AAO, which is dominated by internal atmospheric variability. In comparison with previous observed results, the AO/AAO in internal atmospheric variability bear some similar characteristics, but exhibit a much clearer spatial structure: significant correlation between the North Pacific and North Atlantic centers of action, much stronger and more significant associated precipitation anomalies, and the meridional displacement of upper-tropospheric westerly jet streams in the Northern/Southern Hemisphere. In addition, we examined the relationship between the North Atlantic Oscillation (NAO)/AO and East Asian winter monsoon (EAWM). It has been shown that in the internal atmospheric variability, the EAWM variation is significantly related to the NAO through upper-tropospheric atmospheric teleconnection patterns.

Latitudinal variation of the effect of aviation NOx emissions on atmospheric ozone and methane and related climate metrics

We evaluate the response to regional and latitudinal changes in aircraft NOx emissions using several climate metrics (radiative forcing (RF), Global Warming Potential (GWP), Global Temperature change Potential (GTP)). Global chemistry transport model integrations were performed with sustained perturbations in regional aircraft and aircraft-like NOx emissions. The RF due to the resulting ozone and methane changes is then calculated. We investigate the impact of emission changes for specific geographical regions (approximating to USA, Europe, India and China) and cruise altitude emission changes in discrete latitude bands covering both hemispheres. We find that lower latitude emission changes (per Tg N) cause ozone and methane RFs that are about a factor of 6 larger than those from higher latitude emission changes. The net RF is positive for all experiments. The meridional extent of the RF is larger for low latitude emissions. GWPs for all emission changes are positive, with tropical emissions having the largest values; the sign of the GTP depends on the choice of time horizon.

Copper deposition on TiO2 from copper(II)hexafluoroacetylacetonate

The authors have studied the adsorption of CuII(hfac)2 on the surface of a model oxide system, TiO2(110), and probed the molecular stability with respect to thermal cycling, using atomic scale imaging by scanning tunneling microscopy supported by x-ray photoemission spectroscopy. They find that at 473 K, the adsorbed metal-organic molecules begin to dissociate and release Cu atoms which aggregate and form Cu nanoparticles. These Cu nanoparticles ripen over time and the size (height) distribution develops into a bimodal distribution. Unlike other organometallic systems, which show a bimodal distribution due to enhanced nucleation or growth at surface step edges, the nanoparticles do not preferentially form at steps. The reduced mobility of the Cu islands may be related to the co-adsorbed ligands that remain in very small clusters on the surface.

Correction of atmospheric delay effects in radar interferometry using a nested mesoscale atmospheric model

The temporal variability of the atmosphere through which radio waves pass in the technique of differential radar interferometry can seriously limit the accuracy with which the method can measure surface motion. A forward, nested mesoscale model of the atmosphere can be used to simulate the variable water content along the radar path and the resultant phase delays. Using this approach we demonstrate how to correct an interferogram of Mount Etna in Sicily associated with an eruption in 2004-5. The regional mesoscale model (Unified Model) used to simulate the atmosphere at higher resolutions consists of four nested domains increasing in resolution (12, 4, 1, 0.3 km), sitting within the analysis version of a global numerical model that is used to initiate the simulation. Using the high resolution 3D model output we compute the surface pressure, temperature and the water vapour, liquid and solid water contents, enabling the dominant hydrostatic and wet delays to be calculated at specific times corresponding to the acquisition of the radar data. We can also simulate the second-order delay effects due to liquid water and ice.

Cloud droplet size and liquid water path retrievals from zenith radiance measurements: examples from the Atmospheric Radiation Measurement Program and the Aerosol Robotic Network

The ground-based Atmospheric Radiation Measurement Program (ARM) and NASA Aerosol Robotic Net- work (AERONET) routinely monitor clouds using zenith ra- diances at visible and near-infrared wavelengths. Using the transmittance calculated from such measurements, we have developed a new retrieval method for cloud effective droplet size and conducted extensive tests for non-precipitating liquid water clouds. The underlying principle is to combine a liquid-water-absorbing wavelength (i.e., 1640 nm) with a non-water-absorbing wavelength for acquiring information on cloud droplet size and optical depth. For simulated stratocumulus clouds with liquid water path less than 300 g m−2 and horizontal resolution of 201 m, the retrieval method underestimates the mean effective radius by 0.8μm, with a root-mean-squared error of 1.7 μm and a relative deviation of 13%. For actual observations with a liquid water path less than 450 g m−2 at the ARM Oklahoma site during 2007– 2008, our 1.5-min-averaged retrievals are generally larger by around 1 μm than those from combined ground-based cloud radar and microwave radiometer at a 5-min temporal resolution. We also compared our retrievals to those from combined shortwave flux and microwave observations for relatively homogeneous clouds, showing that the bias between these two retrieval sets is negligible, but the error of 2.6 μm and the relative deviation of 22 % are larger than those found in our simulation case. Finally, the transmittance-based cloud effective droplet radii agree to better than 11 % with satellite observations and have a negative bias of 1 μm. Overall, the retrieval method provides reasonable cloud effective radius estimates, which can enhance the cloud products of both ARM and AERONET.