Formation compaction associated with thermal cooling in geothermal reservoirs

Concepts from the Theory of Interacting Continua are employed to develop constitutive relations for thermoelastic fluid-saturated porous media; these constitutive relations are shown to be equivalent to those of Biot in isothermal limit. For uniaxial compaction, the constitutive theory is simplified to give formation compaction as a function of changes in fluid pressure and temperature. The formation compaction due to thermal cooling is equal to 3hη[(1 + v)/3(1 − v)]ΔT.     

δCa evolution in a carbonate aquifer and its bearing on the equilibrium isotope fractionation factor for calcite

To improve understanding of Ca isotope transport during water-rock interaction on the continents, we measured dissolved δ⁴⁴Ca values along a 236 km flow path in the Madison aquifer, South Dakota, where fluids have chemically evolved according to dolomite and anhydrite dissolution, calcite precipitation, and Ca-for-Na ion-exchange over a timescale spanning ~ 15 kyr. We used a reactive transport model employing rate data constrained from major ion mass-balances to evaluate the extent to which calcite precipitation and ion-exchange fractionate Ca isotopes. Elevated δ⁴⁴Ca values during the initial and final stages of water transport possibly result from calcite precipitation under supersaturated conditions and Ca-for-Na ion-exchange, respectively. However, for the bulk of the flow path, δ⁴⁴Ca values evolve by mixing between anhydrite and dolomite dissolution, with no fractionation during calcite precipitation under saturated conditions. We attribute the absence of Ca isotope fractionation to the long timescale of water-rock interaction and slow rate of calcite precipitation, which have enabled fluids to chemically and isotopically equilibrate with calcite. We therefore conclude that the equilibrium Ca isotope fractionation factor between calcite and water (Δ_cal–w) is very close to zero. To the extent that the Madison aquifer typifies other groundwater systems where calcite slowly precipitates from solutions at or near chemical equilibrium, this study suggests that groundwater contributions to δ⁴⁴Ca variability on the continents can be modeled according to simple mixing theory without invoking isotope discrimination.     

Estimating the mean speed of laminar overland flow using dye injection‐uncertainty on rough surfaces

A common method for estimating mean flow speeds in studies of surface runoff is to time the travel of a dye cloud across a measured flow path. Motion of the dye front reflects the surface flow speed, and a correction must be employed to derive a value for the profile mean speed, which is always lower. Whilst laminar flow conditions are widespread in the interrill zone, few data are available with which to establish the relationship linking surface and profile mean speeds, and there are virtually none for the flow range 100 < Re < 500 (Re = Reynolds number) which is studied here. In laboratory experiments on a glued sand board, mean flow speeds were estimated from both dye speeds and the volumetric flow relation v = Q/ wd with d measured using a computer‐controlled needle gauge at 64 points. In order to simulate conditions applicable to many dryland soils, the board was also roughened with plant litter and with ceramic tiles (to simulate surface stone cover). Results demonstrate that in the range 100 < Re < 500, there is no consistent relation between surface flow speeds and the profile mean. The mean relationship is v = 0·56 v _surf, which departs significantly from the theoretical smooth‐surface relation v = 0·67 v _surf, and exhibits a considerable scatter of values that show a dependence on flow depth. Given the inapplicability of any fixed conversion factor, and the dependence on flow depth, it is suggested that the use of dye timing as a method for estimating v be abandoned in favour of precision depth measurement and the use of the relation v = Q/ wd , at least within the laminar flow range tested. Copyright © 2001 John Wiley & Sons, Ltd.     

Thermal nature of the source of the Indian gravity anomaly

The non-hydrostatic part of the second zonal harmonic term has been included in the Indian gravity anomaly. The interpretation of the potential, gravity and deflection of the vertical has been produced by a new method. The depth of the source is 580 km and the solution satisfies the temperature gradient of 1°K/km. The ratio Δ V_p/Δp > 0 shows the thermal nature of the density distribution near 670 km.     

