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: The 16S rRNA gene (16S rDNA) codes for RNA that plays a fundamental role during translation in the ribosome and is used extensively as a marker gene to establish relationships among bacteria. However, the complementary non-coding 16S rDNA (nc16S rDNA) has been ignored. An idea emerged in the course of analyzing bacterial 16S rDNA sequences in search for nucleotide composition and substitution patterns: Does the nc16S rDNA code? If so, what does it code for? More importantly: Does 16S rDNA evolution reflect its own evolution or the evolution of its counterpart nc16S rDNA? The objective of this minireview is to discuss these thoughts. nc strands often encode small RNAs (sRNAs), ancient components of gene regulation. nc16S rDNA sequences from different bacterial groups were used to search for possible matches in the Bacterial Small Regulatory RNA Database. Intriguingly, the sequence of one published sRNA obtained from Legionella pneumophila (GenBank: AE017354.1) showed high non-random similarity with nc16S rDNA corresponding in part to the V5 region especially from Legionella and relatives. While the target(s) of this sRNA is unclear at the moment, its mere existence might open up a new chapter in the use of the 16S rDNA to study relationships among bacteria.

Description: Thermotolerance of the fungus Fomes sp. EUM1 was evaluated in solid state fermentation (SSF). This thermotolerant strain improved both hyphal invasiveness (38%) and length (17%) in adverse thermal conditions exceeding 30°C and to a maximum of 40°C. In contrast, hyphal branching decreased by 46% at 45°C. The production of cellulases over corn stover increased 1.6-fold in 30°C culture conditions, xylanases increased 2.8-fold at 40°C, while laccase production improved 2.7-fold at 35°C. Maximum production of lignocellulolytic enzymes was obtained at elevated temperatures in shorter fermentation times (8–6 days), although the proteases appeared as a thermal stress response associated with a drop in lignocellulolytic activities. Novel and multiple isoenzymes of xylanase (four bands) and cellulase (six bands) were secreted in the range of 20–150 kDa during growth in adverse temperature conditions. However, only a single laccase isoenzyme (46 kDa) was detected. This is the first report describing the advantages of a thermotolerant white-rot fungus in SSF. These results have important implications for large-scale SSF, where effects of metabolic heat are detrimental to growth and enzyme production, which are severely affected by the formation of high temperature gradients.

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: The rhizobacterium Serratia plymuthica 4Rx13 emits the novel and unique volatile sodorifen (C 16 H 26 ), which has a polymethylated bicyclic structure. Transcriptome analysis revealed that gene SOD_c20750 (annotated as terpene cyclase) is involved in the biosynthesis of sodorifen. Here we show that this gene is located in a small cluster of four genes ( SOD_c20750 – SOD_c20780 ), and the analysis of the knockout mutants demonstrated that SOD_c20760 (annotated as methyltransferase) and SOD_c20780 (annotated as isopentenyl pyrophosphate (IPP) isomerase) are needed for the biosynthesis of sodorifen, while a sodorifen-negative phenotype was not achieved with the SOD_c20770 (annotated as deoxy-xylulose-5-phosphate (DOXP) synthase) mutant. Altogether, the function of this new gene cluster was assigned to the biosynthesis of this structurally unusual volatile compound sodorifen.

Description: Among staphylococci Staphylococcus saprophyticus is the only species that is typically uropathogenic and an important cause of urinary tract infections in young women. The amino acid D-serine occurs in relatively high concentrations in human urine and has a bacteriostatic or toxic effect on many bacteria. In uropathogenic Escherichia coli and S. saprophyticus , the amino acid regulates the expression of virulence factors and can be used as a nutrient. The ability of uropathogens to respond to or to metabolize D-serine has been suggested as a factor that enables colonization of the urinary tract. Until now nothing is known about D-serine transport in S.   saprophyticus . We generated mutants of putative transporter genes in S.   saprophyticus 7108 that show homology to the D-serine transporter cyc A of E. coli and tested them in a D-serine depletion assay to analyze the D-serine uptake rate of the cells. The mutant of SPP1070 showed a strong decrease in D-serine uptake. Therefore, SSP1070 was identified as a major D-serine transporter in S. saprophyticus 7108 and was named D-serine transporter A (DstA). D-serine caused a prolonged lag phase of S. saprophyticus in a chemically defined medium. This negative effect was dependent on the presence of DstA.

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: Genes encoding fluoride transporters have been identified in bacterial and archaeal species. The genome sequence of the cariogenic  Streptococcus mutans  bacteria suggests the presence of a putative fluoride transporter, which is referred to as a chloride channel permease. Two homologues of this gene (GenBank locus tags SMU_1290c and SMU_1289c) reside in tandem in the genome of  S. mutans . The aim of this study was to determine whether the chloride channel permeases contribute to fluoride resistance. We constructed SMU_1290c- and SMU_1289c-knockout  S. mutans  UA159 strains. We also constructed a double-knockout strain lacking both genes. SMU_1290c or SMU_1289c was transformed into a fluoride transporter- disrupted  Escherichia coli strain. All bacterial strains were cultured under appropriate conditions with or without sodium fluoride, and fluoride resistance was evaluated. All three gene-knockout  S. mutans  strains showed lower resistance to sodium fluoride than did the wild-type strain. No significant changes in resistance to other sodium halides were recognized between the wild-type and double-knockout strains. Both SMU_1290c and SMU_1289c transformation rescued fluoride transporter-disrupted  E. coli  cell from fluoride toxicity. We conclude that the chloride channel permeases contribute to fluoride resistance in  S. mutans .

Description: Clostridium thermocellum and Thermoanaerobacterium saccharolyticum are bacteria under investigation for production of biofuels from plant biomass. Thermoanaerobacterium saccharolyticum has been engineered to produce ethanol at high yield (〉90% of theoretical) and titer (〉70 g/l). Efforts to engineer C. thermocellum have not, to date, been as successful, and efforts are underway to transfer the ethanol production pathway from T. saccharolyticum to C. thermocellum . One potential challenge in transferring metabolic pathways is the possibility of incompatible levels of nicotinamide cofactors. These cofactors (NAD + , NADH, NADP + and NADPH) and their oxidation state are important in the context of microbial redox metabolism. In this study we directly measured the concentrations and reduced oxidized ratios of these cofactors in a number of strains of C. thermocellum and T. saccharolyticum by using acid/base extraction and enzymatic assays. We found that cofactor ratios are maintained in a fairly narrow range, regardless of the metabolic network modifications considered. We have found that the ratios are similar in both organisms, which is a relevant observation in the context of transferring the T. saccharolyticum ethanol production pathway to C. thermocellum .

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: Measurements of ground deformation can be used to identify and interpret geophysical processes occurring at volcanoes. Most studies rely on a single geodetic technique, or fit a geophysical model to the results of multiple geodetic techniques. Here we present a methodology that combines GPS, Total Station measurements and InSAR into a single reference frame to produce an integrated 3-D geodetic velocity surface without any prior geophysical assumptions. The methodology consists of five steps: design of the network, acquisition and processing of the data, spatial integration of the measurements, time series computation and finally the integration of spatial and temporal measurements. The most significant improvements of this method are (1) the reduction of the required field time, (2) the unambiguous detection of outliers, (3) an increased measurement accuracy and (4) the construction of a 3-D geodetic velocity field. We apply this methodology to ongoing motion on Arenal's western flank. Integration of multiple measurement techniques at Arenal volcano revealed a deformation field that is more complex than that described by individual geodetic techniques, yet remains consistent with previous studies. This approach can be applied to volcano monitoring worldwide and has the potential to be extended to incorporate other geodetic techniques and to study transient deformation.

