ALBERT

Description: Seismic waves sensitive to the outermost part of the Earth's liquid core seem to be affected by a stably stratified layer at the core–mantle boundary. Such a layer could have an observable signature in both long-term and short-term variations of the magnetic field of the Earth, which are used to probe the flow at the top of the core. Indeed, with the recent SWARM mission, it seems reasonable to be able to identify waves propagating in the core with period of several months, which may play an important role in the large-scale dynamics. In this paper, we characterize the influence of a stratified layer at the top of the core on deep quasi-geostrophic (Rossby) waves. We compute numerically the quasi-geostrophic eigenmodes of a rapidly rotating spherical shell, with a stably stratified layer near the outer boundary. Two simple models of stratification are taken into account, which are scaled with commonly adopted values of the Brunt–Väisälä frequency in the Earth's core. In the absence of magnetic field, we find that both azimuthal wavelength and frequency of the eigenmodes control their penetration into the stratified layer: the higher the phase speed, the higher the permeability of the stratified layer to the wave motion. We also show that the theory developed by Takehiro & Lister for thermal convection extends to the whole family of Rossby waves in the core. Adding a magnetic field, the penetrative behaviour of the quasi-geostrophic modes (the so-called fast branch) is insensitive to the imposed magnetic field and only weakly sensitive to the precise shape of the stratification. Based on these results, the large-scale and high-frequency modes (1–2 month periods) may be detectable in the geomagnetic data measured at the Earth's surface, especially in the equatorial area where the modes can be trapped.

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: 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: A 3-D magnetotelluric (MT) inversion code using unstructured tetrahedral elements has been developed in order to correct the topographic effect by directly incorporating it into computational grids. The electromagnetic field and response functions get distorted at the observation sites of MT surveys because of the undulating surface topography, and without correcting this distortion, the subsurface structure can be misinterpreted. Of the two methods proposed to correct the topographic effect, the method incorporating topography explicitly in the inversion is applicable to a wider range of surveys. For forward problems, it has been shown that the finite element method using unstructured tetrahedral elements is useful for the incorporation of topography. Therefore, this paper shows the applicability of unstructured tetrahedral elements in MT inversion using the newly developed code. The inversion code is capable of using the impedance tensor, the vertical magnetic transfer function (VMTF), and the phase tensor as observational data, and it estimates the subsurface resistivity values and the distortion tensor of each observation site. The forward part of the code was verified using two test models, one incorporating topographic effect and one without, and the verifications showed that the results were almost the same as those of previous works. The developed inversion code was then applied to synthetic data from a MT survey, and was verified as being able to recover the resistivity structure as well as other inversion codes. Finally, to confirm its applicability to the data affected by topography, inversion was performed using the synthetic data of the model that included two overlapping mountains. In each of the cases using the impedance tensor, the VMTF and the phase tensor, by including the topography in the mesh, the subsurface resistivity was determined more proficiently than in the case using the flat-surface mesh. Although the locations of the anomalies were not accurately estimated by the inversion using distorted impedance tensors due to the slightly undervalued gain, these locations were correctly estimated by using undistorted impedance tensors or adding VMTFs in the data. Therefore, it can be concluded that the inversion using the unstructured tetrahedral element effectively prevents the misinterpretation of subsurface resistivity and recovers subsurface resistivity proficiently by representing the topography in the computational mesh.

Description: Separating the contribution of different hematite coercivity grains to the magnetic fabric is a standing problem in rock magnetism because of the common occurrence of thermochemical alterations when measuring the anisotropy of thermal remanence. A technique that eliminates this bias is presented, which is useful when there is a need to separate the fabric of detrital from pigmentary hematite, for example. The method is based on stepwise thermal demagnetization of saturation isothermal remanent magnetizations (IRMs) applied orthogonally on three sister specimens, allowing calculation of the anisotropy tensor from the three components of each demagnetized IRM vector, avoiding the necessity of having to apply IRMs to thermochemically altered specimens. Vector subtraction allows determining the anisotropy tensor for specific unblocking-temperature ranges. The anisotropies of the pigmentary, specular and total hematite of the Mauch Chunk Formation red beds of Pennsylvania have been measured from an oriented block sample and results are compared to previous anisotropy measurements performed using the high-field anisotropy of isothermal remanence technique (hf-AIR), which measures total undifferentiated hematite. Experiments were conducted using non-saturating 1 T and fully saturating 5.5 T fields: both experimental sets seem capable of measuring the orientation of the specularite anisotropy principal axes, but 5.5 T are needed to capture the orientation of the higher coercivity pigmentary grains. The magnitudes of the principal axes, instead, are only faithfully measured using 5.5 T fields and yield somewhat higher anisotropies than those measured by hf-AIR. The fundamental requirement for this technique is homogeneous material among the three sister specimens, which is a significant limitation; homogeneity tests allow assessment of applicability of the method and reliability of the results.

Description: Archaeomagnetic field models cover longer timescales than historical models and may therefore resolve the motion of geomagnetic features on the core–mantle boundary (CMB) in a more meaningful statistical sense. Here we perform a detailed appraisal of archaeomagnetic field models to infer some aspects of the physics of the outer core. We characterize and compare the identification and tracking of reversed flux patches (RFPs) in order to assess the RFPs robustness. We find similar behaviour within a family of models but differences among different families, suggesting that modelling strategy is more influential than data set. Similarities involve recurrent positions of RFPs, but no preferred direction of motion is found. The tracking of normal flux patches shows similar qualitative behaviour confirming that RFPs identification and tracking is not strongly biased by their relative weakness. We also compare the tracking of RFPs with that of the historical field model gufm1 and with seismic anomalies of the lowermost mantle to explore the possibility that RFPs have preferred locations prescribed by lower mantle lateral heterogeneity. The archaeomagnetic field model that most resembles the historical field is interpreted in terms of core dynamics and core–mantle thermal interactions. This model exhibits correlation between RFPs and low seismic shear velocity in co-latitude and a shift in longitude. These results shed light on core processes, in particular we infer toroidal field lines with azimuthal orientation below the CMB and large fluid upwelling structures with a width of about 80° (Africa) and 110° (Pacific) at the top of the core. Finally, similar preferred locations of RFPs in the past 9 and 3 kyr of the same archaeomagnetic field model suggest that a 3 kyr period is sufficiently long to reliably detect mantle control on core dynamics. This allows estimating an upper bound of 220–310 km for the magnetic boundary layer thickness below the CMB.

Description: Post-depositional reductive diagenesis usually results in partial or entire cleansing of the pristine palaeomagnetic signal, therefore, its intensity is important to be assessed for sediments that are in the purpose of retrieving palaeomagnetic information. Grain size, rock magnetic and geochemical studies on the entire core, along with scanning electron microscope observations and X-ray diffraction analyses for representative samples were carried out on a Holocene sediment core retrieved from the deep water part of Huguangyan maar lake (HGY), southeast China. The pristine magnetic mineral assemblage of the studied core is domianted by superparamagnetic (SP) and stable single domain titanomagnetite, and high coercivity minerals are not detectable. Based on down-core variations of the average grain size ( M Z ), total organic carbon (TOC), detrital elements (Al, Ti, Fe and Mn) and the concentration and mineralogy of magnetic minerals, the studied core could be divided into three subsections. The uppermost subsection is the least affected by diagenesis, with detrital titanomagnetite as the dominant magnetic mineral. This is owing to low TOC contents, but high detrital input generated by weak Asian summer monsoon intensity during the late Holocene. The intermediate subsection shows down-core progressively enhanced dissolution of detrital titanomagnetite, and concomitant formation of authigenic pyrite and siderite, which indicates down-core progressively enhanced diagenesis generated by down-core progressive increasing TOC content, but decreasing detrital input as the result of down-core progressively strengthened Asian summer monsoon intensity. The pristine magnetic mineral assemblage has been profoundly modified in the lowermost subsection. At certain positions of the lowermost subsection, detrital titanomagnetite has been even completely dissolved via diagenesis, giving place to authigenic pyrite and siderite. High TOC content, but low detrital input generated from strong Asian summer monsoon intensity during the early Holocene are accountable for intensive diagenesis in the lowermost subsection. Complete erasing of detrital magnetic input signal at certain positions of the lowermost subsection, and considerable formation of authigenic siderite indicate that palaeomagnetic records of the studied core have been significantly compromised. The studied core has relatively higher TOC content, lower detrital matter content, calmer sedimentary environments, and less DO available at its water–sediment interface than the cores retrieved at relatively shallower water depths, which all contribute to its relatively stronger diagenesis. Progressive thickening of the upper two subsections with increasing water depth is owing to progressive increase in sedimentation rate with increasing water depth, which is the key factor in determining the thickness of each diagenetic subsection of cores from HGY. It would be better that lake sediments for palaeomagnetic investigations collected at a water depth shallower than the depth of its thermocline.

Description: The rocks in the crust and the upper mantle of the Earth are believed to exhibit electrical anisotropy to some extent. It is beneficial to further understand and recognize the propagation of the electromagnetic waves in the Earth by investigating the magnetotelluric (which is one of the main geophysical techniques to probe the deep structures in the Earth) responses of the media with anisotropic conductivity structures. In this study, we examine the magnetotelluric fields over an idealized 2-D model consisting of two segments with axially anisotropic conductivity structures overlying a perfect conductor basement by a quasi-static analytic approach. The resulting analytic solution could not only contribute to the electromagnetic induction theory in the anisotropic Earth but also serve as at least an initial standard solution which could be used to validate the reliability and accuracy of the numerical algorithms developed for modelling the magnetotelluric responses of the 2-D media with much more general anisotropic conductivity.

Description: Several works have reported that haematite has non-linear initial susceptibility at room temperature, like pyrrhotite or titanomagnetite, but there is no explanation for the observed behaviours yet. This study sets out to determine which physical property (grain size, foreign cations content and domain walls displacements) controls the initial susceptibility. The performed measurements include microprobe analysis to determine magnetic phases different to haematite; initial susceptibility (300 K); hysteresis loops, SIRM and backfield curves at 77 and 300 K to calculate magnetic parameters and minor loops at 77 K, to analyse initial susceptibility and magnetization behaviours below Morin transition. The magnetic moment study at low temperature is completed with measurements of zero field cooled–field cooled and AC susceptibility in a range from 5 to 300 K. The minor loops show that the non-linearity of initial susceptibility is closely related to Barkhausen jumps. Because of initial magnetic susceptibility is controlled by domain structure it is difficult to establish a mathematical model to separate magnetic subfabrics in haematite-bearing rocks.

Description: Chromosome-long haplotyping of human genomes is important to identify genetic variants with differing gene expression, in human evolution studies, clinical diagnosis, and other biological and medical fields. Although several methods have realized haplotyping based on sequencing technologies or population statistics, accuracy and cost are factors that prohibit their wide use. Borrowing ideas from group testing theories, we proposed a clone-based haplotyping method by overlapping pool sequencing. The clones from a single individual were pooled combinatorially and then sequenced. According to the distinct pooling pattern for each clone in the overlapping pool sequencing, alleles for the recovered variants could be assigned to their original clones precisely. Subsequently, the clone sequences could be reconstructed by linking these alleles accordingly and assembling them into haplotypes with high accuracy. To verify the utility of our method, we constructed 130 110 clones in silico for the individual NA12878 and simulated the pooling and sequencing process. Ultimately, 99.9% of variants on chromosome 1 that were covered by clones from both parental chromosomes were recovered correctly, and 112 haplotype contigs were assembled with an N50 length of 3.4 Mb and no switch errors. A comparison with current clone-based haplotyping methods indicated our method was more accurate.

Description: Sexual differentiation of malaria parasites into gametocytes in the vertebrate host and subsequent gamete fertilization in mosquitoes is essential for the spreading of the disease. The molecular processes orchestrating these transitions are far from fully understood. Here, we report the first transcriptome analysis of male and female Plasmodium falciparum gametocytes coupled with a comprehensive proteome analysis. In male gametocytes there is an enrichment of proteins involved in the formation of flagellated gametes; proteins involved in DNA replication, chromatin organization and axoneme formation. On the other hand, female gametocytes are enriched in proteins required for zygote formation and functions after fertilization; protein-, lipid- and energy-metabolism. Integration of transcriptome and proteome data revealed 512 highly expressed maternal transcripts without corresponding protein expression indicating large scale translational repression in P. falciparum female gametocytes for the first time. Despite a high degree of conservation between Plasmodium species, 260 of these ‘repressed transcripts’ have not been previously described. Moreover, for some of these genes, protein expression is only reported in oocysts and sporozoites indicating that repressed transcripts can be partitioned into short- and long-term storage. Finally, these data sets provide an essential resource for identification of vaccine/drug targets and for further mechanistic studies.

Description: A series of important geological events occurred in the Tibetan Plateau area during the Jurassic, such as the collision of the Lhasa and Qiangtang Terranes, the closure of the Meso-Tethyan Ocean, the opening of the Neo-Tethyan Ocean and the cessation of the mega-monsoon. The ~3000 m thick Jurassic sedimentary sequence in the Qiangtang Basin on the central Tibetan Plateau, which is called the Yanshiping (YSP) Group, recorded these geological events. However, the chronology of the sequence is surprisingly poorly constrained. Here, we perform a detailed palaeomagnetic analysis on the ~1060 m thick middle and upper portions of the YSP Group (the Xiali and Suowa Formations) in the YSP section of the eastern Qiangtang Basin. Three bivalve zones at stratigraphic intervals of ~40–140, 640–800 and 940–1040 m are identified, which yield a Bathonian–Callovian age for the Lower Xiali Fm., a Callovian–Oxfordian age for the Lower Suowa Fm. and an Oxfordian–Kimmeridgian age for the Upper Suowa Fm. A total of 544 oriented palaeomagnetic samples were collected from the section. By combining thermal and alternating field demagnetizations, clear characteristic remanent magnetization (ChRM) directions are isolated for most of the samples. The robust ChRM directions pass fold and reversals tests, which support the primary nature of the ChRMs and yield a palaeopole at 76.8°N/297.2°E (dp = 2.2°, dm = 3.7°). A total of 27 normal and 26 reversed polarity zones were successfully recorded in the section. Combined with fossil age constraints, results suggest that the section is plausibly composed of a Callovian-Early Kimmeridgian age sedimentary sequence.

Description: DNA methylation is a mechanism for long-term transcriptional regulation and is required for normal cellular differentiation. Failure to properly establish or maintain DNA methylation patterns leads to cell dysfunction and diseases such as cancer. Identifying DNA methylation signatures in complex tissues can be challenging owing to inaccurate cell enrichment methods and low DNA yields. We have developed a technique called laser capture microdissection-reduced representation bisulfite sequencing (LCM-RRBS) for the multiplexed interrogation of the DNA methylation status of cytosine–guanine dinucleotide islands and promoters. LCM-RRBS accurately and reproducibly profiles genome-wide methylation of DNA extracted from microdissected fresh frozen or formalin-fixed paraffin-embedded tissue samples. To demonstrate the utility of LCM-RRBS, we characterized changes in DNA methylation associated with gonadectomy-induced adrenocortical neoplasia in the mouse. Compared with adjacent normal tissue, the adrenocortical tumors showed reproducible gains and losses of DNA methylation at genes involved in cell differentiation and organ development. LCM-RRBS is a rapid, cost-effective, and sensitive technique for analyzing DNA methylation in heterogeneous tissues and will facilitate the investigation of DNA methylation in cancer and organ development.

Description: This study investigates pressure effects on the magnetic properties of non-interacting single-domain (SD) magnetite. Using a high-pressure cell specially designed for a Magnetic Property Measurement System, magnetic hysteresis measurements were conducted under high pressures of up to 1 GPa on natural plagioclase crystals containing much acicular SD magnetite. Coercivity and saturation magnetization were nearly constant with pressure, while saturation remanent magnetization and coercivity of remanence decreased with pressure at moderate rates of –8 per cent GPa –1 and –18 per cent GPa –1 , respectively. These results suggest that temperature effects govern the magnetic behaviour of acicular SD magnetite grains in the middle and lower crusts.

