A coupled model for simulating surface water and groundwater interactions in coastal wetlands

Coastal wetlands are characterized by strong, dynamic interactions between surface water and groundwater. This paper presents a coupled model that simulates interacting surface water and groundwater flow and solute transport processes in these wetlands. The coupled model is based on two existing (sub) models for surface water and groundwater, respectively: ELCIRC (a three‐dimensional (3‐D) finite‐volume/finite‐difference model for simulating shallow water flow and solute transport in rivers, estuaries and coastal seas) and SUTRA (a 3‐D finite‐element/finite‐difference model for simulating variably saturated, variable‐density fluid flow and solute transport in porous media). Both submodels, using compatible unstructured meshes, are coupled spatially at the common interface between the surface water and groundwater bodies. The surface water level and solute concentrations computed by the ELCIRC model are used to determine the boundary conditions of the SUTRA‐based groundwater model at the interface. In turn, the groundwater model provides water and solute fluxes as inputs for the continuity equations of surface water flow and solute transport to account for the mass exchange across the interface. Additionally, flux from the seepage face was routed instantaneously to the nearest surface water cell according to the local sediment surface slope. With an external coupling approach, these two submodels run in parallel using time steps of different sizes. The time step (Δt_g) for the groundwater model is set to be larger than that (Δt_s) used by the surface water model for computational efficiency: Δt_g = M × Δt_s where M is an integer greater than 1. Data exchange takes place between the two submodels through a common database at synchronized times (e.g. end of each Δt_g). The coupled model was validated against two previously reported experiments on surface water and groundwater interactions in coastal lagoons. The results suggest that the model represents well the interacting surface water and groundwater flow and solute transport processes in the lagoons. Copyright © 2011 John Wiley & Sons, Ltd.     

Analysis of geomorphology citations in the last quarter of the 20th century

Three hundred and twenty‐eight geomorphology articles published in the last quarter of the 20th century were cited 20 or more times in Institute for Scientific Information (ISI) indices, as of 15 May 2001. At the close of the 20th century, well‐cited geomorphology is highly multidisciplinary and interdisciplinary with the most dominant fields being in biological, civil engineering, earth science, geography, geological, and soils disciplines. The very strong English‐language bias of well‐cited journal articles creates a geographical bias in study site selection, which may in turn bias geomorphic theory. Water‐based research (fluvial processes and landforms, riparian, drainage basin) dominates well‐cited papers, with the ‘hottest’ subfield in the 1990s being riparian research with a biological emphasis. Over 90 journals publish well‐cited papers, but Earth Surface Processes and Landforms hosts the largest number of well‐cited papers. Copyright © 2002 John Wiley & Sons, Ltd.     

Towards quantifying the negative feedback regulation of peatland evaporation to drought

The frequency and intensity of drought is projected to increase within the boreal region under future climatic conditions. Peatlands are widely considered to regulate water loss under drought conditions, increasing surface resistance (r_s) and reducing evaporative losses. This maintains peat moisture content, increasing the resilience of these globally important carbon stores. However, the magnitude and form of this important negative feedback response remains uncertain. To address this, we monitored the response of r_s to drought within four peat cores under controlled meteorological conditions. When the water‐table was dropped to a depth of 0.30 m and the humidity reduced to ≤40%, a step shift in r_s from ~50 s m^‐1 up to 1000 s m^‐1 was observed within burned and unburned peat, which virtually shuts down evaporation, limiting water loss. We show that measured near‐surface tension cannot be used to directly calculate this transition in peat surface resistance. However, empirical relationships that account for strong vertical variations in tension through the near‐surface and/or disequilibrium between pore air and near‐surface pore water pressure provide the potential to incorporate this negative feedback response into peatland ecohydrological models. Further observations are necessary to examine this response under dynamic atmospheric conditions. We suggest that the link between surface temperature and evaporation provides potential to further examine this feedback in either burned peatlands or peatlands with a low vascular vegetation cover. Copyright © 2013 John Wiley & Sons, Ltd.     

First‐order exchange coefficient coupling for simulating surface water–groundwater interactions: parameter sensitivity and consistency with a physics‐based approach

