ALBERT

Description: Previous formalisms for determining the static perturbation of spherically symmetric self-gravitating elastic Earth models due to displacement dislocations deal with each infinitesimal element of the fault system in its epicentral reference frame. In this work, we overcome this restriction and present novel and compact formulas for obtaining the perturbation due to the whole fault system in an arbitrary and common reference frame. Furthermore, we show that, even in an arbitrary reference frame, it is still possible to discriminate the contributions associated with the polar, bipolar and quadrupolar patterns of the seismic source response, as well as their relation with the along strike, along dip and tensile components of the displacement dislocation. These results allow a better understanding of the relation between the static perturbation and the whole fault system, and find direct applications in geodetic problems, like the modelling of long-wavelength geoid or gravity data from GRACE and GOCE space missions and of the perturbation of the deviatoric inertia tensor of the Earth.

Description: Sensitive instruments like strainmeters and tiltmeters are necessary for measuring slowly varying low amplitude Earth deformations. Nonetheless, laser and fibre interferometers are particularly suitable for interrogating such instruments due to their extreme precision and accuracy. In this paper, a practical design of a simple pendulum borehole tiltmeter based on laser fibre interferometric displacement sensors is presented. A prototype instrument has been constructed using welded borosilicate with a pendulum length of 0.85 m resulting in a main resonance frequency of 0.6 Hz. By implementing three coplanar extrinsic fibre Fabry-Perot interferometric probes and appropriate signal filtering, our instrument provides tilt measurements that are insensitive to parasitic deformations caused by temperature and pressure variations. This prototype has been installed in an underground facility (Rustrel, France) where results show accurate measurements of Earth strains derived from Earth and ocean tides, local hydrologic effects, as well as local and remote earthquakes. The large dynamic range and the high sensitivity of this tiltmeter render it an invaluable tool for numerous geophysical applications such as transient fault motion, volcanic strain and reservoir monitoring.

Description: Four antibiotics (pamamycin, oligomycin A, oligomycin B and echinosporin) were isolated and characterized from the fermentation broth of the marine Streptomyces strains B8496 and B8739. Bioassays revealed that each of these compounds impaired motility and caused subsequent lysis of P. viticola zoospores in a dose- and time-dependent manner. Pamamycin displayed the strongest motility inhibitory and lytic activities (IC 50 0.1 μg mL –1 ) followed by oligomycin B (IC 50 0.15 and 0.2 μg mL –1 ) and oligomycin F (IC 50 0.3 and 0.5 μg mL –1 ). Oligomycin A and echinosporin also showed motility inhibitory activities against the zoospores with IC 50 values of 3.0 and 10.0 μg mL –1 , respectively. This is the first report of motility inhibitory and lytic activities of these antibiotics against zoospores of a phytopathogenic peronosporomycete. Structures of all the isolated compounds were determined based on detailed spectroscopic analysis.

Description: In sulfidic environments, microbes oxidize reduced sulfur compounds via several pathways. We used metagenomics to investigate sulfur metabolic pathways from microbial mat communities in two subterranean sulfidic streams in Lower Kane Cave, WY, USA and from Glenwood Hot Springs, CO, USA. Both unassembled and targeted recA gene assembly analyses revealed that these streams were dominated by Epsilonproteobacteria and Gammaproteobacteria , including groups related to Sulfurovum , Sulfurospirillum , Thiothrix and an epsilonproteobacterial group with no close cultured relatives. Genes encoding sulfide:quinone oxidoreductase (SQR) were abundant at all sites, but the specific SQR type and the taxonomic affiliation of each type differed between sites. The abundance of thiosulfate oxidation pathway genes (Sox) was not consistent between sites, although overall they were less abundant than SQR genes. Furthermore, the Sox pathway appeared to be incomplete in all samples. This work reveals both variations in sulfur metabolism within and between taxonomic groups found in these systems, and the presence of novel epsilonproteobacterial groups.

Description: Pseudomonas aeruginosa is an opportunistic pathogen with high resistance to a wide variety of antimicrobials. The multidrug resistance pump MexAB-OprM promotes the efflux of various antibiotics, mostly when mutations accumulate in the transcriptional regulators MexR, NalC and NalD, thereby causing MexAB-OprM overexpression. In this work, a characterization of 50 P. aeruginosa isolates obtained from Brazilian agricultural soils to determine the reasons of their resistance to aztreonam was done. The majority of the isolates showed higher aztreonam resistance than wild-type strain by MIC method. DNA sequence analysis of mexR , nalC and nalD genes from 13 of these isolates showed the amino acid substitution in NalC for all tested isolates, just one mutation was detected in MexR and none in NalD. Furthermore, an increase in the level of mexA expression by real-time RT-PCR analysis in eight isolates harboring mutations in NalC was found. Although there was not a relationship between MIC of aztreonam and the level of mexA expression, on the other hand, the results presented here suggest that novel mutations in NalC, including Arg 97 -Gly and Ala 186 -Thr, are related to MexAB-OprM overexpression causing aztreonam resistance in P. aeruginosa environmental isolates.

Description: Sedge-dominated wetlands on the Qinghai–Tibetan Plateau are methane emission centers. Methanotrophs at these sites play a role in reducing methane emissions, but relatively little is known about the composition of active methanotrophs in these wetlands. Here, we used DNA stable isotope probing to identify the key active aerobic methanotrophs in three sedge-dominated wetlands on the plateau. We found that Methylocystis species were active in two peatlands, Hongyuan and Dangxiong. Methylobacter species were found to be active only in Dangxiong peat. Hongyuan peat had the highest methane oxidation rate, and cross-feeding of carbon from methanotrophs to methylotrophic Hyphomicrobium species was observed. Owing to a low methane oxidation rate during the incubation, the labeling of methanotrophs in Maduo wetland samples was not detected. Our results indicate that there are large differences in the activity of methanotrophs in the wetlands of this region.

Description: Here we present the generation and function of two sets of bacterial plasmids that harbor fluorescent genes encoding either blue, cyan, yellow or red fluorescent proteins. In the first set, protein expression is controlled by the strong and constitutive nptII promoter whereas in the second set, the strong tac promoter was chosen that underlies LacI q regulation. Furthermore, the plasmids are mobilizable, contain Tn 7 transposons and a temperature-sensitive origin of replication. Using Escherichia coli S17-1 as donor strain, the plasmids allow fast and convenient Tn 7 -transposon delivery into many enterobacterial hosts, such as the here-used E. coli O157:H7. This procedure omits the need of preparing competent recipient cells and antibiotic resistances are only transiently conferred to the recipients. As the fluorescence proteins show little to no overlap in fluorescence emission, the constructs are well suited for the study of multicolored synthetic bacterial communities during biofilm production or in host colonization studies, e.g. of plant surfaces. Furthermore, tac promoter-reporter constructs allow the generation of so-called reproductive success reporters, which allow to estimate past doublings of bacterial individuals after introduction into environments, emphasizing the role of individual cells during colonization.

Description: The geocentre motion is the motion of the centre of mass of the entire Earth, considered an isolated system, in a terrestrial system of reference. We first derive a formula relating the harmonic degree-1 Lagrangian variation of the gravity at a station to both the harmonic degree-1 vertical displacement of the station and the displacement of the whole Earth's centre of mass. The relationship is independent of the nature of the Earth deformation and is valid for any source of deformation. We impose no constraint on the system of reference, except that its origin must initially coincide with the centre of mass of the spherically symmetric Earth model. Next, we consider the geocentre motion caused by surface loading. In a system of reference whose origin is the centre of mass of the solid Earth, we obtain a specific relationship between the gravity variation at the surface, the geocentre displacement and the load Love number $h^{\prime }_1$ , which demands the Earth's structure and rheological behaviour be known. For various networks of real or fictitious stations, we invert synthetic signals of surface gravity variations caused by atmospheric loading to retrieve the degree-1 variation of gravity. We then select six well-distributed stations of the Global Geodynamics Project, which is a world network of superconducting gravimeters, to invert actual gravity data for the degree-1 variations and determine the geocentre displacement between the end of 2004 and the beginning of 2012, assuming it to be due to surface loading. We find annual and semi-annual displacements with amplitude 0.5–2.3 mm.

Description: Spa -typing and microarray techniques were used to study epidemiological changes in methicillin-resistant Staphylococcus aureus (MRSA) in South-East Austria. The population structure of 327 MRSA isolated between 2002 and 2012 was investigated. MRSA was assigned to 58 different spa types and 14 different MLST CC (multilocus sequence type clonal complexes); in particular, between 2007 and 2012, an increasing diversity in MRSA clones could be observed. The most abundant clonal complex was CC5. On the respective SCC mec cassettes, the CC5 isolates differed clearly within this decade and CC5/SCC mec I, the South German MRSA, predominant in 2002, was replaced by CC5/SCC mec II, the Rhine-Hesse MRSA in 2012. Whereas in many European countries MLST CC22-MRSA (EMRSA 15, the Barnim epidemic MRSA) is predominant, this clone occurred in Austria nearly 10 years later than in neighbouring countries. CC45, the Berlin EMRSA, epidemic in Germany, was only sporadically found in South-East Austria. The Irish ST8-MRSA-II represented by spa -type t190 was frequently found in 2002 and 2007, but disappeared in 2012. Our results demonstrate clonal replacement of MRSA clones within the last years in Austria. Ongoing surveillance is warranted for detection of changes within the MRSA population.

Description: This study aimed to investigate the effects of dietary fibre sources on the gut microbiota in suckling piglets, and to test the hypothesis that a moderate increase of dietary fibre may affect the gut microbiota during the suckling period. Suckling piglets were fed different fibre-containing diets or a control diet from postnatal day 7 to 22. Digesta samples from cecum, proximal colon and distal colon were used for Pig Intestinal Tract Chip analysis. The data showed that the effects of fibre-containing diet on the gut microbiota differed in the fibre source and gut location. The alfalfa diet increased Clostridium cluster XIVb and Sporobacter termitidis in the cecum compared to the pure cellulose diet. Compared to the control diet, the alfalfa diet also increased Coprococcus eutactus in the distal colon, while the pure cellulose diet decreased Eubacterium pyruvativorans in the cecum. The pure cellulose diet increased Prevotella ruminicola compared to the wheat bran diet. Interestingly, the alfalfa group had the lowest abundance of the potential pathogen Streptococcus suis in the cecum and distal colon. These results indicated that a moderate increase in dietary fibres affected the microbial composition in suckling piglets, and that the alfalfa inclusion produced some beneficial effects on the microbial communities.

Description: One function of the gut microbiota gaining recent attention, especially in herbivorous mammals and insects, is the metabolism of plant secondary metabolites (PSMs). We investigated whether this function exists within the gut communities of a specialist avian herbivore. We sequenced the cecal metagenome of the Greater Sage-Grouse ( Centrocercus urophasianus ), which specializes on chemically defended sagebrush ( Artemisia spp.). We predicted that the cecal metagenome of the sage-grouse would be enriched in genes associated with the metabolism of PSMs when compared to the metagenome of the domestic chicken. We found that representation of microbial genes associated with ‘xenobiotic degradation and metabolism’ was 3-fold higher in the sage-grouse cecal metagenomes when compared to that of the domestic chicken. Further, we identified a complete metabolic pathway for the degradation of phenol to pyruvate, which was not detected in the metagenomes of the domestic chicken, bovine rumen or 14 species of mammalian herbivores. Evidence of monoterpene degradation (a major class of PSMs in sagebrush) was less definitive, although we did detect genes for several enzymes associated with this process. Overall, our results suggest that the gut microbiota of specialist avian herbivores plays a similar role to the microbiota of mammalian and insect herbivores in degrading PSMs.

Description: Intracellular endosymbiotic bacteria are common and can play a crucial role for insect pathology. Therefore, such bacteria could be a potential key to our understanding of major losses of Western honey bees ( Apis mellifera ) colonies. However, the transmission and potential effects of endosymbiotic bacteria in A. mellifera and other Apis spp. are poorly understood. Here, we explore the prevalence and transmission of the genera Arsenophonus , Wolbachia , Spiroplasma and Rickettsia in Apis spp. Colonies of A. mellifera ( N = 33, with 20 eggs from worker brood cells and 100 adult workers each) as well as mated honey bee queens of A. cerana , A. dorsata and A. florea ( N = 12 each) were screened using PCR. While Wolbachia , Spiroplasma and Rickettsia were not detected, Arsenophonus spp. were found in 24.2% of A. mellifera colonies and respective queens as well as in queens of A. dorsata (8.3%) and A. florea (8.3%), but not in A. cerana . The absence of Arsenophonus spp. from reproductive organs of A. mellifera queens and surface-sterilized eggs does not support transovarial vertical transmission. Instead, horizontal transmission is most likely.

Description: We review the theory of the Earth's elastic and gravitational response to a surface disk load. The solutions for displacement of the surface and the geoid are developed using expansions of Legendre polynomials, their derivatives and the load Love numbers. We provide a matlab  function called diskload that computes the solutions for both uncompensated and compensated disk loads. In order to numerically implement the Legendre expansions, it is necessary to choose a harmonic degree, n max , at which to truncate the series used to construct the solutions. We present a rule of thumb (ROT) for choosing an appropriate value of n max , describe the consequences of truncating the expansions prematurely and provide a means to judiciously violate the ROT when that becomes a practical necessity.

