An efficient numerical model for multicomponent compressible flow in fractured porous media

An efficient and accurate numerical model for multicomponent compressible single-phase flow in fractured media is presented. The discrete-fracture approach is used to model the fractures where the fracture entities are described explicitly in the computational domain. We use the concept of cross flow equilibrium in the fractures. This will allow large matrix elements in the neighborhood of the fractures and considerable speed up of the algorithm. We use an implicit finite volume (FV) scheme to solve the species mass balance equation in the fractures. This step avoids the use of Courant–Freidricks–Levy (CFL) condition and contributes to significant speed up of the code. The hybrid mixed finite element method (MFE) is used to solve for the velocity in both the matrix and the fractures coupled with the discontinuous Galerkin (DG) method to solve the species transport equations in the matrix. Four numerical examples are presented to demonstrate the robustness and efficiency of the proposed model. We show that the combination of the fracture cross-flow equilibrium and the implicit composition calculation in the fractures increase the computational speed 20–130 times in 2D. In 3D, one may expect even a higher computational efficiency.     

Understanding soil‐pipe flow and its role in ephemeral gully erosion

The role of soil pipeflow in ephemeral gully erosion is not well understood. Experiments were conducted on continuous soil pipes to better understand the role of internal erosion of soil pipes and its relation to ephemeral gully development. Soil beds of 140 cm length, 100 cm width and 20 cm depth had a single soil pipe of different initial sizes (2, 4, 6, 8, and 10 mm diameter) extend from a water reservoir to the outlet. Experiments were run on Providence silt loam and Smithdale loam soils under a constant head of 15 cm established for 30 min. Either the tunnel collapsed or the head could not be maintained. Soil pipes that were initially 2 and 4 mm clogged instantaneously at their mouth and did not exhibit flow, whereas, pipes initially ≥ 6 mm enlarged by 268, 397, and 699% on average for the 6, 8, and 10 mm diameters, respectively. Critical shear stress values were found to be essentially zero, and erodibility values gave erosion indexes that were extremely high. The rapid internal erosion resulted in erratic flow and sediment concentrations with periods of no flow as pipes were temporarily clogged followed by surges of high flow and high sediment concentrations. Tensiometers within 6 cm of the soil pipes did not exhibit pressure increases typically associated with pipe clogging. Flow through 10 mm diameter soil pipes exhibited tunnel collapse for both soils tested. Tunnel collapse typically occurred within minutes of flow establishment suggesting that ephemeral gullies could be misinterpreted as being caused by convergent surface flow if observations were made after the runoff event instead of when flow is first established through soil pipes. Copyright © 2011 John Wiley & Sons, Ltd.     

Modeling granular soils liquefaction using coupled lattice Boltzmann method and discrete element method

In this paper, a novel coupled pore-scale model of pore-fluid interacting with discrete particles is presented for modeling liquefaction of saturated granular soil. A microscale idealization of the solid phase is achieved using the discrete element method (DEM) while the fluid phase is modeled at a pore-scale using the lattice Boltzmann method (LBM). The fluid forces applied on the particles are calculated based on the momentum exchange between the fluid and particles. The presented model is based on a first principles formulation in which pore-pressure develops due to actual changes in pore space as particles׳ rearrangement occurs during shaking. The proposed approach is used to model the response of a saturated soil deposit subjected to low and large amplitude seismic excitations. Results of conducted simulations show that at low amplitude shaking, the input motion propagates following the theory of wave propagation in elastic solids. The deposit response to the strong input motion indicates that liquefaction took place and it was due to reduction in void space during shaking that led to buildup in pore-fluid pressure. Soil liquefaction was associated with soil stiffness degradation and significant loss of interparticle contacts. Simulation results also indicate that the level of shaking-induced shear strains and associated volumetric strains play a major role in the onset of liquefaction and the rate of pore-pressure buildup.     

Multifluid flow in bedded porous media: laboratory experiments and numerical simulations

