An experimental study on the ripple–dune transition

Flume experiments were conducted on different bed stages across the ripple–dune transition. As flow velocity increases, an initially flat bed surface (made of fairly uniform sandy material) is gradually transformed into a two‐dimensional rippled bed. With further increase in velocity, two‐dimensional ripples are replaced by irregular, linguoid ripples. As the average velocity necessary for the ripple–dune transition to occur is imposed on the bed surface, these non‐equilibrium linguoid ripples are further transformed into larger, two‐dimensional dunes. For each of these stages across the transition, a concrete mould of the bed was created and the flow structure above each fixed bed surface investigated. An acoustic Doppler velocimeter was used to study the flow characteristics above each bed surface. Detailed profiles were used along a transect located in the middle of the channel. Results are presented in the form of spatially averaged profiles of various flow characteristics and of contour maps of flow fields (section view). They clearly illustrate some important distinctions in the flow structure above the different bedform types associated with different stages during the transition. Turbulence intensity and Reynolds stresses gradually increase throughout the transition. Two‐dimensional ripples present a fairly uniform spatial distribution of turbulent flow characteristics above the bed. Linguoid ripples induce three‐dimensional turbulence structure at greater heights above the bed surface and turbulence intensity tends to increase steadily with height above bed surface in the wake region. A very significant increase in turbulence intensity and momentum exchange occurs during the transition from linguoid ripples to dunes. The turbulent flow field properties above dunes are highly dependent on the position along and above the bed surface and these fields present a very high degree of spatial variability (when compared with the rippled beds). Further investigations under natural conditions emphasizing sediment transport mechanisms and rates during the transition should represent the next step of analysis, together with an emphasis on quadrant analysis. Copyright © 2001 John Wiley & Sons, Ltd.     

Assessment of 2DH and pseudo‐3D modelling platforms in a large saline aquatic system: Lake Urmia,Iran

The main objective of this paper is to provide comparative quantitative examinations on the capabilities of two‐dimensional horizontal and pseudo‐three‐dimensional (3D) modelling approaches for simulating spatial and temporal variability of the flow and salinity in Lake Urmia, Iran. The water quality in the lake has been an environmentally important subject partly because this shallow hypersaline aquatic ecosystem is considered to be one of the largest natural habitats of a unique multicellular organism, Artemia urmiana. This brine shrimp is the major food source for many of the protected and rare shorebirds that visit the lake. A. urmiana can grow and survive in certain ranges of salinity, and their disappearance could lead to an alteration of existing equilibria. The lake has also experienced considerable man‐made changes during the past three decades. A newly built crossing embankment almost divided the lake into two northern and southern halves. A relatively small opening of 1.25 km in the new embankment provides water connections between the two halves. As a result, the flow and salinity regimes have been significantly changed. This might have had adverse serious impacts on the lake ecosystem. In the current study, the two‐dimensional horizontal hydrodynamic model has been found to provide reasonable predictions for the flow regime in the lake, whereas its salinity predictions have not been consistent with the field observations. The pseudo‐3D model has produced results fairly close to the salinity measurements and its temporal and spatial variations. The pseudo‐3D model has been used for evaluating the embankment effects on the lake hydrodynamics and on the salinity conditions. The effectiveness of introducing a different number or length of openings in the embankment for restoring the pre‐embankment conditions has also been examined. These remedy options have been found not to offer substantial improvements to the lake's existing ecosystem. Copyright © 2013 John Wiley & Sons, Ltd.     

Variations in thickness of fault‐controlled dikes and its implications: A case study of Gosung area,South‐East Korea

Dikes are natural records that can be used to understand the way magma flows in the crust. Coastal platform outcrops in Gosung, South Korea, show clear evidences that their intrusion took place along pre‐existing fractures. We analyzed outcropping dikes, measuring variations in dike thickness as well as fracture density (cumulative number of fractures along strike) and geometry around the dikes. The geometry and thickness variations of dikes intruded along pre‐existing fractures can be interpreted to understand the effect of pre‐existing fractures to evolution on magma flow, especially related with fault damage zones. This helps us to gain a better understanding of magma and fluid flow along pre‐existing fractures. Magma flow is greater along planes that strike perpendicular to the direction of least compressive horizontal stress, and along well connected fractures that show a high degree of connectivity. At the fault tip and linkage damage zone, there is a concentration of extensional fractures; in these areas injected dikelets can form. As faults become linked, the fracture density increases, until they become fully linked and act as one through‐going fault plane. As faults evolve, the boundary conditions of the faults vary and this has an impact on dike characteristics. Fracture geometry around dikes that intruded pre‐existing faults can be used as a record of fault evolution and this can give insights into how the maturity of a fault system can affect to the related magma or fluid flow characteristics.     

