Integration of Tracer Test Data to Refine Geostatistical Hydraulic Conductivity Fields Using Sequential Self-Calibration Method

On the basis of local measurements of hydraulic conductivity,geostatistical methods have been found to be useful in heterogeneity characterization of a hydraulic conductivity field on a regional scale. However,the methods are not suited to directly integrate dynamic production data,such as,hydraulic head and solute concentration,into the study of conductivity distribution. These data,which record the flow and transport processes in the medium,are closely related to the spatial distribution of hydraulic conductivity. In this study,a three-dimensional gradient-based inverse method-the sequential self-calibration (SSC) method-is developed to calibrate a hydraulic conductivity field,initially generated by a geostatistical simulation method,conditioned on tracer test results. The SSC method can honor both local hydraulic conductivity measurements and tracer test data. The mismatch between the simulated hydraulic conductivity field and the reference true one,measured by its mean square error (MSE),is reduced through the SSC conditional study. In comparison with the unconditional results,the SSC conditional study creates the mean breakthrough curve much closer to the reference true curve,and significantly reduces the prediction uncertainty of the solute transport in the observed locations. Further,the reduction of uncertainty is spatially dependent,which indicates that good locations,geological structure,and boundary conditions will affect the efficiency of the SSC study results.     

ODSIM: An Object-Distance Simulation Method for Conditioning Complex Natural Structures

Stochastic simulation of categorical objects is traditionally achieved either with object-based or pixel-based methods. Whereas object-based modeling provides realistic results but raises data conditioning problems, pixel-based modeling provides exact data conditioning but may lose some features of the simulated objects such as connectivity. We suggest a hybrid dual-scale approach to combine both shape realism and strict data conditioning. The procedure combines the distance transform to a skeleton object representing coarse-scale structures, plus a classical pixel-based random field and threshold representing fine-scale features. This object-distance simulation method (ODSIM) uses a perturbed distance to objects and is particularly appropriate for modeling structures related to faults or fractures such as karsts, late dolomitized rocks, and mineralized veins. We demonstrate this method to simulate dolomite geometry and discuss strategies to apply this method more generally to simulate binary shapes.     

Anisotropic Mean Traveltime Curves: A Method to Estimate Anisotropic Parameters from 2D Transmission Tomographic Data

We present the mathematical deduction and properties of the mean traveltime curves for homogeneous elliptical anisotropic media. These curves generalize their isotropic counterparts which have been introduced in the past as a simple data quality analysis technique at the pre-inversion stage for 2D transmission experiments, allowing the inference of prior velocity models to gain stability at the tomographic inversion. Also, the anisotropy parameters (maximum velocity, anisotropic direction and ratio) are shown to affect the shape of these curves. The degree of asymmetry of the anisotropic mean traveltime curves (displacement of the mean time and standard deviation minima from the middle of the gathering line) is related to the direction of anisotropy which can then be visually estimated. Least squares’ fitting of the anisotropic theoretical models to their experimental counterparts is an effective method to estimate at the pre-inversion stage a macroscopic elliptical anisotropic velocity model, valid at the scale of the experiment, and able to match the experimental mean traveltime distribution. Sensitivity analysis has shown that the mean curve is less prone to errors than the standard deviation curve. Parameter identification from the standard deviation curve becomes unstable for noise levels higher than 5%; data errors produce smearing of the value of the estimated anisotropy ratio and wrong directions of anisotropy biased towards zero degrees. Also, identification from the mean traveltime curve becomes stable when the maximum velocity is well constrained. Finally, this methodology is illustrated with the application to the Grimsel data set. Performing MTC analysis is always recommended since it does not need high numerical requirements, and as shown in the sensibility analysis section, errors in data can be misinterpreted as geological anisotropies. J.L. Fernández Martínez is a visiting professor at UC Berkeley, Department of Civil and Environmental Eng., CA 94720-1710.     

An Effective Method to Determine the Distribution Boundary of Cobalt-Rich Fe-Mn Crusts on a Guyot： Synchronous Application of Sub-bottom Profiling and Deep-Sea Video Recording

Cobalt richferromanganesecrust,hereaftercalled Fe Mncrust,isoneoftheimportantmarinemineralre sourcesintheinternationalseabed.Fe Mncrustoccurs onthesurfaceofseamounts,whichareenrichedinco balt,nickel,copper,platinumgroupelements,rare earthelementsandothe…     

A GIS and Remote Sensing-based Analysis of Land Use Change Using the Asymmetric Relation Analysis Method: A Case Study from the City of Hangzhou,China

Urbanization has been an important component of land use and land cover change. Quantifying land use patterns and their changes is essential for monitoring and assessing the impact of urbanization. Data on existing and changing land use provide crucial clues for future development. Combining the application of remote sensing with geographic information systems (GIS) data, this study accomplishes (i) the creation of time-series land use maps in the core urban area of the city of Hangzhou, China; (ii) the examination of the land use cross-tabulation matrix to assess the transitions between land use types; and (iii) the introduction of a new asymmetric relation analysis method to explore the dependency relationship of different times in land use. The investigation was based on a 12-year time series dataset (1990–2002), compiled from fine-resolution satellite imagery and GIS-based maps. The results of the study showed that, first, there are strong associations between different stages (1990, 1996, and 2002) in land use development, and, secondly, the dependency between land use states was asymmetric in both stages of development. In addition, urbanization in Hangzhou has resulted in notable increases in the surface area of urban land and sharp decreases in the surface area of cropland.     