Understanding light nonaqueous-phase liquid (LNAPL) movement in heterogeneous vadose environments is important for effective remediation design. We investigated LNAPL movement near a sloping fine- over coarse-grained textural interface, forming a capillary barrier. LNAPL flow experiments were performed in a glass chamber (50 cm×60 cm×1.0 cm) using two silica sands (12/20 and 30/40 sieve sizes). Variable water saturations near the textural interface were generated by applying water uniformly to the sand surface at various flow rates. A model LNAPL (Soltrol® 220) was subsequently released at two locations at the sand surface. Visible light transmission was used to quantitatively determine water saturations prior to LNAPL release and to observe LNAPL flow paths. Numerical simulations were performed using the Subsurface Transport Over Multiple Phases (STOMP) simulator, employing two nonhysteretic relative permeability–saturation–pressure (k–S–P) models. LNAPL movement strongly depended on the water saturation in the fine-grained sand layer above the textural interface. In general, reasonable agreement was found between observed and predicted water saturations near the textural interface and LNAPL flow paths. Discrepancies between predictions based on the van Genuchten/Mualem (VGM) and Brooks–Corey/Burdine (BCB) k–S–P models existed in the migration speed of the simulated LNAPL plume and the LNAPL flow patterns at high water saturation above the textural interface. In both instances, predictions based on the BCB model agreed better with experimental observations than predictions based on the VGM model. The results confirm the critical role water saturation plays in determining LNAPL movement in heterogeneous vadose zone environments and that accurate prediction of LNAPL flow paths depends on the careful selection of an appropriate k–S–P model.     

A modified nearly analytic discrete method and wavefield simulations in transversely isotropic media

Under the general case, rocks under the ground can be approximately considered as an elastic medium. Elastic wave equation is a partial differential equation, which describes the elastic wave propagation in elastic media. Simple elastic wave equation can be solved analytically, however most wave equations are very complex, which can only be solved by using numerical methods. Numerical simulations of seismic wave fields have become an important method for studying seismic wave propagation in co…     

求解3D弹性波方程的并行WNAD方法及其TI介质中的波场模拟

准确模拟TTI介质中弹性波的传播是研究地震各向异性、AVO反演的基础. 在二维加权近似解析离散化(WNAD)算法的基础上, 本文发展的并行WNAD算法是一种研究三维横向各向同性(TI)介质中弹性波传播的、快速高效的数值模拟方法. 我们首先介绍三维WNAD方法的构造过程, 然后与经典的差分格式——交错网格(SG)算法进行了比较. 理论分析和数值算例表明, WNAD算法比交错网格算法更适合在高性能计算机上进行大规模弹性波场模拟. 同时, 本文利用并行的WNAD方法研究了弹性波在TTI介质中的传播规律, 观测了TI介质中弹性波传播的重要特征:横波分离、体波耦合和速度各向异性等. 在TTI介质分界面处, 弹性波产生更加复杂的折射、反射和波型转化, 使得波场非常复杂, 研究和辨别不同类型的波能够加深我们对由裂隙诱导的各向异性介质的认识.     

An efficient numerical model for incompressible two-phase flow in fractured media

Various numerical methods have been used in the literature to simulate single and multiphase flow in fractured media. A promising approach is the use of the discrete-fracture model where the fracture entities in the permeable media are described explicitly in the computational grid. In this work, we present a critical review of the main conventional methods for multiphase flow in fractured media including the finite difference (FD), finite volume (FV), and finite element (FE) methods, that are coupled with the discrete-fracture model. All the conventional methods have inherent limitations in accuracy and applications. The FD method, for example, is restricted to horizontal and vertical fractures. The accuracy of the vertex-centered FV method depends on the size of the matrix gridcells next to the fractures; for an acceptable accuracy the matrix gridcells next to the fractures should be small. The FE method cannot describe properly the saturation discontinuity at the matrix–fracture interface. In this work, we introduce a new approach that is free from the limitations of the conventional methods. Our proposed approach is applicable in 2D and 3D unstructured griddings with low mesh orientation effect; it captures the saturation discontinuity from the contrast in capillary pressure between the rock matrix and fractures. The matrix–fracture and fracture–fracture fluxes are calculated based on powerful features of the mixed finite element (MFE) method which provides, in addition to the gridcell pressures, the pressures at the gridcell interfaces and can readily model the pressure discontinuities at impermeable faults in a simple way. To reduce the numerical dispersion, we use the discontinuous Galerkin (DG) method to approximate the saturation equation. We take advantage of a hybrid time scheme to alleviate the restrictions on the size of the time step in the fracture network. Several numerical examples in 2D and 3D demonstrate the robustness of the proposed model. Results show the significance of capillary pressure and orders of magnitude increase in computational speed compared to previous works.     

The effect of instream logs on river-bank erosion: Field measurements of hydraulics and erosion rates

