Geomodelling of a fluvial system with semi-supervised support vector regression

Fluvial deposits are a challenge for modelling flow in sub-surface reservoirs. Connectivity and continuity of permeable bodies have a major impact on fluid flow in porous media. Contemporary object-based and multipoint statistics methods face a problem of robust representation of connected structures. An alternative approach to model petrophysical properties is based on machine learning algorithm ? Support Vector Regression (SVR). Semi-supervised SVR is able to establish spatial connectivity taking into account the prior knowledge on natural similarities. SVR as a learning algorithm is robust to noise and captures dependencies from all available data. Semi-supervised SVR applied to a synthetic fluvial reservoir demonstrated robust results, which are well matched to the flow performance

Early diagenesis of bone and tooth apatite in fluvial and marine settings: Constraints from combined oxygen isotope, nitrogen and REE analysis

Fossil bones and teeth of Late Pleistocene terrestrial mammals from Rhine River gravels (RS) and the North Sea (NS), that have been exposed to chemically and isotopically distinct diagenetic fluids (fresh water versus seawater), were investigated to study the effects of early diagenesis on biogenic apatite. Changes in phosphate oxygen isotopic composition (delta O-18(PO4)), nitrogen content (wt.% N) and rare earth element (REE) concentrations were measured along profiles within bones that have not been completely fossilized, and in skeletal tissues (bone, dentine, enamel) with different susceptibilities to diagenetic alteration. Early diagenetic changes of elemental and isotopic compositions of apatite in fossil bone are related to the loss of the stabilizing collagen matrix. The REE concentration is negatively correlated with the nitrogen content, and therefore the amount of collagen provides a sensitive proxy for early diagenetic alteration. REE patterns of RS and NS bones indicate initial fossilization in a fresh water fluid with similar REE compositions. Bones from both settings have nearly collagen-free, REE-, U-, F- and Sr-enriched altered outer rims, while the collagen-bearing bone compacta in the central part often display early diagenetic pyrite void-fillings. However, NS bones exposed to Holocene seawater have outer rim delta O-18(PO4) values that are 1.1 to 2.6 parts per thousand higher compared to the central part of the same bones (delta O-18(PO4) = 18.2 +/- 0.9 parts per thousand, n = 19). Surprisingly, even the collagen-rich bone compacta with low REE contents and apatite crystallinity seems altered, as NS tooth enamel (delta O-18(PO4) =15.0 +/- 0.3 parts per thousand, n=4) has about 3%. lower delta O-18(PO4) values, values that are also similar to those of enamel from RS teeth. Therefore, REE concentration, N content and apatite crystallinity are in this case only poor proxies for the alteration of delta O-18(PO4) values. Seawater exposure of a few years up to 8 kyr can change the delta O-18(PO4) values of the bone apatite by > 3 parts per thousand. Therefore, bones fossilized in marine settings must be treated with caution for palaeoclimatic reconstructions. However, enamel seems to preserve pristine delta O-18(PO4) values on this time scale. Using species-specific calibrations for modern mammals, a mean delta O-18(H2O) value can be reconstructed for Late Pleistocene mammalian drinking water of around -9.2 +/- 0.5 parts per thousand, which is similar to that of Late Pleistocene groundwater from central Europe. (c) 2008 Elsevier B.V. All rights reserved.

The link between land-use management and fluvial flood risk : a chaotic conception?

There is much policy interest in the possible linkages that might exist between land use and downstream fluvial flood risk. On the one hand, this position is sustained by observations from plot- and field-scale studies that suggest land management does affect runoff. On the other, upscaling these effects to show that land-management activities impact upon flood risk at larger catchment scales has proved to be elusive. This review considers the reasons for why this upscaling is problematic. We argue that, rather than it reflecting methodological challenges associated with the difficulties of modelling hydrological processes over very large areas and during extreme runoff events, it reflects the fact that any linkage between land management and flood risk cannot be generalized and taken out of its specific spatial (catchment) and temporal (flood event) context. We use Sayer's (1992) notion of a `chaotic conception' to describe the belief that there is a simple and general association between land management and downstream flood risk rather than the impacts of land management being spatially and temporally contingent in relation to the particular geographical location, time period and scale being considered. Our argument has important practical consequences because it implies that land-management activities to reduce downstream flood risk will be different to traditional flood-reduction interventions such as levees. The purpose of demonstration projects then needs careful consideration such that conclusions made for one project are not transferred uncritically to other scales of analysis or geographical locations.

