Evaluation of longitudinal and transverse dispersivities/distance ratios for tracer test in a radially convergent flow field with scale-dependent dispersion

This study presents a novel mathematical model for analysis of non-axisymmetrical solute transport in a radially convergent flow field with scale-dependent dispersion. A two-dimensional, scale-dependent advection–dispersion equation in cylindrical coordinates is derived based on assuming that the longitudinal and transverse dispersivities increase linearly with the distance of the solute transported from its injected source. The Laplace transform finite difference technique is applied to solve the two-dimensional, scale-dependent advection–dispersion equation with variable-dependent coefficients. Concentration contours for different times, breakthrough curves of average concentration over concentric circles with a fixed radial distance, and breakthrough curves of concentration at a fixed observation point obtained using the scale-dependent dispersivity model are compared with those from the constant dispersivity model. The salient features of scale-dependent dispersion are illustrated during the non-axisymmetrical transport from the injection well into extraction well in a convergent flow field. Numerical tests show that the scale-dependent dispersivity model predicts smaller spreading than the constant-dispersivity model near the source. The results also show that the constant dispersivity model can produce breakthrough curves of averaged concentration over concentric circles with the same shape as those from the proposed scale-dependent dispersivity model at observation point near the extraction well. Far from the extracting well, the two models predict concentration contours with significantly different shapes. The breakthrough curves at observation point near the injection well from constant dispersivity model always produce lesser overall transverse dispersion than those from scale-dependent dispersivity model. Erroneous dimensionless transverse/longitudinal dispersivity ratio may result from parametric techniques which assume a constant dispersivity if the dispersion process is characterized by a distance-dependent dispersivity relationship. A curve-fitting method with an example is proposed to evaluate longitudinal and transverse scale-proportional factors of a field with scale-dependent dispersion.     

Simulation of flow and solute coupled 2‐D problems with velocity‐dependent dispersion coefficient based on the network method

The governing equations of coupled density‐driven flow and solute transport problems in porous media, with velocity‐dependent dispersion coefficient, are strongly nonlinear and must be solved numerically. This contribution presents a network model, based on the network simulation method, capable of simulating the transient solution to this kind of problem efficiently and with a relatively low computational time. The mathematical model is formulated using the stream function and concentration variables. Simulation of the network model is carried out in the standard electric circuit simulation code, Pspice. The present model is first applied to simulate the original benchmark Henry problem, and the solution is compared with those obtained by other authors. A study of the grid size is also carried out. In addition, the modified version of Simpson and Clement of the Henry problem, as well as the groundwater flow in the closed desert basin of Pilot Valley, is studied using the proposed model. Copyright © 2011 John Wiley & Sons, Ltd.     

Analysis of solute transport in a divergent flow tracer test with scale‐dependent dispersion

It has been known for many years that dispersivities increase with solute displacement distance in a subsurface. The increase of dispersivities with solute travel distance results from significant variation in hydraulic properties of porous media and was identified in the literature as scale‐dependent dispersion. In this study, Laplace‐transformed analytical solutions to advection‐dispersion equations in cylindrical coordinates are derived for interpreting a divergent flow tracer test with a constant dispersivity and with a linear scale‐dependent dispersivity. Breakthrough curves obtained using the scale‐dependent dispersivity model are compared to breakthrough curves obtained from the constant dispersivity model to illustrate the salient features of scale‐dependent dispersion in a divergent flow tracer test. The analytical results reveal that the breakthrough curves at the specific location for the constant dispersivity model can produce the same shape as those from the scale‐dependent dispersivity model. This correspondence in curve shape between these two models occurs when the local dispersivity at an observation well in the scale‐dependent dispersivity model is 1·3 times greater than the constant dispersivity in the constant dispersivity model. To confirm this finding, a set of previously reported data is interpreted using both the scale‐dependent dispersivity model and the constant dispersivity model to distinguish the differences in scale dependence of estimated dispersivity from these two models. The analytical result reveals that previously reported dispersivity/distance ratios from the constant dispersivity model should be revised by multiplying these values by a factor of 1·3 for the scale‐dependent dispersion model if the dispersion process is more accurately characterized by scale‐dependent dispersion. Copyright © 2007 John Wiley & Sons, Ltd.     

Experimental study of the effect of roughness and Reynolds number on fluid flow in rough‐walled single fractures: a check of local cubic law

