Quantifying urban land cover change between 2001 and 2006 in the Gulf of Mexico region

We estimated urbanization rates (2001–2006) in the Gulf of Mexico region using the National Land Cover Database (NLCD) 2001 and 2006 impervious surface products. An improved method was used to update the NLCD impervious surface product in 2006 and associated land cover transition between 2001 and 2006. Our estimation reveals that impervious surface increased 416 km² with a growth rate of 5.8% between 2001 and 2006. Approximately 1110.1 km² of non-urban lands were converted into urban land, resulting in a 3.2% increase in the region. Hay/pasture, woody wetland, and evergreen forest represented the three most common land cover classes that transitioned to urban. Among these land cover transitions, more than 50% of the urbanization occurred within 50 km of the coast. Our analysis shows that the close-to-coast land cover transition trend, especially within 10 km off the coast, potentially imposes substantial long-term impacts on regional landscape and ecological conditions.     

2D equivalent linear site effect simulation: example applications to two deep valleys

In the framework of the Sismovalp European project, an equivalent linear 2D code was developed to compute the response of a valley to SH waves, using the discrete wave-number method proposed by Aki and Larner (Aki K, Larner KL (1970) J Geophys Res 75:5). To overcome the frequency upper bound limitation, the Aki and Larner’s method is combined with a one-dimensional computation using a classical multi-layer method (Aki K, Richards PG (1980) Quantitative Seismology: Theory and Methods, vols. 1 & 2. W.H. Freeman & Co, San Francisco). The so-called “Aki–Larner extended method” is associated to an iterative algorithm, as proposed by Seed and Idriss (Seed HB, Idriss IM (1969) Report No. EERC 70–10, Earthquake Research Center, University of California, Berkeley, California) which accounts for the modulus and damping degradation using a linear visco-elastic model. A comparison of the results in the linear and the equivalent linear cases, for a magnitude 6.0 earthquake, shows that the account for the equivalent linear behaviour of the soil significantly reduces the amplification level, especially at frequencies higher than the fundamental resonance frequency of the site. In the case of site effects or microzonation studies devoted to produce design spectra for engineering structures, this can have a major impact on the associated results and costs, depending on the frequency of interest for the considered structure. As a first application of the developed technique, 2D equivalent linear Aki–Larner computations are used to perform the seismic microzonation study of the upper Rhone valley, in the Visp area (Switzerland), a typical 2D alpine valley. These investigations made it possible to determine site specific spectra, associated with different zones, to be used instead of the code spectra that do not take into account the local 2D amplification.     

Isotope characterization of lead in galena from ore deposits of the Aysén Region, southern Chile

Lead isotope analyses of galena from five ore deposits and six prospects in the Aysén region of southern Chile are reported. Most of the deposits are either low sulfidation epithermal gold–silver veins or skarn and manto deposits; the majority are either suspected to be, or dated as, Late Jurassic to mid-Cretaceous. Galena lead isotope data for most of the deposits from southern Chile cluster near the “orogene” within a “plumbotectonic” model framework. Average values (²⁰⁶Pb/²⁰⁴Pb=18.53, ²⁰⁷Pb/²⁰⁴Pb=15.63, and ²⁰⁸Pb/²⁰⁴Pb=38.50) are near Jurassic to Cretaceous model ages on the “orogene” curve of Zartman and Doe (1981) and the second-stage curve of Stacey and Kramers (1975) on a ²⁰⁶Pb/²⁰⁴Pb versus ²⁰⁷Pb/²⁰⁴Pb plot. These model ages are compatible with absolute ages as currently known. The elongate trends in the general cluster indicate mainly an orogenic model fit, suggesting variable mixing of lead from different sources, mainly model upper crust and lesser model mantle and lower crust reservoirs. Galena lead associated with one deposit (El Faldeo) is relatively radiogenic, and lies near a Jurassic age on the “upper crustal” curve of Zartman and Doe (1981), which is compatible with the Ar/Ar age of the deposit. Galena lead isotope clusters define three main groups of deposits. These three groups appear to be related to three mineralizing events, dated by K–Ar and Ar/Ar, in the Late Jurassic (group 3), and in the Early and mid-Cretaceous (groups 1 and 2 respectively). Averages for group 1, the northern group including El Toqui and Katerfeld, are ²⁰⁶Pb/²⁰⁴Pb=18.51, ²⁰⁷Pb/²⁰⁴Pb=15.62, ²⁰⁸Pb/²⁰⁴Pb=38.48. Averages for group 2, the southern group with Fachinal and Mina Silva, are ²⁰⁶Pb/²⁰⁴Pb = 18.56, ²⁰⁷Pb/²⁰⁴Pb=15.63, ²⁰⁸Pb/²⁰⁴Pb=38.52. Averages for group 3, the southernmost group with the El Faldeo, Lago Chacabuco and Lago Cochrane prospects, are ²⁰⁶Pb/²⁰⁴Pb=18.83, ²⁰⁷Pb/²⁰⁴Pb=15.65, ²⁰⁸Pb/²⁰⁴Pb=38.63. The Cretaceous deposits (groups 1 and 2) contain orogene-type lead that becomes increasingly radiogenic southward. Lead from the Late Jurassic deposits (group 3) appears to reflect mixing of orogene lead with highly radiogenic lead. The observed linear array of lead in group 3 probably reflects mixing of orogene lead with highly radiogenic lead, which was likely extracted by selective leaching of mineralizing hydrothermal solutions from the metamorphic basement. Received: 10 July 1999 / Accepted: 15 July 2000     

