G23FM: a tool for meshing complex geological media

We present G23FM, a mesh generation tool for discretizing two- and three-dimensional complex fractured geological media. G23FM includes different techniques to generate finite element grids that maintain the geometric integrity of input surfaces, and geologic data and produce optimal triangular/tetrahedral grids for flow and transport simulations. G23FM generates grid for two-dimensional cross-sections, represents faults and fractures, for three-dimensional fractured media, and has the capability of including finer grids. Different examples are presented to illustrate some of the main features of G23FM.     

Chemical and Strontium Isotopic Compositions of the Hanjiang Basin Rivers in China: Anthropogenic Impacts and Chemical Weathering

The Hanjiang River, the largest tributaries of the Changjiang (Yangtze) River, is the water source area of the Middle Route of China’s South-to-North Water Transfer Project. The chemical and strontium isotopic compositions of the river waters are determined with the main purpose of understanding the contribution of chemical weathering processes and anthropogenic inputs on river solutes, as well as the associated CO₂ consumption in the carbonate-dominated basin. The major ion compositions of the Hanjiang River waters are characterized by the dominance of Ca²⁺ and HCO₃ ⁻, followed by Mg²⁺ and SO₄ ²⁻. The increase in TDS and major anions (Cl⁻, NO₃ ⁻, and SO₄ ²⁻) concentrations from upstream to downstream is ascribed to both extensive influences from agriculture and domestic activities over the Hanjiang basin. The chemical and Sr isotopic analyses indicate that three major weathering sources (dolomite, limestone, and silicates) contribute to the total dissolved loads. The contributions of the different end-members to the dissolved load are calculated with the mass balance approach. The calculated results show that the dissolved load is dominated by carbonates weathering, the contribution of which accounts for about 79.4% for the Hanjiang River. The silicate weathering and anthropogenic contributions are approximately 12.3 and 6.87%, respectively. The total TDS fluxes from chemical weathering calculated for the water source area (the upper Hanjiang basin) and the whole Hanjiang basin are approximately 3.8 × 10⁶ and 6.1 × 10⁶ ton/year, respectively. The total chemical weathering (carbonate and silicate) rate for the Hanjiang basin is approximately 38.5 ton/km²/year or 18.6 mm/k year, which is higher than global mean values. The fluxes of CO₂ consumption by carbonate and silicate weathering are estimated to be 56.4 × 10⁹ and 12.9 × 10⁹ mol/year, respectively.     

Prediction of Floor Water Inrush: The Application of GIS-Based AHP Vulnerable Index Method to Donghuantuo Coal Mine,China

Floor water inrush represents a geohazard that can pose significant threat to safe operations for instance in coal mines in China and elsewhere. Its occurrence is controlled by many factors, and the processes are often not amenable to mathematical expressions. To evaluate the water inrush risk, the paper proposes the vulnerability index approach by coupling the analytic hierarchy process (AHP) and geographic information system (GIS). The detailed procedures of using this innovative approach are shown in a case study in China (Donghuantuo Coal Mine). The powerful spatial data analysis functions of GIS was used to establish the thematic layer of each of the six factors that control the water inrush, and the contribution weights of each factor was determined with the AHP method. The established AHP evaluation model was used to determine the threshold value for each risk level with a histogram of the water inrush vulnerability index. As a result, the mine area was divided into five regions with different vulnerability levels which served as general guidelines for the mine operations. The prediction results were further corroborated with the actual mining data, and the evaluation result is satisfactory.     

Characteristics of winter mass balance of Glacier No.1 at the headwaters of the Urumqi River,Tianshan Mountains

This study analyzes the response of glacier to climate change during the past 49 years in Urumqi River source region, the Tianshan Mountains of China. The temporal and spatial variations of winter mass balance (bn-w) at different time scales were analyzed to identify their response to climate change during 1988–2006 (The observation of winter mass balance observation began in 1988) on the Glacier No.1 at the headwaters of the Urumqi River, Tianshan Mountains, China. The winter accumulation shows a significantly decreasing trend. The results show that the cumulative values on Glacier No.1 is 2,202 mm water equivalent during 1988–2006 and the mean values is 116 mm a⁻¹. Furthermore, the trend analysis of the winter mass balance indicates a rapid decrease since 1990, and the mean mass balance is only 79 mm a⁻¹ during 1997–2006. Winter mass balance correlates well negatively with the total evaporation from September to April (r = −0.68, α = 0.01), and positively with the total precipitation from September to April (r = 0.74, α = 0.01). However, winter mass balance shows a weak correlation with mean minimum air temperature during September to April (r = −0.35), and runoff on September (r = −0.13).     