Description: Characterization of geomagnetic field behaviour on timescales of centuries to millennia is necessary to understand the mechanisms that sustain the geodynamo and drive its evolution. As Holocene paleomagnetic and archeomagnetic data have become more abundant, strategies for regularized inversion of modern field data have been adapted to produce numerous time-varying global field models. We evaluate the effectiveness of several approaches to inversion and data handling, by assessing both global and regional properties of the resulting models. Global Holocene field models cannot resolve Southern hemisphere regional field variations without the use of sediments. A standard data set is used to construct multiple models using two different strategies for relative paleointensity calibration and declination orientation and a selection of starting models in the inversion procedure. When data uncertainties are considered, the results are similar overall regardless of whether we use iterative calibration and reorientation, or co-estimation of the calibration and orientation parameters as part of the inversion procedure. In each case the quality of the starting model used for initial relative paleointensity calibration and declination orientation is crucial and must be based on the best absolute information available. Without adequate initial calibration the morphology of dipole moment variations can be recovered but its absolute value will be correlated with the initial intensity calibrations, an effect that might be mitigated by ensuring an appropriate fit to enough high quality absolute intensity data with low uncertainties. The declination reorientation mainly impacts regional field structure and in the presence of non-zonal fields will result in a non-zero local average. The importance of declination orientation is highlighted by inconsistencies in the West Pacific and Australian sediment records in CALS10k.1b model. Great care must also be taken to assess uncertainties associated with both paleomagnetic and age data and to evaluate the effects of poor data distribution. New consistently allocated uncertainty estimates for sediment paleomagnetic records highlight the importance of adequate uncertainties in the inversion process, as they determine the relative weighting among the data and overall normalized misfit levels which in turn influence the complexity of the inferred field models. Residual distributions suggest that the most appropriate misfit measure is the L 1 norm (minimum absolute deviation) rather than L 2 (least squares), but this seems to have relatively minor impact on the overall results. For future Holocene field modelling we see a need for comprehensive methods to assess uncertainty in individual archeomagnetic data so that these data or models derived from them can be used for reliable initial relative paleointensity calibration and declination orientation in sediments. More work will be needed to assess whether co-estimation or an iterative approach to inversion is more efficient overall. This would be facilitated by realistic and globally consistent data and age uncertainties from the paleomagnetic community.

Description: Galvanic distortion of magnetotelluric (MT) data due to small-scale surficial bodies or due to topography is one of the major factors that prevents accurate imaging of the subsurface. We present a 3-D algorithm for joint inversion of MT impedance tensor data and a frequency-independent full distortion matrix that circumvents this problem. We perform several tests of our algorithm on synthetic data affected by different amounts of distortion. These tests show that joint inversion leads to a better conductivity model compared to the inversion of the MT impedance tensor without any correction for distortion effects. For highly distorted data, inversion without any distortion correction results in strong artefacts and we cannot fit the data to the specified noise level. When the distortion is reduced, we can fit the data to an RMS of one, but still observe artefacts in the shallow part of the model. In contrast, in both cases our joint inversion can fit the data within the assumed noise level and the resulting models are comparable to the inversion of undistorted data. In addition, we show that the elements of the full distortion matrix can be well resolved by our algorithm. Finally, when inverting undistorted data, including the distortion matrix in the inversion only results in a minor loss of resolution. We therefore consider our new approach a promising tool for the general analysis of field MT data.

Description: The magnetic field that originates in the earth's core is transformed across the electrically conducting mantle before being observed, at the earth's surface or above. Assuming that the conductivity depends only on radius, it has been customary to treat the mantle as a linear time-invariant filter for the core magnetic field, with properties (as a function of the frequency ) specified by the transfer function (). An high-frequency approximation to (), which is derived from a three terms WKBJ expansion with –1/2 as small parameter, is found here to reproduce adequately, for low harmonic degrees and/or thin conducting layers, the exact solution, which is evaluated numerically. It is contrasted with the low-frequency estimation of , which consists in a perturbation procedure and in writing () as a series in powers of ( -〉 0). The low-frequency theory is applied to the magnetic variations produced by the geostrophic core flows with about 6 yr period as the phase of these flows is independently determined from their effect on the length of the day. Apart from that, the low-frequency approximation overestimates the screening by the mantle of high-frequency signals, especially the low harmonic degree ones. In practice, the attenuating factor defined from the O ( 2 ) term in the expansion of as -〉 0 cannot be retrieved from analyses of geomagnetic time-series. Application of the mantle filter theory hinges on our knowledge about the time spectrum of the magnetic field at the core surface. The low-frequency theory had been previously applied to observatory series on the assumption that geomagnetic jerks occurring in the core are rare and isolated events. Rather than following up these earlier studies, I note that the spectral density function for the second time derivative of the main magnetic field coefficients is approximately independent of in a frequency range for which the mantle has undoubtedly negligible influence. In the absence of any other information, this scaling law is extrapolated to higher frequencies.

Description: Surface nuclear magnetic resonance (surface-NMR) is a promising technique for exploring shallow subsurface aquifer structures. Surface-NMR can be applied in environments that are characterized as a 1-D layered Earth. The technique utilizes a single loop and is referred to as magnetic resonance sounding. The technique referred to as magnetic resonance tomography (MRT) allows complex 2-D aquifer structures to be explored. Currently, MRT requires multiple loops and a roll along measurement scheme, which causes long survey time. We propose a loop layout using an elongated transmitter and an in-loop receiver arrays (ETRA) to conduct a 2-D survey with just one measurement. We present a comprehensive comparison between the new layout and the common approaches based on sensitivity and resolution analyses and show synthetic and field data. The results show that ETRA generates subsurface images at sufficient resolution with significantly lower survey times than other loop layouts.

Description: The advent in high-throughput-sequencing (HTS) technologies has revolutionized conventional biodiversity research by enabling parallel capture of DNA sequences possessing species-level diagnosis. However, polymerase chain reaction (PCR)-based implementation is biased by the efficiency of primer binding across lineages of organisms. A PCR-free HTS approach will alleviate this artefact and significantly improve upon the multi-locus method utilizing full mitogenomes. Here we developed a novel multiplex sequencing and assembly pipeline allowing for simultaneous acquisition of full mitogenomes from pooled animals without DNA enrichment or amplification. By concatenating assemblies from three de novo assemblers, we obtained high-quality mitogenomes for all 49 pooled taxa, with 36 species 〉15 kb and the remaining 〉10 kb, including 20 complete mitogenomes and nearly all protein coding genes (99.6%). The assembly quality was carefully validated with Sanger sequences, reference genomes and conservativeness of protein coding genes across taxa. The new method was effective even for closely related taxa, e.g. three Drosophila spp., demonstrating its broad utility for biodiversity research and mito-phylogenomics. Finally, the in silico simulation showed that by recruiting multiple mito-loci, taxon detection was improved at a fixed sequencing depth. Combined, these results demonstrate the plausibility of a multi-locus mito-metagenomics approach as the next phase of the current single-locus metabarcoding method.

Description: Chromatin modifiers and histone modifications are components of a chromatin-signaling network involved in transcription and its regulation. The interactions between chromatin modifiers and histone modifications are often unknown, are based on the analysis of few genes or are studied in vitro . Here, we apply computational methods to recover interactions between chromatin modifiers and histone modifications from genome-wide ChIP-Seq data. These interactions provide a high-confidence backbone of the chromatin-signaling network. Many recovered interactions have literature support; others provide hypotheses about yet unknown interactions. We experimentally verified two of these predicted interactions, leading to a link between H4K20me1 and members of the Polycomb Repressive Complexes 1 and 2. Our results suggest that our computationally derived interactions are likely to lead to novel biological insights required to establish the connectivity of the chromatin-signaling network involved in transcription and its regulation.

Description: We studied ionospheric responses to the 2012 April 11 M w 8.6 North Sumatra earthquake using total electron content (TEC) measurements with the regional Global Navigation Satellite System network. This earthquake ruptured the oceanic lithosphere off the Indian Ocean coast of North Sumatra, and is known as the largest strike-slip earthquake ever recorded. Coseismic ionospheric disturbances (CIDs) with rapid TEC enhancement of a few TEC units propagated northward with a speed of acoustic waves (~1 km s –1 ). Resonant atmospheric oscillation with a frequency ~4 mHz have been found as monochromatic oscillation of TEC lasting for an hour after the main shock and the largest aftershock. We compared CID amplitudes of 21 earthquakes world-wide with moment magnitudes ( M w ) 6.6–9.2. They roughly obeyed a law such that CID amplitude increases by two orders of magnitude for the M w increase of three. The 2012 North Sumatra earthquakes slightly deviated negatively from the trend possibly reflecting their strike-slip mechanisms, that is small vertical crustal movements for their magnitudes.

Description: With read lengths of currently up to 2 x 300 bp, high throughput and low sequencing costs Illumina's MiSeq is becoming one of the most utilized sequencing platforms worldwide. The platform is manageable and affordable even for smaller labs. This enables quick turnaround on a broad range of applications such as targeted gene sequencing, metagenomics, small genome sequencing and clinical molecular diagnostics. However, Illumina error profiles are still poorly understood and programs are therefore not designed for the idiosyncrasies of Illumina data. A better knowledge of the error patterns is essential for sequence analysis and vital if we are to draw valid conclusions. Studying true genetic variation in a population sample is fundamental for understanding diseases, evolution and origin. We conducted a large study on the error patterns for the MiSeq based on 16S rRNA amplicon sequencing data. We tested state-of-the-art library preparation methods for amplicon sequencing and showed that the library preparation method and the choice of primers are the most significant sources of bias and cause distinct error patterns. Furthermore we tested the efficiency of various error correction strategies and identified quality trimming (Sickle) combined with error correction (BayesHammer) followed by read overlapping (PANDAseq) as the most successful approach, reducing substitution error rates on average by 93%.

Description: RNA-seq is a sensitive and accurate technique to compare steady-state levels of RNA between different cellular states. However, as it does not provide an account of transcriptional activity per se , other technologies are needed to more precisely determine acute transcriptional responses. Here, we have developed an easy, sensitive and accurate novel computational method, iRNA-seq , for genome-wide assessment of transcriptional activity based on analysis of intron coverage from total RNA-seq data. Comparison of the results derived from iRNA-seq analyses with parallel results derived using current methods for genome-wide determination of transcriptional activity, i.e. global run-on (GRO)-seq and RNA polymerase II (RNAPII) ChIP-seq, demonstrate that iRNA-seq provides similar results in terms of number of regulated genes and their fold change. However, unlike the current methods that are all very labor-intensive and demanding in terms of sample material and technologies, iRNA-seq is cheap and easy and requires very little sample material. In conclusion, iRNA-seq offers an attractive novel alternative to current methods for determination of changes in transcriptional activity at a genome-wide level.

Description: The magnetic field of the Earth's lithosphere arises from rock magnetization contrasts that were shaped over geological times. The field can be described mathematically in spherical harmonics or with distributions of magnetization. We exploit this dual representation and assume that the lithospheric field is induced by spatially varying susceptibility values within a shell of constant thickness. By introducing a statistical assumption about the power spectrum of the susceptibility, we then derive a statistical expression for the spatial power spectrum of the crustal magnetic field for the spatial scales ranging from 60 to 2500 km. This expression depends on the mean induced magnetization, the thickness of the shell, and a power law exponent for the power spectrum of the susceptibility. We test the relevance of this form with a misfit analysis to the observational NGDC-720 lithospheric magnetic field model power spectrum. This allows us to estimate a mean global apparent induced magnetization value between 0.3 and 0.6 A m –1 , a mean magnetic crustal thickness value between 23 and 30 km, and a root mean square for the field value between 190 and 205 nT at 95 per cent. These estimates are in good agreement with independent models of the crustal magnetization and of the seismic crustal thickness. We carry out the same analysis in the continental and oceanic domains separately. We complement the misfit analyses with a Kolmogorov–Smirnov goodness-of-fit test and we conclude that the observed power spectrum can be each time a sample of the statistical one.

Description: Data on the evolution of Earth's magnetic field intensity are important for understanding the geodynamo and planetary evolution. However, the paleomagnetic record in rocks may be adversely affected by many physical processes, which must be taken into account when analysing the palaeointensity database. This is especially important in the light of an ongoing debate regarding core thermal conductivity values, and how these relate to the Precambrian geodynamo. Here, we demonstrate that several data sets in the Precambrian palaeointensity database overestimate the true paleofield strength due to the presence of non-ideal carriers of palaeointensity signals and/or viscous re-magnetizations. When the palaeointensity overestimates are removed, the Precambrian database does not indicate a robust change in geomagnetic field intensity during the Mesoproterozoic. These findings call into question the recent claim that the solid inner core formed in the Mesoproterozoic, hence constraining the thermal conductivity in the core to ‘moderate’ values. Instead, our analyses indicate that the presently available palaeointensity data are insufficient in number and quality to constrain the timing of solid inner core formation, or the outstanding problem of core thermal conductivity. Very young or very old inner core ages (and attendant high or low core thermal conductivity values) are consistent with the presently known history of Earth's field strength. More promising available data sets that reflect long-term core structure are geomagnetic reversal rate and field morphology. The latter suggests changes that may reflect differences in Archean to Proterozoic core stratification, whereas the former suggest an interval of geodynamo hyperactivity at ca. 550 Ma.

Description: Identifying large-scale structural variation in cancer genomes continues to be a challenge to researchers. Current methods rely on genome alignments based on a reference that can be a poor fit to highly variant and complex tumor genomes. To address this challenge we developed a method that uses available breakpoint information to generate models of structural variations. We use these models as references to align previously unmapped and discordant reads from a genome. By using these models to align unmapped reads, we show that our method can help to identify large-scale variations that have been previously missed.

Description: Red Clay underlying the loess-palaeosol sequences on the Chinese Loess Plateau is an eolian deposit. There is a controversy over whether magnetic susceptibility ( ) variations in Red Clay sequence can be used as an indicator of summer palaeomonsoon intensity. This study investigates the magnetic mineralogy, magnetic concentration and magnetic grain size distribution of Jiaxian Red Clay with multimagnetic methods. Our results indicate that the magnetic properties of Jiaxian Red Clay are similar to those of the Quaternary loess-palaeosol sequences, and ultrafine ferrimagnetic grains produced during pedogenesis are responsible for an increase in susceptibility, therefore the enhancement mechanism of Red Clay is similar to that of the overlying loess-palaeosol sequences. This paper explores variations in the Red Clay sequence through spatial and temporal analysis. The susceptibility variation of six sites along a NNE to SSW transect correlate to palaeoclimatic cycles, so can be used to trace the summer palaeomonsoon intensity from a spatial perspective. However, a simple loess-derived calibration function cannot be used to quantitative reconstruct the palaeomonsoon intensity variations thought time. An adjusted calibration function for palaeosols from Red Clay sequence needs to be developed, so that can be used to quantitative reconstruct palaeomonsoon intensity. Further study is necessary to develop such a transfer function.