Sedimentology, geochemistry and mineralogy of the Paleocene-Eocene thermal maximum (PETM) : sediment records from Egypt, India and Spain

A gradual increase in Earth's surface temperatures marking the transition from the late Paleocene to early Eocene (55.8±0.2Ma), represents an extraordinary warming event known as Paleocene-Eocene Thermal Maximum (PETM). Both marine and continental sedimentary records during this period reveal evidences for the massive injection of isotopically light carbon. The carbon dioxide injection from multiple potential sources may have triggered the global warming. The importance of the PETM studies is due to the fact that the PETM bears some striking resemblances to the human-caused climate change unfolding today. Most notably, the culprit behind it was a massive injection of heat-trapping greenhouse gases into the atmosphere and oceans, comparable in volume to what our persistent burning of fossil fuels could deliver in coming centuries. The exact knowledge of what went on during the PETM could help us to foresee the future climate change. The response of the oceanic and continental environments to the PETM is different. Many factors might control the response of the environments to the PETM such as paleogeography, paleotopography, paleoenvironment, and paleodepth. To better understand the mechanisms triggering PETM events, two different environments were studied: 1) shallow marine to inner shelf environment (Wadi Nukhul, Sinai; and the Dababiya GSSP, Luxor, Egypt), and 2) terrestrial environments (northwestern India lignite mines) representing wetland, and fluvial environments (Esplugafreda, Spain) both highlighting the climatic changes observed in continental conditions. In the marine realm, the PETM is characterized by negative ö13Ccar and ô13Corg excursions and shifts in Ô15N to ~0%o values above the P/E boundary and persisting along the interval suggesting a bloom and high production of atmospheric N2-fixers. Decrease in carbonate contents could be due to dissolution and/or dilution by increasing detrital input. High Ti, K and Zr and decreased Si contents at the P/E boundary indicate high weathering index (CIA), which coincides with significant kaolinite input and suggests intense chemical weathering under humid conditions at the beginning of the PETM. Two anoxic intervals are observed along the PETM. The lower one may be linked to methane released from the continental shelf with no change in the redox proxies, where the upper anoxic to euxinic conditions are revealed by increasing U, Mo, V, Fe and the presence of small size pyrite framboids (2-5fim). Productivity sensitive elements (Cu, Ni, and Cd) show their maximum concentrated within the upper anoxic interval suggesting high productivity in surface water. The obtained data highlight that intense weathering and subsequent nutrient inputs are crucial parameters in the chain of the PETM events, triggering productivity during the recovery phase. In the terrestrial environments, the establishment of wetland conditions and consequence continental climatic shift towards more humid conditions led to migration of modern mammals northward following the extension of the tropical belts. Relative ages of this mammal event based on bio-chemo- and paleomagnetic stratigraphy support a migration path originating from Asia into Europe and North America, followed by later migration from Asia into India and suggests a barrier to migration that is likely linked to the timing of the India-Asia collision. In contrast, at Esplugafereda, northeastern Spain, the terrestrial environment reacted differently. Two significant S13C shifts with the lower one linked to the PETM and the upper corresponding to the Early Eocene Thermal Maximum (ETM2); 180/160 paleothermometry performed on two different soil carbonate nodule reveal a temperature increase of around 8°C during the PETM. The prominent increase in kaolinite content within the PETM is linked to increased runoff and/or weathering of adjacent and coeval soils. These results demonstrate that the PETM coincides globally with extreme climatic fluctuations and that terrestrial environments are very likely to record such climatic changes. - La transition Paléocène-Eocène (55,8±0,2 Ma) est marquée par un réchauffement extraordinaire communément appelé « Paleocene-Eocene Thermal Maximum » (PETM). Les données géochimiques caractérisant les sédiments marins et continentaux de cette période indiquent que ce réchauffement a été déclenché par une augmentation massive de CO2 lié à la déstabilisation des hydrates de méthane stockés le long des marges océaniques. L'étude des événements PETM constitue donc un bon analogue avec le réchauffement actuel. Le volume de CO2 émis durant le PETM est comparable avec le CO2 lié à l'activité actuelle humaine. La compréhension des causes du réchauffement du PETM peut être cruciale pour prévoir et évaluer les conséquences du réchauffement anthropogénique, en particulier les répercussions d'un tel réchauffement sur les domaines continentaux et océaniques. De nombreux facteurs entrent en ligne de compte dans le cas du PETM, tels que la paléogéographie, la paléotopographie et les paléoenvironnement. Pour mieux comprendre les réponses environnementales aux événements du PETM, 2 types d'environnements ont été choisis : (1) le domaine marin ouvert mais relativement peu profond (Wadi Nukhul. Sinai, Dababiya, Luxor, Egypte), (2) le milieu continental marécageux humide (mines de lignite, Inde) et fluviatile, semi-aride (Esplugafreda, Pyrénées espagnoles). Dans le domaine marin, le PETM est caractérisé par des excursions négatives du ô13Ccar et ô13Corg et un shift persistant des valeurs de 815N à ~ 0 %o indiquant une forte activité des organismes (bactéries) fixant l'azote. La diminution des carbonates observée durant le PETM peut-être due à des phénomènes de dissolution ou une augmentation des apports terrigènes. Des taux élevés en Ti, K et Zr et une diminution des montants de Si, reflétés par des valeurs des indices d'altération (CIA) qui coïncident avec une augmentation significative des apports de kaolinite impliquent une altération chimique accrue, du fait de conditions plus humides au début du PETM. Deux événements anoxiques globaux ont été mis en évidence durant le PETM. Le premier, situé dans la partie inférieur du PETM, serait lié à la libération des hydrates de méthane stockés le long des talus continentaux et ne correspond pas à des variations significatives des éléments sensibles aux changements de conditions redox. Le second est caractérisé par une augmentation des éléments U, Mo, V et Fe et la présence de petit framboids de pyrite dont la taille varie entre 2 et 5pm. Le second épisode anoxique est caractérisé par une forte augmentation des éléments sensibles aux changements de la productivité (Cu, Ni et Co), indiquant une augmentation de la productivité dans les eaux de surface. Les données obtenues mettent en évidence le rôle crucial joué par l'altération et les apports en nutriments qui en découlent. Ces paramètres sont cruciaux pour la succession des événements qui ont conduit au PETM, et plus particulièrement l'augmentation de la productivité dans la phase de récupération. Durant le PETM, le milieu continental est caractérisé par l'établissement de conditions humides qui ont facilité voir provoqué la migration des mammifères modernes qui ont suivi le déplacement de ces ceintures climatiques. L'âge de cette migration est basé sur des arguments chimiostratigraphiques (isotopes stables), biostratigraphiques et paléomagnétiques. Les données bibliographiques ainsi que celles que nous avons récoltées en Inde, montrent que les mammifères modernes ont d'abord migré depuis l'Asie vers l'Europe, puis dans le continent Nord américain. Ces derniers ne sont arrivés en Inde que plus tardivement, suggérant que le temps de leur migration est lié à la collision Inde-Asie. Dans le Nord-Est de l'Espagne (Esplugafreda), la réponse du milieu continental aux événements PETM est assez différente. Comme en Inde, deux excursions signicatives en ô13C ont été observées. La première correspond au PETM et la seconde est corrélée avec l'optimum thermique de l'Eocène précoce (ETM2). Les isotopes stables de l'oxygène mesurés 2 différents types de nodules calcaires provenant de paléosols suggère une augmentation de 10°C pendant le PETM. Une augmentation simultanée des taux de kaolinite indique une intensification de l'altération chimique et/ou de l'érosion de sols adjacents. Ces résultats démontrent que le PETM coïncide globalement avec des variations climatiques extrêmes qui sont très aisément reconnaissables dans les dépôts continentaux.

