A methodology for regionalizing 3‐D effective porosity at watershed scale in crystalline aquifers

An innovative approach for regionalizing the 3‐D effective porosity field is presented and applied to two large, overexploited, and deeply weathered crystalline aquifers located in southern India. The method derives from earlier work on regionalizing a 2‐D effective porosity field in that part of an aquifer where the water table fluctuates, which is now extended over the entire aquifer using a 3‐D approach. A method based on geological and geophysical surveys has also been developed for mapping the weathering profile layers (saprolite and fractured layers). The method for regionalizing 3‐D effective porosity combines water table fluctuation and groundwater budget techniques at various cell sizes with the use of satellite‐based data (for groundwater abstraction), the structure of the weathering profile, and geostatistical techniques. The approach is presented in detail for the Kudaliar watershed (983 km²) and tested on the 730 km² Anantapur watershed. At watershed scale, the effective porosity of the aquifer ranges from 0.5% to 2% in Kudaliar and between 0.3% and 1% in Anantapur, which agrees with earlier works. Results show that (a) depending on the geology and on the structure of the weathering profile, the vertical distribution of effective porosity can be very different and that the fractured layers in crystalline aquifers are not necessarily characterized by a rapid decrease in effective porosity and (b) that the lateral variations in effective porosity can be larger than the vertical ones. These variations suggest that within a same weathering profile, the density of open fractures and/or degree of weathering in the fractured zone may significantly vary from a place to another. The proposed method provides information on the spatial distribution of effective porosity that is of prime interest in terms of flux and contaminant transport in crystalline aquifers. Implications for mapping groundwater storage and scarcity are also discussed, which should help in improving groundwater resource management strategies.     

A low-dimensional subsurface model for saturated and unsaturated flow processes: ability to address heterogeneity

A low-dimensional model that describes both saturated and unsaturated flow processes in a single equation is presented. Subsurface flow processes in the groundwater, the vadose zone, and the capillary fringe are accounted for through the computation of aggregated hydrodynamic parameters that result from the integration of the governing flow equations from the bedrock to the land surface. The three-dimensional subsurface flow dynamics are thus described by a two-dimensional equation, allowing for a drastic reduction of model unknowns and simplification of the model parameterizations. This approach is compared with a full resolution of the Richards equation in different synthetic test cases. Because the model reduction stems from the vertical integration of the flow equations, the test cases all use different configurations of heterogeneity for vertical cross-sections of a soil-aquifer system. The low-dimensional flow model shows strong consistency with results from a complete resolution of the Richards equation for both the water table and fluxes. The proposed approach is therefore well suited to the accurate reproduction of complex subsurface flow processes.     

Eddy-permitting simulations of the sub-polar North Atlantic: impact of the model bias on water mass properties and circulation

Some previous studies demonstrated that model bias has a strong impact on the quality of long-term prognostic model simulations of the sub-polar North Atlantic Ocean. Relatively strong bias of water mass characteristics is observed in both eddy-permitting and eddy-resolving simulations, suggesting that an increase of model resolution does not reduce significantly the model bias. This study is an attempt to quantify the impact of model bias on the simulated water mass and circulation characteristics in an eddy-permitting model of the sub-polar ocean. This is done through comparison of eddy-permitting prognostic model simulations with the results from two other runs in which the bias is constrained by using spectral nudging. In the first run, the temperature and salinity are nudged towards climatology in the whole column. In the second run, the spectral nudging is applied in the surface 30 m layer and at depths below 560 m only. The biases of the model characteristics of the unconstrained run are similar to those reported in previous eddy-permitting and eddy-resolving studies. The salinity in the surface and intermediate waters of the Labrador Sea waters increases with respect to the climatology, which reduces the stability of the water column. The deep convection in the unconstrained run is artificially intensified and the transport in the sub-polar gyre stronger than in the observations. In particular, the transport of relatively salty and warm Irminger waters into the Labrador Sea is unrealistically high. While the water mass temperature and salinity in the run with spectral nudging in the whole column are closest to the observations, the depth of the winter convection is underestimated in the model. The water mass characteristics and water transport in the run with spectral nudging in the surface and deep layers only are close to observations and at the same time represent well the deep convection in terms of its intensity and position. The source of the bias in the prognostic model run is discussed.     

Short-term impact of reservoir cleaning on the microhabitat use of three non-salmonid fishes in a piedmont river in south west France

