A three-dimensional,iterative mapping procedure for the implementation of an ionosphere-magnetosphere anisotropic Ohm’s law boundary condition in global magnetohydrodynamic simulations

The mathematical formulation of an iterative procedure for the numerical implementation of an ionosphere-magnetosphere (IM) anisotropic Ohm’s law boundary condition is presented. The procedure may be used in global magnetohydrodynamic (MHD) simulations of the magnetosphere. The basic form of the boundary condition is well known, but a well-defined, simple, explicit method for implementing it in an MHD code has not been presented previously. The boundary condition relates the ionospheric electric field to the magnetic field-aligned current density driven through the ionosphere by the magnetospheric convection electric field, which is orthogonal to the magnetic field B, and maps down into the ionosphere along equipotential magnetic field lines. The source of this electric field is the flow of the solar wind orthogonal to B. The electric field and current density in the ionosphere are connected through an anisotropic conductivity tensor which involves the Hall, Pedersen, and parallel conductivities. Only the height-integrated Hall and Pedersen conductivities (conductances) appear in the final form of the boundary condition, and are assumed to be known functions of position on the spherical surface R=R₁ representing the boundary between the ionosphere and magnetosphere. The implementation presented consists of an iterative mapping of the electrostatic potential , the gradient of which gives the electric field, and the field-aligned current density between the IM boundary at R=R₁ and the inner boundary of an MHD code which is taken to be at R₂>R₁. Given the field-aligned current density on R=R₂, as computed by the MHD simulation, it is mapped down to R=R₁ where it is used to compute by solving the equation that is the IM Ohm’s law boundary condition. Then is mapped out to R=R₂, where it is used to update the electric field and the component of velocity perpendicular to B. The updated electric field and perpendicular velocity serve as new boundary conditions for the MHD simulation which is then used to compute a new field-aligned current density. This process is iterated at each time step. The required Hall and Pedersen conductances may be determined by any method of choice, and may be specified anew at each time step. In this sense the coupling between the ionosphere and magnetosphere may be taken into account in a self-consistent manner.     

Numerical analyses of deformability tests in jointed rock — “Joint Perturbation” and “No Tension” Finite Element solutions

SummaryNumerical Analyses of Deformability Tests in Jointed Rock — Joint Perturbation and No Tension Finite Element Solutions Sound design of structures in or upon rock requires a thorough knowledge of the rock mass deformability. Field test data interpretation has generally relied upon analytical or Finite Element linear elastic solutions. However, when testing in jointed rock, these can no longer be readily used since the medium cannot resist the tension induced upon loading. Accordingly, other solutions are required. If the discontinuities of the rock mass can all be mapped and their surface properties determined, the Finite Element model will use a Joint Perturbation solution where the joint elements have variable stiffnesses to account for movements along the fractures. If the discontinuities cannot be satisfactorily mapped, a global approach is provided by a No Tension analysis of the stress transfer type. The two techniques are compared in the case of bore hole jack deformability tests and found to agree remarkably well, indicating an appreciable difference from results of linear elasticity.

---

ZusammenfassungNumerische Analyse von Verformungsversuchen in geklüfteten Fels — Lösungen mittels Kluft-Perturbation und No Tension Analysis auf Grund der Methode der Finiten Elemente Der Entwurf von Bauwerken in oder auf Fels erfordert eine genaue Kenntnis des Formänderungsverhaltens der Felsmasse. Im allgemeinen basiert die Interpretation der Ergebnisse von Feldversuchen auf den Lösungen des zugehörigen Randwertproblems der linearen Elastostatik mittels klassischer analytischer Methoden oder der Methode der Finiten Elemente. Wenn die Versuche aber in geklüftetem Fels durchgeführt werden, müssen die theoretischen Methoden modifiziert werden, da das Medium keine Zugspannungen aufnehmen kann. Zwei Fälle können unterschieden werden. Erstens: die Kluftflächen der Felsmasse können geometrisch erfaßt und ihre Oberflächeneigenschaften können bestimmt werden; in diesem Falle kann die Finite-Elemente-Lösung eine Kluft-Perturbations-Methode verwenden, in welcher spezielle Kluft-Elemente verschiedener Steifigkeiten normal und tangential zu den Kluftflächen die Bewegungen entlang derselben bestimmen. Zweitens: die Diskontinuitäten der Felsmasse können nicht hinreichend beschrieben werden; für diesen Fall ist vorgekehrt durch die Finite-Element-Analyse vom Typ der Spannungsumlagerung unter Ausschaltung der Zugspannungen (No Tension Analysis), um das globale Verhalten der geklüfteten Felsmasse zu erfassen. Die Lösungen dieser beiden Methoden wurden für den Fall der Bohrloch-Verformungs-Tests verglichen; der Vergleich zeigt ausgezeichnete Übereinstimmung, im Gegensatz zu den Lösungen als lineares Elastizitätsproblem.

---