Distributed hydrologic models capable of simulating fully‐coupled surface water and groundwater flow are increasingly used to examine problems in the hydrologic sciences. Several techniques are currently available to couple the surface and subsurface; the two most frequently employed approaches are first‐order exchange coefficients (a.k.a., the surface conductance method) and enforced continuity of pressure and flux at the surface‐subsurface boundary condition. The effort reported here examines the parameter sensitivity of simulated hydrologic response for the first‐order exchange coefficients at a well‐characterized field site using the fully coupled Integrated Hydrology Model (InHM). This investigation demonstrates that the first‐order exchange coefficients can be selected such that the simulated hydrologic response is insensitive to the parameter choice, while simulation time is considerably reduced. Alternatively, the ability to choose a first‐order exchange coefficient that intentionally decouples the surface and subsurface facilitates concept‐development simulations to examine real‐world situations where the surface‐subsurface exchange is impaired. While the parameters comprising the first‐order exchange coefficient cannot be directly estimated or measured, the insensitivity of the simulated flow system to these parameters (when chosen appropriately) combined with the ability to mimic actual physical processes suggests that the first‐order exchange coefficient approach can be consistent with a physics‐based framework. Copyright © 2009 John Wiley & Sons, Ltd.     

Performance of grass and eucalyptus riparian buffers in a pasture catchment,Western Australia,part 1: riparian hydrology

Rainfall takes many flowpaths to reach a stream, and the success of riparian buffers in water quality management is significantly influenced by riparian hydrology. This paper presents results from hydrometric monitoring of riparian buffer hydrology in a pasture catchment. Runoff processes and riparian flowpaths were investigated on two planar hillslopes with regenerating grass and E. globulus buffers. Surface runoff and subsurface flows (A‐ and B‐horizons) were measured for 3 years using surface runoff collectors, subsurface troughs and piezometers. Water volumes moving through the riparian buffers via the measured flowpaths were ranked B‐horizon ? surface runoff ≈ A‐horizon. Runoff volumes through the B‐horizon troughs were an order of magnitude greater than those recorded for the most productive surface runoff plots or the A‐horizon troughs. Subsurface runoff and saturation‐excess overland flow (SOF) were limited to the winter months, whereas infiltration‐excess overland flow (IEOF) can occur all year round during intense storms. Surface runoff was recorded on 33 occasions, mostly during winter (late May–early October), and total annual surface runoff volumes collected by the 20 unconfined (2 m wide) runoff plots varied between > 80 and < 20 m³. Subsurface flow only occurred in winter, and the 6 m wide B‐horizon subsurface troughs flowed above 1 l s^?1 continuously, whereas the A‐horizon troughs flowed infrequently (<6 days per year). In summer, surface runoff occurred as IEOF during intense storms in the E. globulus buffer, but not in the grass buffer. Observations suggest that surface crusting reduced the soil's infiltration capacity in the E. globulus buffer. During winter, SOF and seepage were observed in both buffers, but subsurface flow through the B‐horizon was the dominant flowpath. Key hydrologic differences between the grass and tree buffers are the generation of IEOF in the E. globulus buffer during intense summer storms, and the smaller subsurface runoff volumes and fewer flow days in the E. globulus buffer. Low surface runoff volumes are likely to limit the potential of these buffers to filter pollutants from surface runoff. High subsurface flow volumes and saturated conductivities are also likely to limit the residence time of water in the subsurface domain. Based on their hydrologic performance, the key roles of riparian buffers in this landscape are likely to be displacing sediment and nutrient‐generating activities away from streams and stabilizing channel morphology. Copyright © 2006 John Wiley & Sons, Ltd.     

HydroBlocks: a field‐scale resolving land surface model for application over continental extents