Local cubic law (LCL) is one of the most commonly applied physical laws for flow in single fractures (SF) and fractured media. The foundation of LCL is Darcian flow. This experimental study examines if LCL is valid for flow in a single rough fracture and how the fracture roughness and Reynolds number (Re) affect flow. Similar to the Moody diagram for flow in pipes, a diagram for flow in a single rough fracture has been generated to relate the friction coefficient with Re and the roughness. Under the experimental condition of this study, flow appears to be substantially different from Darcian flow. The flow law of q ∝ eⁿJ^m appears to be valid for describing the flow scheme where q, e, and J are the unit width flux, the average aperture, and the hydraulic gradient. The value of the power index m is found to be around 0·83 ～ 0·98, less than what has been used in Darcian flow (m = 1). The power index n is around 11·2 and 13·0, much greater than the n value used in the LCL (n = 3), and it increases with the average velocity. The Moody type of diagram shows that the friction factor for flow in SFs is influenced by Re and the roughness. It decreases with Re when Re is small, and becomes less sensitive to Re when Re is large enough. It also increases with the roughness. Copyright © 2010 John Wiley & Sons, Ltd.     

Determination of the conduit structure in a karst aquifer based on tracer data—Lurbach system,Austria

A structure model was used to analyse solute‐transport parameter estimates based on tracer breakthrough curves. In the model system, groundwater flow is envisioned to be organised in a complex conduit network providing a variety of short circuits with relative small carrying capacities along different erosion levels. The discharge through the fully filled conduits is limited owing to void geometries and turbulent flow; thus, a hierarchic overflow system evolves where conduits are (re‐)activated or dried up depending on the flow condition. Exemplified on the Lurbach–Tanneben karst aquifer, the applicability of the model approach was tested. Information derived from multi‐tracer experiments performed at different volumetric flow rates enabled to develop a structural model of the karst network, under constraint of the geomorphological and hydrological evolution of the site. Depending on the flow rate, groundwater is divided into up to eight flow paths. The spatial hierarchy of flow paths controls the sequence of flow path activation. Conduits of the topmost level are strongly influenced by reversible alteration processes. Sedimentation or blocking causes an overflow of water to the next higher conduit. Flow path specific dissolutional denudation rates were estimated using the temporal development of the partial discharge. Copyright © 2012 John Wiley & Sons, Ltd.     

The role of channel morphology on the mobility and dispersion of bed sediment in a small gravel‐bed stream

This study uses a unique 10‐year tracer dataset from a small gravel‐bed stream to examine bed mobility and sediment dispersion over long timescales and at a range of spatial scales. Seasonal tracer data that captured multiple mobilizing events was examined, while the effects of morphology on bed mobility and sediment dispersion were captured at three spatial scales: within morphological units (unit scale), between morphological units (reach scale) and between reaches with different channel morphologies (channel scale). This was achieved by analyzing both reach‐average mobility and travel distance data, as well as the development of ‘mobility maps’ that capture the spatial variability in tracer mobility within the channel. The tracer data suggest that sediment transport in East Creek remains near critical the majority of the time, with only rare large events resulting in high mobility rates and grain travel distances large enough to move sediment past dominant bedforms. While a variable capturing both the magnitude and frequency of flow events within a season yielded a better predictor to sediment mobility and dispersion than peak discharge alone, the distribution of events of different magnitude within the season played a large role in determining tracer mobility rates and travel distances. The effects of morphology differed depending on the analysis scale, demonstrating the importance of scale, and therefore study design, when examining the effect of morphology on sediment transport. Copyright © 2016 John Wiley & Sons, Ltd.     

Simulation of field‐scale water flow and bromide transport in a cracked clay soil

Single domain models may seriously underestimate leaching of nutrients and pesticides to groundwater in clay soils with shrinkage cracks. Various two‐domain models have been developed, either empirical or physically based, which take into account the effects of cracks on water flow and solute transport. This paper presents a model concept that uses the clay shrinkage characteristics to derive crack volume and crack depth under transient field conditions. The concept has been developed to simulate field average behaviour of a field with cracks, rather than flow and transport at a small plot. Water flow and solute transport are described with basic physics, which allow process and scenario analysis. The model concept is part of the more general agrohydrological model SWAP, and is applied to a field experiment on a cracked clay soil, at which water flow and bromide transport were measured during 572 days. A single domain model was not able to mimic the field‐average water flow and solute transport. Incorporation of the crack concept considerably improved the simulation of water content and bromide leaching to the groundwater. Still deviations existed between the measured and simulated bromide concentration profiles. The model did not reproduce the observed bromide retardation in the top layer and the high bromide dispersion resulting from water infiltration at various soil depths. A sensitivity analysis showed that the amounts of bromide leached were especially sensitive to the saturated hydraulic conductivity of the top layer, the solute transfer from the soil matrix to crack water flow and the mean residence time of rapid drainage. The shrinkage characteristic and the soil hydraulic properties of the clay matrix showed a low sensitivity. Copyright © 2000 John Wiley & Sons, Ltd.     

