Archaean high-K granitoids produced by remelting of earlier Tonalite–Trondhjemite–Granodiorite (TTG) in the Sangmelima region of the Ntem complex of the Congo craton,southern Cameroon

We present a geochemical and isotopic study that, consistent with observed field relations, suggest Sangmelima late Archaean high-K granite was derived by partial melting of older Archaean TTG. The TTG formations are sodic-trondhjemitic, showing calcic and calc-alkalic trends and are metaluminous to peraluminous. High-K granites in contrast show a potassic calc-alkaline affinity that spans the calcic, calc-alkalic, alkali-calcic and alkalic compositions. The two rock groups (TTG and high-K granites) on the other hand are both ferroan and magnesian. They have a similar degree of fractionation for LREE but a different one for HREE. Nd model ages and Sr/Y ratios define Mesoarchaean and slab-mantle derived magma compositions respectively, with Nb and Ti anomalies indicating a subduction setting for the TTG. Major and trace element in addition to Sr and Nd isotopic compositions support field observations that indicate the derivation of the high-K granitic group from the partial melting of the older TTG equivalent at depth. Geochemical characteristics of the high-K granitic group are therefore inherited features from the TTG protolith and cannot be used for determining their tectonic setting. The heat budget required for TTG partial melting is ascribed to the upwelling of the mantle marked by a doleritic event of identical age as the generated high-K granite melts. The cause of this upwelling is related to linear delamination along mega-shear zones in an intracontinental setting.     

Environmental impact assessment of risk associated with groundwater overdraft remediation in cone of depression,Jining, China

In this paper, a simple scenario and probabilistic approach is used to assess the potential groundwater risk due to proposed overdraft remedial actions in cone of depression, Jining City, China. Focusing on the concentrations of Chloride ions (Cl⁻) and total hardness (TH), the impact of artificial recharge and reduced pumping on groundwater quality and quantity is analysed by using the three-dimensional finite difference groundwater flow and transport model, Visual MODFLOW, to simulate groundwater flow and transport within the study area based on scenarios, and utilizes SURFER software to map risk levels. Although 5, 10 or 15% reduced pumping with artificial recharge leads to more decrease in Cl⁻ and TH concentrations than the 25%, less volume increase is achieved for the remediation of land subsidence and other environmental problems in the cone of depression. The Cl⁻ concentrations in recovered groundwater are within the desired concentration of 200 mg/l; however, TH in some cases are above the maximum permissible limit of 500 mg/l, with an exceedence probability of about 0.67 for recharge and recharge with reduced pumping at 25%. The presence of fractures and hydrogeological complexity greatly determines impacts of remediation, and the 22% reduced pumping with artificial recharge offers an optimum strategy for overdraft remediation in the Jining cone of depression.     

A coupling of mixed and continuous Galerkin finite element methods for poroelasticity II: the discrete-in-time case

In this paper, we formulate a finite element procedure for approximating the coupled fluid and mechanics in Biot’s consolidation model of poroelasticity. Here, we approximate the pressure by a mixed finite element method and the displacements by a Galerkin method. Theoretical convergence error estimates are derived in a discrete-in-time setting. Of particular interest is the case when the lowest-order Raviart–Thomas approximating space or cell-centered finite differences are used in the mixed formulation and continuous piecewise linear approximations are used for displacements. This approach appears to be the one most frequently applied to existing reservoir engineering simulators.     

Estimating groundwater evapotranspiration rates using diurnal water-table fluctuations in a semi-arid riparian zone

In semi-arid climates, phreatophytes draw on shallow aquifers, and groundwater evapotranspiration (ET_G) is a principal component of groundwater budgets. Diurnal water table fluctuations, which often are a product of ET_G, were monitored in the riparian zone of Red Canyon Creek, Wyoming, USA. These fluctuations were higher in a riparian wetland (2–36 mm) than a grass-covered meadow (1–6 mm). The onset and cessation of water-table fluctuations correspond to daily temperatures relative to freezing. Spatial differences were due to vegetation type and specific yield, while temporal changes were due to vegetation dormancy. Ratios of ET_G to potential evapotranspiration (PET), K _c,GW, were similar to ratios of actual evapotranspiration (ET) to PET, K _c, in semi-arid rangelands. Before vegetation senescence, K _c,GW increased between precipitation events, suggesting phreatophytes pull more water from the saturated zone as soil moisture decreases. In contrast, K _c decreases with soil moisture following precipitation events as ET becomes increasingly water-limited. Error in ET_G is primarily from estimates of specific yield (S _y), which is difficult to quantify in heterogeneous sediments. ET_G values may be more reliable because the range of acceptable S _y is smaller than K _c and S _y does not change with vegetation type or soil moisture.     

Recent plumbing system of the Krakatau volcano revealed by teleseismic earthquake distribution

