Earthquakes of the Nepal Himalaya : towards a physical model of the seismic cycle

Home to hundreds of millions of souls and land of excessiveness, the Himalaya is also the locus of a unique seismicity whose scope and peculiarities still remain to this day somewhat mysterious. Having claimed the lives of kings, or turned ancient timeworn cities into heaps of rubbles and ruins, earthquakes eerily inhabit Nepalese folk tales with the fatalistic message that nothing lasts forever. From a scientific point of view as much as from a human perspective, solving the mysteries of Himalayan seismicity thus represents a challenge of prime importance. Documenting geodetic strain across the Nepal Himalaya with various GPS and leveling data, we show that unlike other subduction zones that exhibit a heterogeneous and patchy coupling pattern along strike, the last hundred kilometers of the Main Himalayan Thrust fault, or MHT, appear to be uniformly locked, devoid of any of the “creeping barriers” that traditionally ward off the propagation of large events. The approximately 20 mm/yr of reckoned convergence across the Himalaya matching previously established estimates of the secular deformation at the front of the arc, the slip accumulated at depth has to somehow elastically propagate all the way to the surface at some point. And yet, neither large events from the past nor currently recorded microseismicity nearly compensate for the massive moment deficit that quietly builds up under the giant mountains. Along with this large unbalanced moment deficit, the uncommonly homogeneous coupling pattern on the MHT raises the question of whether or not the locked portion of the MHT can rupture all at once in a giant earthquake. Univocally answering this question appears contingent on the still elusive estimate of the magnitude of the largest possible earthquake in the Himalaya, and requires tight constraints on local fault properties. What makes the Himalaya enigmatic also makes it the potential source of an incredible wealth of information, and we exploit some of the oddities of Himalayan seismicity in an effort to improve the understanding of earthquake physics and cipher out the properties of the MHT. Thanks to the Himalaya, the Indo-Gangetic plain is deluged each year under a tremendous amount of water during the annual summer monsoon that collects and bears down on the Indian plate enough to pull it away from the Eurasian plate slightly, temporarily relieving a small portion of the stress mounting on the MHT. As the rainwater evaporates in the dry winter season, the plate rebounds and tension is increased back on the fault. Interestingly, the mild waggle of stress induced by the monsoon rains is about the same size as that from solid-Earth tides which gently tug at the planets solid layers, but whereas changes in earthquake frequency correspond with the annually occurring monsoon, there is no such correlation with Earth tides, which oscillate back-and-forth twice a day. We therefore investigate the general response of the creeping and seismogenic parts of MHT to periodic stresses in order to link these observations to physical parameters. First, the response of the creeping part of the MHT is analyzed with a simple spring-and-slider system bearing rate-strengthening rheology, and we show that at the transition with the locked zone, where the friction becomes near velocity neutral, the response of the slip rate may be amplified at some periods, which values are analytically related to the physical parameters of the problem. Such predictions therefore hold the potential of constraining fault properties on the MHT, but still await observational counterparts to be applied, as nothing indicates that the variations of seismicity rate on the locked part of the MHT are the direct expressions of variations of the slip rate on its creeping part, and no variations of the slip rate have been singled out from the GPS measurements to this day. When shifting to the locked seismogenic part of the MHT, spring-and-slider models with rate-weakening rheology are insufficient to explain the contrasted responses of the seismicity to the periodic loads that tides and monsoon both place on the MHT. Instead, we resort to numerical simulations using the Boundary Integral CYCLes of Earthquakes algorithm and examine the response of a 2D finite fault embedded with a rate-weakening patch to harmonic stress perturbations of various periods. We show that such simulations are able to reproduce results consistent with a gradual amplification of sensitivity as the perturbing period get larger, up to a critical period corresponding to the characteristic time of evolution of the seismicity in response to a step-like perturbation of stress. This increase of sensitivity was not reproduced by simple 1D-spring-slider systems, probably because of the complexity of the nucleation process, reproduced only by 2D-fault models. When the nucleation zone is close to its critical unstable size, its growth becomes highly sensitive to any external perturbations and the timings of produced events may therefore find themselves highly affected. A fully analytical framework has yet to be developed and further work is needed to fully describe the behavior of the fault in terms of physical parameters, which will likely provide the keys to deduce constitutive properties of the MHT from seismological observations.

A tale of two cobalts

There are important problems to overcome if solar energy or other renewable energy sources are to be used effectively on a global scale. Solar photons must not only be harvested and converted into a usable form, but they must also be efficiently stored so that energy is available for use on cloudy days and at night. In this work, both the energy conversion and energy storage problems are addressed. Specifically, two cobalt complexes were designed and their reactivity probed for applications in energy conversion and storage. The first chapter describes a cobalt complex that is the first example of a dimeric cobalt compound with two singly proton-bridged cobaloxime units linked by a central BO4--bridge. Using electrochemical methods, the redox properties of the dimer were evaluated and it was found to be an electrocatalyst for proton reduction in acetonitrile.