Description: Palaeointensity experiments were carried out to a sample collection from two sections of basalt lava flow sequences of Pliocene age in north central Iceland (Chron C2An) to further refine the knowledge of the behaviour of the palaeomagnetic field. Selection of samples was mainly based on their stability of remanence to thermal demagnetization as well as good reversibility in variations of magnetic susceptibility and saturation magnetization with temperature, which would indicate the presence of magnetite as a product of deuteric oxidation of titanomagnetite. Among 167 lava flows from two sections, 44 flows were selected for the Königsberger–Thellier–Thellier experiment in vacuum. In spite of careful pre-selection of samples, an Arai plot with two linear segments, or a concave-up appearance, was often encountered during the experiments. This non-ideal behaviour was probably caused by an irreversible change in the domain state of the magnetic grains of the pseudo-single-domain (PSD) range. This is assumed because an ideal linear plot was obtained in the second run of the palaeointensity experiment in which a laboratory thermoremanence acquired after the final step of the first run was used as a natural remanence. This experiment was conducted on six selected samples, and no clear difference between the magnetic grains of the experimented and pristine sister samples was found by scanning electron microscope and hysteresis measurements, that is, no occurrence of notable chemical/mineralogical alteration, suggesting that no change in the grain size distribution had occurred. Hence, the two-segment Arai plot was not caused by the reversible multidomain/PSD effect in which the curvature of the Arai plot is dependent on the grain size. Considering that the irreversible change in domain state must have affected data points at not only high temperatures but also low temperatures, f v ≥ 0.5 was adopted as one of the acceptance criteria where f v is a vectorially defined fraction of the linear segment. A measure of curvature k ' was also used to check the linearity of the selected linear segment. It was avoided, however, to reject the result out of hand by the large curvature k of the entire data points because it might still include a linear segment with a large fraction. Combining with the results of Shaw's experiments, 52 palaeointensities were obtained out of 192 specimens, or 11 flow means were obtained out of the 44 lava flows. Most of the palaeointensities were from the upper part of the lava section (Chron C2An.1n) and ranged between 30 and 66 μT. Including two results from the bottom part of the lava section, the mean virtual dipole moment for 2.5–3.5 Ma is 6.3 ± 1.4  x  10 22 Am 2 ( N  = 11), which is ~19 per cent smaller than the present-day dipole moment.

Description: The sedimentary sequence deposited during the deglaciation phase following the last glacial maximum in the Storfjorden trough, on the northwestern Barents Sea south of Svalbard, was sampled with 10 piston and gravity cores during the SVAIS and EGLACOM cruises. Three cores (SV-02, SV-03 and SV-05) collected on the upper continental slope are characterized by a thin (20–40 cm) Holocene interval and a thick (up to 4.5 m in core SV-03) late Pleistocene sequence of finely laminated fine-grained sediments that have been interpreted as plumites deposited during the Melt Water Pulse 1a (MWP-1a). Radiocarbon ages obtained at the top and bottom of this stratigraphic interval revealed that deposition occurred during less than two centuries at around 15 ka ago, with a very high sedimentary rate exceeding 3 cm a –1 . We studied the palaeomagnetic and rock magnetic properties of this interval, by taking magnetic measurements at 1 cm spacing on u-channel samples collected from the three cores. The data show that this sequence is characterized by good palaeomagnetic properties and the palaeomagnetic and rock magnetic trends may be correlated at high resolution from core to core. The obtained palaeomagnetic data therefore offer the unique opportunity to investigate in detail the rate of geomagnetic palaeosecular variation (PSV) in the high northern latitudes at a decadal scale. Notwithstanding the palaeomagnetic trends of the three cores may be closely matched, the amplitude of directional PSV and the consequent virtual geomagnetic pole (VGP) scatter (S) is distinctly higher in one core (SV-05) than in the other two cores (SV-02 and SV-03). This might result from a variable proportion of two distinct populations of magnetic minerals in core SV-05, as suggested by the variable tendency to acquire a gyromagnetic remanent magnetization at high fields during the AF demagnetization treatment. For the plumite interval of cores SV-02 and SV-03, where the magnetic mineralogy is uniform and magnetite is the main magnetic carrier, a S value of about 9° is obtained. We consider this value as a reliable approximation of palaeomagnetic secular variation at a latitude of 75°N over a time interval spanning a couple of centuries around 15 ka ago.

Description: A recent study of the Matuyama–Brunhes (M-B) geomagnetic field reversal recorded in exposed lacustrine sediments from the Sulmona Basin (Italy) provided a continuous, high-resolution record indicating that the reversal of the field direction at the terminus of the M-B boundary (MBB) occurred in less than a century, about 786 ka ago. In the sediment, thin (4–6 cm) remagnetized horizons were recognized above two distinct tephra layers—SUL2-19 and SUL2-20—that occur ~25 and ~35 cm below the MBB, respectively. Also, a faint, millimetre-thick tephra (SUL2-18) occurs 2–3 cm above the MBB. With the aim of improving the temporal resolution of the previous Sulmona MBB record and understanding the possible influence of cryptotephra on the M-B record in the Sulmona Basin, we performed more detailed sampling and analyses of overlapping standard and smaller samples from a 50 cm-long block that spans the MBB. The new data indicate that (i) the MBB is even sharper than previously reported and occurs ~2.5 cm below tephra SUL2-18, in agreement with the previous study; (ii) the MBB coincides with the rise of an intensity peak of the natural remanent magnetization (NRM) intensity, which extends across SUL2-18; (iii) except for a 2-cm-thick interval just above tephra SUL2-18, the rock magnetic parameters ( k , ARM, M r , M s , B c , B cr ) indicate exactly the same magnetic mineralogy throughout the sampled sequence. We conclude that either SUL2-18 resulted in the remagnetization of an interval of about 6 cm (i.e. during the NRM intensity peak spanning ~260 ± 110 yr, according to the estimated local sedimentation rate), and thus the detailed MBB record is lost because it is overprinted, or the MBB is well recorded, occurred abruptly about 2.5 cm below SUL2-18 and lasted less than 13 ± 6 yr. Both hypotheses challenge our understanding of the geomagnetic field behaviour during a polarity transition and/or of the NRM acquisition process in the Sulmona lacustrine sediment.

Description: We test the ability of large-scale velocity fields inferred from geomagnetic secular variation data to produce the global magnetic field of the Earth. Our kinematic dynamo calculations use quasi-geostrophic (QG) flows inverted from geomagnetic field models which, as such, incorporate flow structures that are Earth-like and may be important for the geodynamo. Furthermore, the QG hypothesis allows straightforward prolongation of the flow from the core surface to the bulk. As expected from previous studies, we check that a simple QG flow is not able to sustain the magnetic field against ohmic decay. Additional complexity is then introduced in the flow, inspired by the action of the Lorentz force. Indeed, on centennial timescales, the Lorentz force can balance the Coriolis force and strict quasi-geostrophy may not be the best ansatz. When our columnar flow is modified to account for the action of the Lorentz force, magnetic field is generated for Elsasser numbers larger than 0.25 and magnetic Reynolds numbers larger than 100. This suggests that our large-scale flow captures the relevant features for the generation of the Earth's magnetic field and that the invisible small-scale flow may not be directly involved in this process. Near the threshold, the resulting magnetic field is dominated by an axial dipole, with some reversed flux patches. Time dependence is also considered, derived from principal component analysis applied to the inverted flows. We find that time periods from 120 to 50 yr do not affect the mean growth rate of the kinematic dynamos. Finally, we note that the footprint of the inner core in the magnetic field generated deep in the bulk of the shell, although we did not include one in our computations.

Description: The Auca Mahuida volcano (2.03–0.88 Ma) located east of the Andean thrust front in the Neuquén basin (Argentina) hosts an oil system of thermogenic origin and is affected by the NW–SE striking-faults. Intrusive bodies and the underlying Jurassic sediments constitute the reservoir rocks. Aeromagnetic data collected in the Auca Mahuida area detected multiple dipolar magnetic anomalies, many of which have reverse polarity. Palaeomagnetic measurements on rock samples collected in the field together with available age determinations indicate that the reversely magnetized sources were mainly emplaced during the Matuyama reverse polarity chron while the normal polarity sources were emplaced during the Olduvai and/or Jaramillo subchrons. The location and geometry of the intrusive bodies is poorly known and the customary magnetic inversion is rendered difficult because of multiple natural remanent magnetization directions. To address these difficulties, a voxel inversion was applied to model the vector residual magnetic intensity (VRMI) transformation of the observed total magnetic intensity data. The modelling showed a 1.5 km deep, subcircular ring-shaped intrusion below the summit of the volcano and a series of NW–SE elongated, fault-controlled intrusive bodies to depths up to 3–4 km. Our results show that magnetic data and VRMI modelling help resolve the geometry of the shallow plumbing system of volcanoes with remanently magnetized sources, and estimate the depth and geometry of potential oil reservoirs in volcanic areas.

Description: Any given human individual carries multiple genetic variants that disrupt protein-coding genes, through structural variation, as well as nucleotide variants and indels. Predicting the phenotypic consequences of a gene disruption remains a significant challenge. Current approaches employ information from a range of biological networks to predict which human genes are haploinsufficient (meaning two copies are required for normal function) or essential (meaning at least one copy is required for viability). Using recently available study gene sets, we show that these approaches are strongly biased towards providing accurate predictions for well-studied genes. By contrast, we derive a haploinsufficiency score from a combination of unbiased large-scale high-throughput datasets, including gene co-expression and genetic variation in over 6000 human exomes. Our approach provides a haploinsufficiency prediction for over twice as many genes currently unassociated with papers listed in Pubmed as three commonly-used approaches, and outperforms these approaches for predicting haploinsufficiency for less-studied genes. We also show that fine-tuning the predictor on a set of well-studied ‘gold standard’ haploinsufficient genes does not improve the prediction for less-studied genes. This new score can readily be used to prioritize gene disruptions resulting from any genetic variant, including copy number variants, indels and single-nucleotide variants.

Description: Was the Yukon-Tanana Terrane (YTT), a California-sized part of south-central Yukon, an autochthonous or para-autochthonous part of northern British Columbia in the Early Cretaceous or was it part of a proposed allochthonous ‘Baja B.C.’ continent offshore of southern California? To answer this fundamental question, a paleomagnetic study has been completed on 347 specimens from 24 sites in the 114.7 ± 1.1 Ma Quiet Lake batholith. This 1300 km 2 pluton is composed mostly of massive medium-to-coarse grained biotite quartz monzonite that exhibits no evidence of either deformation or metamorphism, and that intrudes metamorphosed pre-Cretaceous basement rocks of the YTT in southern Yukon. The paleomagnetic analysis utilized thermal and alternating field step demagnetization, and saturation isothermal remanence methods. A well-defined characteristic remanent magnetization (ChRM) direction was isolated throughout the 500–585 °C temperature range at Decl. = 340.6°, Incl. = 77.4° ( N = 14 sites, k = 51.2, A 95 = 5.6°). The ChRM resides in magnetite with a low titanium content and is interpreted to be a primary thermoremanent magnetization. After correction for 490 km of geologically demonstrable dextral displacement on the inboard Tintina fault zone, the Quiet Lake batholith's paleopole is not significantly different at 95 per cent confidence from the co-eval 115 Ma reference paleopole for North America, giving non-significant translation and rotation estimates of 1.4° ± 5.1° (1) northwestwards and 10° ± 13° (1) clockwise, respectively. Thus, this is the first Early Cretaceous paleopole to show clearly that the YTT in Yukon is a para-autochthon that was part of North America's continental margin at that time. Further, after correction for Tintina fault displacement, the eight available Mesozoic YTT paleopoles agree closely with the North American apparent polar wander path (APWP). In contrast, the 22 paleopoles from the Intermontane Belt show the expected behaviour of an allochthonous thin-skin of terranes (i.e. ‘Baja B.C.’ behaviour) and record ~8° of northward translation and ~50° of clockwise rotation relative to both the YTT and North American APWPs during the Mesozoic and Paleogene.

Description: There exists a fundamental as well as practical interest in being able to accurately forecast the future evolution of Earth's magnetic field at decadal to secular ranges. This work enables such forecasts by combining geomagnetic data with an Earth-like numerical model of a convection-driven fluid dynamo. The underlying data assimilation framework builds on recent progress in inverse geodynamo modelling, a method which estimates an internal dynamic structure for Earth's core from a snapshot of the magnetic field and its instantaneous rate of change at the surface, and takes advantage of linear relationships and long-range correlations between observed and hidden state variables. Here the method is further evolved into a single-epoch ensemble Kalman filter, in order to initialise at a given epoch an ensemble of states compatible with the observations and representative of the uncertainties in the estimation of hidden quantities. The ensemble dynamics, obtained by subsequent numerical integration of the prognostic model equations, are found to be governed by a thermal wind balance or equilibrium between buoyancy forces, the Coriolis force and the pressure gradient. The resulting core fluid flow pattern is a quasi-steady eccentric gyre organised in a column parallel to Earth's rotation axis, in equilibrium with a longitudinal hemispheric convective density anomaly pattern. The flow provides induction for the magnetic field, which also undergoes a realistic amount of diffusion. Predictions of the present magnetic field from data taken within the past century show that the ensemble has an average retaining good consistency with the true geomagnetic evolution and an acceptable spread well representative of prediction errors, up to at least a secular range. The predictability of the geodynamo thus appears to significantly exceed previous theoretical expectations based on the chaotic divergence of ensemble members. The assimilation generally outperforms the linear mathematical extrapolations from a 30-yr prediction range onwards, with a 40 per cent improvement in Earth-surface error at a secular range. The geomagnetic axial dipole decay observed over the past two centuries is predicted to continue at a similar pace in the next century, with a further loss of $1.1\pm 0.3 \hspace{2.84526pt}\mathrm{\mu T}$ by year 2115. The focal (or minimum intensity) point of the South Atlantic geomagnetic anomaly is predicted to enter the South Pacific region in the next century, with the anomaly itself further deepening and widening. By year 2065, the minimum intensity is predicted to decrease by $1.46\pm 0.4 \hspace{2.84526pt}\mathrm{\mu T}$ at the Earth surface and the focal point to move 12.8 ± 1.4 deg westwards with a slight northward component. This corresponds to a drift rate of 0.26 deg yr –1 , similar to the typical geomagnetic westward drift observed over the past four centuries. The same drift rate is also predicted until 2115 with a further (but more uncertain) intensity decrease.

Description: Meta-analysis of gene expression has enabled numerous insights into biological systems, but current methods have several limitations. We developed a method to perform a meta-analysis using the elastic net, a powerful and versatile approach for classification and regression. To demonstrate the utility of our method, we conducted a meta-analysis of lung cancer gene expression based on publicly available data. Using 629 samples from five data sets, we trained a multinomial classifier to distinguish between four lung cancer subtypes. Our meta-analysis-derived classifier included 58 genes and achieved 91% accuracy on leave-one-study-out cross-validation and on three independent data sets. Our method makes meta-analysis of gene expression more systematic and expands the range of questions that a meta-analysis can be used to address. As the amount of publicly available gene expression data continues to grow, our method will be an effective tool to help distill these data into knowledge.

Description: Understanding the contribution of biogenic magnetic particles into sedimentary assemblages is a current challenge in palaeomagnetism. It has been demonstrated recently that magnetic particles produced through biologically controlled mineralization processes, such as magnetosomes from magnetotactic bacteria, contribute to the recording of natural remanent magnetization in marine and lacustrian sediments. Contributions from other, biologically induced, mineralization types, which are known from multiple laboratory experiments to include magnetic minerals, remain largely unknown. Here, we report magnetic properties of iron minerals formed by a community of iron- and manganese-reducing bacteria isolated from a natural groundwater deposit during a 2 yr long incubation experiment. The main iron phases of the biomineralized mass are lepidocrocite, goethite and magnetite, each of which has environmental significance. Unlike the majority of the previous studies that reported superparamagnetic grain size, and thus no remanence carrying capacity of biologically induced magnetite, hysteresis and first-order reversal curves measurements in our study have not detected significant superparamagnetic contribution. The biomineralized mass, instead, contains a mixture of single-domain to pseudo-single-domain and multidomain magnetite particles that are capable of carrying a stable chemical remanent magnetization. Isothermal remanent magnetization acquisition parameters and first-order reversal curves signatures of the biomineralized samples deviate from previously proposed criteria for the distinction of extracellular (biologically induced) magnetic particles in mixtures. Given its potential significance as a carrier of natural remanent magnetization, environmental requirements, distribution in nature and the efficiency in the geomagnetic field recording by biologically induced mineralization need comprehensive investigation.

Description: Advanced sequencing technologies have generated a plethora of data for many chromatin marks in multiple tissues and cell types, yet there is lack of a generalized tool for optimal utility of those data. A major challenge is to quantitatively model the epigenetic dynamics across both the genome and many cell types for understanding their impacts on differential gene regulation and disease. We introduce IDEAS, an i ntegrative and d iscriminative e pigenome a nnotation s ystem, for jointly characterizing epigenetic landscapes in many cell types and detecting differential regulatory regions. A key distinction between our method and existing state-of-the-art algorithms is that IDEAS integrates epigenomes of many cell types simultaneously in a way that preserves the position-dependent and cell type-specific information at fine scales, thereby greatly improving segmentation accuracy and producing comparable annotations across cell types.