A design‐variable‐based inelastic hysteretic model for beam–column connections

This paper presents a design‐variable‐based inelastic hysteretic model for beam–column connections. It has been well known that the load‐carrying capacity of connections heavily depends on the types and design variables even in the same connection type. Although many hysteretic connection models have been proposed, most of them are dependent on the specific connection type with presumed failure mechanisms. The proposed model can be responsive to variations both in design choices and in loading conditions. The proposed model consists of two modules: physical‐principle‐based module and neural network (NN)‐based module in which information flow from design space to response space is formulated in one complete model. Moreover, owing to robust learning capability of a new NN‐based module, the model can also learn complex dynamic evolutions in response space under earthquake loading conditions, such as yielding, post‐buckling and tearing, etc. Performance of the proposed model has been demonstrated with synthetic and experimental data of two connection types: extended‐end‐plate and top‐ and seat‐angle with double‐web‐angle connection. Furthermore, the design‐variable‐based model can be customized to any structural component beyond the application to beam–column connections. Copyright © 2007 John Wiley & Sons, Ltd.     

Out‐of‐plane effects in 2D borehole‐to‐surface resistivity tomography and applications in mineral exploration

In this paper, we discuss the effects of anomalous out‐of‐plane bodies in two‐dimensional (2D) borehole‐to‐surface electrical resistivity tomography with numerical resistivity modelling and synthetic inversion tests. The results of the two groups of synthetic resistivity model tests illustrate that anomalous bodies out of the plane of interest have an effect on two‐dimensional inversion and that the degree of influence of out‐of‐plane body on inverted images varies. The different influences are derived from two cases. One case is different resistivity models with the same electrode array, and the other case is the same resistivity model with different electrode arrays. Qualitative interpretation based on the inversion tests shows that we cannot find a reasonable electrode array to determine the best inverse solution and reveal the subsurface resistivity distribution for all types of geoelectrical models. Because of the three‐dimensional effect arising from neighbouring anomalous bodies, the qualitative interpretation of inverted images from the two‐dimensional inversion of electrical resistivity tomography data without prior information can be misleading. Two‐dimensional inversion with drilling data can decrease the three‐dimensional effect. We employed two‐ and three‐dimensional borehole‐to‐surface electrical resistivity tomography methods with a pole–pole array and a bipole–bipole array for mineral exploration at Abag Banner and Hexigten Banner in Inner Mongolia, China. Different inverse schemes were carried out for different cases. The subsurface resistivity distribution obtained from the two‐dimensional inversion of the field electrical resistivity tomography data with sufficient prior information, such as drilling data and other non‐electrical data, can better describe the actual geological situation. When there is not enough prior information to carry out constrained two‐dimensional inversion, the three‐dimensional electrical resistivity tomography survey is the better choice.     

Coexisting different types of zoned garnet in kyanite‐quartz eclogites from the Sanbagawa metamorphic belt: Evidence of deformation‐induced lithological mixing during prograde metamorphism

Garnet grains in Sanbagawa quartz eclogites from the Besshi region, central Shikoku commonly show a zoning pattern consisting of core and mantle/rim that formed during two prograde stages of eclogite and subsequent epidote–amphibolite facies metamorphism, respectively. Garnet grains in the quartz eclogites are grouped into four types (I, II, III, and IV) according to the compositional trends of their cores. Type I garnet is most common and sometimes coexists with other types of garnet in a thin section. Type I core formed with epidote and kyanite during the prograde eclogite facies stage. The inner cores of types II and III crystallized within different whole‐rock compositions of epidote‐free and kyanite‐bearing eclogite and epidote‐ and kyanite‐free eclogite at the earlier prograde stage, respectively. The inner core of type IV probably formed during the pre‐eclogite facies stage. The inner cores of types II, III, and IV, which formed under different P–T conditions of prograde metamorphism and/or whole‐rock compositions, were juxtaposed with the core of type I, probably due to tectonic mixing of rocks at various points during the prograde eclogite facies stage. After these processes, they have shared the following same growth history: (i) successive crystal growth during the later stage of prograde eclogite facies metamorphism that formed the margin of the type I core and the outer cores of types II, III, and IV; (ii) partial resorption of the core during exhumation and hydration stage; and (iii) subsequent formation of mantle zones during prograde metamorphism of the epidote–amphibolite facies. The prograde metamorphic reactions may not have progressed under an isochemical condition in some Sanbagawa metamorphic rocks, at least at the hand specimen scale. This interpretation suggests that, in some cases, material interaction promoted by mechanical mixing and fluid‐assisted diffusive mass transfer probably influences mineral reactions and paragenesis of high‐pressure metamorphic rocks.     