Description: Two types of signals are clearly visible in continuous GPS (cGPS) time-series in Iceland, in particular in the vertical component. The first one is a yearly seasonal cycle, usually sinusoid-like with a minimum in the spring and a maximum in the fall. The second one is a trend of uplift, with higher values the closer the cGPS stations are to the centre of Iceland and ice caps. Here, we study the seasonal cycle signal by deriving its average at 71 GPS sites in Iceland. We estimate the annual and semi-annual components of the cycle in their horizontal and vertical components using a least-squares adjustment. The peak-to-peak amplitude of the cycle of the vertical component at the studied sites ranges from 4 mm near the coastline up to 27 mm at the centre of the Vatnajökull, the largest ice cap in Iceland. The minimum of the seasonal cycle occurs earlier in low lying areas than in the central part of Iceland, consistent with snow load having a large influence on seasonal deformation. Modelling shows that the seasonal cycle is well explained by accounting for elastically induced surface displacements due to snow, atmosphere, reservoir lake and ocean variations. Model displacement fields are derived considering surface loads on a multilayered isotropic spherical Earth. Through forward and inverse modelling, we were able to reproduce a priori information on the average seasonal cycle of known loads (atmosphere, snow in non-glaciated areas and lake reservoir) and get an estimation of other loads (glacier mass balance and ocean). The seasonal glacier mass balance cycle in glaciated areas and snow load in non-glaciated areas are the main contributions to the seasonal deformation. For these loads, induced seasonal vertical displacements range from a few millimetres far from the loads in Iceland, to more than 20 mm at their centres. Lake reservoir load also has to be taken into account on local scale as it can generate up to 20 mm of vertical deformation. Atmosphere load and ocean load are observable and generate vertical displacements in the order of a few millimetres. Inversion results also shows that the Iceland crust is less rigid than the world average. Interannual deviation from the GPS seasonal cycle can occur and are caused by unusual weather conditions over extended period of time.

Description: Wood-rotting fungi possess remarkably diverse extracellular oxidation mechanisms, including enzymes, such as laccase and peroxidases, and Fenton chemistry. The ability to biologically drive Fenton chemistry by the redox cycling of quinones has previously been reported to be present in both ecologically diverging main groups of wood-rotting basidiomycetes. Therefore, we investigated whether it is even more widespread among fungal organisms. Screening of a diverse selection of a total of 18 ascomycetes and basidiomycetes for reduction of the model compound 2,6-dimethoxy benzoquinone revealed that all investigated strains were capable of reducing it to its corresponding hydroquinone. In a second step, depolymerization of the synthetic polymer polystyrene sulfonate was used as a proxy for quinone-dependent Fenton-based biodegradation capabilities. A diverse subset of the strains, including environmentally ubiquitous molds, white-rot fungi, as well as peatland and aquatic isolates, caused substantial depolymerization indicative for the effective employment of quinone redox cycling as biodegradation tool. Our results may also open up new paths to utilize diverse fungi for the bioremediation of recalcitrant organic pollutants.

Description: Ice-binding proteins (IBPs), such as antifreeze proteins (AFPs) and ice-nucleating proteins (INPs), have been described in diverse cold-adapted organisms, and their potential applications in biotechnology have been recognized in various fields. Currently, both IBPs are being applied to biotechnological processes, primarily in medicine and the food industry. However, our knowledge regarding the diversity of bacterial IBPs is limited; few studies have purified and characterized AFPs and INPs from bacteria. Phenotypically verified IBPs have been described in members belonging to Gammaproteobacteria, Actinobacteria and Flavobacteriia classes, whereas putative IBPs have been found in Gammaproteobacteria, Alphaproteobacteria and Bacilli classes. Thus, the main goal of this minireview is to summarize the current information on bacterial IBPs and their application in biotechnology, emphasizing the potential application in less explored fields such as agriculture. Investigations have suggested the use of INP-producing bacteria antagonists and AFPs-producing bacteria (or their AFPs) as a very attractive strategy to prevent frost damages in crops. UniProt database analyses of reported IBPs (phenotypically verified) and putative IBPs also show the limited information available on bacterial IBPs and indicate that major studies are required.

Description: Triazophos is a broad-spectrum and highly effective insecticide, and the residues of triazophos have been frequently detected in the environment. A triazophos-degrading bacterium, Burkholderia sp. SZL-1, was isolated from a long-term triazophos-polluted soil. Strain SZL-1 could hydrolyze triazophos to 1-phenyl-3-hydroxy-1,2,4-triazole, which was further utilized as the carbon sources for growth. The triazophos hydrolase gene trhA , cloned from strain SZL-1, was expressed and homogenously purified using Ni-nitrilotriacetic acid affinity chromatography. TrhA is 55 kDa and displays maximum activity at 25°C, pH 8.0. This enzyme still has nearly 60% activity at the range of 15°C–50°C for 30 min. TrhA was mutated by sequential error prone PCR and screened for improved activity for triazophos degradation. One purified variant protein (Val89-Gly89) named TrhA-M1 showed up to 3-fold improvement in specific activity against triazophos, and the specificity constants of K cat and K cat / K m for TrhA-M1 were improved up to 2.3- and 8.28-fold, respectively, compared to the wild-type enzyme. The results in this paper provided potential material for the contaminated soil remediation and hydrolase genetic structure research.

Description: The metal mining industry faces many large challenges in future years, among which is the increasing need to process low-grade ores as accessible higher grade ores become depleted. This is against a backdrop of increasing global demands for base and precious metals, and rare earth elements. Typically about 99% of solid material hauled to, and ground at, the land surface currently ends up as waste (rock dumps and mineral tailings). Exposure of these to air and water frequently leads to the formation of acidic, metal-contaminated run-off waters, referred to as acid mine drainage, which constitutes a severe threat to the environment. Formation of acid drainage is a natural phenomenon involving various species of lithotrophic (literally ‘rock-eating’) bacteria and archaea, which oxidize reduced forms of iron and/or sulfur. However, other microorganisms that reduce inorganic sulfur compounds can essentially reverse this process. These microorganisms can be applied on industrial scale to precipitate metals from industrial mineral leachates and acid mine drainage streams, resulting in a net improvement in metal recovery, while minimizing the amounts of leachable metals to the tailings storage dams. Here, we advocate that more extensive exploitation of microorganisms in metal mining operations could be an important way to green up the industry, reducing environmental risks and improving the efficiency and the economy of metal recovery.

Description: Data from the Gravity Recovery and Climate Experiment (GRACE) satellite mission can be used to estimate the mass change rate for separate drainage systems (DSs) of the Greenland Ice Sheet (GrIS). One approach to do so is by inversion of the level-2 spherical harmonic data to surface mass changes in predefined regions, or mascons. However, the inversion can be numerically unstable for some individual DSs. This occurs mainly for DSs with a small mass change signal that are located in the interior region of Greenland. In this study, we present a modified mascon inversion approach with an improved implementation of the constraint equations to obtain better estimates for individual DSs. We use separate constraints for mass change variability in the coastal zone, where run-off takes place, and for the ice sheet interior above 2000 m, where mass changes are smaller. A multi-objective optimization approach is used to find optimal prior variances for these two areas based on a simulation model. Correlations between adjacent DSs are suppressed when our optimized prior variances are used, while the mass balance estimates for the combination of the DSs that make up the GrIS above 2000 m are not affected significantly. The resulting mass balance estimates for some DSs in the interior are significantly improved compared to an inversion with a single constraint, as determined by a comparison with mass balance estimates from surface mass balance modelling and discharge measurements. The rate of mass change of the GrIS for the period of January 2003 to December 2012 is found to be –266.1 ± 17.2 Gt yr –1 in the coastal zone and areas below 2000 m, and +8.2 ± 8.6 Gt yr –1 in the interior region.

Description: The static and transient deformations produced by earthquakes cause density perturbations which, in turn, generate immediate, long-range perturbations of the Earth's gravity field. Here, an analytical solution is derived for gravity perturbations produced by a point double-couple source in homogeneous, infinite, non-self-gravitating elastic media. The solution features transient gravity perturbations that occur at any distance from the source between the rupture onset time and the arrival time of seismic P waves, which are of potential interest for real-time earthquake source studies and early warning. An analytical solution for such prompt gravity perturbations is presented in compact form. We show that it approximates adequately the prompt gravity perturbations generated by strike-slip and dip-slip finite fault ruptures in a half-space obtained by numerical simulations based on the spectral element method. Based on the analytical solution, we estimate that the observability of prompt gravity perturbations within 10 s after rupture onset by current instruments is severely challenged by the background microseism noise but may be achieved by high-precision gravity strainmeters currently under development. Our analytical results facilitate parametric studies of the expected prompt gravity signals that could be recorded by gravity strainmeters.

Description: Seismic waves produced by fault ruptures give rise to gravity perturbations. So far, these perturbations have either been modelled as permanent coseismic gravity change in a half-space or spherical Earth model, or as full time-domain model in infinite space. In this paper, we present the explicit solution of gravity perturbations in time domain produced by a double-couple buried in a homogeneous half-space. This result is especially suited to study gravity perturbations up to a few hundreds of kilometres from the epicentre. It facilitates detailed parametric studies of gravity perturbations from fault rupture, and predicts gravity perturbations of real earthquakes with greatly improved accuracy. The results may serve to develop first designs of gravity-assisted earthquake early-warning systems, made possible by a new generation of ultrasensitive gravity gradiometers, which is currently under development.

Description: Continuous gravimetric observations have been made with three successive generations of superconducting gravimeter over 20 yr at Syowa Station ( $39.6\deg$ E, $69.0\deg$ S), East Antarctica. The third-generation instrument, OSG#058, was installed in January 2010 and was calibrated by an absolute gravimeter during January and February, 2010. The estimated scale factor was –73.823 ± 0.053 μGal V –1 (1 μGal = 10 –8 m s –2 ). The first 5 yr of OSG#058 data from 2010 January 7 to 2015 January 10 were decomposed into tidal waves (M3 to Ssa) and other non-tidal components by applying the Bayesian tidal analysis program BAYTAP. Long-term non-tidal gravity residuals, which were obtained by subtracting annual and 18.6 year tidal waves and the predicted gravity response to the Earth's variable rotation, showed significant correlation with the accumulated snow depth measured at Syowa Station. The greatest correlation occurred when the gravity variations lagged the accumulated snow depth by 21 d. To estimate the gravitational effect of the accumulated snow mass, we inferred a conversion factor of 3.13 ± 0.08 μGal m –1 from this relation. The accumulated snow depth at Syowa Station was found to represent an extensive terrestrial water storage (the snow accumulation) around Syowa Station, which was estimated from the Gravity Recovery and Climate Experiment satellite gravity data. The snow accumulation around Syowa Station was detectable by the superconducting gravimeter.

Description: The new release AIUB-RL02 of monthly gravity models from GRACE GPS and K-Band range-rate data is based on reprocessed satellite orbits referring to the reference frame IGb08. The release is consistent with the IERS2010 conventions. Improvements with respect to its predecessor AIUB-RL01 include the use of reprocessed (RL02) GRACE observations, new atmosphere and ocean dealiasing products (RL05), an upgraded ocean tide model (EOT11A), and the interpolation of shallow ocean tides (admittances). The stochastic parametrization of AIUB-RL02 was adapted to include daily accelerometer scale factors, which drastically reduces spurious signal at the 161 d period in C 20 and at other low degree and order gravity field coefficients. Moreover, the correlation between the noise in the monthly gravity models and solar activity is considerably reduced in the new release. The signal and the noise content of the new AIUB-RL02 monthly gravity fields are studied and calibrated errors are derived from their non-secular and non-seasonal variability. The short-period time-variable signal over the oceans, mostly representing noise, is reduced by 50 per cent with respect to AIUB-RL01. Compared to the official GFZ-RL05a and CSR-RL05 monthly models, the AIUB-RL02 stands out by its low noise at high degrees, a fact emerging from the estimation of seasonal variations for selected river basins and of mass trends in polar regions. Two versions of the monthly AIUB-RL02 gravity models, with spherical harmonics resolution of degree and order 60 and 90, respectively, are available for the time period from March 2003 to March 2014 at the International Center for Global Earth Models or from ftp://ftp.unibe.ch/aiub/GRAVITY/GRACE (last accessed 22 March 2016).

Description: The relative gravimeter is the primary terrestrial instrument for measuring spatially and temporally varying gravitational fields. The background noise of the instrument—that is, non-linear drift and random tares—typically requires some form of least-squares network adjustment to integrate data collected during a campaign that may take several days to weeks. Here, we present an approach to remove the change in the observed relative-gravity differences caused by hydrologic or other transient processes during a single campaign, so that the adjusted gravity values can be referenced to a single epoch. The conceptual approach is an example of coupled hydrogeophysical inversion, by which a hydrologic model is used to inform and constrain the geophysical forward model. The hydrologic model simulates the spatial variation of the rate of change of gravity as either a linear function of distance from an infiltration source, or using a 3-D numerical groundwater model. The linear function can be included in and solved for as part of the network adjustment. Alternatively, the groundwater model is used to predict the change of gravity at each station through time, from which the accumulated gravity change is calculated and removed from the data prior to the network adjustment. Data from a field experiment conducted at an artificial-recharge facility are used to verify our approach. Maximum gravity change due to hydrology (observed using a superconducting gravimeter) during the relative-gravity field campaigns was up to 2.6 μGal d –1 , each campaign was between 4 and 6 d and one month elapsed between campaigns. The maximum absolute difference in the estimated gravity change between two campaigns, two months apart, using the standard network adjustment method and the new approach, was 5.5 μGal. The maximum gravity change between the same two campaigns was 148 μGal, and spatial variation in gravity change revealed zones of preferential infiltration and areas of relatively high groundwater storage. The accommodation for spatially varying gravity change would be most important for long-duration campaigns, campaigns with very rapid changes in gravity and (or) campaigns where especially precise observed relative-gravity differences are used in the network adjustment.

Description: We document two kinds of traveling ionospheric disturbances, namely, CTIDs (Co-tsunami-Traveling-Ionospheric-disturbances) and ATIDs (Ahead-of-Tsunami-Traveling-Ionospheric-disturbances) related to the Tohoku-Oki tsunami of 2011 March 11. They are referred to the disturbances that remain behind and ahead of the principal tsunami wave front, respectively. We first note their presence in a numerical experiment performed using a simulation code coupling the tsunami, atmosphere and ionosphere. This code uses the tsunami wavefield as an input and simulates acoustic-gravity waves (AGWs) in the atmosphere and TIDs, in the form of total electron content (TEC) disturbance, in the ionosphere. The simulated TEC reveals the excitation of CTIDs (at about 2 TECU) and ATIDs (at about 1 TECU), representing up to 5 per cent disturbance over the ambient electron density, and they arise from the dissipation of AGWs in the thermosphere. A novel outcome is that during the tsunami passage between ~6° and 12° of epicentral distance, strong ATIDs arrive ~20–60 min ahead of the tsunami wave front covering ~3°–10° of distance from the tsunami location. Simulation results are compared with the far-field observations using GNSS satellites and confirm that ATIDs are the first detected TEC maximum, occurring 20–60 min ahead of the tsunami arrival. Our simulation also confirms the presence of largest TEC maximum representing CTIDs, 10–20 min after the first tsunami wave. ATIDs reported in this study have characteristics that can be potentially used for the early warning of the tsunami.

