Seismogenic basin and range and intra‐arc normal faulting in the central Mexican Volcanic Belt,Querétaro,México

North‐northwest normal faults intersect ENE normal faults in the vicinity of Querétaro City, in central México, affecting the Miocene–Pliocene northern‐central sector of the Mexican Volcanic Belt province. This intersection produced an orthogonal arrangement of grabens, half‐grabens and horsts that include the Querétaro graben. The NNW faults are part of the Taxco–San Miguel de Allende fault system, which is proposed here as part of the southernmost Basin and Range province in México. The ENE to E–W faults are part of the E–W oriented Chapala–Tula fault zone, which has been interpreted as an active intra‐arc fault system of the Mexican Volcanic Belt. Seventy‐four normal faults were mapped, of which the NNW faults are the largest and have the best morphological expression in the region. More numerous, although shorter, are the ENE faults. Total length of the ENE faults is greater than the total length of the NNW faults. Both sets are dominantly normal faults, indicating ENE extension for the NNW set and NNW extension for the ENE set. Field data indicate that displacement on the two fault sets has overlapped in time, as some NNW faults are younger than some ENE faults, which are supposed to be the younger ones. Seismicity in 1998 on a NNW fault indicates ENE active extension on the NNW faults. These observations support our interpretation that the northern Mexican Volcanic Belt lies on the boundary between the Basin and Range province, which is undergoing ENE extension, and the central Mexican Volcanic Belt province, which is undergoing northerly extension. The apparent overlap in space and time of displacements on the two fault sets reflects the difference in stress regime between the two provinces. Copyright © 2005 John Wiley & Sons, Ltd.     

Flexural and gravity modelling of the Mérida Andes and Barinas–Apure Basin, Western Venezuela

The kinematic evolution of the Barinas–Apure Basin and the southern Mérida Andes from Lower Miocene to the Present is numerically modelled using flexural isostatic theory and geophysical and geological data. Two published regional transects are used to build up a reference section, which is then used to constrain important parameters (e.g. shortenings and sedimentary thicknesses) for the flexural modelling. To control the location of the main fault system in the flexural model earthquake information is also used. The estimated flexural elastic thickness of the South American lithosphere beneath the Barinas–Apure Basin and the Mérida Andes Range is 25 km. The value for the final total shortening is 60 km. The flexural isostatic model shows that the Andean uplift has caused the South American lithosphere subsidence and the development of the Barinas–Apure Basin.In addition, gravity modelling was used to understand deep crustal features that could not be predicted by flexural theory. Consequently, the best-fit flexural model is used to build a gravity model across the Mérida Andes and the Barinas–Apure Basin preserving the best-controlled structural features from the flexural modelling (e.g. basin wavelength and depth) and slightly changing the main bodies density values and deep crustal structures. The final gravity model is intended to be representative of the major features affecting the gravity field in the study area. The predicted morphology in the lower crustal level of the final gravity model favours the hypothesis of a present delamination or megathrust of the Maracaibo crust over the South American Shield. This process would use the Conrad discontinuity as a main detachment surface within an incipient NW dipping continental subduction.     

Monitoring ground subsidence induced by salt mining in the city of Tuzla (Bosnia and Herzegovina)

Salt mining induced ground subsidence is a major hazard in the city of Tuzla (Northeastern sector of Bosnia and Herzegovina) and its surroundings since 1950, when solution mining of salt deposits by boreholes began. An analysis of the large (and never before processed) amount of topographical data collected during two periods: from 1956 to the Balkan War, and from 1992 to 2003 has been made. The analysis reveals a cumulative subsidence as high as 12 m during the whole period, causing damage to buildings and infrastructures within an area that includes a large portion of the historical town. Human-induced subsidence, (with rates up to 40 cm/year in the most developed area), has been investigated to recognize the areas affected by the sinking phenomenon and to produce a subsidence hazard. The time series of topographical observations have been enlarged by conducting new surveys in the urban area by modern space-geodesy methodologies, such as static relative GPS (Global Positioning System) and high resolution satellite imageries. The GPS monitoring started in 2004 and detected a decrease in the subsidence rates to 20 cm/year related to the reduction of salt exploitation. There is close correlation between the average subsidence rate and the annual amount of salt extracted.     

