Effects of the geometry of two‐dimensional fractures on their hydraulic aperture and on the validity of the local cubic law

Flow through rough fractures is investigated numerically in order to assess the validity of the local cubic law for different fracture geometries. Two‐dimensional channels with sinusoidal walls having different geometrical properties defined by the aperture, the amplitude, and the wavelength of the walls' corrugations, the corrugations asymmetry, and the phase shift between the two walls are considered to represent different fracture geometries. First, it is analytically shown that the hydraulic aperture clearly deviates from the mean aperture when the walls' roughness, the phase shift, and/or the asymmetry between the fracture walls are relatively high. The continuity and the Navier–Stokes equations are then solved by means of the finite element method and the numerical solutions compared to the theoretical predictions of the local cubic law. Reynolds numbers ranging from 0.066 to 66.66 are investigated so as to focus more particularly on the effect of flow inertial effects on the validity of the local cubic law. For low Reynolds number, typically less than 15, the local cubic law properly describes the fracture flow, especially when the fracture walls have small corrugation amplitudes. For Reynolds numbers higher than 15, the local cubic law is valid under the conditions that the fracture presents a low aspect ratio, small corrugation amplitudes, and a moderate phase lag between its walls.     

Tidal effects on temperature front in the Yellow Sea

Temperature front (TF) is one of the important features in the Yellow Sea, which forms in spring, thrives in summer, and fades in autumn as thermocline declines. TF intensity ⋎S _T ⋎ is defined to describe the distribution of TF. Based on the MASNUM wave-tide-circulation coupled model, temperature distribution in the Yellow Sea was simulated with and without tidal effects. Along 36°N, distribution of TF from the simulated results are compared with the observations, and a quantitative analysis is introduced to evaluate the tidal effects on the forming and maintaining processes of the TF. Tidal mixing and the circulation structure adapting to it are the main causes of the TF. Supported by the National Basic Research Program of China (No. G1999043809) and the National Science Foundation of China (No. 49736190).     

The hydrological influence of black locust plantations in the loess area of northwest China

Black locust (Robina pseudoacacia) has become one of the most important shelter species in the loess area of northwest China. This paper summarizes recent research concerning its hydrological influence, including canopy interception, litter absorption capacity, its effect on rainfall kinetic energy, infiltration rates, surface runoff, soil moisture, and evapotranspiration, and its role in soil conservation. Several predictive models are listed. on the basis of existing results, optimum characteristics for an effective plantation are defined, and problems requiring further research are identified.     

Alluvial architecture of the Holocene Rhine–Meuse delta (the Netherlands)

Ancient fluvial successions often act as hydrocarbon reservoirs. Sub‐surface data on the alluvial architecture of fluvial successions are often incomplete and modelling is performed to reconstruct the stratigraphy. However, all alluvial architecture models suffer from the scarcity of field data to test and calibrate them. The purposes of this study were to quantify the alluvial architecture of the Holocene Rhine–Meuse delta (the Netherlands) and to determine spatio‐temporal trends in the architecture. Five north–south orientated cross‐sections, perpendicular to the general flow direction, were compiled for the fluvial‐dominated part of the delta. These sections were used to calculate the width/thickness ratios of fluvial sandbodies (SBW/SBT) and the proportions of channel‐belt deposits (CDP), clastic overbank deposits (ODP) and organic material (OP) in the succession. Furthermore, the connectedness ratio (CR) between channel belts was calculated for each cross‐section. Distinct spatial and temporal trends in the alluvial architecture were found. SBW/SBT ratios decrease by a factor of ca 4 in a downstream direction. CDP decreases from ca 0·7 (upstream) to ca 0·3 (downstream). OP increases from less than 0·05 in the upstream part of the delta to more than 0·25 in the downstream delta. ODP is approximately constant (0·4). CR is ca 0·25 upstream, which is approximately two times larger than in the downstream part of the delta. Furthermore, CDP in the downstream Rhine–Meuse delta increases after 3000 cal yr BP. These trends are attributed to variations in available accommodation space, floodplain geometry and channel‐belt size. For instance, channel belts tend to narrow in a downstream direction, which reduces SBW/SBT, CDP and CR. Tectonics cause local deviations in the general architectural trends. In addition, the positive correlation between avulsion frequency and the ratio of local to regional aggradation rate probably influenced alluvial architecture in the Rhine–Meuse delta. The Rhine–Meuse data set can be a great resource when developing more sophisticated models for alluvial architecture simulation, which eventually could lead to better characterizations of hydrocarbon reservoirs. To aid such usage of the Rhine–Meuse data set, constraints for relevant parameters are provided at the end of the paper.     