Uncertainty Quantification Of A Semi-Supervised Support Vector Regression Reservoir Model

Uncertainty quantification of petroleum reservoir models is one of the present challenges, which is usually approached with a wide range of geostatistical tools linked with statistical optimisation or/and inference algorithms. Recent advances in machine learning offer a novel approach to model spatial distribution of petrophysical properties in complex reservoirs alternative to geostatistics. The approach is based of semisupervised learning, which handles both ?labelled? observed data and ?unlabelled? data, which have no measured value but describe prior knowledge and other relevant data in forms of manifolds in the input space where the modelled property is continuous. Proposed semi-supervised Support Vector Regression (SVR) model has demonstrated its capability to represent realistic geological features and describe stochastic variability and non-uniqueness of spatial properties. On the other hand, it is able to capture and preserve key spatial dependencies such as connectivity of high permeability geo-bodies, which is often difficult in contemporary petroleum reservoir studies. Semi-supervised SVR as a data driven algorithm is designed to integrate various kind of conditioning information and learn dependences from it. The semi-supervised SVR model is able to balance signal/noise levels and control the prior belief in available data. In this work, stochastic semi-supervised SVR geomodel is integrated into Bayesian framework to quantify uncertainty of reservoir production with multiple models fitted to past dynamic observations (production history). Multiple history matched models are obtained using stochastic sampling and/or MCMC-based inference algorithms, which evaluate posterior probability distribution. Uncertainty of the model is described by posterior probability of the model parameters that represent key geological properties: spatial correlation size, continuity strength, smoothness/variability of spatial property distribution. The developed approach is illustrated with a fluvial reservoir case. The resulting probabilistic production forecasts are described by uncertainty envelopes. The paper compares the performance of the models with different combinations of unknown parameters and discusses sensitivity issues.

