Guided wave propagation in laterally varying media - II. Lg-waves in north-western Europe

Summary. Lg -wave observations at NORSAR from a sequence of explosions straddling the Central Graben in the North Sea, show that the phase is almost extinguished by passing through this structure. For models of the graben structure, based on seismic refraction experiments, with crustal thinning beneath a zone of thickened sediments, numerical modelling of Lg propagation shows very poor transmission. The inclusion of attenuation simulating the effects of scattering from faults and other complex structures further diminishes the transmittivity. The theoretical effect of a crustal pinch is to redistribute S energy over a wide range of groups velocities, a pattern that fits well with the observations.
Three components recordings at NORSAR of explosive charges fired in water covering a wide range of azimuths, show significant transverse components to the seismograms. The transverse energy builds up with travel time and for Lg is of the same order as the vertical component. These results are consistent with the theoretical prediction of progressive net transfer of energy from Rayleigh to Love modes as the wavetrain propagates through a three-dimensionally heterogeneous medium.     

Origin of dacites of Taupo Volcanic Zone, New Zealand

Dacites form a relatively small proportion of lavas in the Taupo Volcanic Zone, New Zealand (5km³), and occur mainly on the eastern side. In this paper their origin is considered in terms of three processes: (a) partial melting of crustal rocks; (b) fractional crystallisation of basalt and andesite; and (c) sub-surface mixing of basic and acid magma. Modelling techniques are used to calculate the most acceptable degree of fractional crystallisation and magma mixing to fit major-element data, and these values are used to compare calculated and observed trace-element values. The success or failure of the model is determined by the closeness of the two sets of values. For partial-melt models, trace-element values alone are calculated by the batch-melting equation.Results indicate that White Island dacite can best be modelled by fractional crystallisation; Manawahe by fractional crystallisation plus limited crustal contamination; Maungaongaonga by partial melting of Western Basement greywacke, and Tauhara by partial melting of this greywacke together with minor mixing with a more basic magma. Results from Parekauau and Horohoro indicate that these lavas are unlikely to have formed by any of the processes examined.     

Distinct metamorphic evolution of alternating silica-saturated and silica-deficient microdomains within garnet in ultrahigh-temperature granulites: An example from Sri Lanka

Here we report the occurrence of garnet porphyroblasts that have overgrown alternating silica-saturated and silica deficient microdomains via different mineral reactions. The samples were collected from ultrahigh-temperature (UHT) metapelites in the Highland Complex, Sri Lanka. In some of the metapelites, garnet crystals have cores formed via a dehydration reaction, which had taken place at silica-saturated microdomains and mantle to rim areas formed via a dehydration reaction at silica-deficient microdomains. In contrast, some other garnets in the same rock cores had formed via a dehydration reaction which occurred at silica-deficient microdomains while mantle to rim areas formed via a dehydration reaction at silica-saturated microdomains. Based on the textural observations, we conclude that the studied garnets have grown across different effective bulk compositional microdomains during the prograde evolution. These microdomains could represent heterogeneous compositional layers (paleobedding/laminations) in the precursor sediments or differentiated crenulation cleavages that existed during prograde metamorphism. UHT metamorphism associated with strong ductile deformation, metamorphic differentiation and crystallization of locally produced melt may have obliterated the evidence for such microdomains in the matrix. The lack of significant compositional zoning in garnet probably due to self-diffusion during UHT metamorphism had left mineral inclusions as the sole evidence for earlier microdomains with contrasting chemistry.     

Source investigation of a small event using empirical Green's functions and simulated annealing

We propose a two-step inversion of three-component seismograms that (1) recovers the far-field source time function at each station and (2) estimates the distribution of co-seismic slip on the fault plane for small earthquakes (magnitude 3 to 4). The empirical Green's function (EGF) method consists of finding a small earthquake located near the one we wish to study and then performing a deconvolution to remove the path, site, and instrumental effects from the main-event signal.
The deconvolution between the two earthquakes is an unstable procedure: we have therefore developed a simulated annealing technique to recover a stable and positive source time function (STF) in the time domain at each station with an estimation of uncertainties. Given a good azimuthal coverage, we can obtain information on the directivity effect as well as on the rupture process. We propose an inversion method by simulated annealing using the STF to recover the distribution of slip on the fault plane with a constant rupture-velocity model. This method permits estimation of physical quantities on the fault plane, as well as possible identification of the real fault plane.
We apply this two-step procedure for an event of magnitude 3 recorded in the Gulf of Corinth in August 1991. A nearby event of magnitude 2 provides us with empirical Green's functions for each station. We estimate an active fault area of 0.02 to 0.15 km² and deduce a stress-drop value of 1 to 30 bar and an average slip of 0.1 to 1.6 cm. The selected fault of the main event is in good agreement with the existence of a detachment surface inferred from the tectonics of this half-graben.     

Formalization of assessing radioactive waste burial site seismic regime parameters from seismological and geological data

A model is proposed that shows the relation of the block structure of the crust and earthquake sources (Sadovskii, 1979; Rodionov, 1979, 1984, 1994; Bugaev, 1999, 2011, 2014, 2015). The model can formalize how to assess the prediction of seismic regime parameters depending on the elastic limit and conditions and rate of deformation of the Earth’s crust. The spent nuclear fuel repository site in Olkiluoto (Finland) and a site in the area of the Krasnoyarsk Mining and Chemical Combine are considered as examples. It is demonstrated that the parameters of the prediction graphs limit the location of the points of magnitude repeatability graphs calculated for a site based on samples of earthquakes in the area according to different authors. This makes it possible to recommend predictive assessment of seismic regime parameters for stability monitoring of the seismic regime and safety analysis of a geological environment’s insulation properties for waste sites from the results of seismological monitoring and high-precision observations of modern movements of the Earth’s crust.     