Hydrochemical and isotopical tracers in the lacustrine aquifer of the Cerro Prieto area, Baja California, México

A series of patterns in the chemical and isotopic composition of the groundwater of the Cerro Prieto area have been recognized. These patterns were used to define four hydrochemical and isotopical zones in the area. One of these areas located at the east side was found to be the source of recharge, while another in the central part is where an older aquifer exist in lacustrine deposits. The low permeability of the area causes the aquifer to act as an aquitard, which at least is partly surrounded by an interface zone where mixing of waters occurs. Tritium and radiocarbon method revealed that the aquitard is an isolate system whose residence time is of 32 ka.     

Construction of three-phase data to model multiphase flow in porous media: Comparing an optimization approach to the finite element approach

Multiphase flow modelling is a major issue in the assessment of groundwater pollution. Three-phase flows are commonly governed by mathematical models that associate a pressure equation with two saturation equations. These equations involve a number of secondary variables that reflect the fluid behaviour in a porous medium. To improve the computational efficiency of multiphase flow simulators, several simplified reformulations of three-phase flow equations have been proposed. However, they require the construction of new secondary variables adapted to the reformulated flow equations. In this article, two different approaches are compared to quantify these variables. A numerical example is given for a typical fine sand.     

Paleosubsidence and active subsidence due to evaporite dissolution in the Zaragoza area (Huerva River valley, NE Spain): processes, spatial distribution and protection measures for transport routes

The lowest 17-km long reach of the Huerva River valley, down to its confluence with the Ebro River in Zaragoza city, flows across salt-bearing evaporites of the Ebro Tertiary Basin (NE Spain). Upstream, the horizontally lying Miocene evaporites are interfingered with non-soluble distal alluvial fan facies (shales and sandstones). The proportion of soluble facies in the Huerva River valley increases in a downstream direction towards the basin depocenter. On the basis of the type and magnitude of the paleosubsidence features, the valley has been divided into four reaches. Along reach I, undeformed terrace deposits less than 4 m thick rest on insoluble detrital bedrock. In reaches II and III, dissolution at the alluvium–bedrock boundary has generated local thickening, deformation and paleocollapse structures, which only affect the alluvial mantle. In reach IV, terrace deposits thicken to over 60 m resulting from a large-scale synsedimentary subsidence. In this sector, subsidence locally affects to both the alluvium and the underlying bedrock. This indicates that dissolution acts at the rockhead beneath the alluvial cover (alluvial karst) and within the evaporitic substratum (interstratal karst). The development of an intraevaporitic karst in reach IV is attributed to gypsum and salt dissolution. Irregular terrace gravel bodies (gravel pockets) embedded in a fine-grained matrix associated with paleocollapse structures have been interpreted as liquefaction–fluidization structures resulting from ground acceleration and suction induced by catastrophic collapses. Subsidence is currently active in the region affecting areas with a thin alluvial cover in reaches III and IV. The low subsidence activity in most of Zaragoza city is explained by the presence of thickened (around 50 m) and indurated alluvial deposits. In the surrounding area, numerous buildings in Cadrete and Santa Fe villages have been severely damaged by subsidence. Natural and human-induced subsidence favours the development of slope movements in the gypsum scarp overlooking Cadrete village. Several transport routes including the Imperial Canal (irrigation canal) and the recently completed Madrid–Barcelona high-speed railway are affected by human-induced sinkholes. The paleocollapse structures exposed in the trenches of this railway and a ring road under construction point to hazardous locations underlain by cavities and collapse structures where special protection measures should be applied. Rigid structures are recommended beneath the high-speed railway with sufficient strength to span the larger sinkholes with no deformation. Electronic monitoring devices linked to a warning system can detect subtle subsidence-induced deformations in carriageways or railways. This research demonstrates that the study of the paleokarst helps to understand the processes involved in the present-day subsidence phenomena and their general spatial distribution.     

Effects of artificial and natural recharge on chemical equilibrium in the Cerro Prieto reservoirs, Baja California, México

