A stability criterion inherent in laws governing alluvial channel flow

The stability criterion of maximum flow efficiency (MFE) has previously been found inherent in typical alluvial channel flow relationships, and this study investigates the general nature of this criterion using a wider range of flow resistance and bedload transport formulae. For straight alluvial channels, in which the effect of sediment sorting is insignificant, our detailed mathematical analysis demonstrates that a flow efficiency factor ε occurs generally as the ratio of sediment (bedload) discharge Q_s to stream power Ω (γQS) in the form of . When ε is maximized (i.e. Q_s is maximized or Ω is minimized), maximally efficient straight channel geometries derived from most flow resistance and bedload transport formulae are found compatible with observed bankfull hydraulic geometry relations. This study provides support for the use of the criteria of MFE, maximum sediment transporting capacity and minimum stream power for understanding the operation of alluvial rivers, and also addresses limitations to the direct application of its findings. Copyright © 2002 John Wiley & Sons, Ltd.     

Field evidence for the influence of weathering on rock erodibility and channel form in bedrock rivers

Erosion processes in bedrock‐floored rivers shape channel cross‐sectional geometry and the broader landscape. However, the influence of weathering on channel slope and geometry is not well understood. Weathering can produce variation in rock erodibility within channel cross‐sections. Recent numerical modeling results suggest that weathering may preferentially weaken rock on channel banks relative to the thalweg, strongly influencing channel form. Here, we present the first quantitative field study of differential weathering across channel cross‐sections. We hypothesize that average cross‐section erosion rate controls the magnitude of this contrast in weathering between the banks and the thalweg. Erosion rate, in turn, is moderated by the extent to which weathering processes increase bedrock erodibility. We test these hypotheses on tributaries to the Potomac River, Virginia, with inferred erosion rates from ~0.1 m/kyr to >0.8 m/kyr, with higher rates in knickpoints spawned by the migratory Great Falls knickzone. We selected nine channel cross‐sections on three tributaries spanning the full range of erosion rates, and at multiple flow heights we measured (1) rock compressive strength using a Schmidt hammer, (2) rock surface roughness using a contour gage combined with automated photograph analysis, and (3) crack density (crack length/area) at three cross‐sections on one channel. All cross‐sections showed significant (p < 0.01 for strength, p < 0.05 for roughness) increases in weathering by at least one metric with height above the thalweg. These results, assuming that the weathered state of rock is a proxy for erodibility, indicate that rock erodibility varies inversely with bedrock inundation frequency. Differences in weathering between the thalweg and the channel margins tend to decrease as inferred erosion rates increase, leading to variations in channel form related to the interplay of weathering and erosion rate. This observation is consistent with numerical modeling that predicts a strong influence of weathering‐related erodibility on channel morphology. Copyright © 2017 John Wiley & Sons, Ltd.     

Substructural first‐ and second‐order model identification for structural damage assessment

This paper presents two methods to perform system identification at the substructural level, taking advantage of reduction in the number of unknowns and degrees of freedom (DOFs) involved, for damage assessment of fairly large structures. The first method is based on first‐order state space formulation of the substructure where the eigensystem realization algorithm (ERA) and the observer/Kalman filter identification (OKID) are used. Identification at the global level is then performed to obtain the second‐order model parameters. In the second method, identification is performed at the substructural level in both the first‐ and second‐order model identification. Both methods are illustrated using numerical simulation studies where results indicate their significantly better performance than identification using the global structure, in terms of efficiency and accuracy. A 12‐DOF system and a fairly large structural system with 50 DOFs are used where the effects of noisy data are considered. In addition to numerical simulation studies, laboratory experiments involving an eight‐storey frame model are carried out to illustrate the performance of the proposed method. The identification results presented in terms of the stiffness integrity index show that the proposed methodology is able to locate and quantify damage fairly accurately. Copyright © 2005 John Wiley & Sons, Ltd.     

Obtaining refined first‐order predictive models of linear structural systems

This study presents an effective method for identifying predictive models and the underlying modal parameters of linear structural systems using only measured output and excitation time histories obtained from dynamic testing. The system under examination is modelled as a first‐order multi‐input multi‐output time‐invariant system, and the structural model is realized using the Eigensystem Realization Algorithm together with the Observer/Kalman filter IDentification algorithm. The identified state‐space model is further refined using a non‐linear optimization technique based on sequential quadratic programming. The numerical examples show that the developed methodology performs very well even in the presence of inadequate instrumentation and measurement noise, and that the methodology is highly capable of creating realistic predictive models of structural systems, as well as estimating their underlying modal parameters. Copyright © 2002 John Wiley & Sons, Ltd.     

