ALBERT

Description: We describe a multihomogeneity theory for source-parameter estimation of potential fields. Similar to what happens for random source models, where the monofractal scaling-law has been generalized into a multifractal law, we propose to generalize the homogeneity law into a multihomogeneity law. This allows a theoretically correct approach to study real-world potential fields, which are inhomogeneous and so do not show scale invariance, except in the asymptotic regions (very near to or very far from their sources). Since the scaling properties of inhomogeneous fields change with the scale of observation, we show that they may be better studied at a set of scales than at a single scale and that a multihomogeneous model is needed to explain its complex scaling behaviour. In order to perform this task, we first introduce fractional-degree homogeneous fields, to show that: (i) homogeneous potential fields may have fractional or integer degree; (ii) the source-distributions for a fractional-degree are not confined in a bounded region, similarly to some integer-degree models, such as the infinite line mass and (iii) differently from the integer-degree case, the fractional-degree source distributions are no longer uniform density functions. Using this enlarged set of homogeneous fields, real-world anomaly fields are studied at different scales, by a simple search, at any local window W , for the best homogeneous field of either integer or fractional-degree, this yielding a multiscale set of local homogeneity-degrees and depth estimations which we call multihomogeneous model. It is so defined a new technique of source parameter estimation (Multi-HOmogeneity Depth Estimation, MHODE), permitting retrieval of the source parameters of complex sources. We test the method with inhomogeneous fields of finite sources, such as faults or cylinders, and show its effectiveness also in a real-case example. These applications show the usefulness of the new concepts, multihomogeneity and fractional homogeneity-degree, to obtain valid estimates of the source parameters in a consistent theoretical framework, so overcoming the limitations imposed by global-homogeneity to widespread methods, such as Euler deconvolution.

Description: In this paper, two separate but related goals are tackled. The first one is to demonstrate that in some saturated rock textures the non-linear behaviour of induced polarization (IP) and the violation of Ohm's law not only are real phenomena, but they can also be satisfactorily predicted by a suitable physical-mathematical model, which is our second goal. This model is based on Fick's second law. As the model links the specific dependence of resistivity and chargeability of a laboratory sample to the injected current and this in turn to its pore size distribution, it is able to predict pore size distribution from laboratory measurements, in good agreement with mercury injection capillary pressure test results. This fact opens up the possibility for hydrogeophysical applications on a macro scale. Mathematical modelling shows that the chargeability acquired in the field under normal conditions, that is at low current, will always be very small and approximately proportional to the applied current. A suitable field test site for demonstrating the possible reliance of both resistivity and chargeability on current was selected and a specific measuring strategy was established. Two data sets were acquired using different injected current strengths, while keeping the charging time constant. Observed variations of resistivity and chargeability are in agreement with those predicted by the mathematical model. These field test data should however be considered preliminary. If confirmed by further evidence, these facts may lead to changing the procedure of acquiring field measurements in future, and perhaps may encourage the design and building of a new specific geo-resistivity meter. This paper also shows that the well-known Marshall and Madden's equations based on Fick's law cannot be solved without specific boundary conditions.

Description: We present a new 3-D traveltime tomography code (TOMO3D) for the modelling of active-source seismic data that uses the arrival times of both refracted and reflected seismic phases to derive the velocity distribution and the geometry of reflecting boundaries in the subsurface. This code is based on its popular 2-D version TOMO2D from which it inherited the methods to solve the forward and inverse problems. The traveltime calculations are done using a hybrid ray-tracing technique combining the graph and bending methods. The LSQR algorithm is used to perform the iterative regularized inversion to improve the initial velocity and depth models. In order to cope with an increased computational demand due to the incorporation of the third dimension, the forward problem solver, which takes most of the run time (~90 per cent in the test presented here), has been parallelized with a combination of multi-processing and message passing interface standards. This parallelization distributes the ray-tracing and traveltime calculations among available computational resources. The code's performance is illustrated with a realistic synthetic example, including a checkerboard anomaly and two reflectors, which simulates the geometry of a subduction zone. The code is designed to invert for a single reflector at a time. A data-driven layer-stripping strategy is proposed for cases involving multiple reflectors, and it is tested for the successive inversion of the two reflectors. Layers are bound by consecutive reflectors, and an initial velocity model for each inversion step incorporates the results from previous steps. This strategy poses simpler inversion problems at each step, allowing the recovery of strong velocity discontinuities that would otherwise be smoothened.

