ALBERT

Description: Siwi caldera, in the Vanuatu arc (Tanna island), is a rare volcanic complex where both persistent eruptive activity (Yasur volcano)and rapid block resurgence (Yenkahe horst) can be investigated simultaneously during a post-caldera stage. Here we provide new constraints on the feeding system of this volcanic complex, based on a detailed study of the petrology, geochemistry and volatile content of Yasur^Siwi bulk-rocks and melt inclusions, combined with measurements of the chemical composition and mass fluxes of Yasur volcanic gases. Major and trace element analyses of Yasur^ Siwi volcanic rocks, together with literature data for other volcanic centers, point to a single magmatic series and possibly long-lived feeding of Tanna volcanism by a homogeneous arc basalt. Olivine-hosted melt inclusions show that the parental basaltic magma, which produces basaltic-trachyandesites to trachyandesites by 50^70% crystal fractionation, is moderately enriched in volatiles ( 1wt % H2O, 0·1wt % S and 0·055 wt % Cl). The basaltic-trachyandesite magma, emplaced at between 4^5 km depth and the surface, preserves a high temperature (1107 158C) and constant H2O content ( 1wt %) until very shallow depths, where it degasses extensively and crystallizes. These conditions, maintained over the past 1400 years of Yasur activity, require early water loss during basalt differentiation, prevalent open-system degassing, and a relatively high heat flow ( 109W). Yasur volcano releases on average 13·4 103 tons d 1 of H2O and 680 tons d 1 of SO2, but moderate amounts of CO2 (840 tons d 1), HCl (165 tons d 1), and HF (23 tons d 1). Combined with melt inclusion data, these gas outputs constrain a bulk magma degassing rate of 5 107 m3 a 1, about a half of which is due to degassing of the basaltic-trachyandesite. We compute that 25 km3 of this magma have degassed without erupting and have accumulated beneath Siwi caldera over the past 1000 years, which is one order of magnitude larger than the accumulated volume uplift of the Yenkahe resurgent block. Hence, basalt supply and gradual storage of unerupted degassed basaltictrachyandesite could easily account for (or contribute to) the Yenkahe block resurgence.

Description: Published

Description: 1077-1105

Description: 2.3. TTC - Laboratori di chimica e fisica delle rocce

Description: 2.4. TTC - Laboratori di geochimica dei fluidi

Description: JCR Journal

Description: restricted

Description: Stromboli is known for its persistent degassing and rhythmic strombolian activity occasionally punctuated by paroxysmal eruptions. The basaltic pumice and scoria emitted during paroxysms and strombolian activity, respectively, differ in their textures, crystal contents and glass matrix compositions, which testify to distinct conditions of crystallization, degassing and magma ascent. We present here an extensive dataset on major elements and volatiles (CO2, H2O, S and Cl) in olivine-hosted melt inclusions and embayments from pyroclasts emplaced during explosive eruptions of variable magnitude. Magma saturation pressures were assessed from the dissolved amounts of H2O and CO2 taking into account the melt composition evolution. Both pressures and melt inclusion compositions indicate that (1) Ca-basaltic melts entrapped in high-Mg olivines (Fo89–90) generate Stromboli basalts through crystal fractionation, and (2) the Stromboli plumbing system can be imaged as a succession of magma ponding zones connected by dikes. The 7–10 km interval, where magmas are stored and differentiate, is periodically recharged by new magma batches, possibly ranging from Ca-basalts to basalts, with a CO2-rich gas phase. These deep recharges promote the formation of bubbly basalt blobs, which are able to intrude the shallow plumbing system (2–4 km), where CO2 gas fluxing enhances H2O loss, crystallization and generation of crystal-rich, dense, degassed magma. Chlorine partitioning into the H2O–CO2-bearing gas phase accounts for its efficient degassing (≥69%) under the open-system conditions of strombolian activity. Paroxysms, however, are generated through predominantly closed-system ascent of basaltic magma batches from the deep storage zone. In this situation crystallization is negligible and sulfur exsolution starts at ≤170 MPa. Chlorine remains dissolved in the melt until lower pressures, only 16% being lost upon eruption. Finally, we propose a continuum in explosive eruption energy, from strombolian activity to large paroxysmal events, ultimately controlled by variable pressurization of the deep feeding system associated with magma and gas recharges.

Description: Published

Description: 603-626

Description: 2.3. TTC - Laboratori di chimica e fisica delle rocce

Description: 3.5. Geologia e storia dei vulcani ed evoluzione dei magmi

Description: JCR Journal

Description: reserved

Description: SUMMARY In his comment, De Basabe criticises our paper ignoring the advantages of unstructured element mesh used in the finite element method, and argues that the Gaussian quadrature grid (GQG) approach is limited to a homogenous or constant layered geological model and does not have the spectral accuracy. In this reply, we give our response to his criticism and comment, and further clarify the accuracy and capability of the GQG approach.

Description: SUMMARY High-porosity channels are important pathways for melt migration in the mantle. To better understand the lithology and geometry of high-porosity melt channels, we conduct high-order accurate numerical simulations of reactive dissolution along a solubility gradient in an upwelling and viscously deformable porous column. In contrast to earlier studies, we assume the dissolution reaction to be at equilibrium, consider a finite soluble mineral abundance, and employ a high-order accurate numerical scheme. Using sustained perturbations in porosity and soluble mineral abundance at the inflow boundary, we explore the structure of steady-state high-porosity melt channels and their associated lithologies over a range of parameters, including soluble mineral abundance, solubility gradient, amplitude and lateral variation in inflow melt flux, melt fraction and upwelling rate. In general, high-porosity dunite channels are transient and shallow parts of pathways for melt migration in the mantle. The lower parts of a high-porosity channel are orthopyroxene-bearing dunite, harzburgite and possibly lherzolite. A wide orthopyroxene-free dunite channel may contain two or three high-porosity melt channels. The depth of dunite channel initiation depends on the solubility gradient, soluble mineral abundance, inflow melt flux, melt suction rate and upwelling rate. The amplitude and length scale of lateral variation in porosity and orthopyroxene abundance at the base of the upwelling column are important in determining the size and dimension of dunite channels, the strength of melt focusing and the melt suction rate, the presence of compacting boundary layer, as well as the number of high-porosity melt branches within a dunite channel. The spatial relations among the high-porosity melt channels, dunite and harzburgite channels documented in our numerical simulations may shed new light on a number of field, petrological and geochemical observations related to melt migration in the mantle.

Description: SUMMARY Geothermal energy is becoming an important clean energy source, however, the stimulation of a reservoir for an Enhanced Geothermal System (EGS) is associated with seismic risk due to induced seismicity. Seismicity occurring due to the water injection at depth have to be well recorded and monitored. To mitigate the seismic risk of a damaging event, an appropriate alarm system needs to be in place for each individual experiment. In recent experiments, the so-called traffic-light alarm system, based on public response, local magnitude and peak ground velocity, was used. We aim to improve the pre-defined alarm system by introducing a probability-based approach; we retrospectively model the ongoing seismicity in real time with multiple statistical forecast models and then translate the forecast to seismic hazard in terms of probabilities of exceeding a ground motion intensity level. One class of models accounts for the water injection rate, the main parameter that can be controlled by the operators during an experiment. By translating the models into time-varying probabilities of exceeding various intensity levels, we provide tools which are well understood by the decision makers and can be used to determine thresholds non-exceedance during a reservoir stimulation; this, however, remains an entrepreneurial or political decision of the responsible project coordinators. We introduce forecast models based on the data set of an EGS experiment in the city of Basel. Between 2006 December 2 and 8, approximately 11 500 m 3 of water was injected into a 5-km-deep well at high pressures. A six-sensor borehole array, was installed by the company Geothermal Explorers Limited (GEL) at depths between 300 and 2700 m around the well to monitor the induced seismicity. The network recorded approximately 11 200 events during the injection phase, more than 3500 of which were located. With the traffic-light system, actions where implemented after an M L 2.7 event, the water injection was reduced and then stopped after another M L 2.5 event. A few hours later, an earthquake with M L 3.4, felt within the city, occurred, which led to bleed-off of the well. A risk study was later issued with the outcome that the experiment could not be resumed. We analyse the statistical features of the sequence and show that the sequence is well modelled with the Omori–Utsu law following the termination of water injection. Based on this model, the sequence will last 31+29/−14 years to reach the background level. We introduce statistical models based on Reasenberg and Jones and Epidemic Type Aftershock Sequence (ETAS) models, commonly used to model aftershock sequences. We compare and test different model setups to simulate the sequences, varying the number of fixed and free parameters. For one class of the ETAS models, we account for the flow rate at the injection borehole. We test the models against the observed data with standard likelihood tests and find the ETAS model accounting for the on flow rate to perform best. Such a model may in future serve as a valuable tool for designing probabilistic alarm systems for EGS experiments.

Description: SUMMARY In principle, we can never know the true Green's function, which is a major error source in seismic waveform inversion. So far, many studies have devoted their efforts to obtain a Green's function as precise as possible. In this study, we propose a new strategy to cope with this problem. That is to say, we introduce uncertainty of Green's function into waveform inversion analyses. Due to the propagation law of errors, the uncertainty of Green's function results in a data covariance matrix with significant off-diagonal components, which naturally reduce the weight of observed data in later phases. Because the data covariance matrix depends on the model parameters that express slip distribution, the inverse problem to be solved becomes non-linear. Applying the developed inverse method to the teleseismic P -wave data of the 2006 Java, Indonesia, tsunami earthquake, we obtained a reasonable slip-rate distribution and moment-rate function without the non-negative slip constraint. The solution was independent of the initial values of the model parameters. If we neglect the modelling errors due to Green's function as in the conventional formulation, the total slip distribution is much rougher with significant opposite slip components, whereas the moment-rate function is much smoother. If we use a stronger smoothing constraint, more plausible slip distribution can be obtained, but then the moment-rate function becomes even smoother. By comparing the observed waveforms with the synthetic waveforms, we found that high-frequency components were well reproduced only by the new formulation. The modelling errors are essentially important in waveform inversion analyses, although they have been commonly neglected.

Description: SUMMARY Geological and geodetic records show that motion is distributed over a broad and complex system of faults in southern California, with many large earthquakes occurring off the main San Andreas fault. This unusually complex plate boundary geometry is imposed by the big bend of the San Andreas Fault and the dynamics of the underlying lithosphere, but lithosphere rheology and distribution of applied forces are still poorly constrained. Geodynamic modelling is thus very uncertain, and distribution of lithospheric deformations highly controversial. Here, exploiting the property that deformation orients olivine crystals, we show that deep lithospheric deformations can be mapped directly with a new imaging approach of seismic anisotropy, with a resolution fine enough for detailed geodynamic interpretation. This method relies entirely on finite-frequency effects in the splitting of SKS waves. We build an extensive data set of more than 3400 SKS splitting measurements from the analysis of SKS waves recorded by all the broad-band stations in southern California. Resolution tests demonstrate that, using this data set, we are able to resolve anisotropic structures smaller than the size of the first Fresnel zone of SKS waves. The good vertical resolution in our tomographic images results from the short inter-station spacing in southern California. In our 3-D anisotropy model, we find no evidence for localized lithospheric shear deformation beneath the San Andreas Fault. Instead, the lithospheric plate boundary is localized along a broad shear zone beneath the East California Shear Zone. Therefore, surface and deep deformation patterns are poorly correlated and most likely decoupled. Active N–S convergence in the Transverse Ranges results in strongly coherent E–W alignment of olivine fast axes in the shallow lithosphere. At the top of the asthenosphere, the observed low level of anisotropy suggests weak lithospheric basal tractions and thus a strong decoupling between the lithosphere and the asthenosphere.

Description: SUMMARY Earthquakes often occur on faults separating materials with different elastic properties. On theoretical grounds, it is expected that earthquakes on such bimaterial interfaces might have a preferred rupture propagation direction, that being the direction of motion of the more compliant material. The goal of this paper is to determine whether a large sample of natural earthquakes on a bimaterial interface exhibits this tendency. Since the creeping section of San Andreas Fault has a large across-fault velocity contrast and has produced thousands of microearthquakes over the last few decades, the rupture directions of ∼3100 magnitude 0.5–3.0 earthquakes were studied using spectral ratio analysis. The spectral ratios of all earthquake pairs in spatially defined clusters were fitted with synthetic spectral ratios at qualified stations. The synthetics were computed from a simple moving point source model in which each modelled earthquake has four parameters: two rupture lengths (one to the SE and one to the NW) and their propagation velocities. The resolution of rupture directivity increases with event size, such that nearly 900 events, mostly those larger than ∼70 m, appear reasonably well resolved. The inversion results suggest that ∼40 per cent of the well-resolved events are roughly bilateral, although more than ∼80 per cent of the 144 events classified as "strongly unilateral" rupture to the SE, consistent with the theoretical prediction. For those rupture halves that were large enough for the propagation speed to be somewhat resolved, that speed was greater by roughly 10 per cent for those halves propagating to the SE, qualitatively consistent with numerical and laboratory experiments. We also found that events with nearby foreshocks within several hours tend to rupture away from those foreshocks, whether to the NW or to the SE, indicating that asymmetry of prior stressing history can exert a stronger influence on rupture directivity than the material contrast.

Description: SUMMARY We present results on the composition, structure and particle size distribution (PSD) of pulverized and damaged granitic rocks in a 42-m-deep core adjacent to the San Andreas Fault near Littlerock, CA. The cored section is composed of pulverized granites and granodiorites, and is cut by numerous mesoscopic secondary shears. The analysis employs XRD, XRF, thin sections and semi-automated particle size analyser methods, including a novel calibration method. The mean particle size for the majority of samples falls between 50 and 470 μm. The PSDs can be fitted by a power law, with D -values ranging between 2.5 and 3.1, as well as by a superposition of two Gaussians. Fracture surface energy calculations based on the observed particle distributions provide very low values, indicating that the part of the total earthquake energy budget expended for breaking or shattering rocks is small. Shear deformation is likely to dominate near secondary faults. The most pronounced fault-related alteration occurs along the secondary shears, and is a function of both composition and depth. The alteration to clay appears to be the result of fluid–rock interaction and brittle deformation under low temperature conditions, rather than of surface-related weathering. The particle size reduction noted in the core reflects multiple mechanisms of comminution. The zones of pulverization that lack significant weathering likely result from repeating episodes of dynamic dilation and contraction.

Description: SUMMARY Broad-band seismic data from the southern African seismic experiment and the AfricaArray network are used to investigate the seismic velocity structure of the upper mantle beneath southern Africa, and in particular beneath the Kaapvaal Craton. A two-plane approximation method that includes a finite frequency sensitivity kernel is employed to measure Rayleigh wave phase velocities, which are inverted to obtain a quasi-3-D shear wave velocity model of the upper mantle. We find phase velocities for the Kaapvaal Craton and surrounding mobile belts that are comparable to those reported by previous studies, and we find little evidence for variation from east to west across the Namaqua–Natal Belt, a region not well imaged in previous studies. A high-velocity upper-mantle lid is found beneath the Kaapvaal Craton and most of southern Africa. For the Kaapvaal Craton, the thickness of the lid (∼150–200 km) is consistent with the lid thicknesses reported in many previous studies. The cratonic lid is underlain by a ∼100-km thick low-velocity zone with a 3.9 per cent maximum velocity reduction. By comparing the velocity model to those published for other Archean cratons, we find few differences, and therefore conclude that there is little evidence in the shear wave velocity structure of the mantle to indicate that the southern African plateau is supported by an upper-mantle thermal anomaly.