The role of soil in vegetated gravelly river braid plains: more than just a passive response?

This paper reviews the role of alluvial soils in vegetated gravelly river braid plains. When considering decadal time scales of river evolution, we argue that it becomes vital to consider soil development as an emergent property of the developing ecosystem. Soil processes have been relatively overlooked in accounts of the interactions between braided river processes and vegetation, although soils have been observed on vegetated fluvial landforms. We hypothesise that soil development plays a major role in the transition (speed and pathway) from a fresh sediment deposit to a vegetated soil-covered landform. Disturbance (erosion and/or deposition), vertical sediment structure (process history), vegetation succession, biological activity and water table fluctuation are seen as the main controls on early alluvial soil evolution. Erosion and deposition processes may not only act as soil disturbing agents, but also as suppliers of ecosystem resources, because of their role in delivering and changing access (e.g. through avulsion) to fluxes of water, fine sediments and organic matter. In turn, the associated initial ecosystem may influence further fluvial landform development, such as through the trapping of fine-grained sediments (e.g. sand) by the engineering action of vegetation and the deposit stabilisation by the developing above and belowground biomass. This may create a strong feedback between geomorphological processes, vegetation succession and soil evolution which we summarise in a conceptual model. We illustrate this model by an example from the Allondon River (CH) and identify the research questions that follow.

Splitting rivers at their seams: bifurcations and avulsion

River bifurcations are critical but poorly understood elements of many geomorphological systems. They are integral elements of alluvial fans, braided rivers, fluvial lowland plains, and deltas and control the partitioning of water and sediment through these systems. Bifurcations are commonly unstable but their lifespan varies greatly. In braided rivers bars and channels migrate, split and merge at annual or shorter timescales, thereby creating and abandoning bifurcations. This behaviour has been studied mainly by geomorphologists and fluid dynamicists. Bifurcations also exist during avulsion, the process of a river changing course on a floodplain or in a delta, which may take 102103 years and has been studied mainly by sedimentologists. This review synthesizes our current understanding of bifurcations and brings together insights from different research communities and different environmental settings. We consider the causes and initiation of bifurcations and avulsion, the physical mechanisms controlling bifurcation and avulsion evolution, mathematical and numerical modelling of these processes, and the possibility of stable bifurcations. We end the review with some open questions. Copyright (c) 2012 John Wiley & Sons, Ltd.

La responsabilité du transporteur pour perte, avarie et/ou livraison tardive de la marchandise : étude de droit suisse avec un aperçu de droit romain

Introduction 1-La notion de transport a l'avantage de se laisser aborder sous des angles aussi variés que passionnants. Elle est en effet ancrée dans les fondements de notre société moderne à tous les échelons et a de tout temps constitué un sujet d'étude et de recherche de premier plan. 2-Si les échanges ont connu une véritable "explosion" avec le développement des moyens modernes que l'on sait (moteur, électricité), les historiens voient dans le Néolithique leur véritable point de départ'. S'en suivent de longs siècles qui verront le développement d'un réseau routier de plus en plus dense et de moyens de transport plus ou moins rapides et sophistiqués. 3-Il va sans dire que la notion de transport est intimement liée à celle de commerce. L'on peut même douter de l'existence de ce dernier sans le premier et rejoindre ainsi l'opinion qui fixe leur apparition simultanément, à savoir trois mille ans avant notre ère. 4-Au fil du temps, le commerce donnera naissance à un nombre important d'us et coutumes formant le terreau de notre droit des contrats ; logiquement, le monde des transports n'a pas échappé à cette ébauche de réglementation, dont est issue une législation extraordinairement étayée. 5-Ce domaine constitue un sujet d'étude des plus intéressants en raison de la diversité des questions qu'il soulève, d'une part, et de son rôle dans notre société, d'autre part. La matière ne fait étonnamment l'objet que de peu d'attention en droit suisses, et ce malgré les changements importants intervenus ces derniers temps en transports ferroviaire, aérien, ou encore fluvial. 6-Nous proposons dès lors de présenter l'état du droit suisse sur une question précise, à savoir la responsabilité du transporteur en cas de perte, avarie et/ou livraison tardive de la marchandise. Nous dresserons un tableau de la situation dans les divers modes de transports et comparerons les résultats obtenus, en nous interrogeant sur les éventuelles différences et leur pertinence. Ces développements seront précédés d'un chapitre de droit romain abordant en particulier la question de la responsabilité du transporteur.

