Iterative tomographic analysis based on automatic refined picking

The ever-growing size of data sets for active and passive seismic imaging makes the availability of automatic procedures for rapid analysis more and more valuable. Such procedures are especially important for time-critical applications like emergency decisions or re-orienting of ongoing seismic surveys. In this paper a new, iterative scheme for 3D traveltime tomography is presented. The technique, based on a tool originally developed for earthquake data, uses cross-correlation to examine waveform similarity and to adjust arrival times on seismic sections. A preliminary set of reference arrival times is first corrected by the cross-correlation lag and then used to build an initial 3D tomographic velocity model through a standard inversion code; traveltimes calculated from this model are then taken as new reference arrivals and the process of pick adjustment is repeated. The result is a tomographic image, upgraded and refined at each iteration of the procedure. The test performed on the waveform data set recorded during the 2001 SERAPIS active seismic survey in the gulfs of Naples and Pozzuoli (Southern Italy) shows that the 3D iterative tomography scheme produces a velocity image of the structure of the Campi Flegrei caldera which is consistent with the results from previous studies, employing just a fraction of the time needed by a human analyst to identify first breaks. We believe that this technique can be effectively employed for rapid analysis of large data-sets within time-critical or time-dependent tasks and for automatic 4D tomographic investigations.     

Summary findings of the fourth international radiocarbon intercomparison (FIRI)(1998–2001)

Interlaboratory comparisons have been widely used in applied radiocarbon science. These are an important part of ongoing quality assurance (QA) programmes, which are vital to the appropriate interpretation of the evidence provided by the ¹⁴C record in Quaternary applications (including climate change and environmental reconstruction). International comparisons of laboratory performance are an essential component of the quality assurance process in radiocarbon dating. If the user community is to have confidence in radiocarbon results, it needs to be assured that laboratories world wide are producing measurements that are reliable and in accordance with ‘good practice’. The findings from the most recent (completed in 2001) and extensive (more than 90 participating laboratories) radiocarbon intercomparison (FIRI) are reported here. This study was designed (i) to assess comparability, or otherwise, of the results from different laboratories and (ii) to quantify the extent and possible causes of any interlaboratory variation. The results demonstrate that there are no significant differences amongst the main measurement techniques (gas proportional counting, liquid scintillation counting and accelerator mass spectrometry (AMS)) but there is evidence of small laboratory offsets relative to known age samples for some laboratories. There is also evidence in some cases of underestimation of measurement precision. Approximately 10% of all results were classified as extreme (outliers) and these results were generated by 14% of the laboratories. Overall, the evidence supports the fact that radiocarbon laboratories are generally accurate and precise but that, notwithstanding internal QA procedures, some problems still occur, which can best be detected by participation in independent intercomparisons such as FIRI, where the results allow individual laboratories to assess their performance and to take remedial measures where necessary. The results from FIRI are significant in that they show a broad measure of agreement between measurements made in different laboratories on a wide range of materials and they also demonstrate no statistically significant difference between measurements made by radiometric or AMS techniques. Copyright © 2002 John Wiley & Sons, Ltd.     

Regional Fluxes Of Momentum And Sensible Heat Over A Sub-Arctic Landscape During Late Winter

Based on measurements at Sodankylä Meteorological Observatory the regional (aggregated) momentum and sensible heat fluxes are estimated for two days over a site in Finnish Lapland during late winter. The forest covers 49% of the area. The study shows that the forest dominates and controls the regional fluxes of momentum and sensible heat in different ways. The regional momentum flux is found to be 10–20% smaller than the measured momentum flux over the forest, and the regional sensible heat flux is estimated to be 30–50% of the values measured over a coniferous forest.The regional momentum flux is determined in two ways, both based on blending height theory. One is a parameterised method, the other represents a numerical solution of an aggregation model. The regional sensible heat flux is determined from the theory of mixed-layer growth. At near neutral conditions the regional momentum flux can be determined independently of the regional sensible heat flux. At unstable conditions the two models become coupled.The information that is needed by the parameterised blending height method and by the mixed-layer evolution method in order to derive the regional fluxes of momentum and sensible heat can be obtained from radiosonde profiles of wind speed and temperature.     

Current and future atmospheric circulation at 500 hPa over Greenland simulated by the CMIP3 and CMIP5 global models

The Greenland ice sheet is projected to be strongly affected by global warming. These projections are either issued from downscaling methods (such as Regional Climate Models) or they come directly from General Circulation Models (GCMs). In this context, it is necessary to evaluate the accuracy of the daily atmospheric circulation simulated by the GCMs, since it is used as forcing for downscaling methods. Thus, we use an automatic circulation type classification based on two indices (Euclidean distance and Spearman rank correlation using the daily 500 hPa geopotential height) to evaluate the ability of the GCMs from both CMIP3 and CMIP5 databases to simulate the main circulation types over Greenland during summer. For each circulation type, the GCMs are compared to three reanalysis datasets on the basis of their frequency and persistence differences. For the current climate (1961–1990), we show that most of the GCMs do not reproduce the expected frequency and the persistence of the circulation types and that they simulate poorly the observed daily variability of the general circulation. Only a few GCMs can be used as reliable forcings for downscaling methods over Greenland. Finally, when applying the same approach to the future projections of the GCMs, no significant change in the atmospheric circulation over Greenland is detected, besides a generalised increase of the geopotential height due to a uniform warming of the atmosphere.     

