DIFFERENTIATING PATH AND SOURCE PROCESSES IN COMPLEX GEOLOGIC STRUCTURE THROUGH PASSIVE SEISMIC STUDIES: BERING GLACIER, ALASKA AND VILLARRICA VOLCANO, CHILE

Measuring shallow seismic sources provides a way to reveal processes that cannot be directly observed, but the correct interpretation and value of these signals depend on the ability to distinguish source from propagation effects. Furthermore, seismic signals produced by a resonating source can look almost identical to those produced by impulsive sources, but modified along the path. Distinguishing these two phenomena can be accomplished by examining the wavefield with small aperture arrays or by recording seismicity near to the source when possible. We examine source and path effects in two different environments: Bering Glacier, Alaska and Villarrica Volcano, Chile. Using three 3-element seismic arrays near the terminus of the Bering Glacier, we have identified and located both terminus calving and iceberg breakup events. We show that automated array analysis provided a robust way to locate icequake events using P waves. This analysis also showed that arrivals within the long-period codas were incoherent within the small aperture arrays, demonstrating that these codas previously attributed to crack resonance were in fact a result of a complicated path rather than a source effect. At Villarrica Volcano, seismometers deployed from near the vent to ~10 km revealed that a several cycle long-period source signal recorded at the vent appeared elongated in the far-field. We used data collected from the stations nearest to the vent to invert for the repetitive seismic source, and found it corresponded to a shallow force within the lava lake oriented N75°E and dipping 7° from horizontal. We also used this repetitive signal to search the data for additional seismic and infrasonic properties which included calculating seismic-acoustic delay times, volcano acoustic-seismic ratios and energies, event frequency, and real-time seismic amplitude measurements. These calculations revealed lava lake level and activity fluctuations consistent with lava lake level changes inferred from the persistent infrasonic tremor.

Seismic Serviceability Assessment of Water Distribution Systems

Water distribution systems are important for life saving facilities especially in the recovery after earthquakes. In this paper, a framework is discussed about seismic serviceability of water systems that includes the fragility evaluation of water sources of water distribution networks. Also, a case study is brought about the performance of a water system under different levels of seismic hazard. The seismic serviceability of a water supply system provided by EPANET is evaluated under various levels of seismic hazard. Basically, the assessment process is based on hydraulic analysis and Monte Carlo simulations, implemented with empirical fragility data provided by the American Lifeline Alliance (ALA, 2001) for both pipelines and water facilities. Represented by the Seismic Serviceability Index (Cornell University, 2008), the serviceability of the water distribution system is evaluated under each level of earthquakes with return periods of 72 years, 475 years, and 2475 years. The system serviceability under levels of earthquake hazard are compared with and without considering the seismic fragility of the water source. The results show that the seismic serviceability of the water system decreases with the growing of the return period of seismic hazard, and after considering the seismic fragility of the water source, the seismic serviceability decreases. The results reveal the importance of considering the seismic fragility of water sources, and the growing dependence of the system performance of water system on the seismic resilience of water source under severe earthquakes.

TIME LAPSE SEISMIC OBSERVATIONS AND EFFECTS OF RESERVOIR COMPRESSIBILITY AT TEAL SOUTH OIL FIELD

