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: Invertebrate animal species that can withstand temperatures as high as 37°C, the human body temperature, are limited. In the present study, we utilized the two-spotted cricket, Gryllus bimaculatus , which lives in tropical and subtropical regions, as an animal model of human pathogenic bacterial infection. Injection of Pseudomonas aeruginosa or Staphylococcus aureus into the hemolymph killed crickets. Injected P. aeruginosa or S. aureus proliferated in the hemolymph until the cricket died. The ability of these pathogenic bacteria to kill the crickets was blocked by the administration of antibiotics. S. aureus gene-knockout mutants of virulence factors, including cvfA, agr and srtA , exhibited decreased killing ability compared with the parent strain. The dose at which 50% of crickets were killed by P. aeruginosa or S. aureus was not decreased at 37°C compared with that at 27°C. Injection of Listeria monocytogenes , which upregulates toxin expression at 37°C, killed crickets, and the dose at which 50% of crickets were killed was decreased at 37°C compared with that at 27°C. These findings suggest that the two-spotted cricket is a useful model animal for evaluating the virulence properties of various human pathogenic bacteria at variable temperature including 37°C.

Description: The OmpA-like protein domain has been associated with peptidoglycan-binding proteins, and is often found in virulence factors of bacterial pathogens. The intracellular pathogen Legionella pneumophila encodes for six proteins that contain the OmpA-like domain, among them the highly conserved uncharacterized protein we named CmpA. Here we set out to characterize the CmpA protein and determine its contribution to intracellular survival of L. pneumophila . Secondary structure analysis suggests that CmpA is an inner membrane protein with a peptidoglycan-binding domain at the C-teminus. A cmpA mutant was able to replicate normally in broth, but failed to compete with an isogenic wild-type strain in an intracellular growth competition assay. The cmpA mutant also displayed significant intracellular growth defects in both the protozoan host Acanthamoeba castellanii and in primary bone marrow-derived macrophages, where uptake into the cells was also impaired. The cmpA phenotypes were completely restored upon expression of CmpA in trans . The data presented here establish CmpA as a novel virulence factor of L. pneumophila that is required for efficient intracellular replication in both mammalian and protozoan hosts.

Description: Shiga toxin-encoding Escherichia coli (STEC) regroup strains that carry genes encoding Shiga toxin (Stx). Among intestinal pathogenic E. coli , enterohaemorrhagic E. coli (EHEC) constitute the major subgroup of virulent STEC. EHEC cause serious human disease such as haemorrhagic colitis and haemolytic-uremic syndrome. While EHEC have evolved from enteropathogenic E. coli , hybrids with enteroaggregative E. coli have recently emerged. Of note, some enteroinvasive E. coli also belong to the STEC group. While the LEE (locus of enterocyte effacement) is a key and prominent molecular determinant in the pathogenicity, neither all EHEC nor STEC contain the LEE, suggesting that they possess additional virulence and colonisation factors. Currently, nine protein secretion systems have been described in diderm-lipopolysaccharide bacteria (archetypal Gram-negative) and can be involved in the secretion of extracellular effectors, cell-surface proteins or assembly of cell-surface organelles, such as flagella or pili. In this review, we focus on the secretome of STEC and related enteropathotypes, which are relevant to the colonisation of biotic and abiotic surfaces. Considering the wealth of potential protein trafficking mechanisms, the different combinations of colonisation factors and modulation of their expression is further emphasised with regard to the ecophysiology of STEC.

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 activity of levofloxacin against planktonic and biofilm Stenotrophomonas maltophilia cells and the role played by the multidrug efflux pump SmeDEF were evaluated under conditions relevant to the cystic fibrosis (CF) lung. MIC, MBC and MBEC of levofloxacin were assessed, against five CF strains, under ‘standard’ (CLSI-recommended) and ‘CF-like’ (pH 6.8, 5% CO 2 , in a synthetic CF sputum) conditions. Levofloxacin was tested against biofilms at concentrations (10, 50 and 100 μg mL –1 ) corresponding to achievable serum levels and sputum levels by aerosolisation. smeD expression was evaluated, under both conditions, in planktonic and biofilm cells by RT-PCR. The bactericidal effect of levofloxacin was decreased, in three out of five strains tested, under ‘CF-like’ conditions (MBC: 2–4 vs 8–16 μg mL –1 , under ‘standard’ and ‘CF-like’ conditions, respectively). Biofilm was intrinsically resistant to levofloxacin, regardless of conditions tested (MBECs ≥ 100 μg mL –1 for all strains). Only under ‘CF-like’ conditions, smeD expression increased during planktonic-to-biofilm transition, and in biofilm cells compared to stationary planktonic cells. Our findings confirmed that S. maltophilia biofilm is intrinsically resistant to therapeutic concentrations of levofloxacin. Under conditions relevant to CF, smeD overexpression could contribute to levofloxacin resistance. Further studies are warranted to define the clinical relevance of our findings .