Description: In autumn 2012, the new release 05 (RL05) of monthly geopotencial spherical harmonics Stokes coefficients (SC) from Gravity Recovery and Climate Experiment (GRACE) mission was published. This release reduces the noise in high degree and order SC, but they still need to be filtered. One of the most common filtering processing is the combination of decorrelation and Gaussian filters. Both of them are parameters dependent and must be tuned by the users. Previous studies have analyzed the parameters choice for the RL05 GRACE data for oceanic applications, and for RL04 data for global application. This study updates the latter for RL05 data extending the statistics analysis. The choice of the parameters of the decorrelation filter has been optimized to: (1) balance the noise reduction and the geophysical signal attenuation produced by the filtering process; (2) minimize the differences between GRACE and model-based data and (3) maximize the ratio of variability between continents and oceans. The Gaussian filter has been optimized following the latter criteria. Besides, an anisotropic filter, the fan filter, has been analyzed as an alternative to the Gauss filter, producing better statistics.

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: Cadmium is a widespread environmental pollutant and poses some potential risks to human health. However, the signaling events controlling cadmium toxicity are not fully understood. In this study, we examined the effect of cadmium chloride on cell viability and the intracellular nitric oxide (NO) level in yeast cells. The results showed that exposure of yeast cells to cadmium (0–100 μM) could induce cell killing with significantly increased intracellular NO levels. Morphological analysis of the nuclei with 4 ' ,6-diamidino-2-phenylindole staining and DNA strand breaks analysis showed that cadmium at 50 μM can induce cell apoptosis in yeast cells. Treatment of yeast cells with cadmium (50 μM) and the nitric oxide scavenger c-PTIO [2-(4-carboxyphenyl)-4,4,5,5-teramethylimidazoline-1-oxyl-3-oxide; 0.2 mM] showed that c-PTIO attenuated the cadmium-induced cell killing. Our findings indicated that cadmium-induced yeast cell killing is mediated by a directly increased intracellular NO level.

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: Effective uptake of fermentable substrates is a fundamentally important aspect of any fermentation process. The solventogenic bacterium Clostridium beijerinckii is noted for its ability to ferment a wide range of carbohydrates, yet few of its sugar transport systems have been characterized. In common with other anaerobes, C. beijerinckii shows a marked dependence on the PEP-dependent phosphotransferase system (PTS) for sugar accumulation. In this study, the gene cbe0751 encoding the sugar-specific domains of a phosphotransferase belonging to the glucose family was cloned into an Escherichia coli strain lacking the ability to take up and phosphorylate glucose. Transformants gained ability to ferment glucose, and also mannose, and further analysis of a selected transformant demonstrated that it could take up and phosphorylate glucose, confirming that cbe0751 encodes a glucose PTS which also recognizes mannose as a substrate. RT-PCR analysis showed that cbe0751 was expressed in cultures grown on both substrates, but also to varying extents during growth on some other carbon sources. Although analogue inhibition studies suggested that Cbe0751 is not the only glucose PTS in C. beijerinckii , this system should nevertheless be regarded as a potential target for metabolic engineering to generate a strain showing improved sugar fermentation properties.

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: Heavy metals are dense chemicals with dual biological role as micronutrients and intoxicants. A few hypersaline environmental systems are naturally enriched with heavy metals, while most metal-contaminated sites are a consequence of human activities. Numerous halotolerant and moderately halophilic Bacteria possess metal tolerance, whereas a few archaeal counterparts share similar features. The main mechanisms underlying heavy metal resistance in halophilic Bacteria and Archaea include extracellular metal sequestration by biopolymers, metal efflux mediated by specific transporters and enzymatic detoxification. Biotransformation of metals by halophiles has implications both for trace metal turnover in natural saline ecosystems and for development of novel bioremediation strategies.

Description: We screened for a gene that inhibits streptomycin production in Streptomyces griseus when it is introduced on a high-copy-number plasmid pIJ702, and obtained a plasmid pKM545. The introduction of pKM545 abolished streptomycin production on all media tested including YMP-sugar and Nutrient broth. S1 protection analysis demonstrated that the introduction of this plasmid downregulated the transcriptional activity of the promoter preceding strR , the pathway-specific transcriptional regulator for streptomycin biosynthesis. The 2.8-kb Bam HI fragment cloned onto pKM545 contained two coding sequences SGR_5442 and 5443. These coding sequences and the two downstream ones (SGR_5444 and 5445) constituted a possible operon structure designated to be rspABCD (regulation of streptomycin production). RspB and RspC exhibited a marked similarity with an ATP-binding domain and a membrane-associating domain of an ABC-2 type transporter, respectively, suggesting that the Rsp proteins comprise a membrane exporter. The gene cluster consisting of the rsp operon and the upstream divergent small coding sequence (SGR_5441) was widely distributed to Streptomyces genome. An rspB mutant of S. griseus produced 3-fold streptomycin of the parental strain in YMP liquid medium. The evidence implies that the Rsp translocator is involved in the export of a substance that specifies the expression level of streptomycin biosynthesis genes in S. griseus .

Description: We previously demonstrated that in Streptomyces coelicolor two-component system AfsQ1/Q2 activates the production of the yellow-colored coelimycin P2 (also named as yCPK) on glutamate-supplemented minimal medium, and the response regulator AfsQ1 could specifically bind to the intergenic region between two structural genes, cpkA and cpkD . Here, a more in-depth investigation was performed to elucidate the mechanism underlying the role of AfsQ1/Q2 in regulating coelimycin P2 biosynthesis. Deletion of afsQ1/Q2 resulted in markedly decreased expression of the whole coelimycin P2 biosynthetic gene cluster. Electrophoretic mobility shift assays revealed that AfsQ1 bound only to the target site identified previously, but not to any other promoters in the gene cluster. Mutations of AfsQ1-binding motif only resulted in drastically reduced transcription of the cpkA/B/C operon (encoding three type I polyketide synthases) and intriguingly, led to enhanced expression of some coelimcyin P2 genes, particularly accA1 and scF . These results suggested the direct role of AfsQ1/Q2 in regulating coelimycin production, which is directly mediated by the structural genes, but not the cluster-situated regulatory genes, and also implied that other unknown mechanisms may be involved in AfsQ1/Q2-mediated regulation of coelimycin P2 biosynthesis.