Macroturbulent coherent structures in an ice‐covered river flow using a pulse‐coherent acoustic Doppler profiler

The aim of this work is to compare macroturbulent coherent structures (MCS) geometry and organization between ice covered and open channel flow conditions. Velocity profiles were obtained using a Pulse‐Coherent Acoustic Doppler Profiler in both open channel and ice‐covered conditions. The friction imposed by the ice cover results in parabolic shaped velocity profiles. Reynolds stresses in the streamwise (u) and vertical (v) components of the flow show positive values near the channel bed and negative values near the ice cover, with two distinctive boundary layers with specific turbulent signatures. Vertically aligned stripes of coherent flow motions were revealed from statistics applied to space‐time matrices of flow velocities. In open channel conditions, the macroturbulent structures extended over the entire depth of the flow whereas they were discontinued and nested close to the boundary walls in ice‐covered conditions. The size of MCS is consequently reduced in scale under an ice cover. The average streamwise length scale is reduced from 2.5 to 0.4Y (u) and from 1.5 to 0.4Y (v) where Y is the flow depth. In open channel conditions, the vertical extent of MCS covers the entire flow depth, whereas the vertical extent was in the range 0.58Y–1Y (u) and 0.81Y–1Y (v) in ice‐covered conditions. Under an ice cover, each boundary wall generates its own set of MCS that compete with each other in the outer region of the flow, enhancing mixing and promoting the dissipation of coherent structures. Copyright © 2012 John Wiley & Sons, Ltd.     

Two‐dimensional physically based finite element runoff model for small agricultural watersheds: I. Model development

Soil erosion by water is the root cause of ecological degradation in the Shiwalik foothills of Northern India. Simulation of runoff and its component processes is a pre‐requisite to develop the management strategies to tackle the problem, successfully. A two‐dimensional physically based distributed numerical model, ROMO2D has been developed to simulate runoff from small agricultural watersheds on an event basis. The model employs the 2‐D Richards equation with sink term to simulate infiltration and soil moisture dynamics in the vadoze zone under variable rainfall conditions, and 2‐D Saint‐Venant equations under the kinematic wave approximation along with Manning's equation as the stage‐discharge equation for runoff routing. The various flow‐governing equations have been solved numerically by employing a Galerkin finite element method for spatial discretization using quadrilateral elements and finite difference techniques for temporal solutions. The ROMO2D computer program has been developed as a class‐based program, coded in C + + in such a way that with minor modifications, the model can be used to simulate runoff on a continuous basis. The model writes output for a runoff hydrograph of each storm. Model development is described in this paper and the results of model testing and field application are to be presented in a subsequent paper. Copyright © 2008 John Wiley & Sons, Ltd.     

Infiltration into frozen soil: From core‐scale dynamics to hillslope‐scale connectivity

