Research Note: Experimental measurements of Q‐contrast reflections

While seismic reflection amplitudes are generally determined by real acoustical impedance contrasts, there has been recent interest in reflections due to contrasts in seismic‐Q. Herein we compare theoretical and modelled seismic reflection amplitudes for two different cases of material contrasts. In case A, we examine reflections from material interfaces that have a large contrast in real‐valued impedance () with virtually no contrast in seismic‐Q. In case B, we examine reflections from material interfaces that have virtually no contrast in but that have very large seismic‐Q contrasts. The complex‐valued reflection coefficient formula predicts non‐zero seismic reflection amplitudes for both cases. We choose physical materials that typify the physics of both case A and case B. Physical modelling experiments show significantly large reflections for both cases – with the reflections in the two cases being phase shifted with respect to each other, as predicted theoretically. While these modelling experiments show the existence of reflections that are predicted by theory, there are still intriguing questions regarding the size of the Q‐contrast reflections, the existence of large Q‐contrast reflections in reservoir rocks and the possible application of Q‐reflection analysis to viscosity estimation in heavy oilfields.     

New Method for Continuous Transmissivity Profiling in Fractured Rock

A new method is presented to search for hydraulically transmissive features in open boreholes in bedrock. A flexible borehole liner made of a watertight, nylon fabric is filled with water to create a constant driving head to evert (reverse of invert) the liner down the hole so that the liner pushes the borehole water out into transmissive fractures or other permeable features. The descent rate is governed by the bulk transmissivity of the remaining permeable features below the liner. Initially, the liner descent rate or velocity is a measure of transmissivity (T) of the entire hole. As the everting liner passes and seals each permeable feature, changes in the liner velocity indicate the position of each feature and an estimate of T using the Thiem equation for steady radial flow. This method has been performed in boreholes with diameters ranging from 96 to 330 mm. Profiling commonly takes a few hours in holes 200‐ to 300‐m long. After arrival of the liner at the bottom of the hole, the liner acts as a seal preventing borehole cross connection between transmissive features at different depths. Liner removal allows the hole to be used for other purposes. The T values determined using this method in a dolostone aquifer were found to be similar to the values from injection tests using conventional straddle packers. This method is not a replacement for straddle‐packer hydraulic testing of specific zones where greater accuracy is desired; however, it is effective and efficient for scanning entire holes for transmissive features.     

Numerical Modeling of Ice Fog in Interior Alaska Using the Weather Research and Forecasting Model

An ice microphysics parameterization scheme has been modified to better describe and understand ice fog formation. The modeling effort is based on observations in the Sub-Arctic Region of Interior Alaska, where ice fog occurs frequently during the cold season due to abundant water vapor sources and strong inversions existing near the surface at extremely low air temperatures. The microphysical characteristics of ice fog are different from those of other ice clouds, implying that the microphysical processes of ice should be changed in order to generate ice fog particles. Ice fog microphysical characteristics were derived with the NCAR Video Ice Particle Sampler during strong ice fog cases in the vicinity of Fairbanks, Alaska, in January and February 2012. To improve the prediction of ice fog in the Weather Research and Forecasting model, observational data were used to change particle size distribution properties and gravitational settling rates, as well as to implement a homogeneous freezing process. The newly implemented homogeneous freezing process compliments the existing heterogeneous freezing scheme and generates a higher number concentration of ice crystals than the original Thompson scheme. The size distribution of ice crystals is changed into a Gamma distribution with the shape factor of 2.0, using the observed size distribution. Furthermore, gravitational settling rates are reduced for the ice crystals since the crystals in ice fog do not precipitate in a similar manner when compared to the ice crystals of cirrus clouds. The slow terminal velocity plays a role in increasing the time scale for the ice crystals to settle to the surface. Sensitivity tests contribute to understanding the effects of water vapor emissions as an anthropogenic source on the formation of ice fog.     

Impact of reach geometry on stream channel sensitivity to extreme floods