Land surface spatial heterogeneity plays a significant role in the water, energy, and carbon cycles over a range of temporal and spatial scales. Until now, the representation of this spatial heterogeneity in land surface models has been limited to over simplistic schemes because of computation and environmental data limitations. This study introduces HydroBlocks – a novel land surface model that represents field‐scale spatial heterogeneity of land surface processes through interacting hydrologic response units (HRUs). HydroBlocks is a coupling between the Noah‐MP land surface model and the Dynamic TOPMODEL hydrologic model. The HRUs are defined by clustering proxies of the drivers of spatial heterogeneity using high‐resolution land data. The clustering mechanism allows for each HRU's results to be mapped out in space, facilitating field‐scale application and validation. The Little Washita watershed in the USA is used to assess HydroBlocks' performance and added benefit from traditional land surface models. A comparison between the semi‐distributed and fully distributed versions of the model suggests that using 1000 HRUs is sufficient to accurately approximate the fully distributed solution. A preliminary evaluation of model performance using available in situ soil moisture observations suggests that HydroBlocks is generally able to reproduce the observed spatial and temporal dynamics of soil moisture. Model performance deficiencies can be primarily attributed to parameter uncertainty. HydroBlocks' ability to explicitly resolve field‐scale spatial heterogeneity while only requiring an increase in computation of one to two orders of magnitude when compared with existing land surface models is encouraging – ensemble field‐scale land surface modelling over continental extents is now possible. Copyright © 2016 John Wiley & Sons, Ltd.     

Determination of bare surface soil moisture from combined temporal evolution of land surface temperature and net surface shortwave radiation

Land surface soil moisture (SSM) is an important variable for hydrological, ecological, and meteorological applications. A multi‐linear model has recently been proposed to determine the SSM content from the combined diurnal evolution of both land surface temperature (LST) and net surface shortwave radiation (NSSR) with the parameters TN (the LST mid‐morning rising rate divided by the NSSR rising rate during the same period) and t_d (the time of daily maximum temperature). However, in addition to the problem that all the coefficients of the multi‐linear model depend on the atmospheric conditions, the model also suffers from the problems of the nonlinearity of TN as a function of the SSM content and the uncertainty of determining the t_d from the diurnal evolution of the LST. To address these problems, a modified multi‐linear model was developed using the logarithm of TN and normalizing t_d by the mid‐morning temperature difference instead of using the TN and t_d. Except for the constant term, the coefficients of all other variables in the modified multi‐linear model proved to be independent of the atmospheric conditions. Using the relevant simulation data, results from the modified multi‐linear model show that the SSM content can be determined with a root mean square error (RMSE) of 0.030m³/m³, provided that the constant term is known or estimated day to day. The validation of the model was conducted using the field measurements at the Langfang site in 2008 in China. A higher correlation is achieved (coefficient of determination: R² = 0.624, RMSE = 0.107m³/m³) between the measured SSM content and the SSM content estimated using the modified multi‐linear model with the coefficients determined from the simulation data. Another experiment is also conducted to estimate the SSM content using the modified model with the constant term calibrated each day by one‐spot measurements at the site. The estimation result has a relatively larger error (RMSE = 0.125m³/m³). Additionally, the uncertainty of the determination of the coefficients is analysed using the field measurements, and the results indicate that the SSM content obtained using the modified model accurately characterizes the surface soil moisture condition. Copyright © 2012 John Wiley & Sons, Ltd.     

Spatio‐temporal distribution of near‐surface and root zone soil moisture at the catchment scale

Soil moisture is highly variable both spatially and temporally. It is widely recognized that improving the knowledge and understanding of soil moisture and the processes underpinning its spatial and temporal distribution is critical. This paper addresses the relationship between near‐surface and root zone soil moisture, the way in which they vary spatially and temporally, and the effect of sampling design for determining catchment scale soil moisture dynamics. In this study, catchment scale near‐surface (0–50 mm) and root zone (0–300 mm) soil moisture were monitored over a four‐week period. Measurements of near‐surface soil moisture were recorded at various resolutions, and near‐surface and root zone soil moisture data were also monitored continuously within a network of recording sensors. Catchment average near‐surface soil moisture derived from detailed spatial measurements and continuous observations at fixed points were found to be significantly correlated (r² = 0·96; P = 0·0063; n = 4). Root zone soil moisture was also found to be highly correlated with catchment average near‐surface, continuously monitored (r² = 0·81; P < 0·0001; n = 26) and with detailed spatial measurements of near‐surface soil moisture (r² = 0·84). The weaker relationship observed between near‐surface and root zone soil moisture is considered to be caused by the different responses to rainfall and the different factors controlling soil moisture for the soil depths of 0–50 mm and 0–300 mm. Aspect is considered to be the main factor influencing the spatial and temporal distribution of near‐surface soil moisture, while topography and soil type are considered important for root zone soil moisture. The ability of a limited number of monitoring stations to provide accurate estimates of catchment scale average soil moisture for both near‐surface and root zone is thus demonstrated, as opposed to high resolution spatial measurements. Similarly, the use of near‐surface soil moisture measurements to obtain a reliable estimate of deeper soil moisture levels at the small catchment scale was demonstrated. Copyright © 2007 John Wiley & Sons, Ltd.     