Assessing the effect of vegetation‐related bank strength on channel morphology and stability in gravel‐bed streams using numerical models

Bank strength due to vegetation dominates the geometry of small stream channels, but has virtually no effect on the geometry of larger ones. The dependence of bank strength on channel scale affects the form of downstream hydraulic geometry relations and the meandering‐braiding threshold. It is also associated with a lateral migration threshold discharge, below which channels do not migrate appreciably across their floodplains. A rational regime model is used to explore these scale effects: it parameterizes vegetation‐related bank strength using a dimensionless effective cohesion, C_r*. The scale effects are explored primarily using an alluvial state space defined by the dimensionless formative discharge, Q*, and channel slope, S, which is analogous to the Q–S diagrams originally used to explore meandering‐braiding thresholds. The analyses show that the effect of vegetation on both downstream hydraulic geometry and the meandering‐braiding threshold is strongest for the smallest streams in a watershed, but that the effect disappears for Q* > 10⁶. The analysis of the migration threshold suggests that the critical discharge ranges from about 5 m³/s to 50 m³/s, depending on the characteristic rooting depth for the vegetation. The analysis also suggests that, where fires frequently affect riparian forests, channels may alternate between laterally stable gravel plane‐bed channels and laterally active riffle‐pool channels. These channels likely do not exhibit the classic dynamic equilibrium associated with alluvial streams, but instead exhibit a cyclical morphologic evolution, oscillating between laterally stable and laterally unstable end‐members with a frequency determined by the forest fire recurrence interval. Copyright © 2008 John Wiley & Sons, Ltd.     

Time‐integrated sampling of fluvial suspended sediment: a simple methodology for small catchments

Fine‐grained (<62·5 µm) suspended sediment transport is a key component of the geochemical flux in most fluvial systems. The highly episodic nature of suspended sediment transport imposes a significant constraint on the design of sampling strategies aimed at characterizing the biogeochemical properties of such sediment. A simple sediment sampler, utilizing ambient flow to induce sedimentation by settling, is described. The sampler can be deployed unattended in small streams to collect time‐integrated suspended sediment samples. In laboratory tests involving chemically dispersed sediment, the sampler collected a maximum of 71% of the input sample mass. However, under natural conditions, the existence of composite particles or flocs can be expected to increase significantly the trapping efficiency. Field trials confirmed that the particle size composition and total carbon content of the sediment collected by the sampler were representative statistically of the ambient suspended sediment. Copyright © 2000 John Wiley & Sons, Ltd.     

Three‐dimensional flow structure and channel change in an asymmetrical compound meander loop,Embarras River,Illinois

The planform dynamics of meandering rivers produce a complex array of meander forms, including elongated meander loops. Thus far, few studies have examined in detail the flow structure within meander loops and the relation of flow structure to patterns of planform change. This field‐based investigation examines relations between three‐dimensional fluid motion and channel change within an elongated, asymmetrical meander loop containing multiple pool–riffle structures. The downstream velocity field is characterized by a high‐velocity core that shifts slightly outward as flow moves through individual lobes of the loop. For some of the measured flows this core becomes submerged below the water surface downstream of the lobe apexes. Vectors of cross‐stream/vertical velocities indicate that skew‐induced helical motion develops within the pools near lobe apexes and decays over riffles where channel curvature is less pronounced. Maximum rates of bank retreat generally occur near lobe apexes where impingement of the flow on the outer channel bank is greatest. However, maximum rates and loci of bank retreat differ for upstream and downstream lobes of the loop, leading to increasing asymmetry of loop geometry over time—a finding consistent with experimental investigations of loop evolution. Copyright © 2003 John Wiley & Sons, Ltd.     

A conceptual model for meander initiation in bedload‐dominated streams