Noble gas tracers for characterisation of flow dynamics and origin of groundwater: A case study in Switzerland

The Grenchen aquifer system in the Swiss Plateau was extensively investigated in order to determine the extent of groundwater contamination and to assess the natural attenuation capacity. Environmental tracer data were applied to estimate groundwater travel times, mixing ratios, and evaluate groundwater origin. Recharge is basically possible in two distinct topographical areas, the immediate vicinity of the town of Grenchen and the elevated plateau of the first Jura Mountain ridge. Groundwater dating was performed with the ³H/³He dating method and supplemented by ⁸⁵Kr measurements. Stable isotope data (δ¹⁸O, δ²H) and dissolved noble gas concentrations allow the determination of the recharge temperature, which is correlated to the recharge elevation. Noble gas temperatures (NGT) decrease in the direction of groundwater flow and range from 10 to 13 °C in the upstream area of the town to 7–9 °C in the downstream river plain. This trend could suggest the admixture of water from the underlying limestone aquifer recharged under cooler infiltration conditions, e.g. at higher recharge elevations. However, it is shown in this study that the difference in NGT does not require such a recharge. Rather, increasing air temperatures over the last 40 years and the urban heat island effect could possibly explain most of the observed temperature shift. Furthermore, it is concluded that the downstream river plain is hydrologically disconnected from the upstream town area. Consequently most water from the town area is drained by the creek Witibach and recharge in the river plain is higher than previously assumed.     

A combined study of macroseismic data and focal mechanisms applied to the West-Bohemian earthquake,Czechoslovakia

The main shock of the West-Bohemian earthquake swarm, Czechoslovakia, (magnitudem=4.5, depthh=10 km) exhibits an irregular areal distribution of macroseismic intensities 6° to 7° MSK-64. Four lobes of the 6° isoseismal are found and the maximum observed intensity is located at a distance of 8 km from the instrumentally determined epicentre. This distribution can be explained by the energy flux of the directS wave generated by a circular source, the hypocentral location and focal mechanism of which are taken from independent instrumental studies. The theoretical intensity, which is assumed to be logarithmically proportional to the integrated squared ground-motion velocity (i.e.,I=const+log v ² (t)dt), fits the observed intensity with an overall root-mean-square error less than 0.5°. It is important that the present intensity data can also be equally well explained by the isotropic source. The fit was attained by means of a horizontally layered model though large fault zones and an extended sedimentary basin suggest a significant lateral heterogeneity of the epicentral region. The results encourage a broader application of the simple modelling technique used.     

Eddy energetics in the Faroe-Shetland channel: a model resolution study

Observations indicate that, perhaps as a result of intense shears produced by the topographic steering of the slope current and merging of the warm North Atlantic water with cold Nordic waters, the Faroe-Shetland Channel is a fertile ground for energetic eddies and meanders. How these meso-scale ( 10 km) features are resolved in a three-dimensional numerical model is studied here with grid-nesting, in which the convergence of the simulated energetics is examined as a function of grid sizes Δ =10, 10/3, 10/5 and 10/7 km. Large-scale forcing to the fine-grid nest consists of steady slope-current transport (2 Sv) and north/south climatological density difference only. The mean fields were found to converge at Δ = 10/3 km, but the eddy fields required Δ = 2 km or less. At the finest grids, meso-scale features more in conformity with those observed were found: meanders with amplitudes of O(10 km) and along-slope propagation speeds of about 0.14 m s- 1 cross-channel vacillation of the slope current, with amplitude of about 10 km and periods of 5–10 days, and vertical excursions of isotherms of O(100 m). The meanders were highly ageostrophic (ζ/f> 1, ζ=relative vorticity, f = Coriolis parameter), and developed as a result of conditions favorable to flow instability originated from south of the channel. When adequately resolved, the meander energetics are comparable to those produced by (large-scale) winds, which suggests that a high-resolution grid (Δ < 2 km) will be required for hindcast and forecast studies.     