Spatial and temporal analysis of global seismological data 1964–2005 reveals a distinct teleseismic earthquake activity producing a columnar-like formation in the continental wedge between the Krakatau volcano at the surface and the subducting slab of the Indo-Australian plate. These earthquakes occur continuously in time, are in the body-wave (m _b) magnitude range 4.5–5.3 and in the depth range 1–100 km. The Krakatau earthquake cluster is vertical and elongated in the azimuth N30°E, suggesting existence of a deep-rooted fault zone cutting the Sunda Strait in the SSW-NNE direction. Possible continuation of the fault zone in the SW direction was activated by an intensive 2002/2003 aftershock sequence, elongated in the azimuth of N55°E. Beneath the Krakatau earthquake cluster, an aseismic gap exists in the Wadati-Benioff zone of the subducting plate at the depths 100–120 km. We interpret this aseismic gap as a consequence of partial melting inhibiting stress concentration necessary to generate stronger earthquakes, whereas the numerous earthquakes observed in the overlying lithospheric wedge beneath the volcano probably reflect magma ascent in the recent plumbing system of the Krakatau volcano. Focal depth of the deepest events (~100 km) of the Krakatau cluster constrains the location of the primary magma generation to greater depths. The ascending magmatic fluids stress fault segments within the Sunda Strait fault zone and change their friction parameters inducing the observed tectonic earthquakes beneath Krakatau.     

Determining the optimum cell size of digital elevation model for hydrologic application

Scale is one of the most important but unsolved issues in various scientific disciplines that deal with spatial data. The arbitrary choice of grid cell size for contour interpolated digital elevation models (DEM) is one of the major sources of uncertainty in the hydrologic modelling process. In this paper, an attempt was made to identify methods for determining an optimum cell size for a contour interpolated DEM in prior to hydrologic modelling. Twenty-meter interval contour lines were used to generate DEMs of five different resolutions, viz., 30, 45, 60, 75, and 90 m using TOPOGRID algorithm. The obtained DEMs were explored for their intrinsic quality using four different methods, i.e., sink analysis, fractal dimension of derived stream network, entropy measurement and semivariogram modelling. These methods were applied to determine the level artifacts (interpolation error) in DEM surface as well as derived stream network, spatial information content and spatial variability respectively. The results indicated that a 90 m cell size is sufficient to capture the terrain variability for subsequent hydrologic modelling in the study area. The significance of this research work is that it provides methods which DEM users can apply to select an appropriate DEM cell size in prior to detailed hydrologic modelling.     

Changes in nutrient uptake of phytoplankton under the interaction between sunlight and phosphate in the Changjiang （Yangtze） River Estuary

We conducted ship-board incubation experiments to investigate changes in nutrient uptake of phytoplankton under different phosphate concentrations and irradiances in the Changjiang River Estuary and its adjacent waters in China. Under 100% natural irradiance the uptake rates of phosphate, silicate, and nitrate were accelerated at high phosphate levels （1.84 μM）, while under low irradiance （about 50% natural irradiance） their uptake rates were restrained at the high but stimulated greatly at the intermediate phosphate concentrations （1.26 μM）, as the growth of phytoplankton, changes in nitrite and ammonium uptake didn＇t follow an obvious pattern. Our results also showed that there were linear relationships between nitrate, silicate and phosphate uptake at different phosphate concentrations under low and high irradiances, and the growth period of phytoplankton was prolonged both at the high phosphate concentrations under high irradiance and at the intermediate concentrations under low irradiance, suggesting that the limitation of phytoplankton growth mainly reflected changes in its growth period, and because no such environment （low irradiance and low phosphate concentrations） actually existed in a high turbidity zone, phytoplankton blooms hardly occurred there. In the absence of irradiance, denitrification occurred readily and phytoplankton was kept decreasing, which resulted in phosphate regeneration.     

Compressibility of titanosilicate melts

The effect of composition on the relaxed adiabatic bulk modulus (K₀) of a range of alkali- and alkaline earth-titanosilicate [X ₂ ^n/n+ TiSiO₅ (X=Li, Na, K, Rb, Cs, Ca, Sr, Ba)] melts has been investigated. The relaxed bulk moduli of these melts have been measured using ultrasonic interferometric methods at frequencies of 3, 5 and 7 MHz in the temperature range of 950 to 1600°C (0.02 Pa s < _s < 5 Pa s). The bulk moduli of these melts decrease with increasing cation size from Li to Cs and Ca to Ba, and with increasing temperature. The bulk moduli of the Li-, Na-, Ca- and Ba-bearing metasilicate melts decrease with the addition of both TiO₂ and SiO₂ whereas those of the K-, Rb- and Cs-bearing melts increase. Linear fits to the bulk modulus versus volume fraction of TiO₂ do not converge to a common compressibility of the TiO₂ component, indicating that the structural role of TiO₂ in these melts is dependent on the identity of the cation. This proposition is supported by a number of other property data for these and related melt compositions including heat capacity and density, as well as structural inferences from X-ray absorption spectroscopy (XANES). The compositional dependence of the compressibility of the TiO₂ component in these melts explains the difficulty incurred in previous attempts to incorporate TiO₂ in calculation schemes for melt compressibility. The empirical relationship KV^-4/3 for isostructural materials has been used to evaluate the compressibility-related structural changes occurring in these melts. The alkali metasilicate and disilicate melts are isostructural, independent of the cation. The addition of Ti to the metasilicate composition (i.e. X₂TiSiO₅), however, results in a series of melts which are not isostructural. The alkaline-earth metasilicate and disilicate compositions are not isostructural, but the addition of Ti to the metasilicate compositions (i.e. XTiSiO₅) would appear, on the basis of modulus-volume systematics, to result in the melts becoming isostructural with respect to compressibility.