Description: Monthly means of the magnetic field measurements at ground observatories are a key data source for studying temporal changes of the core magnetic field. However, when they are calculated in the usual way, contributions of external (magnetospheric and ionospheric) origin may remain, which make them less favourable for studying the field generated by dynamo action in the core. We remove external field predictions, including a new way of characterizing the magnetospheric ring current, from the data and then calculate revised monthly means using robust methods. The geomagnetic secular variation (SV) is calculated as the first annual differences of these monthly means, which also removes the static crustal field. SV time-series based on revised monthly means are much less scattered than those calculated from ordinary monthly means, and their variances and correlations between components are smaller. On the annual to decadal timescale, the SV is generated primarily by advection in the fluid outer core. We demonstrate the utility of the revised monthly means by calculating models of the core surface advective flow between 1997 and 2013 directly from the SV data. One set of models assumes flow that is constant over three months; such models exhibit large and rapid temporal variations. For models of this type, less complex flows achieve the same fit to the SV derived from revised monthly means than those from ordinary monthly means. However, those obtained from ordinary monthly means are able to follow excursions in SV that are likely to be external field contamination rather than core signals. Having established that we can find models that fit the data adequately, we then assess how much temporal variability is required. Previous studies have suggested that the flow is consistent with torsional oscillations (TO), solid body-like oscillations of fluid on concentric cylinders with axes aligned along the Earth's rotation axis. TO have been proposed to explain decadal timescale changes in the length-of-day. We invert for flow models where the only temporal changes are consistent with TO, but such models have an unacceptably large data misfit. However, if we relax the TO constraint to allow a little more temporal variability, we can fit the data as well as with flows assumed constant over three months, demonstrating that rapid SV changes can be reproduced by rather small flow changes. Although the flow itself changes slowly, its time derivative can be locally (temporally and spatially) large, in particular when and where core surface secular acceleration peaks. Spherical harmonic expansion coefficients of the flows are not well resolved, and many of them are strongly correlated. Averaging functions, a measure of our ability to determine the flow at a given location from the data distribution available, are poor approximations to the ideal, even when centred on points of the core surface below areas of high observatory density. Both resolution and averaging functions are noticeably worse for the toroidal flow component, which dominates the flow, than the poloidal flow component, except around the magnetic equator where averaging functions for both components are poor.

Description: The electrical rock conductivity is a sensitive indicator for carbon dioxide (CO 2 ) injection and migration processes. For a reliable balancing of the free CO 2 in pore space with petrophysical models such as Archie's law or for the detection of migrating CO 2 , detailed knowledge of the pore water conductivity during interaction with CO 2 is essential but not available yet. Contrary to common assumptions, pore water conductivity cannot be assumed constant since CO 2 is a reactive gas that dissolves into the pore water in large amounts and provides additional charge carriers due to the dissociation of carbonic acid. We consequently carried out systematic laboratory experiments to quantify and analyse the changes in saline pore water conductivity caused by CO 2 at thermodynamic equilibrium. Electrical conductivity is measured on pore water samples for pressures up to 30 MPa and temperatures up to 80 °C. The parameter range covers the gaseous, liquid and supercritical state of the CO 2 involved. Pore water salinities from 0.006 up to 57.27 g L –1 sodium chloride were investigated as well as selective other ion species. At the same time, the CO 2 concentration in the salt solution was determined by a wet-chemical procedure. A two-regime behaviour appears: for small salinities, we observe an increase of up to more than factor 3 in the electrical pore water conductivity, which strongly depends on the solution salinity (low-salinity regime). This is an expected behaviour, since the additional ions originating from the dissociation of carbonic acid positively contribute to the solution conductivity. However, when increasing salinities are considered this effect is completely diminished. For highly saline solutions, the increased mutual impeding causes the mobility of all ions to decrease, which may result in a significant reduction of conductivity by up to 15 per cent despite the added CO 2 (high-salinity regime). We present the data set covering the pressure, temperature, salinity and ion species dependence of the CO 2 effect. Furthermore, the observations are analysed and predicted with a semi-analytical formulation for the electrical pore water conductivity taking into account the species’ interactions. For the applicability of our results in practice of exploration and monitoring, we additionally provide a purely empirical formulation to compute the impact of CO 2 on pore water conductivity at equilibrium which only requires the input of pressure, temperature and salinity information.

Description: Although the Matuyama–Brunhes boundary (MBB) in the Chinese Loess Plateau (CLP) is very important in reliably correlating Quaternary loess with other sediments in the world, particularly with marine and polar ice cores, its exact stratigraphic position remains controversial. Previous investigations usually placed the MBB between paleosol unit S8 and loess unit L8 in various locations. To better understand the spatial differences in the MBB position, a high-resolution paleomagnetic study was conducted in a loess section of the Lantian Basin at the southern margin of the CLP. The results show that the MBB is situated in the middle of the relatively weak paleosol unit S7, consistent with a recent report on the MBB based on a 10 Be study from the Xifeng and Luochuan loess sections of the central CLP. However, the regional anomalously low magnetic susceptibility in paleosols S7 and S8 indicates that it is more reliable to determine the paleoclimate boundaries between loess and paleosol horizons of this segment with median grain size. Then, the MBB in the Yushan section can be correlated with the bottom of paleosol S7, corresponding to the older part of interglacial marine isotope stage 19. This result temporally reconciles the striking discrepancy of the position of the MBB recorded in between loess and other typical sedimentary sequences, and further confirms that the stratigraphic position of the MBB could spatially vary to a certain extent due to regional sedimentary or paleoclimatic conditions in the marginal areas of the CLP. In the Yushan section, the high-frequency variations of paleomagnetic directions during a long period of ~31 ka before the MBB, however, could not be attributed to a genuine response to the true geomagnetic behaviour. Moreover, the climate offset defined by the magnetic susceptibility and median grain size of the section can be preliminarily attributed to the regional geology and paleoenvironment background. A multiproxy-based stratigraphic division is considered very necessary when paleomagnetic and climatic boundaries are defined exactly in a specific area of the southern CLP.

Description: We study rotating thermal convection in spherical shells as prototype for flow in the cores of terrestrial planets, gas planets or in stars. We base our analysis on a set of about 450 direct numerical simulations of the (magneto)hydrodynamic equations under the Boussinesq approximation. The Ekman number ranges from 10 –3 to 10 –5 . The supercriticality of the convection reaches about 1000 in some models. Four sets of simulations are considered: non-magnetic simulations and dynamo simulations with either free-slip or no-slip flow boundary conditions. The non-magnetic setup with free-slip boundaries generates the strongest zonal flows. Both non-magnetic simulations with no-slip flow boundary conditions and self-consistent dynamos with free-slip boundaries have drastically reduced zonal-flows. Suppression of shear leads to a substantial gain in heat-transfer efficiency, increasing by a factor of 3 in some cases. Such efficiency enhancement occurs as long as the convection is significantly influenced by rotation. At higher convective driving the heat-transfer efficiency tends towards that of the classical non-rotating Rayleigh–Bénard system. Analysis of the latitudinal distribution of heat flow at the outer boundary reveals that the shear is most effective at suppressing heat-transfer in the equatorial regions. Simulations with convection zones of different thickness show that the zonal flows become less energetic in thicker shells, and, therefore, their effect on heat-transfer efficiency decreases. Furthermore, we explore the influence of the magnetic field on the non-zonal flow components of the convection. For this we compare the heat-transfer efficiency of no-slip non-magnetic cases with that of the no-slip dynamo simulations. We find that at E = 10 –5 magnetic field significantly affects the convection and a maximum gain of about 30 per cent (as compared to the non-magnetic case) in heat-transfer efficiency is obtained for an Elsasser number of about 3. Our analysis motivates us to speculate that convection in the polar regions in dynamos at E = 10 –5 is probably in a ‘magnetostrophic’ regime.

Description: ‘Anomalous’ magnetic fabrics in dikes that appear to indicate flow into the wall confound many workers. Here, we present extensive magnetic data on five dikes from Tenerife, Canary Islands, and use these to interpret the causes of the anomalous fabrics. Comparison of the anisotropy of magnetic susceptibility (AMS) and anhysteretic magnetization (AARM) results show that, in some cases, the anomalous fabrics are caused by single-domain grains, which produce AMS fabrics perpendicular to the grain elongation, whereas AARM fabrics are parallel. To check this, hysteresis experiments were used to characterize the domain state. These show most are mixtures of pseudo-single-domain or single-domain plus multi-domain particles, but many have wasp-waisted hysteresis loops, likely indicating mixed populations of stable single-domain and superparamagnetic grains. First-order reversal curves were used to better characterize this and show mixtures of stable single-domain and superparamagnetic grains dominate the magnetic signal. Magnetic particles at the stable single-domain/superparamagnetic threshold are unstable at timespans relevant to the analytical techniques, so they produce complicated results. This suggests that anomalous AMS fabrics in dikes cannot simply be attributed to elongated stable single-domain particles and that mixtures of the different grain types can produce hybrid fabrics, in which the fabrics are neither perpendicular or parallel to the dike plane, that are difficult to interpret without extensive magnetic analysis.

Description: We derive scaling relationships for planetary dynamos based on a balance between energy production and Joule dissipation, and between the curl of the buoyancy and Coriolis forces. These scaling relationships are deduced for the particular case of dynamos driven by helical waves, but are shown to have a much broader applicability. They are consistent with the evidence of the numerical dynamos, yielding predictions consistent with published empirical scaling laws and also with the observed transition from dipolar to multipolar dynamos. A direct comparison with the observational evidence for the planets is hampered by the fact that we do not know what sets the smallest scale of the motion in the planets. Nevertheless, we use our scaling relationships to show that the traditional assumption that the Elsasser number is of order unity is inconsistent with the observation that the gas-giant dynamos are dipolar dynamos, as is the more recent suggestion that the strength of the dipole is independent of rotation rate and controlled by the buoyancy flux alone. On the other hand, we show that the observational data is consistent with the hypothesis that a dipolar dynamo saturates at the lowest permissible magnetic energy compatible with a given buoyancy flux.

Description: A majority of large-scale bacterial genome rearrangements involve mobile genetic elements such as insertion sequence (IS) elements. Here we report novel insertions and excisions of IS elements and recombination between homologous IS elements identified in a large collection of Escherichia coli mutation accumulation lines by analysis of whole genome shotgun sequencing data. Based on 857 identified events (758 IS insertions, 98 recombinations and 1 excision), we estimate that the rate of IS insertion is 3.5 x 10 –4 insertions per genome per generation and the rate of IS homologous recombination is 4.5 x 10 –5 recombinations per genome per generation. These events are mostly contributed by the IS elements IS 1 , IS 2 , IS 5 and IS 186 . Spatial analysis of new insertions suggest that transposition is biased to proximal insertions, and the length spectrum of IS-caused deletions is largely explained by local hopping. For any of the ISs studied there is no region of the circular genome that is favored or disfavored for new insertions but there are notable hotspots for deletions. Some elements have preferences for non-coding sequence or for the beginning and end of coding regions, largely explained by target site motifs. Interestingly, transposition and deletion rates remain constant across the wild-type and 12 mutant E. coli lines, each deficient in a distinct DNA repair pathway. Finally, we characterized the target sites of four IS families, confirming previous results and characterizing a highly specific pattern at IS 186 target-sites, 5'-GGGG(N6/N7)CCCC-3'. We also detected 48 long deletions not involving IS elements.

Description: We present a detailed palaeomagnetic study from 35 sites on Holocene lava flows of the Tongariro Volcanic Centre, central North Island, New Zealand. Prior to the study the eruption ages of these flows were constrained to within a few thousand years by recently published high-precision 40 Ar/ 39 Ar geochronological data and tephrostratigraphic controls. Correlation of flow mean palaeomagnetic directions with a recently published continuous sediment record from Lake Mavora, Fiordland, allows us to reduce the age uncertainty to 300–500 yr in some cases. Our refined ages significantly improve the chronology of Holocene effusive eruptions of the volcanoes of the Tongariro Volcanic Centre. For instance, differences in the palaeomagnetic directions recorded by lavas from the voluminous Iwikau and Rangataua members suggest that individual effusive periods lasted up to thousands of years and that these bursts have been irregularly spaced over time. While over the last few millennia the effusive eruptive activity from Mt Ruapehu has been relatively quiet, the very young age (200–500 BP) of a Red Crater sourced flow suggests that effusive activity around Mt Tongariro lasted into the past few centuries. This adds an important hazard context to the historical record, which has otherwise comprised frequent relatively small, tephra producing, explosive eruptions without the production of lava flows.

Description: The inversion of magnetotelluric data into subsurface electrical conductivity poses an ill-posed problem. Smoothing constraints are widely employed to estimate a regularized solution. Here, we present an alternative inversion scheme that estimates a sparse representation of the model in a wavelet basis. The objective of the inversion is to determine the few non-zero wavelet coefficients which are required to fit the data. This approach falls into the class of sparsity constrained inversion schemes and minimizes the combination of the data misfit in a least-squares 2 sense and of a model coefficient norm in an 1 sense ( 2 - 1 minimization). The 1 coefficient norm renders the solution sparse in a suitable representation such as the multiresolution wavelet basis, but does not impose explicit structural penalties on the model as it is the case for 2 regularization. The presented numerical algorithm solves the mixed 2 - 1 norm minimization problem for the nonlinear magnetotelluric inverse problem. We demonstrate the feasibility of our algorithm on synthetic 2-D MT data as well as on a real data example. We found that sparse models can be estimated by inversion and that the spatial distribution of non-vanishing coefficients indicates regions in the model which are resolved.

Description: We examined the mathematical and computational aspects of the magnetic potential, vector and gradient tensor fields of a tesseroid in a geocentric spherical coordinate system (SCS). This work is relevant for 3-D modelling that is performed with lithospheric vertical scales and global, continent or large regional horizontal scales. The curvature of the Earth is significant at these scales and hence, a SCS is more appropriate than the usual Cartesian coordinate system (CCS). The 3-D arrays of spherical prisms (SP; ‘tesseroids’) can be used to model the response of volumes with variable magnetic properties. Analytical solutions do not exist for these model elements and numerical or mixed numerical and analytical solutions must be employed. We compared various methods for calculating the response in terms of accuracy and computational efficiency. The methods were (1) the spherical coordinate magnetic dipole method (MD), (2) variants of the 3-D Gauss–Legendre quadrature integration method (3-D GLQI) with (i) different numbers of nodes in each of the three directions, and (ii) models where we subdivided each SP into a number of smaller tesseroid volume elements, (3) a procedure that we term revised Gauss–Legendre quadrature integration (3-D RGLQI) where the magnetization direction which is constant in a SCS is assumed to be constant in a CCS and equal to the direction at the geometric centre of each tesseroid, (4) the Taylor's series expansion method (TSE) and (5) the rectangular prism method (RP). In any realistic application, both the accuracy and the computational efficiency factors must be considered to determine the optimum approach to employ. In all instances, accuracy improves with increasing distance from the source. It is higher in the percentage terms for potential than the vector or tensor response. The tensor errors are the largest, but they decrease more quickly with distance from the source. In our comparisons of relative computational efficiency, we found that the magnetic potential takes less time to compute than the vector response, which in turn takes less time to compute than the tensor gradient response. The MD method takes less time to compute than either the TSE or RP methods. The efficiency of the (GLQI and) RGLQI methods depends on the polynomial order, but the response typically takes longer to compute than it does for the other methods. The optimum method is a complex function of the desired accuracy, the size of the volume elements, the element latitude and the distance between the source and the observation. For a model of global extent with typical model element size (e.g. 1 degree horizontally and 10 km radially) and observations at altitudes of 10s to 100s of km, a mixture of methods based on the horizontal separation of the source and observation separation would be the optimum approach. To demonstrate the RGLQI method described within this paper, we applied it to the computation of the response for a global magnetization model for observations at 300 and 30 km altitude.

