Sustainable groundwater development in the coastal Tra Vinh province in Vietnam under saltwater intrusion and climate change

Three-dimensional transient groundwater flow and saltwater transport models were constructed to assess the impacts of groundwater abstraction and climate change on the coastal aquifer of Tra Vinh province (Vietnam). The groundwater flow model was calibrated with groundwater levels (2007–2016) measured in 13 observation wells. The saltwater transport model was compared with the spatial distribution of total dissolved solids. Model performance was evaluated by comparing observed and simulated groundwater levels. The projected rainfalls from two climate models (MIROC5 and CRISO Mk3.6) were subsequently used to simulate possible effects of climate changes. The simulation revealed that groundwater is currently depleted due to overabstraction. Towards the future, groundwater storage will continue to be depleted with the current abstraction regime, further worsening in the north due to saltwater intrusion from inland trapped saltwater and on the coast due to seawater intrusion. Notwithstanding, the impact from climate change may be limited, with the computed groundwater recharge from the two climate models revealing no significant change from 2017 to 2066. Three feasible mitigation scenarios were analyzed: (1) reduced groundwater abstraction by 25, 35 and 50%, (2) increased groundwater recharge by 1.5 and 2 times in the sand dunes through managed aquifer recharge (reduced abstraction will stop groundwater-level decline, while increased recharge will restore depleted storage), and (3) combining 50% abstraction reduction and 1.5 times recharge increase in sand dune areas. The results show that combined interventions of reducing abstraction and increasing recharge are necessary for sustainable groundwater resources development in Tra Vinh province.

     

The Role of Lithospheric Mantle Heterogeneity in the Generation of Plio-Pleistocene Alkali Basaltic Suites from NW Harrat Ash Shaam (Israel)

Plio-Pleistocene volcanism in the Golan and Galilee (northeasternIsrael) shows systematic variability with time and location:alkali basalts were erupted in the south during the Early Pliocene,whereas enriched basanitic lavas erupted in the north duringthe Late Pliocene (Galilee) and Pleistocene (Golan). The basaltsshow positive correlations in plots of ratios of highly to moderatelyincompatible elements versus the concentration of the highlyincompatible element (e.g. Nb/Zr vs Nb, La/Sm vs La) and indiagrams of REE/HFSE (rare earth elements/high field strengthelements) vs REE concentration (e.g. La/Nb vs La). Some of thesecorrelations are not linear but upward convex. ⁸⁷Sr/⁸⁶Sr ratiosvary between 0·7031 and 0·7034 and correlate negativelywith incompatible element concentrations and positively withRb/Sr ratios. We interpret these observations as an indicationthat the main control on magma composition is binary mixingof melts derived from two end-member mantle source components.Based on the high Sr/Ba ratios and negative Rb anomalies inprimitive mantle normalized trace element diagrams and the moderateslopes of MREE–HREE (middle REE–heavy REE) in chondrite-normalizeddiagrams, we suggest that the source for the alkali basalticend-member was a garnet-bearing amphibole peridotite that hadexperienced partial dehydration. The very high incompatibleelement concentrations, low K content, very low Rb contentsand steep MREE–HREE patterns in the basanites are attributedto derivation from amphibole- and garnet-bearing pyroxeniteveins. It is suggested that the veins were produced via partialmelting of amphibole peridotites, followed by complete solidificationand dehydration that effectively removed Rb and K. The requirementfor the presence of amphibole limits both sources to lithosphericdepths. The spatial geochemical variability of the basalts indicatesthat the lithosphere beneath the region is heterogeneous, composedof vein-rich and vein-poor domains. The relatively uniform ¹⁴³Nd/¹⁴⁴Nd(Nd = 4·0–5·2) suggests that the two mantlesources were formed by dehydration and partial melting of anoriginally isotopically uniform reservoir, probably as a resultof a Paleozoic thermal event. KEY WORDS: basanites; lithospheric heterogeneity; magma mixing; amphibole peridotite; pyroxenites     

The timing of ultrahigh-temperature metamorphism in Southern India: U–Th–Pb electron microprobe ages from zircon and monazite in sapphirine-bearing granulites