Long term evolution of distant retrograde orbits in the Earth-Moon system

This work studies the evolution of several Distant Retrograde Orbits (DROs) of varying size in the Earth-Moon system over durations up to tens of millennia. This analysis is relevant for missions requiring a completely hands off, long duration quarantine orbit, such as a Mars Sample Return mission or the Asteroid Redirect Mission. Four DROs are selected from four stable size regions and are propagated for up to 30,000 years with an integrator that uses extended precision arithmetic techniques and a high fidelity dynamical model. The evolution of the orbit’s size, shape, orientation, period, out-of-plane amplitude, and Jacobi constant are tracked. It has been found that small DROs, with minor axis amplitudes of approximately 45,000 km or less decay in size and period largely due to the Moon’s solid tides. Larger DROs (62,000 km and up) are more influenced by the gravity of bodies external to the Earth-Moon system, and remain bound to the Moon for significantly less time.     

Use of Sonication for Enhanced Sampling of Attached Microbes from Groundwater Systems

The vast majority of microorganisms in aquifers live as biofilms on sediment surfaces, which presents significant challenges for sampling as only the suspended microbes will be sampled through normal pumping. The use of a down-well low frequency sonicator has been suggested as a method of detaching microbes from the biofilm and allowing rapid sampling of this community. We developed a portable, easy to use, low-frequency electric sonicator and evaluated its performance for a range of well depths (tested up to 42 m below ground level) and casing types. Three sonicators were characterized in laboratory experiments using a 1 m long tank filled with pea gravel. These included a commercially available pneumatic sonicator, a rotating flexible shaft sonicator, and the prototype electric sonicator. The electric sonicator detached between 56 and 74% of microbes grown on gravel-containing biobags at distances ranging between 2 and 50 cm from the sonicator. The field testing comprises of a total of 55 sampling events from 48 wells located in 4 regions throughout New Zealand. Pre- and post-sonication samples showed an average 33 times increase in bacterial counts. Microbial sequence data showed that the same classes are present in pre- and post-sonicated samples and only slight differences were seen in the proportions present. The sampling process was rapid and the significant increases in bacterial counts mean that microbial samples can be quickly obtained from wells, which permits more detailed analysis than previously possible.     

Modelling chemo‐hydro‐mechanical behaviour of unsaturated clays: a feasibility study

