Bridging the Transition from Mesoscale to Microscale Turbulence in Numerical Weather Prediction Models

With a focus towards developing multiscale capabilities in numerical weather prediction models, the specific problem of the transition from the mesoscale to the microscale is investigated. For that purpose, idealized one-way nested mesoscale to large-eddy simulation (LES) experiments were carried out using the Weather Research and Forecasting model framework. It is demonstrated that switching from one-dimensional turbulent diffusion in the mesoscale model to three-dimensional LES mixing does not necessarily result in an instantaneous development of turbulence in the LES domain. On the contrary, very large fetches are needed for the natural transition to turbulence to occur. The computational burden imposed by these long fetches necessitates the development of methods to accelerate the generation of turbulence on a nested LES domain forced by a smooth mesoscale inflow. To that end, four new methods based upon finite amplitude perturbations of the potential temperature field along the LES inflow boundaries are developed, and investigated under convective conditions. Each method accelerated the development of turbulence within the LES domain, with two of the methods resulting in a rapid generation of production and inertial range energy content associated to microscales that is consistent with non-nested simulations using periodic boundary conditions. The cell perturbation approach, the simplest and most efficient of the best performing methods, was investigated further under neutral and stable conditions. Successful results were obtained in all the regimes, where satisfactory agreement of mean velocity, variances and turbulent fluxes, as well as velocity and temperature spectra, was achieved with reference non-nested simulations. In contrast, the non-perturbed LES solution exhibited important energy deficits associated to a delayed establishment of fully-developed turbulence. The cell perturbation method has negligible computational cost, significantly accelerates the generation of realistic turbulence, and requires minimal parameter tuning, with the necessary information relatable to mean inflow conditions provided by the mesoscale solution.     

On the possibility of a kinetic fractionation of nitrogen stable isotopes during natural diamond growth

In a diamond from New South Wales (Australia), cubic and octahedral growth sectors, as identified by cathodoluminescence (CL), show slight differences in N-contents of 29 and 42 ppm respectively but no significant differences in either δ¹³C, δ¹⁵N and nitrogen aggregation state with values at +1.96‰, +19.4‰, and 25% Type IaAB aggregation, respectively.Two gem cubes from the Orapa kimberlite (Botswana) were studied by CL revealing a nonfaceted cubic growth. Accordingly, nine other gem cubes were combusted and yielded δ¹³C-values from -5.33‰ to -6.63‰, δ¹⁵N from -1.0‰ to -5.5‰, and nitrogen contents from 914 to 1168 ppm, with nitrogen aggregation state being only Type IaA (zero % B). The gem cubes show striking similarities to fibrous/coated diamonds, not only in both δ¹³C ranges (less than 3‰ from -5 to -8‰), but also in the high levels of nitrogen (≈ 1000 ppm), suggesting that the two diamond types are related. Additionally, no δ¹⁵N variation was detected between the cube and octahedral growth sectors of the Australian diamond, in the cube sectors of the nine gem cubes from Botswana, nor in fibrous/coated diamonds previously studied. These analyses contrast with an earlier study on a synthetic diamond, which reported a strong kinetic fractionation of N-isotopes of about 40‰ between cube and octahedral growth. The present evidence, therefore, suggests that kinetic fractionation of N-isotopes does not operate during natural diamond formation.     

Ferropericlase—a lower mantle phase in the upper mantle

Experiments on compositions along the join MgO–NaA³⁺Si₂O₆ (A=Al, Cr, Fe³⁺) show that sodium can be incorporated into ferropericlase at upper mantle pressures in amounts commonly found in natural diamond inclusions. These results, combined with the observed mineral parageneses of several diamond inclusion suites, establish firmly that ferropericlase exists in the upper mantle in regions with low silica activity. Such regions may be carbonated dunite or stalled and degassed carbonatitic melts. Ferropericlase as an inclusion in diamond on its own is not indicative of a lower mantle origin or of a deep mantle plume. Coexisting phases have to be taken into consideration to decide on the depth of origin. The composition of olivine will indicate an origin from the upper mantle or border of the transition zone to the lower mantle and whether it coexisted with ferropericlase in the upper mantle or as ringwoodite. The narrow and flat three phase loop at the border transition zone—lower mantle together with hybrid peridotite plus eclogite/sediments provides an explanation for the varying and Fe-rich nature of the diamond inclusion suite from Sao Luiz, Brazil.     

Effects of topography on lofting gravity flows: Implications for the deposition of deep-water massive sands

