Land cover effects on runoff patterns in eastern Piedmont (USA) watersheds

Physiography and land cover determine the hydrologic response of watersheds to climatic events. However, vast differences in climate regimes and variation of landscape attributes among watersheds (including size) have prevented the establishment of general relationships between land cover and runoff patterns across broad scales. This paper addresses these difficulties by using power spectral analysis to characterize area‐normalized runoff patterns and then compare these patterns with landscape features among watersheds within the same physiographic region. We assembled long‐term precipitation and runoff data for 87 watersheds (first to seventh order) within the eastern Piedmont (USA) that contained a wide variety of land cover types, collected environmental data for each watershed, and compared the datasets using a variety of statistical measures. The effect of land cover on runoff patterns was confirmed. Urban‐dominated watersheds were flashier and had less hydrologic memory compared with forest‐dominated watersheds, whereas watersheds with high wetland coverage had greater hydrologic memory. We also detected a 10–15% urban threshold above which urban coverage became the dominant control on runoff patterns. When spectral properties of runoff were compared across stream orders, a threshold after the third order was detected at which watershed processes became dominant over precipitation regime in determining runoff patterns. Finally, we present a matrix that characterizes the hydrologic signatures of rivers based on precipitation versus landscape effects and low‐frequency versus high‐frequency events. The concepts and methods presented can be generally applied to all river systems to characterize multiscale patterns of watershed runoff. Copyright © 2013 John Wiley & Sons, Ltd.     

Spatially distributed long‐term hydrologic simulation using a continuous SCS CN method‐based hybrid hydrologic model

A continuous Soil Conservation Service (SCS) curve number (CN) method that considers time‐varied SCS CN values was developed based on the original SCS CN method with a revised soil moisture accounting approach to estimate run‐off depth for long‐term discontinuous storm events. The method was applied to spatially distributed long‐term hydrologic simulation of rainfall‐run‐off flow with an underlying assumption for its spatial variability using a geographic information systems‐based spatially distributed Clark's unit hydrograph method (Distributed‐Clark; hybrid hydrologic model), which is a simple few parameter run‐off routing method for input of spatiotemporally varied run‐off depth, incorporating conditional unit hydrograph adoption for different run‐off precipitation depth‐based direct run‐off flow convolution. Case studies of spatially distributed long‐term (total of 6 years) hydrologic simulation for four river basins using daily NEXRAD quantitative precipitation estimations demonstrate overall performances of Nash–Sutcliffe efficiency (E_NS) 0.62, coefficient of determination (R²) 0.64, and percent bias 0.33% in direct run‐off and E_NS 0.71, R² 0.72, and percent bias 0.15% in total streamflow for model result comparison against observed streamflow. These results show better fit (improvement in E_NS of 42.0% and R² of 33.3% for total streamflow) than the same model using spatially averaged gauged rainfall. Incorporation of logic for conditional initial abstraction in a continuous SCS CN method, which can accommodate initial run‐off loss amounts based on previous rainfall, slightly enhances model simulation performance; both E_NS and R² increased by 1.4% for total streamflow in a 4‐year calibration period. A continuous SCS CN method‐based hybrid hydrologic model presented in this study is, therefore, potentially significant to improved implementation of long‐term hydrologic applications for spatially distributed rainfall‐run‐off generation and routing, as a relatively simple hydrologic modelling approach for the use of more reliable gridded types of quantitative precipitation estimations.     

西洞庭湖沉水植物分布格局对环境因子及水文情势差异的响应