Description: Auranofin is an FDA-approved gold-containing compound used for the treatment of rheumatoid arthritis. Recent reports of antimicrobial activity against protozoa and bacteria indicate that auranofin targets the reductive enzyme thioredoxin reductase (TrxR). We evaluated auranofin as well as five auranofin analogs containing N- heterocyclic carbenes (instead of the triethylphosphane present in auranofin) and five gold-carbene controls for their ability to inhibit or kill Helicobacter pylori in vitro . Auranofin completely inhibited bacterial growth at 1.2 μM. Purified H. pylori TrxR was inhibited by auranofin in a cell-free assay (IC 50 ~88 nM). The most active gold(I)- N- heterocyclic carbene compounds exhibited MICs comparable to auranofin against H. pylori (2 μM), while also exhibiting lower toxicities for human embryonic kidney cells (HEK-293T cells). Median toxic concentrations (TC 50 ) were 13–20-fold higher compared to auranofin indicating that they were less cytotoxic. The N- heterocyclic carbene analogs maybe well tolerated, but further evaluation is needed in vivo . Finally, auranofin was synergistic with the antibiotic amoxicillin, suggesting that targeting both the reductive enzyme TrxR and cell wall synthesis may be effective against H. pylori infections.

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: Legionella feeleii is a Gram-negative pathogenic bacterium that causes Pontiac fever and pneumonia in humans. When L. feeleii serogroup 1 (ATCC 35072) was cultured on BCYE agar plates, two types of colonies were observed and exhibited differences in color, opacity and morphology. Since the two colony types are white rugose and brown translucent, they were termed as white rugose L. feeleii (WRLf) and brown translucent L. feeleii (BTLf), respectively. They exhibited different growth capacities in BYE broth in vitro , and it was found that WRLf could transform to BTLf. Under the electron microscope, it was observed that WRLf secreted materials which could be stained with ruthenium red, which was absent in BTLf. When U937 macrophages and HeLa cells were infected with the bacteria, WRLf manifested stronger internalization ability than BTLf. Intracellular growth in murine macrophages and Acanthamoeba cells was affected by the level of initial phagocytosis. WRLf was more resistant to human serum bactericidal action than BTLf. After being inoculated to guinea pigs, both organisms caused fever in the animals. These results suggest that ruthenium red-stained materials secreted in the surroundings may play a crucial role in determining L. feeleii colony morphology and virulence traits.

Description: Piscirickettsia salmonis is a fastidious intracellular pathogen responsible for high mortality rates in farmed salmonids, with serious economic consequences for the Chilean aquaculture industry. Oxytetracycline and florfenicol are the most frequently used antibiotics against P. salmonis , but routine use could contribute to drug resistance. This study identified differentiated florfenicol susceptibilities in two P. salmonis strains, LF-89 and AUSTRAL-005. The less susceptible isolate, AUSTRAL-005, also showed a high ethidium bromide efflux rate, indicating a higher activity of general efflux pump genes than LF-89. The P. salmonis genome presented resistance nodulation division (RND) family members, a family containing typical multidrug resistance-related efflux pumps in Gram-negative bacteria. Additionally, efflux pump acrAB genes were overexpressed in AUSTRAL-005 following exposure to the tolerated maximal concentration of florfenicol, in contrast to LF-89. These results indicate that tolerated maximum concentrations of florfenicol can modulate RND gene expression and increase efflux pump activity. We propose that the acrAB efflux pump is essential for P. salmonis survival at critical florfenicol concentrations and for the generation of antibiotic-resistant bacterial strains.

