Particulate contribution to extinction of visible radiation: Pollution, haze, and fog

A data set acquired by eight particle-dedicated instruments set up on the SIRTA (Site Instrumental de Recherche par Télédétection Atmosphérique, which is French for Instrumented Site for Atmospheric Remote Sensing Research) during the ParisFog field campaign are exploited to document microphysical properties of particles contributing to extinction of visible radiation in variable situations. The study focuses on a 48-hour period when atmospheric conditions are highly variable: relative humidity changes between 50 and 100%, visibility ranges between 65 and 35 000 m, the site is either downwind the Paris area either under maritime influence. A dense and homogeneous fog formed during the night by radiative cooling. In 6 h, visibility decreased down from 30 000 m in the clear-sky regime to 65 m within the fog, because of advected urban pollution (factor 3 to 4 in visibility reduction), aerosol hydration (factor 20) and aerosol activation (factor 6). Computations of aerosol optical properties, based on Mie theory, show that extinction in clear-sky regime is due equally to the ultrafine modes and to the accumulation mode. Extinction by haze is due to hydrated aerosol particles distributed in the accumulation mode, defined by a geometric mean diameter of 0.6 μm and a geometric standard deviation of 1.4. These hydrated aerosol particles still contribute by 20 ± 10% to extinction in the fog. The complementary extinction is due to fog droplets distributed around the geometric mean diameter of 3.2 μm with a geometric standard deviation of 1.5 during the first fog development stage. The study also shows that the experimental set-up could not count all fog droplets during the second and third fog development stages.     

13C/12C composition, a novel parameter to study the downward migration of paper sludge in soils?

δ ¹³C values of crop and forest soils were measured 8 years after disposal of paper sewage sludge. The carbon transfer from paper sludge downward to the first humic layer is evidenced by a ¹³C-enrichnient of up to + 5.6‰ due to the input of ¹³C-enriched sludge carbonates. ¹³C/¹²C composition is thus a novel, sensitive parameter to follow the downward transfer of paper sludge carbon.     

Potential Flank-Collapse of Soufrière Volcano, Guadeloupe, Lesser Antilles? Numerical Simulation and Hazards

The past history of recurrent flank collapses of la Soufrière volcano of Guadeloupe, its structure, its well-developed hydrothermal system and the current activity constitute factors that could promote a future flank collapse, particularly in the case of a significant increase of activity, with or without shallow magmatic input. To address the hazards associated with such a collapse, we model the emplacement of the debris avalanche generated by a flank-collapse event in 1,250 BC (3,100 years B.P.). We use a finite-difference grain-flow model solving mass and momentum conservation equations that are depth-averaged over the slide thickness, and a Coulomb-type friction law with a variable basal (minimum) friction angle. Using the parameter values determined from this simulation, we then simulate the debris avalanche which could be generated by a potential collapse of the present lava dome. We then discuss the region which could be affected by such a future collapse, and additional associated hazards of concern.     

Bathymay : la structure sous-marine de Mayotte révélée par l'imagerie multifaisceaux

The BATHYMAY marine survey was carried out in January 2004 around the French volcanic island of Mayotte, onboard the research vessel Marion Dufresne 2. Multibeam bathymetry revealed for the first time the morphology of the whole outer slope. Preliminary interpretations show broad canyons that deeply incise the slope, numerous volcanoes and huge landslides. Large submarine plateaus are also highlighted, bounded by steep cliffs and active normal faults. A model for Mayotte is suggested, in which eruptive activity, sometimes influenced by normal faulting, was concentrated along a N140° regional line (the Comoros archipelago axis), with landslide activity concentrated on the northeastern and southwestern flanks of this line. To cite this article: J.-C. Audru et al., C. R. Geoscience 338 (2006).     

Environmental flows in hydro-economic models

The protection of environmental flows, as a management objective for a regulating agency, needs to be consistent with the aquifer water balance and the degree of resource renewability. A stylized hydro-economic model is used where natural recharge, which sustains environmental flows, is considered both in the aquifer water budget and in the welfare function as ecosystem damage. Groundwater recharge and the associated natural drainage may be neglected for aquifers containing fossil water, where the groundwater is mined. However, when dealing with an aquifer that constitutes a renewable resource, for which recharge is not negligible, natural drainage should explicitly appear in the water budget. In doing so, the optimum path of net extraction rate does not necessarily converge to the recharge rate, but depends on the costs associated with ecosystem damages. The optimal paths and equilibrium values for the water volume and water extraction are analytically derived, and numerical simulations based on the Western La Mancha aquifer (southwest Spain) illustrate the theoretical results of the study.     