Groundwater recharge estimation in arid hardrock‐alluvium aquifers using combined water‐table fluctuation and groundwater balance approaches

This paper proposes an approach to estimate groundwater recharge using an optimization‐based water‐table fluctuation method combined with a groundwater balance model in an arid hardrock‐alluvium region, located at the Oman–United Arab Emirates border. We introduce an “effective hardrock thickness” term to identify the percentage of the considered hardrock thickness in which effective groundwater flow takes place. The proposed method is based upon a Thiessen polygon zoning approach. The method includes subpolygons to represent specific geologic units and to enhance the confidence of the estimated groundwater recharge. Two linear and 1 nonlinear submodels were developed to evaluate the model components for the calibration (October 1996 to September 2008) and validation (October 2008 to September 2013) periods. Long‐term annual groundwater recharge from rainfall and return flow over the model domain are estimated as 24.62 and 5.71 Mm³, respectively, while the effective groundwater flow circulation is found to occur in the upper 7% of the known hardrock thickness (42 m), confirming conclusions of previous field studies. Considering a total difference in groundwater levels between eastern and western points of the study area of the order of 220 m and a 12‐year monthly calibration period, a weighted root mean squared error in predicted groundwater elevation of 2.75 m is considered quite reasonable for the study area characterized by remarkable geological and hydrogeological diversity. The proposed approach provides an efficient and robust method to estimate groundwater recharge in regions with a complex geological setting in which interaction between fractured and porous media cannot be easily assessed.     

Simulation of the formation of Ushant thermal front

The formation of the Ushant thermal front off northwestern France is simulated numerically with a one-dimensional model of the vertical thermal structure of the sea. This simulation is derived from a Niiler and Kraus (1977, Modelling and prediction of the upper layers of the ocean, pp. 143–172) mixed layer model in which bottom friction is introduced. The model inputs are the meteorological parameters measured at Ushant and the barotropic tidal currents computed with a numerical model (Mariette et al., 1982, Oceanologica Acta, 5, 149–159). Both very fine time steps (Δt = 1h) and spatial steps (horizontal mesh size: Δx = Δy = 2nmi; vertical step: Δz = 1m) are used. This study is intended to be mainly concerned with the processes acting during the formation of the Ushant front. Comparison of the model results with satellite pictures and in situ measurements shows good agreement. We point out that biological blooms could occur according to a spring-neap cycle, during the thermal front formation. The prevailing effect of local processes such as tidally induced bottom friction and energy budget at the sea surface on the thermal evolution of the Mer d'Iroise is then emphasized. During the formation of the front, these effects seem to be more important than the advective processes which are not considered in this work.     

Higher and lowest order mixed finite element approximation of subsurface flow problems with solutions of low regularity

In this work we study mixed finite element approximations of Richards’ equation for simulating variably saturated subsurface flow and simultaneous reactive solute transport. Whereas higher order schemes have proved their ability to approximate reliably reactive solute transport (cf., e.g. [Bause M, Knabner P. Numerical simulation of contaminant biodegradation by higher order methods and adaptive time stepping. Comput Visual Sci 7;2004:61–78]), the Raviart–Thomas mixed finite element method (RT₀) with a first order accurate flux approximation is popular for computing the underlying water flow field (cf. [Bause M, Knabner P. Computation of variably saturated subsurface flow by adaptive mixed hybrid finite element methods. Adv Water Resour 27;2004:565–581, Farthing MW, Kees CE, Miller CT. Mixed finite element methods and higher order temporal approximations for variably saturated groundwater flow. Adv Water Resour 26;2003:373–394, Starke G. Least-squares mixed finite element solution of variably saturated subsurface flow problems. SIAM J Sci Comput 21;2000:1869–1885, Younes A, Mosé R, Ackerer P, Chavent G. A new formulation of the mixed finite element method for solving elliptic and parabolic PDE with triangular elements. J Comp Phys 149;1999:148–167, Woodward CS, Dawson CN. Analysis of expanded mixed finite element methods for a nonlinear parabolic equation modeling flow into variably saturated porous media. SIAM J Numer Anal 37;2000:701–724]). This combination might be non-optimal. Higher order techniques could increase the accuracy of the flow field calculation and thereby improve the prediction of the solute transport. Here, we analyse the application of the Brezzi-Douglas-Marini element (BDM₁) with a second order accurate flux approximation to elliptic, parabolic and degenerate problems whose solutions lack the regularity that is assumed in optimal order error analyses. For the flow field calculation a superiority of the BDM₁ approach to the RT₀ one is observed, which however is less significant for the accompanying solute transport.     

