Numerical modelling of post-seismic rupture propagation after the Sumatra 26.12.2004 earthquake constrained by GRACE gravity data

International audienceIn the last decades, the development of the surface and satellite geodetic and geophysical observations brought a new insights into the seismic cycle, documenting new features of inter-, co-, and post-seismic processes. In particular since 2002 satellite mission GRACE provides monthly models of the global gravity field with unprecedented accuracy showing temporal variations of the Earth's gravity field, including those caused by mass redistribution associated with earthquake processes. When combined with GPS measurements, these new data have allowed to assess the relative importance of afterslip and viscoelastic relaxation after the Sumatra 26.12.2004 earthquake. Indeed the observed post-seismic crustal displacements were fitted well by a viscoelastic relaxation model assuming Burgers body rheology for the asthenosphere (60–220 km deep) with a transient viscosity as low as 4 × 10^17 Pas and constant ~ 10^19 Pas steady state viscosity in the 60–660-km depth range. However, even the low-viscosity asthenosphere provides the amplitude of strain which gravity effect does not exceed 50 per cent of the GRACE gravity variations, thus additional localized slip of about 1 m was suggested at downdip extension of the coseismic rupture. Post-seismic slip at coseismic rupture or its downdip extension has been suggested by several authors but the mechanism of the post-seismic fault propagation has never been investigated numerically. Depth and size of localized slip area as well as rate and time decay during the post-seismic stage were either assigned a priory or estimated by fitting real geodesy or gravity data. In this paper we investigate post-seismic rupture propagation by modelling two consequent stages. First, we run a long-term, geodynamic simulation to self-consistently produce the initial stress and temperature distribution. At the second stage, we simulate a seismic cycle using results of the first step as initial conditions. The second short-term simulation involves three substeps, including additional stress accumulation after part of the subduction channel was locked; spontaneous coseismic slip; formation and development of damage zones producing afterslip. During the last substep post-seismic stress leads to gradual ~1 m slip localized at three faults around ~100-km downdip extension of the coseismic rupture. We used the displacement field caused by the slip to calculate pressure and density variations and to simulate gravity field variations. Wavelength of calculated gravity anomaly fits well to that of the real data and its amplitude provides about 60 per cent of the observed GRACE anomaly. Importantly, the surface displacements caused by the estimated afterslip are much smaller than those registered by GPS networks. As a result cumulative effect of Burgers rheology viscoelastic relaxation (which explains measured GPS displacements and about a half of gravity variations) plus post-seismic slip predicted by damage rheology model (which causes much smaller surface displacements but provides another half of the GRACE gravity variations) fits well to both sets of the real data. Hence, the presented numerical modelling based on damage rheology supports the process of post-seismic downdip rupture propagation previously hypothesized from the GRACE gravity data

Decay of heavy-light hybrids in HQET sum rules

The decay widths of the 0^{++} and 1^{-+} heavy-light hybrids to B(D) and pion are calculated by using the QCD sum rules. The interpolated current of the hybrid is chosen as $g\bar q\gamma_{\alpha}G_{\alpha\mu}^aT^ah_{\it v}(x)$. In order to simplify the calculation and avoid the ambiguity of three-point correlation function, a two-point correlation function between the pion and vacuum is used instead. The decay width of the $0^{++}\to B(D)$ is about 12(16) MeV while the $1^{-+}\to B(D)$ is around 0.4(1.8) MeV. We keep the leading order of 1/M_Q expansion in our calculation for convenience.Comment: 14 pages, latex file, 4 ps figs, Published version, some numerical results change

Decays of the ρ^(1−+)\hat\rho(1^{-+}) Exotic Hybrid and η\eta-η′\eta' Mixing

QCD sum-rules are used to calculate the $\hat\rho(1^{-+})\to\pi\eta, \pi\eta'$ decay widths of the exotic hybrid in two different $\eta-\eta'$ mixing schemes. In the conventional flavour octet-singlet mixing scheme, the decay widths are both found to be small, while in the recently-proposed quark mixing scheme, the decay width $\Gamma_{\hat\rho\to\eta\pi}\approx 250 MeV$ is large compared with the decay width $\Gamma_{\hat\rho\to\eta^\prime\pi}\approx 20 MeV$. These results provide some insight into $\eta$-$\eta'$ mixing and hybrid decay features.Comment: latex2e, 11 pages with 4 embedded eps figures. v 2 corrects reference [5] and minor error in equation (11

A Study of the \eta \pi^{0} Spectrum and Search for a J^{PC} = 1^{-+} Exotic Meson

A partial wave analysis (PWA) of the of the $\eta \pi ^0$ system (where $\eta \to \gamma \gamma$) produced in the charge exchange reaction $\pi ^-p\to \eta \pi ^0n$ at an incident momentum of 18 GeV$/c$ is presented as a function of ${\eta \pi ^0}$ invariant mass, $m_{\eta\pi^0}$, and momentum transfer squared, $t_{\pi^{-}\to\eta\pi}$, from the incident $\pi^-$ to the outgoing ${\eta\pi ^0}$ system. $S$, $P$ and $D$ waves were included in the PWA. The $a_0(980)$ and $a_2(1320)$ states are clearly observed in the overall ${\eta\pi ^0}$ effective mass distribution as well as in the amplitudes associated with $S$ wave and $D$ waves respectively after partial wave decomposition. The observed distributions in moments (averages of spherical harmonics) were compared to the results from the PWA and the two are consistent. The distribution in $t_{\pi^{-}\to\eta\pi}$ for individual $D$ waves associated with natural and unnatural parity exchange in the $t$-channel are consistent with Regge phenomenology. Of particular interest in this study is the $P$ wave since this leads to an exotic $J^{PC}=1^{-+}$ for the $\eta \pi$ system. A $P$ wave is present in the data, however attempts to describe the mass dependence of the amplitude and phase motion with respect to the $D$ wave as a Breit-Wigner resonance are problematic. This has implications regarding the existence of a reported exotic $J^{PC} = 1^{-+}$ meson decaying into $\eta \pi^0$ with a mass near 1.4 GeV$/c^2$.Comment: 19 pages, 29 figures, to appear in Phys. Rev.