Petrology of the Izmir-Karaburun volcanic area (West Turkey)

The Izmir-Karaburun region is located on the West coast of Turkey. In this area volcanic rocks of the late Miocene-Pliocene age outcrop. On the basis of the collected petrographic and geochemical data it has been possible to subdivide these rocks in to three series:a) calc-alkaline series of Karaburun-Koca dag-Izmir (quantitatively the most important). This series is formed by latite-andesites-dacites-rhyodacites.b) Silicic series of Izmir-Lebedos, mainly constituted by alkali rhyolitic rocks.c) Urla series, formed by alkali trachytes and alkali rhyolites, associated with scarce basic lavas of hawaiitic type. A different genesis is assumed for these series. In a first phase the latite-andesitic magma was formed by a partial melting in the lower crust or in the upper mantle. Afterwards a subcrustal magma with alkali basaltic affinity rose slowly through the crust forming an intermediate reservoir and differentiating predominantly towards alkali trachytic terms. Finally silicic magma of Izmir-Lebedos was formed by an anatectic process. It is possible that the fusion has been favoured by the presence of basic magma in the upper crust.     

Experiments on scalar dispersion within a model plant canopy,part III: An elevated line source

An experiment is reported in which heat was released as a passive tracer from an elevated lateral line source within a model plant canopy, with h _s = 0.85 h _c (h _s and h _c being the source and canopy heights, respectively). A sensor assembly consisting of three coplanar hot wires and one cold wire was used to measure profiles of mean temperature % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaiikamaana% aabaGaeqiUdehaaiaacMcaaaa!390C!\[(\overline \theta )\], temperature variance (Σ_θ ²), vertical and streamwise turbulent heat fluxes, and third moments of wind and temperature fluctuations. Conclusions were:

1. Despite the very heterogeneous flow within the canopy, the observed dispersive heat flux (due to spatial correlation between time-averaged temperature and vertical velocity) was small. However, there is evidence from the plume centroid (which was lower than h _s at the source) of systematic recirculating motions within the canopy.
2. The ratio % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr% 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9% vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x% fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeq4Wdm3aaS% baaSqaaiabeI7aXjaab2gacaqGHbGaaeiEaaqabaGccaGGVaWaa0aa% aeaacqaH4oqCaaWaaSbaaSqaaiaab2gacaqGHbGaaeiEaaqabaaaaa!41DF!\[\sigma _{\theta {\text{max}}} /\overline \theta _{{\text{max}}} \] (of maximum values on vertical profiles) decreased from 1 near the source to an asymptotic value of 0.4 far downstream, in good agreement with previous experimental and theoretical work for concentration fluctuations in the surface layer well above the canopy.
3. The eddy diffusivity for heat from the line source (K _HL ) increased, downstream of the source, to a nearly constant ‘far-field’ vertical profile. Within the canopy, the far-field K _HL was an order of magnitude larger than K _HP , the equivalent diffusivity for a plane source; well above the canopy, the two were equal. The time scale defined by (far-field K _HL )/(vertical velocity variance) was independent of height within the canopy.
4. Budgets for temperature variance, vertical heat flux and streamwise heat flux are remarkably similar to the equivalent budgets for an elevated line source in the surface layer well above the canopy, except in the lower part of the canopy in the far field, where vertical transport is much more important than in the surface layer.
5. A random flight simulation of the mean height and depth of the temperature plume was generally in good agreement with experiment. However, details of the temperature and streamwise turbulent heat flux profiles were not correct, suggesting that the model formulation needs to be improved.

     

LEVINSON INVERSION OF EARTHQUAKE GEOMETRY SH-TRANSMISSION SEISMOGRAMS IN THE PRESENCE OF NOISE*

The Kunetz-Claerbout equation for the acoustic transmission problem in a layered medium in its original form establishes the relation between the transmission and the reflec tion response for P-waves in an horizontally layered medium and with vertical incidence. It states that the reflection seismogram due to an impulsive source at the surface, is one side of the autocorrelation of the seismogram due to an impulsive source at depth and a surface receiver. By adapting Claerbout's formulation to the transmission of SH-waves, the Kunetz-Claerbout equation also holds for reflection and transmission coefficients dependent on the incident angle. Thus, earthquake geometry SH-transmission seismograms can be used to caculate corresponding pseudoreflection seismograms which, in turn, can be inverted for the impedance structure using the Levinson algorithm. If the average incidence angle is known, a geometrical correction on the resulting impedance model can improve the resolution of layer thicknesses. In contrast to the inversion of reflection seismograms, the Levinson algorithm is shown to yield stable results for the inversion of transmission seismograms even in the presence of additive noise. This noise stabilization is inherent to the Kunetz-Claerbout equation. Results of inverted SH-wave microearthquake seismograms from the Swabian Jura, SW Germany, seismic zone obtained at recording site Hausen im Tal have been compared with sonic-log data from nearby exploration drilling at Trochtelfingen. The agreement of the main structural elements is fair to a depth of several hundred metres.