Tectonic significance of an earthquake sequence in the Zacoalco half-graben, Jalisco, Mexico

We studied a sequence of small earthquakes that occurred during the months of April and May of 1997, in Jalisco, southwestern Mexico. The earthquakes were located along a set of active faults that form the Zacoalco half-graben (La Lima fault system), west of Lake Chapala, within the rift–rift–rift triple junction. A total of 33 events were located, with magnitudes ranging from 1.5 to 3.5, recorded by a portable array of broadband seismographs. We identified two groups of events: one corresponding to a shallow normal fault, synthetic to La Lima fault system, and another group associated with a deeper fault. The events that occurred on the synthetic fault show normal faulting oriented on a NW–SE plane, dipping shallowly towards the SW. The other group of mechanisms showed either a normal fault oriented NW–SE and dipping steeply to the NE, or a very shallow-dipping normal fault, dipping to the SW. Earthquake distribution and fault plane solutions suggest that the Zacoalco half-graben developed from blocks that rotate as slip occurs on listric faults. These mechanisms could represent the type of motion expected for larger earthquakes in the area, like the one that occurred in 1568.     

Tectonic implications of the earthquake swarm of 1997 in the Michoacan Triangle, Mexico

An earthquake swarm occurred during February and March 1997 in the vicinity of the Tancitaro Volcano, in the southern part of the tectonically complex Michoacan Triangle. A study of these events provides an opportunity to map the active faults in the area and to learn if the orientation and the sense of motion on these faults are consistent with the mapped faults and the alignment of cinder cones in the region. The foci of 230 earthquakes, which could be located, are distributed between 10 and 18 km depth, and show an alignment in, roughly, a NE direction. The focal mechanisms and seismic moments of the 27 best-recorded events were determined by waveform modeling of P and S waves. These mechanisms show two distinct patterns. More than 50% of the solutions are left-lateral strike–slip mechanisms with a normal component. The preferred fault plane strikes NE. Another group of events, probably caused by triggered seismicity on the Chapala–Oaxaca fault zone, shows left-lateral strike–slip mechanisms with a large-thrust component on NW-trending faults. S wave splitting shows 1–2.5% crustal-anisotropy. The direction of the anisotropy coincides with the NE alignment of events, and the preferred nodal plane. This is also the alignment of cinder cones, suggesting that preexisting fractures and cracks are responsible for the seismicity and anisotropic behavior of the crust. The resulting stress orientation, NE compression, is the one expected for the fore-arc region. We conclude that although Michoacan Triangle lies in the Trans-Mexican Volcanic Belt, it does not form part of this stress province where the stress orientation is NS extension.     

Numerical modeling of storm surges in the coast of Mozambique: the cases of tropical cyclones Bonita (1996) and Lisette (1997)

The coast of Mozambique is often affected by storms, particularly tropical cyclones during summer or sometimes midlatitude systems in the southern part. Storm surges combined with high freshwater discharge can drive huge coastal floods, affecting both urban and rural areas. To improve the knowledge about the impact of storm surges in the coast of Mozambique, this study presents the first attempt to model this phenomenon through the implementation of the Princeton Ocean Model (POM) in the Southwestern Indian Ocean domain (SWIO; 2–32°S, 28–85°E) using a regular grid with 1/6° of spatial resolution and 36 sigma levels. The simulation was performed for the period 1979–2010, and the most interesting events of surges were related to tropical cyclones Bonita (1996) and Lisette (1997) that occurred in the Mozambique Channel. The results showed that the model represented well the amplitude and phase of principal lunar and solar tidal constituents, as well as it captured the spatial pattern and magnitudes of SST with slight positive bias in summer and negative bias in winter months. In terms of SSH, the model underestimated the presence of mesoscale eddies, mainly in the Mozambique Channel. Our results also showed that the atmospheric sea level pressure had a significant contribution to storm heights during the landfall of the tropical cyclones Bonita (1996) and Lisette (1997) in the coast of Mozambique contributing with about 20 and 16% of the total surge height for each case, respectively, surpassing the contribution of the tide-surge nonlinear interactions by a factor of 2.     

Characterisation and multifaceted anisotropy assessment of Corvio sandstone for geological CO2 storage studies