Description: SUMMARY To better constrain Baltica's position within Pangea, we conducted a palaeomagnetic study of Permo-Triassic dykes from the Oslo Graben, as a follow-up to an initial, but rather limited, study by Torsvik and colleagues in 1998. The age of these so-called Lunner dykes had previously been determined as ∼240 Ma in that study, but details in their analyses and new 40 Ar/ 39 Ar ages reveal that there may have been some argon loss in the initially dated dyke minerals and that a combined (weighted mean) age of 271 ± 2.7 (2σ) Myr for the dykes is preferable. We find two major components of magnetization in our samples: one carried by an Fe-sulphide (likely pyrrhotite) and the other carried by low-Ti magnetite; these magnetization components may be found together (superposed) in a given sample or they may occur apart. Micronmetre-sized crystals of Ti-Fe oxides, observed with a scanning electron microscope (SEM) show exsolution lamellae, formed upon cooling from intrusion temperatures. Assuming that the submicronmetre-sized (Ti)-magnetite grains that carry a stable remanence are of the same generation as the observed larger grains, we interpret the magnetite remanence in the dykes as of primary, thermoremanent origin. The sulphide remanence appears to be slightly younger, as seen by the SEM observations of pyrite framboids and a Fe-sulphide grain invading a Ti-magnetite grain. Moreover, the sulphide mineralization is likely of region-wide hydrothermal origin. The magnetizations carried by the pyrrhotite and magnetite have nearly identical directions and so, must be nearly of the same age. For this study, we sampled 56 sites including 39 dykes, 10 baked-contact rocks and 7 host rocks removed from the immediate dyke contacts. The dykes and the contact rocks have the same SW and up directions of magnetization, and contain the Fe-sulphide or the magnetite magnetization or both, as diagnosed by their relative unblocking temperatures. However, all the sampled carbonate and igneous host rocks far away from the dykes also have the same directions. Thus, all of the 10 originally planned contact tests are inconclusive. The new palaeopoles of this study are a few degrees apart; the magnetite pole (from dykes only, N = 25) is located at 51°N, 164°E, K = 69, A 95 = 3.5°, whereas the pole calculated from iron sulphide magnetic directions (all rock types, N = 20) is at 54°N, 166°E, K = 112, A 95 = 3.1°. All directions are of reversed polarity, suggesting that the magnetization was acquired during the Kiaman Reversed Superchron. The palaeomagnetic mean result from the magnetite-bearing sites implies a palaeolatitude of Oslo of 23°N, whereas the palaeolatitude calculated from the pyrrhotite magnetizations is 25–27°N, depending on choice of host lithologies. As noted in many previous publications, the palaeomagnetic poles for the late Palaeozoic and Early-Middle Triassic are in conflict with classical Pangea reconstructions. The poles with ages of 250 ± 10 Ma, in particular, previously showed a discrepancy of some 25° or more, when the Gondwana and Laurussia continents are restored to their juxtapositions in the Pangea-A fit, before the opening of the Atlantic Ocean. Proposed solutions to this conundrum have been controversial, involving doubts about (1) the geocentric coaxial dipole field model, (2) the reliability of the palaeomagnetic results or their ages, or (3) the validity of the Pangea-A reconstruction, leading to proposals of a Pangea B reconstruction in which Gondwana is displaced some 3500 km to the east with respect to Laurussia. The significance of our new result for this Pangea controversy resides in its improved age within an early Guadelupian (mid to late Permian) time interval where few results exist from well-dated igneous rocks in either Baltica or Laurentia. There are quite a few results from sedimentary rocks, but these may be suspected to suffer inclination shallowing, and are therefore less suitable to settle a palaeolatitudinal argument. Our new result of the magnetite magnetization, granted it is primary and acquired at about 270 Ma, combined with a new ∼265 Ma result from Argentina and selected other poles from igneous rocks, leaves enough room for the north–south configuration of Pangea A at 270 Ma and avoids the overlap between Baltica and Gondwana that necessitated Pangea B, at least for the Late Permian.

Description: SUMMARY We investigate 3-D basin structures and site resonance frequencies in the İzmit Bay area of Turkey by new geophysical surveys based on 239 single station microtremor and 405-point gravity measurements. A fundamental resonance frequency map of the İzmit Bay was produced from the main peak in the horizontal-to-vertical component spectral ratio (HVSR) of microtremors. The HVSR analysis of the microtremor data reveals single, double, broad peaked or no peak type HVSR curves varying in accordance with the surface geology and spatial extent of the three basins present in the İzmit Bay area. In the deepest sections of the İzmit, Gölcük and Derince basins of the İzmit Bay, the fundamental resonance frequencies are dominantly 0.24–0.30 Hz. These resonance frequencies should be taken into consideration along with higher mode frequencies to construct earthquake resistant structures in the İzmit Bay area. The 3-D gravimetric bedrock depth map of the İzmit basin shows that the basin has an asymmetric shape with its deepest section coinciding with the surface trace of the North Anatolian Fault. The northern shoulder of the basin has a gentle dip on the Kocaeli Peneplain side and the southern shoulder is much steeper, and it is bounded by the Samanlı Mountains. We derive a power-law relationship that correlates the fundamental site resonance frequencies with the sedimentary cover thickness obtained from the gravity and shear wave velocity data in the İzmit Bay. We use this relationship to estimate bedrock depths beneath Gölcük and Derince basins. Our estimation of maximum basin depths is 1400 m for the İzmit and is 800 m for the Gölcük and Derince basins. Finally, we have analysed a converted Sp phase from a local earthquake recording made at site CMP to calculate and verify the sediment thickness estimations obtained from our gravimetric and microtremor analyses. This calculation shows close agreement with that of the gravimetric and microtremor results. Our results show that the basins in the İzmit Bay area have a very thick sedimentary cover with very low shear velocities underlined by hard bedrock, forming a sharp impedance contrast. We anticipate that these results will be a key contribution to the quantitative assessment of seismic hazard for the İzmit Bay area before the occurrence of strong earthquakes in the Marmara region.

Description: SUMMARY To constrain the forecast horizon of geomagnetic data assimilation, it is of interest to quantify the range of predictability of the geodynamo. Following earlier work in the field of dynamic meteorology, we investigate the sensitivity of numerical dynamos to various perturbations applied to the magnetic, velocity and temperature fields. These perturbations result in some errors, which affect all fields in the same relative way, and grow at the same exponential rate  λ=τ −1 e , independent of the type and the amplitude of perturbation. Errors produced by the limited resolution of numerical dynamos are also shown to produce a similar amplification, with the same exponential rate. Exploring various possible scaling laws, we demonstrate that the growth rate is mainly proportional to an advection timescale. To better understand the mechanism responsible for the error amplification, we next compare these growth rates with two other dynamo outputs which display a similar dependence on advection: the inverse  τ −1 SV  of the secular-variation timescale, characterizing the secular variation of the observable field produced by these dynamos; and the inverse  (τ mag diss ) −1  of the magnetic dissipation time, characterizing the rate at which magnetic energy is produced to compensate for Ohmic dissipation in these dynamos. The possible role of viscous dissipation is also discussed via the inverse  (τ kin diss ) −1  of the analogous viscous dissipation time, characterizing the rate at which kinetic energy is produced to compensate for viscous dissipation. We conclude that  τ e  tends to equate  τ mag diss  for dynamos operating in a turbulent regime with low enough Ekman number, and such that  τ mag diss 〈 τ kin diss . As these conditions are met in the Earth's outer core, we suggest that  τ e  is controlled by magnetic dissipation, leading to a value  τ e =τ mag diss ≈ 30  yr . We finally discuss the consequences of our results for the practical limit of predictability of the geodynamo.

Description: SUMMARY A closed-form solution is given for a 2-D, transverse electric mode, magnetotelluric (MT) problem. The model system consists of a finite vertical thin conductor with variable integrated conductivity over a perfectly conducting base. A notable property of the solution is that the frequency response possesses a single pole in the complex plane. Systems with finitely many resonances play a central role in the 1-D MT inverse problem based on finite data sets, but until now, no 2-D system of this kind was known. The particular model is shown to be just one of a large class of thin conductors with same the property, and further examples are given. The solutions of the induction problem for members of this family can often be written in compact closed form, making them the simplest known solutions to the 2-D MT problem.

Description: SUMMARY The eikonal equation in an attenuating medium has the form of a complex-valued Hamilton–Jacobi equation and must be solved in terms of the complex-valued traveltime (complex-valued action function). A very suitable approximate method for calculating the complex-valued traveltime right in real space is represented by the perturbation from the reference traveltime calculated along real-valued reference rays to the complex-valued traveltime defined by the complex-valued Hamilton–Jacobi equation. The real-valued reference rays are calculated using the reference Hamiltonian function. The perturbation Hamiltonian function is parametrized by one or more perturbation parameters, and smoothly connects the reference Hamiltonian function with the Hamiltonian function corresponding to a given complex-valued Hamilton–Jacobi equation. Both the reference Hamiltonian function and the perturbation Hamiltonian function may be constructed in different ways, yielding differently accurate perturbation expansions of traveltime. All present perturbation methods use reference rays calculated in a reference anisotropic non-attenuating medium. In this paper, the reference Hamiltonian function is constructed directly using the Hamiltonian function corresponding to a given complex-valued Hamilton–Jacobi equation, and the perturbation Hamiltonian function is linear with respect to the perturbation parameter. The direct construction of the reference Hamiltonian function from the given complex-valued Hamilton–Jacobi equation is very general and accurate, especially for homogeneous Hamiltonian functions of degree N =−1 with respect to the slowness vector.

Description: SUMMARY The lithospheric structure in Oman has been determined by analyzing teleseismic P -receiver functions recorded at broad-band and short-period seismic stations of the Oman Seismological Network. Lithospheric structure is obtained by jointly inverting receiver functions and Rayleigh wave group velocities derived from continental-scale surface wave tomography of Pasyanos. We observe relatively thick crust (40–48 km) within the ophiolite formed mountains in northern Oman. The crustal thickness is about 35 km within the passive continental margin of the southern Oman region. Uppermost (〈5 km) crustal shear wave velocities are faster in the northern ophiolite regions compared to the southern Oman region, while shear velocities in the middle crust are faster in the Southern Oman region compared to the ophiolite region. This observation coincides well with cretaceous to Eocene marine platform sequences overlying Precambrian to Cambrian basement of the southern part. Joint inversion analysis shows that the Moho depth of Oman varies from 34 km in the southern region to 48 km in the northern part.

Description: We present experimental data on the partitioning of Li, Be, B, K, Mg, Sr, Ga, Rb, Cs, Ba, La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, Lu, Th, U, Hf, Zr, Nb and Ta between lawsonite and fluid, and zoisite and fluid at 3·0–3·5 GPa and 650–850°C. The aim is to provide data bearing on the trace element contents of fluids released during dehydration of subducting oceanic crust. Experimental trace element partition coefficients for lawsonite indicate a preference for the light rare earth elements (LREE) over the heavy REE (HREE) and for Be. These characteristics are consistent with the chemical composition of lawsonite in natural rocks. Experimental trace element partition coefficients for zoisite indicate a preference for HREE relative to LREE. This observation, consistent with earlier experimental data, is the reverse of the observed trace element compositions of natural zoisites, indicating the influence of other factors on the trace element contents of this phase. Lattice strain theory explains well the experimentally derived partitioning of divalent cations in the Ca-site between lawsonite and fluid. However, the weak relative fractionation of REE between lawsonite and fluid cannot be explained by lattice strain theory, as previously observed for zoisite–fluid REE partitioning. We combine our experimental data with thermodynamic models of mineral stability to model the compositions of fluids released during subduction of altered normal mid-ocean ridge basalt. The low La/Sm ratio associated with very high Ba/Th in arc magmas can be explained only if allanite is stable in the subducting oceanic crust. This suggests that the crustal fluid component involved in arc magma petrogenesis results from processes occurring in the warm, top part of the subducting slab. Decreasing lawsonite modal proportion with depth is associated with a large release of fluid characterized by low B/Be ratios that could explain the decreasing B/Be ratios in arc magmas with increasing distance from the trench. This implies that an important Be input in arc magma originates from the fluid generated during oceanic crust dehydration.

Description: To contribute to our understanding of the mechanisms and pathways of fluid movement through deeply subducted crust, we investigate high-pressure veins cutting eclogite-facies (~2·0 GPa and ~600°C) metagabbros of the Monviso Ophiolite, Italian Western Alps. The veins consist mainly of omphacite with minor garnet, rutile, talc and accessory zircon. Most of the vein minerals have major and trace element compositions that are comparable with the host-rock minerals, and vein and host-rock zircons have similar Hf isotopic compositions. These observations support the conclusions of previous studies that these veins largely formed from a locally sourced hydrous fluid during prograde or peak metamorphism. However, the bulk-rock Cr and Ni contents of the veins are significantly higher than those of the surrounding host eclogites. We also document distinct Cr-rich (up to weight per cent levels) zones in omphacite, garnet and rutile in some vein samples. Vein garnet and talc also have relatively high MgO and Ni contents. X-ray maps of vein garnet and rutile grains reveal complex internal zoning features, which are largely defined by micrometre-scale variations in Cr content. Some grains have concentric and oscillatory zoning in Cr, whereas others feature a chaotic fracture-like pattern. These Cr-rich zones are associated with high concentrations of Ni, B, As, Sb, Nb, Zr and high ratios of light rare earth elements (LREE) to middle REE (MREE) compared with low-Cr vein and host-rock minerals. Petrological and mass-balance constraints verify that the Cr-rich zones in the veins were not derived from internally sourced fluids, but represent precipitates from an external fluid. The external source that is consistent with the distinctive trace element characteristics of the vein components is antigorite serpentinite, which forms the structural basement of the high-pressure metagabbros. We propose at least two separate growth mechanisms for the Monviso veins. Most vein infillings were formed during progressive prograde metamorphism from locally derived fluid. Influx of the serpentinite-derived or other external fluid was transient and episodic and was probably achieved via brittle fractures, which preferentially formed along the pre-existing vein structures. The dehydration of serpentinite at high pressures in subduction zones may provide crucial volatiles and trace elements for arc magmas. Our results indicate that the movement of these fluids through subducted oceanic crust is likely to be highly channeled and transient so the progressive development of vein systems in mafic rocks may also be crucial for forming channelways for long-distance fluid flow at depth in subduction zones.

Description: Siwi caldera, in the Vanuatu arc (Tanna island), is a rare volcanic complex where both persistent eruptive activity (Yasur volcano) and rapid block resurgence (Yenkahe horst) can be investigated simultaneously during a post-caldera stage. Here we provide new constraints on the feeding system of this volcanic complex, based on a detailed study of the petrology, geochemistry and volatile content of Yasur–Siwi bulk-rocks and melt inclusions, combined with measurements of the chemical composition and mass fluxes of Yasur volcanic gases. Major and trace element analyses of Yasur–Siwi volcanic rocks, together with literature data for other volcanic centers, point to a single magmatic series and possibly long-lived feeding of Tanna volcanism by a homogeneous arc basalt. Olivine-hosted melt inclusions show that the parental basaltic magma, which produces basaltic-trachyandesites to trachyandesites by ~50–70% crystal fractionation, is moderately enriched in volatiles (~1 wt % H 2 O, 0·1 wt % S and 0·055 wt % Cl). The basaltic-trachyandesite magma, emplaced at between 4–5 km depth and the surface, preserves a high temperature (1107 ± 15°C) and constant H 2 O content (~1 wt %) until very shallow depths, where it degasses extensively and crystallizes. These conditions, maintained over the past 1400 years of Yasur activity, require early water loss during basalt differentiation, prevalent open-system degassing, and a relatively high heat flow (~10 9 W). Yasur volcano releases on average ≥ 13·4 x 10 3 tons d –1 of H 2 O and 680 tons d –1 of SO 2 , but moderate amounts of CO 2 (840 tons d –1 ), HCl (165 tons d –1 ), and HF (23 tons d –1 ). Combined with melt inclusion data, these gas outputs constrain a bulk magma degassing rate of ~5 x 10 7 m 3 a –1 , about a half of which is due to degassing of the basaltic-trachyandesite. We compute that 25 km 3 of this magma have degassed without erupting and have accumulated beneath Siwi caldera over the past 1000 years, which is one order of magnitude larger than the accumulated volume uplift of the Yenkahe resurgent block. Hence, basalt supply and gradual storage of unerupted degassed basaltic-trachyandesite could easily account for (or contribute to) the Yenkahe block resurgence.