Description: Differential inhibitors are important for measuring the relative contributions of microbial groups, such as ammonia-oxidizing bacteria (AOB) and ammonia-oxidizing archaea (AOA), to biogeochemical processes in environmental samples. In particular, 2-phenyl-4,4,5,5-tetramethylimidazoline-1-oxyl 3-oxide (PTIO) represents a nitric oxide scavenger used for the specific inhibition of AOA, implicating nitric oxide as an intermediate of thaumarchaeotal ammonia oxidation. This study investigated four alternative nitric oxide scavengers for their ability to differentially inhibit AOA and AOB in comparison to PTIO. Caffeic acid, curcumin, methylene blue hydrate and trolox were tested on Nitrosopumilus maritimus , two unpublished AOA representatives (AOA-6f and AOA-G6) as well as the AOB representative Nitrosomonas europaea . All four scavengers inhibited ammonia oxidation by AOA at lower concentrations than for AOB. In particular, differential inhibition of AOA and AOB by caffeic acid (100 μM) and methylene blue hydrate (3 μM) was comparable to carboxy-PTIO (100 μM) in pure and enrichment culture incubations. However, when added to aquarium sponge biofilm microcosms, both scavengers were unable to inhibit ammonia oxidation consistently, likely due to degradation of the inhibitors themselves. This study provides evidence that a variety of nitric oxide scavengers result in differential inhibition of ammonia oxidation in AOA and AOB, and provides support to the proposed role of nitric oxide as a key intermediate in the thaumarchaeotal ammonia oxidation pathway.

Description: In the present work we illustrate a new local inversion algorithm to retrieve the Moho depth from GOCE (Gravity field and steady-state Ocean Circulation Explorer) gravity field. In details the proposed procedure can be divided into two main steps: the first one consists in recognizing and isolating the different geological provinces in the study area by exploiting information coming from the GOCE global gravity field model. Once the main geological provinces are defined, a function relating the crust density of each province with depth is built and used to reduce the data. The gravitational effects of sediments, topography, bathymetry and upper mantle are also removed. In the second step the residual gravitational field is inverted to retrieve the Moho depth and some information on the crustal density. In particular, the clustering of geological province is performed by means of an automatic Bayesian classification algorithm while the inversion of GOCE residual field is performed by adapting the global algorithm developed in the framework of the GEMMA project to the local scale. The procedure, based on an iterative Wiener filter, allows to compute the Moho depth considering lateral as well as radial variations of crustal density. The algorithm has been applied to the fifth release of GOCE time-wise global gravity field model to infer information on the crustal structure in the Western Balkan area, that is, the region laying between Bulgaria and the Adriatic Sea. This region is one of the most complex and active, from the tectonic point of view, in the whole Europe and it is characterized by the presence of the Alpine-Himalayan orogenic belt, formed by the collision between the African and Eurasian plates, and by the opening of the Pannonian Basin. Results show a good agreement between the obtained geological provinces with the actual knowledge on the region. The resulting Moho depth ranges between about 20 km beneath the Adriatic Sea and 45 km in the Dinarides. Comparisons with available seismic data show differences smaller than 1 km (standard deviation).

Description: Gravity variations associated with Earth's oblateness ( J 2 ) have been observed by satellite laser ranging (SLR) since 1976. The J 2 time-series has been used to measure and help understand many geophysical processes within the Earth system ranging from the mantle to the atmosphere. While post glacial rebound and the Earth climate system are believed to be the primary driving forces of long-term and seasonal J 2 variations, the physical cause of decadal and longer timescale J 2 variations has remained uncertain, although recent evidence indicates that polar ice mass changes are important. In this study, we estimate a variety of climate contributions to J 2 over the period 1979–2010, and find that ice mass variations in Greenland and Antarctica are the dominant cause of observed decadal and longer J 2 variations. Residual variations at periods near 10–11 years may reflect limitations of numerical climate models in estimating mass change variability at long periods, but are also suggestive of potential contribution related to variable solar activity.

Description: Interferometric synthetic aperture radar (InSAR) technology provides a valuable tool for obtaining Earth surface deformation and topography at high spatial resolution for crustal deformation studies. Similar to global positioning system (GPS), InSAR measurements are affected by the Earth's ionospheric and tropospheric layers as the electromagnetic signals significantly refract while propagating through the different layers. While GPS signals propagating through the neutral atmosphere are affected primarily by the distribution, pressure and temperature of atmospheric gases, including water vapour, the propagation through the ionosphere is mainly affected by the number of free electrons along the signal path. Here, we present the use of dense regional GPS networks for extracting tropospheric zenith delays and ionospheric total electron content (TEC) maps in order to reduce the noise levels in InSAR images. The results show significant reduction in the root mean square (RMS) values when simultaneously combining the two corrections, both at short time periods where no surface deformation is expected, and at longer periods, where imaging of localized subsidence and fault creep is enhanced.

Description: The main objective of the Gravity Recovery and Climate Experiment (GRACE) Atmospheric and Oceanic De-Aliasing Level-1B product (AOD1B) is the removal of high-frequency non-tidal mass variations due to sub-monthly mass transport in the atmosphere and oceans. Application of AOD1B shall avoid aliasing of these high-frequency signals into monthly gravity models derived from modern gravity missions and shall help to derive consistent orbit solutions for altimetry and Satellite Laser Ranging missions. The AOD1B 6-h series of spherical harmonic coefficients up to degree and order 100 are routinely generated at the German Research Centre for Geoscience and distributed to the GRACE Science Data System and the user community. Inputs for this product are acquired from numerical weather prediction models which are regularly revised and consequently not stable in time. The latest AOD1B release 5 (RL05) is based, as all other releases, on input from ECMWF and does not resolve this problem of discontinuities present in the surface pressure and surface geopotential input data. This might contaminate the gravity field variations derived from atmospheric mass variations. In this paper we present a method to overcome this problem during future AOD1B product generation, as well as two new Level-2 products (GAE and GAF) that, over land, fix a posteriori the two jumps present in the already distributed Level-2 RL05 monthly gravity models which were based on AOD1B RL05. The impact of the proposed correction on the variations and long-term trend of the total mass of the atmosphere and on the ice mass balance over Antarctica and over Greenland is also illustrated. We found that the GAE/GAF-corrected trend of the global atmospheric mass over the GRACE mission lifetime significantly decreased from –0.05 to –0.02 mm yr –1 in terms of geoid height. A considerable effect (33 per cent) was also found in the quadratic term of ice mass loss over Antarctica which results in an acceleration of 3.2 Gt yr –1  yr –1 smaller than without applying this correction.

Description: Previous studies of Earth rotation perturbations due to ice-age loading have predicted a slow secular drift of the rotation axis relative to the surface geography (i.e. true polar wander, TPW) of order of several degrees over the Plio-Pleistocene. It has been argued that this drift and the change in the geographic distribution of solar insolation that it implies may have been responsible for important transitions in ice-age climate, including the termination of ice-age cycles.We use a revised rotational stability theory that incorporates a more accurate treatment of the Earth's background ellipticity to reconsider this issue, and demonstrate that the net displacement of the pole predicted in earlier studies disappears. This more muted polar motion is due to two factors: first, the revised theory no longer predicts the permanent shift in the rotation axis, or the so-called ‘unidirectional TPW’, that appears in the traditional stability theory; and, second, the increased background ellipticity incorporated in the revised predictions acts to reduce the normal mode amplitudes governing the motion of the pole. We conclude that ice-age-induced TPW was not responsible for the termination of the ice age. This does not preclude the possibility that TPW induced by mantle convective flow may have played a role in major Plio-Pleistocene climate transitions, including the onset of Northern Hemisphere glaciation.

Description: Peatlands of all latitudes play an integral role in global climate change by serving as a carbon sink and a primary source of atmospheric methane; however, the microbial ecology of mid-latitude peatlands is vastly understudied. Herein, next generation Illumina amplicon sequencing of small subunit rRNA genes was utilized to elucidate the microbial communities in three southern Appalachian peatlands. In contrast to northern peatlands, Proteobacteria dominated over Acidobacteria in all three sites. An average of 11 bacterial phyla was detected at relative abundance values 〉1%, with three candidate divisions (OP3, WS3 and NC10) represented, indicating high phylogenetic diversity. Physiological traits of isolates within the candidate alphaproteobacterial order, Ellin 329, obtained here and in previous studies indicate that bacteria of this order may be involved in hydrolysis of poly-, di- and monosaccharides. Community analyses indicate that Ellin 329 is the third most abundant order and is most abundant near the surface layers where plant litter decomposition should be primarily occurring. In sum, members of Ellin 329 likely play important roles in organic matter decomposition, in southern Appalachian peatlands and should be investigated further in other peatlands and ecosystem types.

Description: Marine viruses are the most abundant biological entity in the oceans, the majority of which infect bacteria and are known as bacteriophages. Yet, the bulk of bacteriophages form part of the vast uncultured dark matter of the microbial biosphere. In spite of the paucity of cultured marine bacteriophages, it is known that marine bacteriophages have major impacts on microbial population structure and the biogeochemical cycling of key elements. Despite the ecological relevance of marine bacteriophages, there are relatively few isolates with complete genome sequences. This minireview focuses on knowledge gathered from these genomes put in the context of viral metagenomic data and highlights key advances in the field, particularly focusing on genome structure and auxiliary metabolic genes.

Description: The fynbos biome in South Africa is globally recognised as a plant biodiversity hotspot. However, very little is known about the bacterial communities associated with fynbos plants, despite interactions between primary producers and bacteria having an impact on the physiology of both partners and shaping ecosystem diversity. This study reports on the structure, phylogenetic composition and potential roles of the endophytic bacterial communities located in the stems of three fynbos plants ( Erepsia anceps , Phaenocoma prolifera and Leucadendron laureolum ). Using Illumina MiSeq 16S rRNA sequencing we found that different subpopulations of Deinococcus-Thermus, Alphaproteobacteria, Acidobacteria and Firmicutes dominated the endophytic bacterial communities. Alphaproteobacteria and Actinobacteria were prevalent in P. prolifera , whereas Deinococcus-Thermus dominated in L. laureolum , revealing species-specific host–bacteria associations. Although a high degree of variability in the endophytic bacterial communities within hosts was observed, we also detected a core microbiome across the stems of the three plant species, which accounted for 72% of the sequences. Altogether, it seems that both deterministic and stochastic processes shaped microbial communities. Endophytic bacterial communities harboured putative plant growth-promoting bacteria, thus having the potential to influence host health and growth.

Description: Apparent acceleration in Gravity Recovery and Climate Experiment (GRACE) Antarctic ice mass time-series may reflect both ice discharge and surface mass balance contributions. However, a recent study suggests there is also contamination from errors in atmospheric pressure de-aliasing fields [European Center for Medium-Range Weather Forecast (ECMWF) operational products] used during GRACE data processing. To further examine this question, we compare GRACE atmospheric pressure de-aliasing (GAA) fields with in situ surface pressure data from coastal and inland stations. Differences between the two are likely due to GAA errors, and provide a measure of error in GRACE solutions. Time-series of differences at individual weather stations are fit to four presumed error components: annual sinusoids, a linear trend, an acceleration term and jumps at times of known ECMWF model changes. Using data from inland stations, we estimate that atmospheric pressure error causes an acceleration error of about +7.0 Gt yr –2 , which is large relative to prior GRACE estimates of Antarctic ice mass acceleration in the range of –12 to –14 Gt yr –2 . We also estimate apparent acceleration rates from other barometric pressure (reanalysis) fields, including ERA-Interim, MERRA and NCEP/DOE. When integrated over East Antarctica, the four mass acceleration estimates (from GAA and the three reanalysis fields) vary considerably (by ~2–16 Gt yr –2 ). This shows the need for further effort to improve atmospheric mass estimates in this region of sparse in situ observations, in order to use GRACE observations to measure ice mass acceleration and related sea level change.

Description: The long-wavelength gravity field contains information about processes in the sublithospheric mantle. As satellite-derived gravity models now provide the long to medium-wavelength gravity field at unprecedented accuracy, techniques used to process gravity data need to be updated. We show that when determining these long-wavelengths, the treatment of topographic-isostatic effect (TIE) and isostatic effects (IE) is a likely source of error. We constructed a global isostatic model and calculated global TIE and IE. These calculations were done for ground stations as well as stations at satellite height. We considered both gravity and gravity gradients. Using these results, we determined how much of the gravity signal comes from distant sources. We find that a significant long-wavelength bias is introduced if far-field effects on the topographic effect are neglected. However, due to isostatic compensation far-field effects of the topographic effect are to a large degree compensated by the far-field IE. This means that far-field effects can be reduced effectively by always considering topographic masses together with their compensating isostatic masses. We show that to correctly represent the ultra-long wavelengths, a global background model should be used. This is demonstrated both globally and for a continental-scale case area in North America. In the case of regional modelling, where the ultra-long wavelengths are not of prime importance, gravity gradients can be used to help minimize correction errors caused by far-field effects.

Description: The functioning of many natural and engineered environments is dependent on long distance electron transfer mediated through electrical currents. These currents have been observed in exoelectrogenic biofilms and it has been proposed that microbial biofilms can mediate electron transfer via electrical currents on the centimeter scale. However, direct evidence to confirm this hypothesis has not been demonstrated and the longest known electrical transfer distance for single species exoelectrogenic biofilms is limited to 100 μm. In the present study, biofilms were developed on electrodes with electrically non-conductive gaps from 50 μm to 1 mm and the in situ conductance of biofilms was evaluated over time. Results demonstrated that the exoelectrogenic mixed species biofilms in the present study possess the ability to transfer electrons through electrical currents over a distance of up to 1 mm, 10 times further than previously observed. Results indicate the possibility of interspecies interactions playing an important role in the spatial development of exoelectrogenic biofilms, suggesting that these biological networks might remain conductive even at longer distance. These findings have significant implications in regards to future optimization of microbial electrochemical systems.