RésuméLa méthode des Eléments Finis appliquée aux résultats d'essais de vérins cylindriques dans les roches fissurées — solutions dites Joint Pertubé et Sans Traction La construction de structures sur ou dans des roches nécessite une connaissance approfendie de la déformabilité du massif rocheux. L'analyse des résultats d'essais en place se fait généralement à l'aide de solutions analytiques ou par Eléments Finis, qui supposent une élasticité linéaire du milieu. De fait, dans les roches fissurées, ces solutions sont inadéquates puisque le milieu ne peut résister aux sollicitations de traction. Il s'agit donc d'utiliser d'autres méthodes d'analyse. Lors qu'il est possible de déterminer la géométrie et les propriétés de toutes les discontinuités, l'on pourra utiliser une solution par Eléments Finis, dite Joint Perturbé dans laquelle la rigidité de surface du joint (ou fissure) peut varier lorsque des movements prennent place au niveau du joint. Si cela est impossible, il faudra utiliser une approche globale dite Sans Traction qui redistribue les contraintes de traction et les dissipe. On compare ces deux types de solutions dans le cas des essais de vérins cylindriques unidirectionnels. Leurs résultats sont identiques et ont nettement différents de ceux obtenus avec une solution élastique linéaire.

---

---

Research supported by the National Aeronautics and Space Administration — Contract NAS 8-21432.

---

With 8 Figures     

Evolving metasomatic agent in the Northern Andean subduction zone, deduced from magma composition of the long-lived Pichincha volcanic complex (Ecuador)

Geochemical studies of long-lived volcanic complexes are crucial for the understanding of the nature and composition of the subduction component of arc magmatism. The Pichincha Volcanic Complex (Northern Andean Volcanic Zone) consists of: (1) an old, highly eroded edifice, the Rucu Pichincha, whose lavas are mostly andesites, erupted from 1,100 to 150 ka; and (2) a younger, essentially dacitic, Guagua Pichincha composite edifice, with three main construction phases (Basal Guagua Pichincha, Toaza, and Cristal) which developed over the last 60 ka. This structural evolution was accompanied by a progressive increase of most incompatible trace element abundances and ratios, as well as by a sharp decrease of fluid-mobile to fluid-immobile element ratios. Geochemical data indicate that fractional crystallization of an amphibole-rich cumulate may account for the evolution from the Guagua Pichincha andesites to dacites. However, in order to explain the transition between the Rucu Pichincha andesites and Guagua Pichincha dacites, the mineralogical and geochemical data indicate the predominance of magma mixing processes between a mafic, trace-element depleted, mantle-derived end-member, and a siliceous, trace-element enriched, adakitic end-member. The systematic variation of trace element abundances and ratios in primitive samples leads us to propose that the Rucu Pichincha magmas came from a hydrous-fluid metasomatized mantle wedge, whereas Guagua Pichincha magmas are related to partial melting of a siliceous-melt metasomatized mantle. This temporal evolution implies a change from dehydration to partial melting of the slab, which may be associated with an increase in the geothermal gradient along the slab due to the presence of the subducted Carnegie Ridge at the subduction system. This work emphasizes the importance of studying arc-magma systems over long periods of time (of at least 1 million of years), in order to evaluate the potential variations of the slab contribution into the mantle source of the arc magmatism.     

Space- and ground-based measurements of sulphur dioxide emissions from Turrialba Volcano (Costa Rica)

Remotely sensed measurements of sulphur dioxide (SO₂) emitted by Turrialba Volcano (Costa Rica) are reported for the period September 2009–January 2011. These measurements were obtained using images from Advanced Spaceborne Thermal Emission and Reflexion radiometer, Ozone Monitoring Instrument and a ground-based UV camera. These three very different instruments provide flux measurements in good agreement with each other, which demonstrate that they can be integrated for monitoring SO₂ fluxes. Fluxes from Turrialba increased fourfold in January 2010, following a phreatic explosion that formed a degassing vent in the W crater of Turrialba. Since then, the SO₂ flux has remained high (30–50?kg/s) but seems to be showing a slowly decreasing trend. We interpret this evolution as the start of open vent degassing from a recently intruded magma body. The opening of the degassing vent decreased the confining pressure of the magma body and allowed the gases to bypass the hydrothermal system.     

Transient climate changes in a perturbed parameter ensemble of emissions-driven earth system model simulations