Description: Sulfate-reducing bacteria (SRBs) gain their energy by coupling the oxidation of organic substrate to the reduction of sulfate to sulfide. Several SRBs are able to use alternative terminal electron acceptors to sulfate such as nitrate. Nitrate-reducing SRBs have been isolated from a diverse range of environments. In order to be able to understand the significance of nitrate reduction in SRBs, we need to examine the ecology and physiology of the nitrate-reducing SRB isolates.

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 ADP-ribosylating enzymes are encoded in many pathogenic bacteria in order to affect essential functions of the host. In this study, we show that Neisseria gonorrhoeae possess a locus that corresponds to the ADP-ribosyltransferase NarE, a previously characterized enzyme in N. meningitidis . The 291 bp coding sequence of gonococcal narE shares 100% identity with part of the coding sequence of the meningococcal narE gene due to a frameshift previously described, thus leading to a 49-amino-acid deletion at the N-terminus of gonococcal NarE protein. However, we found a promoter region and a GTG start codon, which allowed expression of the protein as demonstrated by RT-PCR and western blot analyses. Using a gonococcal NarE–6xHis fusion protein, we demonstrated that the gonococcal enzyme underwent auto-ADP-ribosylation but to a lower extent than meningococcal NarE. We also observed that gonoccocal NarE exhibited ADP-ribosyltransferase activity using agmatine and cell-free host proteins as ADP-ribose acceptors, but its activity was inhibited by human β-defensins. Taken together, our results showed that NarE of Neisseria gonorrhoeae is a functional enzyme that possesses key features of bacterial ADP-ribosylating enzymes.

Description: The effect of fructose 1,6-bisphosphate (Fru 1,6-P 2 ) on the regulatory enzymes of pentose phosphate pathway of Escherichia coli was examined. Fru 1,6-P 2 inhibited E. coli transaldolase (EC 2.2.1.2) competitively against fructose 6-phosphate and uncompetitively against erythrose 4-phosphate, whereas Fru 1,6-P 2 did not affect glucose 6-phosphate dehydrogenase (EC 1.1.1.49) and 6-phosphogluconate dehydrogenase (EC 1.1.1.44). Kinetic results can be explained by assuming that transaldolase has two kinds of binding sites for Fru 1,6-P 2 : a competitive binding site for fructose 6-phosphate and a second binding site on the enzyme-erythrose 4-phosphate complex. Fru 1,6-P 2 increased resulting from the stimulation of glycolysis, can inhibit transaldolase and further participates in the elevation of the concentration of ribose 5-phosphate that can be preferentially utilized for anabolic reaction in exponential phase of E. coli .

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 Hdr (heterodisulfide reductase)-like enzyme is predicted, from gene transcript profiling experiments previously published, to be essential in oxidative sulfur metabolism in a number of bacteria and archaea. Nevertheless, no biochemical and physicochemical data are available so far about this enzyme. Genes coding for it were identified in Aquifex aeolicus , a Gram-negative, hyperthermophilic, chemolithoautotrophic and microaerophilic bacterium that uses inorganic sulfur compounds as electron donor to grow. We provide biochemical evidence that this Hdr-like enzyme is present in this sulfur-oxidizing prokaryote (cultivated with thiosulfate or elemental sulfur). We demonstrate, by immunolocalization and cell fractionation, that Hdr-like enzyme is associated, presumably monotopically, with the membrane fraction. We show by co-immunoprecipitation assay or partial purification, that the Hdr proteins form a stable complex composed of at least five subunits, HdrA, HdrB1, HdrB2, HdrC1 and HdrC2, present in two forms of high molecular mass on native gel (~240 and 450 kDa). These studies allow us to propose a revised model for dissimilatory sulfur oxidation pathways in A. aeolicus , with Hdr predicted to generate sulfite.

Description: Geobacter sulfurreducens is an anaerobic soil bacterium that is involved in biogeochemical cycles of elements such as Fe and Mn. Although significant progress has been made in the understanding of the electron transfer processes in G. sulfurreducens , little is known about the regulatory mechanisms involved in their control. To expand the study of gene regulation in G. sulfurreducens , we carried out a genome-wide identification of transcription start sites (TSS) by 5'RACE and by deep RNA sequencing of primary mRNAs in two growth conditions. TSSs were identified along G. sulfurreducens genome and over 50% of them were located in the upstream region of the associated gene, and in some cases we detected genes with more than one TSS. Our global mapping of TSSs contributes with valuable information, which is needed for the study of transcript structure and transcription regulation signals and can ultimately contribute to the understanding of transcription initiation phenomena in G. sulfurreducens .

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: DnpA, a putative de- N -acetylase of the PIG-L superfamily, is required for antibiotic tolerance in Pseudomonas aeruginosa . Exactly how dnpA (gene locus PA5002) directs the formation of antibiotic-tolerant persister cells is currently unknown. Previous research provided evidence for a role in surface-associated process(es), possibly in lipopolysaccharide biosynthesis. In silico sequence analysis of DnpA predicts a single transmembrane domain and N in /C out orientation of DnpA. In contrast, we here show that DnpA is an integral inner membrane protein containing two transmembrane domains, with the major C-terminal part located at the cytoplasmic face. Correct insertion into the inner membrane is necessary for DnpA to promote fluoroquinolone tolerance. The membrane localization of DnpA further supports its role in cell envelope-associated process(es). In addition to shedding light on the biological role of DnpA, this study highlights the risks of overreliance on the predictive value of bioinformatics tools and the importance of rigorous experimental validation of in silico predictions.

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: 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: Cell wall impermeability and active efflux of drugs are among the primary reasons for drug resistance in mycobacteria. Efflux pumps are tripartite membrane localized transport proteins that expel drug molecules outside the cells. Several of such efflux pumps are annotated in mycobacteria, but few have been characterized, like MSMEG_2991, a putative efflux pump permease of Mycobacterium smegmatis . To substantiate this, we overexpressed MSMEG_2991 protein in Escherichia coli 2443. Expression of MSMEG_2991 elevated the resistance towards structurally unrelated groups of antibiotics. An active antibiotic efflux pump nature of MSMEG_2991 was revealed by assessing the acquisition of ciprofloxacin in the absence and presence of the efflux pump inhibitor, carbonyl cyanide m-chlorophenyl hydrazone, indicating the involvement of proton-motive force (pmf) during the efflux activity. MSMEG_2991 expression elevated biofilm formation in E. coli by 4-fold, keeping parity to some of the earlier reported efflux pumps. In silico analysis suggested the presence of 12 transmembrane helices in MSMEG_2991 resembling EmrD efflux pump of E. coli . Based on in vivo and in silico analyses, MSMEG_2991 may be designated as a pmf-mediated multidrug efflux pump protein that expels diverse groups of antibiotics and might as well be involved in the biofilm enhancement.

Description: Bacterial small RNAs (sRNAs) play essential roles in the post-transcriptional control of gene expression. To improve their detection by conventional microarrays, we designed a custom microarray containing a group of probes targeting known and some putative Escherichia coli sRNAs. To assess its potential in detection of sRNAs, RNA profiling experiments were performed with total RNA extracted from E. coli MG1655 cells exponentially grown in rich (Luria–Bertani) and minimal (M9/glucose) media. We found that many sRNAs could yield reasonably strong and statistically significant signals corresponding to nearly all sRNAs annotated in the EcoCyc database. Besides differential expression of two sRNAs (GcvB and RydB), expression of other sRNAs was less affected by the composition of the growth media. Other examples of the differentially expressed sRNAs were revealed by comparing gene expression of the wild-type strain and its isogenic mutant lacking functional poly(A) polymerase I ( pcnB ). Further, northern blot analysis was employed to validate these data and to assess the existence of new putative sRNAs. Our results suggest that the use of custom microarrays with improved capacities for detection of sRNAs can offer an attractive opportunity for efficient gene expression profiling of sRNAs and their target mRNAs at the whole transcriptome level.