Description: Constraining laterally varying structures in planetary interiors is important for understanding both the composition and the internal dynamics of a planet. Recognizing that seismic imaging technique is currently only viable for studying the Earth's interior structures, methods that can be supported by advanced space geodetic techniques may become alternatives to ‘image’ the interiors of other planets. The method of tidal tomography is one possibility, and it relies on high precision measurement of the response of a planet to its body tide. However, it is essential to develop an efficient analytical tool that computes the dependence of tidal response to 3-D interior structures. In this paper, we present a complete formulation of such an analytical tool, which calculates to high accuracy the tidal response of a terrestrial planet with lateral heterogeneities in its elastic and density structures. We treat the lateral heterogeneities as small perturbations and derive the governing equations based on the perturbation theory. In a spherical harmonic representation, equations at each order of perturbation are reduced into multiple matrix equations at harmonics that are allowed by mode couplings, and the total response equals the sum of all those single-harmonic responses, which can be solved semi-analytically. We test our perturbation method by applying it to the Moon with a harmonic degree-1 mantle structure for which the perturbation solutions of the tidal response are compared with those from a fully numerical method. The remarkable agreement between results from these two methods validates the perturbation method. As an example, we then use the perturbation method to evaluate the impact of lunar crustal thickness variations on tidal response of the Moon. We find that lunar crust produces much smaller degree-3 tidal responses than a relatively weak degree-1 structure in the deep lunar mantle. Our calculations show that degree-3 tidal response measurements may hold key constraints on possible degree-1 mantle structure of the Moon, as suggested from previous modelling results.

Description: Seismic data are primarily used in studies of the Earth's lithospheric structure including the Moho geometry. In regions, where seismic data are sparse or completely absent, gravimetric or combined gravimetric-seismic methods could be applied to determine the Moho depth. In this study, we derive and present generalized expressions for solving the Vening Meinesz–Moritz's (VMM) inverse problem of isostasy for a Moho depth determination from gravity and vertical gravity-gradient data. By solving the (non-linear) Fredholm's integral equation of the first kind, the linearized observation equations, which functionally relate the (given) gravity/gravity-gradient data to the (unknown) Moho depth, are derived in the spectral domain. The VMM gravimetric results are validated by using available seismic and gravimetric Moho models. Our results show that the VMM Moho solutions obtained by solving the VMM problem for gravity and gravity-gradient data are almost the same. This finding indicates that in global applications, using the global gravity/gravity-gradient data coverage, the spherical harmonic expressions for the gravimetric forward and inverse modelling yield (theoretically) the same results. Globally, these gravimetric solutions have also a relatively good agreement with the CRUST1.0 and GEMMA GOCE models in terms of their rms Moho differences (4.7 km and 4.1 km, respectively).

Description: Geophysical techniques are widely used to monitor volcanic unrest. A number of studies have also demonstrated that hydrological processes can produce or trigger geophysical signals. Hydrologically induced gravity signals have previously been recorded by specifically designed gravity surveys as well as, inadvertently, by volcano monitoring studies. Water table corrections of microgravity surveys are commonplace. However, the fluctuations of the water table beneath survey locations are often poorly known, and such a correction fails to account for changes in water-mass storage in the unsaturated zone. Here, we combine 2-D axis-symmetrical numerical fluid-flow models with an axis-symmetric, distributed-mass, gravity calculation to model gravity changes in response to fluctuating hydrological recharge. Flow simulations are based on tropical volcanic settings where high surface permeabilities promote thick unsaturated zones. Our study highlights that mass storage (saturation) changes within the unsaturated zone beneath a survey point can generate recordable gravity changes. We show that for a tropical climate, recharge variations can generate gravity variations of over 150 μGal; although, we demonstrate that for the scenarios investigated here, the probability of recording such large signals is low. Our modelling results indicate that microgravity survey corrections based on water table elevation may result in errors of up to 100 μGal. The effect of inter-annual recharge fluctuations dominate over seasonal cycles which makes prediction and correction of the hydrological contribution more difficult. Spatial hydrogeological heterogeneity can also impact on the accuracy of relative gravity surveys, and can even result in the introduction of additional survey errors. The loading fluctuations associated with saturation variations in the unsaturated zone may also have implications for other geophysical monitoring techniques, such as geodetic monitoring of ground deformation.

Description: Bacteriophages are increasingly being used as water quality indicators. Two groups of phages infecting Escherichia coli , somatic and F-specific coliphages, are being considered as indicators of fecal and viral contamination for several types of water around the world. However, some uncertainties remain regarding which coliphages to assess. Recently, E. coli strain CB390 has been reported to be suitable for simultaneous detection of both groups, which seems to be more informative than determining only one of the groups. Here, a significant number of samples from different settings, mostly those where F-specific phages have been reported to outnumber somatic coliphages, are analyzed for somatic coliphages, F-specific RNA phages by standardized methods and coliphages detected by host strain CB390. The results presented here confirm that the numbers of phages counted using CB390 are equivalent to the sum of the somatic and F-specific coliphages counted independently in all settings. Hence the usefulness of this strain for simultaneous detection of somatic and F-specific coliphages is confirmed. Also, sets of data on the presence of coliphages in reclaimed and groundwater are reported.

Description: A new approach based on energy conservation principle for satellite gravimetry mission has been developed and yields more accurate estimation of in situ geopotential difference observables using K-band ranging (KBR) measurements from the Gravity Recovery and Climate Experiment (GRACE) twin-satellite mission. This new approach preserves more gravity information sensed by KBR range-rate measurements and reduces orbit error as compared to previous energy balance methods. Results from analysis of 11 yr of GRACE data indicated that the resulting geopotential difference estimates agree well with predicted values from official Level 2 solutions: with much higher correlation at 0.9, as compared to 0.5–0.8 reported by previous published energy balance studies. We demonstrate that our approach produced a comparable time-variable gravity solution with the Level 2 solutions. The regional GRACE temporal gravity solutions over Greenland reveals that a substantially higher temporal resolution is achievable at 10-d sampling as compared to the official monthly solutions, but without the compromise of spatial resolution, nor the need to use regularization or post-processing.

Description: Several attempts have been made to obtain a radiographic image inside volcanoes using cosmic-ray muons (muography). Muography is expected to resolve highly heterogeneous density profiles near the surface of volcanoes. However, several prior works have failed to make clear observations due to contamination by background noise. The background contamination leads to an overestimation of the muon flux and consequently a significant underestimation of the density in the target mountains. To investigate the origin of the background noise, we performed a Monte Carlo simulation. The main components of the background noise in muography are found to be low-energy protons, electrons and muons in case of detectors without particle identification and with energy thresholds below 1 GeV. This result was confirmed by comparisons with actual observations of nuclear emulsions. This result will be useful for detector design in future works, and in addition some previous works of muography should be reviewed from the view point of background contamination.

Description: Essential to understanding sea level change and its causes during the last interglacial (LIG) is the quantification of uncertainties. In order to estimate the uncertainties, we develop a statistical framework for the comparison of palaeoclimatic sea level index points and GIA model predictions. For the investigation of uncertainties, as well as to generate better model predictions, we implement a massive ensemble approach by applying a data assimilation scheme based on particle filter methods. The different runs are distinguished through varying ice sheet reconstructions based on oxygen-isotope curves and different parameter selections within the GIA model. This framework has several advantages over earlier work, such as the ability to examine either the contribution of individual observations to the results or the probability of specific input parameters. This exploration of input parameters and data leads to a larger range of estimates than previously published work. We illustrate how the assumptions that enter into the statistical analysis, such as the existence of outliers in the observational database or the initial ice volume history, can introduce large variations to the estimate of the maximum highstand. Thus, caution is required to avoid overinterpreting results. We conclude that there are reasonable doubts whether the data sets previously used in statistical analyses are able to tightly constrain the value of maximum highstand during the LIG.

Description: A 3-D density model of the crust and upper mantle beneath the Karoo basin is presented here. The model is constrained using potential field, borehole and seismic data. Uplift of the basin by the end of the Cretaceous has resulted in an unusually high plateau (〉1000 m) covering a large portion of South Africa. Isostatic studies show the topography is largely compensated by changes in Moho depths (~35 km on-craton and 〉45 km off-craton) and changes in lithospheric mantle densities between the Kaapvaal Craton and surrounding regions (~50 kg m –3 increase from on- to off-craton). This density contrast is determined by inverted satellite gravity and gravity gradient data. The highest topography along the edge of the plateau (〉1200 m) and a strong Bouguer gravity low over Lesotho, however, can only be explained by a buoyant asthenosphere with a density decrease of around 40 kg m –3 .

Description: It is well known that Methylosinus trichosporium OB3b has two forms of methane monooxygenase (MMO) responsible for the initial conversion of methane to methanol, a cytoplasmic (soluble) methane monooxygenase and a membrane-associated (particulate) methane monooxygenase, and that copper strongly regulates expression of these alternative forms of MMO. More recently, it has been discovered that M. trichosporium OB3b has multiple types of the methanol dehydrogenase (MeDH), i.e. the Mxa-type MeDH (Mxa-MeDH) and Xox-type MeDH (Xox-MeDH), and the expression of these two forms is regulated by the availability of the rare earth element (REE), cerium. Here, we extend these studies and show that lanthanum, praseodymium, neodymium and samarium also regulate expression of alternative forms of MeDH. The effect of these REEs on MeDH expression, however, was only observed in the absence of copper. Further, a mutant of M. trichosporium OB3b, where the Mxa-MeDH was knocked out, was able to grow in the presence of lanthanum, praseodymium and neodymium, but was not able to grow in the presence of samarium. Collectively, these data suggest that multiple levels of gene regulation by metals exist in M. trichosporium OB3b, but that copper overrides the effect of other metals by an as yet unknown mechanism.

Description: A method to estimate the rotation change in the orientation of the centre-of-figure (CF) frame caused by earthquakes is proposed for the first time. This method involves using the point dislocation theory based on a spherical, non-rotating, perfectly elastic and isotropic (SNREI) Earth. The rotation change in the orientation is related solely to the toroidal displacements of degree one induced by the vertical dip slip dislocation, and the spheroidal displacements induced by an earthquake have no contribution. The effects of two recent large earthquakes, the 2004 Sumatra and the 2011 Tohoku-Oki, are studied. Results showed that the Sumatra and Tohoku-Oki earthquakes both caused the CF frame to rotate by at least tens of μas (micro-arc-second). Although the visible co-seismic displacements are identified and removed from the coordinate time-series, the rotation change due to the unidentified ones and errors in removal is non-negligible. Therefore, the rotation change in the orientation of the CF frame due to seismic deformation should be taken into account in the future in reference frame and geodesy applications.

Description: In this study, we propose an approach for determining the geopotential difference using high-frequency-stability microwave links between satellite and ground station based on Doppler cancellation system. Suppose a satellite and a ground station are equipped with precise optical-atomic clocks (OACs) and oscillators. The ground oscillator emits a signal with frequency f a towards the satellite and the satellite receiver (connected with the satellite oscillator) receives this signal with frequency f b which contains the gravitational frequency shift effect and other signals and noises. After receiving this signal, the satellite oscillator transmits and emits, respectively, two signals with frequencies f b and f c towards the ground station. Via Doppler cancellation technique, the geopotential difference between the satellite and the ground station can be determined based on gravitational frequency shift equation by a combination of these three frequencies. For arbitrary two stations on ground, based on similar procedures as described above, we may determine the geopotential difference between these two stations via a satellite. Our analysis shows that the accuracy can reach 1 m 2 s – 2 based on the clocks’ inaccuracy of about 10 –17 (s s –1 ) level. Since OACs with instability around 10 –18 in several hours and inaccuracy around 10 –18 level have been generated in laboratory, the proposed approach may have prospective applications in geoscience, and especially, based on this approach a unified world height system could be realized with one-centimetre level accuracy in the near future.

Description: Polysulfides (S x 2– ) are sulfide oxidation intermediates that are important for a variety of environmentally relevant processes including pyrite formation, organic matter sulfidization, isotope exchange among reduced sulfur species, and metal chelation. In addition to their chemical reactivity, laboratory experiments with microbial cultures and enzymes indicate both indirect and direct roles for microorganisms in affecting polysulfide chemistry in natural environments through production and consumption. As polysulfides have been detected in a wide array of natural systems ranging from microbial mats to hydrothermal vents, constraining their biogeochemical cycling has broad impacts. However, many questions remain regarding the processes responsible for polysulfide dynamics in these environments and the precise role that microorganisms play in these processes. This review provides a summary of laboratory experiments investigating the role of polysulfides in microbial metabolism, and observations of polysulfides in the environment in order to provide further insight into and highlight open questions about this significant component of the sulfur cycle.

Description: We interpreted the TRIDENT satellite derived gravity field to provide detailed insights into the spatial distribution of the crustal density structures in the area of the Yellow Sea. We used 3-D forward density modelling for the interpretation that incorporated constraints from existing geological and geophysical information. A gravity stripping method is used to separate out the gravity effects of different geological crustal structures. From this analysis we see that (1) the Gunsan sedimentary basin is isostatically compensated. (2) The satellite-derived Bouguer anomalies ranging from 15 to –30 x 10 –5 m s –2 are linked to basin thicknesses in the Yellow Sea. (3) The calculated Moho depth in the Yellow Sea varies from 27 km beneath the deep sedimentary basin to 34 km in the uplifted zones. (4) Moho depth calculations show two distinct areas, characterized by the deepest Moho depths and the largest crustal thicknesses in the Yellow Sea. The one region extends along the Qianliyan Uplift Zone from Jiaodong to Hongsung while the other area extends from southeastern China to Hongsung in the Korean peninsula. Compared to previous works we suggest that they are the part of the collisions zone between North and South China Blocks extending from China to the Korean peninsula via the Yellow Sea.