沉水植物是湖泊生态系统的重要组成部分,其生产力和分布格局受环境因子特别是水文情势决定.洞庭湖是长江流域重要的大型通江湖泊,近年来受人为干扰和气候变化影响,水文节律与水质等环境因子发生改变,导致沉水植物出现衰退现象,急需开展科学恢复,因此有必要对洞庭湖沉水植物深入研究.本研究选取西洞庭湖为研究区域,于2018年和2019年夏季调查了12处典型生境、98个样点的沉水植物与水深、透明度等11个环境因子,采用独立样本T检验和冗余分析方法对沉水植物与环境因子的关系进行分析,对比年际水文情势变化的影响.结果表明：1）西洞庭湖沉水植物主要在水深较浅、水质更优、水体更为稳定的半阻隔子湖和自由连通的湖湾区分布,在河道及水位波动较大的区域分布较少,有、无沉水植物分布样点间存在显著差异的环境因子为水深、透明度、底泥总磷和pH;2）在有沉水植物分布的样点,沉水植物生物量与pH、水深和水体总磷呈显著相关关系;3）自然连通的季节性淹没湖泊沉水植物生物量在2018年高于2019年,可能与2019年5—8月沉水植物关键生长期出现的涨水过程有关,持续的高水位对沉水植物的生长产生了不利影响.维持自然水文节律、湖泊生境异质性与自由连通性、健康的水质等是恢复西洞庭沉水植物的关键,建议在水深低于3 m、营养盐浓度适中、流速及风浪较小的湖湾区或半阻隔湖泊开展沉水植物恢复.     

Impact of climate change on runoff from a mid‐latitude mountainous catchment in central Japan

Hydrologic balance in high‐altitude, mid‐latitude mountain areas is important in terms of the water resources available to associated lowlands. This study examined how current and historical shifts in precipitation (P) patterns and concurrent increases in temperature (T) affected runoff (Q) and other hydrologic components in a mid‐latitude mountain catchment of central Japan, using a combination of long‐term data and a simplified hydrologic model, along with their stochastic treatment. The availability of intensive meteorological and hydrological data from the period 1997–2001 allowed the derivation of key relationships for the current climate that tie the forcing term to the parameters or state variables. By using the data recorded in the period 1965–2001, the force for driving the historical simulation was generated. Based on this model and historical shifts in P and T, the probability density functions of Q (pdf(Q)) was computed. A main novelty in this study is that such a stochastic representation, which is useful for considering the influence of projected shifts in environmental factors on the hydrologic budget, was provided. Despite the large increase in the rate of T in winter and spring, pdf(Q) in spring and summer varied appreciably during the time studied mainly because of an increase in snowmelt. An interannual change in whole‐year Q was robust to shifts in T because while Q in spring increased, in summer it decreased, implying a crucial effect of global warming on mountain hydrologic regimes is change in the timing of Q. Copyright © 2009 John Wiley & Sons, Ltd.     

Fractal techniques and the surface roughness of talus slopes

Linear plots of log N against log G, where N is the number of steps of length G to span a transit, are conventionally used as evidence that geomorphic surfaces are self-similar fractals (i.e. the surfaces have a constant fractal dimension). In this study 42 transits on talus slope surfaces in Niagara and Letchworth Gorges, western New York, are investigated to ascertain whether they are self-similar. Log N-log G plots, which r² values in excess of 0·99 suggest are linear, are found upon more rigorous testing to be curvilinear. It is concluded that the talus slope surfaces are not self-similar, and that log N-log G plots are relatively insensitive to departures from self-similarity. The curvilinearity of the log N-log G plots is explained with the aid of a randomized square-wave model of the talus slope surfaces. This model is used to extend the range of measurement beyond that which was possible in the empirical analysis. The negative of the gradient of the log N -log G relation at a point is the fractal dimension D. Measurements made upon the randomized square-wave model indicate that the relation between D and scale of measurement has an asymmetrical wave shape with a peak (i.e. maximum D) where the scale of measurement is equal to the characteristic scale of roughness. In other words the value of D for a surface is not absolute but depends on the scale of measurement relative to the scale of roughness. Linear regression analysis reveals that at the scale of measurement employed in this study, D is positively correlated with particle size. This is because the values of D fall on the right-hand tail of the wave-shaped relation between D and scale of measurement. Transects (normal to the direction of slope) are found to have higher values of D than profiles (parallel to the direction of slope), and this is explained in terms of particle orientation, shape, and juxtaposition. Because D varies continuously with scale of measurement, there are considerable difficulties in using it to characterize and compare the surface roughness of talus slopes. Generalizing from talus slopes to other ground surfaces, it is evident that to the extent that any natural ground surface has a characteristic scale of roughness, it will depart from self-similarity, and D should be used with caution in quantifying the roughness of the surface. Geomorphologists are therefore urged to be more rigorous in their testing of self-similarity before employing D to characterize surface roughness.     