Description: Escherichia coli DedA/Tvp38 family proteins YghB and YqjA are putative membrane transporters with 62% amino acid identity and overlapping functions. An E. coli strain (BC202) with nonpolar yghB and yqjA mutations displays cell-division defects and temperature sensitivity and is sensitive to antibiotics and alkaline pH. In this study, we performed site-directed mutagenesis on conserved, charged amino acids of YqjA and YghB. We discovered two conserved predicted membrane-embedded arginines (R130 and R136) that are critical for function in both proteins as defined by their ability to complement BC202 phenotypes, when expressed from a plasmid. Lysine can substitute for arginine at position R130 indicating a charge dependence at this position, but could not substitute at R136. In light of the established role that arginine plays in the translocation mechanism of numerous membrane transporters, we hypothesize that these amino acids play a role in the transport mechanism of these DedA/Tvp38 family proteins.

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: Candida albicans is a dimorphic yeast commonly found on human mucosal membranes that switches from yeast to hyphal morphology in response to environmental factors. The change to hyphal growth requires histone H3 modifications by the yeast-specific histone acetyltransferase Rtt109. In addition to its role in morphogenesis, Rtt109-dependent acetylation of histone H3 lysine residues 9 and 56 has regulatory functions during DNA replication and repair. Boric acid (BA) is a broad-spectrum agent that specifically inhibits C. albicans hyphal growth, locking the fungus in its harmless commensal yeast state. The present study characterizes the effect of BA on C. albicans histone acetylation in respect to specificity, time-course and significance. We demonstrate that sublethal concentrations of BA reduce H3K9/H3K56 acetylation, both on a basal level and in response to genotoxic stress. Acetylation at other selected histone sites were not affected by BA. qRT-PCR expression analysis of the DNA repair gene Rad51 indicated no elevated level of genotoxic stress during BA exposure. A forward-mutation analysis demonstrated the BA does not increase spontaneous or induced mutations . The findings suggest that DNA repair remains effective even when histone H3 acetylation decreases and dispels the notion that BA treatment impairs genome integrity in yeast.

Description: The impact of gene mutations on the growth of the cells can be studied using pure cultures. However, the importance of certain proteins and pathways can be also examined via co-culturing wild type and its mutant derivative. Here, the relative fitness of a mutant strain that lacks the global nitrogen regulator, CodY, was examined in Bacillus cereus , a food poisoning Gram-positive bacterium. Fitness measurements revealed that the codY strain was outcompeted when cocultured with the wild-type ATCC 14579 under various rich laboratory medium, and also when inoculated in certain beverages. In nutrient-poor minimal medium, the codY mutant had comparable fitness to the wild-type strain. Interestingly, the relative fitness of the codY strain was antagonistic when it was cultivated in apple or orange juices due to unknown properties of these beverages, highlighting the importance of chemical composition of the test medium during the bacterial fitness measurements.

Description: Out of 8000 candidates from a genetic screening for Pseudomonas putida KT2442 mutants showing defects in biofilm formation, 40 independent mutants with diminished levels of biofilm were analyzed. Most of these mutants carried insertions in genes of the lap cluster, whose products are responsible for synthesis, export and degradation of the adhesin LapA. All mutants in this class were strongly defective in biofilm formation. Mutants in the flagellar regulatory genes fleQ and flhF showed similar defects to that of the lap mutants . On the contrary, transposon insertions in the flagellar structural genes fliP and flgG , that also impair flagellar motility, had a modest defect in biofilm formation. A mutation in gacS , encoding the sensor element of the GacS/GacA two-component system, also had a moderate effect on biofilm formation. Additional insertions targeted genes involved in cell envelope function: PP3222, encoding the permease element of an ABC-type transporter and tolB , encoding the periplasmic component of the Tol-OprL system required for outer membrane stability. Our results underscore the central role of LapA, suggest cross-regulation between motility and adhesion functions and provide insights on the role of cell envelope trafficking and maintenance for biofilm development in P. putida .

Description: Horizontal gene transfer plays a significant role in bacterial evolution and has major clinical importance. Thus, it is vital to understand the mechanisms and kinetics of genetic transformations. Natural transformation is the driving mechanism for horizontal gene transfer in diverse genera of bacteria. Our study introduces a simple and rapid method for the investigation of natural transformation. This highly sensitive system allows the detection of a transformation event directly from a bacterial population without any separation step or selection of cells. The system is based on the bacterial luciferase operon from Photorhabdus luminescens . The studied molecular tools consist of the functional modules luxCDE and luxAB , which involve a replicative plasmid and an integrative gene cassette. A well-established host for bacterial genetic investigations, Acinetobacter baylyi ADP1, is used as the model bacterium. We show that natural transformation followed by homologous recombination or plasmid recircularization can be readily detected in both actively growing and static biofilm-like cultures, including very rare transformation events. The system allows the detection of natural transformation within 1 h of introducing sample DNA into the culture. The introduced method provides a convenient means to study the kinetics of natural transformation under variable conditions and perturbations.

Description: The legume–rhizobial interaction results in the formation of symbiotic nodules in which rhizobia fix nitrogen. During the process of symbiosis, reactive oxygen species (ROS) are generated. Thus, the response of rhizobia to ROS is important for successful nodulation and nitrogen fixation. In this study, we investigated how Azorhizobium caulinodans , a rhizobium that forms both root and stem nodules on its host plant, regulates ROS resistance. We found that in-frame deletions of a gene encoding the putative catalase-peroxidase katG or a gene encoding a LysR-family regulatory protein, oxyR , exhibited increased sensitivity to H 2 O 2 . We then showed that OxyR positively regulated katG expression in an H 2 O 2 -independent fashion. Furthermore, we found that deletion of katG or oxyR led to significant reduction in the number of stem nodules and decrease of nitrogen fixation capacities in symbiosis. Our results revealed that KatG and OxyR are not only critical for antioxidant defense in vitro , but also important for nodule formation and nitrogen fixation during interaction with plant hosts.