Investigating soil moisture feedbacks on precipitation with tests of Granger causality

Granger causality (GC) is used in the econometrics literature to identify the presence of one- and two-way coupling between terms in noisy multivariate dynamical systems. Here we test for the presence of GC to identify a soil moisture (S) feedback on precipitation (P) using data from Illinois. In this framework S is said to Granger cause P if F(P_t|Ω_t−Δt)≠F(P_t|Ω_t−Δt−S_t−Δt) where F denotes the conditional distribution of P, Ω_t−Δt represents the set of all knowledge available at time t−Δt, and Ω_t−Δt−S_t−Δt represents all knowledge except S. Critical for land–atmosphere interaction research is that Ω_t−Δt includes all past information on P as well as S. Therefore that part of the relation between past soil moisture and current precipitation which results from precipitation autocorrelation and soil water balance will be accounted for and not attributed to causality. Tests for GC usually specify all relevant variables in a coupled vector autoregressive (VAR) model and then calculate the significance level of decreased predictability as various coupling coefficients are omitted. But because the data (daily precipitation and soil moisture) are distinctly non-Gaussian, we avoid using a VAR and instead express the daily precipitation events as a Markov model. We then test whether the probability of storm occurrence, conditioned on past information on precipitation, changes with information on soil moisture. Past information on precipitation is expressed both as the occurrence of previous day precipitation (to account for storm-scale persistence) and as a simple soil moisture-like precipitation-wetness index derived solely from precipitation (to account for seasonal-scale persistence). In this way only those fluctuations in moisture not attributable to past fluctuations in precipitation (e.g., those due to temperature) can influence the outcome of the test. The null hypothesis (no moisture influence) is evaluated by comparing observed changes in storm probability to Monte-Carlo simulated differences generated with unconditional occurrence probabilities. The null hypothesis is not rejected (p>0.5) suggesting that contrary to recently published results, insufficient evidence exists to support an influence of soil moisture on precipitation in Illinois.     

Variability in the vertical hyporheic water exchange affected by hydraulic conductivity and river morphology at a natural confluent meander bend

River confluences and their associated tributaries are key morphodynamic nodes that play important roles in controlling hydraulic geometry and hyporheic water exchange in fluvial networks. However, the existing knowledge regarding hyporheic water exchange associated with river confluence morphology is relatively scarce. On January 14 and 15, 2016, the general hydraulic and morphological characteristics of the confluent meander bend (CMB) between the Juehe River and the Haohe River in the southern region of Xi'an City, Shaanxi Province, China, were investigated. The patterns and magnitudes of vertical hyporheic water exchange (VHWE) were estimated based on a one‐dimensional heat steady‐state model, whereas the sediment vertical hydraulic conductivity (K_v) was calculated via in situ permeameter tests. The results demonstrated that 6 hydrodynamic zones and their extensions were observed at the CMB during the test period. These zones were likely controlled by the obtuse junction angle and low momentum flux ratio, influencing the sediment grain size distribution of the CMB. The VHWE patterns at the test site during the test period mostly showed upwelling flow dominated by regional groundwater discharging into the river. The occurrence of longitudinal downwelling and upwelling patterns along the meander bend at the CMB was likely subjected to the comprehensive influences of the local sinuosity of the meander bend and regional groundwater discharge and finally formed regional and local flow paths. Additionally, in dominated upwelling areas, the change in VHWE magnitudes was nearly consistent with that in K_v values, and higher values of both variables generally occurred in erosional zones near the thalweg paths of the CMB, which were mostly made up of sand and gravel. This was potentially caused by the erosional and depositional processes subjected to confluence morphology. Furthermore, lower K_v values observed in downwelling areas at the CMB were attributed to sediment clogging caused by local downwelling flow. The confluence morphology and sediment K_v are thus likely the driving factors that cause local variations in the VHWE of fluvial systems.     