Comparison of polynomial chaos and Gaussian process surrogates for uncertainty quantification and correlation estimation of spatially distributed open-channel steady flows

Data assimilation is widely used to improve flood forecasting capability, especially through parameter inference requiring statistical information on the uncertain input parameters (upstream discharge, friction coefficient) as well as on the variability of the water level and its sensitivity with respect to the inputs. For particle filter or ensemble Kalman filter, stochastically estimating probability density function and covariance matrices from a Monte Carlo random sampling requires a large ensemble of model evaluations, limiting their use in real-time application. To tackle this issue, fast surrogate models based on polynomial chaos and Gaussian process can be used to represent the spatially distributed water level in place of solving the shallow water equations. This study investigates the use of these surrogates to estimate probability density functions and covariance matrices at a reduced computational cost and without the loss of accuracy, in the perspective of ensemble-based data assimilation. This study focuses on 1-D steady state flow simulated with MASCARET over the Garonne River (South-West France). Results show that both surrogates feature similar performance to the Monte-Carlo random sampling, but for a much smaller computational budget; a few MASCARET simulations (on the order of 10–100) are sufficient to accurately retrieve covariance matrices and probability density functions all along the river, even where the flow dynamic is more complex due to heterogeneous bathymetry. This paves the way for the design of surrogate strategies suitable for representing unsteady open-channel flows in data assimilation.     

Using stable water isotopes to identify spatio-temporal controls on groundwater recharge in two contrasting East African aquifer systems

Understanding the spatio-temporal variability in groundwater recharge is a prerequisite to sustainable management of aquifers. Spatial analysis of groundwater stable isotopes uncovered predominant controls on groundwater recharge in the Nairobi aquifer system (NAS) and the South Coast aquifer (SC), two exemplar East African aquifers relied upon by 7 million people. We analysed 368 samples for stable isotopes and basic physico-chemical parameters. The NAS groundwater isotopes are controlled by precipitation orographic effects and enriched recharge from impounded lakes/wetlands; the SC isotopes are correlated with water-table depth influencing evapotranspiration. Global Network of Isotopes in Precipitation (GNIP) data revealed groundwater recharge during months of heavy rains in the NAS, whilst the SC experiences spatio-temporally diffuse recharge. Inferred “isoscapes” show: in the NAS, (1) direct, rapid recharge favoured by faults, well-drained soils and ample rainfall in uplands, (2) delayed recharge from impounded lakes and wetlands in mid-lands, and (3) focused, event-based recharge in floodplains; and in the SC, diffuse recharge complicated by significant water-table evapotranspiration processes.     

Pre-Permian uplift and diffuse extensional deformation in the High Zagros Belt (Iran): integration in the geodynamic evolution of the Arabian plate

The High Zagros Belt includes exposures of Lower Palaeozoic rocks in the core of several thrusted anticlines developed during the Cenozoic Zagros orogeny. The structural style and tectonic evolution of this area during the Palaeozoic remain poorly understood due to the complexity of the subsequent deformation. We present the preliminary results of a field study focusing on the structural geology of Palaeozoic rocks in this area. We confirm the existence of an angular unconformity below the Lower Permian Faraghan formation. In the geological literature, this unconformity is reported as the “Hercynian Unconformity” suggesting a relationship with the Hercynian (Variscan) orogeny, which affected Western Europe and westernmost Africa during the Carboniferous. Surprisingly the only observable structures sealed by this unconformity are N to NE trending normal faults and tilted blocks without any evidence of compressional deformation. This pre-Permian extensional deformation, which is general at the scale of the HZB, raises questions about the geodynamic significance of the “Hercynian unconformity” in the study area and, more generally, in the Arabian plate.     

Formation and preservation of pedogenic carbonates in South India, links with paleo-monsoon and pedological conditions: Clues from Sr isotopes, U-Th series and REEs