Hyperpycnal flows are generated in the marine environment by sediment-laden fresh water discharge into the ocean. They frequently form at river mouths and are also generated in proximal ice-melting settings and are thought to be responsible for transporting a large proportion of suspended river sediment onto and off the continental shelf. Hyperpycnal flows are an example of gravity currents that display reversing buoyancy. This phenomenon is generated by the fresh water interstitial fluid being less dense than that of the ambient seawater. Thus after sufficient particles are sedimented the flow can become positively buoyant and loft, forming a rising plume. Here we present results from physical scale-modelling experiments of lofting gravity currents upon interaction with topography. Topography, in the form of a vertical obstacle, triggered a localised lofting zone on its upstream side. This lofting zone was maintained in a fixed position until the bulk density of the flow had reduced enough to allow lofting along its entire length. The obstructed lofting zone is associated with a sharp increase in deposit thickness. By inference these experimentally established lofting dynamics are applied to improve understanding of the potential for hyperpycnal flows to deposit deep-water massive sands. This study provides a depositional mechanism by which large volumes of sand can be deposited in the absence of traction and the fines removed, leaving thick deposits of structureless sand with a low percentage of mud. This conceptual model for the first time provides a framework by which the geometries of certain deep-water massive sands may be predicted within specific depositional and basinal settings. This is crucial to our understanding of massive sand deposits in modern and ancient turbiditic systems and in the commercial evaluation of hydrocarbon potential of such sedimentary successions.     

Cumulative impact mapping: Advances,relevance and limitations to marine management and conservation,using Canada's Pacific waters as a case study

Analysis of cumulative human impacts in the marine environment is still in its infancy but developing rapidly. In this study, existing approaches were expanded upon, aiming for a realistic consideration of cumulative impacts at a regional scale. Thirty-eight human activities were considered, with each broken down according to stressor types and a range of spatial influences. To add to the policy relevance, existing stressors within and outside of conservation areas were compared. Results indicate the entire continental shelf of Canada's Pacific marine waters is affected by multiple human activities at some level. Commercial fishing, land-based activities and marine transportation accounted for 57.0%, 19.1%, and 17.7% of total cumulative impacts, respectively. Surprisingly, most areas with conservation designations contained higher impact scores than the mean values of their corresponding ecoregions. Despite recent advances in mapping cumulative impacts, many limitations remain. Nonetheless, preliminary analyses such as these can provide information relevant to precautionary management and conservation efforts.     

Bounds on the complex conductivity of geophysical mixtures

Following previous work on bounds for complex dielectrics, bounds on the complex conductivity of a mixture of two isotropic components can be developed which are independent of any special assumption concerning the geometry of the mixture. If certain broad restrictions are assumed, such as isotropy of the mixture, then the bounds can be made more restrictive. These bounds reveal the range of the induced polarization response which can be caused by a mixture of two materials of known complex conductivity. The bounds can also be generalized for spectral responses. The bounds are conservative lithologically in the sense that many of the special models corresponding to boundary responses have lithological counterparts. The chief use for the given bounds is to gain insight into the nature of the induced polarization response. It is also possible to use the bounds to estimate the volume fractions of the components. We illustrate how this is done for the case of a general anisotropic medium.     

Minimum work, fault activity and the growth of critical wedges in fold and thrust belts

Many studies of critical wedges treat the interior of the wedge as continuous and do not address the manner in which it grows from the undeformed state to a typical imbricate wedge. In this paper we present a 2D kinematic–mechanical model which attempts to explain the development of a critical wedge in a fold and thrust belt in terms of both gravitational and frictional work. In the undeformed model a series of thrust faults are defined which have the potential to take up an external displacement. The active fault at a given time is that which minimizes gravitational and frictional work as a result of displacement. Displacement on the active fault causes a change in topography and deformation of other faults which may favour an alternative fault at the next time step. The model is a mixed Lagrangian–Eulerian scheme in which the upper surface, in addition to being deformed, is also subject to erosion, transport and sedimentation. The model predicts propagation of thrust fault activity towards the foreland through time as a result of increasing topographic (gravitational) loads and frictional work on deformed hinterland faults. As the zone of fault activity progresses through the developing critical wedge several faults are active over time-scales of ≈1 Myr. However, a simple chronology or sequence of fault activity cannot be assumed as out-of-sequence thrusting occurs during this overall foreland propagation. The detailed spatial and temporal activity of faults is complex and reflects the interaction between the development of topography, the contrast between basal (décollement) and internal coefficients of friction and the effects of erosion and sedimentation. In particular, rates of erosion and sedimentation are found to be important controls on fault activity both spatially and temporally. Erosion, by locally removing topography above a fault, reduces gravitational and frictional work enabling continued fault activity or reactivation. Sedimentation, conversely, acts to increase gravitational and frictional work on a fault, and therefore has the potential to blanket faults and render them inactive. Model results illustrate the complex feedbacks that can exist between tectonic and surficial mass transport processes.     

Silicon control of strontium and cesium partitioning in hydroxide-weathered sediments