Description: Primitive basalts are rarely found in arcs. The active NW Rota-1 volcano in the Mariana arc has erupted near-primitive lavas, which we have sampled with ROV Hyper-Dolphin (HPD). Samples from the summit (HPD480) and eastern flank (HPD488) include 17 magnesian basalts (51–52 wt % SiO 2 ) with 7·5–9·5 wt % MgO and Mg-number of 61–67, indicating little fractionation. Olivine phenocrysts are as magnesian as Fo 93 and contain 0·4 wt % NiO; the Cr/(Cr + Al) values of spinels are mostly 0·5–0·8, indicating equilibrium with depleted mantle. There are three petrographic groups, based on phenocryst populations: (1) cpx–olivine basalt (COB); (2) plagioclase–olivine basalt (POB); (3) porphyritic basalt. Zr/Y and Nb/Yb are higher in POB (3·1–3·2 and 1·2–1·5, respectively) than in COB (Zr/Y = 2·8–3·0 and Nb/Yb = 0·7–0·9), suggesting that POB formed from lower degrees of mantle melting, or that the COB mantle source was more depleted. On the other hand, COB have Ba/Nb (70–80) and Th/Nb (0·4–0·5) that are higher than for POB (Ba/Nb = 30–35 and Th/Nb = 0·1–0·2), and also have steeper light rare earth element (LREE)-enriched patterns. Moreover, COB have enriched 87 Sr/ 86 Sr and 143 Nd/ 144 Nd, and higher Pb isotope values, suggesting that COB has a greater subduction component than POB. 176 Hf/ 177 Hf between COB and POB are similar and Hf behavior in COB and POB is similar to that of Zr, Y and HREE, suggesting that Hf is not included in the subduction component, which produced the differences between COB and POB. The calculated primary basaltic magmas of NW Rota-1 volcano (primary COB and POB magmas) indicate segregation pressures of 2–1·5 GPa (equivalent to 65–50 km depth). These magmas formed by 24–18% melting of mantle peridotite having Mg-number ~89·5. Diapiric ascent of hydrous peridotite mixed heterogeneously with sediment melts may be responsible for the NW Rota-1 basalts. These two basalt magma types are similar to those found at Sumisu and Torishima volcanoes in the Izu–Bonin arc, with COB representing wetter and POB representing drier magmas, where subduction zone-derived melt components are coupled with the water contents.

Description: The origin of andesite is an important issue in petrology because andesite is the main eruptive product at convergent margins, corresponds to the average crustal composition and is often associated with major Cu–Au mineralization. In this study we present petrographic, mineralogical, geochemical and isotopic data for basaltic andesites of the latest Pleistocene Pilavo volcano, one of the most frontal volcanoes of the Ecuadorian Quaternary arc, situated upon thick (30–50 km) mafic crust composed of accreted Cretaceous oceanic plateau rocks and overlying mafic to intermediate Late Cretaceous–Late Tertiary magmatic arcs. The Pilavo rocks are basaltic andesites (54–57·5 wt % SiO 2 ) with a tholeiitic affinity as opposed to the typical calc-alkaline high-silica andesites and dacites (SiO 2 59–66 wt %) of other frontal arc volcanoes of Ecuador (e.g. Pichincha, Pululahua). They have much higher incompatible element contents (e.g. Sr 650–1350 ppm, Ba 650–1800 ppm, Zr 100–225 ppm, Th 5–25 ppm, La 15–65 ppm) and Th/La ratios (0·28–0·36) than Pichincha and Pululahua, and more primitive Sr ( 87 Sr/ 86 Sr ~0·7038–0·7039) and Nd ( Nd ~ +5·5 to +6·1) isotopic signatures. Pilavo andesites have geochemical affinities with modern and recent high-MgO andesites (e.g. low-silica adakites, Setouchi sanukites) and, especially, with Archean sanukitoids, for both of which incompatible element enrichments are believed to result from interactions of slab melts with peridotitic mantle. Petrographic, mineral chemistry, bulk-rock geochemical and isotopic data indicate that the Pilavo magmatic rocks have evolved through three main stages: (1) generation of a basaltic magma in the mantle wedge region by flux melting induced by slab-derived fluids (aqueous, supercritical or melts); (2) high-pressure differentiation of the basaltic melt (at the mantle–crust boundary or at lower crustal levels) through sustained fractionation of olivine and clinopyroxene, leading to hydrous, high-alumina basaltic andesite melts with a tholeiitic affinity, enriched in incompatible elements and strongly impoverished in Ni and Cr; (3) establishment of one or more mid-crustal magma storage reservoirs in which the magmas evolved through dominant amphibole and clinopyroxene (but no plagioclase) fractionation accompanied by assimilation of the modified plutonic roots of the arc and recharge by incoming batches of more primitive magma from depth. The latter process has resulted in strongly increasing incompatible element concentrations in the Pilavo basaltic andesites, coupled with slightly increasing crustal isotopic signatures and a shift towards a more calc-alkaline affinity. Our data show that, although ultimately originating from the slab, incompatible element abundances in arc andesites with primitive isotopic signatures can be significantly enhanced by intra-crustal processes within a thick juvenile mafic crust, thus providing an additional process for the generation of enriched andesites.

Description: SUMMARY We analyse the spatial distribution of the intensity data points surveyed after the M w 6.3, 2009 L’Aquila (central Italy) earthquake, with the aim to recognize and quantify finite-fault and directivity effects. The study is based on the analysis of the residuals, evaluated with respect to attenuation-with-distance models, calibrated for L’Aquila earthquake. We apply a non-parametric approach considering both the epicentral and the rupture distance, which accounts for the finite extension of the source. Then, starting from a simplified kinematic rupture model of the L’Aquila fault, we compute four directivity predictors proposed in literature, and assess their correlation with intensity residuals. We derive a so-called Intensity Directivity Factor by the correlation between theoretical predictors and observed residuals that allows us to identify and quantify the intensity data points affected by forward and backward directivity during L’Aquila earthquake. We find that the effects are more pronounced in the forward directivity direction and increments up to 1 MCS intensity unit are expected. Moreover, the directivity predictor that accounts for radiation pattern poorly correlates with residuals. These results show that the spatial distribution of the L’Aquila macroseismic field is affected by source effects and in particular that directivity-induced amplification effects can be recognized. We show that the quasi-unilateral rupture propagation along the fault can explain the high-intensity patterns observed along specific direction at relatively large distance from the instrumental epicentre, in accordance with the seismological source models derived from the analysis of instrumental observations.

Description: SUMMARY Temporal variations in the elastic behaviour of the Earth’s crust can be monitored through the analysis of the Earth’s seismic response and its evolution with time. This kind of analysis is particularly interesting when combined with the reconstruction of seismic Green’s functions from the cross-correlation of ambient seismic noise, which circumvents the limitations imposed by a dependence on the occurrence of seismic events. In fact, because seismic noise is recorded continuously and does not depend on earthquake sources, these cross-correlation functions can be considered analogously to records from continuously repeating doublet sources placed at each station, and can be used to extract observations of variations in seismic velocities. These variations, however, are typically very small: of the order of 0.1 per cent. Such accuracy can be only achieved through the analysis of the full reconstructed waveforms, including later scattered arrivals. We focus on the method known as Moving-Window Cross-Spectral Analysis that has the advantage of operating in the frequency domain, where the bandwidth of coherent signal in the correlation function can be clearly defined. We investigate the sensitivity of this method by applying it to microseismic noise cross-correlations which have been perturbed by small synthetic velocity variations and which have been randomly contaminated. We propose threshold signal-to-noise ratios above which these perturbations can be reliably observed. Such values are a proxy for cross-correlation convergence, and so can be used as a guideline when determining the length of microseismic noise records that are required before they can be used for monitoring with the moving-window cross-spectral technique.

Description: SUMMARY The young (〈 5 Ma) Corinth Rift is an ideal natural laboratory to investigate rift deformation mechanisms by comparing extension rates determined by various methods spanning different depth and time ranges. Corinth Rift geodetic extension rates averaged over 5–100 yr have been interpreted to increase from ∼5 mm yr –1 or less in the east to 〉10–15 mm yr –1 in the west. We quantify total upper-crust and whole-crust extension on three profiles across the Corinth Rift. Whole-crust extension is greater across the central rift (∼11–21 km) than across the western part of the rift (∼5–13 km). This correlates with the overall rift morphology, which shows maximum basement subsidence, sediment accumulation, rift width and greatest summed Late Quaternary fault displacements in the central basin, but contrasts with the pattern of geodetic extension rates which are greater to the west of the central basin. The E–W increase in strain rates interpreted from geodetic data cannot have persisted over rift history to produce the observed rift morphology. We suggest the discrepancy between short-term and long-term extension patterns is related to shifts in the loci of maximum extension due to fault growth and linkage during Corinth Rift history, and is likely a characteristic of rift development in general. Total upper-crust and whole-crust extension estimates in the western rift, where extension estimates are best constrained, are within error. We propose that uniform pure-shear extension is a viable extension mechanism in the western rift and crustal extension estimates do not require the existence of a major active N–S dipping detachment fault.

Description: SUMMARY This paper describes an innovative method of volcano deformation measurements, applied to camera images taken from the 2004–2008 eruption period at Mount St. Helens. Dome growth was thought to be characterized by sustained, near-linear rates of a solid dacite plug. Through spatial digital image correlation (DIC) analysis of the camera images, new evidences arise that the deformation and strain rate of the spine was more complex. DIC yielded cumulative and incremental displacements, strain and shear planes at decimetre resolution. It was found that dome extrusion rates are highly non-linear, decelerating prior to partial collapse, followed by a pronounced dome extrusion increase and direction change. Associated processes have been identified through DIC, such as shallow landslides and reworking of talus apron material. The work highlights the strengths of camera strain monitoring, and illustrates that dome growth and collapse is a very dynamic process complexly interplaying with the surrounding.

Description: SUMMARY Direct detection of Rayleigh-wave azimuthal anisotropy is reported by applying an array analysis to broad-band seismic data in Southern California, USA. Our approach has excellent resolution for frequencies between 30 and 60 mHz and good resolution between 10 and 30 mHz. Limitation from array size limits accuracy below 10 mHz and complicated wave propagation effects lead to difficulty above 60 mHz. Between 30 and 60 mHz, azimuthal anisotropy of Rayleigh-wave phase velocity is detected cleanly. Phase velocity results show that the 2θ components dominate the 4θ components. Strength of anisotropy hovers around 1 per cent for the 2θ components with standard error ( 1 σ ) of about 0.4 per cent. Fast axes orientations show stable orientation in the direction 110° clockwise from north and its opposite direction 290° . Depth sensitivity kernels for P -wave and S -wave anisotropy indicate that anisotropy is confined to depths in the upper 100 km and may even be to shallower depths. We speculate that the fast axes alignment may be associated with a strong shearing process over the thickness of lithosphere due to obliquely collisional motions between the Pacific Plate and the North America Plate in the region.

Description: SUMMARY Zhou & Greenhalgh have recently presented an application of the Gaussian quadrature grid to seismic modelling in which the authors propose a meshing scheme that partitions the domain independently of the discontinuities in the media parameters. This comment aims to clarify the implications that this strategy has on the accuracy.

Description: SUMMARY We obtained the crustal thickness ( H ) and average Vp / Vs ratio (κ) for the central (CNCC) and western North China Craton (WNCC) by H –κ stacking of receiver functions. Our results show that H and κ varies significantly but similarly in the northern and southern CNCC, reflecting possible effects of the Phanerozoic cratonotic reactivation on the region. Of the WNCC, the south is featured by κ values typical for Precambrian shields (∼1.77) and variations of H mirroring surface topography; the north presents large fluctuations of κ (1.68–1.93) and a crustal thickening (up to 50 km) uncorrelated with topography. Such N–S differences were thought to be associated with a ∼1.95-Ga suturing event occurred in the north with little effects on the south. Our observations suggest that distinct structural heterogeneities of the NCC not only reflect diverse regional tectonics in Phanerozoic, but also provide evidence for long-term survival of crustal fabrics in stable continental interior.

Description: SUMMARY The aim of this paper is to present an alternative method for the analysis of resistivity data. The methodology was developed during a study to evaluate if electrical resistivity can be used as a tool for analysing subsurface gas dynamics and gas emissions from landfills. The main assumption of this study was that variations in time of resistivity data correspond to variations in the relative amount of gas and water in the soil pores. Field measurements of electrical resistivity, static chamber gas flux and weather data were collected at a landfill in Helsingborg, Sweden. The resistivity survey arrangement consisted of nine lines each with 21 electrodes in an investigation area of 16 ×20 m. The ABEM Lund Imaging System provided vertical and horizontal resistivity profiles every second hour. The data were inverted in Res3Dinv using L 1 -norm–based optimization method with a standard least-squares formulation. Each horizontal soil layer was then represented as a linear interpolated raster model. Different areas underneath the gas flux measurement points were defined in the resistivity model of the uppermost soil layer, and the vertical extension of the zones could be followed at greater depths in deeper layer models. The average resistivity values of the defined areas were calculated and plotted on a time axis, to provide graphs of the variation in resistivity with time in a specific section of the ground. Residual variation of resistivity was calculated by subtracting the resistivity variations caused by the diurnal temperature variations from the measured resistivity data. The resulting residual resistivity graphs were compared with field data of soil moisture, precipitation, soil temperature and methane flux. The results of the study were qualitative, but promising indications of relationships between electrical resistivity and variations in the relative amount of gas and water in the soil pores were found. Even though more research and better data quality is necessary for verification of the results presented here, we conclude that this alternative methodology of working with resistivity data seems to be a valuable and flexible tool for this application.

Description: SUMMARY We analysed high-frequency body waves of local earthquakes to image the damage zone of the Calico fault in the eastern California shear zone. We used generalized ray theory and finite difference methods to compute synthetic seismograms for a low-velocity fault zone (FZ) to model the direct and FZ-reflected P and S traveltimes of local earthquakes recorded by a temporary array across the fault. The low velocity zone boundaries were determined by apparent traveltime delays across the fault. The velocity contrast between the fault zone and host rock was constrained by the traveltime delays of P and S waves and differential traveltimes between the direct and FZ-reflected waves. The dip and depth extent of the low velocity zone were constrained by a systematic analysis of direct P traveltimes of events on both sides of the fault. We found that the Calico fault has a ∼1.3-km-wide low velocity zone in which the P - and S -wave velocity decreased 40 and 50 per cent, respectively, with respect to the host rock. The low velocity zone dips 70° northeast and extends 3 km in depth.

Description: SUMMARY High-grade metamorphic rocks of the upper amphibolite to granulite facies in the Higher Himalayan crystalline (HHC) were studied for their appropriateness for palaeomagnetic investigations. Sampling and fieldwork was accomplished on a south–north transect at ∼87°E in the Khumbu region south of Everest. Paragenesis of iron-sulphides and oxides was found to be variable due to different source rocks and metamorphic conditions in this area. Unblocking temperature ranges, isothermal remanence acquisition, low-temperature magnetic properties and microscopic observations confirm the presence of pyrrhotite and magnetite. Two remanence components were separated with unblocking temperature intervals of 250–330 °C and 370–550 °C, which can be attributed to characteristic remanent magnetizations (ChRM) of monoclinic pyrrhotite (ChRM pyr ) and magnetite (ChRM mag ), respectively. Mainly the ChRM pyr showed a consistent demagnetization behaviour, probably residing in pyrrhotite particles within the single domain to pseudo-single domain range. Normal and reverse polarities were observed, partly co-existing within the unblocking spectrum of single specimens; this (i) supports a thermoremanent (TRM) origin of the ChRM pyr and (ii) indicates that the time span of remanence acquisition was long enough to average out palaeosecular variation. The age of the ChRM pyr can be related to the last metamorphic cooling event at the Early to Late Miocene boundary (16 Ma). The ChRM mag is probably a secondary chemical or thermochemical remanence acquired by hydrothermal processes during retrogression. ChRM pyr and ChRM mag directions are similar suggesting that pyrrhotite and magnetite were magnetized approximately contemporaneously. Separation of the sampling area into three sectors reveals a south–north trend with (i) variable block rotations in the southernmost part close to the Main Central Thrust (five sites), (ii) no rotation in the area between Lukla and Namche (10 sites) and (iii) a significant clockwise rotation ( N = 5 sites, D / I = 42.3°/29.4°, k = 20.7, α 95 = 17.2°) along the Imja Khola north of Namche. We interpret the clockwise rotation as an apparent rotation due to tilting by crustal doming in the Everest region. This doming obviously occurred north of a fault zone located at around Namche and was not affecting the region south of it. The clockwise sense of the apparent vertical-axis rotation indicates that the area north of Namche is located at the western side of the dome structure where westward tilting around a subhorizontal axis occurred (best fit: tilt direction 277°, tilt angle 51°).

Description: SUMMARY The deep structure of the Upper Rhine Graben (URG), a continental rift in SW Germany and E France, is still poorly known. This deficit impedes a full understanding of the geodynamic evolution of this prominent rift. We study the lithosphere–asthenosphere structure using teleseismic waveforms obtained from the passive broad-band TIMO project across the central URG. The recovered, crust-corrected traveltime residuals relative to the iasp91 earth model are tiny (mostly less than 0.2–0.3 s). The average measured slowness (〈1 s deg –1 ) and backazimuth (〈5°) deviations are also very small and do not show any systematic wave front anomalies. These observed perturbation values are smaller than expected ones from synthetic 3-D ray tracing modelling with anomalies exceeding 2–3 per cent seismic velocity in the mantle. Thus there is no significant hint for any deep-seated anomaly such as a mantle cushion, etc. This result means that the rifting process did not leave behind a lower lithospheric signature, which could be clearly verified with high-resolution teleseismic experiments. The only significant traveltime perturbation at the central URG is located at its western side in the upper crust around a known geothermal anomaly. The upper crustal seismic anomaly with traveltime delays of 0.2–0.3 s cannot be explained with increased temperature alone. It is possibly related to a zone of highly altered granite. In the west of our network a traveltime anomaly (0.6–0.7 s delay) related with the Eifel plume is confirmed by the TIMO data set.