Magnetic resonance sounding applied to aquifer characterization

Magnetic resonance sounding (MRS) is distinguished from other geophysical tools used for ground water investigation by the fact that it measures a magnetic resonance signal generated directly from subsurface water molecules. An alternating current pulse energizes a wire loop on the ground surface and the MRS signal is generated; subsurface water is indicated, with a high degree of reliability, by nonzero amplitude readings. Measurements with varied pulse magnitudes then reveal the depth and thickness of water saturated layers. The hydraulic conductivity of aquifers can also be estimated using boreholes for calibration. MRS can be used for both predicting the yield of water supply wells and for interpolation between boreholes, thereby reducing the number of holes required for hydrogeological modeling. An example of the practical application of MRS combined with two-dimensional electrical imaging, in the Kerbernez and Kerien catchments area of France, demonstrates the efficiency of the technique.     

Chalcophile element geochemistry of the Boggy Plain zoned pluton, southeastern Australia: a S-saturated barren compositionally diverse magmatic system

The behavior of the platinum group elements (PGE) and Re in felsic magmas is poorly understood due to scarcity of data. We report the concentrations of Ni, Cu, Re, and PGE in the compositionally diverse Boggy Plain zoned pluton (BPZP), which shows a variation of rock type from gabbro through granodiorite and granite to aplite with a SiO₂ range from 52 to 74 wt %. In addition, major silicate and oxide minerals were analyzed for Ni, Cu, and Re, and a systematic sulfide study was carried out to investigate the role of silicate, oxide, and sulfide minerals on chalcophile element geochemistry of the BPZP. Mass balance calculation shows that the whole rock Cu budget hosted by silicate and oxide minerals is <13 wt % and that Cu is dominantly located in sulfide phases, whereas most of the whole rock Ni budget (>70 wt %) is held in major silicate and oxide minerals. Rhenium is dominantly hosted by magnetite and ilmenite. Ovoid-shaped sulfide blebs occur at the boundary between pyroxene phenocrysts and neighboring interstitial phases or within interstitial minerals in the gabbro and the granodiorite. The blebs are composed of pyrrhotite, pyrite, chalcopyrite, and S-bearing Fe-oxide, which contain total trace metals (Co, Ni, Cu, Ag, Pb) up to ~16 wt %. The mineral assemblage, occurrence, shape, and composition of the sulfide blebs are a typical of magmatic sulfide. PGE concentrations in the BPZP vary by more than two orders of magnitude from gabbro (2.7–7.8 ppb Pd, 0.025–0.116 ppb Ir) to aplite (0.05 ppb Pd, 0.001 ppb Ir). Nickel, Cu, Re, and PGE concentrations are positively correlated with MgO in all the rock types although there is a clear discontinuity between the granodiorite and the granite in the trends for Ni, Rh, and Ir when plotted against MgO. Cu/Pd values gradually increase from 6,100 to 52,600 as the MgO content decreases. The sulfide petrology and chalcophile element geochemistry of the BPZP show that sulfide saturation occurred in the late gabbroic stage of magma differentiation. Segregation and distribution of these sulfide blebs controlled Cu and PGE variations within the BPZP rocks although the magma of each rock type may have experienced a different magma evolution history in terms of crustal assimilation and crystal fractionation. The sulfide melt locked in the cumulate rocks must have sequestered a significant portion of the chalcophile elements, which restricted the availability of these metals to magmatic-hydrothermal ore fluids. Therefore, we suggest that the roof rocks that overlay the BPZP were not prospective for magmatic-hydrothermal Cu, Au, or Cu–Au deposits.     

An integrated zircon geochronological and geochemical investigation into the Miocene plutonic evolution of the Cyclades, Aegean Sea, Greece: Part 1: Geochronology