Description: Extensive and multi-dimensional data sets generated from recent cancer omics profiling projects have presented new challenges and opportunities for unraveling the complexity of cancer genome landscapes. In particular, distinguishing the unique complement of genes that drive tumorigenesis in each patient from a sea of passenger mutations is necessary for translating the full benefit of cancer genome sequencing into the clinic. We address this need by presenting a data integration framework (OncoIMPACT) to nominate patient-specific driver genes based on their phenotypic impact. Extensive in silico and in vitro validation helped establish OncoIMPACT's robustness, improved precision over competing approaches and verifiable patient and cell line specific predictions (2/2 and 6/7 true positives and negatives, respectively). In particular, we computationally predicted and experimentally validated the gene TRIM24 as a putative novel amplified driver in a melanoma patient. Applying OncoIMPACT to more than 1000 tumor samples, we generated patient-specific driver gene lists in five different cancer types to identify modes of synergistic action. We also provide the first demonstration that computationally derived driver mutation signatures can be overall superior to single gene and gene expression based signatures in enabling patient stratification and prognostication. Source code and executables for OncoIMPACT are freely available from http://sourceforge.net/projects/oncoimpact .

Description: Kamacite particles (Fe–Ni, Ni 〈 5 per cent), are very common in extra-terrestrial materials, such as meteorites. It is normally assumed that for kamacite particles to be reliable recorders of magnetic fields, they need to be magnetically uniform (single domain, SD) and thermally stable. Larger particles subdivide into non-uniform multidomain (MD) magnetic structures that produce weaker magnetic signals, while small SD particles become magnetically unstable due to thermal fluctuations and exhibit superparamagnetic behaviour. In this paper we determine the first micromagnetic calculation of the stable SD range domain-state phase diagram for metallic iron; previous calculations were analytical. There is a significant increase in the critical size for the SD/MD threshold size, for example, for cube-shaped iron particles, the critical SD/MD threshold has now been estimated to be 25 nm, compared to 17 nm for previous estimates. The larger critical SD/MD threshold size for iron, agrees better with previously published nanometric observations of domain state for FeNi particles, then early analytical models.

Description: Analysis of rewired upstream subnetworks impacting downstream differential gene expression aids the delineation of evolving molecular mechanisms. Cumulative statistics based on conventional differential correlation are limited for subnetwork rewiring analysis since rewiring is not necessarily equivalent to change in correlation coefficients. Here we present a computational method ChiNet to quantify subnetwork rewiring by statistical heterogeneity that enables detection of potential genotype changes causing altered transcription regulation in evolving organisms. Given a differentially expressed downstream gene set, ChiNet backtracks a rewired upstream subnetwork from a super-network including gene interactions known to occur under various molecular contexts. We benchmarked ChiNet for its high accuracy in distinguishing rewired artificial subnetworks, in silico yeast transcription-metabolic subnetworks, and rewired transcription subnetworks for Candida albicans versus Saccharomyces cerevisiae , against two differential-correlation based subnetwork rewiring approaches. Then, using transcriptome data from tolerant S. cerevisiae strain NRRL Y-50049 and a wild-type intolerant strain, ChiNet identified 44 metabolic pathways affected by rewired transcription subnetworks anchored to major adaptively activated transcription factor genes YAP1 , RPN4 , SFP1 and ROX1 , in response to toxic chemical challenges involved in lignocellulose-to-biofuels conversion. These findings support the use of ChiNet in rewiring analysis of subnetworks where differential interaction patterns resulting from divergent nonlinear dynamics abound.

Description: This paper presents a distributed magnetotelluric inversion scheme based on adaptive finite-element method (FEM). The key novel aspect of the introduced algorithm is the use of automatic mesh refinement techniques for both forward and inverse modelling. These techniques alleviate tedious and subjective procedure of choosing a suitable model parametrization. To avoid overparametrization, meshes for forward and inverse problems were decoupled. For calculation of accurate electromagnetic (EM) responses, automatic mesh refinement algorithm based on a goal-oriented error estimator has been adopted. For further efficiency gain, EM fields for each frequency were calculated using independent meshes in order to account for substantially different spatial behaviour of the fields over a wide range of frequencies. An automatic approach for efficient initial mesh design in inverse problems based on linearized model resolution matrix was developed. To make this algorithm suitable for large-scale problems, it was proposed to use a low-rank approximation of the linearized model resolution matrix. In order to fill a gap between initial and true model complexities and resolve emerging 3-D structures better, an algorithm for adaptive inverse mesh refinement was derived. Within this algorithm, spatial variations of the imaged parameter are calculated and mesh is refined in the neighborhoods of points with the largest variations. A series of numerical tests were performed to demonstrate the utility of the presented algorithms. Adaptive mesh refinement based on the model resolution estimates provides an efficient tool to derive initial meshes which account for arbitrary survey layouts, data types, frequency content and measurement uncertainties. Furthermore, the algorithm is capable to deliver meshes suitable to resolve features on multiple scales while keeping number of unknowns low. However, such meshes exhibit dependency on an initial model guess. Additionally, it is demonstrated that the adaptive mesh refinement can be particularly efficient in resolving complex shapes. The implemented inversion scheme was able to resolve a hemisphere object with sufficient resolution starting from a coarse discretization and refining mesh adaptively in a fully automatic process. The code is able to harness the computational power of modern distributed platforms and is shown to work with models consisting of millions of degrees of freedom. Significant computational savings were achieved by using locally refined decoupled meshes.

Description: A spatially averaged mean-field model for fully or partially ordered members of the ilmenite–hematite solid solution series is rigorously derived from the Heisenberg Hamiltonian by first assuming no temporal correlation of atomic spins, and then by spatially averaging over spins at equivalent atomic positions. The model is based on the geometry of exchange interactions between nearest and next-nearest neighbours and predicts magnetization curves in homogenous solid solutions with variable degree of order. While the general framework presented can also be applied to atomic scale models, and to other solid solution series, here the symmetries of the ilmenite–hematite lattice are exploited to show that four different sublattice magnetizations and six independent combinations of exchange constants determine the temperature variation of the magnetization curves. Comparing measured Curie temperatures T C and M s ( T ) curves to model predictions results in accurate constraints for these combinations. It is also possible to calculate predictions for high-field magnetization slopes HF , which not only improve accurate experimental determination of the Curie temperature but also provide a new magnetic method to estimate the order parameter for ilmenite–hematite solid solution samples.

Description: Ribosome biogenesis, a central and essential cellular process, occurs through sequential association and mutual co-folding of protein–RNA constituents in a well-defined assembly pathway. Here, we construct a network of co-evolving nucleotide/amino acid residues within the ribosome and demonstrate that assembly constraints are strong predictors of co-evolutionary patterns. Predictors of co-evolution include a wide spectrum of structural reconstitution events, such as cooperativity phenomenon, protein-induced rRNA reconstitutions, molecular packing of different rRNA domains, protein–rRNA recognition, etc. A correlation between folding rate of small globular proteins and their topological features is known. We have introduced an analogous topological characteristic for co-evolutionary network of ribosome, which allows us to differentiate between rRNA regions subjected to rapid reconstitutions from those hindered by kinetic traps. Furthermore, co-evolutionary patterns provide a biological basis for deleterious mutation sites and further allow prediction of potential antibiotic targeting sites. Understanding assembly pathways of multicomponent macromolecules remains a key challenge in biophysics. Our study provides a ‘proof of concept’ that directly relates co-evolution to biophysical interactions during multicomponent assembly and suggests predictive power to identify candidates for critical functional interactions as well as for assembly-blocking antibiotic target sites.

Description: Tibetan Plateau uplift and global cooling have generally been thought to have caused the drying of the Asian inland, but how and when these factors drove the aridification is unknown. The Xining Basin at the NE Tibetan Plateau received continuous Eocene–Miocene fine-grained sediments, providing an excellent opportunity to address this question. Here we present detailed rock magnetic and diffuse reflectance spectroscopy (DRS) analyses for a well-dated Cenozoic sedimentary record from the Xiejia section in the basin. Magnetic susceptibility ( ), saturation magnetization (Ms) and saturation isothermal remanent magnetization (SIRM) in this section show a long-term decreasing trend from ~52 to ~25 Ma, well coinciding with global cooling and drying in the region, and an increasing trend since ~25 Ma, which is in contrast to the further progressing aridification of the basin. Thermomagnetic results and DRS-determined hematite contents suggest that the relative content of magnetite and hematite is the main control on the , Ms and SIRM values. We argue that the long-term Eocene–Oligocene global cooling increased the drying of the Asian inland, lowering the lake level and exposing a larger area to low temperature oxidation for longer times, thus producing more hematite and leading to the decreasing trend of magnetic concentration parameters from ~52 to ~25 Ma. An intensive uplift of the NE Tibetan Plateau since ~25 Ma, associated with a change in the sedimentary source, might be responsible for the increase of , Ms and SIRM after 25 Ma.

Description: The emergence of new sequencing technologies has facilitated the use of bacterial whole genome alignments for evolutionary studies and outbreak analyses. These datasets, of increasing size, often include examples of multiple different mechanisms of horizontal sequence transfer resulting in substantial alterations to prokaryotic chromosomes. The impact of these processes demands rapid and flexible approaches able to account for recombination when reconstructing isolates’ recent diversification. Gubbins is an iterative algorithm that uses spatial scanning statistics to identify loci containing elevated densities of base substitutions suggestive of horizontal sequence transfer while concurrently constructing a maximum likelihood phylogeny based on the putative point mutations outside these regions of high sequence diversity. Simulations demonstrate the algorithm generates highly accurate reconstructions under realistically parameterized models of bacterial evolution, and achieves convergence in only a few hours on alignments of hundreds of bacterial genome sequences. Gubbins is appropriate for reconstructing the recent evolutionary history of a variety of haploid genotype alignments, as it makes no assumptions about the underlying mechanism of recombination. The software is freely available for download at github.com/sanger-pathogens/Gubbins , implemented in Python and C and supported on Linux and Mac OS X.

Description: Genomic structural variation (SV), a common hallmark of cancer, has important predictive and therapeutic implications. However, accurately detecting SV using high-throughput sequencing data remains challenging, especially for ‘targeted’ resequencing efforts. This is critically important in the clinical setting where targeted resequencing is frequently being applied to rapidly assess clinically actionable mutations in tumor biopsies in a cost-effective manner. We present BreaKmer, a novel approach that uses a ‘kmer’ strategy to assemble misaligned sequence reads for predicting insertions, deletions, inversions, tandem duplications and translocations at base-pair resolution in targeted resequencing data. Variants are predicted by realigning an assembled consensus sequence created from sequence reads that were abnormally aligned to the reference genome. Using targeted resequencing data from tumor specimens with orthogonally validated SV, non-tumor samples and whole-genome sequencing data, BreaKmer had a 97.4% overall sensitivity for known events and predicted 17 positively validated, novel variants. Relative to four publically available algorithms, BreaKmer detected SV with increased sensitivity and limited calls in non-tumor samples, key features for variant analysis of tumor specimens in both the clinical and research settings.

Description: The Mt Amram igneous complex (AIC) represents northern tip of the Neoproterozoic Arabian Nubian Shield (ANS). For the first time the AIC deep structure was studied using the gravity, aero and ground magnetic, magnetic susceptibility and density measurements and geological data. Analysing all available data at the Amram area we concluded what only monzonite body can be reason for gravity high and coinciding reduced to pole (RTP) maximum. Geological knowledge allowed suggesting its intrusive character and compact body form. Cluster of inverse solutions (Werner deconvolution) localized this body as initial model for forward modelling. Further iterations (2 3/4 -D forward modelling) clarified the monzonite geometry and properties; the modelling allowed also to investigate the non-uniqueness and estimate also the confident intervals for final solution. The research consists three interconnected stages. At the detailed scale, ground magnetic data suggested three magmatic blocks of few hundred meters shifted dextral about 100 m along the Zefunut fault. Estimated accuracy for geometry of the magnetic bodies is a few tens metres. At the middle scale, quantitative gravity and magnetic interpretations provide model of the monzonite body, which is an order of magnitude more than the volume of the felsic rhyolites and granite rocks. Boundary of the whole monzonite body was estimated with accuracy as a hundred meters. As a result we suggest that the parent magma for the AIC is the monzonite, similar to the model suggested for the Timna Igneous Complex 12 km north of the AIC. The model developed can be applied to evaluate the subsurface volumes of the mafic magmatic rocks in adjacent locations. At the regional scale for exposed the Sinai and Arab Saudi Precambrian crystalline shield our approach allows to understand the apparent contradiction between geological predominantly granite composition (low magnetic rocks) and magnetic data. The aeromagnetic data show number strong magnetic anomalies suggesting the presence large volume of high magnetic (mainly basic) rocks at the depth. This problem is proposed for future research.

Description: 3-D electromagnetic (EM) studies of the Earth have advanced significantly over the past decade. Despite a certain success of the 3-D EM inversions of real data sets, the quantitative assessment of the recovered models is still a challenging problem. It is known that one can gain valuable information about model uncertainties from the analysis of Hessian matrix. However, even with modern computational capabilities the calculation of the Hessian matrix based on numerical differentiation is extremely time consuming. Much more efficient way to compute the Hessian matrix is provided by an ‘adjoint sources’ methodology. The computation of Hessian matrix (and Hessian-vector products) using adjoint formulation is now well-established approach, especially in seismic inverse modelling. As for EM inverse modelling we did not find in the literature a description of the approach, which would allow EM researchers to apply this methodology in a straightforward manner to their scenario of interest. In the paper, we present formalism for the efficient calculation of the Hessian matrix using adjoint sources approach. We also show how this technique can be implemented to calculate multiple Hessian-vector products very efficiently. The formalism is general in the sense that it allows to work with responses that arise in EM problem set-ups either with natural- or controlled-source excitations. The formalism allows for various types of parametrization of the 3-D conductivity distribution. Using this methodology one can readily obtain appropriate formulae for the specific sounding methods. To illustrate the concept we provide such formulae for two EM techniques: magnetotellurics and controlled-source sounding with vertical magnetic dipole as a source.

Description: Geomagnetic jerks are sudden changes of trend in the geomagnetic secular variation. The Earth's mantle behaves as a filter for the jerks, causing a delayed and a smoothed signal at the Earth's surface. Backus' mantle filter theory relies on approximating the impulse response function (IRF) of the mantle by a Gaussian. The advantage of this theory is the linear relation between jerks’ delay times and the mantle electrical conductivity, as expressed by kernels. However, the limitations of this theory arise when negative delay and/or smoothing times occur. The applicability of the mantle filter theory is examined by analysing the validity of the Gaussian as an approximation for the composite IRF (CIRF) at a given location. We show that the electrical conductivity of the lower mantle is mostly responsible for the jerk delay time. Alternating sign CIRFs might cause negative delay and/or smoothing times which prevents the use of the mantle filter theory. Adequate/inadequate Gaussian approximations to the CIRFs give small/large differences in the convolved jerk occurrence times. Most observatories yield positive time constants, but in most cases the difference in the jerk occurrence times exceeds 0.5 yr.

Description: A new method of estimating fluid flows near the core–mantle boundary is presented. Both the geomagnetic diffusion and dynamics in a boundary layer are taken into account under the no-slip condition at the core–mantle boundary. Below the boundary layer, the tangentially geostrophic constraint is imposed on the flow to be estimated. The radial component of the geomagnetic field just within the core is calculated from the radial component and its partial derivatives with respect to the radius, first, derived from continuity of the solenoidal field across the finite jump in electrical conductivity and, second, the geomagnetic diffusion. Below the boundary layer, fluid motions thus derived show some features common to those estimated on the basis of the frozen-flux approximation, and their accompanying vortices at mid and high latitudes are well correlated with upwellings or downwellings inside the boundary layer. The positional relation of upwellings and downwellings in equatorial regions inside and below the boundary layer suggests existence of columnar convective cells there. Intense magnetic flux spots seen in equatorial regions are likely to be produced not by downwellings at tops of cyclonic columns but by flux expulsion due to upwellings between adjacent columnar flows. The present new method leads to a new insight on core surface flows that have not been noted so far.

