In situ estimation of the effective chemical diffusion coefficient of a rock matrix in a fractured aquifer

An in situ method of estimating the effective diffusion coefficient for a chemical constituent that diffuses into the primary porosity of a rock is developed by abruptly changing the concentration of the dissolved constituent in a borehole in contact with the rock matrix and monitoring the time-varying concentration. The experiment was conducted in a borehole completed in mudstone on the campus of the University of the Free State in Bloemfontein, South Africa. Numerous tracer tests were conducted at this site, which left a residual concentration of sodium chloride in boreholes that diffused into the rock matrix over a period of years. Fresh water was introduced into a borehole in contact with the mudstone, and the time-varying increase of chloride was observed by monitoring the electrical conductivity (EC) at various depths in the borehole. Estimates of the effective diffusion coefficient were obtained by interpreting measurements of EC over 34 d. The effective diffusion coefficient at a depth of 36 m was approximately 7.8×10^?6 m²/d, but was sensitive to the assumed matrix porosity. The formation factor and mass flux for the mudstone were also estimated from the experiment.     

A continuous/discontinuous Galerkin framework for modeling coupled subsurface and surface water flow

We consider conjunctive surface-subsurface flow modeling, where surface water flow is described by the shallow water equations and ground water flow by Richards’ equation for the vadose zone. Coupling between the models is based on the continuity of flux and water pressure. Numerical approximation of the coupled model using the framework of discontinuous Galerkin (DG) methods is formulated. In the subsurface, the local discontinuous Galerkin (LDG) method is used to approximate ground water velocity and hydraulic head; a DG method is also used to approximate surface water velocity and elevation. This approach allows for a weak coupling of the models and the use of different approximating spaces and/or meshes within each regime. A simplified LDG method based on continuous approximations to water head is also described. Numerical results that investigate physical and numerical aspects of surface–subsurface flow modeling are presented. This work was supported by National Science Foundation grant DMS-0411413.     

Hydrogeochemistry and isotope analyses used to determine groundwater recharge and flow in low-gradient catchments of the province of Buenos Aires,Argentina

The region known as Pampa Plain, in Argentina, is a vast area characterized by slopes of less than 0.05%. The surface sediment is silty sand, mainly Aeolian, referred to as Pampean Loess, which is a phreatic aquifer unit of utmost importance for the water supply of the region. On account of the slight gradient, hydrogeological analyses using only hydraulic measurements are difficult to perform, often leading to confusing results. Thus, the study presented relies on hydrochemical modeling and isotopic determinations as well. The study area comprises three catchments in the inter-mountainous area corresponding to El Moro, Tamangueyú and Seco creeks, and covering 2,570 km² in the province of Buenos Aires. Measurements of piezometric levels and samples of groundwater, surface water and rainwater were carried out between January and March 2005. Major ion results were analyzed by means of hydrochemical graphs and hydrogeochemical modelling using NETPATH. The resulting data show that it is possible to identify local changes in recharge, flow direction and stream/groundwater relationship by using hydrochemical and isotopic information, which may become a useful and more precise tool for the study of particularly flat landscapes.     

High heat flow anomalies on an old oceanic plate observed seaward of the Japan Trench

Thirty-three new measurements on the seaward slope and outer rise of the Japan Trench along a parallel of 38°45′N revealed the existence of high heat flow anomalies on the subducting Pacific plate, where the seafloor age is about 135 m.y.. The most prominent anomaly with the highest value of 114 mW/m² is associated with a small mound on the outer rise, which was reported to be a kind of mud volcano. On the seaward slope of the trench, heat flow is variable: high (70–90 mW/m²) at some locations and normal for the seafloor age (about 50 mW/m²) at others. The spatial variation of heat flow may be related to development of normal faults and horst/graben structures due to bending of the Pacific plate before subduction, with fluid flow along the fault zones enhancing the vertical heat transfer. Possible heat sources of the high heat flow anomalies are intra-plate volcanism in the last several million years like that discovered recently on the Pacific plate east of the Japan Trench.     

Estimates of heat flow and heat production and a thermal model of the São Francisco craton

