Dimensionless measures of turbulent magnetohydrodynamic dissipation rates

The magnetic Reynolds number, R _M, is defined as the product of a characteristic scale and associated flow speed divided by the microphysical magnetic diffusivity. For laminar flows, R _M also approximates the ratio of advective to dissipative terms in the total magnetic energy equation, but for turbulent flows this latter ratio depends on the energy spectra and approaches unity in a steady state. To generalize for flows of arbitrary spectra we define an effective magnetic dissipation number,   R _M,e  , as the ratio of the advection to microphysical dissipation terms in the total magnetic energy equation, incorporating the full spectrum of scales, arbitrary magnetic Prandtl numbers, and distinct pairs of inner and outer scales for magnetic and kinetic spectra. As expected, for a substantial parameter range   R _M,e∼ O (1) ≪ R _M  . We also distinguish   R _M,e  from     where the latter is an effective magnetic Reynolds number for the mean magnetic field equation when a turbulent diffusivity is explicitly imposed as a closure. That   R _M,e  and     approach unity even if   R _M≫ 1  highlights that, just as in hydrodynamic turbulence, energy dissipation of large-scale structures in turbulent flows via a cascade can be much faster than the dissipation of large-scale structures in laminar flows. This illustrates that the rate of energy dissipation by magnetic reconnection is much faster in turbulent flows, and much less sensitive to microphysical reconnection rates compared to laminar flows.     

Regional climate model projections for the State of Washington

Global climate models do not have sufficient spatial resolution to represent the atmospheric and land surface processes that determine the unique regional climate of the State of Washington. Regional climate models explicitly simulate the interactions between the large-scale weather patterns simulated by a global model and the local terrain. We have performed two 100-year regional climate simulations using the Weather Research and Forecasting (WRF) model developed at the National Center for Atmospheric Research (NCAR). One simulation is forced by the NCAR Community Climate System Model version 3 (CCSM3) and the second is forced by a simulation of the Max Plank Institute, Hamburg, global model (ECHAM5). The mesoscale simulations produce regional changes in snow cover, cloudiness, and circulation patterns associated with interactions between the large-scale climate change and the regional topography and land-water contrasts. These changes substantially alter the temperature and precipitation trends over the region relative to the global model result or statistical downscaling. To illustrate this effect, we analyze the changes from the current climate (1970–1999) to the mid twenty-first century (2030–2059). Changes in seasonal-mean temperature, precipitation, and snowpack are presented. Several climatological indices of extreme daily weather are also presented: precipitation intensity, fraction of precipitation occurring in extreme daily events, heat wave frequency, growing season length, and frequency of warm nights. Despite somewhat different changes in seasonal precipitation and temperature from the two regional simulations, consistent results for changes in snowpack and extreme precipitation are found in both simulations.     

Autoxidation of N(III), S(IV), and other Species in Frozen Solution – A Possible Pathway for Enhanced Chemical Transformation in Freezing Systems

The chemistry of glycolaldehyde (hydroxyacetaldehyde) relevant to the troposphere has been investigated using UV absorption spectrometry and FTIR absorption spectrometry in an environmental chamber. Quantitative UV absorption spectra have been obtained for the first time. The UV spectrum peaks at 277 nm with a maximum cross section of (5.5± 0.7)×10^–20 cm² molecule^–1. Studies of the ultraviolet photolysis of glycolaldehyde ( = 285 ± 25 nm) indicated that the overall quantum yield is > 0.5 in one bar of air, with the major products being CH₂OH and HCO radicals. Rate coefficients for the reactions of Cl atoms and OH radicals with glycolaldehyde have been determined to be (7.6± 1.5)×10^–11 and (1.1± 0.3)×10^–11 cm³ molecule^–1 s^–1, respectively, in good agreement with the only previous study. The lifetime of glycolaldehyde in the atmosphere is about 1.0 day for reaction with OH, and > 2.5 days for photolysis, although both wet and dry deposition should also be considered in future modeling studies.     