Infiltration into frozen soil is a key hydrological process in cold regions. Although the mechanisms behind point‐scale infiltration into frozen soil are relatively well understood, questions remain about upscaling point‐scale results to estimate hillslope‐scale run‐off generation. Here, we tackle this question by combining laboratory, field, and modelling experiments. Six large (0.30‐m diameter by 0.35‐m deep) soil cores were extracted from an experimental hillslope on the Canadian Prairies. In the laboratory, we measured run‐off and infiltration rates of the cores for two antecedent moisture conditions under snowmelt rates and diurnal freeze–thaw conditions observed on the same hillslope. We combined the infiltration data with spatially variable data from the hillslope, to parameterise a surface run‐off redistribution model. We used the model to determine how spatial patterns of soil water content, snowpack water equivalent (SWE), and snowmelt rates affect the spatial variability of infiltration and hydrological connectivity over frozen soil. Our experiments showed that antecedent moisture conditions of the frozen soil affected infiltration rates by limiting the initial soil storage capacity and infiltration front penetration depth. However, shallow depths of infiltration and refreezing created saturated conditions at the surface for dry and wet antecedent conditions, resulting in similar final infiltration rates (0.3 mm hr^?1). On the hillslope‐scale, the spatial variability of snowmelt rates controlled the development of hydrological connectivity during the 2014 spring melt, whereas SWE and antecedent soil moisture were unimportant. Geostatistical analysis showed that this was because SWE variability and antecedent moisture variability occurred at distances shorter than that of topographic variability, whereas melt variability occurred at distances longer than that of topographic variability. The importance of spatial controls will shift for differing locations and winter conditions. Overall, our results suggest that run‐off connectivity is determined by (a) a pre‐fill phase, during which a thin surface soil layer wets up, refreezes, and saturates, before infiltration excess run‐off is generated and (b) a subsequent fill‐and‐spill phase on the surface that drives hillslope‐scale run‐off.     

A new passive rolling‐pendulum vibration absorber using a non‐axial‐symmetrical guide to achieve bidirectional tuning

A new, passive, vibroprotective device of the rolling‐pendulum tuned mass damper type is presented that, relying on a proper three‐dimensional guiding surface, can simultaneously control the response of the supporting structure in two mutually orthogonal horizontal directions. Unlike existing examples of ball vibration absorbers, mounted on spherical recesses and effective for axial‐symmetrical structures, the new device is bidirectionally tuneable, by virtue of the optimum shape of the rolling cavity, to both fundamental structural modes, even when the corresponding natural frequencies are different, in such a case recurring to an innovative non‐axial‐symmetrical rolling guide. Based on Appell's non‐holonomic mechanics, a non‐linear dynamic model is first derived for the bidirectional absorber mounted on a 1‐storey 3‐degrees‐of‐freedom linear structure translating under the effect of both imparted base motion and applied dynamic forces. A laboratory‐scaled prototype of the device is then tested to experimentally demonstrate the bidirectional tuning capability and to validate the mathematical model. The design procedure and the seismic performance of the absorber are finally exemplified through numerical simulation. Copyright © 2009 John Wiley & Sons, Ltd.     

Scaling and self‐similarity in one‐dimensional unsteady open channel flow

The conditions under which the Saint Venant equations system for unsteady open channel flow, as an initial–boundary value problem, becomes self‐similar are investigated by utilizing one‐parameter Lie group of point scaling transformations. One of the advantages of this methodology is that the self‐similarity conditions due to the initial and boundary conditions can also be investigated thoroughly in addition to the conditions due to the governing equation. The obtained self‐similarity conditions are compared with the scaling relationships that are derived through the Froude similitude. It is shown that the initial–boundary value problem of a one‐dimensional unsteady open channel flow process in a prototype domain can be self‐similar with that of several different scaled domains. However, the values of all the flow variables (at specified time and space) under different scaled domains can be upscaled to the same values in the prototype domain (at the corresponding time and space), as shown in this study. Distortion in scales of different space dimensions has been implemented extensively in physical hydraulic modelling, mainly because of cost, space and time limitations. Unlike the traditional approach, the distinction is made between the longitudinal–horizontal and transverse–horizontal length scales in this study. The scaled domain obtained by the proposed approach, when scaling ratios of channel width and water depth are equal, is particularly important for the similarity of flow characteristics in a cross‐section because the width‐to‐depth ratio and the inclination angles of the banks are conserved in a cross‐section. It is also shown that the scaling ratio of the roughness coefficient under distorted channel conditions depends on that of hydraulic radius and longitudinal length. The proposed scaling relations obtained by the Lie group scaling approach may provide additional spatial, temporal and economical flexibility in setting up physical hydraulic models. Copyright © 2013 John Wiley & Sons, Ltd.     

Closed‐Form Approximations for Two‐Dimensional Groundwater Age Patterns in a Fresh Water Lens

