Recursive seismic ray modelling: applications in inversion and VSP

The recursive nature of rays in blocky models can be exploited to solve some difficult problems in seismic modelling. Each segment of a ray travels from an initial point up to a reflecting interface, where it is split into reflected and transmitted ray segments, which each continue in a similar way. The tree structure that thus emanates is conveniently handled by a recursive scheme. Recursion allows an automatic generation of all phases on a seismogram, together with all information necessary to analyse or select them. By operating recursively with a ray cell, bounded by a pair of vicinal rays in 2D, or a triplet of vicinal rays in 3D, and two successive isochrons, the two-point ray-tracing problem is reduced to a simple interpolation. Also, the cellular approach allows for a stable and robust evaluation of dynamic ray quantities without any paraxial tracing, which is cumbersome in blocky models of realistic complexity. Geometric shadows are filled by recursively generated diffractions. The recursive ray tracer has found applications in the fast computation of Green's functions in target-oriented inversion and in phase identification in VSP.     

Heterogeneous distributions of CO2 may be more important for dissolution and karstification in coastal eogenetic limestone than mixing dissolution

Mixing dissolution, a process whereby mixtures of two waters with different chemical compositions drive undersaturation with respect to carbonate minerals, is commonly considered to form cavernous macroporosity (e.g. flank margin caves and banana holes) in eogenetic karst aquifers. On small islands, macroporosity commonly originates when focused dissolution forms globular chambers lacking entrances to the surface, suggesting that dissolution processes are decoupled from surface hydrology. Mixing dissolution has been thought to be the primary dissolution process because meteoric water would equilibrate rapidly with calcium carbonate as it infiltrates through matrix porosity and because pCO₂ was assumed to be homogeneously distributed within the phreatic zone. Here, we report data from two abandoned well fields in an eogenetic karst aquifer on San Salvador Island, Bahamas, that demonstrate pCO₂ in the phreatic zone is distributed heterogeneously. The pCO₂ varied from less than log ?2.0 to more than log ?1.0 atm over distances of less than 30 m, generating dissolution in the subsurface where water flows from regions of low to high pCO₂ and cementation where water flows from regions of high to low pCO_2. Using simple geochemical models, we show dissolution caused by heterogeneously distributed pCO₂ can dissolve 2.5 to 10 times more calcite than the maximum amount possible by mixing of freshwater and seawater. Dissolution resulting from spatial variability in pCO₂ forms isolated, globular chambers lacking initial entrances to the surface, a morphology that is characteristic of flank margin caves and banana holes, both of which have entrances that form by erosion or collapse after cave formation. Our results indicate that heterogeneous pCO₂, rather than mixing dissolution, may be the dominant mechanism for observed spatial distribution of dissolution, cementation and macroporosity generation in eogenetic karst aquifers and for landscape development in these settings. Copyright © 2015 John Wiley & Sons, Ltd.     

Field Measurements and Numerical Simulations of the 2004 Tsunami Impact on the East Coast of Sri Lanka

This paper employs a numerical model of tsunami propagation together with documented observations and field measurements of the evidence left behind by the tsunami in December 2004, to identify and interpret the factors that have contributed to the significant spatial variability of the level of tsunami impact along the coastal belt of the eastern province of Sri Lanka. The model results considered in the present analysis include the distribution of the amplitude of the tsunami and the pattern of wave propagation over the continental shelf off the east coast, while the field data examined comprise the maximum water levels measured at or near the shoreline, the horizontal inundation distances and the number of housing and other buildings damaged. The computed maximum amplitude of the tsunami at water points nearest the shoreline along the east coast shows considerable variation ranging from 2.2 m to 11.4 m with a mean value of 5.7 m; moreover, the computed amplitudes agree well with the available field measurements. We also show that the shelf bathymetry off the east coast, particularly the submarine canyons at several locations, significantly influences the near-shore transformation of tsunami waves, and consequently, the spatial variation of the maximum water levels along the coastline. The measured values of inundation also show significant variation along the east coast and range from 70 m to 4560 m with a median value of 700 m. Our analyses of field data also show the dominant influence of the coastal topography and geomorphology on the extent of tsunami inundation. Furthermore, the measured inundation distances indicate no apparent correlation with the computed tsunami heights at the respective locations. We also show that both the computed tsunami heights and the measured inundation distances for the east coast closely follow the log-normal statistical distribution.     

