Stress‐dependent elastic wave velocity of microfractured sandstone

A model for the stress‐dependent elastic wave velocity response of fractured rock mass is proposed based on experimental evidence of stress‐dependent fracture normal and shear stiffness. Previously proposed models and previous experimental studies on stress‐dependent fracture stiffness have been reviewed to provide a basis for the new model. Most of the existing stress‐dependent elastic wave velocity models are empirical, with model parameters that do not have clear physical meanings. To propose the new model, the rock mass is assumed to have randomly oriented microscopic fractures. In addition, the characteristic length of microfractures is assumed to be sufficiently short compared to the rock mass dimensions. The macroscopic stress‐dependent elastic wave velocity response is assumed to be attributed to the stress dependency of fracture stiffness. The stress‐dependent fracture normal stiffness is defined as a generalized power law function of effective normal stress, which is a modification of the Goodman's model. On the other hand, the stress dependency of fracture shear stiffness is modeled as a linear function of normal stress based on experimental data. Ultrasonic wave velocity responses of a dry core sample of Berea sandstone were tested at effective stresses ranging from 2 to 55 MPa. Visual observation of thin sections obtained from the Berea sandstone confirms that the assumptions made for microstructure of rock mass model are appropriate. It is shown that the model can describe the stress‐dependent ultrasonic wave velocity responses of dry Berea sandstone with a set of reasonable material parameter values. Published 2013. This article is a U.S. Government work and is in the public domain in the USA.     

Modelling spatial variations in dissolved oxygen in fine‐grained streams under uncertainty

This paper presents a novel triple‐layer model, called VART DO‐3L, for simulation of spatial variations in dissolved oxygen (DO) in fine‐grained streams, characterized by a fluid mud (fluff or flocculent) layer (an advection‐dominated storage zone) as the interface between overlying stream water and relatively consolidated streambed sediment (a diffusion‐dominated storage zone). A global sensitivity analysis is conducted to investigate the sensitivity of VART DO‐3L model input parameters. Results of the sensitivity analysis indicate that the most sensitive parameter is the relative size of the advection‐dominated storage zones (A_s/A), followed by a lumped reaction term (R) for the flocculent layer, biological reaction rate (μ_o) in diffusive layer and biochemical oxygen demand concentration (L) in water column. In order to address uncertainty in model input parameters, Monte Carlo simulations are performed to sample parameter values and to produce various parameter combinations or cases. The VART DO‐3L model is applied to the Lower Amite River in Louisiana, USA, to simulate vertical and longitudinal variations in DO under the cases. In terms of longitudinal variation, the DO level decreases from 7.9 mg l at the Denham Springs station to about 2.89 mg l^?1 at the Port Vincent station. In terms of vertical variation, the DO level drops rapidly from the overlying water column to the advection‐dominated storage zone and further to the diffusive layer. The DO level (C_F) in the advective layer (flocculent layer) can reach as high as 40% of DO concentration (C) in the water column. The VART DO‐3L model may be applied to similar rivers for simulation of spatial variations in DO level. Copyright © 2014 John Wiley & Sons, Ltd.     

Soil pipe flow tracer experiments: 1. Connectivity and transport characteristics

Much debate has occurred in catchment hydrology regarding the connectivity of flow paths from upslope areas to catchment outlets. This study was conducted in two catchments, one with three upper branches, in a loess soil with a fragipan that fosters lateral flow and exhibits an extensive distribution of soil pipe collapse features. The study aimed to determine the connectivity of multiple soil pipe networks as well as determine pipe flow velocities during storm events. Fluorescein dye was injected directly into soil pipes at the upper most pipe collapse feature of four different hillslopes. Breakthrough curves (BTC) were determined by sampling multiple pipe collapse features downslope. The BTCs were used to determine the ‘average’ (centre of mass) and ‘maximum’ (first arrival) flow velocities. This study confirmed that these catchments contain individual continuous soil pipe networks that extend over 190 m and connect the upper most hillslopes areas with the catchment outlet. While the flow paths are continuous, the individual pipe networks consist of alternating reaches of subsurface flow through soil pipes and reaches of surface flow through gullies formed by pipe collapses. In addition, flow can be occurring both through the subsurface soil pipes simultaneous with surface flow generated by artesian flow from the soil pipes. The pipe flow velocities were as high as 0.3 m/s, which was in the range of streamflow velocities. These pipe flow velocities were also in the range of velocities observed in pinhole erosion tests suggesting that these large, mature soil pipes are still actively eroding. Copyright © 2015 John Wiley & Sons, Ltd.     

How preferential flow delivers pre‐event groundwater rapidly to streams

Groundwater often accounts for a substantial fraction of flood hydrographs, but the processes responsible for this have been unclear. However, many aquifers have preferential flow and this explains how aquifers can be so responsive. In bedrock aquifers, weathering enhances the connectivity and apertures along the most efficient flow paths and hence enhances the permeability. This results in celerities and velocities of the preferential flow in these dual‐porosity aquifers that are two to three orders of magnitude higher than if the aquifers behaved as single‐porosity media. The celerities have been determined from artificial and natural flood pulses, from tidal lags, and from pumping tests. Preferential‐flow velocities have been calculated from tests using applied tracers. Celerities in bedrock aquifers are typically one to two orders of magnitude faster than velocities. The ubiquitous preferential flow in aquifers provides an additional explanation, besides groundwater ridging, for the rapid release of groundwater to streams during storm events.     

Steady‐state analytical solutions of flow and deformation coupling due to a point sink within a finite fluid‐saturated poroelastic layer