Description: Peptidoglycan (also known as murein) is an important envelope component of bacteria, and its turnover usually takes place at considerable levels during normal growth. Amino sugars and murein tripeptide resulting from murein degradation are used for resynthesis of peptidoglycan or as self-generated nutrients or energy sources for cell growth. PgrR (regulator of peptide glycan recycling; formerly YcjZ) was recently identified as a repressor of several genes participating in uptake and degradation of murein tripeptide. In this study, we identified the ycjG gene involved in murein tripeptide degradation as a new direct target of PgrR. The expression of PgrR-regulated genes including ycjY , mppA , mpaA and ycjG was repressed in the presence of a good nitrogen source, but their expression increased under poor nitrogen conditions. Under nitrogen starvation, the pgrR mutant cells exhibited faster growth than wild-type cells, implying that derepression of genes under the control of PgrR may help cells overcome nitrogen limitation. Therefore, these results suggest that nitrogen starvation induces derepression of PgrR-controlled genes involved in uptake and degradation of murein tripeptide, and this may stimulate the utilization of murein tripeptide as a nitrogen source.

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: 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 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: The central metabolite acetyl phosphate (acP) has long been proposed to influence transcription regulation by directly transferring its phosphoryl group to a number of response regulators in many bacterial species. Here, we provide in vitro evidence for this proposition and demonstrate, using an in vitro transcription system, that acP-dependent phosphorylation of aspartate 51 of CpxR induces transcription of one of its regulon members in E. coli , cpxP . We also used this in vitro transcription system to extend our previously reported in vivo data that hypothesized that acetylation of RNA polymerase (RNAP) influences acP-dependent cpxP transcription, using glutamine as a genetic mimic for acetylated arginine 291 of the carboxy-terminal domain of RNAP α subunit. The data we present here lend strong support to the hypothesis that acP has a direct effect on transcription regulation in E. coli via phosphorylation of CpxR, and that RNAP acetylation can modulate this response.

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: Chaperones are a diverse class of ubiquitous proteins that assist other cellular proteins in folding correctly and maintaining their native structure. Many different chaperones cooperate to constitute the ‘proteostasis’ machinery in the cells. It has been proposed earlier that archaeal organisms could be ideal model systems for deciphering the basic functioning of the ‘protein folding machinery’ in higher eukaryotes. Several chaperone families have been characterized in archaea over the years but mostly one protein at a time, making it difficult to decipher the composition and mechanistics of the protein folding system as a whole. In order to deal with these lacunae, we have developed a database of all archaeal chaperone proteins, CrAgDb ( C haperone r epertoire in A rchaeal g enomes). The data have been presented in a systematic way with intuitive browse and search facilities for easy retrieval of information. Access to these curated datasets should expedite large-scale analysis of archaeal chaperone networks and significantly advance our understanding of operation and regulation of the protein folding machinery in archaea. Researchers could then translate this knowledge to comprehend the more complex protein folding pathways in eukaryotic systems. The database is freely available at http://14.139.227.92/mkumar/cragdb/ .

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: In bacteria, cysteine can be synthesized from serine by two steps involving an L-serine O-acetyltransferase (SAT) and a cysteine synthase (CysK). While CysK is found in the publicly available annotated genome from Lactobacillus casei ATCC 334, a gene encoding SAT ( cysE ) is missing. In this study, we found that various strains of L. casei grew in a chemically defined medium containing sulfide as the sole sulfur source, indicating the presence of a serine O-acetyltransferase. The gene lying upstream of cysK is predicted to encode a homoserine trans-succinylase ( metA ). To study the function of this gene, it was cloned from L. casei FAM18110. The purified, recombinant protein did not acylate L-homoserine in vitro . Instead, it catalyzed the formation of O-acetyl serine from L-serine and acetyl-CoA. Furthermore, the plasmid expressing the L. casei gene complemented an Escherichia coli cysE mutant strain but not an E. coli metA mutant. This clearly demonstrated that the gene annotated as metA in fact encodes the SAT function and should be annotated as cysE .

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: The paper in question by Van Camp and co-authors [MVC] challenges previous work showing that ground gravity data arising from hydrology can provide a consistent signal for the comparison with satellite gravity data. The data sets used are similar to those used previously, that is, the gravity field as measured by the GRACE satellites versus ground-based data from superconducting gravimeters (SGs) over the same continental area, in this case Central Europe. One of the main impediments in this paper is the presentation that is frequently confusing and misleading as to what the data analysis really shows, for example, the irregular treatment of annual components that are first subtracted then reappear in the analysis. More importantly, we disagree on specific points. Two calculations are included in our comment to illustrate where we believe that the processing in [MVC] paper is deficient. The first deals with their erroneous treatment of the global hydrology using a truncated spherical harmonic approach which explains almost a factor 2 error in their computation of the loading. The second shows the effect of making the wrong assumption in the GRACE/hydrology/surface gravity comparison by inverting the whole of the hydrology loading for underground stations. We also challenge their claims that empirical orthogonal function techniques cannot be done in the presence of periodic components, and that SG data cannot be corrected for comparisons with GRACE data. The main conclusion of their paper, that there is little coherence between ground gravity stations and this invalidates GRACE comparisons, is therefore questionable. There is nothing in [MVC] that contradicts any of the previous papers that have shown clearly a strong relation between seasonal signals obtained from both ground gravity and GRACE satellite data.

Description: The influence of changes in surface ice-mass redistribution and associated viscoelastic response of the Earth, known as glacial isostatic adjustment (GIA), on the Earth's rotational dynamics has long been known. Equally important is the effect of the changes in the rotational dynamics on the viscoelastic deformation of the Earth. This signal, known as the rotational feedback, or more precisely, the rotational feedback on the sea level equation, has been mathematically described by the sea level equation extended for the term that is proportional to perturbation in the centrifugal potential and the second-degree tidal Love number. The perturbation in the centrifugal force due to changes in the Earth's rotational dynamics enters not only into the sea level equation, but also into the conservation law of linear momentum such that the internal viscoelastic force, the perturbation in the gravitational force and the perturbation in the centrifugal force are in balance. Adding the centrifugal-force perturbation to the linear-momentum balance creates an additional rotational feedback on the viscoelastic deformations of the Earth. We term this feedback mechanism, which is studied in this paper, as the rotational feedback on the linear-momentum balance. We extend both the time-domain method for modelling the GIA response of laterally heterogeneous earth models developed by Martinec and the traditional Laplace-domain method for modelling the GIA-induced rotational response to surface loading by considering the rotational feedback on linear-momentum balance. The correctness of the mathematical extensions of the methods is validated numerically by comparing the polar-motion response to the GIA process and the rotationally induced degree 2 and order 1 spherical harmonic component of the surface vertical displacement and gravity field. We present the difference between the case where the rotational feedback on linear-momentum balance is considered against that where it is not. Numerical simulations show that the resulting difference in radial displacement and sea level change between these situations since the Last Glacial Maximum reaches values of ±25 and ±1.8 m, respectively. Furthermore, the surface deformation pattern is modified by up to 10 per cent in areas of former or ongoing glaciation, but by up to 50 per cent at the bottom of the southern Indian ocean. This also results in the movement of coastlines during the last deglaciation to differ between the two cases due to the difference in the ocean loading, which is seen for instance in the area around Hudson Bay, Canada and along the Chinese, Australian or Argentinian coastlines.

