Application of MODFLOW and geographic information system to groundwater flow simulation in North China Plain, China

MODFLOW is a groundwater modeling program. It can be compiled and remedied according to the practical applications. Because of its structure and fixed data format, MODFLOW can be integrated with Geographic Information Systems (GIS) technology for water resource management. The North China Plain (NCP), which is the politic, economic and cultural center of China, is facing with water resources shortage and water pollution. Groundwater is the main water resource for industrial, agricultural and domestic usage. It is necessary to evaluate the groundwater resources of the NCP as an entire aquifer system. With the development of computer and internet information technology it is also necessary to integrate the groundwater model with the GIS technology. Because the geological and hydrogeological data in the NCP was mainly in MAPGIS format, the powerful function of GIS of disposing of and analyzing spatial data and computer languages such as Visual C and Visual Basic were used to define the relationship between the original data and model data. After analyzing the geological and hydrogeological conditions of the NCP, the groundwater flow numerical simulation modeling was constructed with MODFLOW. On the basis of GIS, a dynamic evaluation system for groundwater resources under the internet circumstance was completed. During the process of constructing the groundwater model, a water budget was analyzed, which showed a negative budget in the NCP. The simulation period was from 1 January 2002 to 31 December 2003. During this period, the total recharge of the groundwater system was 49,374 × 10⁶ m³ and the total discharge was 56,530 × 10⁶ m³ the budget deficit was −7,156 × 10⁶ m³. In this integrated system, the original data including graphs and attribution data could be stored in the database. When the process of evaluating and predicting groundwater flow was started, these data were transformed into files that the core program of MODFLOW could read. The calculated water level and drawdown could be displayed and reviewed online.     

Numerical modeling of rainstorm-induced shallow landslides in saturated and unsaturated soils

For the assessment of shallow landslides triggered by rainfall, the physically based model coupling the infinite slope stability analysis with the hydrological modeling in nearly saturated soil has commonly been used due to its simplicity. However, in that model the rainfall infiltration in unsaturated soil could not be reliably simulated because a linear diffusion-type Richards’ equation rather than the complete Richards’ equation was used. In addition, the effect of matric suction on the shear strength of soil was not actually considered. Therefore, except the shallow landslide in saturated soil due to groundwater table rise, the shallow landslide induced by the loss in unsaturated shear strength due to the dissipation of matric suction could not be reliably assessed. In this study, a physically based model capable of assessing shallow landslides in variably saturated soils is developed by adopting the complete Richards’ equation with the effect of slope angle in the rainfall infiltration modeling and using the extended Mohr–Coulomb failure criterion to describe the unsaturated shear strength in the soil failure modeling. The influence of rainfall intensity and duration on shallow landslide is investigated using the developed model. The result shows that the rainfall intensity and duration seem to have similar influence on shallow landslides respectively triggered by the increase of positive pore water pressure in saturated soil and induced by the dissipation of matric suction in unsaturated soil. The rainfall duration threshold decreases with the increase in rainfall intensity, but remains constant for large rainfall intensity.     

Modeling of water–rock interactions in an aquitard of sandy clay in the coastal area near Beihai,China

Unconsolidated sand, gravel and clay deposits near Beihai and in the Leizhou Peninsula in southern China form an unconfined aquifer, aquitard and a confined aquifer. Water and soil samples were collected from the two aquifers in the coastal Beihai area for the determination of chemical compositions, minerals and soluble ions. Hydrogeochemical modeling of three flow paths through the aquitard are carried out using PHREEQC to determine water–rock interactions along the flow paths. The results indicate that the dissolution of anorthite, fluorite, halite, rhodochrosite and CO₂, and precipitation of potash feldspar and kaolinite may be occurring when groundwater leaks through the aquitard from the unconfined aquifer to the confined aquifer. Cation exchanges between Na and Ca can also happen along the flow paths.     

Shear strength characterization of municipal solid waste at the Suzhou landfill, China

The current practice of slope stability analysis for a municipal solid waste (MSW) landfill usually overlooks the dependence of waste properties on the fill age or embedment depth. Changes in shear strength of MSW as a function of fill age were investigated by performing field and laboratory studies on the Suzhou landfill in China. The field study included sampling from five boreholes advanced to the bottom of the landfill, cone penetration tests and monitoring of pore fluid pressures. Twenty-six borehole samples representative of different fill ages (0 to 13 years) were used to perform drained triaxial compression tests. The field and laboratory study showed that the waste body in the landfill can be sub-divided into several strata corresponding to different ranges of fill age. Each of the waste strata has individual composition and shear strength characteristics. The triaxial test results showed that the MSW samples exhibited a strain-hardening and contractive behavior. As the fill age of the waste increased from 1.7 years to 11 years, the cohesion mobilized at a strain level of 10% was found to decrease from 23.3 kPa to 0 kPa, and the mobilized friction angle at the same strain level increasing from 9.9° to 26°. For a confinement stress level greater than 50 kPa, the shear strength of the recently-placed MSW seemed to be lower than that of the older MSW. This behavior was consistent with the cone penetration test results. The field measurement of pore pressures revealed a perched leachate mound above an intermediate cover of soils and a substantial leachate mound near the bottom of the landfill. The measurements of shear strength properties and pore pressures were utilized to assess the slope stability of the Suzhou landfill.     