Simple closed‐form approximations are presented for calculating the steady‐state groundwater age distribution in two‐dimensional vertical cross sections of idealized fresh water lenses overlying salt water, for aquifers that are vertically semi‐infinite and of finite thickness. The approximations are developed on the basis of existing one‐dimensional analytical solutions for travel‐time calculation in fresh water lenses and approximate streamline formulations. The two‐dimensional age distributions based on the closed‐form solutions match convincingly with numerical simulations. As expected, notable deviations from the numerical solution are encountered at the groundwater flow divide and when submarine groundwater discharge occurs. Ratios of recharge over hydraulic conductivities are varied to explore how the magnitude of the deviations changes, and it is found that the approximate closed‐form solutions perform well over a range of conditions found in natural systems.     

An analytical approach to ascertain saturation‐excess versus infiltration‐excess overland flow in urban and reference landscapes

Uncontrolled overland flow drives flooding, erosion, and contaminant transport, with the severity of these outcomes often amplified in urban areas. In pervious media such as urban soils, overland flow is initiated via either infiltration‐excess (where precipitation rate exceeds infiltration capacity) or saturation‐excess (when precipitation volume exceeds soil profile storage) mechanisms. These processes call for different management strategies, making it important for municipalities to discern between them. In this study, we derived a generalized one‐dimensional model that distinguishes between infiltration‐excess overland flow (IEOF) and saturation‐excess overland flow (SEOF) using Green–Ampt infiltration concepts. Next, we applied this model to estimate overland flow generation from pervious areas in 11 U.S. cities. We used rainfall forcing that represented low‐ and high‐intensity events and compared responses among measured urban versus predevelopment reference soil hydraulic properties. The derivation showed that the propensity for IEOF versus SEOF is related to the equivalence between two nondimensional ratios: (a) precipitation rate to depth‐weighted hydraulic conductivity and (b) depth of soil profile restrictive layer to soil capillary potential. Across all cities, reference soil profiles were associated with greater IEOF for the high‐intensity set of storms, and urbanized soil profiles tended towards production of SEOF during the lower intensity set of storms. Urban soils produced more cumulative overland flow as a fraction of cumulative precipitation than did reference soils, particularly under conditions associated with SEOF. These results will assist cities in identifying the type and extent of interventions needed to manage storm water produced from pervious areas.     

Modeling of a large‐scale magneto‐rheological damper for seismic hazard mitigation. Part I: Passive mode

Magneto‐rheological (MR) dampers are a promising device for seismic hazard mitigation because their damping characteristics can be varied adaptively using an appropriate control law. During the last few decades researchers have investigated the behavior of MR dampers and semi‐active control laws associated with these types of dampers for earthquake hazard mitigation. A majority of this research has involved small‐scale MR dampers. To investigate the dynamic behavior of a large‐scale MR damper, characterization tests were conducted at the Lehigh Network for Earthquake Engineering Simulation equipment site on large‐scale MR dampers. A new MR damper model, called the Maxwell Nonlinear Slider (MNS) model, is developed based on the characterization tests and is reported in this paper. The MNS model can independently describe the pre‐yield and post‐yield behavior of an MR damper, which makes it easy to identify the model parameters. The MNS model utilizes Hershel–Bulkley visco‐plasticity to describe the post‐yield non‐Newtonian fluid behavior, that is, shear thinning and thickening behavior, of the MR fluid that occurs in the dampers. The predicted response of a large‐scale damper from the MNS model along with that from existing Bouc–Wen and hyperbolic tangent models, are compared with measured response from various experiments. The comparisons show that the MNS model achieves better accuracy than the existing models in predicting damper response under cyclic loading. Copyright © 2012 John Wiley & Sons, Ltd.     

A two‐dimensional analytical solution for groundwater flow in a leaky confined aquifer system near open tidal water

Groundwater in coastal areas is commonly disturbed by tidal fluctuations. A two‐dimensional analytical solution is derived to describe the groundwater fluctuation in a leaky confined aquifer system near open tidal water under the assumption that the groundwater head in the confined aquifer fluctuates in response to sea tide whereas that of the overlying unconfined aquifer remains constant. The analytical solution presented here is an extension of the solution by Sun for two‐dimensional groundwater flow in a confined aquifer and the solution by Jiao and Tang for one‐dimensional groundwater flow in a leaky confined aquifer. The analytical solution is compared with a two‐dimensional finite difference solution. On the basis of the analytical solution, the groundwater head distribution in a leaky confined aquifer in response to tidal boundaries is examined and the influence of leakage on groundwater fluctuation is discussed. Copyright © 2001 John Wiley & Sons, Ltd.     