Gratkorn: A benchmark locality for the continental Sarmatian s.str. of the Central Paratethys

This paper presents one of the richest and most complete vertebrate faunas of the late Middle Miocene (~12 Ma) of Central Europe. Up to now, sixty-two vertebrate taxa, comprising all major groups (fishes, amphibians, reptiles, birds, mammals), have been recorded. Based on sedimentological and palaeobiological evidences, this Fossillagerstätte is assumed to originate from a floodplain paleosol formed on top of a braided river sequence. The fauna points to a highly structured, somewhat vegetated landscape with a wide array of habitats (e.g., fluvial channels, sporadically moist floodplains, short-lived ponds, savannah-like open areas and screes). It was preserved due to a rapid drowning and the switch to a freshwater lake environment. Palaeoclimatological data, derived from pedogenic features as well as from biota, indicate an overall semi-arid, subtropical climate with distinct seasonality (mean annual precipitation 486 ± 252 mm, mean annual temperature ~15°C). This underlines the late Middle/early Late Miocene dry-spell in Central Europe. From taphonomical point of view, the irregularly distributed but roughly associated larger vertebrate remains refer to an in situ accumulation of the bone bed. Splintered bones, gnawing marks as well as rhizoconcretions and root corrosion structures record some pre- and post-burial modification of the taphocoenose. However, the findings of pellet remains argue for a very fast burial and thus to a low degree of time-averaging. For this reason, the fossil fauna reflects the original vertebrate community rather well and is a cornerstone for the understanding of late Middle Miocene terrestrial ecosystems in this region. Certainly, Gratkorn will be one of the key faunas for a high-resolution continental biostratigraphy and the comprehension of Europe’s faunal interchanges near the Middle/Late Miocene transition.     

Uncertainty analysis of urban sewer system using spatial simulation of radar rainfall fields: New York City case study

The goal of this study is to investigate the uncertainty of an urban sewer system’s response under various rainfall and infrastructure scenarios by applying a recently developed nonparametric copula-based simulation approach to extreme rainfall fields. The approach allows for Monte Carlo simulation of multiple variables with differing marginal distributions and arbitrary dependence structure. The independent and identically distributed daily extreme rainfall events of the corresponding urban area, extracted from nationwide high resolution radar data stage IV, are the inputs of the spatial simulator. The simulated extreme rainfall fields were used to calculate excess runoff using the Natural Resources Conservation Service’s approach. New York City is selected as a case study and the results highlight the importance of preserving the spatial dependence of rainfall fields between the grids, even for simplified hydrologic models. This study estimates the probability of combined sewer overflows under extreme rainfall events and identifies the most effective locations in New York City to install green infrastructure for detaining excess stormwater runoff. The results of this study are beneficial for planners working on stormwater management and the approach is broadly applicable because it does not rely on extensive sewer system information.     

A rifted inside corner massif on the Mid-Atlantic Ridge at 5°S