Description: In most numerical simulations of the Earth's core the dynamo resides outside the tangent cylinder and may be crudely classified as being of the α 2 type. In this region the flow comprises a sea of thin columnar vortices aligned with the rotation axis, taking the form of alternating cyclones and anticyclones. The dynamo is thought to be driven by these columnar vortices within which the flow is observed to be highly helical, helicity being a crucial ingredient of planetary dynamos. As noted in Davidson, one of the mysteries of this dynamo cartoon is the origin of the helicity, which is observed to be positive in the south and negative in the north. While Ekman pumping at the mantle can induce helicity in some of the overly viscous numerical simulations, it is extremely unlikely to be a significant source within planets. In this paper we return to the suggestion of Davidson that the helicity observed in the less viscous simulations owes its existence to helical wave packets, launched in and around the equatorial plane where the buoyancy flux is observed to be strong. Here we show that such wave packets act as a potent source of planetary helicity, constituting a simple, robust mechanism that yields the correct sign for h north and south of the equator. Since such a mechanism does not rely on the presence of a mantle, it can operate within both the Earth and the gas giants. Moreover, our numerical simulations show that helical wave packets dispersing from the equator produce a random sea of thin, columnar cyclone/anticyclone pairs, very like those observed in the more strongly forced dynamo simulations. We examine the local dynamics of helical wave packets dispersing from the equatorial regions, as well as the overall nature of an α 2 -dynamo driven by such wave packets. Our local analysis predicts the mean emf induced by helical waves, an analysis that rests on a number of simple approximations which are consistent with our numerical experiments, while our global analysis yields exact integral relationships between the mean emf induced by the wave packets and the large-scale dipole and azimuthal field. Combining these local and integral equations yields a kinematic model for an α 2 -dynamo driven by helical waves. Order-of-magnitude estimates based on these equations suggest that such a dynamo is indeed feasible.

Description: The Lower Danube Basin is one of the most important loess regions from Europe, which have provided excellent archives for long-term high-resolution palaeoclimate studies. The aim of this paper is to derive new information on the Middle–Late Pleistocene palaeoenvironment from a high resolution multiproxy assessment of the iron mineralogical composition at the Costinesti loess-palaeosol sequence located on the western Black Sea shore. It is the easternmost loess section in the Romanian loess region studied and its distinct pattern of the proxy records can be used to correlate the lower Danube loess to other key sites of the Moldavia and Ukraine loess regions. To investigate the climatic control on soft and hard ferromagnetic minerals we used several types of rock magnetic properties: magnetic susceptibility and its frequency dependence, anhysteretic remanent magnetization, isothermal remanent magnetization, hysteresis properties and FORC distributions, an unmixing model for isothermal remanent magnetization curves and high field (up to 8 T) isothermal remanence measurements. Our results show that the palaeosol horizons, formed during interglacials and climatically more favored periods of the Pleistocene, experienced pedogenic alteration, resulting in high amounts of superparamagnetic, single domain and pseudosingle domain magnetite/maghemite grains and hematite. The loess layers, formed during glacial periods, are mainly dominated by multidomain and/or pseudosingle domain oxidized magnetite and some hematite, all probably of aeolian origin. Goethite contribution is probably minor and constant both in loess and palaeosol horizons. We review the correlation of the loess sections from the lower Danube basin concluding that the new results can be interpreted as a support for the transition of a Mediterranean type climate to a steppe type climate in the last two interglacial periods in the western Black Sea. Because the pattern of magnetic susceptibility data from the lower Danube basin is changing relative fast with distance from the Black Sea shore, it probably reflects the local influence of the Black Sea on continental scale climatic oscillations during the last 600 ka. The values of background magnetic susceptibility of the Romanian loess-palaeosol sections indicate that the main source area of the dust changed during this climatic transition. Our analysis also shows that the age of the loess-palaeosol sections from the Eastern European low lands (Moldavia and Ukraine) must be revised to be in agreement with the chronostratigraphy of the sections from the Lower Danube Basin loess area.

Description: Statistical models of geomagnetic secular variations in terms of a Giant Gaussian Process (GGP) are compared with palaeontensity databases (Brunhes epoch, effusive rocks). We tested two GGP models TK and QC (suggested by Tauxe-Kent and Quidelleur-Courtillot) against 392 volcanic palaeointensity data, and found both models to be incompatible with the data. Even though at present there is no accepted GGP model of the Brunhes chron, we found interesting to clarify the details of this negative test. While such an incompatibility is clear enough for the TK model, because of its too small dipole coefficient, the incompatibility of the QC model was due the 12 per cent of the data with very small intensities; remaining 88 per cent of data is fully compatible with QC model QC. The troublesome, low-palaeointensity records belong to different sites and ages so their excess presence is not an effect of particular location or time. We found no evidence that the marked discrepancy between the empirical data and theoretical calculations on the QC model can be explained by the possibility that the low-intensity values belong to known excursions of the geomagnetic field. As an alternative, we conjecture that selected low-intensity values can possibly stem from acquisition of chemical remanent magnetization instead of the thermoremanence.

Description: A new method of estimating fluid flows near the core–mantle boundary is presented. Both the geomagnetic diffusion and dynamics in a boundary layer are taken into account under the no-slip condition at the core–mantle boundary. Below the boundary layer, the tangentially geostrophic constraint is imposed on the flow to be estimated. The radial component of the geomagnetic field just within the core is calculated from the radial component and its partial derivatives with respect to the radius, first, derived from continuity of the solenoidal field across the finite jump in electrical conductivity and, second, the geomagnetic diffusion. Below the boundary layer, fluid motions thus derived show some features common to those estimated on the basis of the frozen-flux approximation, and their accompanying vortices at mid and high latitudes are well correlated with upwellings or downwellings inside the boundary layer. The positional relation of upwellings and downwellings in equatorial regions inside and below the boundary layer suggests existence of columnar convective cells there. Intense magnetic flux spots seen in equatorial regions are likely to be produced not by downwellings at tops of cyclonic columns but by flux expulsion due to upwellings between adjacent columnar flows. The present new method leads to a new insight on core surface flows that have not been noted so far.

Description: Practical applications of controlled-source electromagnetic (EM) modelling require solutions for multiple sources at several frequencies, thus leading to a dramatic increase of the computational cost. In this paper, we present an approach using block Krylov subspace solvers that are iterative methods especially designed for problems with multiple right-hand sides (RHS). Their main advantage is the shared subspace for approximate solutions, hence, these methods are expected to converge in less iterations than the corresponding standard solver applied to each linear system. Block solvers also share the same preconditioner, which is constructed only once. Simultaneously computed block operations have better utilization of cache due to the less frequent access to the system matrix. In this paper, we implement two different block solvers for sparse matrices resulting from the finite-difference and the finite-element discretizations, discuss the computational cost of the algorithms and study their dependence on the number of RHS given at once. The effectiveness of the proposed methods is demonstrated on two EM survey scenarios, including a large marine model. As the results of the simulations show, when a powerful preconditioning is employed, block methods are faster than standard iterative techniques in terms of both iterations and time.

Description: Loess deposits in the arid Central Asia contain valuable information on the evolution of local aridification and dust sources in the Northern Hemisphere. Xinjiang is located in the eastern part of Central Asia and previous researches have revealed the complex enhancement of magnetic susceptibility in loess–paleosol sequences. However, systematic magnetic archives of loess deposit in this arid Asian interior are still far from adequate. In this study, magnetic parameters combined with nonmagnetic properties (granulometry and chromaticity) were analysed on a loess section in Shawan (SW), northwestern China. The section shares a similar magnetic composition with those in the Chinese Loess Plateau (CLP) as well as other sites in Xinjiang. Ferrimagnetic components (magnetite and maghemite) dominate the magnetic signal while the contribution of antiferromagnetic phases (like hematite and goethite) and paramagnetic portions are relatively low. There is no specific correlation between magnetic concentration and pedogenic intensity in the SW section. In general, magnetic enhancement was largely influenced by the paleowind intensity. However, a positive correlation between magnetic susceptibility and pedogenesis is observed in the upper part (0–3.5 m depths), which is characterized by a moderate wind intensity. Moreover, pedogenesis might be responsible for the enhancement of fine magnetic particles in paleosols. Magnetic properties are controlled by coarse magnetic particles in the pseudo-single domain state, but a coarse stable single domain phase was found in certain paleosol samples. The input of detrital fractions from a nearby dust source probably controlled the magnetic properties while a superparamagnetic fraction, which has been deemed as a product of pedogenesis in the CLP, is limited in the SW section. Caution is needed to employ magnetic susceptibility directly for paleoclimatic assessment because of its uncertainty in the Xinjiang loess. However, the ARM /SIRM ratio enhanced in paleosols and has more potential as an alternative index for the weakly pedogenic degree in this area.

Description: Mantle control on the geodynamo is often simulated using numerical dynamos with imposed outer boundary heat flux inferred from lower mantle tomography, assuming that seismic and thermal anomalies in the lowermost mantle are highly correlated. However, non-thermal effects might perturb this idealized linear seismic–thermal mapping. Here we use a probabilistic tomography model to isolate the thermal part of the seismic anomaly in order to impose a more realistic core–mantle boundary heat flux pattern on the outer boundary of numerical dynamo simulations. We demonstrate that on time average these dynamo models have more low-latitude convective and magnetic activity than corresponding models with conventional tomographic heat flux. In addition, the low-latitude magnetic flux and kinetic energy contributions are more time-dependent in the dynamo models with a probabilistic tomography heat flux, and thus may recover the observed latitudinal distribution of geomagnetic flux on the core–mantle boundary, which we propose as a morphological criterion for Earth-like dynamo models.

Description: The statistics of directional data on a sphere can be modelled either using the Fisher distribution that is conditioned on the magnitude being unity, in which case the sample space is confined to the unit sphere, or using the latitude–longitude marginal distribution derived from a trivariate Gaussian model that places no constraint on the magnitude. These two distributions are derived from first principles and compared. The Fisher distribution more closely approximates the uniform distribution on a sphere for a given small value of the concentration parameter, while the latitude–longitude marginal distribution is always slightly larger than the Fisher distribution at small off-axis angles for large values of the concentration parameter. Asymptotic analysis shows that the two distributions only become equivalent in the limit of large concentration parameter and very small off-axis angle.

Description: With advancing computational resources, 3-D inversion techniques have become feasible in recent years and are now a more widely used tool for magnetotelluric (MT) data interpretation. Galvanic distortion caused by small-scale near-surface inhomogeneities remains an obstacle for 3-D MT inversion which so far has experienced little attention. If not considered properly, the effect on 3-D inversion can be immense and result in erroneous subsurface models and interpretations. To tackle the problem we implemented inversion of the distortion-free phase tensor into the ModEM inversion package. The dimensionless phase tensor components describe only variations of the conductivity structure. When inverting these data, particular care has to be taken of the conductivity structure in the a priori model, which provides the reference frame when transferring the information from phase tensors into absolute conductivity values. Our results obtained with synthetic data show that phase tensor inversion can recover the regional conductivity structure in presence of galvanic distortion if the a priori model provides a reasonable assumption for the regional resistivity average. Joint inversion of phase tensor data and vertical magnetic transfer functions improves recovery of the absolute resistivity structure and is less dependent on the prior model. We also used phase tensor inversion for a data set of more than 250 MT sites from the central San Andreas fault, California, where a number of sites showed significant galvanic distortion. We find the regional structure of the phase tensor inversion results compatible with previously obtained models from impedance inversion. In the vicinity of distorted sites, phase tensor inversion models exhibit more homogeneous/smoother conductivity structures.

Description: Earth's liquid core hosts a diverse set of waves with periods ranging from days to thousands of years. One class of waves with periods of several decades is known to arise from an interplay between magnetic, Archimedes and Coriolis forces. These so-called MAC waves are thought to be relevant for interpreting historical fluctuations in the geomagnetic field. In this study, we show that MAC waves provide a good description of time-dependent zonal flow at the top of the core. The same collection of waves also offers a simple explanation for observed fluctuations in the dipole field. Both of these predictions require a stratified layer at the top of the core with a thickness of 130–140 km and a buoyancy frequency comparable to Earth's rotation rate. We extend these predictions to include changes in the length of day (LOD) and find that MAC waves can account for about half of the observed fluctuation at decadal periods. Larger fluctuations are possible when electromagnetic stresses couple MAC waves to flow in the interior of the core. In fact, an idealized model for the coupled motion overestimates the LOD fluctuations, probably reflecting limitations in this idealized model. Our results offer support for stable stratification at the top of the core and suggest a common origin for decadal fluctuations in the geomagnetic field and the LOD.

Description: A detailed magnetic analysis has been conducted on the borehole CJ-1 (172.3 m in length) from the Changjiang (Yangtze) River delta. Results show that the dominant magnetic carriers are magnetite and hematite. Palaeomagnetic results reveal that high-frequency changes in palaeomagnetic inclinations are tied significantly to abnormal anisotropy of magnetic susceptibility (AMS) patterns due to effects of the high-energy depositional environment in this region. On the basis of AMS patterns, doubtful palaeomagnetic directional anomalies can be distinguished from authentic palaeomagnetic excursions. Magnetostratigraphic results indicate that the Matuyama–Brunhes boundary (MBB) was recorded at a depth of ~152.5 m. The presence of several short-lived inclination anomalies implied that the sedimentation could be continuous even at the millennial timescale at certain depth intervals bracketing these fast geomagnetic events. In summary, our study provides new insights into constructing reliable magnetostratigraphy in the delta region. Moreover, our new magnetostratigraphy of the Changjiang River delta deposits will facilitate studies on the relevant long-term palaeoenvironmental evolution of the delta.

Description: Non-stoichiometric monoclinic 4C pyrrhotite (Fe 7 S 8 ) is a major magnetic remanence carrier in the Earth's crust and in extraterrestrial materials. Because of its low-temperature magnetic transition around 30 K also known as Besnus transition, which is considered to be an intrinsic property, this mineral phase is easily detectable in natural samples. Although the physical properties of pyrrhotite have intensively been studied, the mechanism behind the pronounced change in magnetization at the low-temperature transition is still debated. Here we report magnetization experiments on a pyrrhotite crystal (Fe 6.6 S 8 ) that consists of a 4C and an incommensurate 5C* superstructure that are different in their defect structure. The occurrence of two superstructures is magnetically confirmed by symmetric inflection points in hysteresis measurements above the transition at about 30 K. The disappearance of the inflection points and the associated change of the hysteresis parameters indicate that the two superstructures become strongly coupled to form a unitary magnetic anisotropy system at the transition. From this it follows that the Besnus transition in monoclinic pyrrhotite is an extrinsic magnetic phenomenon with respect to the 4C superstructure and therefore the physics behind it is in fact different from that of the well-known Verwey transition.

Description: Genetic variants in or near miRNA genes can have profound effects on miRNA expression and targeting. As user-friendly software for the impact prediction of miRNA variants on a large scale is still lacking, we created a tool called miRVaS. miRVaS automates this prediction by annotating the location of the variant relative to functional regions within the miRNA hairpin (seed, mature, loop, hairpin arm, flanks) and by annotating all predicted structural changes within the miRNA due to the variant. In addition, the tool defines the most important region that is predicted to have structural changes and calculates a conservation score that is indicative of the reliability of the structure prediction. The output is presented in a tab-separated file, which enables fast screening, and in an html file, which allows visual comparison between wild-type and variant structures. All separate images are provided for downstream use. Finally, we tested two different approaches on a small test set of published functionally validated genetic variants for their capacity to predict the impact of variants on miRNA expression.