Chemical equilibrium modeling studies of two aquifers in the western part of the Cerro Prieto (CP-I) geothermal field were performed to obtain the activities of the main ions and complexes by using the EQQYAC program. The results provided a pH of 6.44 for the deep aquifer, which is slightly higher than the value of 5.5 found for the shallower reservoir. In the northern part of the field, the reservoir pH has decreased up to 1.5 units compared to that modeled for 1990. Variations of pH with respect to time were also found to occur in the southern part of the field. In spite of the large changes in hydrodynamic conditions occurring in the reservoir due to exploitation, the equilibrium between both K- and Na-feldspars and between K-feldspar and K-mica seems to control the chemical composition of the fluids, in both reservoirs.     

有限单元法与有限差分法在MT一维正演模拟中的对比分析

首先从电磁场所满足的麦克斯韦方程组出发,介绍了大地电磁测深正演的基本理论,并针对一维大地电磁模型加以讨论.运用有限单元法及有限差分法分别推导了大地电磁测深一维正演算法,并运用Matlab 7.0软件编写了相应的程序.为了检验这两种一维正演算法的准确性,设计了均匀半空间模型和层状介质模型,并给出了由本文程序、解析解得到的相应结果和图件,从而对正演结果进行对比分析.结果表明:两种方法的正演结果均真实地反映了模型的地电参数.     

ULF geomagnetic anomalies of possible seismogenic origin observed at Teoloyucan station, México, in 1999–2001: Intermediate and short-time analysis

The analysis of ULF geomagnetic field measured at Teoloyucan station (Central Mexico, 11′35.735″ W, 19 44′45.100″ N, 2280 m height) is presented in an intermediate (± 15 days) and short time scale (the day of the EQ occurrence) in relation to 7 major earthquakes occurred in Mexico in 1999–2001. Local changes in the fractal dynamics of the magnetic field are found to be important: a pronounced fall of the fractal index is frequently observed prior to the main shock. The study of the ULF resonant structure recently discovered in the frequencies f_R1 = 10.2−11.1 mHz and f_R2 = 13.6−14.5 mHz reveals changes in their character probably related to the processes of the earthquakes preparation. The success of the observation of the mentioned anomalies (specially the fractal index decrease) strongly depends on how close is the station from the epicenter, and what is the magnitude of the earthquake.     

Distribution of orogenic gold deposits in relation to fault zones and gravity gradients: targeting tools applied to the Eastern Goldfields, Yilgarn Craton, Western Australia

A quantitative spatial analysis of mineral deposit distributions in relation to their proximity to a variety of structural elements is used to define parameters that can influence metal endowment, deposit location and the resource potential of a region. Using orogenic gold deposits as an example, geostatistical techniques are applied in a geographic-information-systems-based regional-scale analysis in the high-data-density Yilgarn Craton of Western Australia. Metal endowment (gold production and gold ‘rank’ per square kilometer) is measured in incremental buffer regions created in relation to vector lines, such as faults. The greatest metal tonnages are related to intersections of major faults with regional anticlines and to fault jogs, particularly those of dilatant geometry. Using fault length in parameter search, there is a strong association between crustal-scale shear zones/faults and deposits. Nonetheless, it is the small-scale faults that are marginal or peripheral to the larger-scale features that are more prospective. Gravity gradients (depicted as multiscale edges or gravity ‘worms’) show a clear association to faults that host gold deposits. Long wavelength/long strikelength edges, interpreted as dominantly fault-related, have greater metal endowment and provide a first-order area selection filter for exploration, particularly in areas of poor exposure. Statistical analysis of fault, fold and gravity gradient patterns mainly affirms empirical exploration criteria for orogenic gold deposits, such as associations with crustal-scale faults, anticlinal hinge zones, dilational jogs, elevated fault roughness, strong rheological contrasts and medium metamorphic grade rocks. The presence and concurrence of these parameters determine the metallogenic endowment of a given fault system and segments within the system. By quantifying such parameters, the search area for exploration can be significantly reduced by an order of magnitude, while increasing the chance of discovery.     

A non‐linear coupled finite element–boundary element model for the prediction of vibrations due to vibratory and impact pile driving