Description: Pseudomonas aeruginosa is an opportunistic pathogen, known to develop robust biofilms. Its biofilm development increases when antibiotics are presented at subminimal inhibitory concentrations (MICs) for reasons that remain unclear. In order to identify genes that affect biofilm development under such a sublethal antibiotic stress condition, we screened a transposon (Tn) mutant library of PAO1, a prototype P. aeruginosa strain. Among ~5000 mutants, a fiuA gene mutant was verified to form very defective biofilms in the presence of sub-MIC carbenicillin. The fiuA gene encodes ferrichrome receptor A, involved in the iron acquisition process. Of note, biofilm formation was not decreased in the pchpvd mutant defective in the production of pyochelin and pyoverdine, two well-characterized P. aeruginosa siderophore molecules. Moreover, fiuA , a non-polar fiuA deletion mutant, produced a significantly decreased level of elastase, a major virulence determinant. Mouse airway infection experiments revealed that the mutant expressed significantly less pathogenicity. Our results suggest that the fiuA gene has pleiotropic functions that affect P. aeruginosa biofilm development and virulence. The targeting of FiuA could enable the attenuation of P. aeruginosa virulence and may be suitable for the development of a drug that specifically controls the virulence of this important pathogen.

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: 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: Activating transcription factor 3 (ATF3) is a stress-induced transcriptional regulator in eukaryote. The role of ATF3 in cancer has been well defined, but how ATF3 functions in bacterial infection is not well understood. Pneumococcal infection has been shown to induce ATF3 expression, which subsequently enhances cytokine production and provides protection from lethal Streptococcus pneumoniae infection, but the role of ATF3 in other Gram-positive (G + ) infections remains unclear. Here, we report that infection with other G + bacteria ( Staphylococcus aureus and Listeria monocytogenes ) and with G – bacteria (uropathogenic Escherichia coli ) also significantly induced ATF3 expression. Moreover, the production of cytokines (tumor necrosis factor alpha [TNF]-α, interleukin [IL]-1β, IL-6 and interferon [IFN]-) was enhanced by ATF3 in S. aureus and L. monocytogenes infection, but decreased in uropathogenic E. coli (UPEC) infection. In addition, in S. aureus and L. monocytogenes infections, ATF3 WT mice cleared bacteria more efficiently and had higher survival rates than ATF3 knockout mice. However, in UPEC infection, no significant difference was found in survival rate. Taken together, these data suggest that ATF3 provides protection from S. aureus and L. monocytogenes infections; however, the role of ATF3 in UPEC infection is more complicated and should be further elucidated.

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

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

Description: The presence of carbapenemase gene bla KPC-2 in a wide variety of plasmids, especially conjugative plasmids, is key to the rapid, worldwide spread of carbapenemase enzymes. Thirty-eight, non-duplicated, carbapenem-resistant, clinical Klebsiella pneumoniae isolates were collected, all carrying bla KPC-2 -bearing plasmids. Relaxase analysis was used to classify these plasmids; 8 and 30 plasmids belonged to the MOB P3 and MOB F12 subfamilies, respectively. Phylogenetic analysis revealed two genetic subclades in the MOB F12 subfamily and suggested that these subclades might not have originated from the same ancestor. Crossing PCR, used to sequence fully the type IV secretion system (T4SS, essential structures for conjugative plasmids) of the MOB F12 plasmids, found that T4SSs were distinctively different in certain functional genes, e.g. traS and traG. In conclusion, this study delineated the evolution of bla KPC-2 -bearing plasmids at Huashan Hospital, Shanghai, China. The plasmids bearing bla KPC-2 were diverse and the MOB F12 plasmids were dominant in clinical K. pneumoniae isolates.

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

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

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

Description: Enteroaggregative Escherichia coli (EAEC) is an important diarrhoeal pathogen causing diseases in multiple epidemiological and clinical settings. In developing countries like India, diarrhoeal diseases are one of the major killers among paediatric population and oddly, few studies are available from Indian paediatric population on the variability of EAEC virulence genes. In this study, we examined the distribution of plasmid and chromosomal-encoded virulence determinants in EAEC isolates, and analysed cytokines response generated against EAEC with specific aggregative adherence fimbriae (AAF) type in duodenal biopsies using in vitro organ culture (IVOC) mimicking in vivo conditions. Different virulence marker combinations among strains were reflected as a function of specific adhesins signifying EAEC heterogeneity. fis gene emerged as an important genetic marker apart from aggA and aap . Further, EAEC infection in IVOC showed upregulation of IL-8, IL-1β, IL-6, TNF-α and TLR-5 expression. EAEC with AAFII induced significant TLR-5 and IL-8 response, conceivably owing to more pathogenicity markers. This study sheds light on the pattern of EAEC pathotypes prevalent in North Indian paediatric population and highlights the presence of unique virulence combinations in pathogenic strains. Thus, evident diversity in EAEC virulence and multifaceted bacteria-host crosstalk can provide useful insights for the strategic management of diarrhoeal diseases in India, where diarrhoeal outbreaks are more frequent.

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

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

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

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

