Intercomparison of measurements of bulk snow density and water equivalent of snow cover with snow core samplers: Instrumental bias and variability induced by observers

Manually collected snow data are often considered as ground truth for many applications such as climatological or hydrological studies. However, there are many sources of uncertainty that are not quantified in detail. For the determination of water equivalent of snow cover (SWE), different snow core samplers and scales are used, but they are all based on the same measurement principle. We conducted two field campaigns with 9 samplers commonly used in observational measurements and research in Europe and northern America to better quantify uncertainties when measuring depth, density and SWE with core samplers. During the first campaign, as a first approach to distinguish snow variability measured at the plot and at the point scale, repeated measurements were taken along two 20 m long snow pits. The results revealed a much higher variability of SWE at the plot scale (resulting from both natural variability and instrumental bias) compared to repeated measurements at the same spot (resulting mostly from error induced by observers or very small scale variability of snow depth). The exceptionally homogeneous snowpack found in the second campaign permitted to almost neglect the natural variability of the snowpack properties and focus on the separation between instrumental bias and error induced by observers. Reported uncertainties refer to a shallow, homogeneous tundra-taiga snowpack less than 1 m deep (loose, mostly recrystallised snow and no wind impact). Under such measurement conditions, the uncertainty in bulk snow density estimation is about 5% for an individual instrument and is close to 10% among different instruments. Results confirmed that instrumental bias exceeded both the natural variability and the error induced by observers, even in the case when observers were not familiar with a given snow core sampler.     

Interaction between lateral sorting in river bends and vertical sorting in dunes

Sediment is sorted in river bends under the influence of gravity that pulls the heavier grains downslope and secondary flow that drags the finer grains upslope. Furthermore, when dunes are present, sediment is also sorted vertically at the dune lee side. However, sorting functions are poorly defined, since the relation to transverse bed slope and the interaction between lateral and vertical sorting is not yet understood for lack of data under controlled conditions. The objective of this study is to describe lateral sorting as a function of transverse bed slope and to gain an understanding of the interaction between lateral and vertical sorting in river bends. To this end, experiments were conducted with a poorly sorted sediment mixture in a rotating annular flume in which secondary flow intensity can be controlled separately from the main flow velocity, and therefore transverse bed slope towards the inner bend and dune dimensions can be systematically varied. Sediment samples were taken along cross-sections at the surface of dune troughs and dune crests, and over the entire depth at the location of dune crests (bulk samples), which enabled comparison of the relative contribution of vertical sorting by dunes to lateral sorting by the transverse bed slope. The data show that lateral sorting is always the dominant sorting mechanism in bends, and bulk samples showed minor effects of vertical sorting by dunes as long as all grain-size fractions are mobile. An empirical bend sorting model was fitted that redistributes the available sediment fractions over the cross-section as a function of transverse bed slope. Comparison with field data showed that the model accurately reproduces spatially-averaged trends in sorting at the bend apex in single-thread channels. The bend sorting model therefore provides a better definition of bend sorting with conservation of mass by size fraction and adds to current understanding of bend sorting. The implication for numerical modelling is that bend sorting mechanisms can be modelled independently of dunes, allowing the application of the active layer concept.     

An occurrence of grandidierite,kornerupine, and tourmaline in southeastern Ontario,Canada

Grandidierite, kornerupine, and tourmaline occur in high-grade pelitic gneisses from southeastern Ontario, Canada. The kornerupine occurs in quartz-bearing layers associated with biotite, cordierite, garnet, ilmenite, K-feldspar, magnetite, quartz, and, less commonly, sillimanite. Grandidierite is found in quartz-poor, cordierite+sillimanite segregations in contact with biotite, cordierite, ilmenite, K-feldspar, magnetite, sillimanite, and, more rarely, garnet. Tourmaline is sporadically distributed in all compositional layers, but is not in contact with the other borosilicates. There is no textural evidence for a reactive relationship among the three borosilicates. Neither chemical or textural equilibrium has been achieved on the scale of a thin section.It is proposed that the granite, K-feldspar-rich leucosomes, and different borosilicate assemblages in adjacent compositional layers evolved along a path of decreasing pressure and increasing temperature. The P-T path intersected a series of dehydration and melting reactions. This P-T path indicates that uplift had occurred before cooling had started and before the maximum temperature was reached. Corona and symplectite textures developed at various times during uplift both before and after cooling had started.     

Experimental investigation and application of the equilibrium rutile + orthopyroxene = quartz + ilmenite