Description: Analysis of RNA-seq data often detects numerous ‘non-co-linear’ (NCL) transcripts, which comprised sequence segments that are topologically inconsistent with their corresponding DNA sequences in the reference genome. However, detection of NCL transcripts involves two major challenges: removal of false positives arising from alignment artifacts and discrimination between different types of NCL transcripts ( trans -spliced, circular or fusion transcripts). Here, we developed a new NCL-transcript-detecting method (‘NCLscan’), which utilized a stepwise alignment strategy to almost completely eliminate false calls (〉98% precision) without sacrificing true positives, enabling NCLscan outperform 18 other publicly-available tools (including fusion- and circular-RNA-detecting tools) in terms of sensitivity and precision, regardless of the generation strategy of simulated dataset, type of intragenic or intergenic NCL event, read depth of coverage, read length or expression level of NCL transcript. With the high accuracy, NCLscan was applied to distinguishing between trans -spliced, circular and fusion transcripts on the basis of poly(A)- and nonpoly(A)-selected RNA-seq data. We showed that circular RNAs were expressed more ubiquitously, more abundantly and less cell type-specifically than trans -spliced and fusion transcripts. Our study thus describes a robust pipeline for the discovery of NCL transcripts, and sheds light on the fundamental biology of these non-canonical RNA events in human transcriptome.

Description: Surface nuclear magnetic resonance technique, also called magnetic resonance sounding (MRS), is an emerging geophysical method that can detect the presence and spatial variations of the subsurface water content directly. In this paper, we introduce the MRS central loop geometry, in which the receiver loop is smaller than the transmitter loop and placed in its centre. In addition, using a shielded receiver coil we show how this configuration greatly increases signal-to-noise ratio and improves the resolution of the subsurface layers compared to the typically used coincident loop configuration. We compare sensitivity kernels for different loop configurations and describe advantages of the MRS central loop geometry in terms of superior behaviour of the sensitivity function, increased sensitivity values, reduced noise level of the shielded receiver coil, improved resolution matrix and reduced instrument dead time. With no extra time and effort in the field, central-loop MRS makes it possible to reduce measurement time and to measure data in areas with high anthropogenic noise. The results of our field example agree well with the complementary data, namely airborne electromagnetics, borehole data, and the hydrologic model of the area.

Description: Frictional heating during earthquake rupture reveals important information on earthquake mechanisms and energy dissipation. The amount of annealing varies widely and is, as yet, poorly constrained. Here we use magnetic susceptibility versus temperature measurements during cycling to increasingly elevated temperatures to constrain the maximum temperature a slip zone has experienced. The case study comprises sheared clay cored from the Japan Trench subduction plate-boundary fault zone (décollement), which accommodated the large slip of the 2011 M w 9.0 Tohoku-oki earthquake. The décollement was cored during the Integrated Ocean Drilling Program (IODP) Expedition 343, the Japan Trench Fast Drilling Project (JFAST). Heating signatures with estimated maximum temperatures ranging from ~300 to over 500 °C are determined close to the multiple slip surfaces within the décollement. Since it is impossible to tie a specific slip surface to a certain earthquake, thermal evidence for the cumulative effect of several earthquakes is unveiled. This as yet preliminary rock magnetic ‘geothermometer’ would be a useful tool to detect seismic heating along faults that experienced medium temperature rise, a range which is difficult to assess with other approaches.

Description: Feldspars are the most abundant rock-forming minerals in the Earth's crust, but their magnetic properties have not been rigorously studied. This work focuses on the intrinsic magnetic anisotropy of 31 feldspar samples with various chemical compositions. Because feldspar is often twinned or shows exsolution textures, measurements were performed on twinned and exsolved samples as well as single crystals. The anisotropy is controlled by the diamagnetic susceptibility and displays a consistent orientation of principal susceptibility axes; the most negative or minimum susceptibility is parallel to [010], and the maximum (least negative) is close to the crystallographic [001] axis. However, the magnetic anisotropy is weak when compared to other rock-forming minerals, 1.53 x 10 –9 m 3 kg –1 at maximum. Therefore, lower abundance minerals, such as augite, hornblende or biotite, often dominate the bulk paramagnetic anisotropy of a rock. Ferromagnetic anisotropy is not significant in most samples. In the few samples that do show ferromagnetic anisotropy, the principal susceptibility directions of the ferromagnetic subfabric do not display a systematic orientation with respect to the feldspar lattice. These results suggest that palaeointensity estimates of the geomagnetic field made on single crystals of feldspar will not be affected by a systematic orientation of the ferromagnetic inclusions within the feldspar lattice.

Description: Feldspars are the most abundant rock-forming minerals in the Earth's crust, but their magnetic properties have not been rigorously studied. This work focuses on the intrinsic magnetic anisotropy of 31 feldspar samples with various chemical compositions. Because feldspar is often twinned or shows exsolution textures, measurements were performed on twinned and exsolved samples as well as single crystals. The anisotropy is controlled by the diamagnetic susceptibility and displays a consistent orientation of principal susceptibility axes; the most negative or minimum susceptibility is parallel to [010], and the maximum (least negative) is close to the crystallographic [001] axis. However, the magnetic anisotropy is weak when compared to other rock-forming minerals, 1.53 x 10 –9 m 3 kg –1 at maximum. Therefore, lower abundance minerals, such as augite, hornblende or biotite, often dominate the bulk paramagnetic anisotropy of a rock. Ferromagnetic anisotropy is not significant in most samples. In the few samples that do show ferromagnetic anisotropy, the principal susceptibility directions of the ferromagnetic subfabric do not display a systematic orientation with respect to the feldspar lattice. These results suggest that palaeointensity estimates of the geomagnetic field made on single crystals of feldspar will not be affected by a systematic orientation of the ferromagnetic inclusions within the feldspar lattice.

Description: Palaeomagnetic data obtained from archaeological materials are used for reconstructions of the Earth's magnetic field of the past millennia. While many studies tested the reliability of this recorder for palaeointensity only a few studies did this for direction. The study presents an archaeomagnetic and rock magnetic investigation applied to an experimental pottery kiln, which was operated in 2003 to produce stone ware. This kind of high-quality pottery needs a temperature of at least 1160 °C. Shortly before heating of the kiln direct absolute measurements of the absolute geomagnetic field vector have been carried out close to it. After cooling of the kiln 24 oriented palaeomagnetic samples have been taken. Although Curie temperatures are about 580 °C, that is the typical temperature for magnetite, thermal as well as alternating field demagnetisations reveal also a considerable amount of hematite as magnetic carrier. This mixture of magnetite and hematite is dominated by pseudo-single domain grains. Demagnetisation removed in some cases weak secondary components, but in most cases the specimens carried a single component thermoremanent magnetisation. The mean characteristic remanent magnetisation direction agrees on 95 per cent confidence level with the directly measured field direction. Archaeointensity was obtained from five specimens with the Thellier–Coe method and with the multiple-specimen palaeointensity domain-state corrected method. Six of these specimens also provided a result of the Dekkers–Böhnel method, which overestimated the archaeointensity by about 9 per cent compared to the direct value, while after correction for fraction the value agrees very well. For the multiple-specimen palaeointensity domain-state corrected method only fractions between 25 and 75 per cent have been used and specimens showing alteration have been excluded. Above 450 °C many specimens showed alteration of the magnetic grains. Because median destructive temperatures were often above this value in most cases the fraction was less than 50 per cent. Nevertheless the obtained intensity (48.48 ± 0.24 μ) is on 95 per cent confidence level in agreement with the direct observation. Behaviour of the specimens during the Thellier-experiments was not ideal because of narrow unblocking temperature spectra and alteration. Nevertheless, the obtained mean archaeointensity is also in agreement with the direct field observation. Here the relative palaeointensity error is about 6 per cent and very high compared the multiple-specimen palaeointensity domain-state corrected method. The investigation demonstrates that a pottery kiln can provide a very precise estimate of the ancient geomagnetic field vector.

Description: We present a new model of the radial (1-D) conductivity structure of Earth's mantle. This model is derived from more than 10 yr of magnetic measurements from the satellites Ørsted, CHAMP, SAC-C and the Swarm trio as well as the global network of geomagnetic observatories. After removal of core and crustal field as predicted by a recent field model, we fit the magnetic data with spherical harmonic coefficients describing ring current activity and associated induction effects and estimate global C -responses at periods between 1.5 and 150 d. The C -responses are corrected for 3-D effects due to induction in the oceans and inverted for a 1-D model of mantle conductivity using both probabilistic and deterministic methods. Very similar results are obtained, consisting of a highly resistive upper mantle, an increase in conductivity in and beneath the transition zone and a conductive lower mantle. Analysis of the Hessian of the cost function reveals that the data are most sensitive to structures at depths between 800 and 1200 km, in agreement with the results obtained from the probabilistic approach. Preliminary interpretation of the inverted conductivity structure based on laboratory-based conductivity profiles shows that the recovered structure in the lower mantle either requires higher temperatures or the presence of material of high conductivity related to ponding of carbonate melts below the transition zone.

Description: A prominent gravity and geoid low lies just south of the Lofoten peninsula in northern Norway, partly coinciding with the location of Proterozoic granites of the Transscandinavian Igneous Belt and being offset by ca . 100 km to the highest topography of northern Norway. The study area extends both onshore and offshore and lies at the transition between Archaean and Proterozoic lithosphere. The Palaeoproterozoic basement has been overthrusted by the Palaeozoic nappes of the Caledonian orogen and now forms the passive margin of the NE Atlantic. We investigate the gravity anomaly performing combined 3-D geophysical-petrological forward modelling of the lithosphere and sublithospheric upper mantle using the interactive modelling program LitMod3D. We include variations in thickness and composition of the lithospheric mantle in order to include the effects on the rifted margin adjoining the Baltic craton. We compare three possible origins of the anomaly: (i) a low-density upper crust, representing the northward extension of the Transscandinavian Igneous Belt, (ii) a lower crustal source formed by a Moho depression and (iii) a thick, depleted lithospheric mantle of possibly Archaean origin. A similar, yet wider and stronger gravity anomaly is found on the conjugate margin in northeastern Greenland. A shallow crustal source is most consistent with the geophysical data sets. A respective source of the granitic belt, however, is difficult to reconcile with the regional geology both in Fennoscandia and Greenland. An additional contribution form a deeper source is suggested.

Description: Detecting allelic biases from high-throughput sequencing data requires an approach that maximises sensitivity while minimizing false positives. Here, we present Allelome.PRO, an automated user-friendly bioinformatics pipeline, which uses high-throughput sequencing data from reciprocal crosses of two genetically distinct mouse strains to detect allele-specific expression and chromatin modifications. Allelome.PRO extends approaches used in previous studies that exclusively analyzed imprinted expression to give a complete picture of the ‘allelome’ by automatically categorising the allelic expression of all genes in a given cell type into imprinted, strain-biased, biallelic or non-informative. Allelome.PRO offers increased sensitivity to analyze lowly expressed transcripts, together with a robust false discovery rate empirically calculated from variation in the sequencing data. We used RNA-seq data from mouse embryonic fibroblasts from F1 reciprocal crosses to determine a biologically relevant allelic ratio cutoff, and define for the first time an entire allelome. Furthermore, we show that Allelome.PRO detects differential enrichment of H3K4me3 over promoters from ChIP-seq data validating the RNA-seq results. This approach can be easily extended to analyze histone marks of active enhancers, or transcription factor binding sites and therefore provides a powerful tool to identify candidate cis regulatory elements genome wide.

Description: Analysis of more than 10 yr of vertical magnetic transfer function (VTF) estimates obtained at 12 mid-latitude sites, located in different continents and tectonic settings, reveals significant temporal variations for a period range between approximately 250 and 2000 s. The most ubiquitous pattern is a seasonal modulation of the VTF element that relates the vertical to the horizontal north–south magnetic components ( Tx ), which shows a high peak around the June solstice (and a low peak around the December solstice) regardless of the location of the site. To quantify the influence of this source effect on the amplitude of VTFs, we modelled the temporal variations of VTFs using a function with dependence on season and magnetic activity indexes. The model shows that differences between VTF estimates obtained at seasonal peaks can reach 0.08 of Tx absolute values and that the effect increases with latitude and period. Seasonal variations are observed also in the VTF component relating vertical to horizontal east–west magnetic components ( Ty ), but here the pattern with respect to the geographic distribution is less clear. In addition to seasonal trends, we observe long-term modulations correlating with the 11-yr solar cycle at some sites. The influence of these external source effects should be taken into account, before attempting a geological interpretation of the VTFs. It can be misleading, for example, to combine or compare VTFs obtained from long-period geomagnetic data acquired at different seasons or years. An effective method to estimate and remove these source effects from VTFs is by comparison with temporal variations of VTFs from synchronously recorded data at sites located at similar latitude (〈5° of difference) and longitude (〈10° of difference). Source effects in temporal variations of VTFs can be identified as those patterns that exhibit similar amplitudes and significant correlation with the geomagnetic activity at all compared sites. We also provide a second-order polynomial which can be used to estimate the amplitude of the seasonal variations in the Tx component globally as a function of latitude.

Description: Directions recovered from palaeomagnetic samples are usually archived with some quantitative information about their precision, most often in the form of a so-called α 95 angle. Such angles are classically co-estimated with the recovered palaeomagnetic direction from a collection of samples providing individual estimates of this direction. In some instances, however, palaeomagnetic directions have to be inferred from a single sample in which case no α 95 angle can be recovered in this way. Fortunately, the progressive demagnetization techniques and principal component analysis universally used to recover directional information from single samples provide alternative measures of the error affecting the recovered direction, known as Maximum Angular Deviation (MAD) angles. These have so far only been considered as rough quality indicators. Here, however, we show that directions recovered in this way can be assumed to satisfy a Fisher distribution, and that the corresponding MAD angles can be rescaled into α 95 estimates by multiplying it by an appropriate factor, which only depends on the number of demagnetization steps used in the principal component analysis and on whether one relies on a standard or a so-called ‘anchored’ principal component analysis. These coefficients have been tabulated and practical recommendations for taking advantage of them outlined in the final section of the text. They provide simple means for users to produce much needed error bars on declination and inclination time series recovered from sedimentary long sequences.

Description: Identity by descent (IBD) can be reliably detected for long shared DNA segments, which are found in related individuals. However, many studies contain cohorts of unrelated individuals that share only short IBD segments. New sequencing technologies facilitate identification of short IBD segments through rare variants, which convey more information on IBD than common variants. Current IBD detection methods, however, are not designed to use rare variants for the detection of short IBD segments. Short IBD segments reveal genetic structures at high resolution. Therefore, they can help to improve imputation and phasing, to increase genotyping accuracy for low-coverage sequencing and to increase the power of association studies. Since short IBD segments are further assumed to be old, they can shed light on the evolutionary history of humans. We propose HapFABIA, a computational method that applies biclustering to identify very short IBD segments characterized by rare variants. HapFABIA is designed to detect short IBD segments in genotype data that were obtained from next-generation sequencing, but can also be applied to DNA microarray data. Especially in next-generation sequencing data, HapFABIA exploits rare variants for IBD detection. HapFABIA significantly outperformed competing algorithms at detecting short IBD segments on artificial and simulated data with rare variants. HapFABIA identified 160 588 different short IBD segments characterized by rare variants with a median length of 23 kb (mean 24 kb) in data for chromosome 1 of the 1000 Genomes Project. These short IBD segments contain 752 000 single nucleotide variants (SNVs), which account for 39% of the rare variants and 23.5% of all variants. The vast majority—152 000 IBD segments—are shared by Africans, while only 19 000 and 11 000 are shared by Europeans and Asians, respectively. IBD segments that match the Denisova or the Neandertal genome are found significantly more often in Asians and Europeans but also, in some cases exclusively, in Africans. The lengths of IBD segments and their sharing between continental populations indicate that many short IBD segments from chromosome 1 existed before humans migrated out of Africa. Thus, rare variants that tag these short IBD segments predate human migration from Africa. The software package HapFABIA is available from Bioconductor. All data sets, result files and programs for data simulation, preprocessing and evaluation are supplied at http://www.bioinf.jku.at/research/short-IBD .