A simulation of large‐scale groundwater flow and travel time in a fractured rock environment for waste disposal purposes

Two‐dimensional regional groundwater flow was simulated based on a conceptual model of low‐permeability crystalline rocks of the Whiteshell Research Area (WRA) in south‐eastern Manitoba. The conceptual model consists of fracture zones that strike in different directions and dip at various angles in the background rock mass. The thickness and hydraulic properties of the fracture zones in the conceptual model were varied as were the fluid properties and the boundary conditions of the groundwater flow system. The effects of these variations on the groundwater flow pattern and on the convective travel time along pathways from a hypothetical disposal vault at 500 m depth to discharge locations at the ground surface were evaluated. The vault was located in the regional discharge area of the groundwater system. A homogeneous conceptual model of the WRA, having only freshwater flow, formed a groundwater flow pattern with a regional flow system. Local flow systems developed increasingly with the introduction of fracture zones 20 m and 3 m thick, and depth‐dependent fluid density. This indicates a reduction in groundwater residence time by fracture zones and fluid density. Flow pathways were analysed using both a stream‐function and a particle‐tracking technique. The pathways and their lengths from the location of the vault to the surface varied spatially according to the flow patterns. The minimum travel time along these pathways was less than 150 000 and greater than 4 000 000 years in models with and without fracture zones, respectively, indicating that the presence of fracture zones was the major controlling factor. A precise knowledge and refinement of conceptual model parameters is necessary during site selection for waste disposal purposes. Copyright © 2004 John Wiley & Sons, Ltd.     

Scaling properties and spatial organization of snow depth fields in sub‐alpine forest and alpine tundra

This paper reports on a study analysing the spatial distribution functions, the correlation structures, and the power spectral densities of high‐resolution LIDAR snow depths (～1 m) in two adjacent 500 m × 500 m areas in the Colorado Rocky Mountains, one a sub‐alpine forest the other an alpine tundra. It is shown how and why differences in the controlling physical processes induced by variations in vegetation cover and wind patterns lead to the observed differences in spatial organization between the snow depth fields of these environments. In the sub‐alpine forest area, the mean of snow depth increases with elevation, while its standard deviation remains uniform. In the tundra subarea, the mean of snow depth decreases with elevation, while its standard deviation varies over a wide range. The two‐dimensional correlations of snow depth in the forested area indicate little spatial memory and isotropic conditions, while in the tundra they indicate a marked directional bias that is consistent with the predominant wind directions and the location of topographic ridges and depressions. The power spectral densities exhibit a power law behaviour in two frequency intervals separated by a break located at a scale of around 12 m in the forested subarea, and 65 m in the tundra subarea. The spectral exponents obtained indicate that the snow depth fields are highly variable over scales larger than the scale break, while highly correlated below. Based on the observations and on synthetic snow depth fields generated with one‐ and two‐dimensional spectral techniques, it is shown that the scale at which the break occurs increases with the separation distance between snow depth maxima. In addition, the breaks in the forested area coincide with those of the corresponding vegetation height field, while in the tundra subarea they are displaced towards larger scales than those observed in the corresponding vegetation height field. Copyright © 2009 John Wiley & Sons, Ltd.     

Dynamics of suspended sediment transport at field‐scale drain channels of irrigation‐dominated watersheds in the Sonoran Desert,southeastern California

Suspended sediment is a major source of pollution in irrigation‐dominated watersheds. However, little is known about the process and mechanisms of suspended sediment transport in drain channels directly connected to agricultural fields. This paper explains sediment dynamics using averaged 5 min flow discharge Q (m³ s^?1) and suspended sediment concentration C (mg l^?1) collected during one crop season in a small catchment containing a first‐order drain channel and its connected six agricultural fields within the Salton Sea watershed. The statistical properties and average trends of Q and C were investigated for both early (i.e. November) and late (i.e. January) stages of a crop season. Further in‐depth analysis on sediment dynamics was performed by selecting two typical single‐field irrigation events and two multiple‐field irrigation events. For each set of irrigation events, the process of suspended sediment transport was revealed by examining hydrograph and sediment graph responses. The mechanisms underlying suspended sediment transport were investigated by analysing the types of corresponding hysteresis loop. Finally, sediment rating curves for both hourly and daily data at early and late stages and for the entire crop season were established to seek possible sediment‐transport predictive model(s). The study suggests that the complicated processes of suspended sediment transport in irrigation‐dominated watersheds require stochastic rather than deterministic forecasting. Copyright © 2007 John Wiley & Sons, Ltd.     