Modelling snowpack surface temperature in the Canadian Prairies using simplified heat flow models

Three practical schemes for computing the snow surface temperature T_s, i.e. the force–restore method (FRM), the surface conductance method (SCM), and the Kondo and Yamazaki method (KYM), were assessed with respect to T_s retrieved from cloud‐free, NOAA‐AVHRR satellite data for three land‐cover types of the Paddle River basin of central Alberta. In terms of R², the mean T_s, the t‐test and F‐test, the FRM generally simulated more accurate T_s than the SCM and KYM. The bias in simulated T_s is usually within several degrees Celsius of the NOAA‐AVHRR T_s for both the calibration and validation periods, but larger errors are encountered occasionally, especially when T_s is substantially above 0 °C. Results show that the simulated T_s of the FRM is more consistent than that of the SCM, which in turn was more consistent than that of the KYM. This is partly because the FRM considers two aspects of heat conduction into snow, a stationary‐mean diurnal (sinusoidal) temperature variation at the surface coupled to a near steady‐state ground heat flux, whereas the SCM assumes a near steady‐state, simple heat conduction, and other simplifying assumptions, and the KYM does not balance the snowpack heat fluxes by assuming the snowpack having a vertical temperature profile that is linear. Copyright © 2005 John Wiley & Sons, Ltd.     

An approach to using snow areal depletion curves inferred from MODIS and its application to land surface modelling in Alaska

Snowcover areal depletion curves inferred from the moderate resolution imaging spectroradiometer (MODIS) are validated and then applied in NASA's catchment‐based land surface model (CLSM) for numerical simulations of hydrometeorological processes in the Kuparuk River basin (KRB) of Alaska. The results demonstrate that the MODIS snowcover fraction f derived from a simple relationship in terms of the normalized difference snow index compares well with Landsat values over the range 20 ≤ f ≤ 100%. For f < 20%, however, MODIS 500 m subpixel data underestimate the amount of snow by up to 13% compared with Landsat at spatial resolutions of 30 m binned to equivalent 500 m pixels. After a bias correction, MODIS snow areal depletion curves during the spring transition period of 2002 for the KRB exhibit similar features to those derived from surface‐based observations. These results are applied in the CLSM subgrid‐scale snow parameterization that includes a deep and a shallow snowcover fraction. Simulations of the evolution of the snowpack and of freshwater discharge rates for the KRB over a period of 11 years are then analysed with the inclusion of this feature. It is shown that persistent snowdrifts on the arctic landscape, associated with a secondary plateau in the snow areal depletion curves, are hydrologically important. An automated method is developed to generate the shallow and deep snowcover fractions from MODIS snow areal depletion curves. This provides the means to apply the CLSM subgrid‐scale snow parameterization in all watersheds subject to seasonal snowcovers. Improved simulations and predictions of the global surface energy and water budgets are expected with the incorporation of the MODIS snow data into the CLSM. Copyright © 2005 John Wiley & Sons, Ltd.     

Comparing state-of the art near-surface models of a seismic test line from Saudi Arabia