Description: Regional refinement of the gravity field models from satellite data using spherical radial base functions (SRBF) is an ill-posed problem. This is mainly due to the regional confinement of the data and the base functions, which leads to severe instabilities in the solutions. Here, this ill-posedness as well as the related regularization process are investigated. We compare three methods for the choice of the regularization parameter, which have been frequently used in gravity modelling. These methods are (1) the variance component estimation (VCE), (2) the generalized cross validation (GCV) and (3) the L-curve criterion. A particular emphasis is put on the impact of the SRBF type on the regularization parameter. To do this, we include two types of SRBF which are often used for regional gravity field modelling. These are the Shannon SRBF or the reproducing kernel and the Spline SRBF. The investigations are performed on two months of the real GOCE ultrasensitive gravity gradients over Central Africa and Amazon. The solutions are validated against a state-of-the-art global gravity solution. We conclude that if a proper regularization method is applied, both SRBF deliver more or less the same accuracy. We show that when the Shannon wavelet is used, the L-curve method gives the best results, while with the Spline kernel, the GCV outperforms the other two methods. Moreover, we observe that the estimated coefficients for the Spline kernel cannot be spatially interpreted. In contrast, the coefficients obtained for the Shannon wavelet reflect the energy of the recovered gravity field with a correlation factor of above 95 per cent. Therefore, when combined with the L-curve method, the Shannon SRBF is advantageous for regional gravity field estimation, since it is one of the simplest band-limited SRBF. In addition, it delivers promising solutions and the estimated coefficients represent the characteristics of the gravity field within the target region.

Description: This paper describes an alternative acceleration approach for determining GRACE monthly gravity field models. The main differences compared to the traditional acceleration approach can be summarized as: (1) The position errors of GRACE orbits in the functional model are taken into account; (2) The range ambiguity is eliminated via the difference of the range measurements and (3) The mean acceleration equation is formed based on Cowell integration. Using this developed approach, a new time-series of GRACE monthly solution spanning the period January 2003 to December 2010, called Tongji_Acc RL01, has been derived. The annual signals from the Tongji_Acc RL01 time-series agree well with those from the GLDAS model. The performance of Tongji_Acc RL01 shows that this new model is comparable with the RL05 models released by CSR and JPL as well as with the RL05a model released by GFZ.

Description: In this paper, we present a method for incorporating prior geological information into potential field data inversion problem. As opposed to the traditional inverse algorithm, our proposed method takes full advantage of prior geological information as a constraint and thus obtains a new objective function for inversion by adding Lagrangian multipliers and slack variables to the traditional inversion method. These additional parameters can be easily solved during iterations. We used both synthetic and observed data sets to test the stability and validity of the proposed method. Our results using synthetic gravity data show that our new method predicts depth and density anomalies more efficiently and accurately than the traditional inversion method that does not include prior geological constraints. Then using observed gravity data in the Three Gorges area and geological constraint information, we obtained the density distribution of the upper and middle crust in this area thus revealing its geological structure. These results confirm the proposed method's validity and indicate its potential application for magnetism data inversion and exploration of geological structures.

Description: It remains enigmatic how slow slip events (SSEs) interact with other slow seismic events and large distant earthquakes at many subduction zones. Here we model the spatiotemporal slip evolution of the most recent long-term SSE in 2009–2011 in the Bungo Channel region, southwest Japan using GEONET GPS position time-series and a Kalman filter-based, time-dependent slip inversion method. We examine the space-time relationship between the geodetically determined slow slip transient and seismically observed low frequency earthquakes (LFEs) and very-low frequency earthquakes (V-LFEs) near the Nankai trough. We find a strong but distinct temporal correlation between transient slip and LFEs and V-LFEs, suggesting a different relationship to the SSE. We also find the great Tohoku-Oki earthquake appears to disrupt the normal source process of the SSE, probably reflecting large-scale stress redistribution caused by the earthquake. Comparison of the 2009–2011 SSE with others in the same region shows much similarity in slip and moment release, confirming its recurrent nature. Comparison of transient slip with plate coupling shows that slip transients mainly concentrate on the transition zone from strong coupling region to downdip LFEs with transient slip relieving elastic strain accumulation at transitional depth. The less consistent spatial correlation between the long-term SSE and seismic slow earthquakes, and susceptibility of these slow earthquakes to various triggering sources including long-term slow slip, suggests caution in using the seismically determined slow earthquakes as a proxy for slow slip.

Description: Normal mode treatments of the Earth's body tide response were developed in the 1980s to account for the effects of Earth rotation, ellipticity, anelasticity and resonant excitation within the diurnal band. Recent space-geodetic measurements of the Earth's crustal displacement in response to luni-solar tidal forcings have revealed geographical variations that are indicative of aspherical deep mantle structure, thus providing a novel data set for constraining deep mantle elastic and density structure. In light of this, we make use of advances in seismic free oscillation literature to develop a new, generalized normal mode theory for the tidal response within the semi-diurnal and long-period tidal band. Our theory involves a perturbation method that permits an efficient calculation of the impact of aspherical structure on the tidal response. In addition, we introduce a normal mode treatment of anelasticity that is distinct from both earlier work in body tides and the approach adopted in free oscillation seismology. We present several simple numerical applications of the new theory. First, we compute the tidal response of a spherically symmetric, non-rotating, elastic and isotropic Earth model and demonstrate that our predictions match those based on standard Love number theory. Second, we compute perturbations to this response associated with mantle anelasticity and demonstrate that the usual set of seismic modes adopted for this purpose must be augmented by a family of relaxation modes to accurately capture the full effect of anelasticity on the body tide response. Finally, we explore aspherical effects including rotation and we benchmark results from several illustrative case studies of aspherical Earth structure against independent finite-volume numerical calculations of the semi-diurnal body tide response. These tests confirm the accuracy of the normal mode methodology to at least the level of numerical error in the finite-volume predictions. They also demonstrate that full coupling of normal modes, rather than group coupling, is necessary for accurate predictions of the body tide response.

Description: A total of 65 spore-forming mercury-resistant bacteria were isolated from natural environments worldwide in order to understand the acquisition of additional genes by and dissemination of mercury resistance transposons across related Bacilli genera by horizontal gene movement. PCR amplification using a single primer complementary to the inverted repeat sequence of Tn MERI1 -like transposons showed that 12 of 65 isolates had a transposon-like structure. There were four types of amplified fragments: Tn 5084 , Tn 5085 , Tn d MER3 (a newly identified deleted transposon-like fragment) and Tn 6294 (a newly identified transposon). Tn d MER3 is a 3.5-kb sequence that carries a merRETPA operon with no merB or transposase genes. It is related to the mer operon of Bacillus licheniformis strain FA6-12 from Russia. DNA homology analysis shows that Tn 6294 is an 8.5-kb sequence that is possibly derived from Tn d MER3 by integration of a Tn MERI1 -type transposase and resolvase genes and in addition the merR2 and merB1 genes. Bacteria harboring Tn 6294 exhibited broad-spectrum mercury resistance to organomercurial compounds, although Tn 6294 had only merB1 and did not have the merB2 and merB3 sequences for organomercurial lyases found in Tn 5084 of B. cereus strain RC607. Strains with Tn 6294 encode mercuric reductase (MerA) of less than 600 amino acids in length with a single N-terminal mercury-binding domain, whereas MerA encoded by strains MB1 and RC607 has two tandem domains. Thus, Tn d MER3 and Tn 6294 are shorter prototypes for Tn MERI1 -like transposons. Identification of Tn 6294 in Bacillus sp. from Taiwan and in Paenibacillus sp. from Antarctica indicates the wide horizontal dissemination of Tn MERI1 -like transposons across bacterial species and geographical barriers.

Description: A sequence of large earthquakes occurred along the North Anatolian fault in the 20th century. These earthquakes, including the 1999 Izmit/Düzce earthquakes, generally propagated westward towards the Marmara Sea, defining the Main Marmara fault as a potential seismic gap. It is important to conduct a detailed assessment of the seismic hazards along the main Marmara fault because the megacity Istanbul lies only approximately 10 km north of the eastern segment of the Main Marmara fault, which is referred to as the Princes’ Islands Fault segment (PIF). Here, we study the locking status of this fault segment to evaluate the seismic hazard potential. For the first time, combined ascending and descending Interferometric Synthetic Aperture Radar and Global Positioning System observations were used to investigate the crustal deformation associated with the PIF. After careful corrections of the estimated ground velocity, a deformation pattern relating to fault locking near the Princes’ Islands was identified. The modeling results revealed that the slip rate and locking depth of the fault segment show a clear trade-off, which were estimated as 18.9 ± 7.2 mm yr –1 and 12.1 ± 7.0 km, respectively. With a moment accumulation rate of 1.7 ± 0.4  x  10 17 Nm yr –1 (proportional to the product of slip rate and locking depth), our results imply a build-up of a geodetic moment on the PIF and therefore a potential for earthquake hazards in the vicinity of the Istanbul megacity.

Description: Fungi may play an important role in the production of the greenhouse gas nitrous oxide (N 2 O). Bipolaris sorokiniana is a ubiquitous saprobe found in soils worldwide, yet denitrification by this fungal strain has not previously been reported. We aimed to test if B. sorokiniana would produce N 2 O and CO 2 in the presence of organic and inorganic forms of nitrogen (N) under microaerobic and anaerobic conditions. Nitrogen source (organic-N, inorganic-N, no-N control) significantly affected N 2 O and CO 2 production both in the presence and absence of oxygen, which contrasts with bacterial denitrification. Inorganic N addition increased denitrification of N 2 O (from 0 to 0.3 μg N 2 0-N h –1  g –1 biomass) and reduced respiration of CO 2 (from 0.1 to 0.02 mg CO 2 h –1  g –1 biomass). Isotope analyses indicated that nitrite, rather than ammonium or glutamine, was transformed to N 2 O. Results suggest the source of N may play a larger role in fungal N 2 O production than oxygen status.

Description: Legionella pneumophila is a pathogenic bacterium commonly found in water and responsible for severe pneumonia. Free-living amoebae are protozoa also found in water, which feed on bacteria by phagocytosis. Under favorable conditions, some L. pneumophila are able to resist phagocytic digestion and even multiply within amoebae. However, it is not clear whether L. pneumophila could infect at a same rate a large range of amoebae or if there is some selectivity towards specific amoebal genera or strains. Also, most studies have been performed using collection strains and not with freshly isolated strains. In our study, we assess the permissiveness of freshly isolated environmental strains of amoebae, belonging to three common genera (i.e. Acanthamoeba, Naegleria and Vermamoeba ), for growth of L. pneumophila at three different temperatures. Our results indicated that all the tested strains of amoebae were permissive to L. pneumophila Lens and that there was no significant difference between the strains. Intracellular proliferation was more efficient at a temperature of 40°C. In conclusion, our work suggests that, under favorable conditions, virulent strains of L. pneumophila could equally infect a large number of isolates of common freshwater amoeba genera.

Description: We present a distributed slip model for the 1999 M w 6.3 Chamoli earthquake of north India using interferometric synthetic aperture radar (InSAR) data from both ascending and descending orbits and Bayesian estimation of confidence levels and trade-offs of the model geometry parameters. The results of fault-slip inversion in an elastic half-space show that the earthquake ruptured a $9 _{ - 2.2}^{\circ + 3.4}$ northeast-dipping plane with a maximum slip of ~1 m. The fault plane is located at a depth of ~ $15.9_{ - 3.0}^{ + 1.1}$ km and is ~120 km north of the Main Frontal Thrust, implying that the rupture plane was on the northernmost detachment near the mid-crustal ramp of the Main Himalayan Thrust. The InSAR-determined moment is 3.35 x 10 18 Nm with a shear modulus of 30 GPa, equivalent to M w 6.3, which is smaller than the seismic moment estimates of M w 6.4–6.6. Possible reasons for this discrepancy include the trade-off between moment and depth, uncertainties in seismic moment tensor components for shallow dip-slip earthquakes and the role of earth structure models in the inversions. The released seismic energy from recent earthquakes in the Garhwal region is far less than the accumulated strain energy since the 1803 M s 7.5 earthquake, implying substantial hazard of future great earthquakes.

Description: We use a Bayesian formalism combined with a grid node discretization for the linear inversion of gravimetric data in terms of 3-D density distribution. The forward modelling and the inversion method are derived from seismological inversion techniques in order to facilitate joint inversion or interpretation of density and seismic velocity models. The Bayesian formulation introduces covariance matrices on model parameters to regularize the ill-posed problem and reduce the non-uniqueness of the solution. This formalism favours smooth solutions and allows us to specify a spatial correlation length and to perform inversions at multiple scales. We also extract resolution parameters from the resolution matrix to discuss how well our density models are resolved. This method is applied to the inversion of data from the volcanic island of Basse-Terre in Guadeloupe, Lesser Antilles. A series of synthetic tests are performed to investigate advantages and limitations of the methodology in this context. This study results in the first 3-D density models of the island of Basse-Terre for which we identify: (i) a southward decrease of densities parallel to the migration of volcanic activity within the island, (ii) three dense anomalies beneath Petite Plaine Valley, Beaugendre Valley and the Grande-Découverte-Carmichaël-Soufrière Complex that may reflect the trace of former major volcanic feeding systems, (iii) shallow low-density anomalies in the southern part of Basse-Terre, especially around La Soufrière active volcano, Piton de Bouillante edifice and along the western coast, reflecting the presence of hydrothermal systems and fractured and altered rocks.

Description: We compute the gravimetric factor at the Chandler wobble (CW) frequency using time-series from superconducting gravimeters (SG) longer than a decade. We first individually process the polar motion and data at each individual gravity station to estimate the gravimetric factor amplitude and phase, then we make a global analysis by applying a stacking method to different subsets of up to seven SG stations. The stacking is an efficient way of getting rid of local effects and improving the signal-to-noise ratio of the combined data sets. Using the stacking method, we find a gravimetric factor amplitude and phase of 1.118 ± 0.016 and –0.45 ± 0.66 deg, respectively, which is smaller in amplitude than expected. The sources of error are then carefully considered. For both local and global analyses, the uncertainties on our results are reliably constrained by computing the standard deviation of the estimates of the gravimetric factor amplitude and phase for increasing length of the time-series. Constraints on the CW anelastic dissipation can be set since any departure of the gravimetric factor from its elastic value may provide some insights into the dissipative processes that occur at the CW period. In particular, assuming given rheological models for the Earth's mantle enables us to make the link between the gravimetric factor phase and the CW quality factor.