We use 369 individual U–Pb zircon ages from 14 granitoid samples collected on five islands in the Cyclades in the Aegean Sea, Greece, for constraining the crystallisation history of I- and S-type plutons above the retreating Hellenic subduction zone. Miocene magmatism in the Cyclades extended over a time span from 17 to 11 Ma. The ages for S-type granites are systematically ~2 million years older than those for I-type granites. Considering plutons individually, the zircon data define age spectra ranging from simple and unimodal to complex and multimodal. Seven of the 14 investigated samples yield more than one distinct zircon crystallisation age, with one I-type granodiorite sample from Mykonos Island representing the most complex case with three resolvable age peaks. Two samples from S-type granites on Ikaria appear to have crystallised zircon over 2–3 million years, whereas for the majority of individual samples with multiple zircon age populations the calculated ages deviate by 1–1.5 million years. We interpret our age data to reflect a protracted history involving initial partial melting at deeper lithospheric levels, followed by crystallisation and cooling at shallower crustal levels. Our study corroborates published research arguing that pluton construction is due to incremental emplacement of multiple magma pulses over a few million years. Assuming that multiple age peaks of our 14 samples can indeed serve to quantify time spans for magmatic emplacement, our data suggest that Aegean plutons were constructed over a few million years. Our tectonic interpretation of the U–Pb ages is that the S-type granites resulted from partial melting and migmatisation of the lower crust, possibly starting at ~23 Ma. The I-type granites and associated mafic melts are interpreted to reflect the magmatic arc stage in the Cyclades starting at ~15 Ma.     

Macrofaunal involvement in the sublittoral decay of kelp debris: The sea urchin Psammechinus miliaris (Gmelin) (Echinodermata: Echinoidea)

Psammechinus miliaris occurs in the Clyde Sea area in large numbers (<18 individuals per 100 g^?1 weed dry wt) on sublittoral beds of detached Laminaria saccharina. Its rôle in weed decomposition has been examined by comparing its responses (behavioural choice, growth rate, absorption efficiencies of both carbon and protein, gut retention times and rate of faecal output) to fresh and rotting weed.Younger urchins grew faster than older individuals on a diet of rotting weed but not on fresh weed. Large seasonal variation existed, however, with fast growth occurring in June–August and little, or no, growth in December–February, irrespective of diet. Starved controls did not grow. Correcting for seasonality, rotting kelp still promoted faster growth of young urchins than did fresh weed. Larger (older) individuals showed no difference. Urchins fed fresh weed had significantly longer gut retention times. Protein absorption efficiency was higher on fresh than rotting weed, varying with weed protein content and size of urchin. Very young individuals can only digest high protein weed efficiently, eg. material derived from near the frond meristem. Organic carbon content of rotting weed was significantly lower than fresh weed. Carbon absorption efficiencies were significantly higher on fresh weed which related to organic carbon content. Standard-sized urchins fed rotting weed produced larger dry weights of faeces per day, reflecting increased ingestion rate. In closed-system choice experiments urchins preferred rotting weed kinetically. Size-frequency analysis of field populations suggested that weed beds are principally colonized by larval settlement from the plankton.Mature Psammechinus have evolved different ‘strategies’ for exploiting fresh and rotting weed. Fresh weed is relatively difficult to digest and long gut retention times allow high protein absorption efficiencies to be attained. Rotting weed has microbial protein in quantities and a lower organic carbon fraction. Some bacterial protein is seemingly unavailable though and lower protein absorption efficiencies result. Thus gut retention time is shortened and more food passed through the gut. Growth remains equivalent. Substratum digestion is of paramount importance for Psammechinus feeding on either fresh or rotting weed, cf. the ‘classical’ microbe-stripping detritivore of Fenchel.     

Experimental simulation of plagioclase diagenesis atP-T conditions of 3.5 km burial depth

Dissolution of plagioclase under the physical conditions at shallow to intermediate burial depths is a prime candidate for secondary porosity generation in feldspathic siliciclastic sediments. The diagenetic behavior of granular aggregates of plagioclase feldspar and quartz has been investigated by experimentation performed in a Bridgeman-type pressure vessel. The experiments, each of two weeks duration, simulated pressure-temperature conditions approximating 3.5 km burial depth. By using a double-acting pore-fluid reservoir, solutions of various chemistries were cycled through samples composed of oligoclase or labradorite feldspar and quartz (90:10 wt% respectively).Scanning electron microscope analysis of the post-experiment samples reveals dissolution features and precipitated products. Dissolution voids of 10 microns occur typically in areas of maximum stress such as crack-tips and grain contacts. Dissolution on a larger scale is exemplified by topographical smoothing of grain su faces. The dissolved species are subsequently reprecipitated as Ca-enriched overgrowths (possibly zeolites) and clays. These precipitates are found individually on the scale of 10 microns and collectively as surface coatings on both feldspar and quartz grains. Atomic absorption spectroscopic analyses of the pore fluid suggest that the fluid chemistry is consistent with the observed experimental precipitates.These experiments show that clay coatings are unnecessary precursors to grain surface dissolution and that the diagenetic precipitation is not mineral selective. Also, the mass transfer of the dissolved species appears to be localized because grains displaying both dissolution and precipitation features are commonplace. Volume changes due to mineral transformation/alteration may increase secondary porosity if the dissolved species produced from dissolution are only partially involved in reprecipitation and the remaining dissolved material is flushed out by the pore fluids. However, if the mass transfer is primarily local then permeability would significantly decrease as precipitates may choke the pore throats.