A simple analytic model is presented relating local sediment transport capacity to variance in the transverse shear stress distribution in a stream channel. The model is used to develop a physically based conceptual model for the initiation of meandering in straight, bedload‐dominated streams as a result of a feedback mechanism. The feedback maximizes the cross‐sectional shear stress variance and – in order to achieve stability – ultimately minimizes the energy slope at repeated locations along the channel, subject to steady‐state mass flux and the stability of the channel boundary. These locations develop into pools in a fully developed meandering channel; they represent attractor states wherein sediment continuity is satisfied using the least possible energy expenditure per unit length of channel. However, since the cross‐sectional geometry of a pool (and the adjacent bar) is asymmetric, these attractor states are only conditionally stable, requiring strong, curvature‐induced secondary circulation to maintain their asymmetry. Between two successive pools, a stream occupies a metastable, higher energy state (corresponding to a riffle) that requires greater energy expenditure per unit length of channel to transport the same volume of sediment. The model we present links processes at the scale of a channel width to adjustments of the channel sinuosity and slope at the scale of a channel reach. We argue that the reach‐scale extremal hypotheses employed by rational regime models are mathematical formalisms that permit a one‐dimensional theory to describe the three‐dimensional dynamics producing stream morphology. Our model is consistent with the results from stream table experiments, with respect to both the rate of development of meandering and the characteristics of the equilibrium channel morphology. Copyright © 2006 John Wiley & Sons, Ltd.     

Channel‐planform evolution in four rivers of Olympic National Park,Washington, USA: the roles of physical drivers and trophic cascades

Identifying the relative contributions of physical and ecological processes to channel evolution remains a substantial challenge in fluvial geomorphology. We use a 74‐year aerial photographic record of the Hoh, Queets, Quinault, and Elwha Rivers, Olympic National Park, Washington, USA, to investigate whether physical or trophic‐cascade‐driven ecological factors – excessive elk impacts after wolves were extirpated a century ago – are the dominant drivers of channel planform in these gravel‐bed rivers. We find that channel width and braiding show strong relationships with recent flood history. All four rivers widened significantly after having been relatively narrow in the 1970s, consistent with increased flood activity since then. Channel planform also reflects sediment‐supply changes, evident from landslide response on the Elwha River. We surmise that the Hoh River, which shows a multi‐decadal trend toward greater braiding, is adjusting to increased sediment supply associated with rapid glacial retreat. These rivers demonstrate transmission of climatic signals through relatively short sediment‐routing systems that lack substantial buffering by sediment storage. Legacy effects of anthropogenic modification likely also affect the Quinault River planform. We infer no correspondence between channel evolution and elk abundance, suggesting that trophic‐cascade effects in this setting are subsidiary to physical controls on channel morphology. Our findings differ from previous interpretations of Olympic National Park fluvial dynamics and contrast with the classic example of Yellowstone National Park, where legacy effects of elk overuse are apparent in channel morphology; we attribute these differences to hydrologic regime and large‐wood availability. Published 2016. This article is a U.S. Government work and is in the public domain in the USA     

Linearization and first‐order expansion of the rocking motion of rigid blocks stepping on viscoelastic foundation

In structural mechanics there are several occasions where a linearized formulation of the original non‐linear problem reduces considerably the computational effort for the response analysis. In a broader sense, a linearized formulation can be viewed as a first‐order expansion of the dynamic equilibrium of the system about a ‘static’ configuration; yet caution should be exercised when identifying the ‘correct’ static configuration. This paper uses as a case study the rocking response of a rigid block stepping on viscoelastic supports, whose non‐linear dynamics is the subject of the companion paper, and elaborates on the challenge of identifying the most appropriate static configuration around which a first‐order expansion will produce the most dependable results in each regime of motion. For the regime when the heel of the block separates, a revised set of linearized equations is presented, which is an improvement to the unconservative equations published previously in the literature. The associated eigenvalues demonstrate that the characteristics of the foundation do not affect the rocking motion of the block once the heel separates. Copyright © 2008 John Wiley & Sons, Ltd.     

A response‐based decoupling criterion for multiply‐supported secondary systems

It is often infeasible to carry out coupled analyses of multiply‐supported secondary systems for earthquake excitations. ‘Approximate’ decoupled analyses are then resorted to, unless the response errors due to those are significantly high. This study proposes a decoupling criterion to identify such cases where these errors are likely to be larger than an acceptable level. The proposed criterion is based on the errors in the primary system response due to decoupling and has been obtained by assuming (i) the input excitation to be an ideal white noise process, (ii) cross‐modal correlation to be negligible, and (iii) the combined system to be classically damped. It uses the modal properties of the undamped combined system, and therefore, a perturbation approach has been formulated to determine the combined system properties in case of light to moderately heavy secondary systems. A numerical study has been carried out to illustrate the accuracy achieved with the proposed perturbation formulation. The proposed decoupling criterion has been validated with the help of two example primary‐secondary systems and four example excitation processes. Copyright © 2002 John Wiley & Sons, Ltd.     