Description: During megathrust earthquakes, great ruptures are accompanied by large scale mass redistribution inside the solid Earth and by ocean mass redistribution due to bathymetry changes. These large scale mass displacements can be detected using the monthly gravity maps of the GRACE satellite mission. In recent years it has become increasingly common to use the long wavelength changes in the Earth's gravity field observed by GRACE to infer seismic source properties for large megathrust earthquakes. An important advantage of space gravimetry is that it is independent from the availability of land for its measurements. This is relevant for observation of megathrust earthquakes, which occur mostly offshore, such as the $M_{\text{w}} \sim 9$ 2004 Sumatra–Andaman, 2010 Maule (Chile) and 2011 Tohoku-Oki (Japan) events. In Broerse et al. , we examined the effect of the presence of an ocean above the rupture on long wavelength gravity changes and showed it to be of the first order. Here we revisit the implementation of an ocean layer through the sea level equation and compare the results with approximated methods that have been used in the literature. One of the simplifications usually lies in the assumption of a globally uniform ocean layer. We show that especially in the case of the 2010 Maule earthquake, due to the closeness of the South American continent, the uniform ocean assumption is not valid and causes errors up to 57 per cent for modelled peak geoid height changes (expressed at a spherical harmonic truncation degree of 40). In addition, we show that when a large amount of slip occurs close to the trench, horizontal motions of the ocean floor play a mayor role in the ocean contribution to gravity changes. Using a slip model of the 2011 Tohoku-Oki earthquake that places the majority of slip close to the surface, the peak value in geoid height change increases by 50 per cent due to horizontal ocean floor motion. Furthermore, we test the influence of the maximum spherical harmonic degree at which the sea level equation is performed for sea level changes occurring along coastlines, which shows to be important for relative sea level changes occurring along the shore. Finally, we demonstrate that ocean floor loading, self-gravitation of water and conservation of water mass are of second order importance for coseismic gravity changes. When GRACE observations are used to determine earthquake parameters such as seismic moment or source depth, the uniform ocean layer method introduces large biases, depending on the location of the rupture with respect to the continent. The same holds for interpreting shallow slip when horizontal motions are not properly accounted for in the ocean contribution. In both cases the depth at which slip occurs will be underestimated.

Description: Long-term volcanic subsidence provides insight into intereruptive processes, which comprise the longest portion of the eruptive cycle. Ground-based geodetic surveys of Medicine Lake Volcano (MLV), northern CA, document subsidence at rates of ~–10 mm yr –1 between 1954 and 2004. The long observation period plus the duration and stable magnitude of this signal presents an ideal opportunity to study long-term volcanic deformation, but this first requires accurate knowledge of the geometry and magnitude of the source. Best-fitting analytical source models to past levelling and GPS data sets show conflicting source parameters—primarily the model depth. To overcome this, we combine multiple tracks of InSAR data, each with a different look angle, to improve upon the spatial resolution of ground-based measurements. We compare the results from InSAR to those of past geodetic studies, extending the geodetic record to 2011 and demonstrating that subsidence at MLV continues at ~–10 mm yr –1 . Using geophysical inversions, we obtain the best-fitting analytical source model—a sill located at 9–10 km depth beneath the caldera. This model geometry is similar to those of past studies, providing a good fit to the high spatial density of InSAR measurements, while accounting for the high ratio of vertical to horizontal deformation derived from InSAR and recorded by existing levelling and GPS data sets. We discuss possible causes of subsidence and show that this model supports the hypothesis that deformation at MLV is driven by tectonic extension, gravitational loading, plus a component of volume loss at depth, most likely due to cooling and crystallization within the intrusive complex that underlies the edifice. Past InSAR surveys at MLV, and throughout the Cascades, are of variable success due to dense vegetation, snow cover and atmospheric artefacts. In this study, we demonstrate how InSAR may be successfully used in this setting by applying a suite of multitemporal analysis methods that account for atmospheric and orbital noise sources. These methods include: a stacking strategy based upon the noise characteristics of each data set; pixelwise rate-map formation (-RATE) and persistent scatterer InSAR (StaMPS).

Description: In the literature, the inverted coseismic slip models from seismological and geodetic data for the 2011 Tohoku-Oki earthquake portray significant discrepancies, in particular regarding the intensity and the distribution of the rupture near the trench. For a megathrust earthquake, it is difficult to discern the slip along the shallow part of the fault from the geodetic data, which are often acquired on land. In this paper, we discuss the uncertainties in the slip distribution inversion using the geodetic data for the 2011 Tohoku earthquake and the Fully Bayesian Inversion method. These uncertainties are due to the prior information regarding the boundary conditions at the edges of the fault, the dip subduction angle and the smoothing operator. Using continuous GPS data from the Japan Island, the results for the rigid and free boundary conditions show that they produce remarkably different slip distributions at shallow depths, with the latter producing a large slip exceeding 30 m near the surface. These results indicate that the smoothing operator (gradient or Laplacian schemes) does not severely affect the slip pattern. To better invert the coseismic slip, we then introduce the ocean bottom GPS (OB-GPS) data, which improve the resolution of the shallow part of the fault. We obtain a near-trench slip greater than 40 m that reaches the Earth's surface, regardless of which boundary condition is used. Additionally, we show that using a mean dip angle for the fault as derived from subduction models is adequate if the goal is to invert for the general features of the slip pattern of this megathrust event.

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: Complications arise in the interpretation of gravity fields because of interference from systematic degradations, such as boundary blurring and distortion. The major sources of these degradations are the various systematic errors that inevitably occur during gravity field data acquisition, discretization and geophysical forward modelling. To address this problem, we evaluate deconvolution method that aim to detect the clear horizontal boundaries of anomalous sources by the suppression of systematic errors. A convolution-based multilayer projection model, based on the classical 3-D gravity field forward model, is innovatively derived to model the systematic error degradation. Our deconvolution algorithm is specifically designed based on this multilayer projection model, in which three types of systematic error are defined. The degradations of the different systematic errors are considered in the deconvolution algorithm. As the primary source of degradation, the convolution-based systematic error is the main object of the multilayer projection model. Both the random systematic error and the projection systematic error are shown to form an integral part of the multilayer projection model, and the mixed norm regularization method and the primal-dual optimization method are therefore employed to control these errors and stabilize the deconvolution solution. We herein analyse the parameter identification and convergence of the proposed algorithms, and synthetic and field data sets are both used to illustrate their effectiveness. Additional synthetic examples are specifically designed to analyse the effects of the projection systematic error, which is caused by the uncertainty associated with the estimation of the impulse response function.

Description: We propose to test if gravimetry can prove useful in discriminating different models of long-term deep crustal processes in the case of the Taiwan mountain belt. We discuss two existing tectonic models that differ in the deep processes proposed to sustain the long-term growth of the orogen. One model assumes underplating of the uppermost Eurasian crust with subduction of the deeper part of the crust into the mantle. The other one suggests the accretion of the whole Eurasian crust above crustal-scale ramps, the lower crust being accreted into the collisional orogen. We compute the temporal gravity changes caused only by long-term rock mass transfers at depth for each of them. We show that the underplating model implies a rate of gravity change of –6 x 10 –2 μGal yr –1 , a value that increases to 2 x 10 –2 μGal yr –1 if crustal subduction is neglected. If the accretion of the whole Eurasian crust occurs, a rate of 7 x 10 –2 μGal yr –1 is obtained. The two models tested differ both in signal amplitude and spatial distribution. The yearly gravity changes expected by long-term deep crustal mass processes in Taiwan are two orders of magnitude below the present-day uncertainty of land-based gravity measurements. Assuming that these annually averaged long-term gravity changes will linearly accumulate with ongoing mountain building, multidecadal time-series are needed to identify comparable rates of gravity change. However, as gravity is sensitive to any mass redistribution, effects of short-term processes such as seismicity and surface mass transfers (erosion, sedimentation, ground-water) may prevent from detecting any long-term deep signal. This study indicates that temporal gravity is not appropriate for deciphering the long-term deep crustal processes involved in the Taiwan mountain belt.