Differential unroofing within the central metasedimentary Belt of the Grenville Orogen: constraints from 40Ar/39Ar thermochronology

An ⁴⁰Ar/³⁹Ar thermochronological investigation of upper greenschist to granulite facies gneiss, amphibolite and marble was conducted in the Central Metasedimentary Belt (CMB), Ontario, to constrain its cooling history. Incremental ⁴⁰Ar/³⁹Ar release spectra indicate that substantial differential unroofing occurred in the CMB between 1000 and 600 Ma. A consistent pattern of significantly older hornblende and phlogopite ⁴⁰Ar/³Ar cooling ages on the southeast sides of major northeast striking shear zones is interpreted to reflect late displacement due to extensional deformation. Variations in hornblende ⁴⁰Ar/³⁹Ar age plateaus exceeding 200 Ma occur over distances less than 50 km with major age discontinuities occurring across the Robertson Lake shear zone and the Sharbot Lake mylonite zone which separate the Sharbot Lake terrane from the Elzevir and Frontenac terranes. Extensional displacements of up to 14 km are inferred between the Frontenac and Elzevir terranes of the CMB. No evidence for significant post argon-closure vertical displacement is indicated in the vicinity of the Perth Road mylonite within the Frontenac terrane. Variations of nearly 100 Ma in phlogopite ⁴⁰Ar/³⁹Ar plateau ages occur in undeformed marble on either side of the Bancroft Shear Zone. Phlogopites from sheared and mylonitized marble within the shear zone yield ⁴⁰Ar/³⁹Ar diffusional loss profiles, but have older geologically meaningless ages thought to reflect incorporation of excess argon. By 900 Ma, southeast directed extension was occurring throughout the CMB, possibly initiated along previous zones of compressional shearing. An easterly migration of active zones of extension is inferred, possibly related to an earlier, overall easterly migration of active zones of regional thrusting and easterly migration of an ancient subduction zone. The duration of extensional shearing is not well constrained, but must have ceased before 600 Ma as required by the deposition of overlying undeformed Cambrian and/or Ordovician sedimentary rocks.Contribution No. 481 from the Mineralogical Laboratory, University of Michigan     

On the Importance of High-Resolution in Large-Scale Ocean Models

Eddying global ocean models are now routinely used for ocean prediction, and the value-added of a better representation of the observed ocean variability and western boundary currents at that resolution is currently being evaluated in climate models. This overview article begins with a brief summary of the impact on ocean model biases of resolving eddies in several global ocean–sea ice numerical simulations. Then, a series of North and Equatorial Atlantic configurations are used to show that an increase of the horizontal resolution from eddy-resolving to submesoscale-enabled together with the inclusion of high-resolution bathymetry and tides significantly improve the models' abilities to represent the observed ocean variability and western boundary currents. However, the computational cost of these simulations is extremely large, and for these simulations to become routine, close collaborations with computer scientists are essential to ensure that numerical codes can take full advantage of the latest computing architecture.     

1988 Compilation of Elemental Concentration Data for USGS Basalt BCR-1

Concentration data on up to 82 individual constituents in USGS Basalt BCR-1 have been collected from 1395 journal articles and technical reports. These data are summarized in consensus (mean) values with uncertainties expressed as ± one standard deviation. Mean values are also calculated as a function of analytical procedure and all raw data are given in the tables. Recommended values are proposed based upon data criteria used by NIST (National Institute of Standards and Technology, formerly the National Bureau of Standards or NBS).     

Comparison of garnet-biotite, calcite-graphite, and calcite-dolomite thermometry in the Grenville Orogen; Ontario, Canada