We describe results from a 57-member ensemble of transient climate change simulations, featuring simultaneous perturbations to 54 parameters in the atmosphere, ocean, sulphur cycle and terrestrial ecosystem components of an earth system model (ESM). These emissions-driven simulations are compared against the CMIP3 multi-model ensemble of physical climate system models, used extensively to inform previous assessments of regional climate change, and also against emissions-driven simulations from ESMs contributed to the CMIP5 archive. Members of our earth system perturbed parameter ensemble (ESPPE) are competitive with CMIP3 and CMIP5 models in their simulations of historical climate. In particular, they perform reasonably well in comparison with HadGEM2-ES, a more sophisticated and expensive earth system model contributed to CMIP5. The ESPPE therefore provides a computationally cost-effective tool to explore interactions between earth system processes. In response to a non-intervention emissions scenario, the ESPPE simulates distributions of future regional temperature change characterised by wide ranges, and warm shifts, compared to those of CMIP3 models. These differences partly reflect the uncertain influence of global carbon cycle feedbacks in the ESPPE. In addition, the regional effects of interactions between different earth system feedbacks, particularly involving physical and ecosystem processes, shift and widen the ESPPE spread in normalised patterns of surface temperature and precipitation change in many regions. Significant differences from CMIP3 also arise from the use of parametric perturbations (rather than a multimodel ensemble) to represent model uncertainties, and this is also the case when ESPPE results are compared against parallel emissions-driven simulations from CMIP5 ESMs. When driven by an aggressive mitigation scenario, the ESPPE and HadGEM2-ES reveal significant but uncertain impacts in limiting temperature increases during the second half of the twenty-first century. Emissions-driven simulations create scope for development of errors in properties that were previously prescribed in coupled ocean–atmosphere models, such as historical CO₂ concentrations and vegetation distributions. In this context, historical intra-ensemble variations in the airborne fraction of CO₂ emissions, and in summer soil moisture in northern hemisphere continental regions, are shown to be potentially useful constraints, subject to uncertainties in the relevant observations. Our results suggest that future climate-related risks can be assessed more comprehensively by updating projection methodologies to support formal combination of emissions-driven perturbed parameter and multi-model earth system model simulations with suitable observational constraints. This would provide scenarios underpinned by a more complete representation of the chain of uncertainties from anthropogenic emissions to future climate outcomes.     

The Multiple Sources and Patterns of Methane inNorth Sea Waters

The methane concentration in the atmosphere andsurface water was surveyed along 58° N acrossthe North Sea. In addition, the vertical methanedistribution in the water column was determined at sixstations along the transect. The methane contents ofthe surface water as well as in the water column wereextremely inhomogeneous. Input by freshwater fromriver discharge and injection of methane from thesediment were both observed. The survey continued fromthe western side of the North Sea to the Elbe Riverestuary. The Elbe River appears to have low methaneconcentrations compared to other European rivers, itsaverage input into the North Sea is estimated to be70 nmol s-1 of methane. Near 58° N,1°40' E, an abandoned drill site releases about 25 % ofthe North Sea's emission of methane to the atmosphere.The advective methane transport induced by watercirculation was assessed for May 16, 1994, using a 3-DNorth Sea circulation model. For the period of thissurvey, the North Sea's source strength foratmospheric methane is estimated using in situwind velocities. In comparison to the advectivetransport by the water circulation, the gas flux tothe atmosphere appears to be the dominant sink ofNorth Sea methane. This flux is estimated to bebetween 1500 · 106 mol a-1 and 3100 ·106mol a-1, depending on the relationbetween wind speed and gas transfer velocity.     

Shapes of chondrules determined from the petrofabric of the CR2 chondrite NWA 801

The approximately spherical shapes of chondrules has long been attributed to surface tension acting on ~1 mm melt droplets that formed and cooled in the microgravity field of the solar nebula. However, chondrule shapes commonly depart significantly from spherical. In this study, 109 chondrules in a sample of CR2 chondrite NWA 801 were imaged by X-ray computed tomography and best-fitted to ellipsoids. The analysis confirms that many chondrules are indeed not spherical, and also that the chondrules’ collective shape fabric records a definite 13% compaction in the host meteorite. Dehydration of phyllosilicates within chondrules may account for that strain. However, retro-deforming all chondrules shows that a large majority were already far from spherical prior to accretion. Possible models for these initial shapes include prior deformation of individual chondrules in earlier hosts, and, as suggested by previous authors, rotation of chondrules as they were solidifying, and/or “streaming” of molten chondrules by their differential velocities with their gaseous hosts after melting. More in situ 3-D work such as this study on a variety of unequilibrated chondrites, combined with detailed structural petrography, should help further constrain these models and refine our understanding of chondrite formation.     

Institutional Adaptation of Water Resource Infrastructures to Climate Change in Eastern Ontario

Institutional barriers and bridges to local climate change impacts adaptation affecting small rural municipalities and Conservation Authorities (CAs are watershed agencies) in Eastern Ontario (Canada) are examined, and elements of a community-based adaptation strategy related to water infrastructures are proposed as a case-study in community adaptation to climate change. No general water scarcity is expected for the region even under unusually dry weather scenarios. Localized quantity and quality problems are likely to occur especially in groundwater recharge areas. Some existing institutions can be relied on by municipalities to build an effective adaptation strategy based on a watershed/region perspective, on their credibility, and on their expertise. Windows of opportunity or framing issues are offered at the provincial level, the most relevant one in a federal state, by municipal emergency plan requirements and pending watershed source water protection legislation. Voluntary and soon to be mandated climate change mitigation programs at the federal level are other ones.     

Methodological challenges to confirmatory latent variable models of social vulnerability

Natural Hazards - Socially vulnerable communities experience disproportionately negative outcomes following natural disasters and underscoring a need for well-validated measures to identify those...