Fine‐grained bed sediment storage within the main channel systems of the Frome and Piddle catchments,Dorset, UK

Lowland permeable catchments in the UK are particularly prone to sedimentation problems, on account of the increased fine sediment loadings generated by recent land‐use change and their stable seasonal hydrological regimes, which are frequently depleted by groundwater abstraction. Fine‐grained sediment storage on the bed of the main channel systems of the Frome (437 km²) and Piddle (183 km²) catchments, Dorset, UK, has been examined at 29 sites using a sediment remobilization technique. Measurements encompassed the period February 2003–July 2004. At individual sites in the Frome, average values ranged between 410 and 2630 g m^?2, with an overall mean of 918 g m^?2. In the Piddle, the average values for individual sites varied between 260 and 4340 g m^?2, with an overall mean of 1580 g m^?2. Temporal variations in fine bed sediment storage at each site were appreciable, with the coefficients of variation ranging between 43 and 155% in the Frome and between 33 and 160% in the Piddle. Average reach‐scale specific bed sediment storage increased markedly downstream along each main stem from 2 to 29 t km^?1 (Frome) and from 4 to 19 t km^?1 (Piddle). Total fine sediment storage on the channel bed of the Frome varied between 479 t (5 t km^?1) and 1694 t (17 t km^?1), with a mean of 795 t (7 t km^?1), compared with between 371 t (5 t km^?1) and 1238 t (14 t km^?1) with a mean of 730 t (9 t km^?1) in the Piddle. During the study period, fine bed sediment storage was typically equivalent to 18% (Frome) and 57% (Piddle) of the mean annual suspended sediment flux at the study catchment outlets. Copyright © 2006 John Wiley & Sons, Ltd.     

Measuring channel planform change from image time series: A generalizable,spatially distributed,probabilistic method for quantifying uncertainty

Channels change in response to natural or anthropogenic fluctuations in streamflow and/or sediment supply and measurements of channel change are critical to many river management applications. Whereas repeated field surveys are costly and time-consuming, remote sensing can be used to detect channel change at multiple temporal and spatial scales. Repeat images have been widely used to measure long-term channel change, but these measurements are only significant if the magnitude of change exceeds the uncertainty. Existing methods for characterizing uncertainty have two important limitations. First, while the use of a spatially variable image co-registration error avoids the assumption that errors are spatially uniform, this type of error, as originally formulated, can only be applied to linear channel adjustments, which provide less information on channel change than polygons of erosion and deposition. Second, previous methods use a level-of-detection (LoD) threshold to remove non-significant measurements, which is problematic because real changes that occurred but were smaller than the LoD threshold would be removed. In this study, we present a new method of quantifying uncertainty associated with channel change based on probabilistic, spatially varying estimates of co-registration error and digitization uncertainty that obviates a LoD threshold. The spatially distributed probabilistic (SDP) method can be applied to both linear channel adjustments and polygons of erosion and deposition, making this the first uncertainty method generalizable to all metrics of channel change. Using a case study from the Yampa River, Colorado, we show that the SDP method reduced the magnitude of uncertainty and enabled us to detect smaller channel changes as significant. Additionally, the distributional information provided by the SDP method allowed us to report the magnitude of channel change with an appropriate level of confidence in cases where a simple LoD approach yielded an indeterminate result. © 2020 John Wiley & Sons, Ltd.     

A hillslope hydrology approach for catchment‐scale slope stability analysis

Regional analysis of slope stability is often constrained by availability of data. Model requirements for input data cannot be met at the desired spatial resolution because data are either site‐speci?c or non‐existent. Faced with these dif?culties it has often been the practice to assume that certain parameters are uniform throughout the area of interest. An alternative approach proposed here allows a more detailed discrimination of slope stability conditions. Based on the principles of hillslope hydrology, hydrologic information can be generated at suf?cient resolution to allow higher resolution slope stability analysis. Measurements from an instrumented network in a small area have been used to establish index‐based models for topographic and climate‐related controls of piezometric response. The ability to relate groundwater levels to rainfall and topographic parameters provides a means of up‐scaling to larger catchments and ultimately the opportunity to generate a catchment‐wide prediction of the distribution, magnitude and frequency of rainstorm‐generated groundwater levels. The example provided in this study uses the topography index of TOPMODEL in GIS to predict the spatial patterns of groundwater elevation for seasonal soil moisture conditions and given rainfall inputs. This allows modelling of catchment‐wide response of soil water to rainstorms with different return periods (representing different magnitudes), and is an essential prerequisite for a probabilistic regional slope stability analysis. Copyright © 2004 John Wiley & Sons, Ltd.     