Description: Coagulase-negative staphylococci are thought to act as reservoirs of antibiotic resistance genes that can be transferred to Staphylococcus aureus , thus hindering the combat of this bacterium. In this work, we analyzed the presence of plasmids conferring resistance to the antibiotic mupirocin—widely used to treat and prevent S. aureus infections in hospital environments—in nosocomial S. haemolyticus strains. About 12% of the 75 strains tested were resistant to mupirocin, and this phenotype was correlated with the presence of plasmids. These plasmids were shown to be diverse, being either conjugative or mobilizable, and capable of transferring mupirocin resistance to S. aureus . Our findings reinforce that S. haemolyticus , historically and mistakenly considered as a less important pathogen, is a reservoir of resistance genes which can be transferred to other bacteria, such as S. aureus , emphasizing the necessity of more effective strategies to detect and combat this emergent opportunistic pathogen.

Description: Emergence of races in Fusarium oxysporum f. sp. lycopersici ( Fol ) is caused by loss or mutation of at least one avirulence ( AVR ) gene. The product of AVR1 is a small protein (Avr1) secreted by Fol in tomato xylem sap during infection. This protein triggers Fol race 1 specific resistance (I) in tomato, indicating that AVR1 is an AVR gene. Deletion of AVR1 in race 1 resulted in the emergence of race 2, and an additional mutation in AVR2 generated race 3. Previously, we reported a new biotype of race 3, KoChi-1, in which AVR1 was truncated by a transposon Hormin , which suggested a new route to evolution of races in Fol . However, to date no race 2 isolate carrying Hormin -truncated AVR1 has been reported. In this report, we describe such isolates, represented by Chiba-5, in which Hormin insertion occurred in AVR1 at a position different from that in KoChi-1. AVR1 truncation in both isolates resulted in production of defective Avr1 proteins. Chiba-5 and KoChi-1 belong to different phylogenetic clades, A1 and A2, respectively, suggesting that insertion of Hormin in AVR1 in Chiba-5 and KoChi-1 occurred as independent evolutionary events.

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: The LBIT-1200 strain of Bacillus thuringiensis was recently isolated from soil, and showed a 6.4 and 9.5 increase in toxicity, against Manduca sexta and Trichoplusia ni , respectively, compared to HD-73. However, LBIT-1200 was still highly similar to HD-73, including the production of bipyramidal crystals containing only one protein of ~130 000 kDa, its flagellin gene sequence related to the kurstaki serotype, plasmid and RepPCR patterns similar to HD-73, no production of β-exotoxin and no presence of VIP genes. Sequencing of its cry gene showed the presence of a cry1Ac -type gene with four amino acid differences, including two amino acid replacements in domain III, compared to Cry1Ac1, which may explain its higher toxicity. In conclusion, the LBIT-1200 strain is a variant of the HD-73 strain but shows a much higher toxicity, which makes this new strain an important candidate to be developed as a bioinsecticide, once it passes other tests, throughout its biotechnological development.

Description: Gastrointestinal (GI) leakage in Clostridium difficile -associated diarrhea (CDAD) is well known but is not routinely assessed in clinical practice. Serum (1-〉3)-β-D-glucan (BG), a fungal cell wall component used as a biomarker for invasive fungal disease, was tested in a CDAD mouse model with and without probiotics. Higher serum fluorescein isothiocyanate-dextran (FITC-dextran) and spontaneous gram-negative bacteremia, GI leakage indicators, were frequently found in CDAD mice, which died compared with those which survived. BG, serum macrophage inflammatory protein-2 and FITC-dextran but not quantitative blood bacterial count differentiated the clinical severity. Interestingly, a specific dose of Lactobacillus rhamnosus L34 attenuated CDAD and decreased serum BG and FITC-dextran, but not other parameters. BG also showed a higher area under the receiver operating characteristic curve for 7-day mortality than FITC-dextran. Fifty-five percent of CDAD mice with BG ≥ 60 pg/ml (the human negative cut-off value for invasive fungal disease) at 1 day after C. difficile gavage died within 7 days. In conclusion, s erum BG was elevated in mice with severe CDAD, an established model of GI leakage with a strong association with mortality rate. BG monitoring in patients with CDAD is of interest as both a potential prognostic tool and a therapeutic efficacy indicator.