This paper presents a non‐linear coupled finite element–boundary element approach for the prediction of free field vibrations due to vibratory and impact pile driving. Both the non‐linear constitutive behavior of the soil in the vicinity of the pile and the dynamic interaction between the pile and the soil are accounted for. A subdomain approach is used, defining a generalized structure consisting of the pile and a bounded region of soil around the pile, and an unbounded exterior linear soil domain. The soil around the pile may exhibit non‐linear constitutive behavior and is modelled with a time‐domain finite element method. The dynamic stiffness matrix of the exterior unbounded soil domain is calculated using a boundary element formulation in the frequency domain based on a limited number of modes defined on the interface between the generalized structure and the unbounded soil. The soil–structure interaction forces are evaluated as a convolution of the displacement history and the soil flexibility matrices, which are obtained by an inverse Fourier transformation from the frequency to the time domain. This results in a hybrid frequency–time domain formulation of the non‐linear dynamic soil–structure interaction problem, which is solved in the time domain using Newmark's time integration method; the interaction force time history is evaluated using the θ‐scheme in order to obtain stable solutions. The proposed hybrid formulation is validated for linear problems of vibratory and impact pile driving, showing very good agreement with the results obtained with a frequency‐domain solution. Linear predictions, however, overestimate the free field peak particle velocities as observed in reported field experiments during vibratory and impact pile driving at comparable levels of the transferred energy. This is mainly due to energy dissipation related to plastic deformations in the soil around the pile. Ground vibrations due to vibratory and impact pile driving are, therefore, also computed with a non‐linear model where the soil is modelled as an isotropic elastic, perfectly plastic solid, which yields according to the Drucker–Prager failure criterion. This results in lower predicted free field vibrations with respect to linear predictions, which are also in much better agreement with experimental results recorded during vibratory and impact pile driving. Copyright © 2008 John Wiley & Sons, Ltd.     

A computational framework to predict the water‐leakage pressure of segmental joints in underwater shield tunnels using an advanced finite element method

The construction of shield tunnels under riverbeds and seabeds has considerably increased over the past decades. Due to the ultra‐high water head, water leakage through tunnel joints is a major concern during a tunnel's service life. One practical solution to prevent groundwater penetration is to implement ethylene‐propylene‐diene‐monomer gaskets at the segmental joints. However, numerical simulation of fluid pressure penetration into rubber materials remains a challenging problem in computational mechanics. Severe mesh distortions can occur due to large deformation. Consequently, a convergent solution is difficult to achieve. This paper presents an Abaqus‐based numerical framework to solve the previously mentioned problem using the implicit finite element solver. The key aspects of this framework are twofold: (1) a remesh and re‐map algorithm to overcome the excessive mesh distortion, and (2) simulation of fluid penetration into the contact interface of the gaskets to reproduce the water‐leakage process at the tunnel joints. The proposed framework is first tested to simulate the gasket‐in‐groove mechanical behavior and is then validated using experimental data and the solution produced by an explicit finite element solver. The developed framework is then adopted to predict the water‐leakage pressure at gasketed tunnel joints to illustrate the practical applications. Finally, the numerical results are compared with experimental data to demonstrate the accuracy and robustness of the proposed method and confirm its superiority and effectiveness over existing methods. This novel method can be used by tunnel designers to analyze and estimate the waterproof behavior of gasketed joints in shield tunnels without performing extensive experimental testing works.     

Identification of the origin of salinization in groundwater using multivariate statistical analysis and geochemical modeling: a case study of Kaohsiung,Southwest Taiwan

The study of brine aquifers in southern Taiwan is highly complicated by hybrid geochemical reactions, which obscure important geochemical information. Using multivariate analysis on major and minor ion compositions normalized by Cl⁻ content, chemical constituents were combined into two principal components representing brine mixing and mineral precipitation. Comparing to multivariate analysis on the original data, this procedure reveals more geochemical information. It demonstrates that the brine groundwater of the region is primarily composed of highly evaporated seawater. The evaporation ratio is >70%; a point at which calcite, dolomite and gypsum precipitate. Oxygen and hydrogen isotopic compositions confirm this inference; and further, geochemical modeling quantitatively determined the evaporation ratio to be about 85%. Natural boron contamination is a consequence of brine groundwater. Two evolutionary trends in the plotting of the Cl/B ratio versus Cl⁻ can be identified: (1) Cl/B ratio decreases with boron being released from clay minerals when brine aquifers are flushed with freshwater; and (2) Cl/B ratio increases when seawater of a high Cl/B ratio infiltrates coastal aquifers.     