We report here U–Pb electron microprobe ages from zircon and monazite associated with corundum- and sapphirine-bearing granulite facies rocks of Lachmanapatti, Sengal, Sakkarakkottai and Mettanganam in the Palghat–Cauvery shear zone system and Ganguvarpatti in the northern Madurai Block of southern India. Mineral assemblages and petrologic characteristics of granulite facies assemblages in all these localities indicate extreme crustal metamorphism under ultrahigh-temperature (UHT) conditions. Zircon cores from Lachmanapatti range from 3200 to 2300 Ma with a peak at 2420 Ma, while those from Mettanganam show 2300 Ma peak. Younger zircons with peak ages of 2100 and 830 Ma are displayed by the UHT granulites of Sengal and Ganguvarpatti, although detrital grains with 2000 Ma ages are also present. The Late Archaean-aged cores are mantled by variable rims of Palaeo- to Mesoproterozoic ages in most cases. Zircon cores from Ganguvarpatti range from 2279 to 749 Ma and are interpreted to reflect multiple age sources. The oldest cores are surrounded by Palaeoproterozoic and Mesoproterozoic rims, and finally mantled by Neoproterozoic overgrowths. In contrast, monazites from these localities define peak ages of between 550 and 520 Ma, with an exception of a peak at 590 Ma for the Lachmanapatti rocks. The outermost rims of monazite grains show spot ages in the range of 510–450 Ma.While the zircon populations in these rocks suggest multiple sources of Archaean and Palaeoproterozoic age, the monazite data are interpreted to date the timing of ultrahigh-temperature metamorphism in southern India as latest Neoproterozoic to Cambrian in both the Palghat–Cauvery shear zone system and the northern Madurai Block. The data illustrate the extent of Neoproterozoic/Cambrian metamorphism as India joined the Gondwana amalgam at the dawn of the Cambrian.     

Early Neoproterozoic granulite facies metamorphism of mafic dykes from the Vestfold Block,east Antarctica

Proterozoic mafic dykes from the southwestern Vestfold Block experienced heterogeneous granulite facies metamorphism, characterized by spotted or fractured garnet‐bearing aggregates in garnet‐absent groundmass. The garnet‐absent groundmass typically preserves an ophitic texture composed of lathy plagioclase, intergranular clinopyroxene and Fe–Ti oxides. Garnet‐bearing domains consist mainly of a metamorphic assemblage of garnet, clinopyroxene, orthopyroxene, hornblende, biotite, plagioclase, K‐feldspar, quartz and Fe–Ti oxides. Chemical compositions and textural relationships suggest that these metamorphic minerals reached local equilibrium in the centre of the garnet‐bearing domains. Pseudosection calculations in the model system NCFMASHTO (Na₂O–CaO–FeO–MgO–Al₂O₃–SiO₂–H₂O–TiO₂–Fe₂O₃) yield P–T estimates of 820–870 °C and 8.4–9.7 kbar. Ion microprobe U–Pb zircon dating reveals that the NW‐ and N‐trending mafic dykes were emplaced at 1764 ± 25 and 1232 ± 12 Ma, respectively, whereas their metamorphic ages cluster between 957 ± 7 and 938 ± 9 Ma. The identification of granulite facies mineral inclusions in metamorphic zircon domains is also consistent with early Neoproterozoic metamorphism. Therefore, the southwestern margin of the Vestfold Block is inferred to have been buried to depths of ~30–35 km beneath the Rayner orogen during the late stage of the late Mesoproterozoic/early Neoproterozoic collision between the Indian craton and east Antarctica (i.e. the Lambert Terrane or the Ruker craton including the Lambert Terrane). The lack of penetrative deformation and intensive fluid–rock interaction in the rigid Vestfold Block prevented the nucleation and growth of garnet and resulted in the heterogeneous granulite facies metamorphism of the mafic dykes.     

Emplacement ages and sources of kimberlites and related rocks in southern Africa: U–Pb ages and Sr–Nd isotopes of groundmass perovskite