Cation partitioning and speciation in an aqueous soil suspension may depend on the coupling of reaction time, sorbate amount and mineral weathering reactions. These factors were varied in sediment suspension experiments to identify geochemical processes that affect migration of Sr²⁺ and Cs⁺ introduced to the subsurface by caustic high level radioactive waste (HLRW). Three glacio-fluvial and lacustrine sediments from the Hanford Site (WA, USA) were subjected to hyperalkaline (pH > 13), Na-Al-NO₃-OH solution conditions within a gradient field of (i) sorptive concentration (10⁻⁵-10⁻³ m) and (ii) reaction time (0-365 d). Strontium uptake (q_Sr) exceeded that of cesium at nearly all reaction times. Sorbent affinity for both Cs⁺ and Sr²⁺ increased with clay plus silt content at early times, but a prolonged slow uptake process was observed over the course of sediment weathering that erased the texture effect for Sr²⁺; all sediments showed similar mass normalized uptake after several months of reaction time. Strontium became progressively recalcitrant to desorption after 92 d, with accumulation and aging of neoformed aluminosilicates. Formation of Cs⁺ and Sr²⁺-containing cancrinite and sodalite was observed after 183 d by SEM and synchrotron μ-XRF and μ-XRD. EXAFS data for q_Sr ≈ 40 mmol kg⁻¹ showed incorporation of Sr²⁺ into both feldspathoid and SrCO₃(s) coordination environments after one year. Adsorption was predominant at early times and low sorbate amount, whereas precipitation, controlled largely by sediment Si release, became increasingly important at longer times and higher sorbate amount. Kinetics of contaminant desorption at pH 8 from one year-weathered sediments showed significant dependence on background cation (Ca²⁺ versus K⁺) composition. Results of this study indicate that co-precipitation and ion exchange in neoformed aluminosilicates may be an important mechanism controlling Sr²⁺ and Cs⁺ mobility in siliceous sediments impacted by hyperalkaline HLRW.     

Inter‐comparison of hydro‐climatic regimes across northern catchments: synchronicity,resistance and resilience

The higher mid‐latitudes of the Northern Hemisphere are particularly sensitive to climate change as small differences in temperature determine frozen ground status, precipitation phase, and the magnitude and timing of snow accumulation and melt. An international inter‐catchment comparison program, North‐Watch, seeks to improve our understanding of the sensitivity of northern catchments to climate change by examining their hydrological and biogeochemical responses. The catchments are located in Sweden (Krycklan), Scotland (Mharcaidh, Girnock and Strontian), the United States (Sleepers River, Hubbard Brook and HJ Andrews) and Canada (Catamaran, Dorset and Wolf Creek). This briefing presents the initial stage of the North‐Watch program, which focuses on how these catchments collect, store and release water and identify ‘types’ of hydro‐climatic catchment response. At most sites, a 10‐year data of daily precipitation, discharge and temperature were compiled and evaporation and storage were calculated. Inter‐annual and seasonal patterns of hydrological processes were assessed via normalized fluxes and standard flow metrics. At the annual‐scale, relations between temperature, precipitation and discharge were compared, highlighting the role of seasonality, wetness and snow/frozen ground. The seasonal pattern and synchronicity of fluxes at the monthly scale provided insight into system memory and the role of storage. We identified types of catchments that rapidly translate precipitation into runoff and others that more readily store water for delayed release. Synchronicity and variance of rainfall–runoff patterns were characterized by the coefficient of variation (cv) of monthly fluxes and correlation coefficients. Principal component analysis (PCA) revealed clustering among like catchments in terms of functioning, largely controlled by two components that (i) reflect temperature and precipitation gradients and the correlation of monthly precipitation and discharge and (ii) the seasonality of precipitation and storage. By advancing the ecological concepts of resistance and resilience for catchment functioning, results provided a conceptual framework for understanding susceptibility to hydrological change across northern catchments. Copyright © 2010 John Wiley & Sons, Ltd.     

Collaborative geomatics and the Mushkegowuk Cree First Nations: Fostering adaptive capacity for community-based sub-arctic natural resource management

The remote First Nation (FN) communities of the Mushkegowuk Territory on the west coast of James Bay, Ontario, Canada are currently facing increased development pressures and the imposition of a government land use planning process. The land use planning process is mandated in the Far North Act (received Royal Assent on September 23, 2010). There is a need for capacity enhancement for community-based natural resource planning and management in the Territory. A number of frameworks are emerging for addressing change brought on by resource development and building resilience to such change at the community level. Among these include the concept of adaptive capacity. In collaboration with FN community leaders, we explored the use of “collaborative geomatics” tools to foster adaptive capacity. Our action research suggests that collaborative geomatics technologies should enhance the Mushkegowuk First Nations’ adaptive capacity to address environmental and policy change by allowing them to collect and manage data collaboratively (e.g., traditional environmental knowledge, western science) to create opportunities for innovative community development, including natural resource development and management.