Description: The computation of quasi-static deformation for axisymmetric viscoelastic structures on a gravitating spherical earth is addressed using the spectral element method (SEM). A 2-D spectral element domain is defined with respect to spherical coordinates of radius and angular distance from a pole of symmetry, and 3-D viscoelastic structure is assumed to be azimuthally symmetric with respect to this pole. A point dislocation source that is periodic in azimuth is implemented with a truncated sequence of azimuthal order numbers. Viscoelasticity is limited to linear rheologies and is implemented with the correspondence principle in the Laplace transform domain. This leads to a series of decoupled 2-D problems which are solved with the SEM. Inverse Laplace transform of the independent 2-D solutions leads to the time-domain solution of the 3-D equations of quasi-static equilibrium imposed on a 2-D structure. The numerical procedure is verified through comparison with analytic solutions for finite faults embedded in a laterally homogeneous viscoelastic structure. This methodology is applicable to situations where the predominant structure varies in one horizontal direction, such as a structural contrast across (or parallel to) a long strike-slip fault.

Description: The geodetic rates for the gravity variation and vertical uplift in polar regions subject to past and present-day ice-mass changes (PDIMCs) provide important insight into the rheological structure of the Earth. We provide an update of the rates observed at Ny-Ålesund, Svalbard. To do so, we extract and remove the significant seasonal content from the observations. The rate of gravity variations, derived from absolute and relative gravity measurements, is –1.39 ± 0.11 μGal yr –1 . The rate of vertical displacements is estimated using GPS and tide gauge measurements. We obtain 7.94 ± 0.21 and 8.29 ± 1.60 mm yr –1 , respectively. We compare the extracted signal with that predicted by GLDAS/Noah and ERA-interim hydrology models. We find that the seasonal gravity variations are well-represented by local hydrology changes contained in the ERA-interim model. The phase of seasonal vertical displacements are due to non-local continental hydrology and non-tidal ocean loading. However, a large part of the amplitude of the seasonal vertical displacements remains unexplained. The geodetic rates are used to investigate the asthenosphere viscosity and lithosphere/asthenosphere thicknesses. We first correct the updated geodetic rates for those induced by PDIMCs in Svalbard, using published results, and the sea level change due to the melting of the major ice reservoirs. We show that the latter are at the level of the geodetic rate uncertainties and are responsible for rates of gravity variations and vertical displacements of –0.29 ± 0.03 μGal yr –1 and 1.11 ± 0.10 mm yr –1 , respectively. To account for the late Pleistocene deglaciation, we use the global ice evolution model ICE-3G. The Little Ice Age (LIA) deglaciation in Svalbard is modelled using a disc load model with a simple linear temporal evolution. The geodetic rates at Ny-Ålesund induced by the past deglaciations depend on the viscosity structure of the Earth. We find that viscous relaxation time due to the LIA deglaciation in Svalbard is more than 60 times shorter than that due to the Pleistocene deglaciation. We also find that the response to past and PDIMCs of an Earth model with asthenosphere viscosities ranging between 1.0 and 5.5 x 10 18 Pa s and lithosphere (resp. asthenosphere) thicknesses ranging between 50 and 100 km (resp. 120 and 170 km) can explain the rates derived from geodetic observations.

Description: The 2 principle and the unbiased predictive risk estimator are used to determine optimal regularization parameters in the context of 3-D focusing gravity inversion with the minimum support stabilizer. At each iteration of the focusing inversion the minimum support stabilizer is determined and then the fidelity term is updated using the standard form transformation. Solution of the resulting Tikhonov functional is found efficiently using the singular value decomposition of the transformed model matrix, which also provides for efficient determination of the updated regularization parameter each step. Experimental 3-D simulations using synthetic data of a dipping dike and a cube anomaly demonstrate that both parameter estimation techniques outperform the Morozov discrepancy principle for determining the regularization parameter. Smaller relative errors of the reconstructed models are obtained with fewer iterations. Data acquired over the Gotvand dam site in the south-west of Iran are used to validate use of the methods for inversion of practical data and provide good estimates of anomalous structures within the subsurface.

Description: Global navigation satellite systems (GNSSs) have revealed that a mega-thrust earthquake that occurs in an island-arc trench system causes post-seismic crustal deformation. Such crustal deformation data have been interpreted by combining three mechanisms: afterslip, poroelastic rebound and viscoelastic relaxation. It is seismologically important to determine the contribution of each mechanism because it provides frictional properties between the plate boundaries and viscosity estimates in the asthenosphere which are necessary to evaluate the stress behaviour during earthquake cycles. However, the observation sites of GNSS are mostly deployed over land and can detect only a small part of the large-scale deformation, which precludes a clear separation of the mechanisms. To extend the spatial coverage of the deformation area, recent studies started to use satellite gravity data that can detect long-wavelength deformations over the ocean. To date, compared with theoretical models for calculating the post-seismic crustal deformation, a few models have been proposed to interpret the corresponding gravity variations. Previous approaches have adopted approximations for the effects of compressibility, sphericity and self-gravitation when computing gravity changes. In this study, a new spectral-finite element approach is presented to consider the effects of material compressibility for Burgers viscoelastic earth model with a laterally heterogeneous viscosity distribution. After the basic principles are explained, it is applied to the 2004 Sumatra–Andaman earthquake. For this event, post-seismic deformation mechanisms are still a controversial topic. Using the developed approach, it is shown that the spatial patterns of gravity change generated by the above three mechanisms clearly differ from one another. A comparison of the theoretical simulation results with the satellite gravity data obtained from the Gravity Recovery and Climate Experiment reveals that both afterslip and viscoelastic relaxation are occurring. Considering the spatial patterns in satellite gravity fields is an effective method for investigating post-seismic deformation mechanisms.

Description: Some of the major geothermal anomalies in central Europe are linked to tectonic structures within the top of crystalline basement, which modify strongly the top of this basement. Their assessment is a major challenge in exploration geophysics. Gravity has been proven to be suitable for the detection of mainly large scale lithological and structural inhomogeneities. Indeed, it is well known and proven by different wells that, for example, in northern Switzerland extended negative anomalies are linked to such structures. Due to depth limitation of wells, there vertical extension is often unknown. In this study, we have investigated the potential of gravity for the geometrical characterization of such basement structures. Our approach consists in the combination of the series of Butterworth filters, geological modelling and best-fitting between observed and computed residual anomalies. In this respect, filters of variable wavelength are applied to observed and computed gravity data. The geological model is discretized into a finite element mesh. Near-surface anomalies and the effect of the sedimentary cover were eliminated using cut-off wavelength of 10 km and geological and seismic information. We analysed the potential of preferential Butterworth filtering in a sensitivity study and applied the above mentioned approach to part of the Swiss molasses basin. Sensitivity analyses reveal that such sets of residual anomalies represents a pseudo-tomography revealing the distribution of different structures with depth. This finding allows for interpreting negative anomalies in terms of 3-D volumes. Best-fitting then permits determination of the most likely 3-D geometries of such basement structures. Our model fits both, geological observations and gravity: among 10 deep boreholes in the studied area, six reach the respective units and confirm our distribution of the negative (and positive) anomalies.

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: 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: We purified a chitinase from pilei extractions of Coprinopsis cinerea fruiting bodies by ammonium sulfate precipitation and CM Sepharose cation exchange chromatography. MALDI-TOF/TOF MS analysis characterized this purified chitinase as a putative class V chitinase, ChiB1. ChiB1 hydrolyzed colloidal chitin and chitosan, whereas it did not hydrolyze chitin powder. ChiB1 cleaved only pNP-(GlcNAc) 2 , rather than pNP-GlcNAc or pNP-(Glc-NAc) 3 , to release nitrophenol. ChiB1 preferably and progressively released (GlcNAc)2 from (GlcNAc)6 and digested (GlcNAc)6 to two molecules of (GlcNAc)3 in a small proportion, but did not split (GlcNAc) 2 , so it is an exochitinase. ChiB1 has an optimum temperature range of 35°C to 40°C and an optimum pH of 5.0. ChiB1 exhibited Km and Vmax values of 2.63 mg ml –1 and 2.31 μmol min –1  mg protein –1 for colloidal chitin, respectively. The ChiB1 gene, along with another putative endochitinase (class III chitinase gene), was expressed dominantly among eight predicted chitinase genes in the genome, and its expression level increased with the maturation of fruiting bodies. ChiB1 incubation released a large amount of soluble β-glucan fractions from alkali-insoluble cell wall fractions of C. cinerea fruiting bodies, thereby it may promote the degradation of cell walls in synergy with the β-1,3-glucanases during pileus autolysis.