Detection and Characterization of Preferential Flow Paths in the Downstream Area of a Hazardous Landfill

We have used surface-based electrical resistivity tomography to detect and characterize preferential hydraulic pathways in the immediate downstream area of an abandoned, hazardous landfill. The landfill occupies the void left by a former gravel pit and its base is close to the groundwater table and lacking an engineered barrier. As such, this site is remarkably typical of many small- to medium-sized waste deposits throughout the densely populated and heavily industrialized foreland on both sides of the Alpine arc. Outflows of pollutants lastingly contaminated local drinking water supplies and necessitated a partial remediation in the form of a synthetic cover barrier, which is meant to prevent meteoric water from percolating through the waste before reaching the groundwater table. Any future additional isolation of the landfill in the form of lateral barriers thus requires adequate knowledge of potential preferential hydraulic pathways for outflowing contaminants. Our results, inferred from a suite of tomographically inverted surfaced-based electrical resistivity profiles oriented roughly perpendicular to the local hydraulic gradient, indicate that potential contaminant outflows would predominantly occur along an unexploited lateral extension of the original gravel deposit. This finds its expression as a distinct and laterally continuous high-resistivity anomaly in the resistivity tomograms. This interpretation is ground-truthed through a litholog from a nearby well. Since the probed glacio-fluvial deposits are largely devoid of mineralogical clay, the geometry of hydraulic and electrical pathways across the pore space of a given lithological unit can be assumed to be identical, which allows for an order-of-magnitude estimation of the overall permeability structure. These estimates indicate that the permeability of the imaged extension of the gravel body is at least two to three orders-of-magnitude higher than that of its finer-grained embedding matrix. This corroborates the preeminent role of the high-resistivity anomaly as a potential preferential flow path.

Modelling hydrodynamics in the Rio Parana, Argentina : an evaluation and inter-comparison of reduced-complexity and physics based models applied to a large sand-bed river

Depth-averaged velocities and unit discharges within a 30 km reach of one of the world's largest rivers, the Rio Parana, Argentina, were simulated using three hydrodynamic models with different process representations: a reduced complexity (RC) model that neglects most of the physics governing fluid flow, a two-dimensional model based on the shallow water equations, and a three-dimensional model based on the Reynolds-averaged Navier-Stokes equations. Row characteristics simulated using all three models were compared with data obtained by acoustic Doppler current profiler surveys at four cross sections within the study reach. This analysis demonstrates that, surprisingly, the performance of the RC model is generally equal to, and in some instances better than, that of the physics based models in terms of the statistical agreement between simulated and measured flow properties. In addition, in contrast to previous applications of RC models, the present study demonstrates that the RC model can successfully predict measured flow velocities. The strong performance of the RC model reflects, in part, the simplicity of the depth-averaged mean flow patterns within the study reach and the dominant role of channel-scale topographic features in controlling the flow dynamics. Moreover, the very low water surface slopes that typify large sand-bed rivers enable flow depths to be estimated reliably in the RC model using a simple fixed-lid planar water surface approximation. This approach overcomes a major problem encountered in the application of RC models in environments characterised by shallow flows and steep bed gradients. The RC model is four orders of magnitude faster than the physics based models when performing steady-state hydrodynamic calculations. However, the iterative nature of the RC model calculations implies a reduction in computational efficiency relative to some other RC models. A further implication of this is that, if used to simulate channel morphodynamics, the present RC model may offer only a marginal advantage in terms of computational efficiency over approaches based on the shallow water equations. These observations illustrate the trade off between model realism and efficiency that is a key consideration in RC modelling. Moreover, this outcome highlights a need to rethink the use of RC morphodynamic models in fluvial geomorphology and to move away from existing grid-based approaches, such as the popular cellular automata (CA) models, that remain essentially reductionist in nature. In the case of the world's largest sand-bed rivers, this might be achieved by implementing the RC model outlined here as one element within a hierarchical modelling framework that would enable computationally efficient simulation of the morphodynamics of large rivers over millennial time scales. (C) 2012 Elsevier B.V. All rights reserved.