Ferric/ferrous ratio in silicate melts: a new model for 1 atm data with special emphasis on the effects of melt composition

The effect of MgO and total FeO on ferric/ferrous ratio in model multicomponent silicate melts was investigated experimentally in the temperature range 1300–1500 °C at 1 atm total pressure in air. We demonstrate that the addition of these weak network modifier cations results in an increase of Fe³⁺/Fe²⁺ ratio in both mafic and silicic melts. Based on present and published experimental data, a new empirical equation is proposed to predict the ferric/ferrous ratio as a function of oxygen fugacity, temperature and melt composition. In contrast to previous equations, the compositional effect of melts on the Fe³⁺/Fe²⁺ ratio is not only modeled by the sum of the molar fraction of the individual oxide components. Additional interactions terms have also been incorporated. The main advantage of the proposed model is its applicability for a wide compositional range. However, its application to felsic melts (>?68 wt% SiO₂) is not recommended. Other advantages of this equation and differences when compared with previous models are discussed.     

The International Association of Hydrogeologists (IAH): reflecting on 60 years of contributions to groundwater science and water management

The 60th anniversary of the founding of the International Association of Hydrogeologists (IAH) is an important milestone that allows pause for reflection on how the association has evolved over the years and the contributions it has made to groundwater science and water management. IAH was founded in 1956 at the 20th International Geological Congress and developed rapidly during the 1980s and 1990s in response to a growing global interest in groundwater mapping and in sound approaches to resource protection and sustainable aquifer management. Incorporated in 2000, IAH has now secured its position as the world’s leading international association specialising in groundwater with over 4,100 members in 131 countries. Much credit for this success must go to members, past and present, whose individual efforts and collaboration with sister institutions are documented here. These members have shaped the association’s goals and contributed selflessly to its scientific programmes, publications and educational and charitable activities. Looking ahead to the next 60 years, it is essential that IAH does not rest on past achievements but listens and adjusts to the needs of members while continuing to pursue its mission of furthering the understanding, wise use and protection of groundwater resources throughout the world.     

Fluids in the peridotite–water system up to 6 GPa and 800°C: new experimental constrains on dehydration reactions

The phase assemblages and compositions in a K-free lherzolite + H₂O system were determined between 4 and 6 GPa and 700–800°C, and the dehydration reactions occurring at subarc depth in subduction zones were constrained. Experiments were performed on a rocking multi-anvil apparatus using a diamond-trap setting. The composition of the fluid phase was measured using the recently developed cryogenic LA–ICP–MS technique. Results show that, at 4 GPa, the aqueous fluid coexisting with residual lherzolite (~85 wt% H₂O) doubles its solute load when chlorite transforms to the 10-Å phase between 700 and 750°C. The 10-Å phase breaks down at 4 and 5 GPa between 750 and 800°C and at 6 GPa between 700 and 750°C, leaving a dry lherzolite coexisting with a fluid phase containing 58–67 wt% H₂O, again doubling the total dissolved solute load. The fluid fraction in the system increases from 0.2 when a hydrous mineral is present to 0.4 when coexisting with a dry lherzolite. Our data do not reveal the presence of a hydrous peridotite solidus below 800°C. The directly measured fluid compositions demonstrate a fundamental change in the (MgO + FeO) to SiO₂ mass ratio of fluid solutes occurring at a depth of ca. 120–150 km (in the temperature window of 700–800°C), from (MgO–FeO)-dominated at 4 GPa [with (MgO + FeO)/SiO₂ ratio of 1.41–1.56] to SiO₂-dominated at 5–6 GPa (ratios of 0.61–0.82). The mobility of Al₂O₃ increases by more than one order of magnitude across this P–T interval and demonstrates that Al₂O₃ is compatible in an aqueous fluid coexisting with the anhydrous ol-opx-cpx ± grt assemblage. This shift in the fluid composition correlates with changes in the phase assemblage of the residual silicates. The hitherto unknown fundamental change in (MgO + FeO)/SiO₂ ratio and prominent increase in Al₂O₃ of the aqueous fluid with progressive subduction will likely inspire novel concepts on mantle wedge metasomatism by slab fluids.     

Effect of $$\hbox {CO}_2$$ solubility on dissolution rate of calcite in saline aquifers for temperature range of 50–100 $$^{\circ } \hbox {C}$$∘C and pressures up to 600 bar: alterations of fractures geometry in carbonate rocks by $$\hbox {CO}_2$$-acidified brines

An empirical model is developed to predict the dissolution rate of calcite in saline solutions that are saturated with respect to dissolved \(\hbox {CO}_2\) over a broad range of both subcritical and supercritical conditions. The focus is on determining the rate of calcite dissolution within a temperature range of 50–100 \(^\circ \hbox {C}\) and pressures up to 600 bar, relevant for \(\hbox {CO}_2\) sequestration in saline aquifers. A general reaction kinetic model is used that is based on the extension of the standard Arrhenius equation with an added, solubility-dependent, pH term to account for the saturated concentration of dissolved \(\hbox {CO}_2\). On the basis of this kinetic model, a new rate equation is obtained using multi-parameter, nonlinear regression of experimental data to determine the dissolution of calcite as a function of temperature, pressure and salinity. Different models for the activity coefficient of \(\hbox {CO}_2\) dissolved in saline solutions are accounted for. The new rate equation helps us obtain good agreement with experimental data, and it is applied to study the geochemically induced alterations of fracture geometry due to calcite dissolution.