Description: Mid-ocean ridge basalts (MORB) from the Arctic Ocean have been significantly less studied than those from other oceans. The Arctic ridges (Gakkel Ridge and Lena Trough) are ultraslow-spreading ridges with low melt productivity and are thus the best locations to investigate mantle heterogeneity. We report the major and trace element and Sr–Nd–Pb–Hf isotope compositions of basalts generated along the Lena Trough and the westernmost part of the Gakkel Ridge in the Arctic Ocean. Basalts from the northern Lena Trough and westernmost Gakkel Ridge (NLT–WGR) have compositions close to normal MORB. The geochemical composition of the NLT–WGR lavas confirms a binary mixing model involving melts from a depleted MORB mantle source and a Spitsbergen amphibole-bearing subcontinental lithospheric mantle (SCLM) source. In contrast, in the central part of the Lena Trough (CLT), the basalts are alkalic with relatively high Mg-number (60–65), high SiO 2 (51·0–51·6 wt %), Al 2 O 3 (18·1–18·4 wt %), Na 2 O (4·0–4·2 wt %), K 2 O (1·0–1·6 wt %), K 2 O/TiO 2 (0·6–0·9) and (La/Sm) PM (1·4–1·8), and low FeO (6·5–6·8 wt %) contents. These basalts display isotope variations with 87 Sr/ 86 Sr ranging from 0·70361 to 0·70390, 143 Nd/ 144 Nd from 0·51283 to 0·51290 ( Nd + 3·7 to +5·2), 176 Hf/ 177 Hf from 0·28313 to 0·28322 ( Hf + 11·6 to +14·9) and 206 Pb/ 204 Pb from 17·752 to 17·884, 207 Pb/ 204 Pb from 15·410 to 15·423 and 208 Pb/ 204 Pb from 37·544 to 37·670. These isotope compositions clearly distinguish the CLT lavas from those generated along the Gakkel Ridge. For the CLT lavas, involvement of a phlogopite- or amphibole- and (possibly garnet)-bearing SCLM source component is proposed. Owing to SCLM contamination along the entire length of the Lena Trough, we classify the Lena Trough as an ocean–continent transition boundary. Magmatism similar to that observed in the Lena Trough would be expected to occur wherever ocean spreading initiates.

Description: SUMMARY Recent revisions to the satellite-derived vertical gravity gradient (VGG) data reveal more detail of the ocean bottom and have allowed us to develop a non-linear inversion method to detect seamounts in VGG data. We approximate VGG anomalies over seamounts as sums of individual, partially overlapping, elliptical polynomial functions, which allows us to form a non-linear inverse problem by fitting the polynomial model to the observations. Model parameters for a potential seamount include geographical location, peak VGG amplitude, major and minor axes of the elliptical base, and the azimuth of the major axis. The non-linear inversion is very sensitive to the initial values for the location and amplitude; hence, they are constrained by the centre and amplitude of the uppermost contours obtained with a 1-Eötvös contour interval. With these initial conditions from contouring, we execute a step-wise and fully automated inversion and obtain optimal model estimates for potential seamounts; these are statistically evaluated for significance using the Akaike Information Criterion and  F  tests. A logarithmic barrier technique is applied to ensure positivity of all seamount amplitudes. After automatic and manual inspections of the model parameters we estimate actual heights and basal ellipses of the inspected potential seamounts directly from the predicted bathymetry grid. In this study, we find globally 24 643 potential seamounts ( h  ≥ 0.1 km) that are located away from continental margins; 8458 potential seamounts are taller than 1 km. Although our global estimate is significantly lower than predictions from previous studies, a first-order reconciliation of the size-frequency statistics obtained from those studies reveals that the previous counts are systematically overestimated. Because of the ambiguity of gravity signals due to small seamounts of  h  〈 1 km and the overlap with abyssal hills, we estimate the global seamount census to lie in the 40 000–55 000 range. The seamount data from this study are accessible from http://www.soest.hawaii.edu/PT/SMTS .

Description: SUMMARY The growth of the Earthscope/USArray Transportable Array (TA) has prompted the development of new methods in surface wave tomography that track phase fronts across the array and map the traveltime field for each earthquake or for each station from ambient noise. Directionally dependent phase velocities are determined locally by measuring the gradient of the observed traveltime field without the performance of a formal inversion. This method is based on the eikonal equation and is, therefore, referred to as ‘eikonal tomography’. Eikonal tomography is a bent-ray theoretic method, but does not account for finite frequency effects such as wave interference, wave front healing, or backward scattering. This shortcoming potentially may lead to both systematic bias and random error in the phase velocity measurements, which would be particularly important at the longer periods studied with earthquakes. It is shown here that eikonal tomography can be improved by using amplitude measurements to construct a geographically localized correction via the Helmholtz equation. This procedure should be thought of as a finite-frequency correction that does not require the construction of finite-frequency kernels and is referred to as ‘Helmholtz tomography’. We demonstrate the method with Rayleigh wave measurements following earthquakes between periods of 30 and 100 s in the western US using data from the TA. With Helmholtz tomography at long periods (〉50 s): (1) resolution of small-scale isotropic structures, which correspond to known geological features, is improved, (2) uncertainties in the isotropic phase velocity maps are reduced, (3) the directionally dependent phase velocity measurements are less scattered, (4) spurious 1-psi azimuthal anisotropy near significant isotropic structural contrasts is reduced, and (5) estimates of 2-psi anisotropy are better correlated across periods.

Description: SUMMARY Geomagnetic polarity transitions may be significantly more complex than are currently depicted in many sedimentary and lava-flow records. By splicing together paleomagnetic results from earlier studies at Steens Mountain with those from three newly studied sections of Oregon Plateau flood basalts at Catlow Peak and Poker Jim Ridge 70–90 km to the southeast and west, respectively, we provide support for this interpretation with the most detailed account of a magnetic field reversal yet observed in volcanic rocks. Forty-five new distinguishable transitional (T) directions together with 30 earlier ones reveal a much more complex and detailed record of the 16.7 Ma reversed (R)-to-normal (N) polarity transition that marks the end of Chron C5Cr. Compared to the earlier R-T-N-T-N reversal record, the new record can be described as R-T-N-T-N-T-R-T-N. The composite record confirms earlier features, adds new west and up directions and an entire large N-T-R-T segment to the path, and fills in directions on the path between earlier directional jumps. Persistent virtual geomagnetic pole (VGP) clusters and separate VGPs have a preference for previously described longitudinal bands from transition study compilations, which suggests the presence of features at the core–mantle boundary that influence the flow of core fluid and distribution of magnetic flux. Overall the record is consistent with the generalization that VGP paths vary greatly from reversal to reversal and depend on the location of the observer. Rates of secular variation confirm that the flows comprising these sections were erupted rapidly, with maximum rates estimated to be 85–120 m ka −1 at Catlow and 130–195 m ka −1 at Poker Jim South. Paleomagnetic poles from other studies are combined with 32 non-transitional poles found here to give a clockwise rotation of the Oregon Plateau of 11.4°± 5.6° with respect to the younger Columbia River Basalt Group flows to the north and 14.5°± 4.6° with respect to cratonic North America (95 per cent confidence interval).

Description: SUMMARY The principal aim of performing a survey or experiment is to maximize the desired information within a data set by minimizing the post-survey uncertainty on the ranges of the model parameter values. Using Bayesian, non-linear, statistical experimental design (SED) methods we show how industrial scale amplitude variations with offset (AVO) surveys can be constructed to maximize the information content contained in AVO crossplots, the principal source of petrophysical information from seismic surveys. The design method allows offset dependent errors, previously not allowed in non-linear geoscientific SED methods. The method is applied to a single common-midpoint gather. The results show that the optimal design is highly dependent on the ranges of the model parameter values when a low number of receivers is being used, but that a single optimal design exists for the complete range of parameters once the number of receivers is increased above a threshold value. However, when acquisition and processing costs are considered we find that a design with constant spatial receiver separation survey becomes close to optimal. This explains why regularly-spaced, 2-D seismic surveys have performed so well historically, not only from the point of view of noise attenuation and imaging in which homogeneous data coverage confers distinct advantages, but also to provide data to constrain subsurface petrophysical information.

Description: SUMMARY Tsunami generated by the Late Bronze Age (LBA) eruption of Thera were simulated using synthetic tide records produced for selected nearshore (∼20 m depths) sites of northern Crete, the Cyclades Islands, SW Turkey and Sicily. Inundation distances inland were also calculated along northern Crete. Modelling was performed by incorporating fully non-linear Boussinesq wave theory with two tsunamigenic mechanisms. The first involved the entry of pyroclastic flows into the sea, assuming a thick (55 m; 30 km 3 ) flow entering the sea along the south coast of Thera in three different directions all directed towards northern Crete, then a thin pyroclastic flow (1 m; 1.2 km 3 ) entering the sea along the north coast of Thera directed towards the Cyclades Islands. Flows were modelled as a solid block that slowly decelerates along a horizontal surface. The second mechanism assumed caldera collapse, of 19 km 3 and 34 km 3 modelled as a dynamic landslide producing rapid vertical displacements. Calculated nearshore wave amplitudes varied from a few metres to 28 m along northern Crete from pyroclastic flows, and up to 19 m from caldera collapse (34 km 3 volume). Inundation distances on Crete were 250–450 m. Waves produced by pyroclastic flows were highly focused, however, as a function of sea entry direction. Smaller volume pyroclastic flows produced nearshore wave amplitudes up to 4 m in the Cyclades islands north of Thera. Wave amplitudes in the Cyclades from smaller volume caldera collapse (19 km 3 ) were up to 24 m, whereas in SW Turkey were as low as 2.1 and 0.8 m (Didim and Fethye where LBA tsunami deposits have been found). Wave amplitudes for the larger volume caldera collapse (34 km 3 ) were generally 2.5–3 times larger than those generated by the smaller volume collapse (19 km 3 ). These results provide estimates for understanding possible consequences of tsunami impact in LBA coastal zones, thus providing criteria at archaeological sites for detecting inundation damage, as well as for contemporary hazard assessment; they also provide additional criteria for deciphering homogenite layers in the abyssal stratigraphy of the Ionian and eastern Mediterranean Seas.

Description: SUMMARY The first archaeomagnetic secular variation (SV) curves for the whole Southern Balkan Peninsula are presented. These are based on all data within a 700 km circle centred at Thessaloniki (40.60 o N, 23.00 o E). This data set consists of 325 directional and 625 intensity data mainly from Greece, Bulgaria, Serbia and southern Hungary. Some data from southern Italy are also included. The sliding moving window technique, was used to calculate a continuous SV curve for intensity while the directional SV curves were calculated using the bivariate extension of the Fisher statistics. These curves are well constrained and clearly show the main features of the geomagnetic field variation in this region during the last eight millennia. Comparisons with the predictions of the SCHA.DIF.3K and SCHA.DIF.8K regional and the CALS7K.2 and ARCH3K.1 global geomagnetic field models show a good agreement for the last 3000 years but differences for older times. The Balkan SV curves identify several rapid changes of the geomagnetic field in eastern Europe and can be used as reference curves for archaeomagnetic dating in the Balkan Peninsula.

Description: SUMMARY Receiver function analysis reveals a 7 km step-like change in crustal thickness across the Taranaki–Ruapehu Line (TR-line) of North Island, New Zealand. The TR-line runs east–west between the active andesite volcanoes of Mt Taranaki and Ruapehu and marks the southern-most extent of subduction zone volcanism in New Zealand. North of the TR-line receiver functions show a strong and sharp P -to- S ( Ps ) conversion at 25 ± 1.5 km depth, which is interpreted as a shallow Moho. At the TR-line the Moho Ps conversion deepens across a step to ∼32 km depth and weakens. Further south the Moho deepens to 〉35 km. Most of the 7 km step-change in crustal thickness occurs over a lateral distance of ∼8 km, yet there is little surface or topographic manifestation of the line. Given Fresnel zone considerations, the dip of the Moho offset could vary between 90° and 45°. Gravity, seismic attenuation and electrical data all show that the TR-line is not only a step in crustal thickness but also a profound lithospheric boundary as mantle properties, such as attenuation ( Q p −1 ), implied density, and electrical resistivity change abruptly across the line. An east–west oriented cluster of earthquakes with hypocentral depths of 20–40 km is centred on the Moho step. We propose that the Moho step, the earthquakes and the rapid change in mantle properties across the TR boundary are causally related. Processes that could be responsible for the phenomena described here include the rapid removal of mantle lithosphere and lower crust to the north of the TR-line.

Description: SUMMARY We investigate the slip distribution of the 2008 May 12 Wenchuan  M s 7.9 earthquake using GPS data and InSAR measurements under a linear inversion scheme, with emphasis on the effect of three factors, including constraint on rake, different discretizations, and layered elastic model. Within our inversion parameterization context, we find the most influential factor would be constraint on rake. Without constraint on rake, the slip model seems physically wrong under the depth of 15 km, due to the limited depth resolution of the geodetic data used, especially the one orbit of InSAR measurements. Thus it is necessary to add  a priori  to the slip rake to obtain a reasonable fault source model. Different discretizations of the subfault patches have a notable impact on the slip distribution. Also, the layered elastic model predicts more slip at depth than does the half-space model, by about 15–20 per cent. The characteristics of slip distribution established through our inversions include some points as follows: (1) the most reliable results would be expected at depth range of 0–15 km; (2) three peak-slip asperities are inverted, at Yingxiu county (near epicentre), Yuejiashan county, and Beichuan county, respectively; (3) the inverted rake along the Yingxiu-Beichuan fault changes from dominantly thrusting motion at the southwest segments to dominant right-lateral or even pure right-lateral strike slip at the northeast segments, and (4) only thrust slip is occurred on the Guanxian-Jiangyou fault.

Description: SUMMARY Recordings from a temporary seismic network in the Rwenzori region of the East African Rift show seven earthquakes at a depth range between 53 and 60 km. Evidence from receiver functions indicates a crustal thickness of about 32 km in this region suggesting that the earthquakes are located well within the mantle beneath the rift. The occurrence of faulting at this depth is of relevance for models of continental-rift initiation and related deep magmatic processes.

Description: SUMMARY We present an analysis on thermal properties, density and porosity of clastic, chemical/biochemical and intrasedimentary volcanic rocks collected from petroleum exploration wells of the Po Basin (Northern Italy). Moreover, we investigate the applicability of the Hashin–Shtrikman's model for a mineral aggregate in combination with the Zimmerman's model that takes pore shape into account, to calculate the bulk thermal conductivity. In case of macroscopically isotropic rocks, deviations between predicted and measured values range from −2.2 per cent to 6.9 per cent, and significantly decrease if a proper pore aspect ratio is chosen. Regarding the volumetric heat capacity, approximate estimates were obtained by means of a weighted average of the volumetric heat capacity of the mineral grains and that of the pore-filling water. The differences between the computed and measured values range from –6.2 per cent to 4.9 per cent and, on average, the computed volumetric heat capacity is lower by 1.6 per cent. The water loss during compaction and the temperature increase with depth are main factors controlling thermal properties. An anisotropy effect occurs in the case of rocks rich in sheet silicates. Due to rotation of these minerals, the vertical thermal conductivity of sheet silicates decreases exponentially with burial depth from 2.13 W m −1 K −1 at the surface to 0.52 W m −1 K −1 at 4.5 km. The laboratory data allow the formulation of compaction curves for the different sedimentary rock types. Examples of estimations of in situ vertical thermal conductivity and heat flow are finally given for two petroleum wells at which lithostratigraphic information is known in good detail.

Description: SUMMARY We present forward and adjoint spectral-element simulations of coupled acoustic and (an)elastic seismic wave propagation on fully unstructured hexahedral meshes. Simulations benefit from recent advances in hexahedral meshing, load balancing and software optimization. Meshing may be accomplished using a mesh generation tool kit such as CUBIT, and load balancing is facilitated by graph partitioning based on the SCOTCH library. Coupling between fluid and solid regions is incorporated in a straightforward fashion using domain decomposition. Topography, bathymetry and Moho undulations may be readily included in the mesh, and physical dispersion and attenuation associated with anelasticity are accounted for using a series of standard linear solids. Finite-frequency Fréchet derivatives are calculated using adjoint methods in both fluid and solid domains. The software is benchmarked for a layercake model. We present various examples of fully unstructured meshes, snapshots of wavefields and finite-frequency kernels generated by Version 2.0 ‘Sesame’ of our widely used open source spectral-element package SPECFEM3D.