Constraining 3-D electrical resistance tomography with GPR reflection data forimproved aquifer characterization

Surface-based ground penetrating radar (GPR) and electrical resistance tomography (ERT) are common tools for aquifer characterization, because both methods provide data that are sensitive to hydrogeologically relevant quantities. To retrieve bulk subsurface properties at high resolution, we suggest incorporating structural information derived from GPR reflection data when inverting surface ERT data. This reduces resolution limitations, which might hinder quantitative interpretations. Surface-based GPR reflection and ERT data have been recorded on an exposed gravel bar within a restored section of a previously channelized river in northeastern Switzerland to characterize an underlying gravel aquifer. The GPR reflection data acquired over an area of 240×40 m map the aquifer's thickness and two internal sub-horizontal regions with different depositional patterns. The interface between these two regions and the boundary of the aquifer with then underlying clay are incorporated in an unstructured ERT mesh. Subsequent inversions are performed without applying smoothness constraints across these boundaries. Inversion models obtained by using these structural constraints contain subtle resistivity variations within the aquifer that are hardly visible in standard inversion models as a result of strong vertical smearing in the latter. In the upper aquifer region, with high GPR coherency and horizontal layering, the resistivity is moderately high (N300 Ωm). We suggest that this region consists of sediments that were rearranged during more than a century of channelized flow. In the lower low coherency region, the GPR image reveals fluvial features (e.g., foresets) and generally more heterogeneous deposits. In this region, the resistivity is lower (~200 Ωm), which we attribute to increased amounts of fines in some of the well-sorted fluvial deposits. We also find elongated conductive anomalies that correspond to the location of river embankments that were removed in 2002.

Stochastic inversion of tracer test and electrical geophysical data to estimate hydraulic conductivities.

Quantifying the spatial configuration of hydraulic conductivity (K) in heterogeneous geological environments is essential for accurate predictions of contaminant transport, but is difficult because of the inherent limitations in resolution and coverage associated with traditional hydrological measurements. To address this issue, we consider crosshole and surface-based electrical resistivity geophysical measurements, collected in time during a saline tracer experiment. We use a Bayesian Markov-chain-Monte-Carlo (McMC) methodology to jointly invert the dynamic resistivity data, together with borehole tracer concentration data, to generate multiple posterior realizations of K that are consistent with all available information. We do this within a coupled inversion framework, whereby the geophysical and hydrological forward models are linked through an uncertain relationship between electrical resistivity and concentration. To minimize computational expense, a facies-based subsurface parameterization is developed. The Bayesian-McMC methodology allows us to explore the potential benefits of including the geophysical data into the inverse problem by examining their effect on our ability to identify fast flowpaths in the subsurface, and their impact on hydrological prediction uncertainty. Using a complex, geostatistically generated, two-dimensional numerical example representative of a fluvial environment, we demonstrate that flow model calibration is improved and prediction error is decreased when the electrical resistivity data are included. The worth of the geophysical data is found to be greatest for long spatial correlation lengths of subsurface heterogeneity with respect to wellbore separation, where flow and transport are largely controlled by highly connected flowpaths.

Dynamic drag modeling of submerged aquatic vegetation canopy flows

Vegetation has a profound effect on flow and sediment transport processes in natural rivers, by increasing both skin friction and form drag. The increase in drag introduces a drag discontinuity between the in-canopy flow and the flow above, which leads to the development of an inflection point in the velocity profile, resembling a free shear layer. Therefore, drag acts as the primary driver for the entire canopy system. Most current numerical hydraulic models which incorporate vegetation rely either on simple, static plant forms, or canopy-scaled drag terms. However, it is suggested that these are insufficient as vegetation canopies represent complex, dynamic, porous blockages within the flow, which are subject to spatially and temporally dynamic drag forces. Here we present a dynamic drag methodology within a CFD framework. Preliminary results for a benchmark cylinder case highlight the accuracy of the method, and suggest its applicability to more complex cases.