First‐order exchange coefficient coupling for simulating surface water–groundwater interactions: parameter sensitivity and consistency with a physics‐based approach

Distributed hydrologic models capable of simulating fully‐coupled surface water and groundwater flow are increasingly used to examine problems in the hydrologic sciences. Several techniques are currently available to couple the surface and subsurface; the two most frequently employed approaches are first‐order exchange coefficients (a.k.a., the surface conductance method) and enforced continuity of pressure and flux at the surface‐subsurface boundary condition. The effort reported here examines the parameter sensitivity of simulated hydrologic response for the first‐order exchange coefficients at a well‐characterized field site using the fully coupled Integrated Hydrology Model (InHM). This investigation demonstrates that the first‐order exchange coefficients can be selected such that the simulated hydrologic response is insensitive to the parameter choice, while simulation time is considerably reduced. Alternatively, the ability to choose a first‐order exchange coefficient that intentionally decouples the surface and subsurface facilitates concept‐development simulations to examine real‐world situations where the surface‐subsurface exchange is impaired. While the parameters comprising the first‐order exchange coefficient cannot be directly estimated or measured, the insensitivity of the simulated flow system to these parameters (when chosen appropriately) combined with the ability to mimic actual physical processes suggests that the first‐order exchange coefficient approach can be consistent with a physics‐based framework. Copyright © 2009 John Wiley & Sons, Ltd.     

Topological invariants of field lines rooted to planes

Abstract

Two open curves with fixed endpoints on a boundary surface can be topologically linked. However, the Gauss linkage integral applies only to closed curves and cannot measure their linkage. Here we employ the concept of relative helicity in order to define a linkage for open curves. For a magnetic field consisting of closed field lines, the magnetic helicity integral can be expressed as the sum of Gauss linkage integrals over pairs of lines. Relative helicity extends the helicity integral to volumes where field lines may cross the boundary surface. By analogy, linkages can be defined for open lines by requiring that their sum equal the relative helicity.

With this definition, the linkage of two lines which extend between two parallel planes simply equals the number of turns the lines take about each other. We obtain this result by first defining a gauge-invariant, one-dimensional helicity density, i.e. the relative helicity of an infinitesimally thin plane slab. This quantity has a physical interpretation in terms of the rate at which field lines lines wind about each other in the direction normal to the plane. A different method is employed for lines with both endpoints on one plane; this method expresses linkages in terms of a certain Gauss linkage integral plus a correction term. In general, the linkage number of two curves can be put in the form L=r + n, |r|≦1J2, where r depends only on the positions of the endpoints, and n is an integer which reflects the order of braiding of the curves.

Given fixed endpoints, the linkage numbers of a magnetic field are ideal magneto-hydrodynamic invariants. These numbers may be useful in the analysis of magnetic structures not bounded by magnetic surfaces, for example solar coronal fields rooted in the photosphere. Unfortunately, the set of linkage numbers for a field does not uniquely determine the field line topology. We briefly discuss the problem of providing a complete and economical classification of field topologies, using concepts from the theory of braid equivalence classes.     

Soil hydraulic conductivities of latosols under pasture,forest and teak in Rondonia,Brazil