Local‐ and regional‐scale biomorphodynamics due to tree uprooting in semi‐natural and managed montane forests of the Sudetes Mountains,Central Europe

In this work, direct and indirect geomorphic consequences of wind‐related tree uprooting are examined, using an extensive dataset from the mountain range of the Sudetes, Poland. The role of local conditions in influencing the geomorphic efficacy of tree uprooting is examined, as well as issues of upscaling individual observations from experimental sites. This problem is approached at a range of spatial and observational scales, from monitoring of root plate degradation over time through to examination of wind effects at a slope scale and region‐wide analysis. In our study area the mean root plate volume is between 0.4 and 4.2 m³ for spruce and 2.4 m³ for beech, and their degradation may last tens of years. The density of relict pit‐and‐mound microtopography varies from 2.7 up to 40 pairs per hectare and the maximum coverage of terrain is 4.7%. The volume of treethrow mounds varies from 0.5 to 3.1 m³ and mounds seem to outlive the pits formed in the same episode of disturbance. However, in specific lithological and topographic conditions, pit‐and‐mound topography does not form. The maximum biogenic transport attributable to a single windstorm event is c. 80 m³ ha^?1, while soil turnover times are calculated in the order of 1000–10 000 years. Rock fragment ‘mining’ is an important biogeomorphic process, both in terms of impact on hillslope surfaces and on soil properties. Gravel armours and small‐scale stepped topography may form instead of typical pit–mound associations in specific circumstances. Managed forests appear more prone to wind damage and associated geomorphic consequences. In the Sudetes Mountains, the variable role of tree uprooting in local and regional hillslope denudation is governed by forest stand structure, topography and regolith properties, with the former significantly influenced by human activity. Copyright © 2016 John Wiley & Sons, Ltd.     

Catchment‐scale conceptual modelling of water and solute transport in the dual flow system of the karst critical zone

Hydrological and hydrochemical processes in the critical zone of karst environments are controlled by the fracture‐conduit network. Modelling hydrological and hydrochemical dynamics in such heterogeneous hydrogeological settings remains a research challenge. In this study, water and solute transport in the dual flow system of the karst critical zone were investigated in a 73.5‐km² catchment in southwest China. We developed a dual reservoir conceptual run‐off model combined with an autoregressive and moving average model with algorithms to assess dissolution rates in the “fast flow” and “slow flow” systems. This model was applied to 3 catchments with typical karst critical zone architectures, to show how flow exchange between fracture and conduit networks changes in relation to catchment storage dynamics. The flux of bidirectional water and solute exchange between the fissure and conduit system increases from the headwaters to the outfall due to the large area of the developed conduits and low hydraulic gradient in the lower catchment. Rainfall amounts have a significant influence on partitioning the relative proportions of flow and solutes derived from different sources reaching the underground outlet. The effect of rainfall on catchment function is modulated by the structure of the karst critical zone (e.g., epikarst and sinkholes). Thin epikarst and well‐developed sinkholes in the headwaters divert more surface water (younger water) into the underground channel network, leading to a higher fraction of rainfall recharge into the fast flow system and total outflow. Also, the contribution of carbonate weathering to mass export is also higher in the headwaters due to the infiltration of younger water with low solute concentrations through sinkholes.     

Investigation of the relationship between runoff and atmospheric oscillations,sea surface temperature,and local‐scale climate variables in the Yellow River headwaters region