Description: Acidithiobacillus ferrivorans is an acidophile implicated in low-temperature biomining for the recovery of metals from sulfide minerals. Acidithiobacillus ferrivorans obtains its energy from the oxidation of inorganic sulfur compounds, and genes encoding several alternative pathways have been identified. Next-generation sequencing of At. ferrivorans RNA transcripts identified the genes coding for metabolic and electron transport proteins for energy conservation from tetrathionate as electron donor. RNA transcripts suggested that tetrathionate was hydrolyzed by the tetH1 gene product to form thiosulfate, elemental sulfur and sulfate. Despite two of the genes being truncated, RNA transcripts for the SoxXYZAB complex had higher levels than for thiosulfate quinone oxidoreductase ( doxDA genes). However, a lack of heme-binding sites in soxX suggested that DoxDA was responsible for thiosulfate metabolism. Higher RNA transcript counts also suggested that elemental sulfur was metabolized by heterodisulfide reductase ( hdr genes) rather than sulfur oxygenase reductase ( sor ). The sulfite produced as a product of heterodisulfide reductase was suggested to be oxidized by a pathway involving the sat gene product or abiotically react with elemental sulfur to form thiosulfate. Finally, several electron transport complexes were involved in energy conservation. This study has elucidated the previously unknown At. ferrivorans tetrathionate metabolic pathway that is important in biomining.

Description: TRP1 is a frequently used auxotrophic marker for genetic modifications and selections in trp – budding yeast strains, including the commonly used wild-type strain W303a. However, we found that introduction of the TRP1 gene into a trp – strain significantly affected vegetative growth at low and high temperatures. Therefore, caution should be needed when working in a trp – background strain and using the TRP1 marker to study stress response phenotypes, particularly when analyzing temperature sensitivities.

Description: The factors driving the dominance of the Calvin–Benson–Bassham cycle (CBB) or reductive citric acid cycle (rCAC) in autotrophic microorganisms in different habitats are debated. Based on costs for synthesizing a few metabolic intermediates, it has been suggested that the CBB poses a disadvantage due to higher metabolic cost. The purpose of this study was to extend this estimate of cost from metabolite synthesis to biomass synthesis. For 12 gammaproteobacteria (CBB) and five epsilonproteobacteria (rCAC), the amount of ATP to synthesize a gram of biomass from CO 2 was calculated from genome sequences via metabolic maps. The eleven central carbon metabolites needed to synthesize biomass were all less expensive to synthesize via the rCAC (66%–89% of the ATP needed to synthesize them via CBB). Differences in cell compositions did result in differing demands for metabolites among the organisms, but the differences in cost to synthesize biomass were small among organisms that used a particular pathway (e.g. rCAC), compared to the difference between pathways (rCAC versus CBB). The rCAC autotrophs averaged 0.195 moles ATP per g biomass, while their CBB counterparts averaged 0.238. This is the first in silico estimate of the relative expense of both pathways to generate biomass.

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: Bacteria carry a number of genes that cause cell growth arrest or cell lysis upon expression. Notably, defective prophages retain many lysis proteins. Here, we identified a novel lytic gene, ydfD , on the Qin prophage segment of the Escherichia coli genome. YdfD lyses 99.9% of cells within 2 h of its induction. The co-expression of the upstream gene, dicB , encoding a cell division inhibitor, as well as sulA , encoding another cell division inhibitor, abolished YdfD-induced cell lysis. These results imply that YdfD-induced lysis is a cell division-dependent event. We further found that by deleting the hydrophobic 22-residue N-terminal domain, the resulting 42-residue C-terminal domain was still toxic to cause cell lysis. We propose that YdfD, associated with the cytoplasmic membrane, inhibits an essential cellular process(s).

Description: Many bacteria are capable of accumulating intracellular granules of polyhydroxyalkanoates (PHA). In this work, we developed confocal microscopy analysis of bacterial cells to study changes in the diameters of cells as well as PHA granules during growth and PHA accumulation in the bacterium Cupriavidus necator H16 (formerly Ralstonia eutropha ). The cell envelope was stained by DiD ® fluorescent probe and PHA granules by Nile Red. Signals from both probes were separated based on their spectral and fluorescence life-time properties. During growth and PHA accumulation, bacterial cells increased their length but the width of the cells remained constant. The volume fraction of PHA granules in cells increased during PHA accumulation, nevertheless, its value did not exceed 40 vol. % regardless of the PHA weight content. It seems that bacterial cultures lengthen the cells in order to control the PHA volume portion. However, since similar changes in cell length were also observed in a PHA non-accumulating mutant, it seems that there is no direct control mechanism, which regulates the prolongation of the cells with respect to PHA granules volume. It is more likely that PHA biosynthesis and the length of cells are influenced by the same external stimuli such as nutrient limitation.

Description: Three putative regulatory genes, namely ttmRI, ttmRII and ttmRIII , which are present in the tetramycin ( ttm ) biosynthetic gene cluster, were found in Streptomyces ahygroscopicus . Disruption of ttmRI, ttmRII or ttmRIII reduced tetramycin production, and their complementation restored production to varying degrees. Gene expression analysis of the wild-type (WT) and mutant strains through reverse transcriptase–polymerase chain reaction (RT-PCR) of the ttm gene cluster showed that the expression levels of most of the biosynthetic genes were reduced in ttmRI , ttmRII and ttmRIII . Electrophoretic mobility shift assays demonstrated that TtmRI, TtmRII and TtmRIII bound the promoters of several genes in the ttm gene cluster. This study found that these three proteins are a group of positive regulators that activate the transcription of the ttm gene cluster in S. ahygroscopicus . In addition, ttmRII had a reduced degree of grey pigmentation. Thus, TtmRII has a pleiotropic regulatory function in the tetramycin biosynthetic pathway and in development.