Assessment of disaster resilience capacity of hillslope communities with high risk for geological hazards

This study presents a novel preparedness assessment method for assessing hazard mitigation and environmental planning of hillslope communities. A professional questionnaire was utilized to weight each indicator. Communities in Hsinchu, Taichung and Nantou counties with debris flow hazards were taken as study samples. Debris flow risk and landslide susceptibility for each community were determined using Geographic Information System (GIS) technology and logistic regression analysis. Thus, a novel risk assessment method for evaluating disaster resilience capacity of hillslope communities was established. This method was then applied to assess casualties caused by Typhoon Herb in 1996 and Typhoon Mindulle in 2004. Additionally, the analytical results generated by this assessment method were discussed with the aim of developing references for implementation of risk analysis, increasing the effectiveness of disaster mitigation, and reducing future loss of life and property.     

Silicon isotope fractionation in bamboo and its significance to the biogeochemical cycle of silicon

A systematic investigation on silica contents and silicon isotope compositions of bamboos was undertaken. Seven bamboo plants and related soils were collected from seven locations in China. The roots, stem, branch and leaves for each plant were sampled and their silica contents and silicon isotope compositions were determined. The silica contents and silicon isotope compositions of bulk and water-soluble fraction of soils were also measured. The silica contents of studied bamboo organs vary from 0.30% to 9.95%. Within bamboo plant the silica contents show an increasing trend from stem, through branch, to leaves. In bamboo roots the silica is exclusively in the endodermis cells, but in stem, branch and leaves, the silica is accumulated mainly in epidermal cells. The silicon isotope compositions of bamboos exhibit significant variation, from −2.3‰ to 1.8‰, and large and systematic silicon isotope fractionation was observed within each bamboo. The δ³⁰Si values decrease from roots to stem, but then increase from stem, through branch, to leaves. The ranges of δ³⁰Si values within each bamboo vary from 1.0‰ to 3.3‰. Considering the total range of silicon isotope composition in terrestrial samples is only 7‰, the observed silicon isotope variation in single bamboo is significant and remarkable. This kind of silicon isotope variation might be caused by isotope fractionation in a Rayleigh process when SiO₂ precipitated in stem, branches and leaves gradually from plant fluid. In this process the Si isotope fractionation factor between dissolved Si and precipitated Si in bamboo (α_pre-sol) is estimated to be 0.9981. However, other factors should be considered to explain the decrease of δ³⁰Si value from roots to stem, including larger ratio of dissolved H₄SiO₄ to precipitated SiO₂ in roots than in stem. There is a positive correlation between the δ³⁰Si values of water-soluble fractions in soils and those of bulk bamboos, indicating that the dissolved silicon in pore water and phytoliths in soil is the direct sources of silicon taken up by bamboo roots. A biochemical silicon isotope fractionation exists in process of silicon uptake by bamboo roots. Its silicon isotope fractionation factor (α_bam-wa) is estimated to be 0.9988. Considering the distribution patterns of SiO₂ contents and δ³⁰Si values among different bamboo organs, evapotranspiration may be the driving force for an upward flow of a silicon-bearing fluid and silica precipitation. Passive silicon uptake and transportation may be important for bamboo, although the role of active uptake of silicic acid by roots may not be neglected. The samples with relatively high δ³⁰Si values all grew in soils showing high content of organic materials. In contrast, the samples with relatively low δ³⁰Si values all grew in soil showing low content of organic materials. The silicon isotope composition of bamboo may reflect the local soil type and growth conditions. Our study suggests that bamboos may play an important role in global silicon cycle.     

Soil n-alkane δD vs. altitude gradients along Mount Gongga, China

The altitude effect on the isotopic composition of precipitation and its application to paleoelevation reconstruction using authigenic or pedogenic minerals have been intensively studied. However, there are still no studies on variations in biomarker δD along altitude transects to investigate its potential as a paleoelevation indicator, although it has been observed that δD of higher plant lipid may record changes in precipitation δD (δD_p). Here, we present δD values of higher plant-derived C₂₇, C₂₉, and C₃₁n-alkanes from surface soil along the eastern slope of Mount Gongga, China with great changes in physical variables and vegetation over a range from 1000 to 4000 m above sea level. The weighted-mean δD values of these n-alkanes (δD_wax) show significant linear correlations with predicted δD_p values (R² = 0.76) with an apparent isotopic enrichment (ε_wax-p) of −137 ± 9‰, indicating that soil δD_wax values track overall δD_p variation along the entire altitudinal transect. Leaf δD_wax is also highly correlated with mountain altitude by a significant quadratic relationship (R² = 0.80). Evapotranspiration is found declining with altitude, potentially lowering δD_wax values at higher elevations. However, this evapotranspiration effect is believed to be largely compensated by the opposing effect of vegetation changes, resulting in less varied ε_wax-p values over the slope transect. This study therefore confirms the potential of using leaf δD_wax to infer paleoelevations, and more generally, to infer the δD of precipitation.     