Because hydrogen gas is difficult to handle and store, the hydrogenation of CO2 and later dehydrogenation of the liquid product, formic acid, has been proposed as a hydrogen storage system. Thus, a second complex, described in chapter two, supported by a triphosphine ligand framework was used as a catalyst precursor for this key dehydrogenation step. The studies here demonstrate the efficacy of the complex as a precatalyst for the desired reaction, with good conversion of starting formic acid to CO2 and H2. In order to better understand the properties of the triphosphine cobalt complex, a synthetic procedure for substituting electron donating groups (e.g., methoxy groups) onto the ligand was investigated, yielding a novel diphosphine cobalt(II) complex.

DNA-mediated oxidation of transcription factor p53

Transcription factor p53 is the most commonly altered gene in human cancer. As a redox-active protein in direct contact with DNA, p53 can directly sense oxidative stress through DNA-mediated charge transport. Electron hole transport occurs with a shallow distance dependence over long distances through the π-stacked DNA bases, leading to the oxidation and dissociation of DNA-bound p53. The extent of p53 dissociation depends upon the redox potential of the response element DNA in direct contact with each p53 monomer. The DNA sequence dependence of p53 oxidative dissociation was examined by electrophoretic mobility shift assays using radiolabeled oligonucleotides containing both synthetic and human p53 response elements with an appended anthraquinone photooxidant. Greater p53 dissociation is observed from DNA sequences containing low redox potential purine regions, particularly guanine triplets, within the p53 response element. Using denaturing polyacrylamide gel electrophoresis of irradiated anthraquinone-modified DNA, the DNA damage sites, which correspond to locations of preferred electron hole localization, were determined. The resulting DNA damage preferentially localizes to guanine doublets and triplets within the response element. Oxidative DNA damage is inhibited in the presence of p53, however, only at DNA sites within the response element, and therefore in direct contact with p53. From these data, predictions about the sensitivity of human p53-binding sites to oxidative stress, as well as possible biological implications, have been made. On the basis of our data, the guanine pattern within the purine region of each p53-binding site determines the response of p53 to DNA-mediated oxidation, yielding for some sequences the oxidative dissociation of p53 from a distance and thereby providing another potential role for DNA charge transport chemistry within the cell.

To determine whether the change in p53 response element occupancy observed in vitro also correlates in cellulo, chromatin immunoprecipition (ChIP) and quantitative PCR (qPCR) were used to directly quantify p53 binding to certain response elements in HCT116N cells. The HCT116N cells containing a wild type p53 were treated with the photooxidant [Rh(phi)2bpy]³⁺, Nutlin-3 to upregulate p53, and subsequently irradiated to induce oxidative genomic stress. To covalently tether p53 interacting with DNA, the cells were fixed with disuccinimidyl glutarate and formaldehyde. The nuclei of the harvested cells were isolated, sonicated, and immunoprecipitated using magnetic beads conjugated with a monoclonal p53 antibody. The purified immounoprecipiated DNA was then quantified via qPCR and genomic sequencing. Overall, the ChIP results were significantly varied over ten experimental trials, but one trend is observed overall: greater variation of p53 occupancy is observed in response elements from which oxidative dissociation would be expected, while significantly less change in p53 occupancy occurs for response elements from which oxidative dissociation would not be anticipated.

The chemical oxidation of transcription factor p53 via DNA CT was also investigated with respect to the protein at the amino acid level. Transcription factor p53 plays a critical role in the cellular response to stress stimuli, which may be modulated through the redox modulation of conserved cysteine residues within the DNA-binding domain. Residues within p53 that enable oxidative dissociation are herein investigated. Of the 8 mutants studied by electrophoretic mobility shift assay (EMSA), only the C275S mutation significantly decreased the protein affinity (KD) for the Gadd45 response element. EMSA assays of p53 oxidative dissociation promoted by photoexcitation of anthraquinone-tethered Gadd45 oligonucleotides were used to determine the influence of p53 mutations on oxidative dissociation; mutation to C275S severely attenuates oxidative dissociation while C277S substantially attenuates dissociation. Differential thiol labeling was used to determine the oxidation states of cysteine residues within p53 after DNA-mediated oxidation. Reduced cysteines were iodoacetamide labeled, while oxidized cysteines participating in disulfide bonds were ¹³C₂D₂-iodoacetamide labeled. Intensities of respective iodoacetamide-modified peptide fragments were analyzed using a QTRAP 6500 LC-MS/MS system, quantified with Skyline, and directly compared. A distinct shift in peptide labeling toward ¹³C₂D₂-iodoacetamide labeled cysteines is observed in oxidized samples as compared to the respective controls. All of the observable cysteine residues trend toward the heavy label under conditions of DNA CT, indicating the formation of multiple disulfide bonds potentially among the C124, C135, C141, C182, C275, and C277. Based on these data it is proposed that disulfide formation involving C275 is critical for inducing oxidative dissociation of p53 from DNA.