The impact of exceptional events on erosion,bedload transport and channel stability in a step‐pool channel

Sediment transport in the Erlenbach, a small stream with step‐pool morphology in the canton of Schwyz, Switzerland, has been monitored for more than 20 years. During this time three exceptional events (events with high sediment yield and long return times that have a large effect on channel morphology) have impacted the stream and partly or completely rearranged the existing step‐pool morphology. In the aftermath of the events, sediment transport rates at a given discharge and total sediment yield remained elevated for about a year or longer. For the last event, dated on the 20 June 2007, observations of boulder mobility and step destruction were used to interpret channel stability. Boulders with median diameters of up to 135 cm and estimated weights of more than 2·5 tons have moved during the 2007 event. Using hydraulic observations and shear stress calculations boulders up to 65 cm in diameter were predicted to have been fully mobile in peak conditions, even if form resistance and increased critical stresses needed for the initiation of motion in steep streams were taken into account. For two of the events, estimated peak shear stresses at the bed exceeded 1000 Pa, calculated both from observations of the flow hydraulics and from boulder mobility. This suggests that highly energetic flows occur relatively frequently in small, steep streams and that large boulders can be transported by fluvial processes in such streams. The observations have potential significance for hazard risk mitigation, stream engineering and restoration. Copyright © 2009 John Wiley & Sons, Ltd.     

A real‐space cellular automaton laboratory

Geomorphic investigations may benefit from computer modelling approaches that rely entirely on self‐organization principles. In the vast majority of numerical models, instead, points in space are characterized by a variety of physical variables (e.g. sediment transport rate, velocity, temperature) recalculated over time according to some predetermined set of laws. However, there is not always a satisfactory theoretical framework from which we can quantify the overall dynamics of the system. For these reasons, we prefer to concentrate on interaction patterns using a basic cellular automaton modelling framework. Here we present the Real‐Space Cellular Automaton Laboratory (ReSCAL), a powerful and versatile generator of 3D stochastic models. The objective of this software suite, released under a GNU licence, is to develop interdisciplinary research collaboration to investigate the dynamics of complex systems. The models in ReSCAL are essentially constructed from a small number of discrete states distributed on a cellular grid. An elementary cell is a real‐space representation of the physical environment and pairs of nearest‐neighbour cells are called doublets. Each individual physical process is associated with a set of doublet transitions and characteristic transition rates. Using a modular approach, we can simulate and combine a wide range of physical processes. We then describe different ingredients of ReSCAL leading to applications in geomorphology: dune morphodynamics and landscape evolution. We also discuss how ReSCAL can be applied and developed across many disciplines in natural and human sciences. Copyright © 2013 John Wiley & Sons, Ltd.     

Enhanced application of root‐reinforcement algorithms for bank‐stability modeling

Riparian vegetation is known to exert a number of mechanical and hydrologic controls on bank stability. In particular, plant roots provide mechanical reinforcement to a soil matrix due to the different responses of soils and roots to stress. Root reinforcement is largely a function of the strength of the roots crossing potential shear planes, and the number and diameter of such roots. However, previous bank stability models have been constrained by limited field data pertaining to the spatial and temporal variability of root networks within stream banks. In this paper, a method is developed to use root‐architecture data to derive parameters required for modeling temporal and spatial changes in root reinforcement. Changes in root numbers over time were assumed to follow a sigmoidal curve, which commonly represents the growth rates of organisms. Regressions for numbers of roots crossing potential shear planes over time showed small variations between species during the juvenile growth phase, but extrapolation led to large variations in root numbers by the time the senescent phase of the sigmoidal growth curve had been reached. In light of potential variability in the field data, the mean number of roots crossing a potential shear plane at each year of tree growth was also calculated using data from all species and an additional sigmoidal regression was run. After 30 years the mean number of roots predicted to cross a 1 m shear plane was 484, compared with species‐specific curves whose values ranged from 240 roots for black willow trees to 890 roots for western cottonwood trees. In addition, the effect of spatial variations in rooting density with depth on stream‐bank stability was modeled using the bank stability and toe erosion model (BSTEM). Three root distributions, all approximating the same average root reinforcement (5 kPa) over the top 1 m of the bank profile, were modeled, but with differing vertical distributions (concentrated near surface, non‐linear decline with depth, uniform over top meter). It was found that stream‐bank F_S varied the most when the proportion of the failure plane length to the depth of the rooting zone was greatest. Copyright © 2008 John Wiley & Sons, Ltd.     