This is the first substantial field measurement of river-bank erosion around fallen logs in rivers. Whilst numerous studies have established that living trees can stabilize river banks, and that fallen trees can cause scour of the river bed, knowledge of bank erosion effects from logs is largely restricted to qualitative observations. Recent flume studies suggest that a single log can increase near-bank velocity (and thus erosion) and this increase is related to the blockage ratio of the log and the distance between the log and bank. However, hydraulic interactions between logs can reduce this increase or even decrease the near-bank velocity. These theories, developed in a straight flume, have not been tested in the field. We measured erosion rates (relative to controls) on river banks adjacent to 35 large logs for 2 years, and velocity distributions around 11 logs during a near-bankfull flow in anabranching channels of the River Murray, SE Australia. These channels have abundant large instream logs, consistent bank material, and consistent regulated high flows. The field results generally supported the velocity changes caused by single and multiple logs in the flume studies, with single logs increasing near-bank velocity, but with the hydraulic interactions between successive logs tending to reduce this increase. Flow patterns caused by logs adjacent to curved banks were more complicated as the local effects of logs reinforced or weakened recirculating flows. Instream logs did not change overall, average, bank erosion rates, but they tended to shift the erosion from bank top to bank toe. However, individual logs increased or decreased bank erosion rates in patterns that generally concur with the near-bank velocity changes predicted in flume studies: that isolated logs increased erosion rates whilst hydraulically interacting logs did not increase erosion rates. © 2020 John Wiley & Sons, Ltd.     

Multi-block lattice Boltzmann simulations of solute transport in shallow water flows

We report a two-dimensional multi-block lattice Boltzmann model for solute transport in shallow water flows, which is developed based on the advection–diffusion equation for mass transport and the shallow water equations for the flows. A weighting factor is included in the centered scheme for improved accuracy. The model is firstly verified by simulating three benchmark tests: wind-driven circulation in a dish-shaped lake, jet-forced flow in a circular basin, and flow formed by two parallel streams containing different uniform concentrations at the same constant velocity; and then it is applied to a practical wind-induced flow, Baiyangdian Lake, which is characterized by irregular geometries and complex bathymetries. The numerical results have shown that the model is able to produce accurate and detailed results for both water flows and solute transport, which is attractive, especially for flows in narrow zones of practical terrains and certain areas with largely varying pollutant concentrations.     

2D numerical simulations of earthquake ground motion: examples from the Marche Region,Italy

A detailed numerical simulation of the ground motion and a site response analysis for two towns in the Marche Region (Treia and Cagli) is carried out on the basis of structural models deduced from available geological and geophysical data. In both cases, the reference event is an M = 5.7 earthquake associated with a normal fault located beneath each town. The ground motion is computed using the 2D spectral element method (SPEM 2D). The method solves the propagation of the seismic field through complex geological structures and enables an estimate of the effects of deep crustal structure, superficial geology, and topography on ground motion. Numerical simulations of the seismic field are performed along 2D vertical planes containing the seismic source. Strong ground motion has not been yet recorded in the two towns; therefore, the numerical simulation of ground motion represents a way to overcome the lack of instrumental data. The simulations carried out for Treia show that ground motion is influenced by both source mechanism and effects due to propagation through the geological structure, while ground motion in Cagli features strong local effects, caused by the presence of alluvial deposits under a large area of the town.     

Laboratory soil piping and internal erosion experiments: evaluation of a soil piping model for low‐compacted soils

Mechanistic models have been proposed for soil piping and internal erosion on well‐compacted levees and dams, but limited research has evaluated these models in less compacted (more erodible) soils typical of hillslopes and streambanks. This study utilized a soil box (50 cm long, 50 cm wide and 20 cm tall) to conduct constant‐head, soil pipe and internal erosion experiments for two soils (clay loam from Dry Creek and sandy loam from Cow Creek streambanks) packed at uniform bulk densities. Initial gravimetric moisture contents prior to packing were 10, 12 and 14% for Dry Creek soil and 8, 12, and 14% for Cow Creek soil. A 1‐cm diameter rod was placed horizontally along the length of the soil bed during packing and carefully removed after packing to create a continuous soil pipe. A constant head was maintained at the inflow end. Flow rates and sediment concentrations were measured from the pipe outlet. Replicate submerged jet erosion tests (JETs) were conducted to derive erodibility parameters for repacked samples at the same moisture contents. Flow rates from the box experiments were used to calibrate the mechanistic model. The influence of the initial moisture content was apparent, with some pipes (8% moisture content) expanding so fast that limited data was collected. The mechanistic model was able to estimate equivalent flow rates to those observed in the experiments, but had difficulty matching observed sediment concentrations when the pipes rapidly expanded. The JETs predicted similar erodibility coefficients compared to the mechanistic model for the more erodible cases but not for the less erodible cases (14% moisture content). Improved models are needed that better define the changing soil pipe cross‐section during supply‐ and transport‐limited internal erosion, especially for piping through lower compacted (more erodible) soils as opposed to more well‐compacted soils resulting from constructing levees and dams. Copyright © 2013 John Wiley & Sons, Ltd.     