Description: Chlamydial species are common intracellular parasites that cause various diseases, mainly characterized by persistent infection, which lead to inflammatory responses modulated by pattern recognition receptors (PRRs). The best understood PRRs are the extracellular Toll-like receptors, but recent significant advances have focused on two important proteins, NOD1 and NOD2, which are members of the intracellular nucleotide-binding oligomerization domain receptor family and are capable of triggering the host innate immune signaling pathways. This results in the production of pro-inflammatory cytokines, which is vital for an adequate host defense against intracellular chlamydial infection. NOD1/2 ligands are known to derive from peptidoglycan, and the latest research has resolved the paradox of whether chlamydial species possess this bacterial cell wall component; this finding is likely to promote in-depth investigations into the interaction between the NOD proteins and chlamydial pathogens. In this review, we summarize the basic characteristics and signal transduction functions of NOD1 and NOD2 and highlight the new research on the roles of NOD1 and NOD2 in the host defense against chlamydial infection.

Description: Two strains of Aeromonas salmonicida , YK and BG, were isolated from largemouth bronze gudgeon and northern whitefish in China, and identified as A. salmonicida subsp. salmonicida based on phylogenetic analysis of vapA and 16S rRNA gene sequences. YK and BG originated from freshwater fish, one of which belonged to the cyprinid family, and the strains showed a difference in virulence. Subsequently, we performed whole genome sequencing of the strains, and comparison of their genomic sequences to the genome of the A449 reference strain revealed various genomic rearrangements, including a new variant of the genomic island AsaGEI in BG, designated as AsaGEI2c . This is the first report on a GEI of A. salmonicida strain from China. Furthermore, both YK and BG strains contained a Tn7 transposon inserted at the same position in the chromosome. Finally, IS-dependent rearrangements on pAsa5 are deemed likely to have occurred, with omission of the resD gene in both strains as well as omission of genes related to the IncF conjugal transfer system in the YK isolate. This study demonstrates that A. salmonicida subsp. salmonicida can infect non-salmonids (cyprinids) in addition to salmonids, and that AsaGEI2c might be useful as a geographical indicator of Chinese A. salmonicida subsp. salmonicida isolates.

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: Burkholderia pseudomallei is an intracellular Gram-negative bacterial pathogen and the causative agent of melioidosis, a widespread disease in Southeast Asia. Reactive nitrogen, in an intermediate form of nitric oxide (NO), is one of the first lines of defense used by host cells to eliminate intracellular pathogens, through the stimulation of inducible nitric oxide synthase (iNOS). Studies in phagocytotic cells have shown that the iNOS response is muted in B. pseudomallei infection, and implicated the rpoS sigma factor as a key regulatory factor mediating suppression. The liver is a main visceral organ affected by B. pseudomallei , and there is little knowledge about the interaction of liver cells and B. pseudomallei . This study investigated the induction of iNOS, as well as autophagic flux and light-chain 3 (LC3) localization in human liver (HC04) cells in response to infection with B. pseudomallei and its rpoS deficient mutant. Results showed that the rpoS mutant was unable to suppress iNOS induction and that the mutant showed less induction of autophagy and lower co-localization with LC3, and this was coupled with a lower intracellular growth rate. Combining these results suggest that B. pseudomallei rpoS is an important factor in establishing infection in liver cells.

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: Bacteriophages are viruses that infect bacteria. There are an estimated 10 31 phage on the planet, making them the most abundant form of life. We are rapidly approaching the centenary of their identification, and yet still have only a limited understanding of their role in the ecology and evolution of bacterial populations. Temperate prophage carriage is often associated with increased bacterial virulence. The rise in use of technologies, such as genome sequencing and transcriptomics, has highlighted more subtle ways in which prophages contribute to pathogenicity. This review discusses the current knowledge of the multifaceted effects that phage can exert on their hosts and how this may contribute to bacterial adaptation during infection.

Description: Phthiocerol dimycocerosates (PDIMs) and phenolic glycolipids (PGLs) contribute to the pathogenicity of several mycobacteria. Biosynthesis of these virulence factors requires polyketide synthases and other enzymes that represent potential targets for the development of adjuvant antivirulence drugs. We used six isogenic Mycobacterium marinum mutants, each with a different gene knockout in the PDIM/PGL biosynthetic pathway, to probe the pleiotropy of mutations leading to PDIM – PGL – , PDIM + PGL – or PDIM – PGL + phenotypes. We evaluated the M. marinum mutants for changes in antibiotic susceptibility, cell envelope permeability, biofilm formation, surface properties, sliding motility and virulence in an amoeba model. The analysis also permitted us to begin exploring the hypothesis that different gene knockouts rendering the same PDIM and/or PGL deficiency phenotypes lead to M. marinum mutants with equivalent pleiotropic profiles. Overall, the results of our study revealed a complex picture of pleiotropic patterns emerging from different gene knockouts, uncovered unexpected phenotypic inequalities between mutants, and provided new insight into the phenotypic consequences of gene knockouts in the PDIM/PGL biosynthetic pathway.