Collision frequencies in the D-region

Partial reflection differential absorption and differential phase observations have been used to systematically study collision frequencies (ν) in D-region of the lower ionosphere. The observations made with the large MF radar located at Buckland Park (35°S, 138°E) near Adelaide in the period September 1996–December 2000 show ν values larger than predicted by previously used models. The new estimates are compared with values calculated using new collision frequency momentum cross-sections for N₂(σ) measured in the laboratory. The two types of measurements are found to be in good agreement, with a moderate seasonal variation of ν at a constant height. No change due to solar cycle variations is found.     

The effect of expansion ratio on the critical Reynolds number in single fracture flow with sudden expansion

Flow in a single fracture (SF) is an important research subject in groundwater hydrology, hydraulic engineering, radioactive nuclear waste repository and geotechnical engineering. An abruptly changing aperture is a unique type of SF. This study discusses the relation between the values of the critical Reynolds number (Re_c) for the onset of symmetry breaking of flow and the expansion ratio (E) of SF, which is defined as the ratio between the outlet (D) and inlet (d) apertures. This study also investigates the effect of inlet aperture d on Re_c for flow in an SF with abruptly changing apertures (SF‐ACA) using the finite volume method. Earlier numerical and experimental results showed that flow is symmetric in respect to the central plane of the SF‐ACA at small Reynolds number (Re) but becomes asymmetric when Re is sufficiently large. Our simulations show that the value of Re_c decreases with the increasing E, and the relationship between the logarithm of Re_c and E can be described accurately using either a quadratic polynomial function or a logarithmic function. However, the relationship of Re_c and d for a given E value is vague, and Re_c becomes even less sensitive to d when E increases. This study also reveals that the hydraulic gradient (J) and flow velocity (v) follow a super‐linear relationship that can be fitted almost perfectly by the Forchheimer equation. The inertial component (J_i) of J increases monotonically with Re, whereas the viscous component (J_v) of J decreases monotonically with Re. The Re value corresponding to equal inertial and viscous components of J (named as the transitional point Re) decreases when E increases, and such a transitional point Re should be closely related to the critical Reynolds number Re_c, although a rigorous theoretical proof is not yet available. Copyright © 2015 John Wiley & Sons, Ltd.     

Ionospheric index of solar activity based on the data of measurements of the spacecraft signals characteristics

The method for estimating an ionospheric index of solar-activity (IISA) based on the processing of spacecraft radio signals is suggested. The IISA values have been obtained by comparison between the measured and calculated variations of radio-signal characteristics. To calculate the variations of radio-signal characteristics, the straight rays approximation and the solar-activity index (Wolf numbers W and/or values of F_10.7 solar flux) as a control parameter of the ionospheric model have been used. The suggested method was tested using spacecraft radio signals from the radio-navigation system “CIKADA”. The reduced differences of phases (ΔΦ) for frequencies 150 and 400 MHz were measured and the same characteristics were calculated by integration along the ray of radio-wave propagation between the receiver and the satellite. The IRI-95 has been used as a background ionospheric model. The satellite co-ordinates were determined using the orbit parameters recorded in the navigation messages. Minimization of the difference measured and calculated ΔΦ using arbitrary time steps, or during whole time intervals of observation, gives the IISA corresponding the satellite pass. Daily IISA values were obtained by averaging over all communication contacts during a day (20–30 passes). Testing this approach based on the measurements during March/April 1997, 1998, shows that on magnetically quiet days differences between IISA and the primary solar activity indices are about 5%.     