Description: SUMMARY Allègre et al. recently presented new experimental data regarding the dependence of the streaming potential coupling coefficient with the saturation of the water phase. Such experiments are important to model the self-potential response associated with the flow of water in the vadose zone and the electroseismic/seismoelectric conversions in unsaturated porous media. However, the approach used to interpret the data is questionable and the conclusions reached by Allègre et al. likely incorrect

Description: SUMMARY Revil and Linde recently commented our paper concerning streaming potential (SP) measurements in unsaturated sand during drainage experiments. This comment suggests that the approach used to infer SP coefficients was inappropriate for unsaturated conditions, and therefore yielded wrong conclusions and ‘unphysical’ results regarding the behaviour of the relative SP coefficient. This reply argues that even if in Allègre et al . we neglected some secondary electrokinetic sources, the resulting conclusions are still representative of the behaviour of the true SP coefficient, and that the remarks of Revil & Linde arose from a misunderstanding of the drainage experiment conditions. We also find support for our results from a comparison between our observations and previous experimental studies.

Description: SUMMARY A small-aperture (40 m) short-period seismic array was installed during four months on the Vercors massif (Western French Alps) at the top of a limestone column which collapsed one month and a half later. During this monitoring period, 193 seismic events were recorded by the seven seismometers of the array. Signal analysis yielded three main types of nearby seismic events to be identified from temporal and spectral characteristics: microearthquakes (single or multiple events), individual block falls and rock falls. It turned out that 60 per cent of the 193 events were classified as microearthquakes, exhibiting distinct  P  and  S  waves, while 17 per cent of these events remained unclassified. Out of the microearthquakes, 40 events with a good signal-to-noise ratio were selected and processed.  P-  and  S -wave traveltimes were picked on the records and the inferred hypocentral distances agree with the two zones of the scarp exhibiting fresh ruptures after the fall. Polarization analysis of the 3-C records, along with numerical simulations, allowed discriminating between the two possible rupture mechanisms (toppling and sliding). Shear rupture (sliding) was the predominant mode in the lower part of the column whereas traction rupture (toppling) affected the upper part. Finally, the comparison between the ground motions recorded on the column and on the rock mass showed a systematic amplification on the column. Signal processing and numerical modelling both suggest that this amplification resulted from the excitation of the natural frequencies of the column and is particularly high (〉3) for microearthquakes occurring at the column-to-mass interface.

Description: SUMMARY The kinematic fractal source model presented in this study is able to simulate broad-band accelerograms with spectral amplitudes proportional to a fraction of the directivity coefficient C d in the far-field approximation. This approach is based on a composite source description, where subevents are generated using a fractal distribution of sizes, and, by summation, produces k -square space distribution of the slip. Each elementary source is described as a crack-type slip model growing circularly from a nucleation point when the rupture front reaches it. In order to better control the directivity effect, the location of the nucleation point for an elementary source is assumed to be scale-dependent. For the larger sources, the nucleation point is located near the intercept of the crack with the rupture front, whereas for smaller sources, it is randomly chosen within the crack. For simplicity, a constant rupture velocity is assumed. Each subevent is set up with scale-dependent rise-time, assuming a boxcar source–time function, hence filtering out its own high frequency radiation. The resulting mean slip-velocity functions are very similar to the ones derived from dynamic rupture modelling. Ground motion synthetics are computed by convolving the slip-velocity functions with the Green's functions. It is demonstrated that, in the far-field approximation, accelerogram spectra follow the ω 2 model with amplitudes controlled by the frequency-dependent directivity effects. In particular, spectral amplitudes at high-frequencies are proportional to a fraction of C d . These results were verified for few earthquake magnitudes. In addition, a validation exercise was made in the near-fault region by modelling the complete wavefield of strong ground motion at a few receivers and for several rupture scenarios. The synthetic strong-motion parameters are compared to the ones predicted by empirical attenuation relationships. It is shown that calculated standard deviations are in good agreement with the empirical ones, as well as the ground-motion parameter amplitudes predicted as a function of distance for whole interval of source distance considered in modelling. Minor differences were found in peak ground-accelerations computed at large distance from the fault, a problem related to the simplified response of the medium.

Description: SUMMARY Geodynamic models incorporating metamorphic phase transformations almost invariably assume the validity of the Boussinesq approximation that violates conservation of mass. In such models metamorphic density changes take place without volumetric effects. We assess the impact of the Boussinesq approximation by comparing models of orogeny accompanied by lower crustal eclogitization both with and without the approximation. Our results demonstrate that the approximation may cause errors approaching 100 per cent in characteristic measures of orogenic shape. Mass conservation errors in Boussinesq models amplify with model time. Mass conservative models of metamorphism are therefore essential to understand long-term tectonic evolution and to assess the importance of the different geodynamic processes.

Description: The Western Alpine Sesia–Lanzo Zone (SLZ) is a sliver of eclogite-facies continental crust exhumed from mantle depths in the hanging wall of a subducted oceanic slab. Eclogite-facies felsic and basic rocks sampled across the internal SLZ show different degrees of retrograde metamorphic overprint associated with fluid influx. The weakly deformed samples preserve relict eclogite-facies mineral assemblages that show partial fluid-induced compositional re-equilibration along grain boundaries, brittle fractures and other fluid pathways. Multiple fluid influx stages are indicated by replacement of primary omphacite by phengite, albitic plagioclase and epidote as well as partial re-equilibration and/or overgrowths in phengite and sodic amphibole, producing characteristic step-like compositional zoning patterns. The observed textures, together with the map-scale distribution of the samples, suggest open-system, pervasive and reactive fluid flux across large rock volumes above the subducted slab. Thermodynamic modelling indicates a minimum amount of fluid of 0·1–0·5 wt % interacting with the wall-rocks. Phase relations and reaction textures indicate mobility of K, Ca, Fe and Mg, whereas Al is relatively immobile in these medium-temperature–high-pressure fluids. Furthermore, the thermodynamic models show that recycling of previously fractionated material, such as in the cores of garnet porphyroblasts, largely controls the compositional re-equilibration of the exhumed rock body.

Description: SUMMARY We present a new 3-D transversely isotropic shear wave velocity model of the European and Mediterranean upper mantle obtained by analysis of surface waves. Data used are fundamental-mode Rayleigh and Love group velocity measurements in the period range 35–170 s, taken on seismograms recorded by European stations for regional earthquakes. The tomographic inversion to map the 3-D earth structure is split into two steps. First, we regionalize the group velocity dispersion measurements, obtaining distinct geographical group velocity maps at different periods; then, each local dispersion curve is inverted separately to find the shear wave velocity structure at depth. The inversion benefits from using a priori information from a 3-D global mantle model (S20RTS) and a new detailed European crustal model (EPcrust) to constrain the shallower layers. The inversion scheme follows a non-linear iterative algorithm by which Rayleigh and Love group slowness are inverted simultaneously for the best-fitting isotropic Voigt shear wave speed and radial anisotropy parameter ( v SH − v SV ) . Final merging of the v S profiles results in a new higher resolution 3-D model of European upper mantle. We find that Western Europe and Mediterranean Sea are mainly characterized by relatively low velocities, strongly contrasting with the fast roots of the Eastern European Craton. Many regional scale structures are also evident in the model, thus providing insights into the complex geodynamic framework of the European continent. Most prominent are the low-velocity West Mediterranean spreading basins and European Cenozoic rift system, and seismically fast features connected to subduction of Adria microplate, Hellenic Arc and Calabrian Arc. Radial anisotropy does not vary very significantly with respect to the PREM profile, as available data only resolve lateral variations to a limited degree due to trade-off with velocity. EPmantle has the potential to provide a reliable seismological reference for the upper-mantle structure in the broad European region.

Description: SUMMARY The main objective of this paper is to present a methodology for predicting spectrum compatible seismograms at sites where main shock is not recorded by other sites from which the main event has been observed using recorded aftershocks. The genetic algorithm and multitaper method are employed to minimize differences between the calculated synthetic seismogram and observed data. The multitaper spectrum approach being used permits the reduction of spectral leakage providing the possibility of matching the Fourier spectral amplitudes corresponding to synthetic and observed waveforms. The proposed technique is applied to the 2006 Silakhor earthquake ( M w = 6.1, Iran) as a case study. The six component waveforms (L, T and V) at two stations were simulated by minimizing errors between the synthetic response spectra, the Fourier spectral amplitudes and those of the observed. The results are shown to be in a good agreement with those of the observed data. The three components of waveforms at the other three stations were predicted incorporating the calculated parameters. The validity of the technique is proven through demonstrating good agreement of the response spectra and Fourier spectral amplitudes corresponding to the predicted and observed seismograms at other three stations. It is concluded that, the methodology is applicable for predicting realistic acceleration time histories to be used in seismic performance evaluation of existing structures in the region under study.

Description: SUMMARY Crustal density models derived from seismic velocity models by means of velocity–density conversions typically reproduce the main features of the observed gravity anomaly over the area but often show significant misfits. Given the uncertainty in the relationship between velocity and density, seismically derived density models should be regarded as an initial estimate of the true subsurface density structure. In this paper, we present a method for estimating the adjustments necessary to a seismically derived density model to improve the fit to gravity data. The method combines the Genetic Algorithm paradigm with linear inversion as a way to approach the non-linear and linear aspects of the problem. The models are divided into three layers representing the sedimentary column, the crystalline crust and the lithospheric mantle; the depths of these layers are determined from the seismic velocity model. Each of the layers is divided into a number of provinces and a density adjustment (Δρ) value is found for each province so that the residual gravity (difference between the observed gravity anomaly and the anomaly calculated for the seismically derived model) is minimized while keeping Δρ between predefined bounds. The preferred position of the province boundaries is found through the artificial evolution of a population of solutions. Given the stochastic nature of the algorithm and the non-uniqueness of the problem, different realizations can yield different solutions. By performing multiple realizations we can analyse a set of solutions by taking their mean and standard deviation, providing not only an estimate of the Δρ distribution in the subsurface but also an estimate of the associated uncertainty. Synthetic tests prove the ability of the algorithm to accurately recover the location of province boundaries and the Δρ values for a known model when using noise-free synthetic data. When noise is added to the data, the algorithm broadly recovers the features that define the known model despite greater standard deviations of the solutions and the occurrence of artefacts in the mean solutions. The algorithm was applied to four profiles across the Caribbean–South America Plate boundary. Some general patterns in the distribution of Δρ were observed consistently in the profiles and are correlated with the interpretations of the velocity models. Positive Δρ values in the sedimentary layer, negative Δρ values for island arc and extended island arc crust with an abrupt change to positive values in South American crust, and positive values in the mantle under the continent and island arc with a transition to negative values under the Caribbean oceanic crust.

Description: SUMMARY Active high-temperature (〉150 °C) geothermal areas like the Krafla caldera, NE-Iceland, often show distinct magnetic lows in aeromagnetic anomaly maps suggesting a destruction of magnetic minerals by hydrothermal activity. The main alteration processes in such an environment are low-temperature oxidation (〈350 °C, maghemitization) and fluid–rock interactions. We investigated the rock magnetic properties [natural remanent magnetization (NRM) magnetic susceptibility and their temperature and field variation] and the mineralogy, using X-ray diffraction, microscopic methods and electron microprobe analyses, of two drill cores (KH1 and KH3) from the rim of the Krafla caldera. The drill cores have distinctly lower NRM values (average 〈3 A m –1 ) compared to younger surface basalts (average 20 A m –1 ) along with a large variation in magnetic susceptibility (1.3 × 10 −7 – 4.9 × 10 −5 m 3 kg –1 ). The secondary mineral assemblage (sulphides, sphene, rutile and chlorite) indicates an alteration within the chlorite–smectite zone for both cores without depth zoning. Optical miscroscopy in combination with the Bitter technique and backscatter electron microscopy along with the thermomagnetic analyses allow distinguishing two different magnetomineralogical groups of titanomaghemite: (1) titanomaghemite with intermediate titanium concentration and probably high vacancy concentration, and (2) titanomaghemite with low titanium concentration and low vacancy concentration. The mineral assemblages, textures and magnetic properties deduced from the mentioned magnetic measurements indicate two-stage transformation mechanism: (1) Dissolution of titanium at low pH under oxidizing conditions. The ulvöspinel component of titanomagnetite and ilmenite forms rutile or sphene, and Fe 2 + migrates out of the spinel lattice forming titanomaghemite. (2) Formation of pyrite and dissolution of remaining titanomaghemite under reducing and acidic conditions. The latter mechanism produces ghost textures (all titanomaghemite is transformed and only their former grain shapes are preserved), with only paramagnetic minerals left and ferrimagnetic minerals nearly dissolved. This mechanism could explain the significant magnetization loss, which is seen in many local magnetic anomaly lows within the oceanic crust and volcanic islands like Iceland or Hawaii. The production of nanoporous textures in titanomaghemites is suggested as a mechanism for the enhancement of the magnetic susceptibility values related to the hydrothermal alteration of Krafla.

Description: SUMMARY We determine centroid-moment tensor (CMT) solutions by minimizing waveform differences between observed and simulated seismograms based on an adjoint method. Synthetic seismograms and Fréchet derivatives are calculated based on a spectral-element method. The non-linear adjoint CMT inversion algorithm requires three simulations for each iteration: one ‘forward’ simulation to obtain synthetics for the current source parameters, one ‘adjoint’ simulation which involves injecting time-reversed differences between observed and simulated seismograms as simultaneous virtual sources at each of the receivers, and an extra forward simulation to compute the step length in the conjugate-gradient direction. Whereas the vertical component of the adjoint wavefield reflects the radiation pattern near the centroid location, the components of the adjoint strain tensor capture the elements of the moment tensor. We use the method to determine adjoint CMT solutions for two representative southern California earthquakes using recent 3-D crustal model CVM-6.2. The adjoint CMT solutions are in good agreement with classical Hessian-based CMT solutions involving 3-D Green's functions. In general, adjoint CMT inversions require fewer numerical simulations than traditional Hessian-based inversions. This faster convergence holds promise for multiple moment-tensor and kinematic rupture inversions in 3-D earth models.

Description: SUMMARY We present a parallel goal-oriented adaptive finite element method that can be used to rapidly compute highly accurate solutions for 2.5-D controlled-source electromagnetic (CSEM) and 2-D magnetotelluric (MT) modelling problems. We employ unstructured triangular grids that permit efficient discretization of complex modelling domains such as those containing topography, dipping layers and multiple scale structures. Iterative mesh refinement is guided by a goal-oriented error estimator that considers the relative error in the strike aligned fields and their spatial gradients, resulting in a more efficient mesh refinement than possible with a previous approach based on the absolute errors. Reliable error estimation is accomplished by a dual weighted residual method that is carried out via hierarchical basis computations. Our algorithm is parallelized over frequencies, wavenumbers, transmitters and receivers, where adaptive refinement is performed in parallel on subsets of these parameters. Mesh sharing allows an adapted mesh generated for a particular frequency and wavenumber to be shared with nearby frequencies and wavenumbers, thereby efficiently reducing the parallel load of the adaptive refinement calculations. We demonstrate the performance of our algorithm on a large cluster computer through scaling tests for a complex model that includes strong seafloor topography variations and multiple thin stacked hydrocarbon reservoirs. In tests using up to 800 processors and a realistic suite of CSEM data parameters, our algorithm obtained run-times as short as a few seconds to tens of seconds.