An updated analysis of geothermal data from the highland area of eastern Brazil has been carried out and the characteristics of regional variations in geothermal gradients and heat flow examined. The database employed includes results of geothermal measurements at 45 localities. The results indicate that the Salvador craton and the adjacent metamorphic fold belts northeastern parts of the study area are characterized by geothermal gradients in the range of 6–17°C/km. The estimated heat flow values fall in the range of 28–53 mW/m², with low values in the cratonic area relative to the fold belts. On the other hand, the São Francisco craton and the intracratonic São Francisco sedimentary basin in the southwestern parts are characterized by relatively higher gradient values, in the range of 14–42°C/km, with the corresponding heat flow values falling in the range of 36–89 mW/m². Maps of regional variations indicate that high heat flow anomaly in the São Francisco craton is limited to areas of sedimentary cover, to the west of the Espinhaço mountain belt. Crustal thermal models have been developed to examine the implications of the observed intracratonic variations in heat flow. The thermal models take into consideration variation of thermal conductivity with temperature as well as change of radiogenic heat generation with depth. Vertical distributions of seismic velocities were used in obtaining estimates of radiogenic heat production in crustal layers. Crustal temperatures are calculated based on a procedure that makes simultaneous use of the Kirchoff and Generalized Integral Transforms, providing thereby analytical solutions in 2D and 3D geometry. The results point to temperature variations of up to 300°C at the Moho depth, between the northern Salvador and southern São Francisco cratons. There are indications that differences in rheological properties, related to thermal field, are responsible for the contrasting styles of deformation patterns in the adjacent metamorphic fold belts.     

Application of FLATModel, a 2D finite volume code, to debris flows in the northeastern part of the Iberian Peninsula

FLATModel is a two-dimensional shallow-water approximation code with corrections and modifications that create a simulation tool adapted to debris-flows behaviour. FLATModel uses the finite volume method with the numerical implementation of the Godunov scheme and includes correction terms regarding the effect of flow over high slopes and curvature. Additionally, the stop-and-go phenomenon, the basal entrainment and a correction regarding the front inclination of the final deposit are incorporated into FLATModel. In addition, different flow resistance laws were integrated in the numerical code including Bingham, Herschel–Bulkley and Voellmy fluid model. Firstly, our numerical model was validated using analytical solutions of a dam-break scenario and published data on a laboratory experiment. Secondly, three real events, which occurred in the northeastern part of the Iberian Peninsula, were back-calculated. Although field observations of the three events are not very detailed, the back-analyses revealed interesting patterns on the flow dynamics, and the numerical results generally showed good agreement with field data. Comparing the different flow resistance laws, the Voellmy fluid model presents the best behaviour regarding both the flow behaviour and the deposit characteristics. Preliminary simulation runs incorporating the effect of basal entrainment offered satisfactory results, although the final volume is rather sensitive on the selected friction angle of channel-bed material. The outcomes regarding the correction of the calculated front inclination of the final deposit showed that this implementation strongly improves the simulation results and better represents steep fronts of final deposits.     

Inversely-Mapped Analytical Solutions for Flow Patterns around and within Inclined Elliptic Inclusions in Fluid-Saturated Rocks

In this paper, an inverse mapping is used to transform the previously-derived analytical solutions from a local elliptical coordinate system into a conventional Cartesian coordinate system. This enables a complete set of exact analytical solutions to be derived rigorously for the pore-fluid velocity, stream function, and excess pore-fluid pressure around and within buried inclined elliptic inclusions in pore-fluid-saturated porous rocks. To maximize the application range of the derived analytical solutions, the focal distance of an ellipse is used to represent the size of the ellipse, while the length ratio of the long axis to the short one is used to represent the geometrical shape of the ellipse. Since the present analytical solutions are expressed in a conventional Cartesian coordinate system, it is convenient to investigate, both qualitatively and quantitatively, the distribution patterns of the pore-fluid flow and excess pressure around and within many different families of buried inclined elliptic inclusions. The major advantage in using the present analytical solution is that they can be conveniently computed in a global Cartesian coordinate system, which is widely used in many scientific and engineering computations. As an application example, the present analytical solutions have been used to investigate how the dip angle of an inclined elliptic inclusion affects the distribution patterns of the pore-fluid flow and excess pore-fluid pressure when the permeability ratio of the elliptic inclusion is of finite but nonzero values.     