Constraints on element mobility associated with the conversion of granulite to eclogite along fractures in an anorthositic complex on Holsnøy, Norway

On Holsnøy, an island off the coast of Western Norway, an anorthositic complex metamorphosed to granulite facies was partially overprinted by a later eclogite facies metamorphism. Eclogite facies rocks (containing omphacite, garnet, kyanite and hydrous phases such as mica and zoisite) occur in shear zones of various scales and adjacent to veins. Previous studies of shear zones on Holsnøy reported evidence for substantial element mobility (Jamtveit et al ., 1990; Mattey et al ., 1994). In this work, we compare chemical compositions of granulite and its undeformed eclogitized equivalent adjacent to veins in locations where a single band of granulite can be traced and sampled as it approaches the vein. This tracing is crucial because the pre-granulite rocks cover a substantial compositional range, indicative of a petrologically variable protolith consisting of anorthosite, gabbro and jotunite. We analysed multiple core samples collected across nine separate granulite-eclogite transition zones located at veins in anorthositic, jotunitic and gabbroic protoliths for major and trace elements. For each transition, no compositional difference between the average granulite and average eclogite composition was found at the 90% confidence level except for LOI (loss on ignition), which was consistently significantly higher in the eclogite samples. Although not significant at the 90% confidence level for any single traverse, the average eclogite concentrations of SiO₂ , Na₂O, Cs, As and Br exceed the average granulite concentrations for eight or all nine of the traverses. For most traverses, statistical analysis of the data limits any gain of SiO₂ in the eclogites to no more than a few relative per cent. Other than the introduction of volatile substances, presumably an H₂O-rich fluid, eclogitization associated with vein formation was essentially isochemical.     

Seismic and geological structure of the crust in the transition from Baltica to Palaeozoic Europe in SE Poland—CELEBRATION 2000 experiment, profile CEL02

The CELEBRATION 2000 together with the earlier POLONAISE'97 deep seismic sounding experiments was aimed at the recognition of crustal structure in the border zone between the Precambrian East European Craton (Baltica) and Palaeozoic Europe. The CEL02 profile of the CELEBRATION family is a 400-km long SW–NE transect, running in Poland from the Upper Silesia Block (USB), across the Małopolska Block (MB) and the Trans-European Suture Zone (TESZ) to the East European Craton (EEC). The structure along CEL02 was interpreted using both 2D tomography and forward ray-tracing techniques as well as 2D gravity modelling.The crustal thickness along CEL02 varies from 32–35 km in the USB to 45–47 km beneath the TESZ and the EEC. The USB is a clearly distinctive crustal block with the characteristic high velocity lower crust (7.1–7.2 km/s), interpreted as a fragment of Gondwana. The Kraków–Lubliniec Fault is a terrane boundary produced by soft docking of the USB with the MB. The Małopolska crust fundamentally differs from the USB and has a strong connection with Baltica. It is a transitional, 150- to 200-km wide unit composed of the extended Baltican lower crust and the overlying low velocity (5.15–5.9 km/s) Neoproterozoic metasediments in the up to 18-km thick upper crust. The Łysogóry Unit has its crustal structure identical with that of Małopolska, thus it is connected with Baltica and cannot be interpreted as a Gondwana-derived terrane. Higher velocity and density bodies found below the Mazovia–Lublin Graben at a depth of 12 km and at the base of the lower crust, might be a result of mantle-derived mafic intrusions accompanying the extension of Baltica. By the preliminary 2D gravity modelling, we have reconfirmed the need for considering the increased TESZ mantle density in comparison to the EEC and USB mantle.     

The role of small‐scale fold and fault development in seismogenic zones: example of the western Huércal‐Overa Basin (eastern Betic Cordillera,Spain)