We investigated the changes of saturated hydraulic conductivity, K_sat , with depth of latosols developed on Precambrian basement rocks under primary rainforest, pasture and teak. In all cases, K_sat decreased with depth, with most of the decrease occurring between the surface and a depth of 30 cm. In conjunction with prevailing rainfall intensities and frequencies, this anisotropy supports a pronounced lateral component of hillslope flow paths, and also of overland flow under pasture. Our results are at variance with data from other latosols where K_sat tends to increase with depth, and hence suggest that considerable restraint is needed in generalization and extrapolation until results from a co‐ordinated effort at hydrology‐oriented data collection become available. Copyright © 1999 John Wiley & Sons, Ltd.     

Jet instability and the stabilizing effect of topography on jets in two-layer rotating systems

Abstract

Laboratory experiments concerning azimuthal jets in two-layer rotating systems in the absence and presence of bottom topography aligned along the jets have been conducted. The jets were forced by the selective withdrawal of fluid from the upper layer of a two-fluid system contained in a circular dishpan geometry. The principal parameters measured in the experiments were the jet Rossby number, Ro, and a stratification parameter F = r ₁/(λ₁λ₂)^1/2 where r ₁ is the radius of the circular disc used for the selective withdrawal (i.e., r ₁ is the approximate radius of curvature of the jet) and λ₁,λ₂ are the internal Rossby radii of deformation in the upper and lower fluids, respectively.

The no-topography experiments show that for a sufficiently small F, the particular value depending on Ro, the jet is stable for the duration of the experiment. For sufficiently large F, again as a function of Ro, the jet becomes unstable, exhibiting horizontal wave disturbances from modes three to seven. An Ro against F flow regime diagram is presented.

Experiments are then conducted in the presence of a bottom topography having constant cross-section and extending around a mid-radius of the dishpan. The axis of the topography is in the vicinity of the jet axis forced in the no-topography experiments and the crest of the topography is in the vicinity of the interface between the two fluids (i.e., the front associated with the jet). The experiments show that in all cases investigated the jet tends to be stabilized by the bottom topography. Experiments with the topography in place, but with the interface between the fluids being above the topography crest, are shown to be unstable but more irregular than their no-topography counterparts.

Various quantitative measurements of the jet are presented. It is shown, for example, that the jet Rossby number defined in terms of the fluid withdrawal rate from the tank. Q, can be well correlated with a dimensionless vorticity gradient, V_G , across the upper layer jet. This allows for an assessment of the stability characteristics of a jet based on a knowledge of V_G (which can be estimated given a jet profile) and F.     

On the Behavior of the Stable Boundary Layer and the Role of Initial Conditions

Previous studies of the stable atmospheric boundary layer using techniques of nonlinear dynamical systems (MCNIDER et al., 1995) have shown that the equations support multiple solutions in certain parameter spaces. When geostrophic speed is used as a bifurcation parameter, two stable equilibria are found—a warm solution corresponding to the high-wind regime where the surface layer of the atmosphere stays coupled to the outer layer, and a cold solution corresponding to the low-wind, decoupled case. Between the stable equilibria is an unstable region where multiple solutions exist. The bifurcation diagram is a classic S shape with the foldback region showing the multiple solutions. These studies were carried out using a simple two-layer model of the atmosphere with a fairly complete surface energy budget. This allowed the dynamical analysis to be carried out on a coupled set of four ordinary differential equations. The present paper extends this work by examining additional bifurcation parameters and, more importantly, analyzing a set of partial differential equations with full vertical dependence. Simple mathematical representations of classical problems in dynamical analysis often exhibit interesting behavior, such as multiple solutions, that is not retained in the behavior of more complete representations. In the present case the S-shaped bifurcation diagram remains with only slight variations from the two-layer model. For the parameter space in the foldback region, the evolution of the boundary layer may be dramatically affected by the initial conditions at sunset. An eigenvalue analysis carried out to determine whether the system might support pure limit-cycle behavior showed that purely complex eigenvalues are not found. Thus, any cyclic behavior is likely to be transient.     