Description: Analysing independent 1-yr data sets of 10 European superconducting gravimeters (SG) reveals statistically significant temporal variations of M2 tidal parameters. Both common short-term (〈2 yr) and long-term (〉2 yr) features are identified in all SG time-series but one. The averaged variations of the amplitude factor are about 0.2. The path of load vector variations equivalent to the temporal changes of tidal parameters suggests the presence of an 8.85 yr modulation (lunar perigee). The tidal waves having the potential to modulate M2 with this period belong to the 3rd degree constituents. Their amplitude factors turn out to be much closer to body tide model predictions than that of the main 2nd degree M2, which indicates ocean loading for 3rd degree waves to be less prominent than for 2nd degree waves within the M2 group. These two different responses to the loading suggest that the observed modulation is more due to insufficient frequency resolution of limited time-series rather than to time variable loading. Presently, SG gravity time-series are still too short to prove if time variable loading processes are involved too as in case of the annual M2 modulation known to appear for analysis intervals of less than 1 yr. Whatever the variations are caused by, they provide the upper accuracy limit for earth model validation and permit estimating the temporal stability of SG scale factors and assessing the quality of gravity time-series.

Description: A sequence of large earthquakes occurred along the North Anatolian fault in the 20th century. These earthquakes, including the 1999 Izmit/Düzce earthquakes, generally propagated westward towards the Marmara Sea, defining the Main Marmara fault as a potential seismic gap. It is important to conduct a detailed assessment of the seismic hazards along the main Marmara fault because the megacity Istanbul lies only approximately 10 km north of the eastern segment of the Main Marmara fault, which is referred to as the Princes’ Islands Fault segment (PIF). Here, we study the locking status of this fault segment to evaluate the seismic hazard potential. For the first time, combined ascending and descending Interferometric Synthetic Aperture Radar and Global Positioning System observations were used to investigate the crustal deformation associated with the PIF. After careful corrections of the estimated ground velocity, a deformation pattern relating to fault locking near the Princes’ Islands was identified. The modeling results revealed that the slip rate and locking depth of the fault segment show a clear trade-off, which were estimated as 18.9 ± 7.2 mm yr –1 and 12.1 ± 7.0 km, respectively. With a moment accumulation rate of 1.7 ± 0.4  x  10 17 Nm yr –1 (proportional to the product of slip rate and locking depth), our results imply a build-up of a geodetic moment on the PIF and therefore a potential for earthquake hazards in the vicinity of the Istanbul megacity.

Description: Deformation analysis in general and strain analysis in particular using permanent GPS networks require proper analysis of time-series in which all functional effects are taken into consideration and all stochastic effects are captured using an appropriate noise model. This contribution addresses both issues when considering the strain parameters of a GPS network. Estimates of spatial correlation, time correlated noise, and multivariate power spectrum for daily position time-series of the Southern California Integrated GPS Network (SCIGN) stations collected between 1996 and 2011 are obtained. Significant signals with periods of 13.63 d and those related to the GPS draconitic year are identified in these time-series. We aim to assess the effect of a realistic noise model of the series on the uncertainties of the strain parameters including displacements, normal and shear strains, and rotations. For the SCIGN network considered, the following results are highlighted. Contrary to the common belief, the uncertainties of the displacements parameters become smaller when taking a realistic noise model into account. This however was not the case when assessing the noise characteristics of the normal and shear strain, and rotation parameters. The uncertainties increase nearly by a factor of two, in agreement to what is expected. Some of the significant deformation parameters of the white noise model become less significant in case of the realistic noise model.

Description: Crustal vertical deformation (CVD) observed by continuous GPS height time-series can be explained largely by surface loading effects recovered from both Gravity Recover and Climate Experiment (GRACE) and General Circulation Models (GCMs) data. We first show that lower degree CVD spatial spectrum due to the Earth's elastic response to a uniform surface loading plays more important roles than that of high-degree case. We then demonstrate that GRACE data with 300–400 km spatial resolution have the ability to detect 99 per cent power of global and regional CVD in spatial spectrum domain using a global frequency–wavenumber spectrum method. We can just use either GRACE or GCMs 36 degree/order (d/o) spherical harmonic coefficients (SHCs) which correspond to 500 km spatial resolution to acquire more than 90 per cent variance of total CVD modeled by up to 180 d/o SHCs at 98 per cent global gridpoints. Globally, CVD modeled by GRACE loading can explain 72 per cent annual amplitude and 69 per cent variance of GPS observed height time-series, which is better than the GCMs results of 64 per cent for annual amplitude and 41 per cent for variance. Using a three cornered hat method, we also show that the noise level of monthly averaged CVD is about 3 mm for both GPS height time-series and GRACE loading result, while that of GCMs result is only 1.3 mm.

Description: We present a new method to derive 3-D surface deformation from an integration of interferometric synthetic aperture radar (InSAR) images and Global Navigation Satellite System (GNSS) observations based on Akaike's Bayesian Information Criterion (ABIC), considering relationship between deformations of neighbouring locations. This method avoids interpolated errors by excluding the interpolation of GNSS into the same spatial resolution as InSAR images and harnesses the data sets and the prior smooth constraints of surface deformation objectively and simultaneously by using ABIC, which were inherently unresolved in previous studies. In particular, we define surface roughness measuring smoothing degree to evaluate the performance of the prior constraints and deduce the formula of the covariance for the estimation errors to estimate the uncertainty of modelled solution. We validate this method using synthetic tests and the 2008 M w 7.9 Wenchuan earthquake. We find that the optimal weights associated with ABIC minimum are generally at trade-off locations that balance contributions from InSAR, GNSS data sets and the prior constraints. We use this method to evaluate the influence of the interpolated errors from the Ordinary Kriging algorithm on the derivation of surface deformation. Tests show that the interpolated errors may contribute to biasing very large weights imposed on Kriged GNSS data, suggesting that fixing the relative weights is required in this case. We also make a comparison with SISTEM method, indicating that our method allows obtaining better estimations even with sparse GNSS observations. In addition, this method can be generalized to provide a solution for situations where some types of data sets are lacking and can be exploited further to account for data sets such as the integration of displacements along radar lines and offsets along satellite tracks.

Description: While it has been known for some time that offsets in the time-series of Global Navigation Satellite System (GNSS) position estimates degrade station velocity determinations, the magnitude of the effect has not been clear. Using products of the International GNSS Service (IGS), we assess the impact empirically by injecting progressively larger numbers of artificial offsets and solving for a series of long-term secular GNSS frames. Our results show that the stability of the IGS global frame datum is fairly robust, with significant effects at the formal error level only for the R x (and Y-pole) and R z rotational orientations. On the other hand, station velocity estimates are more seriously affected, especially the vertical component. For the typical IGS station, the mean vertical rate uncertainty is already limited to 0.34 mm yr –1 for the current set of position discontinuities. If the number of breaks doubles, which might occur using newer detection schemes, then that uncertainty will worsen by ~40 per cent to 0.48 mm yr –1 . This error source is generally a more important component of realistic velocity uncertainties than any other, including accounting for temporal correlations in the GNSS data. The only way to improve future GNSS velocity estimates is to severely limit manmade displacements at the tracking stations.

Description: We have extended backwards from 2001 to 1979 the current release 05 (RL05) of the Gravity Recovery and Climate Experiment (GRACE) Atmospheric and Oceanic De-aliasing Level-1B (AOD1B) product and studied the impact of this and a previous release 04 (RL04) of the AOD1B product on precise orbits of five altimetry satellites (ERS-1, ERS-2, TOPEX/Poseidon, Envisat and Jason-1) for the time span 1991–2012, as compared to the case when no AOD1B product is used. We have found that using AOD1B RL05 product reduces root mean square (RMS) fits of satellite laser ranging (SLR) observations by about 1.0–6.4 per cent, 2-d arc overlaps in radial, cross-track and along-track directions by about 1.3–12.0, 0.3–10.0 and 2.0–10.0 per cent, respectively, for various satellites tested, as compared to the case without AOD1B product. Using AOD1B RL05 product instead of RL04 one reduces SLR RMS fits by 0.1–0.7 per cent, 2-d arc overlaps in radial, cross-track and along-track directions by 0.1–0.6, 0.1–1.3 and 0.2–1.2 per cent, respectively, for the satellite orbits tested. The multi-mission crossover analysis shows that the application of an AOD1B product reduces the scatter of radial errors by 0.4–2.8 per cent for the satellite missions studied. At the regions with the most pronounced changes the use of the AOD1B products improves the consistency between the sea level as measured by the TOPEX and ERS-2 missions and by the Jason-1 and Envisat missions by 5 to 10 per cent (globally by about 2 per cent). The results of our study show that extended AOD1B RL05 product performs better than the AOD1B RL04 and improves orbits of altimetry satellites and consistency of sea level products.

Description: In general, observations are normally considered to refer to an epoch in time, however, observations take time. During this time span temporal variations of the observable alias the measurement. Similar phenomenon can be defined in the space domain as well: data treated to refer to a geographical location often contains integrated information of the surroundings. In each case the appropriate signal content can partially be recovered by desmoothing the averaged data. The present study delivers the theoretical foundation of a desmoothing method, and suggests its use on different applications in geodesy. The theoretical formulation of the desmoothing has been derived for 1-D and 2-D signals, the latter is interpreted on a plain and also on a sphere. The presented case studies are less elaborated, but intended to demonstrate the need and usefulness of the desmoothing tool.

Description: Extracting geophysical signals from Global Positioning System (GPS) coordinate time-series is a well-established practice that has led to great insights into how the Earth deforms. Often small discontinuities are found in such time-series and are traceable to either broad-scale deformation (i.e. earthquakes) or discontinuities due to equipment changes and/or failures. Estimating these offsets accurately enables the identification of coseismic deformation estimates in the former case, and the removal of unwanted signals in the latter case which then allows tectonic rates to be estimated more accurately. We develop a method to estimate accurately discontinuities in time series of GPS positions at specified epochs, based on a so-called ‘offset series’. The offset series are obtained by varying the amount of GPS data before and after an event while estimating the offset. Two methods, a mean and a weighted mean method, are then investigated to produce the estimated discontinuity from the offset series. The mean method estimates coseismic offsets without making assumptions about geophysical processes that may be present in the data (i.e. tectonic rate, seasonal variations), whereas the weighted mean method includes estimating coseismic offsets with a model of these processes. We investigate which approach is the most appropriate given certain lengths of available data and noise within the time-series themselves. For the Sumatra–Andaman event, with 4.5 yr of pre-event data, we show that between 2 and 3 yr of post-event data are required to produce accurate offset estimates with the weighted mean method. With less data, the mean method should be used, but the uncertainties of the estimated discontinuity are larger.

Description: A new analytical method for the computation of a truncated series of solid spherical harmonic coefficients (HCs) from data on a spheroid (i.e. an oblate ellipsoid of revolution) is derived, using a transformation between surface and solid spherical HCs. A two-step procedure is derived to extend this transformation beyond degree and order (d/o) 520. The method is compared to the Hotine–Jekeli transformation in a numerical study based on the EGM2008 global gravity model. Both methods are shown to achieve submicrometre precision in terms of height anomalies for a model to d/o 2239. However, both methods result in spherical harmonic models that are different by up to 7.6 mm in height anomalies and 2.5 mGal in gravity disturbances due to the different coordinate system used. While the Hotine–Jekeli transformation requires the use of an ellipsoidal coordinate system, the new method uses only spherical polar coordinates. The Hotine–Jekeli transformation is numerically more efficient, but the new method can more easily be extended to cases where (a linear combination of) normal derivatives of the function under consideration are given on the surface of the spheroid. It therefore provides a solution to many types of ellipsoidal boundary-value problems in the spectral domain.

Description: Growth media have been developed to facilitate the enrichment and isolation of acidophilic and acid-tolerant sulfate-reducing bacteria (aSRB) from environmental and industrial samples, and to allow their cultivation in vitro . The main features of the ‘standard’ solid and liquid devised media are as follows: (i) use of glycerol rather than an aliphatic acid as electron donor; (ii) inclusion of stoichiometric concentrations of zinc ions to both buffer pH and to convert potentially harmful hydrogen sulphide produced by the aSRB to insoluble zinc sulphide; (iii) inclusion of Acidocella aromatica (an heterotrophic acidophile that does not metabolize glycerol or yeast extract) in the gel underlayer of double layered (overlay) solid media, to remove acetic acid produced by aSRB that incompletely oxidize glycerol and also aliphatic acids (mostly pyruvic) released by acid hydrolysis of the gelling agent used (agarose). Colonies of aSRB are readily distinguished from those of other anaerobes due to their deposition and accumulation of metal sulphide precipitates. Data presented illustrate the effectiveness of the overlay solid media described for isolating aSRB from acidic anaerobic sediments and low pH sulfidogenic bioreactors.

Description: We explore Earth's elastic deformation response to ocean tidal loading (OTL) using kinematic Global Positioning System (GPS) observations and forward-modelled predictions across South America. Harmonic coefficients are extracted from up to 14 yr of GPS-inferred receiver locations, which we estimate at 5 min intervals using precise point positioning. We compare the observed OTL-induced surface displacements against predictions derived from spherically symmetric, non-rotating, elastic and isotropic (SNREI) Earth models. We also compare sets of modelled predictions directly for various ocean-tide and Earth-model combinations. The vector differences between predicted displacements computed using separate ocean-tide models reveal uniform-displacement components common to all stations in the South America network. Removal of the network-mean OTL-induced displacements from each site substantially reduces the vector differences between observed and predicted displacements. We focus on the dominant astronomical tidal harmonics from three distinct frequency bands: semidiurnal (M 2 ), diurnal (O 1 ) and fortnightly (M f ). In each band, the observed OTL-induced surface displacements strongly resemble the modelled displacement-response patterns, and the residuals agree to about 0.3 mm or better. Even with the submillimetre correspondence between observations and predictions, we detect regional-scale spatial coherency in the final set of residuals, most notably for the M 2 harmonic. The spatial coherency appears relatively insensitive to the specific choice of ocean-tide or SNREI-Earth model. Varying the load model or 1-D elastic structure yields predicted OTL-induced displacement differences of order 0.1 mm or less for the network. Furthermore, estimates of the observational uncertainty place the noise level below the magnitude of the residual displacements for most stations, supporting our interpretation that random errors cannot account for the entire misfit. Therefore, the spatially coherent residuals may reveal deficiencies in the a priori SNREI Earth models. In particular, the residuals may indicate sensitivity to regional deviations from standard globally averaged Earth structure due to the presence of the South American craton.