Contaminant transport in one‐dimensional single fractured media: semi‐analytical solution for three‐member decay chain with pulse and Heaviside input sources

A semi‐analytical solution of the one‐dimensional transport for considering a three‐member decay chain in a single fracture with pulse and Heaviside input sources has been studied using the Laplace transform and its numerical inversion. The results reveal that breakthrough curves of dimensionless concentration for the decay chain of Np‐237, U‐233, and Th‐229 in the fracture can be well demonstrated in the temporal and spatial domains. The conditions with and without retardation effects are also compared. During the preliminary screening phase the solutions are suitable for performance assessment on radioactive waste disposal sites under a one‐dimensional single fracture condition. Copyright © 2007 John Wiley & Sons, Ltd.     

Analysis of benzene transport in a two‐dimensional aquifer model

In this study, the fate and transport of aqueous benzene was investigated in a laboratory‐scale homogeneous aquifer by conducting a two‐dimensional plume test. Benzene solution was introduced as a pulse type along the width of the aquifer model through a recharge zone situated at the upper‐left part of the model and followed by a steady state flow. Solution samples were collected at various locations on the front side of the model to capture two‐dimensional plumes at discrete time intervals. The benzene plumes showed a moderate retardation relative to chloride plumes observed from the previous study conducted for the same aquifer model. The retardation factor was obtained from the ratio of travel distances of benzene peaks to chloride peaks from the injection point, computed using a line integral method. Mass recovery of aqueous benzene revealed that there was a significant reduction of benzene mass, indicating the occurrence of volatilization and/or irreversible sorption during transport. Thus, retardation along with volatilization and/or irreversible sorption may be important processes affecting the fate and transport of aqueous benzene in the aquifer model. Copyright © 2005 John Wiley & Sons, Ltd.     

Estimation of Salt Water Upconing Using a Steady‐State Solution for Partial Completion of a Pumped Well

A new steady‐state analytical solution to the two‐dimensional radial‐flow equation was developed for drawdown (head) conditions in an aquifer with constant transmissivity, no‐flow conditions at the top and bottom, constant head conditions at a known radial distance, and a partially completed pumping well. The solution was evaluated for accuracy by comparison to numerical simulations using MODFLOW. The solution was then used to estimate the rise of the salt water‐fresh water interface (upconing) that occurs under a pumping well, and to calculate the critical pumping rate at which the interface becomes unstable, allowing salt water to enter the pumping well. The analysis of salt water‐fresh water interface rise assumed no significant effect on upconing by recharge; this assumption was tested and supported using results from a new steady‐state analytical solution developed for recharge under two‐dimensional radial‐flow conditions. The upconing analysis results were evaluated for accuracy by comparison to those from numerical simulations using SEAWAT for salt water‐fresh water interface positions under mild pumping conditions. The results from the equation were also compared with those of a published numerical sharp‐interface model applied to a case on Cape Cod, Massachusetts. This comparison indicates that estimating the interface rise and maximum allowable pumping rate using the analytical method will likely be less conservative than the maximum allowable pumping rate and maximum stable interface rise from a numerical sharp‐interface model.     

Modelling two‐dimensional steady‐state groundwater flow and flow sensitivity to boundary conditions in blanket peat complexes

This study used a two‐dimensional steady‐state finite‐element groundwater flow model to simulate groundwater flow in two Newfoundland blanket peat complexes and to examine flow system sensitivity to changes in water table recharge and aquifer properties. The modelling results were examined within the context of peat‐forming processes in the two complexes. Modelled flow compared favourably with observed flow. The sensitivity analyses suggested that more highly decomposed bog peat along bog margins probably has/had a positive impact on net peat accumulation within bog interiors. Peat with lower hydraulic conductivity along bog margins effectively impedes lateral drainage, localizes water table drawdown to extreme bog margins, and elevates water tables along bog interiors. Peat formation and elevated water tables in adjacent poor fens/laggs currently rely on placic and ortstein horizons impeding vertical drainage and water flow inputs from adjacent bogs. Modest reductions in atmospheric recharge were found to govern bog‐flow‐system geometries in a way that would adversely affect paludification processes in adjacent fens/laggs. Copyright © 2004 John Wiley & Sons, Ltd.