Experimental study on transient and steady‐state dynamics of bedforms in supply limited configuration

The purpose of the present study is to investigate experimentally the development of bedforms in a configuration where the sediment supply is limited. The experimental setup is a rectangular closed duct combining an innovative system to control the rate of sediment supply Q_in , and a digitizing system to measure in real time the 3D bedform topography. We carried out different sets of experiments with two sediment sizes (100 µm and 500 µm) varying both the sediment supply and the water flow rate to obtain a total of 46 different configurations. After a transient phase, steady sub‐centimeter bedforms of various shapes have been observed: barchans dunes, straight transverse dune, linguoid transverse dunes and bedload sheets. Height, spacing, migration speed, and mean bed elevation of the equilibrium bedforms were measured. For a given flow rate, two regimes were identified with fine sediment: (i) a monotonic increasing regime where the equilibrium bedform height and velocity increase with the sediment supply rate Q_in and (ii) an invariant regime for which both parameters are almost independent of Q_in. For coarse sediment, only the first regime is observed. We interpret the saturation of height and velocity for fine sediment bedforms as the transition from a supply‐limited regime to a transport‐limited regime in which the bedload flux has reached its maximum value under the prevailing flow conditions. We also demonstrate that all experiments can be rescaled if the migration speed and height of the bedforms are, respectively, divided and multiplied by the cube of the shear velocity. This normalization is independent of grain size and of bedform morphology. These experimental results provide a new quantification of the factors controlling equilibrium height and migration speed of bedforms in supply‐limited conditions against which theoretical and numerical models can be tested.     

The influence of flow discharge variations on the morphodynamics of a diffluence–confluence unit on a large river

Bifurcations are key geomorphological nodes in anabranching and braided fluvial channels, controlling local bed morphology, the routing of sediment and water, and ultimately defining the stability of their associated diffluence–confluence unit. Recently, numerical modelling of bifurcations has focused on the relationship between flow conditions and the partitioning of sediment between the bifurcate channels. Herein, we report on field observations spanning September 2013 to July 2014 of the three‐dimensional flow structure, bed morphological change and partitioning of both flow discharge and suspended sediment through a large diffluence–confluence unit on the Mekong River, Cambodia, across a range of flow stages (from 13 500 to 27 000 m³ s^?1). Analysis of discharge and sediment load throughout the diffluence–confluence unit reveals that during the highest flows (Q = 27 000 m³ s^?1), the downstream island complex is a net sink of sediment (losing 2600 ± 2000 kg s^?1 between the diffluence and confluence), whereas during the rising limb (Q = 19 500 m³ s^?1) and falling limb flows (Q = 13 500 m³ s^?1) the sediment balance is in quasi‐equilibrium. We show that the discharge asymmetry of the bifurcation varies with discharge and highlight that the influence of upstream curvature‐induced water surface slope and bed morphological change may be first‐order controls on bifurcation configuration. Comparison of our field data to existing bifurcation stability diagrams reveals that during lower (rising and falling limb) flow the bifurcation may be classified as unstable, yet transitions to a stable condition at high flows. However, over the long term (1959–2013) aerial imagery reveals the diffluence–confluence unit to be fairly stable. We propose, therefore, that the long‐term stability of the bifurcation, as well as the larger channel planform and morphology of the diffluence–confluence unit, may be controlled by the dominant sediment transport regime of the system. © 2017 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Microtopographic evolution of mineral surfaces as a tool to identify and date young fault scarps in bedrock