Equilibria in the Sirf (Silica-Ilmenite-Rutile-Ferrosilite) system: $${\text{SiO}}_{\text{2}} + ({\text{Mg,Fe}}){\text{TiO}}_{\text{3}} {\text{ + (Mg,Fe)SiO}}_{\text{3}} $$ have been calibrated in the range 800–1100° C and 12–26 kbar using a piston-cylinder apparatus to assess the potential of the equilibria for geobarometry in granulite facies assemblages that lack garnet. Thermodynamic calculations indicate that the two end-member equilibria involving quartz + geikielite = rutile + enstatite, and quartz + ilmenite = rutile + ferrosilite, are metastable. We therefore reversed equilibria over the compositional range Fs_40–70, using Ag₈₀Pd₂₀ capsules with \(f_{{\text{O}}_{\text{2}} } \) buffered at or near iron-wüstite. Ilmenite compositions coexisting with orthopyroxene are \(X_{{\text{MgTiO}}_{\text{3}} }^{{\text{Ilm}}} \) of 0.06 to 0.15 and \(X_{{\text{Fe}}_{\text{2}} {\text{O}}_{\text{3}} }^{{\text{Ilm}}} \) of 0.00 to 0.01, corresponding toK _D values of 13.3, 10.2, 9.0 and 8.0 (±0.5) at 800, 900, 1000 and 1100° C, respectively, whereK _D =(XMg/XFe)^Opx/(XMg/XFe)^Ilm. Pressures have been calculated using equilibria in the Sirf system for granulites from the Grenville Province of Ontario and for granulite facies xenoliths from central Mexico. Pressures are consistent with other well-calibrated geobarometers for orthopyroxeneilmenite pairs from two Mexican samples in which oxide textures appear to represent equilibrium. Geologically unreasonable pressures are obtained, however, where oxide textures are complex. Application of data from this study on the equilibrium distribution of iron and magnesium between ilmenite and orthopyroxene suggests that some ilmenite in deep crustal xenoliths is not equilibrated with coexisting pyroxene, while assemblages from exposed granulite terranes have reequilibrated during retrogression. The Sirf equilibria are sensitive to small changes in composition and may be used for determination of activity/composition (a/X) relations of orthopyroxene if an ilmenite model is specified. A symmetric regular solution model has been used for orthopyroxene in conjunction with activity models for ilmenite available from the literature to calculatea/X relations in orthopyroxene of intermediate composition. Data from this study indicate that FeSiO₃?MgSiO₃ orthopyroxene exhibits small, positive deviations from ideality over the range 800–1100°C.     

Correction: Hydraulic tomography analysis of municipal-well operation data with geology-based groundwater models

Hydrogeology Journal - A Correction to this paper has been published: https://doi.org/10.1007/s10040-021-02351-x     

Relationship among titanium,rare earth elements,U–Pb ages and deformation microstructures in zircon: Implications for Ti-in-zircon thermometry

A zircon grain in an orthopyroxene–garnet–phlogopite–zircon–rutile-bearing xenolith from Udachnaya, Siberia, preserves a pattern of crystallographic misorientation and subgrain microstructure associated with crystal–plastic deformation. The zircon grain records significant variations in titanium (Ti) from 2.6 to 30 ppm that corresponds to a difference in calculated Ti-in-zircon temperatures of over several hundred degrees Celsius. The highest Ti concentration is measured at subgrain centres (30 ppm), and Ti is variably depleted at low-angle boundaries (down to 2.6 ppm). Variations in cathodoluminescence coincide with the deformation microstructure and indicate localised, differential enrichment of rare earth elements (REE) at low-angle boundaries. Variable enrichment of U and Th and systematic increase of Th/U from 1.61 to 3.52 occurs at low-angle boundaries. Individual SHRIMP-derived U–Pb ages from more deformed zones (mean age of 1799 ± 40, n = 22) are systematically younger than subgrain cores (mean age of 1851 ± 65 Ma, n = 7), and indicate that open system behaviour of Ti–Th–U occurred shortly after zircon growth, prior to the accumulation of significant radiogenic Pb. Modelling of trace-element diffusion distances for geologically reasonable thermal histories indicates that the observed variations are ~ 5 orders of magnitude greater than can be accounted for by volume diffusion. The data are best explained by enhanced diffusion of U, Th and Ti along deformation-related fast-diffusion pathways, such as dislocations and low-angle (< 5°) boundaries. These results indicate chemical exchange between zircon and the surrounding matrix and show that Ti-in-zircon thermometry and U–Pb geochronology from deformed zircon may not yield information relating to the conditions and timing of primary crystallisation.     

Flow,stress and sediment resuspension in a shallow tidal channel

