Development of a 3D GIS and its application to karst areas

There is a growing interest in modeling and analyzing karst phenomena in three dimensions. This paper integrates geology, groundwater hydrology, geographic information system (GIS), database management system (DBMS), visualization and data mining to study karst features in Huaibei, China. The 3D geo-objects retrieved from the karst area are analyzed and mapped into different abstract levels. The spatial relationships among the objects are constructed by a dual-linker. The shapes of the 3D objects and the topological models with attributes are stored and maintained in the DBMS. Spatial analysis was then used to integrate the data in the DBMS and the 3D model to form a virtual reality (VR) to provide analytical functions such as distribution analysis, correlation query, and probability assessment. The research successfully implements 3D modeling and analyses in the karst area, and meanwhile provides an efficient tool for government policy-makers to set out restrictions on water resource development in the area.     

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.     

The double solid reactant method: II. An application to the shallow groundwaters of the Porto Plain,Vulcano Island (Italy)

This paper presents an example of application of the double solid reactant method (DSRM) of Accornero and Marini (Environmental Geology, 2007a), an effective way for modeling the fate of several dissolved trace elements during water–rock interaction. The EQ3/6 software package was used for simulating the irreversible water–rock mass transfer accompanying the generation of the groundwaters of the Porto Plain shallow aquifer, starting from a degassed diluted crateric steam condensate. Reaction path modeling was performed in reaction progress mode and under closed-system conditions. The simulations assumed: (1) bulk dissolution (i.e., without any constraint on the kinetics of dissolution/precipitation reactions) of a single solid phase, a leucite-latitic glass, and (2) precipitation of amorphous silica, barite, alunite, jarosite, anhydrite, kaolinite, a solid mixture of smectites, fluorite, a solid mixture of hydroxides, illite-K, a solid mixture of saponites, a solid mixture of trigonal carbonates and a solid mixture of orthorhombic carbonates. Analytical concentrations of major chemical elements and several trace elements (Cr, Mn, Fe, Ni, Cu, Zn, As, Sr and Ba) in groundwaters were satisfactorily reproduced. In addition to these simulations, similar runs for a rhyolite, a latite and a trachyte permitted to calculate major oxide contents for the authigenic paragenesis which are comparable, to a first approximation, with the corresponding data measured for local altered rocks belonging to the silicic, advanced argillic and intermediate argillic alteration facies. The important role played by both the solid mixture of trigonal carbonates as sequestrator of Mn, Zn, Cu and Ni and the solid mixture of orthorhombic carbonates as scavenger of Sr and Ba is emphasized.

Flood hazard delineation combining geomorphological and hydrological methods: an example in the Northern Iberian Peninsula

Flood mapping requires the combination and integration of geomorphological and hydrological-hydraulic methods; however, despite this, there is very little scientific literature that compares and validates both methods. Two types of analysis are addressed in the present article. On the one hand, maps of flood plains have been elaborated using geomorphological evidence and historical flood data in the mountainous area of northwestern Spain, covering an area of more then 232 km² of floodplains. On the other hand, a hydrometeorological model has been developed (Clark semidistributed unit hydrograph) in the Sarria River basin (155 km², NW Spain). This basin is not gauged, hence the model was subjected to a goodness-of-fit test of its parameter (curve number) by means of Monte Carlo simulation. The peak flows obtained by means of the hydrological model were used for hydraulic modeling (one-phase, one-dimensional and steady flow) in a 4 km² urban stretch of the river bed. The delineation of surface areas affected by floods since 1918, as well as those analyzed subsequent to the geomorphological study, reveals a high degree of reliability in the delineation of the flooded areas with frequent recurrence intervals (<50 years). If we compare these flooded surface areas with the estimate obtained by the hydrological-hydraulic method we can see that the latter method overestimates the extent of the surface water by 144% for very frequent recurrence intervals (>10 years) and underestimates it as the recurrence interval increases, by up to 80% less floodplain for exceptional events (>500 years). Finally, a management map is put forth combining the most reliable results available by integrating both methods. Originally presented at the Sixth International Conference on Geomorphology.     

Numerical modeling of tide-surge interaction along Orissa coast of India

