Discontinuous Galerkin methods for modeling Hurricane storm surge

Storm surge due to hurricanes and tropical storms can result in significant loss of life, property damage, and long-term damage to coastal ecosystems and landscapes. Computer modeling of storm surge can be used for two primary purposes: forecasting of surge as storms approach land for emergency planning and evacuation of coastal populations, and hindcasting of storms for determining risk, development of mitigation strategies, coastal restoration and sustainability.Storm surge is modeled using the shallow water equations, coupled with wind forcing and in some events, models of wave energy. In this paper, we will describe a depth-averaged (2D) model of circulation in spherical coordinates. Tides, riverine forcing, atmospheric pressure, bottom friction, the Coriolis effect and wind stress are all important for characterizing the inundation due to surge. The problem is inherently multi-scale, both in space and time. To model these problems accurately requires significant investments in acquiring high-fidelity input (bathymetry, bottom friction characteristics, land cover data, river flow rates, levees, raised roads and railways, etc.), accurate discretization of the computational domain using unstructured finite element meshes, and numerical methods capable of capturing highly advective flows, wetting and drying, and multi-scale features of the solution.The discontinuous Galerkin (DG) method appears to allow for many of the features necessary to accurately capture storm surge physics. The DG method was developed for modeling shocks and advection-dominated flows on unstructured finite element meshes. It easily allows for adaptivity in both mesh (h) and polynomial order (p) for capturing multi-scale spatial events. Mass conservative wetting and drying algorithms can be formulated within the DG method.In this paper, we will describe the application of the DG method to hurricane storm surge. We discuss the general formulation, and new features which have been added to the model to better capture surge in complex coastal environments. These features include modifications to the method to handle spherical coordinates and maintain still flows, improvements in the stability post-processing (i.e. slope-limiting), and the modeling of internal barriers for capturing overtopping of levees and other structures. We will focus on applications of the model to recent events in the Gulf of Mexico, including Hurricane Ike.     

Bimodal wave spectra in lower Chesapeake Bay, sea bed energetics and sediment transport during winter storms

The transition zone separating estuarine environments from the coastal ocean is characterized not only by distinctive morphological and sedimentary trends but by unique hydrodynamic forces as well. Lower Chesapeake Bay, a large coastal estuary within the Mid-Atlantic Bight of the U.S. East Coast, experiences complex wave and current-induced forces produced during winter storms. Wave and current measurements made near Thimble Shoal Light over five winter seasons show that most storms simultaneously produce both ocean and bay-generated wave trains that appear as distinct bimodal peaks in directional spectra. Analysis of selected storm wave records reveal that lower-frequency ocean waves, although nominally lower in amplitude than higher-frequency bay waves, are roughly equivalent to bay waves in terms of energy expended on beds of fine- to medium-grained sand at either end of the Thimble Shoal Channel. Grain-friction energy dissipation estimates calculated for waves and currents suggest that waves provide more net energy capable of transporting bottom sediment than currents, although strong barotropic flows briefly encountered during a major storm on 13–14 March 1993, exceeded wave energy expended at the bed by almost an order of magnitude. From analyses of wave orbital velocity spectra, it is shown that dual wave trains characterized by differences in peak frequency and direction may assist each other through interactions that increase their combined contribution to frictional energy dissipation and inferred sediment transport at the bed.     

A comparison of biohydrodynamic interaction within mangrove and saltmarsh boundaries

Mangrove forests and saltmarshes are recognized for their roles in wave and current attenuation, although a comparison of in situ observations between woody and herbaceous plants is needed in order to understand the different mechanisms of bio‐physical interaction within coastal wetlands. The aim of our study was to compare the mechanisms of flow reduction and energy dissipation by mangrove trees and saltmarsh grass in a subtropical area where tidal currents dominate. Fieldwork was conducted to measure the hydrodynamic processes occurring at the boundaries between a bare mudflat and vegetated tidal flat, as the flow transitions from a bare mudflat to either mangrove or saltmarsh. Synchronous acoustic Doppler velocimetry (ADV) measurements at three sites revealed that the mangrove was more effective than the saltmarsh grass at flow reduction. In addition, a considerable rotation in flow direction was observed as the flow entered the mangrove trees, while rotation was considerably less pronounced within the saltmarsh edge. The mechanism for this difference was explained through a combination of changes in drag force and eddy viscosity over the two vegetation types. Although overall the mangrove was observed to dissipate energy more effectively than the saltmarsh, the relative efficiency of the vegetation at dissipating turbulent energy was found to vary with the maximum water level of tidal cycle. When the maximum water level remained below the mangrove canopy bottom (‘bio‐line’), the energy dissipation ability of the mangrove was relatively low, as a result of the presence of rigid, sparse trunks rather than denser saltmarsh grass found near the bed; when the maximum water level was sufficiently high to reach the mangrove canopy, the ability of the mangrove to dissipate energy was significantly increased, becoming more effective than the saltmarsh grass. Copyright © 2016 John Wiley & Sons, Ltd.     