Salinograph trends as indicators of the recession characteristics of stream components

Compared to hydrograph recession analysis, which is widely applied in engineering hydrology, the quantitative assessment of stream salinity with time (i.e. the salinograph) has received significantly less attention. In particular, while in many previous hydrological studies an inverse relationship between hydrograph and salinograph responses is apparent, the concept of salinity accession (the inversely related salinity counterpart to hydrograph recession) has not been introduced nor quantitatively evaluated in previous literature. In this study, we conduct a mathematical analysis of salinograph accession, and determine new quantitative relationships between salinity accession and hydrograph recession parameters. An equation is formulated that reproduces the general trend in salinity accession. A salinity accession parameter k_c is then introduced and is shown to be the ratio of direct runoff to total stream flow recession parameters: k_r/k. The groundwater recession parameter k_g was estimated using a simple and rapid method that uses both salinograph and hydrograph data. Salinity accession type‐curves illustrate that under certain conditions, the relative steepness of individual salinographs is dependent upon the ratio of groundwater salinity to direct runoff salinity: C_g/C_r. The salinity accession algorithms are applied to two contrasting field settings: Scott Creek, South Australia and Sandy Creek, northern Queensland, Australia. It was found that k_g > k during periods of obvious stream flow recession, for the events analysed. Salinograph accession behaviour was fairly similar for both sites, despite contrasting environments. Using assumed end‐member salinities for groundwater and direct runoff based upon field observations, the behaviour of k_c from the Scott Creek site was approximately reproduced by varying the initial groundwater to runoff flow ratio: Q_g0/Q_r0, within reasonable parameter ranges. The use of salinograph information when used in addition to standard hydrograph analyses provided useful information on recession characteristics of stream components. Copyright © 2008 John Wiley & Sons, Ltd.     

A stochastic approach to nonlinear unconfined flow subject to multiple random fields

In this study, the KLME approach, a moment-equation approach based on the Karhunen–Loeve decomposition developed by Zhang and Lu (Comput Phys 194(2):773–794, 2004), is applied to unconfined flow with multiple random inputs. The log-transformed hydraulic conductivity F, the recharge R, the Dirichlet boundary condition H, and the Neumann boundary condition Q are assumed to be Gaussian random fields with known means and covariance functions. The F, R, H and Q are first decomposed into finite series in terms of Gaussian standard random variables by the Karhunen–Loeve expansion. The hydraulic head h is then represented by a perturbation expansion, and each term in the perturbation expansion is written as the products of unknown coefficients and Gaussian standard random variables obtained from the Karhunen–Loeve expansions. A series of deterministic partial differential equations are derived from the stochastic partial differential equations. The resulting equations for uncorrelated and perfectly correlated cases are developed. The equations can be solved sequentially from low to high order by the finite element method. We examine the accuracy of the KLME approach for the groundwater flow subject to uncorrelated or perfectly correlated random inputs and study the capability of the KLME method for predicting the head variance in the presence of various spatially variable parameters. It is shown that the proposed numerical model gives accurate results at a much smaller computational cost than the Monte Carlo simulation.     

Performance of the White method for estimating groundwater evapotranspiration under conditions of deep and fluctuating groundwater