The Yellow River headwaters region (YRHR) contributes nearly 40% of total flow in the Yellow River basin, which is suffering from a serious water shortage problem. Investigation of the relationship between runoff and climate variables is important for understanding the variation trend of runoff in the YRHR under global climate change. Global and local climate variables, including the West Pacific subtropical high; northern hemisphere polar vortex (NH); Tibetan Plateau Index B (TPI‐B); southern oscillation index; sea surface temperature; and precipitation, evaporation, and temperature, were fully considered to explore the relationship with runoff at Jimai, Maqu, and Tangnaihai stations from 1956 to 2014. The results reveal that runoff had a decreasing trend, which will likely be maintained in the future, and there was a significant change in runoff around 1995 at all stations. Correlation analysis indicated that runoff was dominated by precipitation, NH, temperature, and TPI‐B, and a substantial correlation was observed with sea surface temperature and evaporation, but there was little correlation with West Pacific subtropical high and southern oscillation index. Furthermore, impacts of climate change on runoff variations were distinctly different at different temporal scales. Three dominant runoff periodicities were identified by a singular spectrum analysis‐multitaper method and continuous wavelet transform, that is, 1.0‐, 6.9‐, and 24.8‐year runoff periodicities. In addition, runoff was positively correlated with temperature at a 1‐year periodicity, negatively correlated with TPI‐B at a 6.9‐year periodicity, and positively correlated with NH at a 24.8‐year periodicity, that is, temperature, TPI‐B, and NH‐controlled runoff at annual, interannual, and interdecadal scales. Further, all analyses of the stations in the YRHR showed excellent consistency. The results will provide valuable information for water resource management in the YRHR.     

Implications of climate change in the twenty‐first century for simulated magnitude and frequency of bed‐material transport in tributaries of the Saint‐Lawrence River

More frequent extreme flood events are likely to occur in many areas in the twenty‐first century due to climate change. The impacts of these changes on sediment transport are examined at the event scale using a 1D morphodynamic model (SEDROUT4‐M) for three tributaries of the Saint‐Lawrence River (Québec, Canada) using daily discharge series generated with a hydrological model (HSAMI) from three global climate models (GCMs). For all tributaries, larger flood events occur in all future scenarios, leading to increases in bed‐material transport rates, number of transport events and number of days in the year where sediment transport occurs. The effective and half‐load discharges increase under all GCM simulations. Differences in flood timing within the tributaries, with a shift of peak annual discharge from the spring towards the winter, compared to the hydrograph of the Saint‐Lawrence River, generate higher sediment transport rates because of increased water surface slope and stream power. Previous research had shown that channel erosion is expected under all GCMs' discharge scenarios. This study shows that, despite lower bed elevations, flood risk is likely to increase as a result of higher flood magnitude, even with falling base level in the Saint‐Lawrence River. Copyright © 2010 John Wiley & Sons, Ltd.     

Discriminating between the roles of late Pleistocene palaeodischarge and geological‐topographic inheritance in fluvial longitudinal profile and channel development

This paper investigates the influences of palaeohydrology and geological‐topographic inheritance in shaping the channel of the lower River Suir, southeast Ireland. Results of acoustic surveys of the lower River Suir and Waterford Harbour reveal two scales of pseudo‐cyclic river bedforms. Longitudinal elevation profiles of the geological topography (undulating bedrock and till‐mantled bedrock) bounding the present floodplain swath reveal pseudo‐cyclicity in that terrain too. Spectral and statistical analyses are used to quantify the cyclicity of the long profile and geological‐topographic series. These methods show that the dominant cyclicity of the long profile reflects autocorrelation more than inheritance of cyclicity from the bounding geological topography. The cyclicity of the long profile mainly reflects a hydraulic control on pool‐spacing, although some cyclicity probably has been inherited from the geological‐topography. Channel‐forming palaeodischarge is estimated based on the dominant pool‐spacing revealed by spectral analysis, validated using relationships between meander wavelength, channel cross‐sectional geometry and hydraulically‐informed discharge reconstruction. The palaeodischarge estimates are in close agreement and are two orders of magnitude greater than present flood maxima. Significantly, these palaeodischarge estimates also agree closely with palaeodischarge calculated for the submerged Pleistocene palaeochannel that extends across the near‐shore continental shelf from Waterford Harbour. The pool‐sequence of the lower Suir and the submerged palaeochannel represent a former land‐system that was active during a period of low relative sea level during the last glacial. More broadly, the paper offers insights into the landscape evolution of formerly glaciated regions that experienced very wide discharge variability during and after the transition from glacial to interglacial regimes, in a context of complex relative sea level change. Copyright © 2017 John Wiley & Sons, Ltd.     

Spatial and temporal analysis of hillslope–channel coupling and implications for the longitudinal profile in a dryland basin