Alteration of bottom roughness by benthic organisms in a sandy coastal environment

This work deals with the potential influence of benthic communities on the sediment dynamics of the coastal zone, and specifically with the modifications to bottom roughness caused by communities and their effects on wave propagation across the coastal profile. Time-series of video observations of the sea bottom on the Ebro delta coast were analysed in order to estimate the bottom roughness associated with physical and biological morphological components and the sediment reworking rates caused by epifaunal organisms. Biological roughness was mainly caused by ophiuroids and tanatocenosis of bivalve and gastropod shells, which changed their abundance during the study period. The total biological roughness (K_bio) ranged between 0.27 and 0.81 cm and represented a significant part (<20%) of the total form drag roughness. Flattening of ripples caused by bioturbation was observed under low-energy conditions. Surface sediment perturbation and bioturbation rates were also estimated. Based on these observations a wave propagation model was applied in order to carry out a sensitivity analysis of the significance of biological roughness on wave dissipation in the study area under different wave conditions.     

Coastal Resiliency Groundwater Considerations

The potential for rising groundwater is an important consideration in any coastal resiliency assessment. Unlike other groundwater modeling that focuses mostly on contaminant tracking, coastal groundwater resiliency assessments are primarily concerned with the potential for groundwater emergence induced by sea level rise. This provides more options for modelers that range from simplified water table elevation models to fully integrated groundwater and storm water models. The selection is dependent on available data and project needs. However, despite the relative simplicity of some of the techniques, all the methods benefit from a professional with hydrogeological training.     

A study of dissipation of wind-waves by mud at Cassino Beach,Brazil: Prediction and inversion

The impact of a non-rigid seafloor on the wave climate at Cassino Beach, Brazil, May–June 2005 is studied using field measurements and a numerical wave model. The measurements consist of wave data at four locations; rheology and mud thickness from grab samples; and an estimate of the horizontal distribution of mud based on echo-soundings. The dissipation of waves by a non-rigid bottom is represented in the wave model by treating the mud layer as a viscous fluid. Applied for 431 time periods, the model without this type of dissipation has a strong tendency to overpredict nearshore wave energy, except during a period of large storm waves. Two model variations which include this dissipation have a modest tendency to underpredict the nearshore wave energy. An inversion methodology is developed and applied to infer an alternate mud distribution which, when used with the wave model, yields the observed waveheights.     

Effect of offshore waves and vegetation on the sediment budget in the Virginia Coast Reserve (VA)

The potential for rapid coastline modification in the face of sea-level rise or other stressors is alarming, since coasts are often densely populated and support valuable infrastructure. In addition to coastal submergence, nutrient-related water pollution is a growing concern for coastal wetlands. Previous studies found that the Suspended Sediment Concentration (SSC) of coastal wetlands acts as a first-order control of their sustainability, but SSC dynamics are poorly understood. Our study focuses on the Virginia Coast Reserve (VCR) Long Term Ecological Research (LTER) site, a shallow multiple tidal inlet system in the USA. We apply numerical modelling (Delft3D-SWAN) and subsequent analyses to determine SSC dynamics within the VCR. In particular, we consider two important controls on SSC in the system: vegetation (seagrass and salt marsh) and offshore waves. Our results show that vegetation colonies and increased wave energy lengthen water residence time. The reduction in the tidal prism decreases SSC export from the bay via tidal inlets, leading to increased sediment retention in the bay. We found that alongshore currents can enhance lagoon SSC by importing fine sediments from an adjacent inlet along the coastline. Our numerical experiments on vegetation seasonality can improve the understanding of wave climate impact on coastal bay sediment budget. Offshore waves increase sediment export from coastal bays, particularly during winter seasons with low vegetation density. Therefore, our study can help managers and stakeholders to understand how to implement restoration strategies for the VCR. © 2020 John Wiley & Sons, Ltd.     

