Buoyant melting instabilities beneath extending lithosphere: 1. Numerical models

Buoyant decompression melting instabilities in regions of partially molten upper mantle have been proposed to be an important process that might account for some characteristics of intraplate volcanism on Earth and other terrestrial planets. The instability is driven by variations in the melting rate within a partially molten layer whenever a relative decrease in density accompanies decompression melting of ascending mantle. Here, the development of buoyant decompression melting instabilities in a plane layer of passively upwelling and partially melting mantle beneath diffusely extending lithosphere is studied using numerical convection models covering a wide range of physical parameters. We find that the occurrence and nature of these instabilities in such a scenario is strongly affected by the rate of extension and melt percolation, as well as depth distribution of solid density variations arising from melt depletion. In some cases, instabilities do not occur during extension, but only develop after extension has slowed or stopped completely. This behavior creates two pulses of magma generation due to passive upwelling accompanying extension followed by the subsequent instability and is favored by a faster rate of extension, higher mantle viscosity, higher rate of melt percolation, and smaller amount of solid residuum depletion‐derived buoyancy. Larger degrees of solid density changes accompanying melt depletion can enhance the instability of partially molten mantle during extension but decrease the cumulative volume of generated melt. This kind of behavior modifies the conventional expectation of spatially and temporally correlated volcanism and extension and may lend insight into the observed increase in localized volcanic activity following Miocene Basin and Range extension in the western United States

Buoyant melting instabilities beneath extending lithosphere: 2. Linear analysis

In a companion paper, numerical models reveal that buoyant melting instabilities can occur beneath extending lithosphere for a sufficiently small mantle viscosity, extension rate, and rate of melt percolation. However, in some cases, instabilities do not develop during extension but only occur after extension slows or stops. These results are suggestive of a critical behavior in the onset of these kinds of instabilities and motivate a linear analysis to study the onset of instability in a partially melting, passively upwelling plane layer of mantle beneath extending lithosphere. The model we employ includes the effects of buoyancy arising from thermal expansion, the presence of a retained fraction of partial melt, and depletion of the solid by melt extraction. We find a critical behavior in the onset of instability controlled by melt retention buoyancy that is characterized by a “Rayleigh” number M, such that M must exceed some critical value M_(crit) which depends on the efficiency of Stokes rise of a partially molten body relative to the rate of background percolation. Comparison of this theory to the numerical results in the companion paper yields a close quantitative agreement. We also find that solid depletion buoyancy can either stabilize or destabilize a partially melting layer, depending upon both the distribution of preexisting depletion and the magnitude of density changes with depth. This theory is compared with previous studies of buoyant melting instabilities beneath mid‐ocean ridges where similar behavior was reported, and it suggests that the stability of passively upwelling, partially melting mantle underlying both narrow and wide rift settings is controlled by similar processes

Buoyant melting instabilities beneath extending lithosphere: 2. Linear analysis

In a companion paper, numerical models reveal that buoyant melting instabilities can occur beneath extending lithosphere for a sufficiently small mantle viscosity, extension rate, and rate of melt percolation. However, in some cases, instabilities do not develop during extension but only occur after extension slows or stops. These results are suggestive of a critical behavior in the onset of these kinds of instabilities and motivate a linear analysis to study the onset of instability in a partially melting, passively upwelling plane layer of mantle beneath extending lithosphere. The model we employ includes the effects of buoyancy arising from thermal expansion, the presence of a retained fraction of partial melt, and depletion of the solid by melt extraction. We find a critical behavior in the onset of instability controlled by melt retention buoyancy that is characterized by a “Rayleigh” number M, such that M must exceed some critical value M_(crit) which depends on the efficiency of Stokes rise of a partially molten body relative to the rate of background percolation. Comparison of this theory to the numerical results in the companion paper yields a close quantitative agreement. We also find that solid depletion buoyancy can either stabilize or destabilize a partially melting layer, depending upon both the distribution of preexisting depletion and the magnitude of density changes with depth. This theory is compared with previous studies of buoyant melting instabilities beneath mid‐ocean ridges where similar behavior was reported, and it suggests that the stability of passively upwelling, partially melting mantle underlying both narrow and wide rift settings is controlled by similar processes

Buoyant melting instabilities beneath extending lithosphere: 1. Numerical models

Buoyant decompression melting instabilities in regions of partially molten upper mantle have been proposed to be an important process that might account for some characteristics of intraplate volcanism on Earth and other terrestrial planets. The instability is driven by variations in the melting rate within a partially molten layer whenever a relative decrease in density accompanies decompression melting of ascending mantle. Here, the development of buoyant decompression melting instabilities in a plane layer of passively upwelling and partially melting mantle beneath diffusely extending lithosphere is studied using numerical convection models covering a wide range of physical parameters. We find that the occurrence and nature of these instabilities in such a scenario is strongly affected by the rate of extension and melt percolation, as well as depth distribution of solid density variations arising from melt depletion. In some cases, instabilities do not occur during extension, but only develop after extension has slowed or stopped completely. This behavior creates two pulses of magma generation due to passive upwelling accompanying extension followed by the subsequent instability and is favored by a faster rate of extension, higher mantle viscosity, higher rate of melt percolation, and smaller amount of solid residuum depletion‐derived buoyancy. Larger degrees of solid density changes accompanying melt depletion can enhance the instability of partially molten mantle during extension but decrease the cumulative volume of generated melt. This kind of behavior modifies the conventional expectation of spatially and temporally correlated volcanism and extension and may lend insight into the observed increase in localized volcanic activity following Miocene Basin and Range extension in the western United States