Description: SUMMARY Geomagnetic induction responses such as geomagnetic depth sounding (GDS), magnetotelluric (MT), and horizontal geomagnetic transfer function (HTF) at long periods are used to estimate the electrical conductivity in the deep mantle. The responses in the period range that are shorter than 10 5 s (about 1 day) are in many cases considered to be local or regional induction problems in which the source field is approximated by plane waves and therefore the sphericity of Earth is not taken into account. In the period range between 10 4 and 10 5 s, the most dominant signature of the magnetic field variation is the solar quiet daily ( Sq ) variation and its higher harmonics. Therefore, when we obtain the responses due to the quasi-white background spectrum composed of plane waves, we regard the Sq field variations as noises in estimating the responses, and line spectra of the variations are removed from the observed time-series before the responses are calculated. However, with this approach, the calculated responses tend to possess a discontinuity at a period of about 10 4 s, and the response functions show common features at longer periods irrespective of the location of the observation site. It is particularly well known that the imaginary part of the induction vector tends to have a significant westward component for periods ranging between 10 4 and 10 5 s. Such features cannot be easily explained by the effect of the electrical conductivity structure alone. Examination of the phase of the HTF implies that the responses in the same period range are affected by the signature of sources of finite wavelength moving westward. Thus, it is suggested that the response functions in this period range were under the effect of the Sq field variations, even though the line spectra of them were removed and the responses were estimated at periods separate from the harmonic periods of Sq field variation. We examined how the influences of the Sq field appear on the responses numerically by a forward modelling. Results show most of the characteristic features in observed response functions can be ascribed to Sq source effects.

Description: SUMMARY Dependent on the ‘intrinsic’ effects on the crystal lattice of the rock constituents and the diminishing ‘extrinsic’ effects of pores and microcracks, elastic wave velocity versus pressure trends in cracked rocks are characterized by non-linear velocity increase at low pressure. At high pressure the ‘extrinsic’ influence vanishes and the velocity increase becomes approximately linear. Usually, the transition between non-linear and linear behaviour, the ‘crack closure pressure’, is not accessible in an experiment, because actual equipment is limited to lower pressure. For this reason, several model functions for describing velocity—pressure trends were proposed in the literature to extrapolate low-pressure P -wave velocity measurements to high pressures and, in part, to evaluate the ‘intrinsic’ velocity—pressure trend from low-pressure data. Knowing the ‘intrinsic’ velocity trend is of particular importance for the quantification of the crack influence at low pressure, at high pressure, the ‘intrinsic’ trend describes the velocity trend as a whole sufficiently well. Checking frequently used model functions for suitability led to the conclusion that all relations are unsuitable for the extrapolation and, if applicable, the estimation of the ‘intrinsic’ velocity trend. However, it can be shown that the ‘intrinsic’ parameters determined by means of a suitable model function, the zero pressure velocity and the pressure gradient depend on maximum experimental pressure in a non-linear way. Our approach intends to obtain better estimates of particular parameters from observed non-linear behaviour. A converging exponential function is used to approximate particular trends, assuming that the point of convergence of the function represents a better estimate of the zero pressure velocity and the pressure gradient, respectively. Whether the refined ‘intrinsic’ velocity trend meets the ‘true intrinsic’ velocity trend within acceptable errors cannot be proven directly due to missing experimental data at very high pressure. We, therefore, conclude that our approach cannot ensure absolutely certain ‘intrinsic’ velocity trends, however, it can be shown that the optimized trends approximate the ‘true intrinsic’ velocity trend better as all the other relations do.

Description: SUMMARY Several papers in recent years have reported instabilities in P -wave modelling, based on an acoustic approximation, for inhomogeneous transversely isotropic media with tilted symmetry axis (TTI media). In particular, instabilities tend to occur if the axis of symmetry varies rapidly in combination with strong contrasts of medium parameters, which is typically the case at the foot of a steeply dipping salt flank. In a recent paper, we have proposed and demonstrated a P -wave modelling approach for TTI media, based on rotated stress and strain tensors, in which the wave equations reduce to a coupled set of two second-order partial differential equations for two scalar stress components: a normal component along the variable axis of symmetry and a lateral component of stress in the plane perpendicular to that axis. Spatially constant density is assumed in this approach. A numerical discretization scheme was proposed which uses discrete second-derivative operators for the non-mixed second-order derivatives in the wave equations, and combined first-derivative operators for the mixed second-order derivatives. This paper provides a complete and rigorous stability analysis, assuming a uniformly sampled grid. Although the spatial discretization operator for the TTI acoustic wave equation is not self-adjoint, this operator still defines a complete basis of eigenfunctions of the solution space, provided that the solution space is somewhat restricted at locations where the medium is elliptically anisotropic. First, a stability analysis is given for a discretization scheme, which is purely based on first-derivative operators. It is shown that the coefficients of the central difference operators should satisfy certain conditions. In view of numerical artefacts, such a discretization scheme is not attractive, and the non-mixed second-order derivatives of the wave equation are discretized directly by second-derivative operators. It is shown that this modification preserves stability, provided that the central difference operators of the second-order derivatives dominate over the twice applied operators of the first-order derivatives. In practice, it turns out that this is almost the case. Stability of the desired discretization scheme is enforced by slightly weighting down the mixed second-order derivatives in the wave equation. This has a minor, practically negligible, effect on the kinematics of wave propagation. Finally, it is shown that non-reflecting boundary conditions, enforced by applying a taper at the boundaries of the grid, do not harm the stability of the discretization scheme.

Description: SUMMARY Broad-band seismograms of teleseismic events recorded at the Croatian Seismological Network were used to compute radial receiver functions (RFs) for eight locations in the External Dinarides. Waveform modelling was performed by a multistep matching of the theoretical RFs computed for horizontally layered 1-D isotropic models with the averaged observed RFs. Constraints from existing deep seismic sounding profiles, traveltime curves of regional crustal seismic phases and intuitive inferences gained from interactive forward modelling were used to construct initial 1-D models of the Earth. A non-linear inversion was performed in two steps—a grid search followed by the Monte Carlo search for the model parameters. Concurrently, RFs from different azimuths were stacked to obtain trade-off estimates of crustal thickness versus V p / V s ratios. The Moho depths were found in the range from around 40 km for Northern Adriatic stations to over 55 km for stations in the central part of the External Dinarides. Comparing our results with recent maps of the Moho topography inferred from seismic and gravimetric data, we find that for some stations the agreement between our results and the existing Moho maps is very good. For the others, we find the Mohorovičić discontinuity to be considerably deeper, indicating some of the thickest crust in Europe. Although it is plausible that such a deep Moho could be a consequence of a complex tectonic setting of the region (e.g. overlapping of two large tectonic units—the Adriatic microplate and the Dinarides), this result will have to be verified in the future studies using various other geophysical techniques.

Description: SUMMARY Gas hydrates are a potential energy resource and hazard for drilling and infrastructure, yet estimates of global volume vary by over three orders of magnitude. Hydrates are electrically resistive compared to water saturated sediment and so electromagnetic methods provide an additional tool to seismic surveys and drilling for determining hydrate saturations. A marine electromagnetic survey was carried out at Hydrate Ridge, Oregon, USA, with the aim of testing the use of controlled source electromagnetic (CSEM) and magnetotelluric (MT) methods to map gas hydrate and free gas below the gas hydrate stability zone. A 2-D CSEM inversion supports the scenario deduced from previous seismic and drilling results, which indicate two mechanisms of hydrate emplacement: a transport-dominated and reaction-dominated regime. A prominent resistive region of 2.5–4 Ωm at a depth of about 130 mbsf, near the seismic bottom simulating reflector (BSR), suggests that 27 to 46 per cent of the bulk volume is filled with hydrate, depending on whether Archie's Law or the Hashin-Strikman bounds are used. This is representative of a reaction-dominated regime for hydrate emplacement, and where a significant low velocity zone exists based on a seismic tomography inversion, suggests large quantities of free gas below the BSR. Electrical resistivity logging while drilling (LWD) data show general agreement with the CSEM inversion model except for a CSEM-derived resistive region at seismic horizon A, known to transport free gas into the gas hydrate stability zone. Inversion of MT data collected simultaneously during the CSEM survey provides a complimentary low-resolution image of the shallow sediments and shows folding in the accretionary complex sediments similar to that imaged by a tomographic seismic velocity model.

Description: SUMMARY We present a new 3-D vector finite element code and demonstrate its strength by modelling a realistic marine CSEM scenario. Unstructured tetrahedral meshes easily allow for the inclusion of arbitrary seafloor bathymetry so that natural environments are mapped into the model in a close-to-reality way. A primary/secondary field approach, an adaptive mesh refinement strategy as well as a higher order polynomial finite element approximation improve the solution accuracy. A convergence study strongly indicates that the use of higher order finite elements is beneficial even if the solution is not globally smooth. The marine CSEM scenario also shows that seafloor topography gives an important response which needs to be reproduced by numerical modelling to avoid the misinterpretation of measurements.

Description: SUMMARY The geographical distribution of the (1981–2005) seismicity in southern California forms a ±150 km broad zone adjacent to the Pacific–North America plate boundary, ranging from depths of ∼1–∼30 km, with the bulk of the focal depths in the range of 2–12 km. The distribution of the seismicity that includes both mainshock–aftershock sequences and background events is affected by both static and kinematic geophysical parameters of the crust. The static parameters include heat flow, topography, crustal density, V p / V s ratio, hypocentral fault-distance and crustal thickness from receiver functions. The tectonic loading is represented by kinematic parameters such as the crustal shear strain rate field, and the dilatational strain rate field. In our analysis, we normalize the seismicity relative to the areal density of the range of values of each of the parameters. Most of the seismicity occurs in areas of average heat flow, low to intermediate topography, average V p / V s and high late Quaternary fault density, and forms seismogenic zones that extend through the brittle crust. The location of late Quaternary faults, often described as zones of weakness, influences the geographical distribution of seismicity more than any other parameter. Although above or below average crustal properties such as high heat flow, thin crust or very low V p / V s values exist, these properties only influence the spatial distribution of seismicity in a minor way. As an example, the Salton Trough area of low topography, high heat flow, high V p / V s , high shear strain rate and thin crust has distributed seismicity within a thin seismogenic zone. Also, somewhat surprisingly, areas of high topography, low heat flow, low V p / V s , low shear strain rate and thick crust have low seismicity rates but a thin seismogenic zone. We determine an empirical relationship between heat flow and crustal thickness to show how the ∼400 °C temperature isotherm gradually deepens with crustal thickness and forms the base of the seismogenic zone for crustal thicknesses from 22 to 37 km. For crustal thickness ranging from 37 to 43 km, the ∼250 °C isotherm forms the base of the seismogenic zone, suggesting that seismic faulting in these regions is confined to the top of the upper crust (12–14 km), and thus does not accommodate plate motion.

Description: SUMMARY In this study the attenuation mechanism of seismic wave energy in north central Italy is estimated using low-magnitude earthquake local data recorded at six stations managed by INGV. Most of the analysed events are located along the Alpine chain in the zone of Iseo and Garda lakes, while a minor part in the Po valley. The zone investigated is characterized by the occurrence of significantly intense earthquakes (magnitude up to 6.6) the most recent occurred in 2004 close to the city of Salò on the coast of the Garda lake ( M w = 5.0). Due to the high population density and presence of industrial activity the investigated area is characterized by a high seismic risk. First, the ordinary Multiple Lapse Time Window Analysis (MLTWA) method is applied in the assumption of uniform velocity and scattering and the couple of B 0 , the seismic albedo and Le −1 , the extinction length inverse (corresponding to the total attenuation coefficient) is calculated in the frequency bands of 1.5, 3, 6 and 12 Hz. To retrieve more realistic estimates, the obtained values of B 0 and Le −1 are corrected taking into account the effects of a depth-dependent earth model, consisting of an earth structure characterized by a transparent upper mantle and a heterogeneous crust. We find that the corrected intrinsic and scattering attenuation parameters (which are proportional to the inverse of the intrinsic/scattering quality factors, Q −1 I and Q −1 S ) are strongly frequency dependent, with a prevalence of scattering attenuation over the intrinsic dissipation. The corrected and uncorrected values of total Q are in agreement with the total Q values obtained with different approaches for the same area.

Description: SUMMARY High-frequency and full waveform synthetic seismograms on a 3-D laterally heterogeneous earth model are simulated using the theory of coupled normal modes. The set of coupled integral equations that describe the 3-D response are simplified into a set of uncoupled integral equations by using the Born approximation to calculate scattered wavefields and the pure-path approximation to modulate the phase of incident and scattered wavefields. This depends upon a decomposition of the aspherical structure into smooth and rough components. The uncoupled integral equations are discretized and solved in the frequency domain, and time domain results are obtained by inverse Fourier transform. Examples show the utility of the normal mode approach to synthesize the seismic wavefields resulting from interaction with a combination of rough and smooth structural heterogeneities. This approach is applied to an ∼4 Hz shallow crustal wave propagation around the site of the San Andreas Fault Observatory at Depth (SAFOD).

Description: SUMMARY We analyse data from seismic stations surrounding the Alboran Sea between Spain and North Africa to constrain variations of the lithosphere–asthenosphere boundary (LAB) in the region. The technique used is the receiver function technique, which uses  S -to- P  converted teleseismic waves at the LAB below the seismic stations. We confirm previous data suggesting a shallow (60–90 km) LAB beneath the Iberian Peninsula and we observe a similarly shallow LAB beneath the Alboran Sea where the lithosphere becomes progressively thinner towards the east. A deeper LAB (90–100 km) is observed beneath the Betics, the south of Portugal and Morocco. The structure of the LAB in the entire region does not seem to show any indication of subduction related features. We also observe good  P  receiver function signals from the seismic discontinuities at 410 and 660 km depth which do not indicate any upper-mantle anomaly beneath the entire region. This is in agreement with the sparse seismic activity in the mantle transition zone suggesting the presence of only weak and regionally confined anomalies.

Description: SUMMARY We perform numerical experiments to investigate the accuracy of broad-band moment magnitude, M wp . We conduct these experiments by measuring M wp from synthetic seismograms and comparing the resulting values to the moment magnitudes used in the calculation of synthetic seismograms. In the numerical experiments using point sources, we have found that there is a significant dependence of M wp on focal mechanisms, and that depths phases have a large impact on M wp estimates, especially for large shallow earthquakes. Numerical experiments using line sources suggest that the effects of source finiteness and rupture propagation on M wp estimates are on the order of 0.2 magnitude units for vertical fault planes with pure dip-slip mechanisms and 45° dipping fault planes with pure dip-slip (thrust) mechanisms, but that the dependence is small for strike-slip events on a vertical fault plane. Numerical experiments for huge thrust faulting earthquakes on a fault plane with a shallow dip angle suggest that the M wp estimates do not saturate in the moment magnitude range between 8 and 9, although they are underestimates. Our results are consistent with previous studies that compared M wp estimates to moment magnitudes calculated from seismic moment tensors obtained by analyses of observed data.

Description: SUMMARY In areas dominated by Glacial Isostatic Adjustment (GIA), the free-air gravity anomaly rate can be converted to uplift rate to good approximation by using a simple spectral relation. We provide quantitative comparisons between gravity rates derived from monthly gravity field solutions (GFZ Potsdam, CSR Texas, IGG Bonn) from the Gravity Recovery and Climate Experiment (GRACE) satellite mission with uplift rates measured by GPS in these areas. The band-limited gravity data from the GRACE satellite mission can be brought to very good agreement with the point data from GPS by using scaling factors derived from a GIA model (the root-mean-square of differences is 0.55 mm yr –1 for a maximum uplift rate signal of 10 mm yr –1 ). The root-mean-square of the differences between GRACE derived uplift rates and GPS derived uplift rates decreases with increasing GRACE time period to a level below the uncertainty that is expected from GRACE observations, GPS measurements and the conversion from gravity rate to uplift rate. With the current length of time-series (more than 8 yr) applying filters and a hydrology correction to the GRACE data does not reduce the root-mean-square of differences significantly. The smallest root-mean-square was obtained with the GFZ solution in Fennoscandia and with the CSR solution in North America. With radial gravity rates in excellent agreement with GPS uplift rates, more information on the GIA process can be extracted from GRACE gravity field solutions in the form of tangential gravity rates, which are equivalent to a rate of change in the deflection of the vertical scaled by the magnitude of gravity rate vector. Tangential gravity rates derived from GRACE point towards the centre of the previously glaciated area, and are largest in a location close to the centre of the former ice sheet. Forward modelling showed that present day tangential gravity rates have maximum sensitivity between the centre and edge of the former ice sheet, while radial gravity rates are most sensitive in the centre of the former ice sheet. As a result, tangential gravity rates offer constraints on a two-layer mantle viscosity profile that are different from radial gravity rates, which can be exploited in future GIA studies.

Description: SUMMARY The study of electromagnetic (EM) field sensitivities is useful for assessing the feasibility of controlled-source electromagnetic (CSEM) surveys. Sensitivity calculations are also a principal building block of EM inversion schemes. Sensitivities are formally given by the derivatives of the EM field components with respect to conductivity. For horizontally layered media, these derivatives can be evaluated analytically, offering advantages in computational efficiency and accuracy over numerical evaluation. We present a complete set of explicit analytic expressions for the EM field sensitivities in 1-D VTI-anisotropic media for horizontal and vertical electric and magnetic dipole sources, and also for finite horizontal electric sources. Since our derivations are based on a formulation for EM fields that is quite general in allowing for sources and receivers at any depth, our sensitivity expressions exhibit the same generality. We verify our expressions by comparison to numerical solutions, and finally present application examples that demonstrate the utility and versatility of these expressions for CSEM feasibility studies.