The structure of the Mid-Atlantic Ridge at 5°S was investigated during a recent cruise with the FS Meteor. A major dextral transform fault (hereafter the 5°S FZ) offsets the ridge left-laterally by 80 km. Just south of the transform and to the west of the median valley, the inside corner (IC – the region bounded by the ridge and the active transform) is marked by a major massif, characterized by a corrugated upper surface. Fossil IC massifs can also be identified further to the west. Unusually, a massif almost as high as the IC massif also characterizes the outside corner (OC) south of the inactive fracture zone and to the east of the median valley. This OC massif has axis-parallel dimensions identical to the IC massif and both are bounded on their sides closest to the spreading axis by abrupt, steep slopes. An axial volcanic ridge is well developed in the median valley both south of the IC/OC massifs and in an abandoned rift valley to the east of the OC massif, but is absent along the new ridge-axis segment between the IC and OC massifs. Wide-angle seismic data show that between the massifs, the crust of the median valley thins markedly towards the FZ. These observations are consistent with the formation of the OC massif by the rifting of an IC core complex and the development of a new spreading centre between the IC and OC massifs. The split IC massif presents an opportunity to study the internal structure of the footwall of a detachment fault, from the corrugated fault surface to deeper beneath the fault, without recourse to drilling. Preliminary dredging recovered gabbros from the scarp slope of the rifted IC massif, and serpentinites and gabbros from the intersection of this scarp with the corrugated surface. This is compatible with a concentration of serpentinites along the detachment surface, even where the massif internally is largely plutonic in nature.     

Stable isotopes and metal contamination in caged marine mussel Mytilus galloprovincialis

Metal concentrations and isotopic composition were measured in different tissues of the mussel Mytilus galloprovincialis in waters of the Balearic Islands (Western Mediterranean) in order to assess pollution levels. The isotopic composition was correlated with lead, cadmium, selenium and nickel obtained from the digestive gland and foot of the mussels. Significant negative correlations were found between cadmium, selenium and zinc and the mussel foot, mainly for ¹³C. Significant correlations were also found between lead and cadmium and the digestive gland. Pearson correlations indicated that the ¹³C isotopic signal in foot is a good proxy for the concentration of metals such as lead, cadmium, selenium and zinc. Similarly, ¹⁵N isotopic signatures in the digestive gland reflected the lead and cadmium concentration.     

Seismic gaps and plate tectonics: Seismic potential for major boundaries

The theory of plate tectonics provides a basic framework for evaluating the potential for future great earthquakes to occur along major plate boundaries. Along most of the transform and convergent plate boundaries considered in this paper, the majority of seismic slip occurs during large earthquakes, i.e., those of magnitude 7 or greater. The concepts that rupture zones, as delineated by aftershocks, tend to abut rather than overlap, and large events occur in regions with histories of both long- and short-term seismic quiescence are used in this paper to delineate major seismic gaps.In detail, however, the distribution of large shallow earthquakes along convergent plate margins is not always consistent with a simple model derived from plate tectonics. Certain plate boundaries, for example, appear in the long term to be nearly aseismic with respect to large earthquakes. The identification of specific tectonic regimes, as defined by dip of the inclined seismic zone, the presence or absence of aseismic ridges and seamounts on the downgoing lithospheric plate, the age contrast between the overthrust and underthrust plates, and the presence or absence of back-arc spreading, have led to a refinement in the application of plate tectonic theory to the evaluation of seismic potential.The term seismic gap is taken to refer to any region along an active plate boundary that has not experienced a large thrust or strike-slip earthquake for more than 30 years. A region of high seismic potential is a seismic gap that, for historic or tectonic reasons, is considered likely to produce a large shock during the next few decades. The seismic gap technique provides estimates of the location, size of future events and origin time to within a few tens of years at best.The accompanying map summarizes six categories of seismic potential for major plate boundaries in and around the margins of the Pacific Ocean and the Caribbean, South Sandwich and Sunda (Indonesia) regions for the next few decades. These categories range from what we consider high to low potential for being the site of large earthquakes during that period of time. Categories 1, 2 and 6 define a time-dependent potential based on the amount of time elapsed since the last large earthquake. The remaining categories, 3, 4, and 5, are used for areas that have ambiguous histories for large earthquakes; their seismic potential is inferred from various tectonic criteria. These six categories are meant to be interpreted as forecasts of the location and size of future large shocks and should not be considered to be predictions in which a precise estimate of the time of occurrence is specified.Several of the segments of major plate boundaries that are assigned the highest potential, i.e., category 1, are located along continental margins, adjacent to centers of population. Some of them are hundreds of kilometers long. High priority should be given to instrumenting and studying several of these major seismic gaps since many are now poorly instrumented. The categories of potential assigned here provide a rationale for assigning prorities for instrumentation, for future studies aimed at predicting large earthquakes and for making estimates of tsunami potential.Lamont-Doherty Geological Observatory Contribution No. 2906.     