Faulting that results in surface ruptures through bedrock can be particularly difficult to date. For example, stratigraphic control on the age of faulting, based on the age of the bedrock, often leaves unacceptably large uncertainty on the age of the faulting. From a paleoseismological perspective, there is a clear need to determine if a bedrock fault scarp is actually a young feature. For young fault ruptures that create fresh mineral surfaces, analysis of microtopography developed by weathering of the mineral surface may provide a quantifiable method for determining the fault age. The direct quantitative measurement of mineral surface microtopography using Atomic Force Microscopy affords a novel method to study the rupture ages of active faults. The method for using microtopographic evolution of mineral surfaces depends on three conditions. The first condition is that freshly exposed mineral cleavage surfaces, which can be described geometrically as planes, are formed during a rupture event. The formation of these fresh surfaces is analogous to the initiation of a weathering ‘clock’ that defines time t=0. Following cleavage formation dissolution of the planar mineral surface occurs. The rate of dissolution for a mineral species under given climatic conditions, governs the rate of mineral surface alteration. Thus as dissolution proceeds, the roughness of the mineral surface increases. We suggest that the progression of microtopographic roughness over time, which can be estimated by computing quantitative statistics derived from digital mineral surface topography, will systematically vary until a steady state surface topography is reached. The fractal dimension, D_f, is one such measure of surface roughness where, D_f at time t=0 is 2. The dissolution of the mineral surface increases the fractal dimension as the removal of material proceeds. We posit that somewhere between D_f=2 and D_f=3, the microtopography reaches a steady state. Therefore, in the pre-steady state stage of surface roughness, the quantitative measure of roughness of the mineral may serve as a measure of time elapsed since faulting. The period of time this initial stage of surface roughening represents is dependent on the mineral and as a consequence, its dissolution rate, in a specific set of environmental conditions. The time elapsed since fault rupture and grain cleavage can also be estimated from the measurement of the volume of material removed through dissolution. If part of the original cleavage surface remains and can be identified then AFM measurements of the surface microtopography can be used to calculate the dissolved volume per unit area.     

Understanding and modeling basin hydrology: interpreting the hydrogeological signature

Basin landscapes possess an identifiable spatial structure, fashioned by climate, geology and land use, that affects their hydrologic response. This structure defines a basin's hydrogeological signature and corresponding patterns of runoff and stream chemistry. Interpreting this signature expresses a fundamental understanding of basin hydrology in terms of the dominant hydrologic components: surface, interflow and groundwater runoff. Using spatial analysis techniques, spatially distributed watershed characteristics and measurements of rainfall and runoff, we present an approach for modelling basin hydrology that integrates hydrogeological interpretation and hydrologic response unit concepts, applicable to both new and existing rainfall‐runoff models. The benefits of our modelling approach are a clearly defined distribution of dominant runoff form and behaviour, which is useful for interpreting functions of runoff in the recruitment and transport of sediment and other contaminants, and limited over‐parameterization. Our methods are illustrated in a case study focused on four watersheds (24 to 50 km²) draining the southern coast of California for the period October 1988 though to September 2002. Based on our hydrogeological interpretation, we present a new rainfall‐runoff model developed to simulate both surface and subsurface runoff, where surface runoff is from either urban or rural surfaces and subsurface runoff is either interflow from steep shallow soils or groundwater from bedrock and coarse‐textured fan deposits. Our assertions and model results are supported using streamflow data from seven US Geological Survey stream gauges and measured stream silica concentrations from two Santa Barbara Channel–Long Term Ecological Research Project sampling sites. Copyright © 2004 John Wiley & Sons, Ltd.     

Nonlinear modal analysis of penetrative convection

Abstract

The modal expansion procedure has been used to analyze penetrative convection that arises when a thin unstable layer is embedded between two stable regions. The Boussinesq approximation is applied in which the effect of compressibility and stratification are neglected. Various calculations have been made, with one and two modes, for Rayleigh numbers ranging from the critical value to more than 10⁵ times critical. The effect of decreasing the Prandtl number has also been investigated.