In October of 2004, a 3-d observational program to measure flow and sediment resuspension within a coastal intertidal salt marsh was conducted in the North Inlet/Winyah Bay National Estuarine Research Reserve located near Georgetown, South Carolina. Current and acoustic backscatter profiles were obtained from a moored acoustic Doppler current profiler (ADCP) deployed in a shallow tidal channel during the spring phase of the tidal cycle under high discharge conditions. The channel serves as a conduit between Winyah Bay, a large brackish estuary, and North Inlet, a saline intertidal coastal salt marsh with little freshwater input. Salinity measurements indicate that the water column is vertically well mixed during flood, but becomes vertically stratified during early ebb. The stratification results from brackish (15 psu) Winyah Bay water entering North Inlet via the tidal channel, suggesting an exchange mechanism that permits North Inlet to receive a fraction of the poor water quality and high discharge flow from upland rivers. Although maximum flood currents exceed maximum ebb currents by 0.2 m s⁻¹, suspended sediment concentrations are highest during the latter ebb phase and persist for a longer fraction of the ebb cycle. Even though the channel is flood-dominated, the higher concentrations occurring over a longer fraction of the ebb phase indicate net particulate transport from Winyah Bay to North Inlet during spring tide accompanied by high discharge. Our evidence suggests that the higher concentrations during ebb result from increased bed friction caused by flow asymmetries and variations in water depth in which the highest stresses occur near the end of ebb near low water despite stronger maximum currents during flood.     

Fluid ascent through the solid lithosphere and its relation to earthquakes

The Earth is continuously expelling gases and liquids from great depths—juvenile volatiles from the mantle and recycled metamorphic products. Some of these fluids ascend through liquid rock in volcanic processes, but others utilize fractures and faults as conduits through the solid lithosphere. The latter process may have a major influence on earthquakes, since fluids at near lithostatic pressures appear to be required to activate deep faults that would otherwise remain locked.Fluids can be driven upward through solid rock by buoyancy, but only if present in sufficient concentration to form large-scale domains occupying interconnected fracture porosity. A growing fluid domain becomes so mobilized only when it attains the critical vertical dimension required for hydrostatic instability. This dimension, depending on the ultimate compressive yield strength of the rock, may be as much as several kilometers.Any column of fluid ascending through fractures in the solid lithosphere from a prolific deep source must become organized into a vertical sequence of discrete domains, separated by fluid-pressure discontinuities. This is required because a continuous hydrostatic-fluid-pressure profile extending from an arbitrarily deep source to the surface cannot be permitted by the finite strength of rock. A vertically stacked sequence of domains allows the internal fluid-pressure profile to approximate the external rock-stress profile in a stepwise fashion. The pressure discontinuity below the base of the uppermost hydrostatic domain may be responsible for some occurrences of so-called anomalous geopressures. An ascending stream of fluid that percolates upward from a deep source through a column of domains must encounter a sequence of abrupt pressure decreases at the transitions between successive domains. If supercritical gases act as solvents, the dissolved substances may drop out of solution at such pressure discontinuities, resulting in a local concentration of minerals and other substances.At great depths, brittle fracture would normally be prevented by high pressure and temperature, with all excessive stress discharged by ductile flow. Rock strata invaded by an ascending fluid domain are weakened, however, because cracks generated or reactivated by the high-pressure fluid can support the overburden, greatly reducing internal friction. This reduction of strength may cause a previously stressed rock to fail, resulting in hydraulic shear fracture. Thus, earthquakes may be triggered by the buoyant migration of deep-source fluids.The actual timing of the failure that leads to such an earthquake may be determined by the relatively rapid inflation of a fluid domain and not by any significant increase in the probably much slower rate of regional tectonic strain. Many earthquake precursory phenomena may be secondary symptoms of an increase in pore-fluid pressure, and certain coseismic phenomena may result from the venting of high-pressure fluids when faults break the surface. Instabilities in the migration of such fluid domains may also contribute to or cause the eruption of mud volcanoes, magma volcanoes, and kimberlite pipes.     

The mid-ocean ridge axial valley as a steady-state neck

In this paper the mid-ocean ridge axial valley is modelled as a steady-state lithospheric neck in which lithospheric stretching balances lithospheric accretion. Conversely, the axial high is a steady-state lithospheric bulge. The lithosphere is modelled as a thin plate with a Newtonian rheology. It is shown that an axial valley will occur if the rate of viscosity increase away from the ridge axis is faster than the rate at which accretion decreases. An axial high will occur if the opposite condition holds. This is consistent with the observation that axial valleys occur at low spreading rates and axial highs at high spreading rates. By fitting our model to profiles across the Mid-Atlantic Ridge and the East Pacific Rise and assuming the lithospheric thickness at the ridge axis to be 5 km, we find accretion widths of 6–8 km. We find the width over which there is a significant increase in lithospheric viscosity to be also 6–8 km.