Morphological analysis of the drainage system in the Eastern Alps

We study the morphology of the major rivers draining the Eastern Alps to test whether the active tectonics of this part of the orogen is reflected in the shape of channel profiles of the river network. In our approach we compare channel profiles measured from digital elevation models with numerically modelled channel profiles using a stream power approach. It is shown that regions of high stream power coincide largely with regions of highest topography and largest uplift rates, while the forelands and the Pannonian Basin are characterised by a significantly lower stream power. From stream power modelling we conclude that there is young uplift at the very east of the Eastern Alps, in the Bohemian Massif and in the Pohorje Range. The impact of the Pleistocene glaciations is explored by comparing properties of rivers that drain in proximal and distal positions relative to the ice sheet during the last glacial maximum. Our analysis shows that most knick points, wind gaps and other non-equilibrium features of catchments covered by ice during the last glaciations (Salzach, Enns) can be correlated with glacial processes. In contrast the ice free catchments of the Mur and Drava are characterized by channels in morphological equilibrium at the first approximation and are showing only weak evidence of the strong tectonic activity within these catchments. Finally, the channel profiles of the Adige and the divide between the upper Rhine and Danube catchments differ significantly from the other catchments. We relate this to the fact that the Adige and the Rhine respond to different base levels from the remainder of the Eastern Alps: The Adige may preserve a record from the Messininan base level change and the Rhine is subject to the base level lowering in the Rhine Graben.     

Migration of the Ganga river and its implication on hydro-geological potential of Varanasi area, U.P., India

Borehole data reveals that during Late Quaternary, the Ganga river was non-existent in its present location near Varanasi. Instead, it was flowing further south towards peripheral craton. Himalayan derived grey micaceous sands were being carried by southward flowing rivers beyond the present day water divide of Ganga and mixed with pink arkosic sand brought by northward flowing peninsular rivers. Subsequently, the Ganga shifted to its present position and got incised. Near Varanasi, the Ganga river is flowing along a NW-SE tectonic lineament. The migration of Ganga river is believed to have been in response to basin expansion caused due to Himalayan tectonics during Middle Pleistocene times. Multi-storied sand bodies generated as a result of channel migration provide excellent aquifers confined by a thick zone of muddy sediments near the surface. Good quality potable water is available at various levels below about 70 m depth in sandy aquifers. Craton derived gravelly coarse-to-medium grained sand forms the main aquifer zones of tens of meter thickness with enormous yield. In contrast, the shallow aquifers made up of recycled interfluve silt and sandy silt occur under unconfined conditions and show water-level fluctuation of a few meters during pre-and post-monsoon periods.     

Application of MODFLOW and geographic information system to groundwater flow simulation in North China Plain, China

MODFLOW is a groundwater modeling program. It can be compiled and remedied according to the practical applications. Because of its structure and fixed data format, MODFLOW can be integrated with Geographic Information Systems (GIS) technology for water resource management. The North China Plain (NCP), which is the politic, economic and cultural center of China, is facing with water resources shortage and water pollution. Groundwater is the main water resource for industrial, agricultural and domestic usage. It is necessary to evaluate the groundwater resources of the NCP as an entire aquifer system. With the development of computer and internet information technology it is also necessary to integrate the groundwater model with the GIS technology. Because the geological and hydrogeological data in the NCP was mainly in MAPGIS format, the powerful function of GIS of disposing of and analyzing spatial data and computer languages such as Visual C and Visual Basic were used to define the relationship between the original data and model data. After analyzing the geological and hydrogeological conditions of the NCP, the groundwater flow numerical simulation modeling was constructed with MODFLOW. On the basis of GIS, a dynamic evaluation system for groundwater resources under the internet circumstance was completed. During the process of constructing the groundwater model, a water budget was analyzed, which showed a negative budget in the NCP. The simulation period was from 1 January 2002 to 31 December 2003. During this period, the total recharge of the groundwater system was 49,374 × 10⁶ m³ and the total discharge was 56,530 × 10⁶ m³ the budget deficit was −7,156 × 10⁶ m³. In this integrated system, the original data including graphs and attribution data could be stored in the database. When the process of evaluating and predicting groundwater flow was started, these data were transformed into files that the core program of MODFLOW could read. The calculated water level and drawdown could be displayed and reviewed online.     