Predicting spatial and temporal variations in bank erosion due to extreme floods presents a long‐standing challenge in geomorphology. We develop two methodologies for rapid, regional‐scale assessments of stream reaches susceptible to channel widening. The first proposes that channel widening occurs when unit stream power exceeds a critical threshold (300 W/m²). The second is motivated by the observation that widening often occurs at channel bends. We introduce a new metric, the bend stress parameter, which is proportional to the centripetal force exerted on a concave bank. We propose that high centripetal forces generate locally high bank shear forces and enhance channel bank erosion. We test both metrics using the geomorphic signature of Tropical Storm Irene (2011) on the White and the Saxtons Rivers, Vermont. Specifically, we test if reaches where significant channel widening occurred during Irene required one or both metrics to exceed threshold values. We observe two distinct styles of channel widening. Where unit stream power and bend stress parameter are high, widening is usually due to bank retreat. Elsewhere widening is usually due to the stripping of the upstream end of mid‐channel islands. Excluding widening associated with the stripping of the heads of mid‐channel islands, almost all the widening (> 98%) occurred along reaches identified as susceptible to widening. The combined metrics identify up to one‐quarter of the reaches lacking susceptibility to channel widening. Copyright © 2014 John Wiley & Sons, Ltd.     

Remote measurement of river morphology via fusion of LiDAR topography and spectrally based bathymetry

This study developed and evaluated a hybrid approach to remote measurement of river morphology that combines LiDAR topography with spectrally based bathymetry. Comparison of filtered LiDAR point clouds with surveyed cross‐sections indicated that subtle features on low‐relief floodplains were accurately resolved by LiDAR but that submerged areas could not be detected due to strong absorption of near‐infrared laser pulses by water. The reduced number of returns made the active channel evident in a LiDAR point density map. A second dataset suggested that pulse intensity also could be used to discriminate land from water via a threshold‐based masking procedure. Fusion of LiDAR and optical data required accurate co‐registration of images to the LiDAR, and we developed an object‐oriented procedure for achieving this alignment. Information on flow depths was derived by correlating pixel values with field measurements of depth. Highly turbid conditions dictated a positive relation between green band radiance and flow depth and contributed to under‐prediction of pool depths. Water surface elevations extracted from the LiDAR along the water's edge were used to produce a continuous water surface that preserved along‐channel variations in slope. Subtracting local flow depths from this surface yielded estimates of the bed elevation that were then combined with LiDAR topography for exposed areas to create a composite representation of the riverine terrain. The accuracy of this terrain model was assessed via comparison with detailed field surveys. A map of elevation residuals showed that the greatest errors were associated with underestimation of pool depths and failure to capture cross‐stream differences in water surface elevation. Nevertheless, fusion of LiDAR and passive optical image data provided an efficient means of characterizing river morphology that would not have been possible if either dataset had been used in isolation. Copyright © 2011 John Wiley & Sons, Ltd.     

Elastic source model of the North Mono eruption (1325–1368 A.D.) based on shoreline deformation

Topographic data from the Shuttle Radar Topography Mission (SRTM) captures the permanent deformation of a prominent highstand of Mono Lake, California USA. Deformation of the Dechambeau Ranch highstand shoreline was measured using the elevation of the beach berm—shoreline bluff break-in-slope. Point source models and a boundary element dike model were used to approximate the source of deformation underneath the northern end of the Mono Craters. The point source model could not adequately explain the observed deformation. The dike model yielded the best results for a NW trending dike dipping 60° NE and inflated to widths greater than 60 m. The results suggest that the geometry of the source is more complex than a simple vertical dike and that the deformation is better explained with a dipping dike following a normal fault, or an elongated cryptodome.     

The role of biophysical interactions within the ijzermonding tidal flat sediment dynamics

This paper focuses on the importance of biophysical interactions on short-term and long-term sediment dynamics. Therefore, various biological (macrobenthos, photopigments, colloidal EPS) and physical parameters (grain size, water content, sediment stability, bed level) were determined (bi)monthly in nine sampling plots on the IJzermonding tidal flat (Belgium, 51°08′N, 2°44′E) during three consecutive years (July 2005–June 2008). Results showed that sediment stability varied on the short timescale and was directly influenced by biota, while bed level varied mainly on the long-term due to interannual variability. The short-term dynamic relationships between mud content, water content, fucoxanthin and macrobenthos density resulted in a seasonal mud deposition and erosion cycle, and directly influenced sediment stability. Moreover, macrobenthos was proven to be the most important parameter determining sediment stability. On the long-term, a shift was observed from high fucoxanthin/chla concentration, high mud content and zero to moderate densities of Corophium volutator towards low fucoxanthin/chl a and mud content and high Corophium densities, which resulted in a transition from net accretion to net erosion. However, most measured variables proved to be poor predictors for these long-term bed level changes, indicating that external physical forces, such as waves and storminess, probably were the most important factors triggering long-term sediment dynamics. Nevertheless, biota indirectly influenced bed level changes by mediating short-term changes in sediment stability, thereby influencing the erodability of the sediment. The macrobenthos, and especially the mud shrimp Corophium, was suggested as the (indirect) driving destabilising factor for the sampling plots in the IIzermonding when considering the long-term evolution.     