Description: Iterative substitution of the coupled Marchenko equations is a novel methodology to retrieve the Green's functions from a source or receiver array at an acquisition surface to an arbitrary location in an acoustic medium. The methodology requires as input the single-sided reflection response at the acquisition surface and an initial focusing function, being the time-reversed direct wavefield from the acquisition surface to a specified location in the subsurface. We express the iterative scheme that is applied by this methodology explicitly as the successive actions of various linear operators, acting on an initial focusing function. These operators involve multidimensional crosscorrelations with the reflection data and truncations in time. We offer physical interpretations of the multidimensional crosscorrelations by subtracting traveltimes along common ray paths at the stationary points of the underlying integrals. This provides a clear understanding of how individual events are retrieved by the scheme. Our interpretation also exposes some of the scheme's limitations in terms of what can be retrieved in case of a finite recording aperture. Green's function retrieval is only successful if the relevant stationary points are sampled. As a consequence, internal multiples can only be retrieved at a subsurface location with a particular ray parameter if this location is illuminated by the direct wavefield with this specific ray parameter. Several assumptions are required to solve the Marchenko equations. We show that these assumptions are not always satisfied in arbitrary heterogeneous media, which can result in incomplete Green's function retrieval and the emergence of artefacts. Despite these limitations, accurate Green's functions can often be retrieved by the iterative scheme, which is highly relevant for seismic imaging and inversion of internal multiple reflections.

Description: Employing dynamic reciprocity can be an effective tool to simplify the calculation of elastic wavefields for borehole problems and to check the results. We analytically obtain the reciprocity relations for the elastodynamic fields generated by multipole sources in a fluid–solid configuration: if the multipole sources are located in the fluid, the particle displacement due to a dipole source is reciprocal to the particle acceleration due to a single force; the fluid pressure due to a dipole source is reciprocal to the particle acceleration due to a monopole source; the particle displacement due to a quadrupole source is reciprocal to the spatial partial derivative of the particle acceleration due to a single force; the fluid pressure due to a quadrupole source is reciprocal to the spatial partial derivative of the particle acceleration due to a monopole source. These relations are tested by numerical experiments for different borehole problems, including acoustic logging, single-well imaging and vertical seismic profiling. A reciprocity test can be used as a quick check of a finite-difference algorithm and the implementation of the sources, although it cannot detect errors due to improper discretization of the interfaces.

Description: Clonal populations accumulate mutations over time, resulting in different haplotypes. Deep sequencing of such a population in principle provides information to reconstruct these haplotypes and the frequency at which the haplotypes occur. However, this reconstruction is technically not trivial, especially not in clonal systems with a relatively low mutation frequency. The low number of segregating sites in those systems adds ambiguity to the haplotype phasing and thus obviates the reconstruction of genome-wide haplotypes based on sequence overlap information. Therefore, we present EVORhA, a haplotype reconstruction method that complements phasing information in the non-empty read overlap with the frequency estimations of inferred local haplotypes. As was shown with simulated data, as soon as read lengths and/or mutation rates become restrictive for state-of-the-art methods, the use of this additional frequency information allows EVORhA to still reliably reconstruct genome-wide haplotypes. On real data, we show the applicability of the method in reconstructing the population composition of evolved bacterial populations and in decomposing mixed bacterial infections from clinical samples.

Description: The area of the 9.1-km-deep Continental Deep Drillhole (KTB) in Germany is used as a case study for a geothermal reservoir situated in folded and faulted metamorphic crystalline crust. The presented approach is based on the analysis of 3-D seismic reflection data combined with borehole data and hydrothermal numerical modelling. The KTB location exemplarily contains all elements that make seismic prospecting in crystalline environment often more difficult than in sedimentary units, basically complicated tectonics and fracturing and low-coherent strata. In a first step major rock units including two known nearly parallel fault zones are identified down to a depth of 12 km. These units form the basis of a gridded 3-D numerical model for investigating temperature and fluid flow. Conductive and advective heat transport takes place mainly in a metamorphic block composed of gneisses and metabasites that show considerable differences in thermal conductivity and heat production. Therefore, in a second step, the structure of this unit is investigated by seismic waveform modelling. The third step of interpretation consists of applying wavenumber filtering and log-Gabor-filtering for locating fractures. Since fracture networks are the major fluid pathways in the crystalline, we associate the fracture density distribution with distributions of relative porosity and permeability that can be calibrated by logging data and forward modelling of the temperature field. The resulting permeability distribution shows values between 10 –16 and 10 –19 m 2 and does not correlate with particular rock units. Once thermohydraulic rock properties are attributed to the numerical model, the differential equations for heat and fluid transport in porous media are solved numerically based on a finite difference approach. The hydraulic potential caused by topography and a heat flux of 54 mW m –2 were applied as boundary conditions at the top and bottom of the model. Fluid flow is generally slow and mainly occurring within the two fault zones. Thus, our model confirms the previous finding that diffusive heat transport is the dominant process at the KTB site. Fitting the observed temperature–depth profile requires a correction for palaeoclimate of about 4 K at 1 km depth. Modelled and observed temperature data fit well within 0.2 °C bounds. Whereas thermal conditions are suitable for geothermal energy production, hydraulic conditions are unfavourable without engineered stimulation.