The Orissa coast of India is one of the most vulnerable regions of extreme sea levels associated with severe tropical cyclones. There was extensive loss of life and property due to the October 1999 super cyclone, which devastated large part of the Orissa coast. The shallow nature of the head bay, presence of a large number of deltas formed by major rivers of Orissa such as Mahanadi and Dhamra, and high tidal range are responsible for storm surge flooding in the region. Specifically, rising and falling tidal phases influence the height, duration, and arrival time of peak surge along the coast. The objective of the present study is to evaluate the tide-surge interaction during the 1999 Orissa cyclone by using nonlinear vertically integrated numerical models. The pure tidal solution for the head bay region of the Bay of Bengal provides the initial condition for the fine resolution nested grid Orissa model. However, the feedback from the Orissa model does not affect the head bay model as the study provides a one-way interaction. Numerical experiments are performed to study the tide-surge interaction by considering various relative phases of the tidal waves with the surge-wave produced by 1999 Orissa cyclone. The comparison, although utilizing only the limited estimates of tidal data, appears adequate to assert that the principal features are reproduced correctly.     

EPR investigation of the methyl radical,the hydrogen atom and carbon oxide radicals in Maxixe-type beryl

The EPR spectra of Maxixe-type beryl contain a large number of overlapping signals. The angular dependence of the 1:3:3:1 signal typical for the CH₃ radical shows that this radical is located at the center of the channel cavity with its symmetry axis parallel to the crystal c-axis and is rotating around this axis. Its EPR spectrum is axially symmetric with g _// = 2.00263, g _⊥ = 2.00249 and A_// = 2.288 mT, A_⊥ = 2.256 mT. These anisotropies have the opposite signs of those found for surface-adsorbed methyl radicals. Hydrogen atoms are located at position 2a at the center of the beryl cavity and the EPR parameters of the narrow doublet signal are A ₀ = 1,407 MHz and g = 2.00230. Another doublet signal, which is broader and has axial symmetry with g _// = 2.00265, g _⊥ = 2.00625 and A_// = 0.895 mT, A_⊥ = 0.885 mT, could come from a HCO₃ radical. One narrow and easily saturated signal with g _// = 2.00227 and g _⊥ = 2.00386 is interpreted to arise from a carbon monoxide radical in the beryl channel, oriented with its axis parallel to the crystal c-axis. Additional weak doublet lines, which have similar g values as the carbon monoxide radical, are created by nearby hydrogens. A powder spectrum with g _// = 2.0017 and g _⊥ = 2.0004 appears upon UV irradiation of the single crystal and is easily saturated. This spectrum is interpreted to arise from a carbon dioxide radical, which rotates around its symmetry axis.     

Influence of Structural Non-Stationarity of Surface Roughness on Morphological Characterization and Mechanical Deformation of Rock Joints

Summary Structural non-stationarity of surface roughness affects accurate morphological characterization as well as mechanical behaviour of rock joints at the laboratory scale using samples with a size below the stationarity threshold. In this paper, the effect of structural non-stationarity of surface roughness is investigated by studying the scale dependence of surface roughness and mechanical behaviour of rock joints. The results show that the structural non-stationarity mainly affects the accurate characterization of the surface roughness of the fracture samples. It also controls the amount and location of the contact areas during shear tests, which in turn affects the mechanical properties and asperity degradation of the samples. It is concluded that for accurate determination of the morphological and mechanical properties of rock joints at laboratory and field scales, samples with size equal to or larger than the stationarity threshold are required. Author’s address: Nader Fardin, Rock Mechanics Group, Department of Mining Engineering, Faculty of Engineering, University of Tehran, P.O. Box: 11365/4563, Tehran, Iran     

Ground motion estimation during the Kashmir earthquake of 8th October 2005

In this article, analytical methods have been used to estimate ground motion during the 8 October 2005, Kashmir earthquake. Peak ground acceleration (PGA) values at several stations in the epicentral region have been estimated by empirical analytical source mechanism models. As an alternate analysis, PGA estimates have also been obtained using the stochastic finite fault seismological model. The estimated PGAs are compared with that obtained from damage values. A PGA contour map in the near-source region is provided. It is found that very near to the epicenter, PGA would have reached more than 1 g. It is demonstrated that empirical analytical models can be effectively used to estimate ground motion due to rupture of active faults.     

TBM在不良地质洞条件下的施工技术

TBM掘进机以其快速、高效、安全、优质等优点越来越广泛地被应用于隧洞开挖施工中,尤其更适用于深埋超长隧洞,然而在不良地质洞段中TBM掘进缓慢,甚至有卡刀可能,反而不如钻爆法灵活,这就需要根据围岩性状采取特殊技术处理措施,辅以监控量测手段对支护方案进行验证、调整支护措施、修正设计参数等。结合辽宁大伙房输水隧洞工程,总结了在不良地质洞条件下的超前地质预报方法、不良地质段处理措施以及围岩变形监测方法。