Micromechanical approach to unsaturated water flow in structured geomaterials by two-scale computations

This article presents a micromechanical approach to the problem of unsaturated water flow in heterogeneous porous media in transient conditions. The numerical formulation is based on the two-scale model obtained previously by periodic homogenization. It allows for a coupled solution of the non-linear flow equations at macroscopic and microscopic scales and takes into account the macroscopic anisotropy of the medium and the local non-equilibrium of the capillary pressure. The model was applied to simulate two-dimensional water infiltration at constant flux into an initially dry medium containing inclusions of square and rectangular shapes. The numerical results showed the influence of the inclusion–matrix conductivity ratio and the local geometry on the macroscopic behavior. The influence of the conductivity ratio manifested itself by the acceleration or retardation of the onset of the macroscopic water flux at the outlet, while the local geometry (anisotropy) significantly affected the macroscopic spatial distribution of the water flux. Such type of approach can be extended to simulate coupled phenomena (for example hydro-mechanical problems) with evolving local geometry.     

The dynamics and thermodynamics of large ash flows

 Ash flow deposits, containing up to 1000 km³ of material, have been produced by some of the largest volcanic eruptions known. Ash flows propagate several tens of kilometres from their source vents, produce extensive blankets of ash and are able to surmount topographic barriers hundreds of metres high. We present and test a new model of the motion of such flows as they propagate over a near horizontal surface from a collapsing fountain above a volcanic vent. The model predicts that for a given eruption rate, either a slow (10–100 m/s) and deep (1000–3000 m) subcritical flow or a fast (100–200 m/s) and shallow (500–1000 m) supercritical flow may develop. Subcritical ash flows propagate with a nearly constant volume flux, whereas supercritical flows entrain air and become progressively more voluminous. The run-out distance of such ash flows is controlled largely by the mass of air mixed into the collapsing fountain, the degree of fragmentation and the associated rate of loss of material into an underlying concentrated depositional system, and the mass eruption rate. However, in supercritical flows, the continued entrainment of air exerts a further important control on the flow evolution. Model predictions show that the run-out distance decreases with the mass of air entrained into the flow. Also, the mass of ash which may ascend from the flow into a buoyant coignimbrite cloud increases as more air is entrained into the flow. As a result, supercritical ash flows typically have shorter runout distances and more ash is elutriated into the associated coignimbrite eruption columns. We also show that one-dimensional, channellized ash flows typically propagate further than their radially spreading counterparts. As a Plinian eruption proceeds, the erupted mass flux often increases, leading to column collapse and the formation of pumiceous ash flows. Near the critical conditions for eruption column collapse, the flows are shed from high fountains which entrain large quantities of air per unit mass. Our model suggests that this will lead to relatively short ash flows with much of the erupted material being elutriated into the coignimbrite column. However, if the mass flux subseqently increases, then less air per unit mass is entrained into the collapsing fountain, and progressively larger flows, which propagate further from the vent, will develop. Our model is consistent with observations of a number of pyroclastic flow deposits, including the 1912 eruption of Katmai and the 1991 eruption of Pinatubo. The model suggests that many extensive flow sheets were emplaced from eruptions with mass fluxes of 10⁹–10¹⁰ kg/s over periods of 10³–10⁵ s, and that some indicators of flow "mobility" may need to be reinterpreted. Furthermore, in accordance with observations, the model predicts that the coignimbrite eruption columns produced from such ash flows rose between 20 and 40 km. Received: 25 August 1995 / Accepted: 3 April 1996     

地下水非稳定流的LTFAM算法

渗流域内应用拉普拉斯变换(LT)建立相应的有限分析(FAM)方程,顾及渗流域内地下水流的初始条件和边界条件,可在LT空间构成一个封闭的以水头像函数为变量的线性方程组。将此方程组所得的解,通过Stehfest数值反演公式,可归化为时间域的解(水头)。由于时间t被隐含在数值方程内,从而克服了传统数值法按时段(△t)逐步迭代的缺陷,提高了计算效率,也为用嵌入法建立地下水流管理模型提供了一条捷径。