As with many rivers of the Pyrenees Mountains, reservoirs on the River Ariège (SW Francc) are used for hydroelectricity generation, irrigation and to maintain summer flows in the River Garonne. Two river reservoirs on the River Ariège were simultaneously emptied in April 1991 to remove accumulated sediment. Fish and ten environmental variables were sampled using Point Abundance Sampling by electrofishing before (February 1991) and one year later (February 1992) at three study sites downstream of the dams. We examined the short-term impact of reservoir cleaning on the relative abundance, habitat availability and habitat use of three species of small fish: minnow Phoxinus phoxinus (L.), gudgeon Gobio gobio (L.) and stone loach Barbatula barbatula (L.).¶After the reservoir cleaning, areas covered by gravel, sand and mud increased significantly as the result of the sediment release. The three species did not react in the same manner to these environmental changes. Stone loach appeared to be the most sensitive species, with its relative abundance decreasing at all sites. Conversely, minnow seemed to be the most resistant species and as such is a leading pioneer species for colonizing new or impacted stretches of rivers channels. Finally, the number of gudgeon decreased particularly at the most downstream site where large deposits of mud were recorded.     

Linking Tarim Basin sea retreat (west China) and Asian aridification in the late Eocene

The Tarim Basin in western China formed the easternmost margin of a shallow epicontinental sea that extended across Eurasia and was well connected to the western Tethys during the Paleogene. Climate modelling studies suggest that the westward retreat of this sea from Central Asia may have been as important as the Tibetan Plateau uplift in forcing aridification and monsoon intensification in the Asian continental interior due to the redistribution of the land‐sea thermal contrast. However, testing of this hypothesis is hindered by poor constraints on the timing and precise palaeogeographic dynamics of the retreat. Here, we present an improved integrated bio‐ and magnetostratigraphic chronological framework of the previously studied marine to continental transition in the southwest Tarim Basin along the Pamir and West Kunlun Shan, allowing us to better constrain its timing, cause and palaeoenvironmental impact. The sea retreat is assigned a latest Lutetian–earliest Bartonian age (ca. 41 Ma; correlation of the last marine sediments to calcareous nannofossil Zone CP14 and correlation of the first continental red beds to the base of magnetochron C18r). Higher up in the continental deposits, a major hiatus includes the Eocene–Oligocene transition (ca. 34 Ma). This suggests the Tarim Basin was hydrologically connected to the Tethyan marine Realm until at least the earliest Oligocene and had not yet been closed by uplift of the Pamir–Kunlun orogenic system. The westward sea retreat at ca. 41 Ma and the disconformity at the Eocene–Oligocene transition are both time‐equivalent with reported Asian aridification steps, suggesting that, consistent with climate modelling results, the sea acted as an important moisture source for the Asian continental interior.     

Momentum and heat transfers in the surface layer over a frozen sea

Wind and temperature profiles from a 10-m mast were measured over the frozen Baltic Sea. The Monin-Obukhov similarity theory applies well to runs carefully selected according to stationarity criteria. This provides relatively low-scattered results with the roughness length z ₀ 0.04 cm, the drag coefficient C _D 1.5 × 10^-3 and the Stanton number C _H 1.00 × 10^-3 under near-neutral conditions. The roughness length, however, behaves in a peculiar way under intense stratification conditions. The reasons proposed for this could lead to an extension of the theory. The mechanisms for momentum and heat transfers are also examined, revealing that there are two different regimes, the smooth and the rough, for the wall friction but not for the heat transfer. Further, the scatter of the data for momentum and heat transfer under the aerodynamically rough regime can be explained by the distinction between the type of thermal stability conditions.     

Tectono‐sedimentary evolution of the western Corinth rift (Central Greece)

The Corinth rift (Greece) is one of the world's most active rifts. The early Plio‐Pleistocene rift is preserved in the northern Peloponnese peninsula, south of the active Corinth rift. Although chronostratigraphic resolution is limited, new structural, stratigraphic and sedimentological data for an area >400 km² record early rift evolution in three phases separated by distinct episodes of extension rate acceleration and northward fault migration associated with major erosion. Minimum total N–S extension is estimated at 6.4–7.7 km. The earliest asymmetrical, broad rift accommodated slow extension (0.6–1 mm a^?1) over >3 Myrs and closed to the west. North‐dipping faults with throws of 1000–2200 m defined narrow blocks (4–7 km) with little footwall relief. A N‐NE flowing antecedent river system infilled significant inherited relief (Lower group). In the earliest Pleistocene, significant fluvial incision coincided with a 15 km northward rift margin migration. Extension rates increased to 2–2.5 mm a^?1. The antecedent rivers then built giant Gilbert‐type fan deltas (Middle group) north into a deepening lacustrine/marine basin. N‐dipping, basin margin faults accommodated throws <1500 m. Delta architecture records initiation, growth and death of this fault system over ca. 800 ka. In the Middle Pleistocene, the rift margin again migrated 5 km north. Extension rate increased to 3.4–4.8 mm a^?1. This transition may correspond to an unconformity in offshore lithostratigraphy. Middle group deltas were uplifted and incised as new hangingwall deltas built into the Gulf (Upper group). A final increase to present‐day extension rates (11–16 mm a^?1) probably occurred in the Holocene. Fault and fault block dimensions did not change significantly with time suggesting control by crustal rheological layering. Extension rate acceleration may be due to strain softening or to regional tectonic factors.     