Landscapes on the edge: examining the role of climatic interactions in shaping coastal watersheds using a coastal–terrestrial landscape evolution model

Incised coastal gullies (ICGs) are dynamic features found at the terrestrial‐coastal interface. Their geomorphic evolution is driven by the interactions between processes of fluvial knickpoint migration and coastal cliff erosion. Under scenarios of future climate change the frequency and magnitude of the climatological drivers of both terrestrial (fluvial and hillslope) and coastal (cliff erosion) processes are likely to change, with an adjunct impact on these types of coastal features. Here we explore the response of an incised coastal gully to changes in both terrestrial and coastal climate in order to elucidate the key process interactions which drive ICG evolution. We modify an extant landscape evolution model, CHILD, to incorporate processes of soft‐cliff erosion. This modified version, termed the Coastal‐Terrestrial‐CHILD (CT‐CHILD) model, is then employed to explore the interactions between changing terrestrial and coastal driving forces on the future evolution of an ICG found on the south‐west Isle of Wight, UK. It was found that the magnitude and frequency of storm events will play a key role in determining the future trajectory of ICGs, highlighting a need to understand the role of event sequencing in future projections of landscape evolution. Furthermore, synergistic (positive) and antagonistic (negative) interactions were identified between coastal and terrestrial parameters, such as wave height intensity and precipitation duration, which act to modulate the impact of changes in any one parameter. Of note was the role played by wave height intensity in driving coastal erosion, which was found to play a more important role than sea‐level rise in determining rates of coastal erosion. This highlights the need for a greater focus on wave height in studies of soft‐cliff erosion. Copyright © 2014 John Wiley & Sons, Ltd.     

The extreme 2013/2014 winter storms: hydrodynamic forcing and coastal response along the southwest coast of England

The southwest coast of England was subjected to an unusually energetic sequence of Atlantic storms during the 2013/2014 winter, with the 8‐week period from mid‐December to mid‐February representing the most energetic period since at least 1953. A regional analysis of the hydrodynamic forcing and morphological response of these storms along the SW coast of England highlighted the importance of both storm‐ and site‐specific conditions. The key factor that controls the Atlantic storm wave conditions along the south coast of southwest England is the storm track. Energetic inshore wave conditions along this coast require a relatively southward storm track which enables offshore waves to propagate up the English Channel relatively unimpeded. The timing of the storm in relation to the tidal stage is also important, and coastal impacts along the macro‐tidal southwest coast of England are maximised when the peak storm waves coincide with spring high tide. The role of storm surge is limited and rarely exceeds 1 m. The geomorphic storm response along the southwest coast of England displayed considerable spatial variability; this is mainly attributed to the embayed nature of the coastline and the associated variability in coastal orientation. On west‐facing beaches typical of the north coast, the westerly Atlantic storm waves approached the coastline shore‐parallel, and the prevailing storm response was offshore sediment transport. Many of these north coast beaches experienced extensive beach and dune erosion, and some of the beaches were completely stripped of sediment, exposing a rocky shore platform. On the south coast, the westerly Atlantic storm waves refract and diffract to become southerly inshore storm waves and for the southeast‐facing beaches this results in large incident wave angles and strong eastward littoral drift. Many south coast beaches exhibited rotation, with the western part of the beaches eroding and the eastern part accreting. © 2015 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

A numerical study on the effects of wave–current–surge interactions on the height and propagation of sea surface waves in Charleston Harbor during Hurricane Hugo 1989

The effects of wave–current interactions on ocean surface waves induced by Hurricane Hugo in and around the Charleston Harbor and its adjacent coastal waters are examined by using a three-dimensional (3D) wave–current coupled modeling system. The 3D storm surge modeling component of the coupled system is based on the Princeton Ocean Model (POM), the wave modeling component is based on the third generation wave model, Simulating WAves Nearshore (SWAN), and the inundation model is adopted from [Xie, L., Pietrafesa, L. J., Peng, M., 2004. Incorporation of a mass-conserving inundation scheme into a three-dimensional storm surge model. J. Coastal Res., 20, 1209–1223]. The results indicate that the change of water level associated with the storm surge is the primary cause for wave height changes due to wave–surge interaction. Meanwhile, waves propagating on top of surge cause a feedback effect on the surge height by modulating the surface wind stress and bottom stress. This effect is significant in shallow coastal waters, but relatively small in offshore deep waters. The influence of wave–current interaction on wave propagation is relatively insignificant, since waves generally propagate in the direction of the surface currents driven by winds. Wave–current interactions also affect the surface waves as a result of inundation and drying induced by the storm. Waves break as waters retreat in regions of drying, whereas waves are generated in flooded regions where no waves would have occurred without the flood water.     