Geodynamo and mantle convection simulations on the Earth Simulator using the Yin-Yang grid

We have developed finite difference codes based on the Yin-Yang grid for the geodynamo simulation and the mantle convection simulation. The Yin-Yang grid is a kind of spherical overset grid that is composed of two identical component grids. The intrinsic simplicity of the mesh configuration of the Yin-Yang grid enables us to develop highly optimized simulation codes on massively parallel supercomputers. The Yin-Yang geodynamo code has achieved 15.2 Tflops with 4096 processors on the Earth Simulator. This represents 46% of the theoretical peak performance. The Yin-Yang mantle code has enabled us to carry out mantle convection simulations in realistic regimes with a Rayleigh number of $10^7$ including strongly temperature-dependent viscosity with spatial contrast up to $10^6$.Comment: Plenary talk at SciDAC 200

Improving gait classification in horses by using inertial measurement unit (IMU) generated data and machine learning

For centuries humans have been fascinated by the natural beauty of horses in motion and their different gaits. Gait classification (GC) is commonly performed through visual assessment and reliable, automated methods for real-time objective GC in horses are warranted. In this study, we used a full body network of wireless, high sampling-rate sensors combined with machine learning to fully automatically classify gait. Using data from 120 horses of four different domestic breeds, equipped with seven motion sensors, we included 7576 strides from eight different gaits. GC was trained using several machine-learning approaches, both from feature-extracted data and from raw sensor data. Our best GC model achieved 97% accuracy. Our technique facilitated accurate, GC that enables in-depth biomechanical studies and allows for highly accurate phenotyping of gait for genetic research and breeding. Our approach lends itself for potential use in other quadrupedal species without the need for developing gait/animal specific algorithms

JC Virus T-Antigen Regulates Glucose Metabolic Pathways in Brain Tumor Cells

Recent studies have reported the detection of the human neurotropic virus, JCV, in a significant population of brain tumors, including medulloblastomas. Accordingly, expression of the JCV early protein, T-antigen, which has transforming activity in cell culture and in transgenic mice, results in the development of a broad range of tumors of neural crest and glial origin. Evidently, the association of T-antigen with a range of tumor-suppressor proteins, including p53 and pRb, and signaling molecules, such as β-catenin and IRS-1, plays a role in the oncogenic function of JCV T-antigen. We demonstrate that T-antigen expression is suppressed by glucose deprivation in medulloblastoma cells and in glioblastoma xenografts that both endogenously express T-antigen. Mechanistic studies indicate that glucose deprivation-mediated suppression of T-antigen is partly influenced by 5′-activated AMP kinase (AMPK), an important sensor of the AMP/ATP ratio in cells. In addition, glucose deprivation-induced cell cycle arrest in the G1 phase is blocked with AMPK inhibition, which also prevents T-antigen downregulation. Furthermore, T-antigen prevents G1 arrest and sustains cells in the G2 phase during glucose deprivation. On a functional level, T-antigen downregulation is partially dependent on reactive oxygen species (ROS) production during glucose deprivation, and T-antigen prevents ROS induction, loss of ATP production, and cytotoxicity induced by glucose deprivation. Additionally, we have found that T-antigen is downregulated by the glycolytic inhibitor, 2-deoxy-D-glucose (2-DG), and the pentose phosphate inhibitors, 6-aminonicotinamide and oxythiamine, and that T-antigen modulates expression of the glycolytic enzyme, hexokinase 2 (HK2), and the pentose phosphate enzyme, transaldolase-1 (TALDO1), indicating a potential link between T-antigen and metabolic regulation. These studies point to the possible involvement of JCV T-antigen in medulloblastoma proliferation and the metabolic phenotype and may enhance our understanding of the role of viral proteins in glycolytic tumor metabolism, thus providing useful targets for the treatment of virus-induced tumors

Persistence of strong silica-enriched domains in the Earth's lower mantle

The composition of the lower mantle—comprising 56% of Earth’s volume—remains poorly constrained. Among the major elements, Mg/Si ratios ranging from ∼0.9–1.1, such as in rocky Solar-System building blocks (or chondrites), to ∼1.2–1.3, such as in upper-mantle rocks (or pyrolite), have been proposed. Geophysical evidence for subducted lithosphere deep in the mantle has been interpreted in terms of efficient mixing, and thus homogenous Mg/Si across most of the mantle. However, previous models did not consider the effects of variable Mg/Si on the viscosity and mixing efficiency of lower-mantle rocks. Here, we use geodynamic models to show that large-scale heterogeneity associated with a 20-fold change in viscosity, such as due to the dominance of intrinsically strong (Mg, Fe)SiO3–bridgmanite in low-Mg/Si domains, is sufficient to prevent efficient mantle mixing, even on large scales. Models predict that intrinsically strong domains stabilize mantle convection patterns, and coherently persist at depths of about 1,000–2,200 km up to the present-day, separated by relatively narrow up-/downwelling conduits of pyrolitic material. The stable manifestation of such bridgmanite-enriched ancient mantle structures (BEAMS) may reconcile the geographical fixity of deep-rooted mantle upwelling centres, and geophysical changes in seismic-tomography patterns, radial viscosity, rising plumes and sinking slabs near 1,000 km depth. Moreover, these ancient structures may provide a reservoir to host primordial geochemical signatures