The Uplift of the Longmenshan Thrust Belt and Subsidence of the West Sichuan Foreland Basin

Based on fission track dating of apatite, and measurement of vitrinite reflectance of rock samples from the Longmenshan (Longmen Mountain)area and the West Sichuan foreland basin and computer modelling it is concluded that (l)the Songpan-Garze fold belt has uplifted at least by 3-4 km with an uplift rate of no less than 0.3-0.4 mm/a since 10 Ma B.P.; (2) the Longmenshan thrust nappe belt has uplifted at least by 5-6 km with an uplift rate of more than 0.5- 0.6 mm /a since 10 Ma B.P.; (3) the Longmenshan detachment belt has uplifted by 1 - 2 km at a rate of 0.016-0.032 mm/a since 60 Ma B.P.; (4) the West Sichuan foreland basin has uplifted by 1.7-3 km at a rate of 0.028-0.05 mm/a since 60 Ma B.P.; (5) the uplift rate of the area on the west side of the Beichuan-Yingxiu-Xiaoguanzi fault for the last 10 Ma is 40 times as much as that on its east side; (6) the uplifting of the the Songpan - Garze fold belt and the subsidence of the West Sichuan foreland basin 60 Ma ago exhibit a mirro-image correlation, i.e     

Dextral‐Slip Thrust Faulting and Seismic Events of the Ms 8.0 Wenchuan Earthquake,Longmenshan Mountains,Eastern Margin of the Tibetan Plateau

Abstract: Dextral-slip thrust movement of the Songpan-Garzê terrain over the Sichuan block caused the Ms 8.0 Wenchuan earthquake of May 12, 2008 and offset the Central Longmenshan Fault (CLF) along a distance of ~250 km. Displacement along the CLF changes from Yingxiu to Qingchuan. The total oblique slip of up to 7.6 m in Yingxiu near the epicenter of the earthquake, decreases northeastward to 5.3 m, 6.6 m, 4.4 m, 2.5 m and 1.1 m in Hongkou, Beichuan, Pingtong, Nanba and Qingchuan, respectively. This offset apparently occurred during a sequence of four reported seismic events, EQ1–EQ4, which were identified by seismic inversion of the source mechanism. These events occurred in rapid succession as the fault break propagated northeastward during the earthquake. Variations in the plunge of slickensides along the CLF appear to match these events. The Mw 7.5 EQ1 event occurred during the first 0–10 s along the Yingxiu-Hongkou section of the CLF and is characterized by 1.7 m vertical slip and vertical slickensides. The Mw 8.0 EQ2 event, which occurred during the next 10–42 s along the Yingxiu-Yanziyan section of the CLF, is marked by major dextral-slip with minor thrust and slickensides plunging 25°–35° southwestward. The Mw 7.5 EQ3 event occurred during the following 42–60 s and resulted in dextral-slip and slickensides plunging 10° southwestward in Beichuan and plunging 73° southwestward in Hongkou. The Mw 7.7 EQ4 event, which occurred during the final 60–95 s along the Beichuan-Qingchuan section of the CLF, is characterized by nearly equal values of dextral and vertical slips with slickensides plunging 45°–50° southwestward. These seismic events match and evidently controlled the concentrations of landslide dams caused by the Wenchuan earthquake in Longmenshan Mountains.