Modeling boulder transport by coastal waves on cliff topography: Case study at Hachijo Island,Japan

Clifftop coastal boulders transported by storm waves or tsunamis have been reported around the world. Although numerical calculation of boulder transport is a strong tool for the identification of tsunami or storm boulders, and for estimation of the wave size emplacing boulders, models which can reasonably solve boulder transport from below a cliff or from a cliff-edge onto a cliff-top do not yet exist. In this study, we developed a new numerical formulation for cliff-top deposition of boulders from the cliff edge or below the cliff, with validation from laboratory tests. We then applied the model using storm and tsunami wave forcing to simulate the observed boulder deposits at the northwest coast of Hachijo Island, Japan. Using the model, the actual distribution of boulders was explained well using a reasonable storm wave height without assumption of anomalously high-water level by storm surge. Results show that boulder transport from the cliff edge or under the cliff onto the cliff-top was possible from a tsunami with periods of 5~10 min or storm waves with no storm surge. However, the actual distribution of boulders on the cliff was explained only by storm waves, but not by tsunami. Therefore, the boulders distributed at this site are likely of storm wave origin. Our developed model for the boulder transport calculation can be useful for identifying a boulder's origin and can reasonably calculate cliff-top deposition of boulders by tsunami and storm waves. © 2019 John Wiley & Sons, Ltd. © 2019 John Wiley & Sons, Ltd.     

Geomorphological control on boulder transport and coastal erosion before,during and after an extreme extra‐tropical cyclone

Extreme wave events in coastal zones are principal drivers of geomorphic change. Evidence of boulder entrainment and erosional impact during storms is increasing. However, there is currently poor time coupling between pre‐ and post‐storm measurements of coastal boulder deposits. Importantly there are no data reporting shore platform erosion, boulder entrainment and/or boulder transport during storm events – rock coast dynamics during storm events are currently unexplored. Here, we use high‐resolution (daily) field data to measure and characterize coastal boulder transport before, during and after the extreme Northeast Atlantic extra‐tropical cyclone Johanna in March 2008. Forty‐eight limestone fine‐medium boulders (n = 46) and coarse cobbles (n = 2) were tracked daily over a 0.1 km² intertidal area during this multi‐day storm. Boulders were repeatedly entrained, transported and deposited, and in some cases broken down (n = 1) or quarried (n = 3), during the most intense days of the storm. Eighty‐one percent (n = 39) of boulders were located at both the start and end of the storm. Of these, 92% were entrained where entrainment patterns were closely aligned to wave parameters. These data firmly demonstrate rock coasts are dynamic and vulnerable under storm conditions. No statistically significant relationship was found between boulder size (mass) and net transport distance. Graphical analyses suggest that boulder size limits the maximum longshore transport distance but that for the majority of boulders lying under this threshold, other factors influence transport distance. Paired analysis of 20 similar sized and shaped boulders in different morphogenic zones demonstrates that geomorphological control affects entrainment and transport distance – where net transport distances were up to 39 times less where geomorphological control was greatest. These results have important implications for understanding and for accurately measuring and modelling boulder entrainment and transport. Coastal managers require these data for assessing erosion risk. © 2016 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Wetland hydroperiod classification in the western prairies using multitemporal synthetic aperture radar