Can the Mechanism for π1→ηπ,η′π\pi_1\to \eta\pi,\eta'\pi Hybrid Decays be Detected?

Two mechanisms for the $\pi_1$ ($J^{PC}=1^{-+}$) hybrid meson decay processes $\pi_1\to\eta\pi,\eta'\pi$ are investigated. These mechanisms are applied to $\phi\to\eta\gamma,\eta'\gamma$ and $J/\psi\to\eta\gamma,\eta'\gamma$ decays to illustrate the validity of the decay mechanisms and to obtain independent information on the coupling of $\eta,\eta'$ to quark and gluonic operators. From this information, we find that $\Gamma(\pi_1\to\eta\pi)/\Gamma(\pi_1\to\eta'\pi)$ is substantially different in the two decay mechanisms, and hence future experimental measurements of this ratio will provide valuable information for substantiating the hybrid nature of these states and for determining the mechanism for these hybrid decays.Comment: 5 pages, revtex, 1 eps figure embedded in manuscript. Analysis and references extended in v

Evidence for Exotic J^{PC}=1^{-+} Meson Production in the Reaction pi- p --> eta pi- p at 18 GeV/c

Details of the analysis of the eta pi- system studied in the reaction pi^{-} p --> eta pi^{-} p at 18 GeV/c are given. Separate analyses for the 2 gamma and pi+ pi- pi0 decay modes of the eta are presented. An amplitude analysis of the data indicates the presence of interference between the a(2)(1320)- and a J^{PC}=1^{-+} wave between 1.2 and 1.6 GeV/c^2. The phase difference between these waves shows phase motion not attributable solely to the a(2)(1320)-. The data can be fitted by interference between the a(2)(1320)- and an exotic 1^{-+} resonance with M = 1370 +-16 +50 -30} MeV/c^2 and Gamma = 385 +- 40 +65 -105 MeV/c^2. Our results are compared with those of other experiments.Comment: 50 pages of text and 34 figure

The masses and decay widths of heavy hybrid mesons

We first derive the mass sum rules for the heavy hybrid mesons to obtain the binding energy and decay constants in the leading order of HQET. The pionic couplings between the lightest $1^{-+}$ hybrid $(Q\bar q g)$ and the lowest three heavy meson doublets are calculated with the light cone QCD sum rules. With $SU_f (3)$ flavor symmetry we calculate the widths for all the possible two-body decay processes with a Goldstone boson in the final state. The total width of the $1^{-+}$ hybrid is estimated to be 300 MeV. We find the dominant decay mode of the $1^{-+}$ hybrid is $1^{-+}\to \pi + 1^+$ where the $1^+$ heavy meson belongs to the $(1^+,2^+)$ doublet. Its branching ratio is about 80% so this mode can be used for the experimental search of the lowest heavy hybrid meson.Comment: 20 pages + 12 PS figures, introduction revised, Fig 7 updated, to appear in Phys. Rev.

Interpretation of Experimental J^PC Exotic Signals

We investigate theoretical interpretations of the 1.4 GeV J^PC exotic resonance reported by the E852 collaboration. It is argued that interpretation in terms of a hybrid meson is untenable. A K-matrix analysis shows that the 1.4 GeV enhancement in the E852 eta pi data can be understood as an interference of a non-resonant Deck-type background and a resonance at 1.6 GeV. A final state rescattering calculation shows that the 1.6 GeV hybrid has a eta pi width which is bounded above by 57 \pm 14 MeV.Comment: 23 pages, LaTeX, 4 encapsulated postscript figures. Accepted for publication by Physical Review

Regulation of the Fruit-Specific PEP Carboxylase SlPPC2 Promoter at Early Stages of Tomato Fruit Development

The SlPPC2 phosphoenolpyruvate carboxylase (PEPC; EC 4.1.1.31) gene from tomato (Solanum lycopersicum) is differentially and specifically expressed in expanding tissues of developing tomato fruit. We recently showed that a 1966 bp DNA fragment located upstream of the ATG codon of the SlPPC2 gene (GenBank AJ313434) confers appropriate fruit-specificity in transgenic tomato. In this study, we further investigated the regulation of the SlPPC2 promoter gene by analysing the SlPPC2 cis-regulating region fused to either the firefly luciferase (LUC) or the β-glucuronidase (GUS) reporter gene, using stable genetic transformation and biolistic transient expression assays in the fruit. Biolistic analyses of 5′ SlPPC2 promoter deletions fused to LUC in fruits at the 8th day after anthesis revealed that positive regulatory regions are mostly located in the distal region of the promoter. In addition, a 5′ UTR leader intron present in the 1966 bp fragment contributes to the proper temporal regulation of LUC activity during fruit development. Interestingly, the SlPPC2 promoter responds to hormones (ethylene) and metabolites (sugars) regulating fruit growth and metabolism. When tested by transient expression assays, the chimeric promoter:LUC fusion constructs allowed gene expression in both fruit and leaf, suggesting that integration into the chromatin is required for fruit-specificity. These results clearly demonstrate that SlPPC2 gene is under tight transcriptional regulation in the developing fruit and that its promoter can be employed to drive transgene expression specifically during the cell expansion stage of tomato fruit. Taken together, the SlPPC2 promoter offers great potential as a candidate for driving transgene expression specifically in developing tomato fruit from various tomato cultivars