The White method is a simple but the most frequently applied approach to estimate groundwater evapotranspiration (ET_g) from groundwater level diurnal signals. Because of a lack of direct measurements of ET_g, it is difficult to evaluate the performance of the White method, particularly in field environments with variable groundwater fluctuations. A 2‐year field observation in a hyper‐arid riparian tamarisk (Tamarix spp.) stand with deep groundwater depth in the lower Tarim River basin of China was conducted to measure the surface evapotranspiration (ET_s) and groundwater table. The performance of the White method and the influences of the variable groundwater fluctuations on the determinations of the specific yield (S_y) and recharge rate of groundwater (r) in the White method were investigated. The results showed that the readily available S_y determined by Loheide's method was feasible but must be finely determined based on the soil textures in the layers in which the groundwater level fluctuated. A constant S_y value for a defined porous medium could be assumed regardless of the discharge or recharge processes of groundwater. The time span of 0000 h to 0600 h for r estimation for the White method worked best than other time spans. A 2‐day moving average of r values further enhance ET_g estimation. Slight effects of environmental or anthropogenic disturbances on the diurnal fluctuations of groundwater level did not influence the ET_g estimations by the White method. Our results provide valuable references to the application of the White method for estimating daily ET_g in desert riparian forests with deep groundwater depth. Copyright © 2015 John Wiley & Sons, Ltd.     

Geochemistry of the acid Kawah Putih lake, Patuha Volcano, West Java, Indonesia

Kawah Putih is a summit crater of Patuha volcano, West Java, Indonesia, which contains a shallow, 300 m-wide lake with strongly mineralized acid–sulfate–chloride water. The lake water has a temperature of 26–34°C, pH=<0.5–1.3, S^tot=2500–4600 ppm and Cl=5300–12 600 ppm, and floating sulfur globules with sulfide inclusions are common. Sulfur oxyanion concentrations are unusually high, with S₄O₆²⁻+S₅O₆²⁻+S₆O₆²⁻=2400 – 4200 ppm. Subaerial fumaroles (<93°C) on the lake shore have low molar SO₂/H₂S ratios (<2), which is a favorable condition to produce the observed distribution of sulfur oxyanion species. Sulfur isotope data of dissolved sulfate and native sulfur show a significant ³⁴S fractionation (ΔSO₄–S_e of 20‰), probably the result of SO₂ disproportionation in or below the lake. The lake waters show strong enrichments in ¹⁸O and D relative to local meteoric waters, a result of the combined effects of mixing between isotopically heavy fluids of deep origin and meteoric water, and evaporation-induced fractionation at the lake surface. The stable-isotope systematics combined with energy-balance considerations support very rapid fluid cycling through the lake system. Lake levels and element concentrations show strong seasonal fluctuations, indicative of a short water residence time in the lake as well.Thermodynamic modeling of the lake fluids indicates that the lake water is saturated with silica phases, barite, pyrite and various Pb, Sb, Cu, As, Bi-bearing sulfides when sulfur saturation is assumed. Precipitating phases predicted by the model calculations are consistent with the bulk chemistry of the sulfur-rich bottom sediments and their identified mineral phases. Much of the lake water chemistry can be explained by congruent rock dissolution in combination with preferential enrichments from entering fumarolic gases or brines and element removal by precipitating mineral phases, as indicated by a comparison of the fluids, volcanic rocks and lake bed sediment.Flank springs on the mountain at different elevations vary in composition, and are consistent with local rock dissolution as a dominant factor and pH-dependent element mobility. Discharges of warm sulfate- and chloride-rich water at the highest elevation and a near-neutral spring at lower level may contain a small contribution of crater-lake water. The acid fluid-induced processes at Patuha have led to the accumulation of elements that are commonly associated with volcano-hosted epithermal ore deposits. The dispersal of heavy metals and other potentially toxic elements from the volcano via the local drainage system is a matter of serious environmental concern.     

Peak surface strains during strong earthquake motion

Peak amplitudes of surface strains during strong earthquake ground motion can be approximated by ε = Aν_max/β₁, where ν_max is the corresponding peak particle velocity, β₁ is the velocity of shear waves in the surface layer, and A is a site specific scaling function. In a 50 m thick layer with shear wave velocity β₁ 300 m/s, A 0·4 for the radial strain ε_rr, A 0·2 for the tangential strain ε_rθ, and A 1·0 for the vertical strain, ε_z. These results are site specific and representative of strike slip faulting and of soil in Westmoreland, in Imperial Valley, California. Similar equations can be derived for other sites with known shear wave velocity profile versus depth.