Description: A common dye of prussian blue (PB) as an indicator was used to develop a colorimetric method for detecting the efficacy of the antibiotics in vitro. Considering the electronic production capacity of microbial respiration, ferricyanide was employed in transferring electrons from target microorganism of Escherichia coli ( E. coli ) to produce ferrocyanide. Subsequently, ferrocyanide reacted with ferric ions to form PB. In view of relationship between the PB yield and the bacterial activity, the efficacy of the antibiotics on E. coli was directly detected at 700 nm of PB absorption. When the 5% activity of antibiotics on 20 isolates of E. coli was quantified as 5% efficacy, the applied concentrations of eight antibiotics, such as cefepime, ceftriaxone sodium, cefoperazone sodium, piperacillin sodium, amoxicillin, gentamicin, amikacin and levofloxacin were 2, 2, 4, 4, 10, 4, 8 and 8 μg mL –1 , respectively. To compare with minimum inhibitory concentration results obtained by Clinical and Laboratory Standards Institute broth macrodilution method, the results of PB methods showed good agreements except with gentamicin. Paired t- test result ( P ) also showed that difference between two methods was statistically significant ( P = 0.006).

Description: Planar faults are widely adopted during inversions to determine slip distributions and fault geometries using geodetic observations; however, little research has been conducted with respect to curved faults. We attribute this to the lack of an appropriate parameterized modelling method. In this paper, we present a curved-fault modelling method (CFMM) that describes a curved fault according to specific parameters, and we also develop a corresponding hybrid iterative inversion algorithm (HIIA) to perform inversions for parametric curved-fault geometries and slips. The results of the strike-component and dip-component synthetic tests show that a complex S-shaped fault surface and a circular slip distribution are successfully recovered, indicating the strong performance of the CFMM and HIIA methods. In addition, we describe and verify a scenario for determining the number of necessary geometrical parameters for the HIIA and examine the case study of the Wenchuan earthquake, which occurred on a complex listric fault surface. During the iteration process of the HIIA, both the fault geometry and slip distribution of the Beichuan and Pengguan faults converge to optimal values, indicating a Beichuan fault (BCF) model with a continuous listric shape and gradual steepening from the southwest to the northeast, which is highly consistent with geological survey results. Both the synthetic and real-world case studies show that the HIIA and the CMFF are superior to the conventional fault modelling method based on rectangular planes and that these models have the potential for use in more integrated research involving inversion studies, such as joint slip/curved-fault-geometry inversions that take into account data resolving power.

Description: Metarhizium acridum is an entomopathogenic fungus commonly used as a bioinsecticide. The conidium is the fungal stage normally employed as field inoculum in biological control programs and must survive under field conditions such as high ultraviolet-B (UV-B) exposure. Light, which is an important stimulus for many fungi, has been shown to induce the production of M. robertsii conidia with increased stress tolerance. Here we show that a two-hour exposure to white or blue/UV-A light of fast-growing mycelium induces tolerance to subsequent UV-B irradiation. Red light, however, does not have the same effect. In addition, we established that this induction can take place with as little as 1 min of white-light exposure. This brief illumination scheme could be relevant in future studies of M. acridum photobiology and for the production of UV-B resistant mycelium used in mycelium-based formulations for biological control.

Description: Photorhabdus (Enterobacteriaceae) bacteria are pathogenic to insects and mutualistic with entomopathogenic Heterorhabditis nematodes . Photorhabdus luminescens subsp. akhurstii LN2, associated with Heterorhabditis indica LN2, shows nematicidal activity against H. bacteriophora H06 infective juveniles (IJs). In the present study, an rpoS mutant of P. luminescens LN2 was generated through allelic exchange to examine the effects of rpoS deletion on the nematicidal activity and nematode development. The results showed that P. luminescens LN2 required rpoS for nematicidal activity against H06 nematodes, normal IJ recovery and development of H. indica LN2, however, not for the bacterial colonization in LN2 and H06 IJs. This provides cues for further understanding the role of rpoS in the mutualistic association between entomopathogenic nematodes and their symbionts.

Description: Catechol 2, 3-dioxygenase (C23O) is the key enzyme for aerobic aromatic degradation. Based on clone libraries and quantitative real-time polymerase chain reaction, we characterized diversity and distribution patterns of C23O genes in surface sediments of the Bohai Sea. The results showed that sediments of the Bohai Sea were dominated by genes related to C23O subfamily I.2.A. The samples from wastewater discharge area (DG) and aquaculture farm (KL) showed distinct composition of C23O genes when compared to the samples from Bohai Bay (BH), and total organic carbon was a crucial determinant accounted for the composition variation. C6BH12-38 and C2BH2-35 displayed the highest gene copies and highest ratios to the 16S rRNA genes in KL, and they might prefer biologically labile aromatic hydrocarbons via aquaculture inputs. Meanwhile, C7BH3-48 showed the highest gene copies and highest ratios to the 16S rRNA genes in DG, and this could be selective effect of organic loadings from wastewater discharge. An evident increase in C6BH12-38 and C7BH3-48 gene copies and reduction in diversity of C23O genes in DG and KL indicated composition perturbations of C23O genes and potential loss in functional redundancy. We suggest that ecological habitat and trophic specificity could shape the distribution of C23O genes in the Bohai Sea sediments.

Description: LysR-type transcriptional regulators (LTTRs) regulate various cellular processes in bacteria. pnpR is an LTTR-encoding gene involved in the regulation of hydroquinone (HQ) degradation, and its effects on the cellular processes of Pseudomonas putida DLL-E4 were investigated at the physiological, biochemical and molecular levels. Reverse transcription polymerase chain reaction revealed that pnpR positively regulated its own expression and that of the pnpC1C2DECX1X2 operon; additionally, pnpR partially regulated the expression of pnpA when P. putida was grown on para -nitrophenol (PNP) or HQ. Strains DLL-E4 and DLL- pnpR exhibited similar cellular morphologies and growth rates. Transcriptome analysis revealed that pnpR regulated the expression of genes in addition to those involved in PNP degradation. A total of 20 genes were upregulated and 19 genes were downregulated by at least 2-fold in strain DLL- pnpR relative to strain DLL-E4. Bioinformatic analysis revealed putative PnpR-binding sites located in the upstream regions of genes involved in PNP degradation, carbon catabolite repression and other cellular processes. The utilization of L-aspartic acid, L-histidine, L-pyroglutamic acid, L-serine, -aminobutyric acid, D,L-lactic acid, D-saccharic acid, succinic acid and L-alaninamide was increased at least 1.3-fold in strain DLL- pnpR as shown by BIOLOG assays, indicating that pnpR plays a potential negative regulation role in the utilization of carbon sources.

Description: Globally gridded estimates of monthly-mean anomalies of terrestrial water storage (TWS) are estimated from the most recent GRACE release 05a of GFZ Potsdam in order to provide non-geodetic users a convenient access to state-of-the-art GRACE monitoring data. We use an ensemble of five global land model simulations with different physics and different atmospheric forcing to obtain reliable gridded scaling factors required to correct for spatial leakage introduced during data processing. To allow for the application of this data-set for large-scale monitoring tasks, model validation efforts, and subsequently also data assimilation experiments, globally gridded estimates of TWS uncertainties that include (i) measurement, (ii) leakage and (iii) re-scaling errors are provided as well. The results are generally consistent with the gridded data provided by Tellus, but deviate in some basins which are largely affected by the uncertainties of the model information required for re-scaling, where the approach based on the median of a small ensemble of global land models introduced in this paper leads to more robust results.

Description: Landfills are significant global sources of atmospheric methane, but little is known about the ecology and community structure of methanogens in these sites. Here, we investigated the methanogen community based on methyl coenzyme M reductase A gene amplicons in the vertical profiles of three different sites at a municipal landfill complex in China. Links between methanogen communities and refuse properties were explored using multivariate analysis. Clone library results showed that most clones (92%) were related to the hydrogenotrophic methanogens, Methanomicrobiales. Almost all of the Methanomicrobiales clones retrieved in this study are members of the genus Methanoculleus . Eight clones were affiliated with the genus Methanofollis . The remaining clones were clustered within the genus Methanosarcina . Terminal restriction fragment length polymorphism profiles showed that the landfill was predominated by 22 taxa, making up 69%–96% of the community. Of these, a single taxon comprised 36%–65% of the communities across all sites and depths. Principal components analysis separated the methanogen community into three groups, irrespective of site or depth. Redundancy analysis suggested that total phosphorus and pH play roles in structuring methanogen communities in landfills.

Description: Regional recovery of the disturbing gravitational potential in the area of Central Europe from satellite gravitational gradients data is discussed in this contribution. The disturbing gravitational potential is obtained by inverting surface integral formulas which transform the disturbing gravitational potential onto disturbing gravitational gradients in the spherical local north-oriented frame. Two numerical approaches that solve the inverse problem are considered. In the first approach, the integral formulas are rigorously decomposed into two parts, that is, the effects of the gradient data within near and distant zones. While the effect of the near zone data is sought as an inverse problem, the effect of the distant zone data is synthesized from the global gravitational model GGM05S using spectral weights given by truncation error coefficients up to the degree 150. In the second approach, a reference gravitational field up to the degree 180 is applied to reduce and smooth measured gravitational gradients. In both cases we recovered the disturbing gravitational potential from each of the four well-measured gravitational gradients of the GOCE satellite separately as well as from their combination. Obtained results are compared with the EGM2008, DIR-r2, TIM-r2 and SPW-r2 global gravitational models. The best fit was achieved for EGM2008 and the second approach combining all four well-measured gravitational gradients with rms of 1.231 m 2  s –2 .

Description: Here we report a newly identified ‘Chalky back’ phenomenon in banana prawns ( Fenneropenaeus merguiensis ) farmed in North Queensland, Australia. This was characterized by localized white discoloured segmentation of the cervical groove, moreover, after cooking the prawns exploded, making them unfit for commercial sale. Histological examination revealed breakdown of gut and abdominal muscle tissue in some moribund specimens. We selectively isolated Vibrio spp., which are known prawn pathogens, from healthy and Chalky back specimens. Isolated bacteria were identified, typed and tested for the presence of eight virulence genes (VGs), biofilm formation, adherence and cytotoxicity to fish cells. In all, 32 isolates were recovered and identified as Vibrio harveyi , V. owensii , V. sinaloensis -like, V. campbellii , V. shilonii , Vibrio sp. and Photobacterium damselae using 16S rRNA gene sequencing. All V. harveyi carried VGs coding for haemolysin, tox R and flagella; formed biofilm; and adhered to both cell lines. This was similar to the V. sinaloensis -like strains that were only isolated from Chalky back specimens. Our data suggest that Vibrio spp. may play a role in the pathogenesis of Chalky back. This study is the first report of Chalky back phenomenon in farmed banana prawns that needs to be closely monitored by the industry.

Description: We present efficient Fourier-domain algorithms for modelling gravity effects due to topographic masses. The well-known Parker's formula originally based on the standard fast Fourier transform (FFT) algorithm is modified by applying the Gauss–FFT method instead. Numerical precision of the forward and inverse Fourier transforms embedded in Parker's formula and its extended forms are significantly improved by the Gauss–FFT method. The topographic model is composed of two major aspects, the geometry and the density. Versatile geometric representations, including the mass line model, the mass prism model, the polyhedron model and smoother topographic models interpolated from discrete data sets using high-order splines or pre-defined by analytical functions, in combination with density distributions that vary both laterally and vertically in rather arbitrary ways following exponential or general polynomial functions, now can be treated in a consistent framework by applying the Gauss–FFT method. The method presented has been numerically checked by space-domain analytical and hybrid analytical/numerical solutions already established in the literature. Synthetic and real model tests show that both the Gauss–FFT method and the standard FFT method run much faster than space-domain solutions, with the Gauss–FFT method being superior in numerical accuracy. When truncation errors are negligible, the Gauss–FFT method can provide forward results almost identical to space-domain analytical or semi-numerical solutions in much less time.

Description: The gravity gradient tensor (GGT) has been increasingly used in practical applications, but the advantages and the disadvantages of the analysis of GGT components versus the analysis of the vertical component of the gravity field are still debated. We analyse the performance of joint inversion of GGT components versus separate inversion of the gravity field alone, or of one tensor component. We perform our analysis by inspection of the Picard Plot, a Singular Value Decomposition tool, and analyse both synthetic data and gradiometer measurements carried out at the Vredefort structure, South Africa. We show that the main factors controlling the reliability of the inversion are algebraic ambiguity (the difference between the number of unknowns and the number of available data points) and signal-to-noise ratio. Provided that algebraic ambiguity is kept low and the noise level is small enough so that a sufficient number of SVD components can be included in the regularized solution, we find that: (i) the choice of tensor components involved in the inversion is not crucial to the overall reliability of the reconstructions; (ii) GGT inversion can yield the same resolution as inversion with a denser distribution of gravity data points, but with the advantage of using fewer measurement stations.

Description: This paper compares GRACE (Gravity Recovery and Climate Experiment) and ICESat (Ice, Cloud and land Elevation Satellite) observations to confirm whether the observed gravity increase in the Tibetan Plateau (TP) was primarily caused by lake storage gain, and comprehensively analyses the changing pattern of lake level over 2003–2009. An improved automated method was used to obtain lake-level changes and the underestimation of lake water storage was considered due to lake area expansion and lake density. The result demonstrates that GRACE recorded a mass gain (16.43 ± 1.65/11.79 ± 1.25 gt a –1 ) in the total/inner TP, of which lake storage increase accounts for (8.78 ± 0.75/7.53 ± 0.56 gt a –1 ) based on ICESat. The northwestern residual may be stored in new lakes and soil moisture as a result of net precipitation gain. According to the character of the lake-level changes, we divide the TP into four subregions. Generally, the changing pattern of lake level concurs with the distribution of precipitation, which is increasing in the inner TP and decreasing in the upstream area of the Indus and Brahmaputra Rivers. An excess of rainfall in the northeastern TP in the summer of 2005 and 2009 caused a simultaneous large increase in water level in many lakes. The correlation of lake changes with precipitation demonstrates that precipitation rather than glacial melt is the main cause of lake-level change in most places. Nonetheless, the meltwater is a considerable supplement for lakes near glaciers such as Selin Co and Nam Co, which partly explains why GRACE indicates a much weaker signal in this region.