One of the original ocean-bottom time-lapse seismic studies was performed at the Teal South oil field in the Gulf of Mexico during the late 1990’s. This work reexamines some aspects of previous work using modern analysis techniques to provide improved quantitative interpretations. Using three-dimensional volume visualization of legacy data and the two phases of post-production time-lapse data, I provide additional insight into the fluid migration pathways and the pressure communication between different reservoirs, separated by faults. This work supports a conclusion from previous studies that production from one reservoir caused regional pressure decline that in turn resulted in liberation of gas from multiple surrounding unproduced reservoirs. I also provide an explanation for unusual time-lapse changes in amplitude-versus-offset (AVO) data related to the compaction of the producing reservoir which, in turn, changed an isotropic medium to an anisotropic medium. In the first part of this work, I examine regional changes in seismic response due to the production of oil and gas from one reservoir. The previous studies primarily used two post-production ocean-bottom surveys (Phase I and Phase II), and not the legacy streamer data, due to the unavailability of legacy prestack data and very different acquisition parameters. In order to incorporate the legacy data in the present study, all three poststack data sets were cross-equalized and examined using instantaneous amplitude and energy volumes. This approach appears quite effective and helps to suppress changes unrelated to production while emphasizing those large-amplitude changes that are related to production in this noisy (by current standards) suite of data. I examine the multiple data sets first by using the instantaneous amplitude and energy attributes, and then also examine specific apparent time-lapse changes through direct comparisons of seismic traces. In so doing, I identify time-delays that, when corrected for, indicate water encroachment at the base of the producing reservoir. I also identify specific sites of leakage from various unproduced reservoirs, the result of regional pressure blowdown as explained in previous studies; those earlier studies, however, were unable to identify direct evidence of fluid movement. Of particular interest is the identification of one site where oil apparently leaked from one reservoir into a “new” reservoir that did not originally contain oil, but was ideally suited as a trap for fluids leaking from the neighboring spill-point. With continued pressure drop, oil in the new reservoir increased as more oil entered into the reservoir and expanded, liberating gas from solution. Because of the limited volume available for oil and gas in that temporary trap, oil and gas also escaped from it into the surrounding formation. I also note that some of the reservoirs demonstrate time-lapse changes only in the “gas cap” and not in the oil zone, even though gas must be coming out of solution everywhere in the reservoir. This is explained by interplay between pore-fluid modulus reduction by gas saturation decrease and dry-frame modulus increase by frame stiffening. In the second part of this work, I examine various rock-physics models in an attempt to quantitatively account for frame-stiffening that results from reduced pore-fluid pressure in the producing reservoir, searching for a model that would predict the unusual AVO features observed in the time-lapse prestack and stacked data at Teal South. While several rock-physics models are successful at predicting the time-lapse response for initial production, most fail to match the observations for continued production between Phase I and Phase II. Because the reservoir was initially overpressured and unconsolidated, reservoir compaction was likely significant, and is probably accomplished largely by uniaxial strain in the vertical direction; this implies that an anisotropic model may be required. Using Walton’s model for anisotropic unconsolidated sand, I successfully model the time-lapse changes for all phases of production. This observation may be of interest for application to other unconsolidated overpressured reservoirs under production.

Stratigraphic analysis of lake level fluctuations in Lake Ohrid: an integration of high resolution hydro-acoustic data and sediment cores

Abstract. Ancient Lake Ohrid is a steep-sided, oligotrophic, karst lake that was tectonically formed most likely within the Pliocene and often referred to as a hotspot of endemic biodiversity. This study aims on tracing significant lake level fluctuations at Lake Ohrid using high-resolution acoustic data in combination with lithological, geochemical, and chronological information from two sediment cores recovered from sub-aquatic terrace levels at ca. 32 and 60m water depth. According to our data, significant lake level fluctuations with prominent lowstands of ca. 60 and 35m below the present water level occurred during Marine Isotope Stage (MIS) 6 and MIS 5, respectively. The effect of these lowstands on biodiversity in most coastal parts of the lake is negligible, due to only small changes in lake surface area, coastline, and habitat. In contrast, biodiversity in shallower areas was more severely affected due to disconnection of today sublacustrine springs from the main water body. Multichannel seismic data from deeper parts of the lake clearly image several clinoform structures stacked on top of each other. These stacked clinoforms indicate significantly lower lake levels prior to MIS 6 and a stepwise rise of water level with intermittent stillstands since its existence as water-filled body, which might have caused enhanced expansion of endemic species within Lake Ohrid.

Quantitative 3D strain analysis in analogue experiments: Integration of X-ray computed tomography and digital volume correlation techniques