Orbital scale vegetation change in Africa

Palynological records of Middle and Late Pleistocene marine sediments off African shores is reviewed in order to reveal long-term patterns of vegetation change during climate cycles. Whether the transport of pollen and spores from the source areas on the continent to the ocean floor is mainly by wind or predominantly by rivers depends on the region. Despite the differences in transportation, accumulation rates in the marine sediments decline exponentially with distance to the shore. The marine sediments provide well-dated records presenting the vegetation history of the main biomes of western and southern Africa. The extent of different biomes varied with the climate changes of the glacial interglacial cycle. The Mediterranean forest area expanded during interglacials, the northern Saharan desert during glacials, and the semi-desert area in between during the transitions. In the sub-Saharan mountains ericaceous scrubland spread mainly during glacials and the mountainous forest area often increased during intermediate periods. Savannahs extended or shifted to lower latitudes during glacials. While the representation of the tropical rain forest fluctuated with summer insolation and precession, that of the subtropical biomes showed more obliquity variability or followed the pattern of glacial and interglacials.     

Influence des rejets atmosphériques d'anhydride sulfureux sur la qualité de l'eau des lacs de la région de Rouyn‐Noranda*

Résumé

En mars 1982, le ministère de l'Environnement du Québec a entrepris un échantillonnage physico‐chimique de 64 lacs de tête de la région de Rouyn‐Noranda. Les objectifs de l'étude étaient de déterminer la sensibilité des lacs de cette région et de délimiter l'aire d'influence de la plus importante source anthropique d'émission d'anhydride sulfureux au Québec.

Le traitement graphique et statistique des résultats, ainsi que l'application de modèles empiriques d'acidification ont permis de dégager les zones les plus vulnérables à l'acidification. L'application de ces techniques d'analyse a montré, entre autres, que les concentrations de sulfates en lac étaient beaucoup plus élevées à proximité de Rouyn‐Noranda où est située l'usine Mines Noranda Itée. De fait, la concentration de sulfates dans les lacs de cette région passe de 16 mg L^?1, aux alentours de Rouyn‐Noranda, à moins de 2 mg L^?1 en direction du nord‐est, alors qu'elle baisse à 8 mg L^?1 en bordure de la rivière des Outaouais. De plus, en faisant abstraction de la distribution en forme de cloche des concentrations de sulfates autour de Rouyn‐Noranda, l'étude a permis de faire ressortir la décroissance générale des sulfates du sud‐ouest vers le nord‐est.

Les points majeurs qui ressortent de cette étude sont la présence de lacs très acides à proximité de Rouyn‐Noranda et V extrême vulnérabilité des plans d'eau situés au sud et sud‐est de la ville. En effet, l'eau des lacs de la région du Témiscamingue est en général moins minéralisée et plus acide que celle des lacs au nord de cette région. Cette distinction est due à la nature géologique des substrats dans les bassins versants. Ainsi, les lacs sensibles à l'acidification sont principalement situés sur des bassins versants associés à des roches granitiques et gneissiques, comparativement aux roches volcaniques et aux dépôts d'argiles lacustres qui caractérisent les lacs aux eaux plus minéralisées.     

Transport and deposition of suspended particles in saturated porous media: hydrodynamic effect

A laboratory study was undertaken to determine the transport and deposition rate of suspended particles in columns of saturated porous media (gravel and glass beads), where the porous media were subjected to steady-state flow. Silt particles with a mode of 14 μm diameter (used as the suspended particles) and fluorescein (as the conservative tracer) were injected into the columns in short pulses. The breakthrough curves were competently described with the analytical solution of a convection–dispersion equation with a first-order deposition rate. The experiments were performed using different flow rates. The suspended particle size distribution, the porous media, and the flow rates themselves were the main factors retained in this study to investigate the mechanisms governing the transport and deposition kinetics in detail. The results showed the existence of a flow rate, beyond which suspended particles travel faster than the conservative tracer. A decrease of the deposition rate of suspended particles beyond a critical flow velocity was also observed. Such behaviour led to consideration of the couple hydrodynamic-gravity forces at high flow rates. As the hydrodynamic force increases, particle deposition rates are reduced due to the effect of hydrodynamic forces inhibiting the deposition.