The influence of the pedogenic and climatic contexts on the formation and preservation of pedogenic carbonates in a climosequence in the Western Ghats (Karnataka Plateau, South West India) has been studied. Along the climosequence, the current mean annual rainfall (MAR) varies within a 80 km transect from 6000 mm at the edge of the Plateau to 500 mm inland. Pedogenic carbonates occur in the MAR range of 500-1200 mm. In the semi-arid zone (MAR: 500-900 mm), carbonates occur (i) as thick hardpan calcretes on pediment slopes and (ii) as nodular horizons in polygenic black soils (i.e. vertisols). In the sub-humid zone (MAR: 900-1500 mm), pedogenic carbonates are disseminated in the black soil matrices either as loose, irregular and friable nodules of millimetric size or as indurated botryoidal nodules of centimetric to pluricentimetric size. They also occur at the top layers of the saprolite either as disseminated pluricentimetric indurated nodules or carbonate-cemented lumps of centimetric to decimetric size.Chemical and isotopic (⁸⁷Sr/⁸⁶Sr) compositions of the carbonate fraction were determined after leaching with 0.25 N HCl. The corresponding residual fractions containing both primary minerals and authigenic clays were digested separately and analyzed. The trend defined by the ⁸⁷Sr/⁸⁶Sr signatures of both labile carbonate fractions and corresponding residual fractions indicates that a part of the labile carbonate fraction is genetically linked to the local soil composition. Considering the residual fraction of each sample as the most likely lithogenic source of Ca in carbonates, it is estimated that from 24% to 82% (55% on average) of Ca is derived from local bedrock weathering, leading to a consumption of an equivalent proportion of atmospheric CO₂. These values indicate that climatic conditions were humid enough to allow silicate weathering: MAR at the time of carbonate formation likely ranged from 400 to 700 mm, which is 2- to 3-fold less than the current MAR at these locations.The Sr, U and Mg contents and the (²³⁴U/²³⁸U) activity ratio in the labile carbonate fraction help to understand the conditions of carbonate formation. The relatively high concentrations of Sr, U and Mg in black soil carbonates may indicate fast growth and accumulation compared to carbonates in saprolite, possibly due to a better confinement of the pore waters which is supported by their high (²³⁴U/²³⁸U) signatures, and/or to higher content of dissolved carbonates in the pore waters. The occurrence of Ce, Mn and Fe oxides in the cracks of carbonate reflects the existence of relatively humid periods after carbonate formation. The carbonate ages determined by the U-Th method range from 1.33 ± 0.84 kyr to 7.5 ± 2.7 kyr and to a cluster of five ages around 20 kyr, i.e. the Last Glacial Maximum period. The young occurrences are only located in the black soils, which therefore constitute sensitive environments for trapping and retaining atmospheric CO₂ even on short time scales. The maximum age of carbonates depends on their location in the climatic gradient: from about 20 kyr for centimetric nodules at Mule Hole (MAR = 1100 mm/yr) to 200 kyr for the calcrete at Gundlupet (MAR = 700 mm/yr, Durand et al., 2007). The intensity of rainfall during wet periods would indeed control the lifetime of pedogenic carbonates and thus the duration of inorganic carbon storage in soils.     

Estimating aquifer thickness using multiple pumping tests

A method to estimate aquifer thickness and hydraulic conductivity has been developed, consisting of multiple pumping tests. The method requires short-duration pumping cycles on an unconfined aquifer with significant seasonal water-table fluctuations. The interpretation of several pumping tests at a site in India under various initial conditions provides information on the change in hydrodynamic parameters in relation to the initial water-table level. The transmissivity linearly decreases compared with the initial water level, suggesting a homogeneous distribution of hydraulic conductivity with depth. The hydraulic conductivity is estimated from the slope of this linear relationship. The extrapolation of the relationship between transmissivity and water level provides an estimate of the aquifer thickness that is in good agreement with geophysical investigations. The hydraulically active part of the aquifer is located in both the shallow weathered and the underlying densely fractured zones of the crystalline basement. However, no significant relationship is found between the aquifer storage coefficient and initial water level. This new method contributes to filling the methodological gap between single pumping tests and hydraulic tomography, in providing information on the variation of the global transmissivity according to depth. It can be applied to any unconfined aquifer experiencing large seasonal water-table fluctuations and short pumping cycles.