The NW–SE shortening between the African and the Eurasian plates is accommodated in the eastern Betic Cordillera along a broad area that includes large N‐vergent folds and kilometric NE–SW sinistral faults with related seismicity. We have selected the best exposed small‐scale tectonic structures located in the western Huércal‐Overa Basin (Betic Cordillera) to discuss the seismotectonic implications of such structures usually developed in seismogenic zones. Subvertical ESE–WNW pure dextral faults and E–W to ENE–ESW dextral‐reverse faults and folds deform the Quaternary sediments. The La Molata structure is the most impressive example, including dextral ESE–WNW Neogene faults, active southward‐dipping reverse faults and associated ENE–WSW folds. A molar M1 assigned to Mimomys savini allows for precise dating of the folded sediments (0.95–0.83 Ma). Strain rates calculated across this structure give ～0.006 mm a^?1 horizontal shortening from the Middle Pleistocene up until now. The widespread active deformations on small‐scale structures contribute to elastic energy dissipation around the large seismogenic zones of the eastern Betics, decreasing the seismic hazard of major fault zones. Copyright © 2009 John Wiley & Sons, Ltd.     

Finite element modelling of elastic intraplate stresses due to heterogeneities in crustal density and mechanical properties for the Jabalpur earthquake region,central India

Deep lower crustal intraplate earthquakes are infrequent and the mechanism of their occurrence is not well understood. The Narmada-Son-lineament region in central India has experienced two such events, the 1938 Satpura earthquake and the 1997 Jabalpur earthquake, having a focal depth of more than 35 km. We have estimated elastic stresses due to the crustal density and mechanical properties heterogeneities along the Hirapur-Mandla profile passing through the Jabalpur earthquake region to analyse conditions suitable for the concentration of shear stresses in the hypocentral region of this earthquake. Elastic stresses have been computed by a finite element method for a range of material parameters. The results indicate that the shear stresses generated by the density heterogeneities alone are not able to locally enhance the stress concentration in the hypocentral region. The role of mechanical properties of various crustal layers is important in achieving this localization of stresses. Among a range of material parameters analysed, the model with a mechanically strong lower crust overlying a relatively weak sub-Moho layer is able to enhance the stress concentration in the hypocentral region, implying a weaker mantle in comparison to the lower crust for this region of central India.     

Early biodégradation of ligneous organic materials and its relation to ore deposition in the tréves zn-pb ore body (gard, france)

Previous petrographical and geochemical studies on the organic content of mineralized and non-mineralized carbonate series in the Trèves Zn-Pb stratiform deposit led to the conclusion that ore deposition occurred during the diagenetic history of the sediments, in a mixed fresh/hydrothermal system (Gauthier and Disnar, 1984; Gauthier et al., 1985a). These circulations of fluids are probably responsible for the dolomitization of the original limestones and for the deposition of the ore through microbial activity (Gauthier, 1984; Gauthier el al., 1985b).The occurrence of this activity is supported by the formation of two kinds of organic components which are attributable to this uncommon bio-alteration process: alphatic-rich bitumen-like products and aromatic polyethers. The chemical composition of these two types of compounds suggests that they were produced from ligneous debris by (sulfate-reducing?) bacteria.Other indications of the occurrence of a biological activity have been found in the nature and relative abundance of microbial amino compounds. However, these constituents are thought to have been introduced only recently into the sediments whereas other evidence supports the idea that ore deposition occurred before diagenesis.     

High resolution mapping,description and interpretation of the foreshore in front of the great undercliff landslides; the Plateau,Dowlands and Bindon within the Axmouth to Lyme Regis Undercliffs National Nature Reserve,south Devon coast,England, and how that might inform understanding of the mechanisms of failure

The Axmouth to Lyme Regis Undercliff National Nature Reserve contains three highly unusual landslides; the Plateau, Dowlands and Bindon. The failure surface of the Plateau landslide is clearly within the Triassic strata rather than on the Jurassic-Cretaceous unconformity, which is the primary failure surface elsewhere along this coast. The nature of the famous Bindon Landslide of 1839, which formed Goat Island and the Chasm, has been the subject of much debate since it was first described by Buckland and Conybeare in 1840. While it is broadly agreed to be a translational failure, the stratigraphic location of the failure surface is still not agreed. This paper examines the Plateau slip in detail and attempts to apply what can be seen there to Dowlands and Bindon. While Pitts undertook a highly detailed mapping exercise of the landslides throughout the undercliff itself, he did not map the foreshore in any detail. Therefore an Unmanned Aerial Vehicle or drone was used to generate a high resolution photomontage of the foreshore in front of all three landslides to enable detailed mapping of the features contained within it. The results are available on the web. The mapping does not resolve the question of the location of the failure surface of both Dowlands and Bindon but it provides further detail for interpretation.