Groundmass perovskite has been dated by LA-ICPMS in 135 kimberlites and related rocks from 110 localities across southern Africa. Sr and/or Nd isotopes have been analysed by LA-MC-ICPMS in a subset of these and integrated with published data. The age distribution shows peaks at 1,600–1,800, 1,000–1,200, 500–800 and 50–130 Ma. The major “bloom” of Group I kimberlites at ca 90 ± 10 Ma was preceded by a slow build-up in magmatic activity from ca 180 Ma. The main pulse of Group II kimberlites at 120–130 Ma was a distinct episode within this build-up. Comparison of the isotopic data with seismic tomography images suggests that metasomatized subcontinental lithospheric mantle (SCLM) with very low ε _Nd and high ⁸⁷Sr/⁸⁶Sr, (the isotopic signature of Group II kimberlites) was focused in low-Vs zones along translithospheric structures. Such metasomatized zones existed as early as 1,800 Ma, but were only sporadically tapped until the magmatic build-up began at ca 180 Ma, and contributed little to the kimberlitic magmas after ca 110 Ma. We suggest that these metasomatized volumes resided in the deep SCLM and that their low-melting point components were “burned off” by rising temperatures, presumably during an asthenospheric upwelling that led to SCLM thinning and a rise in the ambient geotherm between 120 and 90 Ma. The younger Group I kimberlites therefore rarely interacted with such SCLM, but had improved access to shallower volumes of differently metasomatized, ancient SCLM with low ⁸⁷Sr/⁸⁶Sr and intermediate ε _Nd (0–5). The kimberlite compositions therefore reflect the evolution of the SCLM of southern Africa, with metasomatic-enrichment events from as early as 1.8 Ga, through a major thermal and compositional change at ca 110 Ma, and the major kimberlite “bloom” around 90 Ma.     

Fluid inclusion microthermometry and the P-T evolution of gold-bearing hydrothermal fluids in the Niuxinshan gold deposit, eastern Hebei province, NE China

Niuxinshan is a typical example of the numerous mesothermal gold deposits formed during Mesozoic tectono-magmatic reactivation of the Archean North China Craton in eastern Hebei province. Gold occurs in quartz-sulfide lodes in Archean amphibolites and also in greisen zones in the Mesozoic Niuxinshan granite stock. Four mineralization stages can be recognized from early to late: (1) quartz-K-feldspar, (2) quartz-pyrite, (3) quartz-polysulfide, and (4) quartz-carbonate. Gold mineralization mainly occurs in stages 2 and 3. Fluid inclusions in quartz and fluorite from greisen zones in the Niuxinshan granite, and inclusions in vein quartz and sphalerite from stages 1 to 3 in the amphibolites, have been studied by microthermometry. Three compositional types of inclusions are recognized: type 1 (Tp1) are H₂O-CO₂-bearing inclusions and include primary (Tp1-P) and secondary (Tp1-S) inclusions. These are found in quartz and fluorite from the greisen zones as well as in vein quartz and sphalerite from stages 1 to 3. The Tp1-P inclusions are considered to represent the gold-bearing hydrothermal fluids. Type 2 (Tp2-S) are secondary H₂O-CO₂ + solid phase inclusions in fluorite from the greisen zones. Type 3 (Tp3-S) are secondary aqueous inclusions with a solid phase which coexist with the Tp2-S in fluorite from the greisen zones. The Tp1-P inclusions show variable V_CO2 (commonly 0.3 to 0.6) and X_CO2 values (mainly 0.1 to 0.4). The salinities of inclusions cluster around 3 to 11 wt.% NaCl equivalent and their homogenization temperatures to the liquid phase (Th(L)) fall dominantly in the range of 260 to 360 °C. The compositional variations of inclusions in stage 1 probably result from exsolution of magmatic fluids at various stages; immiscibility or boiling of the fluids can be ruled out. The compositional variations of inclusions in the greisen zones and in vein stages 2 and 3 are attributed to cooling, mixing (dilution), and necking-down of the fluids. The Tp1-S and Tp2-S inclusions show salinities of 3 to 6 wt.% NaCl equivalent and X_CO2 values of 0.04 to 0.17. Th(L) clusters at 240 to 260 °C. The Tp3-S inclusions have salinities of 3 to 6 wt.% NaCl equivalent and Th(L) of 170 to 240 °C. Isochoric reconstructions, combined with oxygen and sulfur isotope geothermometry of mineral pairs, give trapping P-T conditions for the gold-bearing fluids. The greisen zones formed at 310 to 460 °C and 1.3 to 3.7 kbar; stage 1 veins at 300 to 430 °C and 1.2 to 3.7 kbar; stage 2 veins at 290 to 380 °C and 1 to 3 kbar; stage 3 veins at 250 to 350 °C and 1 to 3 kbar. H₂O-CO₂ fluids with low to moderate salinities and moderate to high densities (0.66 to 1.01 g/cm³) dominated at early mineralization stages, and evolved towards H₂O-richer and CO₂- and less saline fluids through time. The retrograde P-T evolution probably resulted from regional uplift and cooling of gold-bearing hydrothermal fluids. The gold bisulfide complex was dominant in the fluids during mineralization and gold deposition was mainly induced by decreases of temperature and pressure, as well as destabilization of the bisulfide complex during sulfidization of wall rocks. Received: 16 March 1998 / Accepted: 11 January 1999     