Description: Traditional processing of Global Navigation Satellite System (GNSS) data using dedicated scientific software has provided the highest levels of positional accuracy, and has been used extensively in geophysical deformation studies. To achieve these accuracies a significant level of understanding and training is required, limiting their availability to the general scientific community. Various online GNSS processing services, now freely available, address some of these difficulties and allow users to easily process their own GNSS data and potentially obtain high quality results. Previous research into these services has focused on Continually Operating Reference Station (CORS) GNSS data. Less research exists on the results achievable with these services using large campaign GNSS data sets, which are inherently noisier than CORS data. Even less research exists on the quality of velocity fields derived from campaign GNSS data processed through online precise point positioning services. Particularly, whether they are suitable for geodynamic and deformation studies where precise and reliable velocities are needed. In this research, we process a very large campaign GPS data set (spanning 10 yr) with the online Jet Propulsion Laboratory Automated Precise Positioning Service. This data set is taken from a GNSS network specifically designed and surveyed to measure deformation through the central North Island of New Zealand. This includes regional CORS stations. We then use these coordinates to derive a horizontal and vertical velocity field. This is the first time that a large campaign GPS data set has been processed solely using an online service and the solutions used to determine a horizontal and vertical velocity field. We compared this velocity field to that of another well utilized GNSS scientific software package. The results show a good agreement between the CORS positions and campaign station velocities obtained from the two approaches. We discuss the implications of these results for how future GNSS campaign field surveys might be conducted and how their data might be processed.

Description: We have developed a method to estimate the geometry, location and densities of anomalies coming from 2-D gravity data based on compact gravity inversion technique. Compact gravity inversion is simple, fast and user friendly but severely depends on the number of model parameters, that is, by increasing the model parameters, the anomalies tend to concentrate near the surface. To overcome this ambiguity new weighting functions based on density contrast, depth, and compactness models have been introduced. Variable compactness factors have been defined here to get either a sharp or a smooth model based on the depth of the source or existence of prior information. Depth weighting derived from one station of gravity data whereas the effect of gravity data is 2-D and 3-D. To compensate this limitation an innovating weighting function namely kernel function has been introduced which multiplies with weight and compactness matrixes to yield a general model weighting function. The method is tested using three different sets of synthetic examples: a body at various depths (20, 40, 80 and 140 m), two bodies at the same depth but various distances to estimate lateral resolution and three bodies with negative and positive density contrast in different depths. The method is also applied to three real gravity data of Woodlawn massive sulphide body, sulphides mineralization of British Colombia and iron ore body of Missouri. The method produces solutions consistent with the known geologic attributes of the gravity sources, illustrating its potential practicality.

Description: This paper resurrects a version of Poisson's Partial Differential Equation (PDE) associated with the gravitational field at the Earth's surface and illustrates how the PDE possesses a capability to extract the mass density of Earth's topography from land-based gravity data. Herein, first we propound a theorem which mathematically introduces this version of Poisson's PDE adapted for the Earth's surface and then we use this PDE to develop a method of approximating the terrain mass density. Also, we carry out a real case study showing how the proposed approach is able to be applied to a set of land-based gravity data. In the case study, the method is summarized by an algorithm and applied to a set of gravity stations located along a part of the north coast of the Persian Gulf in the south of Iran. The results were numerically validated via rock-samplings as well as a geological map. Also, the method was compared with two conventional methods of mass density reduction. The numerical experiments indicate that the Poisson PDE at the Earth's surface has the capability to extract the mass density from land-based gravity data and is able to provide an alternative and somewhat more precise method of estimating the terrain mass density.

Description: Plant-growth-promoting bacteria belonging to Azospirillum and Pseudomonas genera are major inhabitants of the rhizosphere. Both are increasingly commercialized as crops inoculants. Interspecific interaction in the rhizosphere is critical for inoculants aptness. The objective of this work was to evaluate Azospirillum and Pseudomonas interaction in mixed biofilms by co-cultivation of the model strains Azospirillum brasilense Sp245 and Pseudomonas protegens CHA0. The results revealed enhanced growth of both strains when co-cultured in static conditions. Moreover, Sp245 biofilm formed in plastic surfaces was increased 2-fold in the presence of CHA0. Confocal microscopy revealed highly structured mixed biofilms showing Sp245 mainly on the bottom and CHA0 towards the biofilm surface. In addition, A. brasilense biofilm was thicker and denser when co-cultured with P. protegens. In a colony–colony interaction assay, Sp245 changed nearby CHA0 producing small colony phenotype, which accounts for a diffusible metabolite mediator; though CHA0 spent medium did not affect Sp245 colony phenotype. Altogether, these results point to a cooperative interaction between A. brasilense Sp245 and P. protegens CHA0 in which both strains increase their static growth and produce structured mixed biofilms with a strain-specific distribution.

Description: Knowledge about the factors shaping the rumen microbiota in wild animals is limited. Therefore, the aim of this study was to compare the microbiota from the three cervid species moose ( Alces alces , n = 5), red deer ( Cervus elaphus , n = 4) and roe deer ( Capreolus capreolus , n = 12), sharing the same habitat. Using deep 16S rRNA gene sequencing, we found that the largest species moose had the highest number of unique operational taxonomic units. Furthermore, red deer and moose shared more of the microbiota, compared with the smallest species, roe deer, with Firmicutes and Euryarchaeota being significantly overrepresented for the shared microbiota. These differences could not be explained by diet or range. The animals largely shared the same range, and there are no systematic differences in diet. We therefore believe rumen physiology can be one of the main contributing factors to the observed distribution of the rumen microbiota in cervid species.

Description: Mesorhizobium loti MAFF303099 has a functional Type III secretion system (T3SS) that is involved in the determination of competitiveness for legume nodulation. Here we demonstrate that the transcriptional factor TtsI, which positively regulates T3SS genes expression, is involved in a negative regulation of M. loti swimming motility in soft-agar. Conditions that induce T3SS expression affect flagella production. The same conditions also affect promoter activity of M. loti visN gene, a homolog to the positive regulator of flagellar genes that has been described in other rhizobia. Defects in T3SS complex assembly at membranes limited the negative regulation of motility by the expression of TtsI.

Description: The Moho surface can be determined according to isostatic theories, and among them, the recent Vening Meinesz-Moritz (VMM) theory of isostasy has been successfully applied for this purpose. In this paper, this method is studied from a theoretical prospective and its connection to the Airy–Heiskanen (AH) and Vening Meinesz original theories are presented. Jeffrey's inverse solution to isostasy is developed according to the recent developments of the VMM method and both are compared in similar situations. It is shown that they are generalizations of the AH model in a global and continuous domain. In the VMM spherical harmonic solution for Moho depth, the mean Moho depth contributes only to the zero-degree term of the series, while in Jeffrey's solution it contributes to all frequencies. In addition, the VMM spherical harmonic series is improved further so that the mean Moho can contribute to all frequencies of the solution. This modification makes the VMM global solution superior to the Jeffrey one, but in a global scale, the difference between both solutions is less than 3 km. Both solutions are asymptotically convergent and we present two methods to obtain smooth solutions for Moho from them.

Description: Frequent burning is commonly undertaken to maintain diversity in temperate grasslands of southern Australia. How burning affects below-ground fungal community diversity remains unknown. We show, using a fungal rDNA metabarcoding approach (Illumina MiSeq), that the fungal community composition was influenced by fire regime (frequency) but not time-since-fire. Fungal community composition was resilient to direct fire effects, most likely because grassland fires transfer little heat to the soil. Differences in the fungal community composition due to fire regime was likely due to associated changes that occur in vegetation with recurrent fire, via the break up of obligate symbiotic relationships. However, fire history only partially explains the observed dissimilarity in composition among the soil samples, suggesting a distinctiveness in composition in each grassland site. The importance of considering changes in soil microbe communities when managing vegetation with fire is highlighted.

Description: Soil is thought to be important both as a source and a sink of carbonyl sulfide (COS) in the troposphere, but the mechanism affecting COS uptake, especially for fungi, remains uncertain. Fungal isolates that were collected randomly from forest soil showed COS-degrading ability at high frequencies: 38 out of 43 isolates grown on potato dextrose agar showed degradation of 30 ppmv COS within 24 h. Of these isolates, eight degraded 30 ppmv of COS to below the detection limit within 2 h. These isolates also showed an ability to degrade COS included in ambient air (around 500 pptv) and highly concentrated (12 500 ppmv) level, even though the latter is higher than the lethal level for mammals. COS-degrading activity was estimated by using ergosterol as a biomass index for fungi. Trichoderma sp. THIF08 had the highest COS-degrading activity of all the isolates. Interestingly, Umbelopsis/Mortierella spp. THIF09 and THIF13 were unable to degrade 30 ppmv COS within 24 h, and actually emitted COS during the cultivation in ambient air. These results indicate a fungal contribution to the flux of COS between the terrestrial and atmospheric environments.

Description: A composite transposon is a mobile genetic element consisting of two insertion sequences (ISs) flanking a segment of cargo DNA often containing antibiotic resistance (AR) genes. Composite transposons can move as a discreet unit. There have been recently several reports on a novel mechanism of movement of an IS 26 -based composite transposon through the formation of a translocatable unit (TU), carrying the internal DNA segment of a composite transposon and one copy of a flanking IS. In this study, we determined the presence of composite transposons and TUs in human oral metagenomic DNA using PCR primers from common IS elements. Analysis of resulting amplicons showed four different IS 1216 composite transposons and one IS 257 composite transposon in our metagenomic sample. As our PCR strategy would also detect TUs, PCR was carried out to detect circular TUs predicted to originate from these composite transposons. We confirmed the presence of two novel TUs, one containing an experimentally proven antiseptic resistance gene and another containing a putative universal stress response protein (UspA) encoding gene. This is the first report of a PCR strategy to amplify the DNA segment on composite transposons and TUs in metagenomic DNA. This can be used to identify AR genes associated with a variety of mobile genetic elements from metagenomes.

Description: Carbonyl sulfide (COS) is an atmospheric trace gas and one of the sources of stratospheric aerosol contributing to climate change. Although one of the major sinks of COS is soil, the distribution of COS degradation ability among bacteria remains unclear. Seventeen out of 20 named bacteria belonging to Actinomycetales had COS degradation activity at mole fractions of 30 parts per million by volume (ppmv) COS. Dietzia maris NBRC 15801 T and Mycobacterium sp. THI405 had the activity comparable to a chemolithoautotroph Thiobacillus thioparus THI115 that degrade COS by COS hydrolase for energy production. Among 12 bacteria manifesting rapid degradation at 30 ppmv COS, D. maris NBRC 15801 T and Streptomyces ambofaciens NBRC 12836 T degraded ambient COS (~500 parts per trillion by volume). Geodermatophilus obscurus NBRC 13315 T and Amycolatopsis orientalis NBRC 12806 T increased COS concentrations. Moreover, six of eight COS-degrading bacteria isolated from soils had partial nucleotide sequences similar to that of the gene encoding clade D of β-class carbonic anhydrase, which included COS hydrolase. These results indicate the potential importance of Actinomycetes in the role of soils as sinks of atmospheric COS.

Description: In a pioneering study, Wahr & Bergen developed the widely adopted, pseudo-normal mode framework for predicting the impact of anelastic effects on the Earth's body tides. Lau et al. have recently derived an extended normal mode treatment of the problem (as well as a minor variant of the theory known as the direct solution method) that makes full use of theoretical developments in free oscillation seismology spanning the last quarter century and that avoids a series of assumptions and approximations adopted in the traditional theory for predicting anelastic effects. There are two noteworthy differences between these two theories: (1) the traditional theory only considers perturbations to the eigenmodes of an elastic Earth, whereas the new theory augments this set of modes to include the relaxation modes that arise in anelastic behaviour; and (2) the traditional theory approximates the complex perturbation to the tidal Love number as a scaled version of the complex perturbation to the elastic moduli, whereas the new theory computes the full complex perturbation to each eigenmode. In this study, we highlight the above differences using a series of synthetic calculations, and demonstrate that the traditional theory can introduce significant error in predictions of the complex perturbation to the Love numbers due to anelasticity and the related predictions of tidal lag angles. For the simplified Earth models we adopt, the computed lag angles differ by ~20 per cent. The assumptions in the traditional theory have important implications for previous studies that use model predictions to correct observables for body tide signals or that analyse observations of body tide deformation to infer mantle anelastic structure. Finally, we also highlight the fundamental difference between apparent attenuation (i.e. attenuation inferred from observations or predicted using the above theories) and intrinsic attenuation (i.e. the material property investigated through experiments), where both are often expressed in terms of lag angles or Q –1 . In particular, we demonstrate the potentially significant (factor of two or more) bias introduced in estimates of Q –1 and its frequency dependence in studies that have treated Q –1 determined from tidal phase lags or measured experimentally as being equal. The observed or theoretically predicted lag angle (or apparent Q –1 ) differs from the intrinsic, material property due to inertia, self-gravity and effects associated with the energy budget. By accounting for these differences we derive, for a special case, an expression that accurately maps apparent attenuation predicted using the extended normal mode formalism of Lau et al. into intrinsic attenuation. The theory allows for more generalized mappings which may be used to robustly connect observations and predictions of tidal lag angles to results from laboratory experiments of mantle materials.

Description: We consider a new approach to both the forward and inverse problems in post-seismic deformation. We present a method for forward modelling post-seismic deformation in a self-gravitating, heterogeneous and compressible earth with a variety of linear and nonlinear rheologies. We further demonstrate how the adjoint method can be applied to the inverse problem both to invert for rheological structure and to calculate the sensitivity of a given surface measurement to changes in rheology or time-dependence of the source. Both the forward and inverse aspects are illustrated with several numerical examples implemented in a spherically symmetric earth model.