Impacts of at-site wastewater disposal systems on the groundwater aquifer in arid regions: case of Sfax City, Southern Tunisia

Groundwater in Sfax City (Tunisia) has been known since the beginning of the century for its deterioration in quality, as a result of wastewater recharge into the aquifer. An average value of 12 × 10⁶ m³ of untreated wastewater reaches the groundwater aquifer each year. This would result not only in a chemical and biological contamination of the groundwater, but also in an increase of the aquifer piezometric level. Quantitative impacts were evaluated by examining the groundwater piezometric level at 57 surface wells and piezometers. The survey showed that, during the last two decades, the groundwater level was ever increasing in the urban area with values reaching 7 m in part; and decreasing in Sidi Abid (agricultural area) with values exceeding −3 m. Groundwater samples for chemical and microbial analysis were collected from 41 wells spread throughout the study area. Results showed significantly elevated levels of sodium, chlorides, nitrates and coliform bacteria all over the urban area. High levels (NO₃: 56–254 mg/l; Na >1,500 mg/l; Coliforms >30/100 ml) can be related to more densely populated areas with a higher density of pit latrine and recharge wells. Alternatively results showed a very variable chemical composition of groundwater, e.g. electrical conductivity ranges from 4,040 to19,620 μs/cm and the dry residual varies between 1.4 and 14 g/l with concentrations increasing downstream. Furthermore a softening of groundwater in Set Ezzit (highly populated sector) was observed.     

Hydrogeochemical investigations and groundwater provinces of the Friuli Venezia Giulia Plain aquifers,northeastern Italy

Water resources are a key factor, particularly for the planning of the sustainable regional development of agriculture, as well as for socio-economic development in general. A hydrochemical investigation was conducted in the Friuli Venezia Giulia aquifer systems to identify groundwater evolution, recharge and extent of pollution. Temperature, pH, electric conductivity, total dissolved solids, alkalinity, total hardness, SAR, Ca²⁺, Na⁺, K⁺, Mg²⁺, Cl⁻, SO₄ ²⁻, NO₃ ⁻, HCO₃ ⁻ _, water quality and type, saturation indexes and the environmental stable isotope δ¹⁸O were determined in 149 sampling stations. The pattern of geochemical and oxygen stable isotope variations suggests that the sub-surface groundwater (from phreatic and shallow confined aquifers) is being recharged by modern precipitations and local river infiltrations. Four hydrogeological provinces have been recognised and mapped in the Friuli Venezia Giulia Plain having similar geochemical signatures. These provinces have different degrees of vulnerability to contamination. The deep confined groundwater samples are significantly less impacted by surface activities; and it appears that these important water resources have very low recharge rates and would, therefore, be severely impacted by overabstraction.     

Glacier change and glacier runoff variation in the Tuotuo River basin,the source region of Yangtze River in western China

Glaciers in the Tuotuo River basin, western China, have been monitored in recent decades by applying topographical maps and high-resolution satellite images. Results indicate that most of glaciers in the Tuotuo River basin have retreated in the period from 1968/1971 to 2001/2002, and their shrinkage area is 3.2% of the total area in the late 1960s. To assess the influence of glacier runoff on river runoff, a modified degree–day model including potential clear-sky direct solar radiation has been applied to the glaciated regions of the river basin over the period 1961–2004. It was found that glacier runoff has increased in the last 44 years, especially in the 1990s when a two-thirds increase in river runoff was derived from the increase in glacier runoff caused by loss of ice mass in the entire Tuotuo River basin.