Church-based social capital, networks and geographical scale: Katrina evacuation, relocation, and recovery in a New Orleans Vietnamese American community

This research examines the role of social capital and networks to explain the evacuation, relocation, and recovery experiences of a Vietnamese American community in New Orleans, Louisiana in the aftermath of Hurricane Katrina. As the single largest community institution, the parish church’s complex bonding and bridging social capital and networks proved particularly critical in part because of its historically based ontological security. The process of evacuation, but especially relocation and recovery, was dependent on deploying co-ethnic social capital and networks at a variety of geographical scales. Beyond the local or community scale, extra-local, regional, and national scales of social capital and networks reproduced a spatially redefined Vietnamese American community. Part of the recovery process included constructing discursive place-based collective-action frames to successfully contest a nearby landfill that in turn engendered social capital and networks crossing ethnic boundaries to include the extra-local African American community. Engaging social capital and networks beyond the local geographical scale cultivated a Vietnamese American community with an emergent post-Katrina cultural and political identity.     

Origin of Fe-Ti Oxide Ores in Mafic Intrusions: Evidence from the Panzhihua Intrusion, SW China

Economic concentrations of Fe–Ti oxides occur as massive,conformable lenses or layers in the lower part of the Panzhihuaintrusion, Emeishan Large Igneous Province, SW China. Mineralchemistry, textures and QUILF equilibria indicate that oxidesin rocks of the intrusion were subjected to extensive subsolidusre-equilibration and exsolution. The primary oxide, reconstructedfrom compositions of titanomagnetite in the ores and associatedintergrowths, is an aluminous titanomagnetite (Usp₄₀) with 40wt % FeO, 34 wt % Fe₂O₃, 16·5 wt % TiO₂, 5·3 wt% Al₂O₃, 3·5 wt % MgO and 0·5 wt % MnO. This compositionis similar to the bulk composition of the oxide ore, as inferredfrom whole-rock data. This similarity strongly suggests thatthe ores formed from accumulation of titanomagnetite crystals,not from immiscible oxide melt as proposed in earlier studies.The occurrence of oxide ores in the lower parts of the Panzhihuaintrusion is best explained by settling and sorting of densetitanomagnetite in the ferrogabbroic parental magma. This magmamust have crystallized Fe–Ti oxides relatively early andabundantly, and is likely to have been enriched in Fe and Tibut poor in SiO₂. These features are consistent with fractionationof mantle-derived melts under relatively high pressures (10kbar), followed by emplacement of the residual magma at 5 kbar.This study provides definitive field and geochemical evidencethat Fe–Ti oxide ores can form by accumulation in ferrogabbro.We suggest that many other massive Fe–Ti oxide depositsmay have formed in a similar fashion and that high concentrationsof phosphorus or carbon, or periodic fluctuation of fO₂ in themagma, are of secondary importance in ore formation. KEY WORDS: ELIP; Fe–Ti oxide ore; layered intrusion; Panzhihua; QUILF     

Sorption of Ni on Na-rectorite studied by batch and spectroscopy methods

Sorption of Ni²⁺ on Na-rectorite as a function of contact time, temperature, pH and fulvic acid (FA)/humic acid (HA) was studied under ambient conditions. A pseudo-second-order rate equation was used to simulate the kinetic sorption. The removal of Ni²⁺ increased with increasing pH. The presence of FA/HA enhanced the sorption of Ni²⁺ at low pH values, whereas no drastic effect of FA/HA on Ni²⁺ uptake to rectorite was found at high pH values. The diffuse layer model (DLM) fitted the experimental data of Ni²⁺ sorption in the absence and presence of FA/HA very well with the aid of FITEQL 3.2. The Langmuir, Freundlich and Dubinin–Radushkevich (D–R) models were used to simulate the sorption isotherms of Ni²⁺ at different temperatures. The thermodynamic data (ΔH⁰, ΔS⁰, ΔG⁰) were calculated from the temperature dependent sorption isotherms and the results suggested that the sorption process of Ni²⁺ on rectorite was spontaneous and endothermic. The sorption and species of Ni²⁺ on rectorite in the presence and absence of FA/HA was also investigated and characterized by XPS. The spectroscopic analysis indicated no drastic structural changes of Na-rectorite and the sorption of Ni²⁺ mainly occurred on the surface and at the edge position of Na-rectorite.