The Elzevir Terrane of the Grenville Orogen in southern Ontario contains metapelites and abundant graphitic marbles that were regionally metamorphosed from the upper greenschist to upper amphibolite facies. Comparative thermometry was undertaken with widely used calibrations for the systems garnet-biotite, calcite-dolomite, and calcite-graphite. Temperatures that are obtained from matrix biotites paired with prograde garnet near-rim analyses are usually consistent with those determined using calcite-graphite thermometry. However, calcite-graphite thermometry occasionally yields low temperatures due to lack of equilibration of anomalously light graphite. Application of calcite-graphite and garnet-biotite systems may yield temperatures up to 70 °C higher than calcite-dolomite in amphibolite facies rocks. Calcite-dolomite temperatures most closely approach those from calcite-graphite and garnet-biotite when the samples contain a single generation of dolomite and calcite grains contain no visible dolomite exsolution lamellae. However, some of these samples yield temperatures considerably lower than temperatures calculated from calcite-graphite and garnet-biotite thermometry, indicating that the calcite-dolomite thermometer may have been partially reset during retrogression. Estimated peak metamorphic temperatures of regional metamorphism between Madoc (upper greenschist facies) and Bancroft (upper amphibolite facies) range from 500 to 650 °C. These results place the chlorite-staurolite isograd at 540 °C, the kyanite-sillimanite isograd at 590 °C, and the sillimanite-K-feldspar isograd at 650 °C. Although each thermometer may have an absolute uncertainty of as much as ±50 °C, the 50 to 60 °C temperature differences between the isograds are probably accurate to 10 to 20 °C. An incomplete picture of the thermal gradients can result from the application of only one thermometer in a given area. Simultaneous application of several systems allows one to recognize and overcome the inherent limitations of each thermometer. Received: 26 March 1997 / Accepted: 15 April 1998     

Impact of climate change on the hydrogeology of two basins in northern France

This study presents an analysis of climate-change impacts on the water resources of two basins located in northern France, by integrating four sources of uncertainty: climate modelling, hydrological modelling, downscaling methods, and emission scenarios. The analysis focused on the evolution of the water budget, the river discharges and piezometric heads. Seven hydrological models were used, from lumped rainfall-discharge to distributed hydrogeological models, and led to quite different estimates of the water-balance components. One of the hydrological models, CLSM, was found to be unable to simulate the increased water stress and was, thus, considered as an outlier even though it gave fair results for the present day compared to observations. Although there were large differences in the results between the models, there was a marked tendency towards a decrease of the water resource in the rivers and aquifers (on average in 2050 about ?14 % and ?2.5 m, respectively), associated with global warming and a reduction in annual precipitation (on average in 2050 +2.1 K and ?3 %, respectively). The uncertainty associated to climate models was shown to clearly dominate, while the three others were about the same order of magnitude and 3–4 times lower. In terms of impact, the results found in this work are rather different from those obtained in a previous study, even though two of the hydrological models and one of the climate models were used in both studies. This emphasizes the need for a survey of the climatic-change impact on the water resource.     

Numerical modeling of nonlinear water waves over heterogeneous porous beds

The transformation of nonlinear water waves over porous beds is studied by applying a numerical model based on Chen's [2006. Fully nonlinear Boussinesq-type equations for waves and currents over porous beds. Journal of Engineering Mechanics, 132:2, 220–230] Boussinesq-type equations for highly nonlinear waves on permeable beds. The numerical model uses a high-order time-marching solution and fourth-order finite-difference schemes for discretization of first-order spatial derivatives to obtain a computational accuracy consistent with the model equations. By forcing the wave celerity and spatial porous-damping rate of the linearized model to match the exact linear theory for horizontal porous bed over a prescribed range of relative depths, the values of the model parameters are optimally determined. Numerical simulations of the damped wave propagation over finite-thickness porous layer demonstrate the accuracy of both the numerical model and governing equations, which have been shown by prior theoretical analyses to be accurate for both nominal and thick porous layers. These simulations also elucidate on the significance of the higher-order porous-damping terms and the influence of the hydraulic parameters. Application of the model to the simulation of the wave field around a laboratory-scale submerged porous mound provides a measure of its capability, as well as useful insight into the scaling of the porous-resistance coefficients. For application to heterogeneous porous beds, the assumption of weak spatial variation of the porous resistance is examined using truncated forms of the governing equations. The results indicate that the complete set of Boussinesq-type equations is applicable to porous beds of nonhomogeneous makeup.