We present a seismic Test Line, provided by Saudi Aramco for various research teams, to highlight a few major challenges in land data processing due to near‐surface anomalies. We discuss state‐of‐the‐art methods used to compensate for shallow distortions, including single‐layer, multilayer, plus/minus, refraction and tomostatics methods. They are a starting point for the new technologies presented in other papers, all dealing with the same challenging data described here. The difficulties on the Test Line are mostly due to the assumption of vertical raypaths, inherent in classical applications of near‐surface correction statics. Even the most detailed velocity/depth model presents difficulties, due to the compleX‐raypath. There is a need for methods which are based on more complex models andtheories.     

Flow threads in surface run‐off: implications for the assessment of flow properties and friction coefficients in soil erosion and hydraulics investigations

Soil surface microtopography produces non‐uniform surface run‐off, in which narrow threads of relatively deep and fast ?ow move within broader, shallower, slower‐moving regions. This kind of ?ow is probably widespread, given that microtopography is itself common. Methods used to record the properties of surface run‐off include grid‐ or transect‐based depth observations, with a single mean ?ow speed derived by calculation from V = Q/WD, and the use of dye timing to estimate velocity, with an effective mean depth calculated from D = Q/WV. Because these methods allow only single, ?ow‐?eld mean values to be derived for V or D, neither is well suited to non‐uniform ?ows. The use of depth data to derive a ?ow‐?eld mean V furthermore implicitly applies area weighting to the depth data; likewise, the use of dye speeds for V inherently overestimates mean V because dye dominantly follows the faster ?ow threads. The associated errors in derived parameters such as friction coef?cients are not readily quanti?ed and appear not to have been addressed previously. New ?eld experiments made on untilled soil surfaces in arid western NSW, Australia, explore these circumstances and the implications for deriving meaningful measures of ?ow properties, including friction coef?cients. On surfaces deliberately chosen for their very subtle microtopography, average thread velocities are shown to be commonly 2·5 times greater than the ?ow‐?eld mean, and locally 6–7 times greater. On the other hand, non‐thread ?ow speeds lie below the ?ow‐?eld mean, on average reaching only 84 per cent of this value, and often considerably less. Flow‐?eld means conceal the existence of regions of the ?ow ?eld whose properties are statistically distinct. Results con?rm that a reliance on ?ow‐?eld average depths yields estimates of friction coef?cients that are biased toward the shallower, high‐roughness parts of the ?ow, while if dye speeds are relied upon the results are biased toward the deeper, smoother threads of ?ow. A new approach to the evaluation of friction coef?cients in non‐uniform ?ows is advanced, involving the determination of separate coef?cients for threads and non‐thread zones of the ?ow ?eld. In contrast, ?ow‐?eld friction coef?cients as they are customarily derived in run‐off plot experiments subsume these distinct coef?cients in proportions that are generally unknown. The value of such coef?cients is therefore questionable. Copyright © 2004 John Wiley & Sons, Ltd.     

Simplified modelling of areal average infiltration at the hillslope scale

Two models for estimating expected areal‐average infiltration rate, ī, at the hillslope scale are presented. The first relies upon the condition of a negligible infiltration of surface water running downslope (run‐on process) into a previous heterogeneous soil. It is an adapted version of an earlier semi‐analytical model. The second incorporates the run‐on process and is based on a lumped approach that uses an effective saturated hydraulic conductivity. This latter was parameterized in terms of the main characteristics of rainfall and soil. Both the models were tested by comparison with the results carried out by Monte‐Carlo simulations over different soil types. It was found that the first model simulated ī with maximum errors in magnitude typically less than 10%. The second model provided similar errors in the total volume of overland flow, and the rising limb of the hydrograph experienced a distortion. Lastly, satisfactory results were obtained by comparing the model without run‐on with an empirical approach particularly accurate for fine‐textured soils. Copyright © 2002 John Wiley & Sons, Ltd.     

Role of fluids in surface deformation caused by the 1999 Chi‐Chi earthquake in Taiwan