We present a comprehensive characterisation of the physical, mineralogical, geomechanical, geophysical, and hydrodynamic properties of Corvio sandstone. This information, together with a detailed assessment of anisotropy, is needed to establish Corvio sandstone as a useful laboratory rock‐testing standard for well‐constrained studies of thermo–hydro–mechanical–chemical coupled phenomena associated with CO₂ storage practices and for geological reservoir studies in general. More than 200 core plugs of Corvio sandstone (38.1 and 50 mm diameters, 2:1 length‐to‐diameter ratio) were used in this characterisation study, with a rock porosity of 21.7 ± 1.2%, dry density 2036 ± 32 kg m^?3, and unconfined compressive and tensile strengths of 41 ± 3.28 and 2.3 ± 0.14 MPa, respectively. Geomechanical tests show that the rock behaves elastically between ～10 and ～18 MPa under unconfined conditions with associated Young's modulus and Poisson's ratio of 11.8 ± 2.8 GPa and 0.34 ± 0.01 GPa, respectively. Permeability abruptly decreases with confining pressure up to ～10 MPa and then stabilises at ～1 mD. Ultrasonic P‐ and S‐wave velocities vary from about 2.8–3.8 km s^?1 and 1.5–2.4 km s^?1, respectively, over confining and differential pressures between 0.1 and 35 MPa, allowing derivation of associated dynamic elastic moduli. Anisotropy was investigated using oriented core plugs for electrical resistivity, elastic wave velocity and attenuation, permeability, and tracer injection tests. Corvio sandstone shows weak transverse isotropy (symmetry axis normal to bedding) of <10% for velocity and <20% for attenuation.     

Evaluation of ASCE 7 equations for designing acceleration‐sensitive nonstructural components using data from instrumented buildings

This study uses instrumented buildings and models of code‐based designed buildings to validate the results of previous studies that highlighted the need to revise the ASCE 7 F_p equation for designing nonstructural components (NSCs) through utilizing oversimplified linear and nonlinear models. The evaluation of floor response spectra of a large number of instrumented buildings illustrates that, unlike the ASCE 7 approach, the in‐structure and the component amplification factors are a function of the ratio of NSC period to the supporting building modal periods, the ground motion intensity, and the NSC location. It is also shown that the recorded ground motions at the base of instrumented buildings in most cases are significantly lower than design earthquake (DE) ground motions. Because ASCE 7 is meant to provide demands at a DE level, for a more reliable evaluation of the F_p equation, 2 representative archetype buildings are designed based on the ASCE 7‐16 seismic provisions and exposed to various ground motion intensity levels (including those consistent with the ones experienced by instrumented buildings and the DE). Simulation results of the archetype buildings, consistent with previous numerical studies, illustrate the tendency of the ASCE 7 in‐structure amplification factor, [1 + 2(z/h)] , to significantly overestimate demands at all floor levels and the ASCE 7 limit of to in many cases underestimate the calculated NSC amplification factors. Furthermore, the product of these 2 amplification factors (that represents the normalized peak NSC acceleration) in some cases exceeds the ASCE 7 equation by a factor up to 1.50.     

SlugIn 1.0: A Free Tool for Automated Slug Test Analysis

The correct characterization of aquifer parameters is essential for water‐supply and water‐quality investigations. Slug tests are widely used for these purposes. While free software is available to interpret slug tests, some codes are not user‐friendly, or do not include a wide range of methods to interpret the results, or do not include automatic, inverse solutions to the test data. The private sector has also generated several good programs to interpret slug test data, but they are not free of charge. The computer program SlugIn 1.0 is available online for free download, and is demonstrated to aid in the analysis of slug tests to estimate hydraulic parameters. The program provides an easy‐to‐use Graphical User Interface. SlugIn 1.0 incorporates automated parameter estimation and facilitates the visualization of several interpretations of the same test. It incorporates solutions for confined and unconfined aquifers, partially penetrating wells, skin effects, shape factor, anisotropy, high hydraulic conductivity formations and the Mace test for large‐diameter wells. It is available in English and Spanish and can be downloaded from the web site of the Geological Survey of Spain. Two field examples are presented to illustrate how the software operates.     