The erosion rates of cohesive sediments in Venice lagoon,Italy

The stability of cohesive sediments from Venice lagoon has been measured in situ using the benthic flume Sea Carousel. Twenty four stations were occupied during summertime, and a sub-set of 13 stations was re-occupied during the following winter. Erosion thresholds and first-order erosion rates were estimated and showed a distinct difference between inter-tidal and sub-tidal stations. The higher values for inter-tidal stations are the result of exposure that influences consolidation, density, and organic adhesion. The thresholds for each state of sediment motion are well established. However, the rate of erosion once the erosion threshold has been exceeded has been poorly treated. This is because normally a time-series of sediment concentration (C) and bed shear stress (τ₀(t)) is used to define threshold stress or cohesion (τ_crit,z) and erosion rate (E). Whilst solution of the onset of erosion, τ_crit,0, is often reported, the evaluation of the erosion threshold variation through the process of erosion (eroded depth) is usually omitted or not estimated. This usually leads to assumptions on the strength profile of the bed which invariably has no credibility within the topmost mm of the bed where most erosion takes place. It is possible to extract this information from a time-series through the addition of a step in data processing. This paper describes how this is done, and the impact of this on the accuracy of estimates of the excess stress (τ₀(t)–τ_crit,z) on E.     

Probabilistic data integration for characterization of spatial distribution of residual LNAPL

The spatial distribution of residual light non-aqueous phase liquid (LNAPL) is an important factor in reactive solute transport modeling studies. There is great uncertainty associated with both the areal limits of LNAPL source zones and smaller scale variability within the areal limits. A statistical approach is proposed to construct a probabilistic model for the spatial distribution of residual NAPL and it is applied to a site characterized by ultra-violet-induced-cone-penetration testing (CPT–UVIF). The uncertainty in areal limits is explicitly addressed by a novel distance function (DF) approach. In modeling the small-scale variability within the areal limits, the CPT–UVIF data are used as primary source of information, while soil texture and distance to water table are treated as secondary data. Two widely used geostatistical techniques are applied for the data integration, namely sequential indicator simulation with locally varying means (SIS–LVM) and Bayesian updating (BU). A close match between the calibrated uncertainty band (UB) and the target probabilities shows the performance of the proposed DF technique in characterization of uncertainty in the areal limits. A cross-validation study also shows that the integration of the secondary data sources substantially improves the prediction of contaminated and uncontaminated locations and that the SIS–LVM algorithm gives a more accurate prediction of residual NAPL contamination. The proposed DF approach is useful in modeling the areal limits of the non-stationary continuous or categorical random variables, and in providing a prior probability map for source zone sizes to be used in Monte Carlo simulations of contaminant transport or Monte Carlo type inverse modeling studies.     

Statistic inversion of multi-zone transition probability models for aquifer characterization in alluvial fans

Understanding the heterogeneity arising from the complex architecture of sedimentary sequences in alluvial fans is challenging. This paper develops a statistical inverse framework in a multi-zone transition probability approach for characterizing the heterogeneity in alluvial fans. An analytical solution of the transition probability matrix is used to define the statistical relationships among different hydrofacies and their mean lengths, integral scales, and volumetric proportions. A statistical inversion is conducted to identify the multi-zone transition probability models and estimate the optimal statistical parameters using the modified Gauss–Newton–Levenberg–Marquardt method. The Jacobian matrix is computed by the sensitivity equation method, which results in an accurate inverse solution with quantification of parameter uncertainty. We use the Chaobai River alluvial fan in the Beijing Plain, China, as an example for elucidating the methodology of alluvial fan characterization. The alluvial fan is divided into three sediment zones. In each zone, the explicit mathematical formulations of the transition probability models are constructed with optimized different integral scales and volumetric proportions. The hydrofacies distributions in the three zones are simulated sequentially by the multi-zone transition probability-based indicator simulations. The result of this study provides the heterogeneous structure of the alluvial fan for further study of flow and transport simulations.