Description: SUMMARY 3-D numerical models describing convection cells that may develop in a buoyant layer over an upper-mantle hotspot are developed. Results are compared with those of an alternative model, which assumes an upper partially molten layer having the same density as the underlying mantle, but of much higher viscosity. Geophysical and geochemical data over Iceland are used to assess the models. For the case of a buoyant upper layer, a major finding of our modelling is the prediction of a cold downward moving sheet in the upper layer situated just above the axis of the hot underlying mantle hotspot. When compared to geophysical data from Iceland, this modelled feature may explain why: (i) along the hotspot axis, seismic tomography reveals a mantle that remains cold down to ∼80 km depth; (ii) the ridge crest systematically avoids the axis of the hotspot and (iii) the ridge axis is displaced by ∼60 km to the east or the west of the hotspot axis, these ridge axes corresponding to the modelled upwellings of hot material in the uppermost layer driven by the cold downwelling axial sheet. Finally, the model also predicts development of gravity swells a few hundred kilometres long and several tens of kilometres wide, with a peak to valley amplitude of 5–10 mGal, lying over ascending sheets of upwelling material that float over the hotspot. Comparison with the residual gravity map over Iceland does indeed show the presence of a long gravity high along the active eastern section of the ridge. Over the AskjaVolcano, the long gravity high intersects an off-axis gravity high oriented northwest, which joins the Kolbeinsey Ridge. It is proposed that this branching may be the cause of the eastward drift of the western section of the Iceland ridge that took place 3 Ma. Further south, parallel to the same ridge-associated positive gravity anomaly, a more recent branching is found with a very long gravity high, oriented in an easterly direction, but that eventually contours the whole eastern coast of Iceland. This latter feature may be the precursor of future eastward drift of the whole eastern section of the ridge. When combined with experimentally derived models of dry peridotite melting, it is found that the temperature field of the models is consistent with the principal constraints derived from compositional characteristics of Icelandic basalts and the 3-D map of the upper mantle below the Icelandic hotspot derived from S -wave seismic tomography. The high viscosity layer model is also relatively successful at explaining geophysical and petrological data, except that it does not predict the development of hot and cold sheets in the partially molten mantle, and it leads to a large and strong negative gravity and bathymetry anomaly over the Iceland Hotspot, at a place where available data show a swell. However, this latter result is found to depend on the lateral boundary conditions used in the numerical modelling. For this reason we conclude that the presence of a very high viscosity layer below the ridge cannot be ruled out on the basis of gravity and bathymetric data alone. On the other hand, the influence of a floating upper layer provides an interesting and plausible alternative to explain available geophysical data, given that experimental data tend to indicate only modest changes in viscosity of the mantle in response to partial melting.

Description: SUMMARY We use a multiscale approach as a semi-automated interpreting tool of potential fields. The depth to the source and the structural index are estimated in two steps: first the depth to the source, as the intersection of the field ridges (lines built joining the extrema of the field at various altitudes) and secondly, the structural index by the scale function. We introduce a new criterion, called ‘ridge consistency’ in this strategy. The criterion is based on the principle that the structural index estimations on all the ridges converging towards the same source should be consistent. If these estimates are significantly different, field differentiation is used to lessen the interference effects from nearby sources or regional fields, to obtain a consistent set of estimates. In our multiscale framework, vertical differentiation is naturally joint to the low-pass filtering properties of the upward continuation, so is a stable process. Before applying our criterion, we studied carefully the errors on upward continuation caused by the finite size of the survey area. To this end, we analysed the complex magnetic synthetic case, known as Bishop model, and evaluated the best extrapolation algorithm and the optimal width of the area extension, needed to obtain accurate upward continuation. Afterwards, we applied the method to the depth estimation of the whole Bishop basement bathymetry. The result is a good reconstruction of the complex basement and of the shape properties of the source at the estimated points.

Description: SUMMARY The estimation of shear wave velocity and attenuation in near-surface geology is of primary importance in engineering seismology. In fact, such knowledge is essential for site response studies when preparing improved seismic hazard scenarios. In this study, we propose a linear inversion of the spectra of a deconvolved wavefield collected by a borehole array in Istanbul, Turkey. The spectra are calculated using as a reference the recordings collected by a sensor at the surface. This allows us to minimize the effect of the deconvolution filter on the peaks’ amplitudes. The feasibility of the proposed inversion scheme and the parametrization of the velocity and attenuation models were tested and assessed using synthetic data. The real data inversion is carried out using observation from three weak motion events. Several starting models are used for each event, changing the values of the quality factor  Qs , whereas the starting  Vs  profiles are fixed to what was obtained by the analysis of the deconvolved wavefield in the time domain. Our results showed that fairly well-constrained  Vs  models can be obtained while the variability of the  Qs  profile, both between different inversions of the same data set, and between the inversions of different data sets, might be large. This indicates that  Qs  variations observed in models derived by strong motion data inversions should be considered with care. However, the high agreement of the estimated soil profiles in the uppermost layer, both when using different events and different starting models, suggests that any variation in the  S -wave velocity that might be estimated when using strong motion data should be reliable.

Description: SUMMARY The rapid melting of the Earth’s ice reservoirs will produce geographically distinct patterns of sea level change that have come to be known as sea level fingerprints. A basic, gravitationally self-consistent theory for computing these patterns appeared in the 1970s; however, recent, highly discrepant fingerprint calculations have led to suggestions that the algorithms and/or theoretical implementation adopted in many previous predictions is not robust. We present a suite of numerical predictions, including benchmark comparisons with analytic results, that counter this argument and demonstrate the accuracy of most published predictions. Moreover, we show that small differences apparent in calculations published by some groups can be accounted for by subtle differences in the underlying physics. The paper concludes with two sensitivity analyses: (1) we present the first-ever calculation of sea level fingerprints on earth models with 3-D variations in elastic structure and density, and conclude that this added complexity has a negligible effect on the predictions; (2) we compare fingerprints of polar ice sheet mass flux computed under the (very common) assumption of a uniform melt distribution to fingerprints calculated using melt geometries constrained by analysing recent trends in GRACE gravity data. Predictions in the near field of the ice sheets are sensitive to the assumed melt geometry; however, this sensitivity also extends to the far field, particularly in the case of Antarctic mass changes, because of the strong dependence of the rotational feedback signal on the melt geometry. We conclude that inferences of ice sheet mass flux based on modern sea level constraints should consider these more realistic melt geometries.

Description: SUMMARY Inverting for the subsurface velocity distribution by refraction traveltime tomography is a well-accepted imaging method by both the exploration and earthquake seismology communities. A significant drawback, however, is that the recorded traces become noisier with increasing offset from the source position, and so accurate picking of traveltimes in far-offset traces is often prevented. To enhance the signal-to-noise ratio (SNR) of the far-offset traces, we present the theory of supervirtual refraction interferometry where the SNR of far-offset head-wave arrivals can be theoretically increased by a factor proportional to ; here, N is the number of receiver or source positions associated with the recording and generation of the head-wave arrival. There are two steps to this methodology: correlation and summation of the data to generate traces with virtual head-wave arrivals, followed by the convolution of the data with the virtual traces to create traces with supervirtual head-wave arrivals. This method is valid for any medium that generates head-wave arrivals recorded by the geophones. Results with both synthetic traces and field data demonstrate the feasibility of this method. There are at least four significant benefits of supervirtual interferometry: (1) an enhanced SNR of far-offset traces so the first-arrival traveltimes of the noisy far-offset traces can be more reliably picked to extend the useful aperture of the data, (2) the SNR of head waves in a trace that arrive later than the first arrival can be enhanced for accurate traveltime picking and subsequent inversion by later-arrival traveltime tomography, (3) common receiver-pair gathers can be analysed to detect the presence of diving waves in the first arrivals, which can be used to assess the nature of the refracting boundary, and (4) the source statics term is eliminated in the correlation operations so that the timing of the virtual traces is independent of the source excitation time. This suggests the possibility of applying this method to earthquake data recorded by receivers that are inline with the refraction paths and source locations.

Description: SUMMARY We present an analytic boundary layer model for thermal convection with a finite-strength plate and depth-dependent viscosity. The model permits solutions in which convective flow rates in the mantle exceed the plate velocity. The energy balance equations for the lithosphere and convective cell are distinct and the model yields the plate velocity, plate thickness and heat flow, as well as the laterally averaged horizontal flow profile for a convective cell. We demonstrate, by inspection of the lithospheric energy balance, that the dominant plate driving force depends not only on the material properties of the mantle and lithosphere, but also on the velocity of the plate. Multiple solutions are possible with three solution branches for the plate velocity representing three distinct modes of thermal convection. The branch of solutions with the largest plate velocity corresponds to the classic boundary layer solution, with a plate velocity approximately equal to the maximum mantle velocity. This branch reproduces the classic convective scaling laws for an isoviscous fluid and its dynamics are controlled by the mantle material properties. The branch of solutions with intermediate plate velocity represents a convective cell with a sluggish-lid and a plate velocity that is less than the underlying mantle velocity. The dynamics for this solution branch depend on the material properties of both the mantle and the lithosphere. Finally, the lower branch of solutions for the plate velocity yields a convective solution in which the plate velocity is determined entirely by the local dynamics of a thick and strong (relative to the mantle) lithosphere. The dynamics of this solution branch are independent of the material properties of the mantle and depend entirely on the properties of the lithosphere. The introduction of an asthenosphere or low viscosity layer (LVL) beneath the plate can significantly alter the dynamics of the system by affecting plate–mantle coupling. Lowering the LVL viscosity, relative to the lower mantle, can promote plate motion by providing a lubricating layer. However, a very low viscosity LVL beneath a strong plate can decouple plate and mantle and inhibit plate motion, producing a solution with a slow moving plate and a channelized flow in the LVL. Thus, a LVL can sometimes inhibit rather than promote plate tectonics.

Description: SUMMARY We developed an integrated technique for estimating crustal anisotropy with a horizontal axis using radial (R) and transverse (T) receiver functions. The technique includes computing three individual and one joint objective function (JOF) and a reliability analysis of the estimated anisotropy. The individual objective functions (IOFs) are designed to: (1) maximize the peak energy of the stacked R receiver function after a cosine moveout correction in the  Ps  arrival time; (2) to maximize the correlation of the radial receiver functions after a full correction of anisotropy or (3) to minimize the total energy of transverse receiver functions stacked after a removal of crustal anisotropy. The JOF was computed by a weighted average of the three IOFs, and the reliability analysis uses the principle that stacking coherent signals can lead to an increase of signal-to-noise ratio. We applied the technique to synthetic receiver functions generated with 30–60 per cent white noise from a variety of anisotropic and heterogeneous models. The synthetic tests indicate that the proposed technique has good capability to recover the input models. Despite the presence of random and other coherent noises, such as those caused by inhomogeneous structures, in the data, the technique can always provide accurate estimates of crustal anisotropy. We applied the technique to two permanent seismic stations in western China and found significant crustal anisotropy beneath one station located at the northern edge of the Tibetan plateau. The observed fast polarization direction at this station follows the direction of the maximum horizontal tensile stress, suggesting that the observed seismic anisotropy is likely caused by mineral alignment in the lower crust. The station situated in the Sichuan basin, on the other hand, shows little to no seismic anisotropy, which may suggest that the crust beneath the basin is nearly rigid with very little deformation. The developed technique can be applied to any broadband seismic stations that have a good backazimuthal coverage of teleseismic events.

Description: SUMMARY The Markov approximation is a powerful stochastic method for the direct synthesis of wave envelopes in random velocity fluctuated media when the wavelength is shorter than their correlation distance. To apply the Markov approximation to realistic cases, we consider horizontal layered random media characterized by an anisotropic autocorrelation function (ACF), where different layers have different randomnesses and different background velocities having step-like changes. Solving the parabolic master equation for the two frequency mutual coherence function (TFMCF) in random elastic media for the vertical incidence of a plane wavelet from the bottom, we calculate the angular spectrum just before the first velocity boundary. Multiplying transmission or PS conversion coefficients of the boundary by the angular spectrum, we calculate the angular spectrum and TFMCF on the other side of the boundary. Then we solve the master equation for the forward propagating wavelet in the second layer. Taking the same procedure for each layer boundary, we finally obtain the mean square (MS) vector envelopes on the free surface on the top layer. For the practical simulation, we use 2-D random media characterized by a Gaussian ACF. We numerically confirm the validity of the envelope synthesis for a specific case of layered random media with anisotropic ACF by comparing with finite difference (FD) simulations of elastic waves. Considering the Earth structure, the horizontal correlation distance is larger than the vertical one and the velocity fractional fluctuation becomes weak as depth increases, the Markov approximation is good for modelling the primary wavelet and also the converted wavelet for the vertical incident wavelet. We derive an applicable range of the Markov approximation for random media with anisotropic ACF by comparing with FD simulations. The results show that the Markov approximation is accurate when the wavelength is comparable or shorter than the both of vertical and horizontal correlation distances and the MS fractional fluctuation is much smaller than the ratio of squared horizontal correlation distance to the product of the vertical correlation distance and propagation distance.

Description: SUMMARY In finite-difference (FD) acoustic forward modelling, the parameter settings for the perfectly matched layer (PML) are case-dependent. There is no explicit PML formula that can be applied for most acoustic models without tuning, especially, for the fourth-order FD scheme. In this paper, we propose an explicit PML formula for the acoustic frequency-domain FD with second-order and fourth-order accuracies, respectively. The fourth-order FD scheme uses a special treatment for the boundary. The number of points in the PML is fixed to be 10 and 15 for the second-order and fourth-order FD schemes, respectively. The maximum artificial attenuation parameter associated with the PML formula is automatically calculated based on the FD grid size and the value of the compressional velocity of the boundary cells of the interior domain. Numerical tests confirm that this empirical formula achieves the desired accuracy for 2-D and 3-D media for grid sizes varying from 1 to 200 m. For the fourth-order FD scheme, the proposed PML formula works effectively up to 25 points per wavelength for both 2-D and 3-D media. Beyond that, the error of the PML discretization becomes larger than the discretization error in the interior domain. For such cases and to keep a fourth-order accuracy, a larger number of points in the PML (thicker PML region) needs to be employed.

Description: SUMMARY A large number of fundamental-mode Love and Rayleigh wave dispersion curves were determined from seismograms for 3330 earthquakes recorded on 258 globally distributed seismographic stations. The dispersion curves were sampled at periods between 25 and 250 s to determine propagation-phase anomalies with respect to a reference earth model. The data set of phase anomalies was first used to construct global isotropic phase-velocity maps at specific frequencies using spherical-spline basis functions with a nominal uniform resolution of 650 km. Azimuthal anisotropy was then included in the parametrization, and its importance for explaining the data explored. Only the addition of 2ζ azimuthal variations for Rayleigh waves was found to be resolved by the data. In the final stage of the analysis, the entire phase-anomaly data set was inverted to determine a global dispersion model for Love and Rayleigh waves parametrized horizontally using a spherical-spline basis, and with a set of B-splines to describe the slowness variations with respect to frequency. The new dispersion model, GDM52, can be used to calculate internally consistent global maps of phase and group velocity, as well as local and path-specific dispersion curves, between 25 and 250 s.

Description: SUMMARY The separation of Northeast Greenland and Svalbard is the result of large-scale strike slip movements during Cenozoic times. Geological evidence for these movements can be found onshore both on North Greenland and Svalbard. However, the role of the submarine Yermak Plateau (YP) in this process is unclear. The compilation of available multichannel reflection and wide-angle seismic data give new insights into the sedimentary and crustal structure and evolution of the plateau. The flat surface of the present-day plateau is a quite young feature. Up to 2 km of Cenozoic sediments cover a rough basement, which show similarities to the rough topography and strike of geological structures of Spitsbergen Island. In some basins more than 4 km of sedimentary rocks could be mapped. The most pronounced structure is the Sverdrup Bank, which appears to be part of a larger crustal block. P -wave velocities of about 4.5 km s –1 derived from sonobuoy data indicate that its uppermost part is most probably composed of sedimentary or volcanic rocks. We have made a correlation of previously defined seismic units across the YP to outline the history of sediment deposition in the area. The existing graben structures on the plateau might have provided early shallow pathways for water exchange between the Arctic and the Atlantic Oceans. A chaotic sedimentary apron east of the Sverdrup Bank and bright reflections near the Mosby Seamount interpreted as magmatic sills suggest tectonic and magmatic events during the Miocene.