Wetlands represent one of the world's most biodiverse and threatened ecosystem types and were diminished globally by about two‐thirds in the 20th century. There is continuing decline in wetland quantity and function due to infilling and other human activities. In addition, with climate change, warmer temperatures and changes in precipitation and evapotranspiration are reducing wetland surface and groundwater supplies, further altering wetland hydrology and vegetation. There is a need to automate inventory and monitoring of wetlands, and as a study system, we investigated the Shepard Slough wetlands complex, which includes numerous wetlands in urban, suburban, and agricultural zones in the prairie pothole region of southern Alberta, Canada. Here, wetlands are generally confined to depressions in the undulating terrain, challenging wetlands inventory and monitoring. This study applied threshold and frequency analysis routines for high‐resolution, single‐polarization (HH) RADARSAT‐2, synthetic aperture radar mapping. This enabled a growing season surface water extent hyroperiod‐based wetland classification, which can support water and wetland resource monitoring. This 3‐year study demonstrated synthetic aperture radar‐derived multitemporal open‐water masks provided an effective index of wetland permanence class, with overall accuracies of 89% to 95% compared with optical validation data, and RMSE between 0.2 and 0.7 m between model and field validation data. This allowed for characterizing the distribution and dynamics of 4 marsh wetlands hydroperiod classes, temporary, seasonal, semipermanent, and permanent, and mapping of the sequential vegetation bands that included emergent, obligate wetland, facultative wetland, and upland plant communities. Hydroperiod variation and surface water extent were found to be influenced by short‐term rainfall events in both wet and dry years. Seasonal hydroperiods in wetlands were particularly variable if there was a decrease in the temporary or semipermanent hydroperiod classes. In years with extreme rain events, the temporary wetlands especially increased relative to longer lasting wetlands (84% in 2015 with significant rainfall events, compared with 42% otherwise).     

Computing near coastal tidal dynamics from observations and a numerical model

A two-dimensional vertically integrated model of the North Sea is used to compute the distribution of M₂ and M₄ tidal elevations and currents over the region. Comparison of computed and observed elevations and currents in the area shows that the model can accurately reproduce the M₂ tide in the North Sea, although there are difficulties with the M₄ tide particularly in the northern North Sea.Comparison between model and a large number of observations collected in a shallow water region off the east coast of England, revealed that the model can accurately reproduce the tides even in near coastal regions, where model resolution problems can occur. Comparisons of computed and observed M₂ tidal energy fluxes in this region, show that model and observations agree to within the order of 10% (the error associated with the necessary interpolation of the observations in order to compute the energy flux).The problem of computing energy dissipation in the area by subtracting the energy fluxes into and out of the region is shown to be ill-conditioned in that the energy dissipation in the area is comparable to the error in the energy flux. Consequently for the sea region considered here it is not meaningful to compare this energy budget with energy dissipation due to bottom friction.Energy dissipation for the whole of the North Sea is computed using the numerical model and the geographical distribution of dissipation due to bottom friction is given for the M₂ tide.     

Marine control over negative power law scaling of mass wasting events in chalk sea cliffs with implications for future recession under the UKCP09 medium emission scenario

Coastal cliff erosion represents a significant geohazard for people and infrastructure. Forecasting future erosion rates is therefore of critical importance to ensuring the resiliency of coastal communities. We use high precision monitoring of chalk cliffs at Telscombe, UK to generate monthly mass movement inventories between August 2016 and July 2017. Frequency–magnitude analysis of our inventories demonstrate negative power law scaling over 7 orders of magnitude and, for the first time, we report statistically significant correlations between significant wave height (H_s) and power law scaling coefficients (r² values of 0.497 and 0.590 for β and s respectively). Applying these relationships allows for a quantitative method to predict erosion at the site based on H_s probabilities and sea level forecasts derived from the UKCP09 medium emission climate model (A1B). Monte‐Carlo simulations indicate a range of possible erosion scenarios over 70 years (2020–2090) and we assess the impact these may have on the A259 coastal road which runs proximal to the cliffs. Results indicate a small acceleration in erosion compared with those based on current conditions with the most likely scenario at the site being 21.7 m of cliff recession by 2090. However, low‐probability events can result in recession an order of magnitude higher in some scenarios. In the absence of negative feedbacks, we estimate an ~11% chance that the A259 will be breached by coastal erosion by 2090. Copyright © 2018 John Wiley & Sons, Ltd.     

Simulating short‐circuiting flow in a constructed wetland: the implications of bathymetry and vegetation effects

Short‐circuiting flow, commonly experienced in many constructed wetlands, reduces hydraulic retention times in unit wetland cells and decreases the treatment efficiency. A two‐dimensional (2‐D), physically based, distributed modelling approach was used to systematically address the effects of bathymetry and vegetation on short‐circuiting flow, which previously have been neglected or lumped in one‐dimensional wetland flow models. In this study, a 2‐D transient hydrodynamics with advection‐dispersion model was developed using MIKE 21 and calibrated with bromide tracer data collected at the Orlando Easterly Wetland Cell 7. The estimated topographic difference between short‐circuiting flow zone and adjacent area ranged from 0·3 to 0·8 m. A range of the Manning roughness coefficient at the short‐circuiting flow zone was estimated (0·022–0·045 s m^?1/3). Sensitivity analysis of topographical and vegetative heterogeneity deduced during model calibration shows that relic ditches or other ditch‐shaped landforms and the associated sparse vegetation along the main flow direction intensify the short‐circuiting pattern, considerably affecting 2‐D solute transport simulation. In terms of hydraulic efficiency, this study indicates that the bathymetry effect on short‐circuiting flow is more important than the vegetation effect. Copyright © 2009 John Wiley & Sons, Ltd.     

