A seismic source zone model for the seismic hazard assessment of the Italian territory

We designed a new seismic source model for Italy to be used as an input for country-wide probabilistic seismic hazard assessment (PSHA) in the frame of the compilation of a new national reference map.

We started off by reviewing existing models available for Italy and for other European countries, then discussed the main open issues in the current practice of seismogenic zoning.

The new model, termed ZS9, is largely based on data collected in the past 10 years, including historical earthquakes and instrumental seismicity, active faults and their seismogenic potential, and seismotectonic evidence from recent earthquakes. This information allowed us to propose new interpretations for poorly understood areas where the new data are in conflict with assumptions made in designing the previous and widely used model ZS4.

ZS9 is made out of 36 zones where earthquakes with M_w > = 5 are expected. It also assumes that earthquakes with M_w up to 5 may occur anywhere outside the seismogenic zones, although the associated probability is rather low. Special care was taken to ensure that each zone sampled a large enough number of earthquakes so that we could compute reliable earthquake production rates.

Although it was drawn following criteria that are standard practice in PSHA, ZS9 is also innovative in that every zone is characterised also by its mean seismogenic depth (the depth of the crustal volume that will presumably release future earthquakes) and predominant focal mechanism (their most likely rupture mechanism). These properties were determined using instrumental data, and only in a limited number of cases we resorted to geologic constraints and expert judgment to cope with lack of data or conflicting indications. These attributes allow ZS9 to be used with more accurate regionalized depth-dependent attenuation relations, and are ultimately expected to increase significantly the reliability of seismic hazard estimates.     

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.