Braiding process and bank erosion in the Brahmaputra River

The present work explores relations between stream power,braiding intensities and bank erosion in certain stretches of the Brahmaputra River.In this paper,an objective approach is presented to enable quantitative assessment of spatio-temporal behaviour of channel braiding process of the Brahmaputra River by using the Plan Form Index and corresponding estimation of stream power to establish a behavioural pattern of variability of potential energy expenditure.The braiding index is compared for discrete years to understand the morphological behaviour.Subsequently,a real time estimation of stream power for certain stretches of Brahmaputra River is done in order to analyse its variability in braiding intensity and bank erosion.The paper presents the dynamic behaviour of the channel pattern of the Brahmaputra River System in Assam valley of India over a time span of 18 years.The procedure addresses the selection of input parameters from digital satellite images,comprising scenes for the years 1990,1997 and 2007 with specific dates,from Dhubri near Indo-Bangladesh Border to Upper Assam.Deployment of GIS technique has been made to extract the required parameters to derive Plan Form Indices for the entire study reach.Stream power estimation is done for corresponding latest floods and for corresponding dates of image scenes.The study indicated that due to consistent aggradation of riverbed inducing temporal declination of stream power,there is an occurrence of wide spread braiding.This in turn incurs substantial yearly land loss due to bank erosion,caused by flow concentrations due to temporal evolution of multiple channels in the Brahmaputra River.     

Monitoring and numerical modelling of riverbank erosion processes: a case study along the Cecina River (central Italy)

Riverbank retreat along a bend of the Cecina River, Tuscany (central Italy) was monitored across a near annual cycle (autumn 2003 to summer 2004) with the aim of better understanding the factors influencing bank changes and processes at a seasonal scale. Seven flow events occurred during the period of investigation, with the largest having an estimated return period of about 1·5 years. Bank simulations were performed by linking hydrodynamic, fluvial erosion, groundwater flow and bank stability models, for the seven flow events, which are representative of the typical range of hydrographs that normally occur during an annual cycle. The simulations allowed identification of (i) the time of onset and cessation of mass failure and fluvial erosion episodes, (ii) the contributions to total bank retreat made by specific fluvial erosion and mass‐wasting processes, and (iii) the causes of retreat. The results show that the occurrence of bank erosion processes (fluvial erosion, slide failure, cantilever failure) and their relative dominance differ significantly for each event, depending on seasonal hydrological conditions and initial bank geometry. Due to the specific planimetric configuration of the study bend, which steers the core of high velocity fluid away from the bank at higher flow discharges, fluvial erosion tends to occur during particular phases of the hydrograph. As a result fluvial erosion is ineffective at higher peak discharges, and depends more on the duration of more moderate discharges. Slide failures appear to be closely related to the magnitude of peak river stages, typically occurring in close proximity to the peak phase (preferentially during the falling limb, but in some cases even before the peak), while cantilever failures more typically occur in the late phase of the flow hydrograph, when they may be induced by the cumulative effects of any fluvial erosion. Copyright © 2008 John Wiley & Sons, Ltd.     

Debris flow and slide deposits at the top of the Internal Liguride ophiolitic sequence, Northern Apennines, Italy: A record of frontal tectonic erosion in a fossil accretionary wedge