Metasomatized lithospheric mantle beneath Turkana depression in southern Ethiopia (the East Africa Rift): geochemical and Sr–Nd–Pb isotopic characteristics

Mantle xenoliths entrained in Quaternary alkaline basalts from the Turkana Depression in southern Ethiopia (the East Africa Rift) were studied for their geochemical and Sr–Nd–Pb isotopic compositions to constrain the evolution of the lithosphere. The investigated mantle xenoliths are spinel lherzolites in composition with a protogranular texture. They can be classified into two types: anhydrous and hydrous spinel lherzolites; the latter group characterized by the occurrences of pargasite and phlogopite. The compositions of whole-rock basaltic component (CaO = 3.8–5.6 wt%, Al₂O₃ = 2.5–4.1 wt%, and MgO = 34.7–38.1 wt%), spinel (Cr# = 0.062–0.117, Al₂O₃ = 59.0–64.4 wt%) and clinopyroxene (Mg# = 88.4–91.7, Al₂O₃ = 5.2–6.7 wt%) indicate that the lherzolites are fertile and have not experienced significant partial melting. Both types are characterized by depleted ⁸⁷Sr/⁸⁶Sr (0.70180–0.70295) and high ¹⁴³Nd/¹⁴⁴Nd (0.51299–0.51348) with wide ranges of ²⁰⁶Pb/²⁰⁴Pb (17.86–19.68) isotopic compositions. The variations of geochemical and isotopic compositions can be explained by silicate metasomatism induced by different degree of magma infiltrations from ascending mantle plume. The thermobarometric estimations suggest that the spinel lherzolites were derived from depths of 50–70 km (15.6–22.2 kb) and entrained in the alkaline magma at 847–1,052°C. Most of the spinel lherzolites from this study record an elevated geotherm (60–90 mW/m²) that is related to the presence of rising mantle plume in an active tectonic setting. Sm–Nd isotopic systematic gives a mean T_DM model age of 0.95 Ga, interpreted as the minimum depletion age of the subcontinental lithosphere beneath the region.     

Threat of land subsidence in and around Kolkata City and East Kolkata Wetlands,West Bengal,India

This paper attempts to estimate the possible rate of land subsidence of Kolkata City including Salt Lake City and the adjoining East Kolkata Wetlands located at the lower part of the deltaic alluvial plain of South Bengal basin. Demand of groundwater for drinking, agricultural and industrial purposes has increased due to rapid urbanization. The subsurface geology consists of Quaternary sediments comprising a succession of clay, silty clay and sand of various grades. Groundwater occurs mostly under confined condition except in those places where the top aquitard has been obliterated due to the scouring action of past channels. Currently, the piezometric head shows a falling trend and it may be accelerated due to further over-withdrawal of groundwater resulting in land subsidence. The estimated mean land subsidence rate is 13.53 mm/year and for 1 m drop in the piezometric head, the mean subsidence is 3.28 cm. The surface expression of the estimated land subsidence is however, cryptic because of a time lag between the settlement of the thick low-permeable aquitard at the top and its surface expression. Therefore, groundwater of the cities and wetland areas should be developed cautiously based on the groundwater potential to minimize the threat of land subsidence.     

Response of xenotime to prograde metamorphism

Xenotime is a widespread accessory mineral in lower greenschist to upper amphibolite facies metasedimentary rocks from the Palaeoproterozoic Mount Barren Group, southwestern Australia. Xenotime is closely associated with detrital zircon, commonly forming syntaxial outgrowths, in samples of sandstone, micaceous quartzite, slate, phyllite, garnet-bearing semi-pelites, and in kyanite-, garnet-, and staurolite-bearing mica schists. In situ geochronology of xenotime from lower greenschist sandstones has previously yielded multiple U–Pb ages with peaks at ~2.0, ~1.7, and ~1.65 Ga, interpreted to represent the age of detritus, early diagenesis, and a later thermal event, respectively. New U–Pb dating of xenotime in slate yields a major population at ~1.7 Ga with a minor population at ~1.2 Ga, reflecting diagenetic and metamorphic growth, respectively, whereas xenotime in phyllite forms a minor age population at ~1.7 Ga and a main peak at ~1.2 Ga. Mid-greenschist facies semi-pelitic schists (quartz-muscovite-garnet) contain xenotime that formed before 1.8 Ga and at 1.2 Ga, representing detrital and peak metamorphic ages, respectively. Xenotime in samples of amphibolite facies schist (650°C and ~8 kbars) yields U–Pb ages of ~1.2 Ga, coinciding with the time of peak metamorphism. A single analysis of a xenotime core from an amphibolite facies schist gave an age of ~1.8 Ga, consistent with the presence of detrital xenotime. Our results suggest that detrital xenotime may be preserved under greenschist facies conditions, but is largely replaced during upper amphibolite facies conditions. Detrital xenotime is replaced through dissolution–reprecipitation reactions forming compositionally distinct rims during greenschist and amphibolite facies metamorphism at 1.2 Ga. Diagenetic xenotime is present in lower greenschist facies samples, but was not observed in metasedimentary rocks that had experienced temperatures above mid-greenschist facies metamorphism (450°C). The apparent disappearance of detrital and diagenetic xenotime and appearance of metamorphic xenotime during prograde metamorphism indicates that some of the yttrium, heavy rare earth elements, and phosphorus needed for metamorphic xenotime growth are probably derived from the replacement of detrital and diagenetic xenotime.     

Barium-rich fluids and melts in a subduction environment (La Corea and Sierra del Convento mélanges,eastern Cuba)

Whole-rock compositions of muscovite-bearing amphibolite, trondhjemite, pegmatite and quartz-muscovite rocks from Sierra del Convento and La Corea mélanges (eastern Cuba), as well as mineralogy, record complex circulation of Ba-rich fluids and melts in the subduction environment. Partial melting of fluid-fluxed, MORB-derived amphibolite produced trondhjemite magmas that crystallized at depth, in some cases evolving into pegmatites. Qtz-Ms rocks probably crystallized from primary fluids derived from subducted sediments. All these rocks have elevated concentrations of large-ion lithophile elements, especially Ba (up to 11,810 ppm), presumably released from slab sediments by fluids and/or melts. Fluid–rock interaction produced crystallization of phengite in parental amphibolites. The phengite crystallized in all types of rocks is rich in Ba, with concentric zoning, characterized by Ba-rich cores and Ba-poor rims, indicating a compatible behaviour of Ba in the studied systems. Zoning in phengite is governed primarily by the celadonite (tschermak) exchange vector ((Mg,Fe)Si-(Al)_-2), with more moderate contributions of the celsian (BaAl-(KSi)_-1) and paragonite (NaK_-1) exchange vectors. Late remobilization of Ba at relatively low temperature formed retrograde celsian. The compatible behaviour of Ba in the studied rocks strengthens the importance of the stability of phengite for the transfer of LIL elements from the subduction to the volcanic arc environments.     

Modeling a major source of economic resilience to disasters: recapturing lost production

Terrorist attacks and natural disasters have potentially severe economic consequences in terms of property damage and business interruption. However, experience from the September 11 World Trade Center attack and other disasters indicates that the economy has a great deal of resilience. This refers to the ability to dampen the maximum potential economic output (business interruption) loss. One of the most prominent sources of resilience is the ability of businesses to reschedule, or recapture, lost production after the event. Although there have been applications of a fixed parameter recapture factor for each of several aggregated sectors of the economy, there has been little formal analysis of this resilience action. This study offers a theoretic framework for analyzing production rescheduling. It distinguishes the major conditions influencing two aspects that have previously been neglected: (1) the maximum time span over which the rescheduling can take place and (2) the likely decline of the maximum recapture as the business interruption increases. We divide the relevant time path into two periods after recovery. One is a function of a recaptured output path after recovery to the status of normal production. The other is a function for the maximum recaptured production, based on the recaptured output path. The recaptured output path function is assumed to follow a normal distribution function, and hence, total recaptured output follows the cumulative normal distribution function over time after productive capacity is restored. Also, we develop a new cumulative normal distribution function for interruption time duration, which is symmetric with respect to the output axis. This recapture function has unknown parameters. Empirical data on the recaptured amounts following an actual disaster can be used to estimate the parameters of this function using simulation methods.     

Characteristics of pegmatoidal granite exposed near Bayalan,Ajmer district,Rajasthan

The study involves the characterization of pegmatoidal granite, southeast of Beawar, Ajmer district, Rajasthan. Earlier researchers had described this granite as part of the BGC, basement to the Bhim Group of the Delhi Super Group rocks. However, the present study indicates that it is younger than the rocks of Bhim Group of South Delhi Fold Belt, into which it is intrusive. The intrusion is structurally controlled and the outcrop pattern is phacolithic. The granite had intruded post-D₂ deformation of the Delhi orogeny along the axial planes of D₂ folds. The intrusion has also resulted in the formation of a contact aureole about the calc gneisses.