The long‐term evolution of channel longitudinal profiles within drainage basins is partly determined by the relative balance of hillslope sediment supply to channels and the evacuation of channel sediment. However, the lack of theoretical understanding of the physical processes of hillslope–channel coupling makes it challenging to determine whether hillslope sediment supply or channel sediment evacuation dominates over different timescales and how this balance affects bed elevation locally along the longitudinal profile. In this paper, we develop a framework for inferring the relative dominance of hillslope sediment supply to the channel versus channel sediment evacuation, over a range of temporal and spatial scales. The framework combines distinct local flow distributions on hillslopes and in the channel with surface grain‐size distributions. We use these to compute local hydraulic stresses at various hillslope‐channel coupling locations within the Walnut Gulch Experimental Watershed (WGEW) in southeast Arizona, USA. These stresses are then assessed as a local net balance of geomorphic work between hillslopes and channel for a range of flow conditions generalizing decadal historical records. Our analysis reveals that, although the magnitude of hydraulic stress in the channel is consistently higher than that on hillslopes, the product of stress magnitude and frequency results in a close balance between hillslope supply and channel evacuation for the most frequent flows. Only at less frequent, high‐magnitude flows do channel hydraulic stresses exceed those on hillslopes, and channel evacuation dominates the net balance. This result suggests that WGEW exists mostly (~50% of the time) in an equilibrium condition of sediment balance between hillslopes and channels, which helps to explain the observed straight longitudinal profile. We illustrate how this balance can be upset by climate changes that differentially affect relative flow regimes on slopes and in channels. Such changes can push the long profile into a convex or concave condition. © 2018 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Large‐scale sediment bodies and superimposed bedforms on the continental shelf close to the Strait of Gibraltar: interplay of complex oceanographic conditions and physiographic constraints

A previously unknown field of large‐scale sedimentary bodies has been mapped and studied on the continental shelf off the Cape Trafalgar near the Strait of Gibraltar with particular emphasis on the relationship between large‐scale sediment bodies and the superimposed bedforms. This study is based on a grid of 975 km of high‐resolution seismic profiles collected at water depths ranging between 15 and 60 m. High variability of large‐scale sedimentary bodies is attributed to the complex interaction of hydrodynamic agents. The most prominent sedimentary features are sand banks and ridges that indicate long‐term southwest and southward‐directed sediment transport patterns, possibly due to the interplay of two dominant current systems flowing southward and westward. These sediment bodies evolve laterally to distinct external geometries, such as sand shoals in shallow water and sand sheets in the vicinity of larger sand banks that indicate moderate current velocities. In addition, pre‐existing physiography is considered to play a role in the generation of certain sediment bodies, developed over inclined surfaces or confined laterally by elevations. Relationships between superimposed bedforms (mostly very large dunes) and underlying sediment bodies vary across the study area. Most superimposed bedforms occur over the complex mosaic of sediment banks and sheets, suggesting the interaction of several high‐energy currents with different directions, such as tidal and/or wind‐driven currents. Copyright © 2010 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.     

Self‐adjustment in rivers: Evidence for least action as the primary control of alluvial‐channel form and process

Least action principle (LAP) in rivers is demonstrated by maximum flow efficiency (MFE) and is the foundation of variational mechanics based on energy and work rather than Newtonian force and momentum. Empirical evidence shows it to be the primary control for the adjustment of alluvial channels. Because most rivers flow with imposed water and sediment loads down valley gradients they have largely inherited, they self‐regulate energy expenditure to match the work they are required to do to remain stable. Overpowered systems develop a variety of channel patterns to expend excess energy and remain stable. Australia offers an opportunity to study low‐energy rivers closely adjusted to very low continental gradients. The anabranching Marshall and single‐thread Plenty Rivers flow down nearly straight channels with average H numbers [ratio between excess bed shear and width/depth (W/D) ratio] close to the optimum of 0.3 for stationary equilibrium. Ridge‐form divisions of the original channel width create anabranches that radically alter W/D ratios relative to bed shear, the same being true for short‐wide islands on the large low‐gradient Yangtze River in China. In contrast, Mount Chambers Creek in Australia's tectonically more active Flinders Ranges is accreting an alluvial fan with unstable distributary channels exhibiting H numbers well below the optimum. LAP also explains profound biases in Earth's stratigraphic record. Because meandering is an energy‐shedding mechanism, sinuous rivers sequester relatively little sediment resulting in all sequences being just a few tens of metres thick. In contrast, low‐energy braided disequilibrium systems can sequester sediment piles over a kilometre in thickness and tens of kilometres wide. LAP provides a new paradigm for river research by identifying the attractor state controlling river channel evolution. It links advances in theoretical physics to fluvial geomorphology, stratigraphy and hydraulic engineering and opens opportunities for diverse investigations in Earth system science. Copyright © 2016 John Wiley & Sons, Ltd.