Description: Culture medium from an isolate of the fungus Aspergillus candidus was extracted, fractionated and examined to discover compounds antagonistic to plant-parasitic nematodes that are important pathogens of agricultural crops. Column, thin layer and preparative chromatographies and spectral and elemental analyses, were used to isolate and identify two major constituents of an active fraction (Fraction F) obtained from the medium. Compound 1 was identified as 2-hydroxypropane-1, 2, 3-tricarboxylic acid (citric acid). Compound 2 was identified as 3-hydroxy-5-methoxy-3-(methoxycarbonyl)-5-oxopentanoic acid, an isomer of 1, 2-dimethyl citrate. Compound 1 and a citric acid standard, each tested at 50 mg mL –1 in water, decreased hatch from eggs of the plant-parasitic nematode Meloidogyne incognita by more than 94%, and completely immobilized second-stage juveniles after 4–6 days exposure. Fraction F and Compounds 1 and 2 decreased the mobility of adults of the plant-parasitic nematode Ditylenchus destructor in vitro . Fraction F (25 mg mL –1 ) inhibited mobility 〉99% at 72 hrs. Compounds 1 and 2 (50 mg mL –1 ) each inhibited mobility more than 25% at 24 hr and more than 50% at 72 hr. This is the first assignment of nematode-antagonistic properties to specifically identified A. candidus metabolites.

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

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

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

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

Description: An RNAi system based on T7 RNA polymerase (TRNAP) was designed and examined in Aspergillus fumigatus . This system consists of two elements; an inducible T7RNAP expressing cassette and an AMA1-based episomal RNAi plasmid. These constructs were transformed into the A. fumigatus protoplasts and the efficiency of this system was tested in downregulation of alb1 gene. Upon the induction of T7RNAP expression, the recombinant T7RNAP was able to recognize T7 promoters, which were located on the episomal plasmid and in opposite direction. As a result, the bidirectional transcription of alb1 fragment led to the silencing of the target gene. However, our results demonstrated that this silencing system is unstable and may not be applicable in preparation of RNAi libraries.

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: A simple diagnosis of the presence or absence of an infection is an uninformative metric when individuals differ considerably in their tolerance to different infection loads or resistance to rates of disease progression. Models that incorporate the relationship between the progression of the infection with the potential alternate outcomes provide a far more powerful predictive tool than diagnosis alone. The global decline of amphibians has been amplified by Batrachochytrium dendrobatidis , a pathogen that can cause the fatal disease chytridiomycosis. We measured the infection load and observed signs of disease in Litoria aurea . Receiver operating characteristic curves were used to quantify the dissimilarity between the infection loads of L. aurea that showed signs associated with chytridiomycosis and those that did not. Litoria aurea had a 78% probability of developing chytridiomycosis past a threshold of 68 zoospore equivalents (ZE) per swab and chytridiomycosis occurred within a variable range of 0.5–490 ZE. Studies should incorporate a species-specific threshold as a predictor of chytridiomycosis, rather than a binary diagnosis. Measures of susceptibility to chytridiomycosis must account not only for the ability of B. dendrobatidis to increase its abundance on the skin of amphibians but also to determine how each species tolerates these infection loads.

Description: The control and prevention of biofilm-related infections is an important public healthcare issue. Given the increasing antibiotic resistance among bacteria and fungi that cause serious infections in humans, promotion of new strategies combating microorganisms has been essential. One attractive approach to inactivate microorganisms is the use of semiconductor photo-catalysis, which has become the subject of extensive research. In this study, the bactericidal properties of four photo-catalysts, TiO 2 , TiO 2 -S, TiO 2 -Eu and TiO 2 -Eu-S, were investigated against established 24, 48, 72 and 96 h biofilms of Enterococcus . The exposure of biofilms to the catalysts induced the production of superoxide radical anions. The best photo-catalytic inactivation was achieved with the TiO 2 -Eu-S and TiO 2 -S nanopowders and 24 h biofilms. Transmission electron microscopy images showed significant changes in the structure of the biofilm cells following photo-inactivation. The results suggest that doping with europium and modifying the surface with sulphate groups enhanced the bactericidal activity of the TiO 2 nanoparticles against enterococcal biofilms.

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: The Pseudomonas aeruginosa Chp chemosensory system regulates twitching motility, intracellular adenosine 3 '' 5 ' -cyclic monophosphate (cAMP) levels and is postulated to be involved in directional twitching towards phosphatidylethanolamine (PE). Because PilJ is the only methyl-accepting chemotaxis protein (MCP) identified in the Chp system, we determined the role of PilJ in mediating signal transduction for the distinct outputs of this system. Mutants that lack the periplasmic domain of PilJ ( pilJ 74-273 ) showed lower levels of cAMP but retained directional twitching towards PE. While initial studies revealed reduced twitching motility by PilJ 74-273 , this was due to decreased cAMP levels. Our data illustrate the importance of the periplasmic domain of PilJ in regulating cAMP. This is the first time a defined domain within PilJ has been identified as having a distinct role in signal transduction.