Abstract In the Northern Apennines, the Internal Liguride units are characterized by an ophiolite sequence that represents the stratigraphic base of a late Jurassic–early Paleocene sedimentary cover. The Bocco Shale represents the youngest deposit recognized in the sedimentary cover of the ophiolite and can be subdivided into two different groups of deep sea sediments. The first group is represented by slide, debris flow and high density turbidity current-derived deposits, whereas the second group consists of thin-bedded turbidites. Facies analysis and provenance studies indicate, for the former group, small and scarcely evoluted flows that rework an oceanic lithosphere and its sedimentary cover. We interpret the Bocco Shale as an ancient example of a deposit related to the frontal tectonic erosion of the accretionary wedge slope. The frontal tectonic erosion resulted in a large removal of materials, from the accretionary wedge front, that was reworked as debris flows and slide deposits sedimented on the lower plate above the trench deposits. The frontal tectonic erosion was probably connected with subduction of oceanic crust characterized by positive topographic relief. This interpretation can be also applied for the origin of analogous deposits of Western Alps and Corsica.     

The effect of wildfire on runoff and erosion in native Eucalyptus forest

Wildfires raise concerns over the risk of accelerated erosion as a result of increased overland flow and decreased protection of the soil by litter and ground vegetation cover. We investigated these issues following the 1994 fires that burnt large areas of native Eucalyptus forest surrounding Sydney, Australia. A review of previous studies identifies the fire and rainfall conditions that are likely to lead to increased runoff and accelerated erosion. We then compare runoff and erosion between burnt and unburnt sites for 10 months after the 1994 fires. At the scale of hillslope plots, the 1994 fire increased runoff by enhancing soil hydrophobicity, and greatly increased sediment transport, mainly through the reduced ground cover, which lowered substantially the threshold for initial sediment movement. However, both runoff and sediment transport were very localized, resulting in little runoff or sediment yield after the fire at the hillslope catchment scale. We identify that after moderately intense fires, rainfall events of greater than one year recurrence interval are required to generate substantial runoff and sediment yield. Such events did not occur during the monitoring period. Past work shows that mild burns have little effect on erosion, and it is only after the most extreme fires that erosion is produced from small, frequent storms. © 1998 John Wiley & Sons, Ltd.     

Soil erosion and overland flow in Japanese cypress plantations: spatio-temporal variations and a sampling strategy

ABSTRACT

Many studies have focused on soil erosion in unmanaged Japanese cypress plantations because the sparse understory vegetation and litter covering the forest ground enhance soil erosion. In this study, soil erosion, litter, and overland flow measurements were conducted over 14 months to identify the spatio-temporal variation and examine the optimal sample size. Fifteen traps (each 0.25 m wide) were installed in line along the bottom of a 15-m-wide slope. Soil erosion and overland flow had large spatial variations as compared to litter. The temporal coefficient of variation of soil erosion and overland flow was highest during dry seasons, while smaller during wet seasons. The random sampling analysis showed that the rate of decrease in spatio-temporal variation became moderate as the sample size increased beyond six. This result indicated that the optimal sample size was five, the total width of which was equivalent to about 8% of the monitored slope width.     

Rapid assessment of field erosion and sediment transport pathways in cultivated catchments after heavy rainfall events

This paper explores a scale‐adapted erosion mapping method which aims at a rapid assessment of field erosion and sediment transport pathways in catchments up to several square kilometres and compares the results with the output of a well‐known erosion model (LISEM). The mapping method is based on an event‐defined classification scheme of erosion intensity (zero, weak, moderate and strong) that is applied to arable fields, in combination with incision measurements of erosion features for each erosion intensity class on a small sample of fields. Sediment deposition is classified on the basis of quantity indicators and abundance. In addition, relevant conditions and erosion factors are determined for each field. The method was applied to an agricultural catchment (4·2 km²) in the Sundgau (Alsace), after a short but violent thunderstorm in May 2001, to illustrate its potential use and its limitations. The rainfall event led to strong erosion on the arable fields and a muddy flow that caused significant damage in the built‐up area. On the basis of the analyses of the incision measurements in combination with the mapping of erosion intensity classes, total erosion for the catchment was estimated as 15 000 t (an average of about 36 t[sol ]ha). Sediment deposition was found to occur in three major locations: (1) in thalwegs at the interface between maize and downslope winter wheat fields, (2) in downslope headlands where the flow direction suddenly changed due to oriented tillage structures in the perpendicular direction, and (3) the lowest corners of fields which collect all the runoff from the field. Preliminary data analyses suggest that erosion intensity is related to field size and[sol ]or tillage direction and to slope morphology. Model output (LISEM) appeared to depend more strongly on slope gradient than the results obtained with the mapping method. The method yields a database, which can be used as a foundation for conservation strategies in small regions with similar land use and geomorphology. The mapping and modelling methods are compared, and their complementary aspects are highlighted. Copyright © 2005 John Wiley & Sons, Ltd.     