The 1999 Chi‐Chi earthquake significantly altered the landscape of central Taiwan. Surface deformation produced by the earthquake along the trace of the Chelungpu thrust can be classified into two styles: (1) uplift without significant surface rupture, and (2) uplift accompanied by surface rupture. Here we examine areas that exhibited the first style of deformation (e.g. Wufeng). Seismic stress at the time of the main shock may have been relieved by high pore‐fluid pressure in a 300‐m‐thick sand and gravel aquifer. Along the thrust fault, frictional heating of these sediments resulted in thermal expansion and an increase in pore‐fluid pressure. High pore‐fluid pressure damped seismic‐wave energy and enhanced intergranular slips of unconsolidated sandy and gravel sediments, which were possibly assisted by sulphuric acid corrosion, leading to a high sulphate content in the groundwater (c. 70 mg L^?1). These changes permitted surface folding and terrace‐style uplifting to occur without significant rupture. In contrast, other areas in which the second style of deformation is dominant (e.g. Fengyuen‐Shihkang) have thin (0–10 m) sand and gravel deposits and lower concentrations of sulphate (c. 30 mg L^?1) in groundwater. In these areas, sediments were heated but not sufficiently to produce significant thermal expansion and increase in pore‐fluid pressure; accumulation of stress in these locations led to rupture at the ground surface, with the formation of steep fault scarps. The areas exhibiting the first deformation style are characterized by the presence of high pore‐fluid pressure, frictional heat conduction, and possibly chemical corrosion related to sulphuric acid attack and formation of sulphate, in contrast to those involving significant uplift and surface rupture. The areal distribution of these two surface deformation styles suggests that the aforementioned fluid‐related subsurface processes may have altered the characteristics of sediments and caused diverse responses to the quake. Copyright © 2002 John Wiley & Sons, Ltd.     

Salt marsh surface survives true‐to‐scale simulated storm surges

A full‐scale controlled experiment was conducted on an excavated and re‐assembled coastal wetland surface, typical of floristically diverse northwest European saltmarsh. The experiment was undertaken with true‐to‐scale water depths and waves in a large wave flume, in order to assess the impact of storm surge conditions on marsh surface soils, initially with three different plant species and then when this marsh canopy had been mowed. The data presented suggests a high bio‐geomorphological resilience of salt marshes to vertical sediment removal, with less than 0.6 cm average vertical lowering in response to a sequence of simulated storm surge conditions. Both organic matter content and plant species exerted an important influence on both the variability and degree of soil surface stability, with surfaces covered by a flattened canopy of the salt marsh grass Puccinellia experiencing a lower and less variable elevation loss than those characterized by Elymus or Atriplex that exhibited considerable physical damage through stem folding and breakage. Copyright © 2015 John Wiley & Sons, Ltd.     

Innovative analysis and use of high‐resolution DTMs for quantitative interrogation of Earth‐surface processes

This is the era of digital landscapes; the widespread availability of powerful sensing technologies has revolutionized the way it is possible to interrogate landscapes in order to understand the processes sculpting them. Vastly greater areas have now been acquired at ‘high resolution’: currently tens of metres globally to millimetric precision and accuracy locally. This permits geomorphic features to be visualized and analysed across the scales at which Earth‐surface processes operate. Especially exciting is the capturing of process dynamics in repeated surveying, which will only become more important with low‐cost accessible data generation through techniques such as Structure from Motion (SfM). But the key challenge remains; to interpret high resolution Digital Terrain Models (DTMs), particularly by extracting geomorphic features in robust and objective ways and then linking the observed features to the underlying physical processes. In response to the new data and challenges, recent years have seen improved processing of raw data into DTMs, development of data fusion techniques, novel quantitative analysis of topographic data, and innovative geomorphological mapping. The twelve papers collected in this volume sample this progress in interrogating Earth‐surface processes through the analysis of DTMs. They cover a wide range of disciplines and spatio‐temporal scales, from landslide prone landscapes, to agriculturally modified regions, to mountainous landscapes, and coastal zones. They all, however, showcase the quantitative exploitation of information contained in high‐resolution topographic data that we believe will underpin the improvement of our understanding of many elements of Earth‐surface processes. Most of the papers introduced here were first presented in a conference session at the European Geosciences Union General Assembly in 2011. Copyright © 2014 John Wiley & Sons, Ltd.