Description: Histone-like nucleoid-structuring (H-NS) proteins, which are conserved in Gram-negative bacteria, bind DNA and act as the global transcriptional repressors. In this study, we identified and characterized the xrvC gene encoding a H-NS protein in Xathomonas oryzae pv. oryzae ( Xoo ) Philippines strain PXO99 A . Compared with the wild type, the xrvC -deficient mutant of PXO99 A (named PXO99 xrvC ) showed a reduced growth rate in both nutrient-rich and nutrient-limited media. Interestingly, PXO99 xrvC exhibited significantly reduced virulence on rice cultivar IRBB214, but its virulence on 31 other rice cultivars was not affected. Transcriptional analysis revealed that the expression of hrpG , hrpX and hpa1 and of 15 out of 18 tested non-TAL (transcription activator-like) effector genes was decreased significantly in the xrvC mutant compared with that in the wild type. In addition, loss of xrvC also impaired the induction of the rice susceptibility gene Os8N3 in IRBB214 by PXO99 A . Our results suggest that the xrvC gene is involved in bacterial growth, and it plays a vital role in virulence by positively regulating the expression of hrp genes and non-TAL effector genes in PXO99 A and the susceptibility gene Os8N3 in rice.

Description: In Xanthomonas oryzae pv. oryzae , the pathogen of bacterial leaf blight of rice, hrp gene expression is regulated by the key hrp regulators HrpG and HrpX. HrpG regulates hrpX and hrpA , and HrpX regulates the other hrp genes on hrpB–hrpF operons. We previously examined the expression of the HrpX-regulated hrp gene hrcU and demonstrated that hrp gene expression is highly induced in a certain nutrient-poor medium containing xylose. In the present study, we found that the induction level of HrpX-regulated hrp genes was higher in medium with xylose than in media with any other sugar sources (glucose, sucrose and fructose), but that expression of hrpG , hrpX and hrpA was independent of the sugar sources. In western blot analysis, the accumulation of HrpX was reduced in media with a sugar other than xylose, probably as a result of proteolysis, but the addition of xylose canceled this reduced accumulation of the protein. The results suggest that proteolysis of HrpX is an important hrp regulatory mechanism and that xylose specifically suppresses this proteolysis, resulting in active hrp gene expression in X. oryzae pv. oryzae .

Description: Antibiotic therapy has long term consequences in the intestinal microbiome. Clostridium difficile has a well-known role in antibiotic-associated diarrhea, but in addition, persistent infection with this organism may increase the risk for developing inflammatory bowel disease. Here, recent literature on how the intestinal microbiome is altered by antibiotic therapy is presented.

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: The aim of the present study was to verify whether penicillin-resistant, ampicillin-susceptible Enterococcus faecalis (PRASEF) occurred in Brazil prior to the beginning of the 21st century, and to verify whether ampicillin susceptibility can predict susceptibility to other β-lactams in E. faecalis with this inconsistent phenotype. The presence of polymorphisms in the pbp4 gene and genetic diversity among the isolates were investigated. Of 21 PRASEF analyzed, 5 (23.8%) and 4 (19.0%) were imipenem and piperacillin resistant simultaneously by disk diffusion and broth dilution respectively, contradicting the current internationally accepted standards of susceptibility testing. Sequencing of pbp4 gene revealed an amino acid substitution (Asp-573-〉Glu) in all PRASEF isolates but not in the penicillin-susceptible, ampicillin-susceptible E. faecalis . Most PRASEF (90.5%) had related pulsed-field gel electrophoresis profiles, but were different from other PRASEF described to date. Results demonstrate that penicillin-resistant, ampicillin-susceptible phenotype was already a reality in the 1990s in E. faecalis isolates in different Brazilian states, and some of these isolates were also imipenem- and piperacillin-resistant; therefore, internationally accepted susceptibility criteria cannot be applied to these isolates. According to pbp4 gene sequencing, this study suggests that a specific amino acid substitution in pbp4 gene found in all PRASEF analyzed is associated with penicillin resistance.

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