A new scheme for numerical modelling of flow and transport processes in 3D fractured porous media

The objective of this work is to develop a new numerical approach for the three-dimensional modelling of flow and transient solute transport in fractured porous media which would provide an accurate and efficient treatment of 3D complex geometries and inhomogeneities. For this reason, and in order to eliminate as much as possible the number of degrees of freedom, the fracture network, fractures and their intersections, are solved with a coupled 2D–1D model while the porous matrix is solved independently with a 3D model. The interaction between both models is accounted for by a coupling iterative technique. In this way it is possible to improve efficiency and reduce CPU usage by avoiding 3D mesh refinements of the fractures. The approach is based on the discrete-fracture model in which the exact geometry and location of each fracture in the network must be provided as an input. The formulation is based on a multidimensional coupling of the boundary element method-multidomain (BEM-MD) scheme for the flow and boundary element dual reciprocity method-multidomain (BE-DRM-MD) scheme for the transport. Accurate results and high efficiency have been obtained and are reported in this paper.     

Meltwater erosion process of frozen soil as affected by thawed depth under concentrated flow in high altitude and cold regions

Changes in thawed depth of frozen soil caused by diurnal and seasonal temperature fluctuations are commonly found in high altitude and latitude regions of the world. These changes significantly influence hydrologic and erosion processes. Experimental data are necessary to improve the understanding and modeling of the phenomenon. Laboratory experiments were conducted in Beijing to assess the impacts of thawed soil depth, slope gradient, and flow rate on soil erosion by concentrated meltwater flow over an underlying frozen soil layer. Soil samples from watershed were filled in flumes, saturated before being frozen. After the soil was completely frozen, flumes were taken out of storage to thaw the frozen soil from top to the designed depths. Meltwater flow was simulated using a tank filled with water and icecubes at approximately 0°C. The erosion experiments involved four thawed soil depths of 1, 2, 5, and 10 cm; three slope gradients of 5°, 10°, and 15°; and three flow rates of 1, 2, and 4 l/min; and seven rill lengths of 0.5, 1, 2, 3, 4, 5, and 6 m. Sediment‐laden water samples were collected at the lower end of the flume for determination of sediment concentration. The results showed that sediment concentration increased exponentially with rill length to approach a maximum value. The sediment concentrations were closely correlated with thawed soil depth, flow rate, and slope gradient. Shallower thawed depths delivered more sediments than deeper thawed depths. Slope gradient was the primary factor responsible for severe erosion. The effect of flow rate on sediment concentration which decreased with increasing slope gradient, was not as significant as that of slope gradient. Results from these experiments are useful for understanding the effect of thawed soil depth on erosion process in thawed soils subject to freezing and for estimating erosion model parameters. Copyright © 2017 John Wiley & Sons, Ltd.     

Changes in the wind-wave field and related salt-marsh lateral erosion: inferences from the evolution of the Venice Lagoon in the last four centuries

Salt marshes are crucially important ecosystems at the boundary between the land and the sea, that are experiencing significant losses worldwide mainly dictated by the erosion of their margins. Improving our understanding of the mechanisms controlling marsh edge erosion is a key step to address conservation issues and salt-marsh response to changes in the environmental forcing. Here we have employed a complete, coupled Wind-Wave Tidal Model (WWTM) to analyse the temporal evolution of the wave field, and in particular of the mean wave-power density, in the Venice Lagoon over the past four centuries (from 1611 to 2012). We have then related wave-field changes to the observed erosion patterns determined by comparing recent aerial photographs (1978–2010) and historical bathymetric data. The results of our analyses from the Venice Lagoon show that, while wave-fields did not significantly change from 1611 to 1901, a rapid increase in wave-power densities occurred in the last century. This is suggested to depend on the positive feedback between relevant morphological evolutions and changes in the wave field, both influenced by natural forcing and anthropogenic pressures. We also emphasize the existence of a strong positive linear relationship between the volumetric marsh erosion rate and mean wave-power density. We thus suggest that relating salt-marsh lateral erosion rates to properly computed mean wave-power densities provides a valuable tool to address long-term tidal morphodynamics. © 2019 John Wiley & Sons, Ltd.