Effective capabilities of combined chemo‐elasto‐plastic and unsaturated soil models to simulate chemo‐hydro‐mechanical (CHM) behaviour of clays are examined in numerical simulations through selected boundary value problems. The objective is to investigate the feasibility of approaching such complex material behaviour numerically by combining two existing models. The chemo‐mechanical effects are described using the concept of chemical softening consisting of reduction of the pre‐consolidation pressure proposed originally by Hueckel (Can. Geotech. J. 1992; 29 :1071–1086; Int. J. Numer. Anal. Methods Geomech. 1997; 21 :43–72). An additional chemical softening mechanism is considered, consisting in a decrease of cohesion with an increase in contaminant concentration. The influence of partial saturation on the constitutive behaviour is modelled following Barcelona basic model (BBM) formulation (Géotech. 1990; 40 (3):405–430; Can. Geotech. J. 1992; 29 :1013–1032). The equilibrium equations combined with the CHM constitutive relations, and the governing equations for flow of fluids and contaminant transport, are solved numerically using finite element. The emphasis is laid on understanding the role that the individual chemical effects such as chemo‐elastic swelling, or chemo‐plastic consolidation, or finally, chemical loss of cohesion have in the overall response of the soil mass. The numerical problems analysed concern the chemical effects in response to wetting of a clay specimen with an organic liquid in rigid wall consolidometer, during biaxial loading up to failure, and in response to fresh water influx during tunnel excavation in swelling clay. Copyright © 2005 John Wiley & Sons, Ltd.     

Computation of amplification factor of earthquake ground motion for a local sedimentary structure

We present methodology of calculating acceleration and corresponding earthquake ground motion characteristics at a site of interest assuming acceleration at a reference site for two basic configurations. In one configuration we assume that the reference ground motion is not affected by the local structure beneath the site of interest. In the other configuration we assume that the reference ground motion is affected by the local structure. Consequently, the two configurations differ from each other by the presence of the reference site within the computational model. For each of the two configurations we assume two wavefield excitations: a vertical plane-wave incidence and a point double-couple source. We illustrate the methodology on the example of the Grenoble valley. The extensive investigation of effects of local surface sedimentary structures based on the developed methodology is presented in the accompanying article by Moczo et al. (Bull Earthq Eng, 2018) in this volume.     

Towards a microzonation of the Greater Beirut area: an instrumental approach combining earthquake and ambient vibration recordings

Lebanon is situated on the 1000 km long Levant transform fault that separates the Arabic from the African tectonic plates. In Lebanon, the Levant fault splits up into a set of ramifications that had, in the past, generated major destructive earthquakes causing a lot of destruction and thousands of casualties. The most devastating one was the 551 A.D. offshore earthquake that destroyed Beirut, the capital of Lebanon. This paper presents a site effect study in Beirut, aimed at proposing a framework for future microzonation works in the city. It includes two complementary parts. A 6-month, temporary seismological experiment was first conducted to estimate the site response at 10 sites sampling the main geological units of Beirut on the basis of local and regional earthquake recordings. This spatially sparse information was then complemented by a large number (615) of microtremor measurements covering the Beirut municipality and part of its suburbs with a 400 m dense grid. The recordings were analysed with the standard site-to-reference and horizontal-to-vertical spectral ratio methods for earthquake recordings, and the horizontal-to-vertical ratio for ambient noise recordings. Significant ground motion amplification effects (up to a factor of 8) are found in a few areas corresponding to recent deposits. The consistency between results from earthquake and microtremor recordings allows proposing a map of the resonance frequencies within the city and its suburbs, with frequencies ranging from 0.5 to 5 Hz for the deepest deposits, and 5–10 Hz for shallow areas. Finally, the results are discussed and a way to combine the results obtained from the temporary stations to the great number of recordings coming from the permanent Lebanese network is proposed.     

Site effects in Mexico City eight years after the September 1985 Michoacan earthquakes

The purpose of this paper is to take a comprehensive look at site effects in Mexico City for the 1985 Michoacan earthquake. We examine, successively, 1D and 2D models. For the latter, we consider in detail both large scale and small scale heterogeneities, using extensively the Aki-Larner wave propagation method, in the version given by Bard and Gariel. In particular, we make a critical review of the different explanations proposed for the large duration of strong ground motion in the lake zone. Our purpose is two-sided. We first outline the difference between what is well established and what remains still unexplained regarding the seismic response of Mexico City basin. On the other hand, we wish to make explicit the conditions that the proposed models require to explain strong motion duration. Our results allow us to qualify the models proposed to date and to point out what could be the experiments and the new data required to find a truly satisfactory explanation of strong ground motion at Mexico City.