The combination of scaled analogue experiments, material mechanics, X-ray computed tomography (XRCT) and Digital Volume Correlation techniques (DVC) is a powerful new tool not only to examine the 3 dimensional structure and kinematic evolution of complex deformation structures in scaled analogue experiments, but also to fully quantify their spatial strain distribution and complete strain history. Digital image correlation (DIC) is an important advance in quantitative physical modelling and helps to understand non-linear deformation processes. Optical non-intrusive (DIC) techniques enable the quantification of localised and distributed deformation in analogue experiments based either on images taken through transparent sidewalls (2D DIC) or on surface views (3D DIC). X-ray computed tomography (XRCT) analysis permits the non-destructive visualisation of the internal structure and kinematic evolution of scaled analogue experiments simulating tectonic evolution of complex geological structures. The combination of XRCT sectional image data of analogue experiments with 2D DIC only allows quantification of 2D displacement and strain components in section direction. This completely omits the potential of CT experiments for full 3D strain analysis of complex, non-cylindrical deformation structures. In this study, we apply digital volume correlation (DVC) techniques on XRCT scan data of “solid” analogue experiments to fully quantify the internal displacement and strain in 3 dimensions over time. Our first results indicate that the application of DVC techniques on XRCT volume data can successfully be used to quantify the 3D spatial and temporal strain patterns inside analogue experiments. We demonstrate the potential of combining DVC techniques and XRCT volume imaging for 3D strain analysis of a contractional experiment simulating the development of a non-cylindrical pop-up structure. Furthermore, we discuss various options for optimisation of granular materials, pattern generation, and data acquisition for increased resolution and accuracy of the strain results. Three-dimensional strain analysis of analogue models is of particular interest for geological and seismic interpretations of complex, non-cylindrical geological structures. The volume strain data enable the analysis of the large-scale and small-scale strain history of geological structures.

Chemical analysis of ODP Hole 104-642E

During Ocean Drilling Program Leg 104 a 900-m-thick sequence of volcanic rocks was drilled at Hole 642E on the Vøring Plateau, Norwegian Sea. This sequence erupted in two series (upper and lower series) upon continental basement. The upper series corresponds to the seaward-dipping seismic reflectors and comprises a succession of about 122 flows of transitional oceanic tholeiite composition. They have been subdivided into several formations consisting of flows related to each other by crystal fractionation processes, magma mixing, or both. Major- and trace-element chemistry indicates affinities to Tertiary plateau lavas of northeast Greenland and to Holocene lavas from shallow transitional segments of the Mid-Atlantic Ridge, such as Reykjanes Ridge. The tholeiitic magmas have been derived from a slightly LREE-depleted mantle source. Two tholeiitic dikes that intruded the lower series derive from an extremely depleted mantle source. Interlayered volcaniclastic sediments are dominantly ferrobasaltic and more differentiated. They appear to come from a LREE-enriched mantle source, and may have been erupted in close vicinity of the Vøring Plateau during hydroclastic eruptions. The two tholeiitic dikes that intruded the lower series as well as some flows at the base of the upper series show evidence of assimilation of continental upper crustal material.

Seismic data from Colle Gnifetti glacie saddle around borehole KCI, Swiss Alps

Two seismic surveys were carried out on the high-altitude glacier saddle, Colle Gnifetti, Monte Rosa, Italy/Switzerland. Explosive and vibroseismic sources were tested to explore the best way to generate seismic waves to deduce shallow and intermediate properties (<100 m) of firn and ice. The explosive source (SISSY) excites strong surface and diving waves, degrading data quality for processing; no englacial reflections besides the noisy bed reflector are visible. However, the strong diving waves are analyzed to derive the density distribution of the firn pack, yielding results similar to a nearby ice core. The vibrator source (ElViS), used in both P- and SH-wave modes, produces detectable laterally coherent reflections within the firn and ice column. We compare these with ice-core and radar data. The SH-wave data are particularly useful in providing detailed, high-resolution information on firn and ice stratigraphy. Our analyses demonstrate the potential of seismic methods to determine physical properties of firn and ice, particularly density and potentially also crystal-orientation fabric.

Sedimentation and mineral analysis on sediment cores from the Lena Delta