Simulation of the effect of wind speedup in the formation of transverse dune fields

A computer simulation model for transverse‐dune‐field dynamics, corresponding to a uni‐directional wind regime, is developed. In a previous formulation, two distinct problems were found regarding the cross‐sectional dune shape, namely the erosion in the lee of dunes and the steepness of the windward slopes. The first problem is solved by introducing no erosion in shadow zones. The second issue is overcome by introducing a wind speedup (shear velocity increase) factor, which can be accounted for by adding a term to the original transport length, which is proportional to the surface height. By incorporating these features we are able to model dunes whose individual shape and collective patterns are similar to those observed in nature. Moreover we show how the introduction of a non‐linear shear‐velocity‐increase term leads to the reduction of dune height, and this may result in an equilibrium dune field configuration. This is thought to be because the non‐linear increase of the transport length makes the sand trapping efficiency lower than unity, even for higher dunes, so that the incoming and the outgoing sand flux are in balance. To fully describe the inter‐dune morphology more precise dynamics in the lee of the dune must be incorporated. Copyright © 2000 John Wiley & Sons, Ltd.     

Effects of the improper modeling of viscous damping on the first-mode and higher-mode dominated responses of base-isolated buildings

In many finite element platforms, a classical global damping matrix based on the elastic stiffness of the system (including isolators) is usually developed as part of the solution to the equations of motion of base-isolated buildings. The conducted analytical and numerical investigations illustrate that this approach can lead to the introduction of unintended damping to the first and higher vibration modes and the spurious suppression of the respective structural responses. A similar shortcoming might be observed even when a nonclassical damping model (ie, an assembly of the superstructure and isolation system damping sub-matrices) is used. For example, the use of Rayleigh damping approach to develop the superstructure damping sub-matrix can lead to the undesired addition of damping to the isolated mode arising from the mass-proportional component of the superstructure damping. On the other hand, the improper use of nonclassical stiffness-proportional damping (eg, determining the proportional damping coefficient, β_k , based on the first mode) can result in assigning significant damping to the higher-modes and the unintended mitigation of the higher-mode responses. Results show that a nonclassical stiffness-proportional model in which β_k is determined based on the second modal period of a base-isolated building can reasonably specify the intended damping to the higher modes without imparting undesirable damping to the first mode. The nonclassical stiffness-proportional damping can be introduced to the numerical model through explicit viscous damper elements attached between adjacent floors. In structural analysis software such as SAP2000®, the desired nonclassical damping can be also modeled through specifying damping solely to the superstructure material.     

A New Appraisal of Sri Lankan BB Zircon as a Reference Material for LA‐ICP‐MS U‐Pb Geochronology and Lu‐Hf Isotope Tracing

A potential zircon reference material (BB zircon) for laser ablation‐inductively coupled plasma‐mass spectrometry (LA‐ICP‐MS) U‐Pb geochronology and Hf isotope geochemistry is described. A batch of twenty zircon megacrysts (0.5–1.5 cm³) from Sri Lanka was studied. Within‐grain rare earth element (REE) compositions are largely homogeneous, albeit with some variation seen between fractured and homogeneous domains. Excluding fractured cathodoluminescence bright domains, the variation in U content for all analysed crystals ranged from 227 to 368 μg g^?1 and the average Th/U ratios were between 0.20 and 0.47. The Hf isotope composition (0.56–0.84 g/100 g Hf) is homogeneous within and between the grains – mean ¹⁷⁶Hf/¹⁷⁷Hf of 0.281674 ± 0.000018 (2s). The calculated alpha dose of 0.59 × 10¹⁸ g^?1 for a number of BB grains falls within the trend of previously studied, untreated zircon samples from Sri Lanka. Aliquots of the same crystal (analysed by ID‐TIMS in four different laboratories) gave consistent U‐Pb ages with excellent measurement reproducibility (0.1–0.4% RSD). Interlaboratory assessment (by LA‐ICP‐MS) from individual crystals returned results that are within uncertainty equivalent to the TIMS ages. Finally, we report on within‐ and between‐grain homogeneity of the oxygen isotope systematic of four BB crystals (13.16‰ VSMOW).