Rainfall thresholds for shallow landsliding derived from pressure‐head monitoring: cases with permeable and impermeable bedrocks in Boso Peninsula,Japan

Rainfall thresholds for shallow landslide initiation were determined for hillslopes with two types of bedrock, permeable sandstone and impermeable mudstone, in the Boso Peninsula, Japan. The pressure‐head response to rainfall was monitored above a slip scarp due to earlier landslides. Multiple regression analysis estimated the rainfall thresholds for landsliding from the relation between the magnitude of the rainfall event and slope instability caused by the increased pressure heads. The thresholds were expressed as critical combinations of rainfall intensity and duration, incorporating the geotechnical properties of the hillslope materials and also the slope hydrological processes. The permeable sandstone hillslope has a greater critical rainfall and hence a longer recurrence interval than the impermeable mudstone hillslope. This implies a lower potential for landsliding in sandstone hillslopes, corresponding to lower landslide activity. Copyright © 2007 John Wiley & Sons, Ltd.     

Capturing the spatiotemporal variability of bedload transport: A time‐lapse imagery technique

Monitoring sediment transport in morphologically complex and labile channels remains a difficult task, even at the laboratory scale. To address this challenge, a fully automated imagery technique for continuously mapping the spatial and temporal variability of bedload transport is proposed. This method uses differentiated time‐lapse imagery taken from a fixed camera to detect bed variations induced by grain displacement. The technique is not based on tracking the individual particles; rather, it evaluates macroscopic colour changes within a region that contains several grains, which depend on the occurrence and intensity of the bedload transport. Image‐derived data were compared with the sediment flux measured during four flume experiments, and produced good correspondence. The method provides continuous tracking of the location of the transporting channels, and enables estimation of local variations in the magnitude of the bedload flux. Moreover, the spatial extent of the monitoring area offers an unprecedented opportunity to aggregate spatially dense and continuous data at the reach scale, as needed to properly capture the full range of variability of morphologically complex and rapidly evolving gravel‐bed rivers. Despite being limited to laboratory‐scale physical experiments, the method provides useful data to investigate fundamental morphodynamic processes such as bar migration, bank erosion, anabranches opening/closure, and the associate spatial and temporal scales. Further, the data obtained have the potential to enhance numerical model calibration and improve our understanding of the complex dynamics of real‐world settings. Copyright © 2017 John Wiley & Sons, Ltd.     

A coupled one‐dimensional viscoelastic–plastic model for aquitard consolidation caused by hydraulic head variations in aquifers

Accurate and practical calculation of aquitard consolidation is required for a reliable analysis of land subsidence caused by groundwater overexploitation in a multilayered aquifer system because aquitards are generally more compressible than aquifers are. This study proposes a coupled one‐dimensional viscoelastic–plastic consolidation model that considers the combined effect of changes in soil parameters and body force to simulate aquitard consolidation caused by hydraulic head variations in neighbouring aquifers. The proposed model uses variable total stress and simultaneously solves hydraulic head and vertical soil displacement. The constitutive relation based on the Voigt model with different elastic moduli of the spring in normally consolidated and overconsolidated soils is used to describe the viscoelastic–plastic deformation mechanism of aquitards. In addition, the proposed model considers the combined effect of variations in hydraulic conductivity, elastic moduli, and body force on the calculation of aquitard consolidation. Three hypothetical scenarios with various hydraulic head variations in aquifers are used to examine the coupled one‐dimensional viscoelastic–plastic consolidation model. The results show that neglecting plasticity and viscosity of soil causes aquitard consolidation to be respectively underestimated and overestimated. In addition, ignoring body force variation underestimates aquitard consolidation, whereas neglecting soil parameters variation overestimates aquitard consolidation. Two real case scenarios are also studied to further demonstrate the applicability of the coupled one‐dimensional viscoelastic–plastic consolidation model. Copyright © 2015 John Wiley & Sons, Ltd.