Core and outcrop analysis from Lena mouth deposits have been used to reconstruct the Late Quaternary sedimentation history of the Lena Delta. Sediment properties (heavy mineral composition, grain size characteristics, organic carbon content) and age determinations (14C AMS and IR-OSL) are applied to discriminate the main sedimentary units of the three major geomorphic terraces, which form the delta. The development of the terraces is controlled by complex interactions among the following four factors: (1) Channel migration. According to the distribution of 14C and IR-OSL age determinations of Lena mouth sediments, the major river runoff direction shifted from the west during marine isotope stages 5-3 (third terrace deposits) towards the northwest during marine isotope stage 2 and transition to stage 1 (second terrace), to the northeast and east during the Holocene (first terrace deposits). (2) Eustasy. Sea level rise from Last Glacial lowstand to the modern sea level position, reached at 6-5 ka BP, resulted in back-filling and flooding of the palaeovalleys. (3) Neotectonics. The extension of the Arctic Mid-Ocean Ridge into the Laptev Sea shelf acted as a halfgraben, showing dilatation movements with different subsidence rates. From the continent side, differential neotectonics with uplift and transpression in the Siberian coast ridges are active. Both likely have influenced river behavior by providing sites for preservation, with uplift, in particular, allowing accumulation of deposits in the second terrace in the western sector. The actual delta setting comprises only the eastern sector of the Lena Delta. (4) Peat formation. Polygenetic formation of ice-rich peaty sand (''Ice Complex'') was most extensive (7-11 m in thickness) in the southern part of the delta area between 43 and 14 ka BP (third terrace deposits). In recent times, alluvial peat (5-6 m in thickness) is accumulated on top of the deltaic sequences in the eastern sector (first terrace).

Structural behaviour and design criteria of under-deck cable-stayed bridges subjected to seismic action

Under-deck cable-stayed bridges are very effective structural systems for which the strong contribution of the stay cables under live loading allows for the design of very slender decks for persistent and transient loading scenarios. Their behaviour when subjected to seismic excitation is investigated herein and a set of design criteria are presented that relate to the type and arrangement of bearings, the number and configuration of struts, and the transverse distribution of stay cables. The nonlinear behaviour of these bridges when subject to both near-field and far-field accelerograms has been thoroughly investigated through the use of incremental dynamic analyses. An intensity measure that reflects the pertinent contributions to response when several vibration modes are activated was proposed and is shown to be effective for the analysis of this structural type. The under-deck cable-stay system contributes in a very positive manner to reducing the response when the bridges are subject to very strong seismic excitation. For such scenarios, the reduction in the stiffness of the deck because of crack formation, when prestressed concrete decks are used, mobilises the cable system and enhances the overall performance of the system. Sets of natural accelerograms that are compliant with the prescriptions of Eurocode 8 were also applied to propose a set of design criteria for this bridge type in areas prone to earthquakes. Particular attention is given to outlining the optimal strategies for the deployment of bearings

The use of direct boundary element method for gaining insight into complex seismic site response

The boundary element method is specially well suited for the analysis of the seismic response of valleys of complicated topography and stratigraphy. In this paper the method’s capabilities are illustrated using as an example an irregularity stratified (test site) sedimentary basin that has been modelled using 2D discretization and the Direct Boundary Element Method (DBEM). Site models displaying different levels of complexity are used in practice. The multi-layered model’s seismic response shows generally good agreement with observed data amplification levels, fundamental frequencies and the high spatial variability. Still important features such as the location of high frequencies peaks are missing. Even 2D simplified models reveal important characteristics of the wave field that 1D modelling does not show up.

Seismic Hazard and Ground Motion Characterization at the Itoiz Dam (Northern Spain).