Evaluation of heavy metal contamination and implication of multiple sources from Hunchun basin, northeastern China

Present concentrations and distributions of heavy metals through profiles, surface soil, and stream sediment samples in the Hunchun area, north-eastern China, were investigated to determine the elemental background values. This study also aims to characterize potentially toxic materials such as pulverized fly ash (PFA) from power stations or ash and slag from coal used domestically in urban areas, agrochemicals applied inappropriately, and urban sewage sludges from Hunchun City, as well as to ascertain the possibility of natural enrichment through site characterization by mineralogical and geochemical investigation. The distribution of contaminants in the alluvial soils (fluvisol) of this area has been influenced by several interacting factors. The parent alluvial materials from weathered products of amphiboles have made coatings such as ferrihydrite, goethite, and hematite. This natural inheritance factor is supported by the fact that the concentrations of weak acid-extractable (plant-available) heavy metals are very low, except for Fe and Mn. However, in agricultural soils and adjacent stream sediments, an anthropogenic input of Cd, Pb, Ni and Cr by agrochemicals is strongly suggested. Also, F contamination by coal combustion and the dissolution of F-bearing minerals could cause some future problems. Wide distribution and significantly high concentrations of Cd, Fe, Mn, and F in soils throughout the combination of pollutants originating from lithogenic and the anthropogenic sources pose potential problems in utilizing water resources. Received: 14 June 1999 · Accepted: 27 December 1999     

Linear and nonlinear vibration of non-uniform beams on two-parameter foundations using p-elements

A hierarchical finite element is presented for the geometrically nonlinear free and forced vibration of a non-uniform Timoshenko beam resting on a two-parameter foundation. Legendre orthogonal polynomials are used as enriching shape functions to avoid the shear-locking problem. With the enriching degrees of freedom, the accuracy of the computed results and the computational efficiency are greatly improved. The arc-length iterative method is used to solve the nonlinear eigenvalue equation. The computed results of linear and nonlinear vibration analyses show that the convergence of the proposed element is very fast with respect to the number of Legendre orthogonal polynomials used. Since the elastic foundation and the axial load applied at both ends of the beam affect the ratios of linear frequencies associated with the internal resonance, they influence the nonlinear vibration characteristics of the beam. The axial tensile stress of the beam in nonlinear vibration is investigated in this paper, and attention should be paid to the geometrically nonlinear vibration resulting in considerably large axial tensile stress in the beam.     

Application of an optimum design technique for determining the coefficient of consolidation by using piezocone test data

For normally consolidated clay, several researchers have developed a number of theoretical time factors to determine the coefficient of consolidation from piezocone test results. However, depending on assumptions and analytical techniques, it could vary considerably, even for a specific degree of consolidation. In this paper a method is proposed to determine a consistent coefficient of consolidation by applying the concept of an optimum design technique over all ranges of the degree of consolidation. Initial excess pore pressure distribution is assumed to be capable of being obtained by the successive spherical cavity expansion theory. The dissipation of pore pressure is simulated by means of a two-dimensional linear-uncoupled axi-symmetric consolidation analysis. The minimization of differences between measured and predicted excess pore pressure was carried out by the BFGS unconstrained optimum design algorithm with a one-dimensional golden section search technique. By analyzing numerical examples and in-situ test results, it was found that the adopted optimum design technique gives consistent and convergent results.