Description: Greenstreet, S. P. R., Fraser, H. M., Rogers, S. I., Trenkel, V. M., Simpson, S. D., and Pinnegar, J. K. 2012. Redundancy in metrics describing the composition, structure, and functioning of the North Sea demersal fish community. – ICES Journal of Marine Science, 69: 8–22. Broader ecosystem management objectives for North Sea demersal fish currently focus on restoring community size structure. However, most policy drivers explicitly concentrate on restoring and conserving biodiversity, and it has not yet been established that simply restoring demersal fish size composition will be sufficient to reverse declines in biodiversity and ensure a generally healthy community. If different aspects of community composition, structure, and function vary independently, then to monitor all aspects of community general health will require application of a suite of metrics. This assumes low redundancy among the metrics used in any such suite and implies that addressing biodiversity issues specifically will require explicit management objectives for particular biodiversity metrics. This issue of metric redundancy is addressed, and 15 metrics covering five main attributes of community composition, structure, and function are applied to groundfish survey data. Factor analysis suggested a new interpretation of the metric information and indicated that a minimum suite of seven metrics was necessary to ensure that all changes in the general health of the North Sea demersal fish community were monitored properly. Covariance among size-based and species-diversity metrics was low, implying that restoration of community size structure would not necessarily reverse declines in species diversity.

Description: Dankel, D. J., Aps, R., Padda, G., Röckmann, C., van der Sluijs, J. P., Wilson, D. C., and Degnbol, P. 2012. Advice under uncertainty in the marine system. – ICES Journal of Marine Science, 69: 3–7. There is some uncertainty in the fisheries science–policy interface. Although progress has been made towards more transparency and participation in fisheries science in ICES Areas, routine use of state-of-the-art quantitative and qualitative tools to address uncertainty systematically is still lacking. Fisheries science that gives advice to policy-making is plagued by uncertainties; the stakes of the policies are high and value-laden and need therefore to be treated as an example of "post-normal science" (PNS). To achieve robust governance, understanding of the characteristics and implications of the scientific uncertainties for management strategies need to come to the centre of the table. This can be achieved using state-of-the-art tools such as pedigree matrices and uncertainty matrices, as developed by PNS scholars and used in similar science–policy arenas on other complex issues. An explicit extension of the peer community within maritime systems will be required to put these new tools in place. These new competences become even more important as many countries within the ICES Area are now embarking on new policies.

Description: Browman, H.I. 2012. Quo Vadimus. – ICES Journal of Marine Science, 69: 1–2.

Description: Hüssy, K., Coad, J. O., Farrell, E. D., Clausen, L. A. W., and Clarke, M. W. 2012. Age verification of boarfish ( Capros aper ) in the Northeast Atlantic. – ICES Journal of Marine Science, 69: 34–40. The boarfish ( Capros aper ) is a pelagic species of recent interest to the fishing industry, with landings increasing by 〉500% over the past 3 years. The objective of the study was to provide a method for age determination based on whole sagittal otoliths, with the results to be used in stock assessment. Translucent zones laid down at regular intervals are identified by marginal increment analysis as seasonally recurring. Translucent zones are formed between September/October and March/April, regardless of fish age. The occurrence of the first annulus is validated by analysis of presumed daily growth increments. Subsequent annulus deposition is homogenous between individuals and allows general guidelines to be derived for interpreting the age of boarfish using their otoliths.

Description: Kim, D., Yang, E. J., Kim, K. H., Shin, C-W., Park, J., Yoo, S., and Hyun, J-H. 2012. Impact of an anticyclonic eddy on the summer nutrient and chlorophyll a distributions in the Ulleung Basin, East Sea (Japan Sea). – ICES Journal of Marine Science, 69: 23–29. The impact of the anticyclonic Ulleung Warm Eddy (UWE) on the vertical distributions of nutrient and chlorophyll a (Chl a ) concentrations in the Ulleung Basin (UB) was investigated during the contrasting summers of 2005 and 2007. The physical structure of the water column was characterized by an intrathermocline eddy (ITE) in 2005, whereas the UWE remained distant from the sampling transect in 2007. Water column structures appeared to be highly stratified, and nutrients in the surface waters were totally depleted at all stations. In 2005, an exceptionally high concentration of Chl a (5.5 mg m –3 ) was measured below the surface mixed layer in the eddy core (station D3), and values of ~2.5 mg m –3 were observed at the eddy edge (stations D2 and D4). Formation of an ITE efficiently mixed surface and deep-ocean waters, the latter supplying sufficient nutrients to generate an extremely high concentration of Chl a at the base of the subsurface layer. Overall, the results indicated that the anticyclonic UWE plays a key ecological role in supporting substantial phytoplankton biomass in the nutrient-depleted surface waters in summer and maintaining high benthic mineralization in the deep-sea sediments of the UB.

Description: Franchini, P., Sola, L., Crosetti, D., Milana, V., and Rossi, A. R. 2012. Low levels of population genetic structure in the gilthead sea bream, Sparus aurata , along the coast of Italy. – ICES Journal of Marine Science, 69: 41–50. The gilthead sea bream, Sparus aurata , is a coastal, commercially important fish. Contrasting results concerning the genetic structure of the species at different geographic scales have been reported. Here, an investigation is made into the population genetic structure of S. aurata along the coast of Italy, using samples analysed previously and material from new sampling sites (12) and using different microsatellite loci (10). One sample from the eastern Atlantic and three temporal replicates from one site were also included. The presence of a weak (overall F ST = 0.0072), but significant, genetic population subdivision was detected by F -statistics. Temporal replicates indicate genetic data consistency over time. Isolation by distance between the Atlantic and the coast of Italy is suggested by a Mantel test. The distributional pattern of genetic variance obtained by analysis of molecular variation reflects the geographic sampling areas, but is only partially congruent with the results obtained with fewer sites and loci. The dispersal of passive eggs/larvae by the main currents appears to contribute to shaping the gene flow. Given the intensity of sea bream aquaculture activities in Italy, the possibility that aquaculture may have partially contributed to the population genetic pattern detected cannot be excluded.

Description: Charrier, F., Mazel, V., Caraguel, J-M., Abdallah, Y., Le Gurun, L. L., Legault, A., and Laffaille, P. 2012. Escapement of silver-phase European eels, Anguilla anguilla , determined from fishing activities in a Mediterranean lagoon (Or, France). – ICES Journal of Marine Science, 69: 30–33. Escapement of silver eels from a Mediterranean lagoon was estimated by a capture–tagging–recapture and automated tag-reading study. The population of silver-phase eels in the lagoon was estimated to be 13.2 kg ha –1 , with an escapement rate from the commercial fishery of 76.8%.

Description: Jennings, S., and Lee, J. 2012. Defining fishing grounds with vessel monitoring system data. – ICES Journal of Marine Science, 69: 51–63. Methods for defining fishing grounds to support marine spatial planning and management are developed, applied, and compared. The methods are broadly applicable and repeatable because they use vessel monitoring system (VMS) data that are archived and increasingly accessible. For several fleets at regional and national scales, an attempt is made to assess how the choice of criteria for defining grounds influences (i) size, shape, and location, (ii) overlap among grounds, and (iii) the extent to which annual and multi-annual patterns of fishing activity describe grounds used seasonally or by individual vessels. The results show that grounds defined by excluding infrequently fished margins (areas with 〈10% of total fishing activity) are typically 50% smaller than total fished area. However, landings weight or value (LWV) per unit activity can be higher at the margins, with 10% of activity usually accounting for 10–20% of LWV. The removal of fishing activity in the margins, as a consequence of regulation or fleet behaviour, would lead to disproportionately greater reductions in interactions with other fisheries, sectors, and the environment. Accessible high-resolution information on the "anatomy" of all fishing grounds would better inform debates on the allocation and the use of marine space and the integration of fisheries and environmental management.

Description: Bentley, N., Kendrick, T. H., Starr, P. J., and Breen, P. A. 2012. Influence plots and metrics: tools for better understanding fisheries catch-per-unit-effort standardizations. – ICES Journal of Marine Science, 69: 84–88. Standardization of catch per unit effort using generalized linear models (GLMs) is a common procedure that attempts to remove the confounding effects of variables other than abundance. Simple plots and metrics are described to assist understanding the standardization effects of explanatory variables included in GLMs, illustrated with an example based on New Zealand trevally ( Caranx lutescens ) data.

Description: Aldrin, M., Mortensen, B., Storvik, G., Nedreaas, K., Aglen, A., and Aanes, S. 2012. Improving management decisions by predicting fish bycatch in the Barents Sea shrimp fishery. – ICES Journal of Marine Science, 69: 64–74. When the bycatch of juvenile fish within the Barents Sea shrimp fishery is too large, the area is closed to fishing for a certain period. Bycatch is estimated from sampled trawl hauls, for which the shrimp yield is recorded, along with the total number of various bycatch fish species. At present, bycatch estimation is based on a simple estimator, the sum of the number of fish caught within the area of interest within a small time window, divided by the corresponding shrimp yield (in weight). No historical data are used. A model-based estimation is proposed in which spatio-temporal models are constructed for the variation in both the yield of shrimp and the amount of bycatch in space and time. The main effects are described through generalized additive models, and local dependence structures are specified through correlated random effects. Model estimation includes historical and recent data. Experiments with both simulated and real data show that the model-based estimator outperforms the present simple estimator when a low or moderate number of samples (e.g. 〈20) is available, whereas the two estimators are equally good when the number of samples is high.

Description: Jiao, Y., Smith, E. P., O'Reilly, R., and Orth, D. J. 2012. Modelling non-stationary natural mortality in catch-at-age models. – ICES Journal of Marine Science, 69: 105–118. Non-stationarity in the natural mortality of many species has been discussed among research scientists, but no generally applicable models/methods have been developed where a statistical catch-at-age (SCA) model framework is used. Using the Atlantic weakfish ( Cynoscion regalis ) fishery as an example, several SCA models are developed to assess the population dynamics, then compared. Models used included (i) an SCA with constant natural mortality, (ii) an SCA with unknown natural mortality, but with a hierarchical prior distribution from a mixture of distributions based on life-history information, (iii) an SCA with age-specific unknown natural mortality, (iv) an SCA with time-varying natural mortality, following a random-walk process, and (v) an SCA with age-specific time-varying natural mortality. The last two models imply that the population dynamics are not stationary. A Bayesian approach was used to estimate parameters, and performance of the models was compared through goodness-of-fit and the retrospective patterns of the models. A simulation study was then conducted to test the uncertainty resulting from model selection, the efficiency of using the best model selected based on deviance information criterion, and whether changes in natural mortality could be detected. An SCA with time-varying natural mortality, following a random-walk process, is proposed for the example fishery here. The estimated non-stationary temporal patterns in natural mortality were linked further to climate–ocean oscillation indices, to diagnose possible mechanisms/linkages to these temporal variations in natural mortality.

Description: Juntunen, T., Vanhatalo, J., Peltonen, H., and Mäntyniemi, S. 2012. Bayesian spatial multispecies modelling to assess pelagic fish stocks from acoustic- and trawl-survey data. – ICES Journal of Marine Science, 69: 95–104. A Bayesian spatial model was constructed to estimate the abundance of multiple fish species in a pelagic environment. Acoustic- and trawl-survey data were combined with environmental data to predict the spatial distribution of (i) the acoustic backscattering of fish, (ii) the relative proportion of each species, and (iii) their mean length in the Gulf of Finland in the northeastern Baltic Sea. By combining the three spatial model layers, the spatial distribution of the biomass of each species was estimated. The model consists of a linear predictor on environmental variables and a spatial random effect given by a Gaussian process. A Bayesian approach is a natural choice for the task because it provides a theoretically justified means of summarizing the uncertainties from various model layers. In the study area, three species dominate pelagic waters: sprat ( Sprattus sprattus ), herring ( Clupea harengus ), and three-spined stickleback ( Gasterosteus aculeatus ). Results are presented for each model layer and for estimated total biomass for each species in 2 x 2 km lattices. The posterior mean and central 95% credible intervals of total biomass were sprat 45.7 kt (27.7–71.6), herring 24.6 kt (9.7–41.3), and three-spined stickleback 1.9 kt (0.9–3.2).

Description: Smith, W. K., and Solow, A. R. 2012. Missing and presumed lost: extinction in the ocean and its inference. – ICES Journal of Marine Science, 69: 89–94. The number of modern extinctions in the ocean is unknown. The actual demise of the last individual of a species is essentially unobservable, so extinction can only be inferred. Statistical methods are described for inferring extinction from sighting records, species–area considerations, and taxonomic samples collected at two different times. The methods are illustrated using a variety of real datasets, including a sighting record of the Caribbean monk seal and results from three surveys of benthic invertebrates.

Description: Delavenne, J., Metcalfe, K., Smith, R. J., Vaz, S., Martin, C. S., Dupuis, L., Coppin, F., and Carpentier, A. 2012. Systematic conservation planning in the eastern English Channel: comparing the Marxan and Zonation decision-support tools. – ICES Journal of Marine Science, 69: 75–83. The systematic conservation approach is now commonly used for the design of efficient marine protected area (MPA) networks, and identifying these priority areas often involves using specific conservation-planning software. Several such software programmes have been developed in recent years, each differing in the underlying algorithms used. Here, an investigation is made into whether the choice of software influences the location of priority areas by comparing outputs from Marxan and Zonation, two widely used conservation-planning, decision-support tools. Using biological and socio-economic data from the eastern English Channel, outputs are compared and it is shown that the two software packages identified similar sets of priority areas, although the relatively wide distribution of habitat types and species considered offered much flexibility. Moreover, the similarity increased with increasing spatial constraint, especially when using real-world cost data, suggesting that the choice of cost metric has a greater influence on conservation-planning analyses than the choice of software. However, Marxan generally produced more efficient results and Zonation produced results with greater connectivity, so the most appropriate software package will depend on the overall goals of the MPA planning process.

Description: Andersen, B. S., Ulrich, C., Eigaard, O. R., and Christensen, A-S. 2012. Short-term choice behaviour in a mixed fishery: investigating métier selection in the Danish gillnet fishery. – ICES Journal of Marine Science, 69: 131–143. The study presents a short-term effort allocation modelling approach based on a discrete choice random utility model combined with a survey questionnaire to examine the selection of métiers (a combination of fishing area and target species) in the Danish North Sea gillnet fishery. Key decision variables were identified from the survey questionnaire, and relevant proxies for the decision function were identified based on available landings and effort information. Additional variables from the survey questionnaire were further used to validate and verify the outcome of the choice model. Commercial fishers in a mixed fishery make use of a number of decision variables used previously in the literature, but also a number of decision parameters rarely explicitly accounted for, such as price, weather, and management regulation. The seasonal availability of individual target species and within-year changes in monthly catch ration were the main explanatory drivers, but gillnetters were also responsive to information on the whole fishery, fish prices, and distance travelled to fishing grounds. Heterogeneous responses were evident from geographic differences in home harbour, which underpins the need to understand alternative fishing strategies among individual gillnetters better.

Description: Calise, L., and Knutsen, T. 2012. Multifrequency target strength of northern krill ( Meganyctiphanes norvegica ) swimming horizontally. – ICES Journal of Marine Science, 69: 119–130. Multifrequency acoustic measurements on ex situ horizontally swimming krill were made in a novel experimental setting. An ensemble of northern krill ( Meganyctiphanes norvegica ) was introduced to a large enclosure (a mesocosm), and acoustic backscatter was sampled using a multifrequency (70, 120, and 200 kHz) echosounder (Simrad EK60). Two submerged lamps were placed at opposite sides of the mesocosm and switched on and off to induce the krill, by light attraction, to swim horizontally through the acoustic beams. By tracking echoes, animal displacement, swimming speed, and target strength ( TS ) by frequency were estimated. The dominant and secondary modes of the total-length distribution were 21.8 ± 3.0 and 27.8 ± 2.7 mm, respectively. Although krill orientation was assumed stable and the ping rate was high, the range and inter-ping variability of the average TS values were large, decreasing and increasing with frequency, respectively. The overall TS frequency response observed and concurrent measurements at 120 and 200 kHz confirm the theoretical expectation that the acoustic backscatter from the investigated organisms were confined to the Rayleigh and Geometric scattering regions, a finding that might both aid acoustic identification and size-group separation of in situ northern krill.