It is found that in the nonlinear regime, the convective motions penetrate substantially into the stable regions. The flux of kinetic energy plays a crucial role in such penetration, and its existence puts some requirements on the motions: in the single-mode case, they need to be three-dimensional. The extent of penetration amounts to about half of the thickness of the unstable layer on each side of it when the degree of instability and that of stability are comparable in the two domains; it increases as the stability of the outer region is lowered. The penetration depth appears to be independent of all other parameters defining the problem.     

Exact Elastodynamic Green Functions for Simple Types of Anisotropy Derived from Higher-Order Ray Theory

Using higher-order ray theory, we derived exact elastodynamic Green functions for three simple types of homogeneous anisotropy. The first type displays an orthorhombic symmetry, the other two types display transverse isotropy. In all cases, the slowness surfaces of waves are either ellipsoids, spheroids or spheres. All three Green functions are expressed by a ray series with a finite number of terms. The Green functions can be written in explicit and elementary form similar to the Stokes solution for isotropy. In two Green functions, the higher-order ray approximations form a near-singularity term, which is significant near a kiss singularity. In the third Green function, the higher-order ray approximations also form a near-field term, which is significant near the point source. No effect connected with the line singularity was observed.     

Simulation and statistical modelling approaches to investigate hydrologic regime transformations following Eastern hemlock decline

Ecohydrological processes occurring at or near the Earth's surface are strongly influenced by Eastern hemlock [EH; Tsuga canadensis (L.) Carrière], a foundation tree species of eastern North American forests. EH populations are currently threatened by the invasive hemlock woolly adelgid (HWA; Adelges tsugae Annand). HWA HWA populations have been expanding rapidly throughout the EH's range. Catchment-scale research examining the hydrological consequences of HWA infestation is lacking, and plot-scale studies remain conflicted in their findings. Given the complex relationships between canopy interception, unsaturated and saturated groundwater storage, and root water uptake, it is not immediately clear how EH loss will affect the hydrologic cycle. We investigated the impact of EH mortality on stream discharge characteristics across a regional sample of catchments utilizing both simulation and statistical modelling approaches. We first examined the relationship between various catchment characteristics, including EH health, and three hydrological variables through regression analysis. We then employed a non-parametric statistical test to evaluate differences in hydrologic regime trends between non-infested and infested catchments. Finally, we calibrated a physically based hydrologic model and considered differences in optimal model parameter values and simulated overland runoff between non-infested and infested catchments. HWA presence modified several ecohydrological characteristics and precipitation partitioning between groundwater flows and surface runoff, potentially driving higher stream flashiness and overland flow, lower baseflow contributions and catchment storage, shorter flow-path lengths, and variable source area dilation at infested sites. Our results suggest that EH decline will augment flooding potential associated with the increasing frequency and intensity of Atlantic Basin tropical cyclone events. Further, our physically based simulation provides more determinate results than regression analysis, indicating that a purely statistical methodology, commonly utilized in studying the relationship between landcover characteristics and hydrologic regime, neglects dynamic physical ecohydrologic relationships.     

Spatial variability of soils developing on basalt flows in the Potrillo volcanic field,southern New Mexico: prelude to a chronosequence study