Description: Glioma is the most common and fatal primary brain tumour with poor prognosis; however, the functional roles of miRNAs in glioma malignant progression are insufficiently understood. Here, we used an integrated approach to identify miRNA functional targets during glioma malignant progression by combining the paired expression profiles of miRNAs and mRNAs across 160 Chinese glioma patients, and further constructed the functional miRNA–mRNA regulatory network. As a result, most tumour-suppressive miRNAs in glioma progression were newly discovered, whose functions were widely involved in gliomagenesis. Moreover, three miRNA signatures, with different combinations of hub miRNAs (regulations≥30) were constructed, which could independently predict the survival of patients with all gliomas, high-grade glioma and glioblastoma. Our network-based method increased the ability to identify the prognostic biomarkers, when compared with the traditional method and random conditions. Hsa-miR-524-5p and hsa-miR-628-5p, shared by these three signatures, acted as protective factors and their expression decreased gradually during glioma progression. Functional analysis of these miRNA signatures highlighted their critical roles in cell cycle and cell proliferation in glioblastoma malignant progression, especially hsa-miR-524-5p and hsa-miR-628-5p exhibited dominant regulatory activities. Therefore, network-based biomarkers are expected to be more effective and provide deep insights into the molecular mechanism of glioma malignant progression.

Description: A major challenge in cancer genomics is uncovering genes with an active role in tumorigenesis from a potentially large pool of mutated genes across patient samples. Here we focus on the interactions that proteins make with nucleic acids, small molecules, ions and peptides, and show that residues within proteins that are involved in these interactions are more frequently affected by mutations observed in large-scale cancer genomic data than are other residues. We leverage this observation to predict genes that play a functionally important role in cancers by introducing a computational pipeline ( http://canbind.princeton.edu ) for mapping large-scale cancer exome data across patients onto protein structures, and automatically extracting proteins with an enriched number of mutations affecting their nucleic acid, small molecule, ion or peptide binding sites. Using this computational approach, we show that many previously known genes implicated in cancers are enriched in mutations within the binding sites of their encoded proteins. By focusing on functionally relevant portions of proteins—specifically those known to be involved in molecular interactions—our approach is particularly well suited to detect infrequent mutations that may nonetheless be important in cancer, and should aid in expanding our functional understanding of the genomic landscape of cancer.

Description: Genome duplication with hybridization, or allopolyploidization, occurs commonly in plants, and is considered to be a strong force for generating new species. However, genome-wide quantification of homeolog expression ratios was technically hindered because of the high homology between homeologous gene pairs. To quantify the homeolog expression ratio using RNA-seq obtained from polyploids, a new method named HomeoRoq was developed, in which the genomic origin of sequencing reads was estimated using mismatches between the read and each parental genome. To verify this method, we first assembled the two diploid parental genomes of Arabidopsis halleri subsp. gemmifera and Arabidopsis lyrata subsp. petraea ( Arabidopsis petraea subsp. umbrosa ), then generated a synthetic allotetraploid, mimicking the natural allopolyploid Arabidopsis kamchatica . The quantified ratios corresponded well to those obtained by Pyrosequencing. We found that the ratios of homeologs before and after cold stress treatment were highly correlated ( r = 0.870). This highlights the presence of nonstochastic polyploid gene regulation despite previous research identifying stochastic variation in expression. Moreover, our new statistical test incorporating overdispersion identified 226 homeologs (1.11% of 20 369 expressed homeologs) with significant ratio changes, many of which were related to stress responses. HomeoRoq would contribute to the study of the genes responsible for polyploid-specific environmental responses.

Description: Broad-band magnetic susceptibility (MS) measurement, a novel magnetic method capable of quantifying a narrow grain size distribution (GSD) of superparamagnetic (SP) particles by measuring low-field MS at a number of frequency steps spanning four orders of magnitude, has been tested in a loess/palaeosol section at Luochuan in the Chinese Loess Plateau. The studied succession consists of sequences from the latest palaeosol unit (S0) to the upper part of the loess unit (L2), spanning the last glacial–interglacial cycle. Reconstructed GSDs consist of volume fractions on the order of 10 –24 m 3 , and the mean GSDs are modal but with distinctive skewness among the loess, the weakly developed palaeosols (weak palaeosols), and the mature palaeosols. This indicates that the mean volume of SP particles in this loess/palaeosol sequence tends to increase during the transition from loess -〉 weak palaeosol -〉 palaeosol, an indication of grain growth as pedogenesis progresses. Total frequency dependence, or TFD(per cent), the difference between 130 at the lowest (130 Hz) and 500k at the highest (500 kHz) frequencies normalized to 130 , is judged to be a more suitable index than previous frequency dependence parameters for the concentration of SP particles. TFD(per cent) has a strong correlation with 130 , showing a continuous ‘growth curve’ with the rate of increase being highest for the loess, moderate for the weak palaeosols, and saturated for the palaeosols. The characteristic curve suggests that smaller SP particles are preferentially formed in the earlier stage of pedogenesis rather than the later phase when even larger particles are formed in mature palaeosols. These results demonstrate that the broad-band MS measurement method will be useful for the quantitative assessment of magnetic nanoparticles in soils and sediments.

Description: Environmental magnetism can provide simple and fast methods to semi-quantify the contamination levels by heavy metals in urban areas from the relationships between the magnetic properties and the heavy metal concentrations. The aim of this study is to explore the interpretation of selected magnetic parameters as proxies for the evaluation of heavy metal contamination in urban dusts. Dust samples were collected from different districts of Nanjing, Southeast China. The magnetic properties of dusts were analysed and heavy metal (Fe, Al, As, Ba, Cd, Cr, Co, Cu, Hg, Mn, Mo, Ni, Pb, Sr, V and Zn) concentrations were measured. Magnetic minerals, mainly from anthropogenic activities, are dominated by coarse stable single domain-grained ferrimagnetic minerals. The – T curves indicated that the main magnetic minerals are magnetite and hematite. Magnetic properties exhibited similar spatial distribution and significant positive correlations with the concentrations of most metals and with the pollution load index. Magnetic properties had a strong link with the abundance of heavy metals derived from industrial emission and traffic activities, whereas there was a poor correlation with the concentrations of metals from multi-sources such as commercial/domestic sources and natural processes. Besides the low-frequency magnetic susceptibility ( LF ) and the saturation isothermal remanent magnetization, other magnetic parameters can be considered as efficient indicators in the assessment of heavy metal contamination. The linkage between the magnetic properties and heavy metal concentrations in street dusts depends on the sources and nature of their magnetic fraction. Environmental magnetism analyses may provide simple and rapid methods for the assessment of heavy metal contamination in urban street dusts.

Description: We present in this paper an alternative approach to modelling and inversion of surface nuclear magnetic resonance (sNMR) data that has numerous attractive characteristics. By considering the forward and inverse problem in the frequency-domain (FD) it is possible to heavily truncate the band-limited data set with no appreciable loss of information. Furthermore, it introduces a natural method of signal demodulating and sidesteps the rotating frame transformation, which can be difficult to apply to large volumes. By instead remaining in the laboratory frame, phase shifts due to off-resonance transmission and static variation in the Larmor frequencies become simpler to account for. For these reasons, we present the first practical scheme for comprehensively inverting the complex sNMR data set. The use of complex data allows for improved depth and decay constant resolution. However, it requires the knowledge of subsurface electrical conductivity, and by corollary the ability to accurately model those effects in the sNMR record. Finally, the complex FD signal is sensitive to dephasing of the sNMR record due to static magnetic field inhomogeneity. We illustrate the ability to detect dephasing in synthetic and field data and demonstrate an effective inversion which is able to solve for static dephasing effects. This results in better estimation of meaningful decay constants as well as reduced error in extrapolated initial signal amplitude. Both of these parameters are critical for using sNMR in applications requiring estimation of hydraulic parameters.

Description: We report a palaeomagnetic investigation of the last full geomagnetic field reversal, the Matuyama-Brunhes (M-B) transition, as preserved in a continuous sequence of exposed lacustrine sediments in the Apennines of Central Italy. The palaeomagnetic record provides the most direct evidence for the tempo of transitional field behaviour yet obtained for the M-B transition. 40 Ar/ 39 Ar dating of tephra layers bracketing the M-B transition provides high-accuracy age constraints and indicates a mean sediment accumulation rate of about 0.2 mm yr –1 during the transition. Two relative palaeointensity (RPI) minima are present in the M-B transition. During the terminus of the upper RPI minimum, a directional change of about 180 ° occurred at an extremely fast rate, estimated to be less than 2 ° per year, with no intermediate virtual geomagnetic poles (VGPs) documented during the transit from the southern to northern hemisphere. Thus, the entry into the Brunhes Normal Chron as represented by the palaeomagnetic directions and VGPs developed in a time interval comparable to the duration of an average human life, which is an order of magnitude more rapid than suggested by current models. The reported investigation therefore provides high-resolution integrated palaeomagnetic and radioisotopic data that document the fine details of the anatomy and tempo of the M-B transition in Central Italy that in turn are crucial for a better understanding of Earth's magnetic field, and for the development of more sophisticated models that are able to describe its global structure and behaviour.

Description: The anisotropy of magnetic susceptibility (AMS) of soft rocks was measured in order to distinguish between the effect of remote and local strain fields, determine the size of the related inelastic damage zone and resolve the fault-plane solutions of past earthquakes. The AMS fabrics were explored next to late Pleistocene syndepositional normal faults (total displacement up to ~3.5 m) that cross soft lacustrine rocks within the seismically active Dead Sea basin. ‘Deposition fabrics’ prevail meters away from the fault planes and are characterized by scattered maximum and intermediate principal AMS axes. ‘Deformation fabrics’ are detected up to tens of centimetres from the fault planes and are characterized by well-grouped AMS axes, in which one of the principal axes is parallel to the strike of the nearby fault. Variations in the AMS fabrics and magnetic lineations define the size of the inelastic damage zone around the faults. The results demonstrate that the deformation-driven magnetic fabrics and the associated inelastic damage zones are compatible with coseismic dynamic faulting and the effects of the local strain field during earthquakes. Most of the AMS fabrics show a conspicuous similarity to that of the fault-plane solutions, i.e., the principal AMS axes and instantaneous strain ellipsoids are coaxial. These results suggest a novel application of the AMS method for defining the shape and size of the damage zones surrounding dynamic faults and determining the full tensor of the local strain field.

Description: Very low frequency (VLF) electromagnetic waves that penetrate conductive magma-filled dykes generate secondary fields on the surface that can be used to invert for dyke properties. The model used for the interpretation calculates currents induced in a conductive strip by an inducing field that decays exponentially with depth due to the conductivity of the surrounding medium. The differential equations are integrated to give an inhomogeneous Fredholm equation of the second kind with a kernel consisting of a modified Bessel function of the second kind. Numerical methods are typically used to solve for the induced currents in the strip. In this paper, we apply a modified Galerkin–Chebyshev method, which involves separating the kernel into source and field spectra and integrating the source terms to obtain a matrix equation for the unknown coefficients. The incident wave is expressed as a Chebyshev series. The modified Bessel function is separated into a logarithmic singularity and a non-singular remainder, both of which are expanded in complex Chebyshev polynomials. The Chebyshev coefficients for the remainder are evaluated using a fast Fourier transform, while the logarithmic term and incident field have analytic series. The deconvolution then involves a matrix inversion. The results depend on the ratio of strip-size to skin-depth. For infinite skin-depth and a singular conductivity distribution given by $\tau _0 a/\sqrt{a^2 - z^2 }$ (where 0 is the conductance, a is the half-length and z the distance from the centre), Parker gives an analytic solution. We present a similar analytic series solution for the finite skin-depth case, where the size to skin depth ratio is small. Results are presented for different ratios of size to skin depth that can be compared with numerical solutions. We compare full-space and half-space solutions. A fit of the model to VLF data taken above a magma filled dykes in Hawaii and Mt Etna demonstrates that while properties such as depth to top, conductivity ratio, tilt and dip can be determined, the depth to bottom is indeterminate due to the exponential decay of the inducing field.

Description: Correlated data errors (noise) are common in magnetotelluric (MT) data, but MT inversions typically neglect error correlations without investigating the impact of this simplification on inversion results. This paper examines effects of neglecting frequency- and spatially correlated noise on MT inversion, based on a nonlinear Bayesian formulation which quantifies the uncertainties of inversion results in terms of marginal posterior probability densities and credibility intervals. To do so, data with frequency- and spatially correlated noise of differing degrees are generated for several layered (1-D) synthetic cases. Bayesian MT inversions are carried out for these data sets with and without accounting for error correlation (i.e. applying full and diagonal covariance matrices, respectively, in the inversion), and the results are compared. For cases with noise that is strongly correlated over frequency or space, parameter uncertainties estimated using the diagonal-covariance simplification (neglecting error correlations) are found to often be significantly underestimated compared with results computed using the full covariance matrix.

Description: We report a new method to infer continuous time-series of the declination, inclination and intensity of the magnetic field from archaeomagnetic data. Adopting a Bayesian perspective, we need to specify a priori knowledge about the time evolution of the magnetic field. It consists in a time correlation function that we choose to be compatible with present knowledge about the geomagnetic time spectra. The results are presented as distributions of possible values for the declination, inclination or intensity. We find that the methodology can be adapted to account for the age uncertainties of archaeological artefacts and we use Markov chain Monte Carlo to explore the possible dates of observations. We apply the method to intensity data sets from Mari, Syria and to intensity and directional data sets from Paris, France. Our reconstructions display more rapid variations than previous studies and we find that the possible values of geomagnetic field elements are not necessarily normally distributed. Another output of the model is better age estimates of archaeological artefacts.

Description: We present a rigorous method for interpolation of electric and magnetic fields close to an interface with a conductivity contrast. The method takes into account not only a well-known discontinuity in the normal electric field, but also discontinuity in all the normal derivatives of electric and magnetic tangential fields. The proposed method is applied to marine 3-D controlled-source electromagnetic modelling (CSEM) where sources and receivers are located close to the seafloor separating conductive seawater and resistive formation. For the finite-difference scheme based on the Yee grid, the new interpolation is demonstrated to be much more accurate than alternative methods (interpolation using nodes on one side of the interface or interpolation using nodes on both sides, but ignoring the derivative jumps). The rigorous interpolation can handle arbitrary orientation of interface with respect to the grid, which is demonstrated on a marine CSEM example with a dipping seafloor. The interpolation coefficients are computed by minimizing a misfit between values at the nearest nodes and linear expansions of the continuous field components in the coordinate system aligned with the interface. The proposed interpolation operators can handle either uniform or non-uniform grids and can be applied to interpolation for both sources and receivers.

Description: According to field observations, electromagnetic (EM) signals accompanying seismic waves can be recorded. The orders of magnitude of observed coseismic electric and magnetic signals associated with earthquakes are usually around 1 to 10 1 mV km –1 and 10 –2 to 1 nT, respectively. In this paper, we carry out numerical simulation of coseismic EM signals associated with seismic waves due to electrokinetic effect and compare with field observations. The seismic source is represented by a finite fault measuring 15 x 9 km 2 with a max slip displacement 1.5 m, corresponding to a M w 5.9 earthquake. While using the EM surface boundary condition of continuous horizontal EM components, the magnetic signals only accompany the late-arriving S waves at receiver near the ground surface. This is obviously different from field observations. Thus, we adopt another EM surface boundary condition which assumes the ground surface carries surface charge. For the used half-space model, a surface-charge density magnitude | Q sc | in excess of 10 –4 C m –2 is sufficient to make horizontal magnetic components clearly show up at the whole time duration of seismic waves. When | Q sc | increases, the contribution of surface-charge density to coseismic EM signals becomes more and more dominant comparing with that of the seismically induced streaming-current. We estimate the Q sc expected at the Earth's surface might be a value between –5 x 10 –4 and –0.1 C m –2 by the comparison between numerical results and field observations. The vertical magnetic signals only accompany the late-arriving seismic waves, because they are theoretically only induced by SH wave. The field observation results of vertical magnetic signals may be resulted from the scattering effect or the seismic dynamo effect. We conclude that electrokinetic effect combined with surface-charge assumption is one possible generation mechanism of the observed coseismic EM signals.