This paper presents a new hazard-consistent ground motion characterization of the Itoiz dam site, located in Northern Spain. Firstly, we propose a methodology with different approximation levels to the expected ground motion at the dam site. Secondly, we apply this methodology taking into account the particular characteristics of the site and of the dam. Hazard calculations were performed following the Probabilistic Seismic Hazard Assessment method using a logic tree, which accounts for different seismic source zonings and different ground-motion attenuation relationships. The study was done in terms of peak ground acceleration and several spectral accelerations of periods coinciding with the fundamental vibration periods of the dam. In order to estimate these ground motions we consider two different dam conditions: when the dam is empty (T = 0.1 s) and when it is filled with water to its maximum capacity (T = 0.22 s). Additionally, seismic hazard analysis is done for two return periods: 975 years, related to the project earthquake, and 4,975 years, identified with an extreme event. Soil conditions were also taken into account at the site of the dam. Through the proposed methodology we deal with different forms of characterizing ground motion at the study site. In a first step, we obtain the uniform hazard response spectra for the two return periods. In a second step, a disaggregation analysis is done in order to obtain the controlling earthquakes that can affect the dam. Subsequently, we characterize the ground motion at the dam site in terms of specific response spectra for target motions defined by the expected values SA (T) of T = 0.1 and 0.22 s for the return periods of 975 and 4,975 years, respectively. Finally, synthetic acceleration time histories for earthquake events matching the controlling parameters are generated using the discrete wave-number method and subsequently analyzed. Because of the short relative distances between the controlling earthquakes and the dam site we considered finite sources in these computations. We conclude that directivity effects should be taken into account as an important variable in this kind of studies for ground motion characteristics.

Modelling of Topographic Irregularities for Seismic Site Response

Seismic evaluation methodology is applied to an existing viaduct in the south of Spain, near Granada, which is a medium seismicity region. The influence of both geology and topography in the spatial variability of ground motion are studied as well as seismic hazard analysis and ground motion characterization. Artificial hazard-consistent ground motion records are synthesised applying seismic hazard analysis and site effects are estimated through a diffraction study. Direct BEM is used to calculate the valley displacement response to vertically propagating SV waves and transfer functions are generated allowing the transformation of free field motion to motion at each support. A closed formulae is used to estimate these transfer function. Finally, the results obtained are compared.

Multi-hazard analysis and identification of priority settlements for land management in Haiti

Seismic hazard study in “La Hispaniola” island in connection with the land tenure situation in the region, in order to define priority areas with a high risk, where some land management recommendations are proposed. The seismic hazard assessment has been carried out following the probabilistic method with a seismogenic zonation and including the major faults of the region as independent units. In order to identify the priority areas, it has taken into account, besides the seismic hazard study, the map of changes of static Coulomb failure stress and the landslide hazard map.

An evaluation of seismic hazard in La Hispaniola, after the 2010 Haiti earthquake

An evaluation of the seismic hazard in La Hispaniola Island has been carried out, as part of the cooperative project SISMO-HAITI, supported by the Technical University of Madrid (UPM) and developed by several Spanish Universities, the National Observatory of Environment and Vulnerability) ONEV of Haiti, and with contributions from the Puerto Rico Seismic Network (PRSN) and University Seismological Institute of Dominican Republic (ISU). The study was aimed at obtaining results suitable for seismic design purposes. It started with the elaboration of a seismic catalogue for the Hispaniola Island, requiring an exhaustive revision of data reported by more than 20 seismic agencies, apart from these from the PRSN and ISU. The final catalogue contains 96 historical earthquakes and 1690 instrumental events, and it was homogenized to moment magnitude, Mw. Seismotectonic models proposed for the region were revised and a new regional zonation was proposed, taking into account geological andtectonic data, seismicity, focal mechanisms, and GPS observations. In parallel, attenuation models for subduction and crustal zones were revised in previous projects and the most suitable for the Caribbean plate were selected. Then, a seismic hazard analysis was developed in terms of peak ground acceleration, PGA, and spectral accelerations, SA (T), for periods of 0.1, 0.2, 0.5, 1 and 2s, using the Probabilistic Seismic Hazard Assessment (PSHA) methodology. As a result, different hazard maps were obtained for the quoted parameters, together with Uniform Hazard Spectra for Port au Prince and the main cities in the country. Hazard deaggregation was also carried out in these towns, for the target motion given by the PGA and SA (1s) obtained for return periods of 475, 975 and 2475 years. Therefore, the controlling earthquakes for short- and long-period target motions were derived. This study was started a few months after the 2010 earthquake, as a response to an aid request from the Haitian government to the UPM, and the results are available for the definition of the first building code in Haiti.