Description: Data on biological mechanisms of aging are mostly obtained from cross-sectional study designs. An inherent disadvantage of this design is that inter-individual differences can mask small but biologically significant age-dependent changes. A serially sampled design (same individual at different time points) would overcome this problem but is often limited by the relatively small numbers of available paired samples and the statistics being used. To overcome these limitations, we have developed a new vector-based approach, termed three-component analysis, which incorporates temporal distance, signal intensity and variance into one single score for gene ranking and is combined with gene set enrichment analysis. We tested our method on a unique age-based sample set of human skin fibroblasts and combined genome-wide transcription, DNA methylation and histone methylation (H3K4me3 and H3K27me3) data. Importantly, our method can now for the first time demonstrate a clear age-dependent decrease in expression of genes coding for proteins involved in translation and ribosome function. Using analogies with data from lower organisms, we propose a model where age-dependent down-regulation of protein translation-related components contributes to extend human lifespan.

Description: A major challenge in the field of shotgun metagenomics is the accurate identification of organisms present within a microbial community, based on classification of short sequence reads. Though existing microbial community profiling methods have attempted to rapidly classify the millions of reads output from modern sequencers, the combination of incomplete databases, similarity among otherwise divergent genomes, errors and biases in sequencing technologies, and the large volumes of sequencing data required for metagenome sequencing has led to unacceptably high false discovery rates (FDR). Here, we present the application of a novel, gene-independent and signature-based metagenomic taxonomic profiling method with significantly and consistently smaller FDR than any other available method. Our algorithm circumvents false positives using a series of non-redundant signature databases and examines G enomic O rigins T hrough T axonomic CHA llenge (GOTTCHA). GOTTCHA was tested and validated on 20 synthetic and mock datasets ranging in community composition and complexity, was applied successfully to data generated from spiked environmental and clinical samples, and robustly demonstrates superior performance compared with other available tools.

Description: This paper presents the results from the Deflo-hydroacoustic experiment in the Southern Indian Ocean using three autonomous underwater hydrophones, complemented by two permanent hydroacoustic stations. The array monitored for 14 months, from November 2006 to December 2007, a 3000 x 3000 km wide area, encompassing large segments of the three Indian spreading ridges that meet at the Indian Triple Junction. A catalogue of 11 105 acoustic events is derived from the recorded data, of which 55 per cent are located from three hydrophones, 38 per cent from 4, 6 per cent from five and less than 1 per cent by six hydrophones. From a comparison with land-based seismic catalogues, the smallest detected earthquakes are m b 2.6 in size, the range of recorded magnitudes is about twice that of land-based networks and the number of detected events is 5–16 times larger. Seismicity patterns vary between the three spreading ridges, with activity mainly focused on transform faults along the fast spreading Southeast Indian Ridge and more evenly distributed along spreading segments and transforms on the slow spreading Central and ultra-slow spreading Southwest Indian ridges; the Central Indian Ridge is the most active of the three with an average of 1.9 events/100 km/month. Along the Sunda Trench, acoustic events mostly radiate from the inner wall of the trench and show a 200-km-long seismic gap between 2 °S and the Equator. The array also detected more than 3600 cryogenic events, with different seasonal trends observed for events from the Antarctic margin, compared to those from drifting icebergs at lower (up to 50°S) latitudes. Vocalizations of five species and subspecies of large baleen whales were also observed and exhibit clear seasonal variability. On the three autonomous hydrophones, whale vocalizations dominate sound levels in the 20–30 and 100 Hz frequency bands, whereas earthquakes and ice tremor are a dominant source of ambient sound at frequencies 〈20 Hz.