Two‐dimensional hydraulic flood modelling using a finite‐element mesh decomposed according to vegetation and topographic features derived from airborne scanning laser altimetry

Airborne scanning laser altimetry (LiDAR) is an important new data source that can provide two‐dimensional river flood models with spatially distributed floodplain topography for model bathymetry, together with vegetation heights for parameterization of model friction. Methods are described for improving such models by decomposing the model's finite‐element mesh to reflect floodplain vegetation features such as hedges and trees having different frictional properties to their surroundings, and significant floodplain topographic features having high height curvatures. The decomposition is achieved using an image segmentation system that converts the LiDAR height image into separate images of surface topography and vegetation height at each point. The vegetation height map is used to estimate a friction factor at each mesh node. The spatially distributed friction model has the advantage that it is physically based, and removes the need for a model calibration exercise in which free parameters specifying friction in the channel and floodplain are adjusted to achieve best fit between modelled and observed flood extents. The scheme was tested in a modelling study of a flood that occurred on the River Severn, UK, in 1998. A satellite synthetic aperture radar image of flood extent was used to validate the model predictions. The simulated hydraulics using the decomposed mesh gave a better representation of the observed flood extent than the more simplistic but computationally efficient approach of sampling topography and vegetation friction factors on to larger floodplain elements in an undecomposed mesh, as well as the traditional approach using no LiDAR‐derived data but simply using a constant floodplain friction factor. Use of the decomposed mesh also allowed velocity variations to be predicted in the neighbourhood of vegetation features such as hedges. These variations could be of use in predicting localized erosion and deposition patterns that might result in the event of a flood. Copyright © 2003 John Wiley & Sons, Ltd.     

Evidence for elevated coastal vulnerability following large‐scale historical oyster bed harvesting

Living coastal barriers, such as coral reefs, tidal marshes, mangroves and shellfish beds are widely recognized for their potential role in mitigating flood risk. Limited data exists, however, for assessing the effectiveness of these natural defenses as forms of flood mitigation. In particular, very few mature shellfish beds exist today for modern study due to their destruction in the past few centuries. As an alternative method of study, we present here sedimentary reconstructions of storm overwash from coastal ponds internal to New York Harbor. We use these reconstructions to show that the initial degradation of oyster beds following European settlement of the area coincides with a significant increase in wave‐derived overwash deposition at all three of our field sites. Numerical simulations of two flood events of record in the harbor (Hurricane Sandy and a severe winter storm in 1992) were run without and with oyster beds of varying heights (1 m above the seafloor‐to‐intertidal). Simulations show that the removal of these oyster beds increases wave energy directly off‐shore of our field sites by between 30% and 200%. Sedimentary reconstructions and wave modeling experiments therefore both support oyster beds serving as a significant form of coastal protection prior to European disturbance. Copyright © 2016 John Wiley & Sons, Ltd.     

How is water availability related to the land use and morphology of an inland valley wetland in Kenya?

Small inland valley wetlands contribute substantially to the livelihoods of rural communities in East Africa. Their conversion into farmland is driven by water availability. We quantified spatial-temporal dynamics of water availability in a headwater wetland in the humid zone of Kenya. Climatic conditions, soil moisture contents, groundwater levels and discharge data were monitored. A land-use map and a digital elevation model of the valley bottom were created to relate variations in soil moisture to dominant land uses and valley morphology.Upland crops occupied about a third of the wetland area, while approximately a quarter of the wet, central part of the valley bottom was designated for flood-tolerant taro, grown either by itself or in association or in rotation with upland crops. Finally, natural vegetation was found in 3% of the mapped area, mainly in sections with nearpermanent soil saturation.The HBV rainfall-runoff model's overestimation of stream discharge during the long dry season of the hydrological year 2010/2011 can be explained by the strong seasonal impact of water abstraction on the wetland's water balance.Our study vividly demonstrates the necessity of multi-method approaches for assessing the impact of management practices on water availability in valley bottom wetlands in East Africa.