If the systematic spatial variability of soils in a chronosequence is identified and accounted for, the accuracy of quantitative data derived from soil chronosequence studies will be increased. A sample design using landscape positions with minimal variability could result in more accurate chronofunctions from these studies. Four basalt flows in the Potrillo volcanic field, southern New Mexico, with ages ranging between 20 ka and 260 ka (⁴⁰Ar/³⁹Ar and/or cosmogenic ³He methods) provide a sound basis for a soil chronosequence study. Basalt flow surface relief in the Potrillos reduces with time as depressions fill with basalt rubble and aeolian dust. Soil variability is primarily a function of landscape position with respect to ridges and swales in the original basalt flow topography. Soils developing over original topographic lows (swale soils) form primarily in aeolian dust, have larger amounts of total carbonate and soluble salts, and display greater variability than soils developing over original topographic highs (ridge soils). It is thus concluded that ridge soils, which have minimal variability, should be employed for a soil chronosequence study of basalt surfaces in the Potrillo volcanic field. The spatial variability of swale soils results in part from the significant hydrologic variability of low-lying landscape positions. Depth profiles of chloride concentrations suggest that hydrologic variability systematically correlates with the size and shape of depressions in which soils are forming. The infilling of depressions with aeolian material results in increasingly arid hydrologic conditions both by increasing the volume of aeolian material that is being drained and by reducing the catchment area for runoff into the depression. Depressions fill at different rates, however, depending on their size, shape and catchment area. Small, narrow depressions fill quickly, and their associated soils form under more arid conditions and have stronger development than soils in large depressions. Therefore, a number of geomorphic surfaces of varying age may develop on a single isochronous basalt flow. Each of these surfaces will have unique hydrologic characteristics and consequently different degrees of soil development. The pre-burial high water flux evident in depressions suggests that basalt flows may play an important role in aquifer recharge in this area of New Mexico. Copyright © 1999 John Wiley & Sons, Ltd.     

Boundary element simulation of scattering of elastic waves by 3-D cracks

Numerical modelling techniques are now becoming common for understanding the complicated nature of seismic wave propagation in fractured rock. Here the Indirect Boundary Element Method (IBEM) is applied to study scattering of elastic waves by cracks. The problem addressed in this paper is the diffraction of P and S waves by open 3-D cracks of arbitrary shape embedded in a homogeneous isotropic medium. The IBEM yields the value of the jump of displacements between opposite surfaces of the crack, often called Crack Opening Displacement (COD). This is used to evaluate the solution away from the crack. We use a multi-regional approach which consists of splitting a surface S into two identical surfaces S⁺ and S⁻ chosen such that the crack lies at the interface. The resulting integral equations are not hyper-singular and wave propagation within media that contain open cracks can be rigorously solved. In order to validate the method, we compare results of displacements of a penny-shaped crack for a vertical incident P-wave with the classic results by Mal (1970) obtaining excellent agreement. This comparison gives us confidence to study cases where no analytic solutions exist. Some examples of incidence of P or S waves upon cracks with various shapes are depicted and the salient aspects of the method are also discussed. Both frequency and time-domain results are included.     

Characteristics of turbulence over the semi-fixed desert area north of Xinjiang,China

As the largest fixed and semi-fixed desert in China, the Gurbantünggüt Desert undergoes a long period of snow cover in the winter and the rapid growth of ephemeral plants in the spring, presenting obvious seasonal changes in the underlying desert surface type, which can lead to variation in the turbulence of the near-surface boundary layer turbulence over the desert. In this study, gradient tower data and eddy covariance data from 2017 were analysed to investigate the turbulence characteristics of the different surface boundary layers in the hinterland of the Gurbantünggüt Desert. The results indicate that stable atmospheric conditions in the desert occur exclusively during the early morning and at night in the desert, and the onset and duration of this stable state varies seasonally. Two regimes of intermittent turbulence occur during the night, a weak turbulent regime that occurs when the wind speed is less than the threshold and a strong turbulent regime when the wind speed exceeds the threshold, and different wind speed thresholds were observed at each level. These parameters follow a seasonal pattern of summer (July) > spring (April) > autumn (October) > winter (January) in terms of magnitude. The mean turbulence intensities of the along-wind, cross-wind and vertical wind are 0.5, 0.47 and 0.14, respectively, with I_u > I_v > I_w. The normalized standard deviation of the wind velocity components (σ_u, σ_v and σ_w) generally satisfies a 1/3 power-law relation. Our results show that the night-time turbulence regime classification for the Gurbantünggüt Desert strongly depends on meteorological and orographic features, and the intermittent turbulent events have the non-stationarity of the flow in common. The results can contribute to the study of land surface processes, climate change and desertification in inland arid desert areas.