Description: The Mycobacterium tuberculosis Rv0679c protein is a surface protein that contributes to host cell invasion. We previously showed that a single nucleotide transition of the Rv0679c gene leads to a single amino acid substitution from asparagine to lysine at codon 142 in the Beijing genotype family. In this study, we examined the immunological effect of this substitution. Several recombinant proteins were expressed in Escherichia coli and Mycobacterium smegmatis and characterized with antisera and two monoclonal antibodies named 5D4-C2 and 8G10-H2. A significant reduction of antibody binding was detected by enzyme-linked immunosorbent assay (ELISA) and western blot analysis in the Lys142-type protein. This reduction of 8G10-H2 binding was more significant, with the disappearance of a signal in the proteins expressed by recombinant mycobacteria in western blot analysis. In addition, epitope mapping analysis of the recombinant proteins showed a linear epitope by 5D4-C2 and a discontinuous epitope by 8G10-H2. The antibody recognizing the conformational epitope detected only mycobacterial Asn142-type recombinant protein. Our results suggest that a single amino acid substitution of Rv0679c has potency for antigenic change in Beijing genotype strains.

Description: Listeria monocytogenes possesses the highest number of leucine-rich repeat (LRR)-containing proteins among all Gram-positive bacteria; these LRR-containing molecules are known as the ‘internalin’ family. To understand the functions of largely uncharacterized LRR-containing molecules, we constructed seven deletion mutants in the L. monocytogenes H7858 strain targeting genes in this family and tested their virulence. Among the seven mutants, the LMOh7858_0369 strain and the LMOh7858_2546 strain showed significantly impaired invasiveness of HepG2 cells. We further tested the virulence of these two strains in the intravascular sepsis model using BALB/c mice. Interestingly, the LMOh7858_0369 strain showed significant reduction in organ colonization, bacteremia and invasion of the brain compared with the parental wild-type strain. Host immune responses to listerial intravascular infection were measured at 24 and 72 h post-infection. Transcript levels of several proinflammatory cytokines and chemokines were significantly lower when induced by the lmOh7858_0369 strain than when induced by the wild type. These results suggest that the putative LRR-containing protein encoded by LMOh7858_0369 might be a novel virulence factor of the L. monocytogenes H7858 strain.

Description: Several Gram-positive pathogens scavenge host-derived heme to satisfy their nutritional iron requirement. However, heme is a toxic molecule capable of damaging the bacterial cell. Gram-positive pathogens within the phylum Firmicutes overcome heme toxicity by sensing heme through HssRS, a two-component system that regulates the heme detoxification transporter HrtAB. Here we show that heme sensing by HssRS and heme detoxification by HrtAB occur in the insect pathogen Bacillus thuringiensis . We find that in B. thuringiensis , HssRS directly regulates an operon, hrmXY , encoding hypothetical membrane proteins that are not found in other Firmicutes with characterized HssRS and HrtAB systems. This novel HssRS-regulated operon or its orthologs BMB171_c3178 and BMB171_c3330 are required for maximal heme resistance. Furthermore, the activity of HrmXY is not dependent on expression of HrtAB. These results suggest that B. thuringiensis senses heme through HssRS and induces expression of separate membrane-localized systems capable of overcoming different aspects of heme toxicity.

Description: The aim of this study was to understand the pathogenesis of motile aeromonas septicemia caused by an emergent, high virulent Aeromonas hydrophila (vAh) in channel catfish, Ictalurus punctatus . Adipose fin clipped catfish were challenged with vAh using a waterborne challenge method, and the distribution of vAh over a time course was detected and quantified using real-time polymerase chain reaction. The results showed that 77.8% of fish died within 48 h post challenge with mean day to death of 1.5 days. At 2 h post challenge, vAh (inferred from genomic DNA copies or genome equivalents) was detected in all external and internal tissues sampled. Gill had the highest vAh cells at 1 h post challenge. Spleen harbored the most vAh cells among internal organs at 4 h post challenge. The tissues/organs with most vAh cells detected at 8 h post challenge were adipose fin, blood, intestine, kidney and skin, while liver showed the highest vAh cells at 24 h post challenge. These results suggest that vAh was able to rapidly proliferate and spread, following wound infection, through the fish blood circulation system and cause mortality within 8–24 h.