An exact steady‐state closed‐form solution is presented for coupled flow and deformation of an axisymmetric isotropic homogeneous fluid‐saturated poroelastic layer with a finite radius due to a point sink. The hydromechanical behavior of the poroelastic layer is governed by Biot's consolidation theory. Boundary conditions on the lateral surface are specifically chosen to match the appropriate finite Hankel transforms and simplify the transforms of the governing equations. Ordinary differential equations in the transformed domain are solved, and then the analytical solutions in the physical space for the pore pressure and the displacements are finally obtained by using finite Hankel inversions. The analytical solutions at some special locations such as the top and bottom surfaces, lateral surface, and the symmetrical axis are given and analyzed. And a case study for the consolidation of a water‐saturated soft clay layer due to pumping is conducted. The analytical solution is verified against the finite element solution. Meanwhile, an analysis of coupled hydromechanical behavior is carried out herein. The presented analytical solution is an exact solution to the practical poroelastic problem within an axisymmetric finite layer. It can provide us a better understanding of the poroelastic behavior of the finite layer due to fluid extraction. Besides, it can be applied to calibrate numerical schemes of axisymmetric poroelasticity within finite domains. Copyright © 2017 John Wiley & Sons, Ltd.     

Displacement velocity effects on rock fracture shear strengths

Triaxial shear tests are performed to assess the effects of displacement velocity and confining pressure on shear strengths and dilations of tension-induced fractures and smooth saw-cut surfaces prepared in granite, sandstone and marl specimens. A polyaxial load frame is used to apply confining pressures between 1 and 18 MPa with displacement velocities ranging from 1.15 × 10⁻⁵ to 1.15 × 10⁻² mm/s. The results indicate that the shearing resistances of smooth saw-cut surfaces tend to be independent of the displacement velocity and confining pressure. Under each confinement the peak and residual shear strengths and dilation rates of rough fractures increase with displacement velocities. The sheared-off areas increase when the confining pressure increases, and the displacement rate decreases. The velocity-dependent shear strengths tend to act more under high confining pressures for the rough fractures in strong rock (granite) than for the smoother fractures in weaker rocks (sandstone and marl). An empirical criterion that explicitly incorporates the effects of shear velocity is proposed to describe the peak and residual shear strengths. The criterion fits well to the test results for the three tested rocks.     

Effects of surface flux parameterization on the numerically simulated intensity and structure of Typhoon Morakot (2009)

The effects of surface flux parameterizations on tropical cyclone(TC) intensity and structure are investigated using the Advanced Research Weather Research and Forecasting(WRF-ARW) modeling system with high-resolution simulations of Typhoon Morakot(2009).Numerical experiments are designed to simulate Typhoon Morakot(2009) with different formulations of surface exchange coefficients for enthalpy(C_K) and momentum(C_D) transfers,including those from recent observational studies based on in situ aircraft data collected in Atlantic hurricanes.The results show that the simulated intensity and structure are sensitive to C_K and C_D,but the simulated track is not.Consistent with previous studies,the simulated storm intensity is found to be more sensitive to the ratio of C_K/C_D than to C_K or C_D alone.The pressure-wind relationship is also found to be influenced by the exchange coefficients,consistent with recent numerical studies.This paper emphasizes the importance of C_D and C_K on TC structure simulations.The results suggest that C_D and C_K have a large impact on surface wind and flux distributions,boundary layer heights,the warm core,and precipitation.Compared to available observations,the experiment with observed C_D and C_K generally simulated better intensity and structure than the other experiments,especially over the ocean.The reasons for the structural differences among the experiments with different C_D and C_K setups are discussed in the context of TC dynamics and thermodynamics.     

Study of high temperature targets identification and temperature retrieval experimental model in SWIR remote sensing based Landsat8

For surface features in short-wave infrared (SWIR, 1.3–3.0 μm) in remote sensing imagery, pixel values depict the total energy including reflection and emission. For surface features at normal temperature in SWIR band, emission energy can be ignored. While for surface features at high temperature in SWIR band, emission energy is equal to or even higher than the reflection energy. So remote sensing imagery of SWIR band can be used to separate emission and reflection energy as well as to realize temperature retrieval of high temperature targets. In this study, the seventh band (SWIR band) of Landsat8 OLI remote sensing imagery is used to perform the theoretical model research for temperature retrieval of high temperature targets. In the meantime, it is also used with the corresponding observation experiment of synchronization satellite to check the theoretical model. The result shows that the radiant flux density for mixed pixels with high temperature targets is higher than adjacent pixels without high temperature targets. Thus, the high temperature pixels can be identified in SWIR band. The retrieval results of temperature and fractional area for high temperature targets are consistent with reality. In the study, the result illustrates that it is effective to identify high temperature targets in remote sensing imagery of SWIR band and the model is appropriate for temperature retrieval use.     

基于抗差估计的实时单站GPS同震速度获取方法

基于实时单站GPS载波相位历元差分测速模型,引入抗差最小二乘估计,并根据IGGⅢ方案选择合适的等价权因子来削弱小周跳和粗差对结果的影响,并采用一组静态数据和2011年日本"3.11"地震期间MIZU站的高频数据对算法进行测试。结果表明,在含有小周跳和粗差的情况下,抗差最小二乘估计能够明显改善解算速度,可以实时获取测站毫米级的同震速度。     

浅水多波束系统SONIC 2024在码头测深中的应用

文中介绍浅水多波束系统SONIC 2024的技术优点,从姿态参数校正、潮汐改正、声速改正、测线滤波及曲面滤波等方面探讨多波束数据处理的关键技术。以码头实测数据验证该系统在水深测量应用中的便捷性和可靠性。