Simulated effect of a forest road on near‐surface hydrologic response and slope stability

Forest management practices often result in significant changes to hydrologic and geomorphic responses at or near the earth's surface. A well‐known, but not fully tested, hypothesis in hillslope hydrology[sol ]geomorphology is that a near‐surface permeability contrast, caused by the surface compaction associated with forest roads, can result in diverted subsurface flow paths that produce increased up‐slope pore pressures and slope failure. The forest road focused on in this study is located in a steep forested, zero‐order catchment within the H. J. Andrews Experimental Forest (Oregon). A three‐phase modelling effort was employed to test the aforementioned hypothesis: (i) two‐dimensional (vertical slice), steady‐state, heterogeneous, saturated subsurface flow simulations at the watershed scale for establishing the boundary conditions for the catchment‐scale boundary‐value problem in (ii); (ii) two‐dimensional (vertical slice), transient, heterogeneous, variably saturated subsurface flow simulations at the catchment scale for estimating near‐surface hydrologic response and pore pressure distributions; and (iii) slope stability analyses, using the infinite slope approach, driven by the pore pressure distributions simulated in (ii), for assessing the impact of the forest road. Both observed and hypothetical rainfall events are used to drive the catchment‐scale simulations. The results reported here support the hypothesis that a forest road can have an effect on slope stability. The permeability contrast associated with the forest road in this study led to a simulated altering of slope‐parallel subsurface flow with increased pore pressures up‐slope of the road and, for a large rainfall event, a slope failure prediction. Copyright © 2005 John Wiley & Sons, Ltd.     

Simulated effect of a forest road on near‐surface hydrologic response: redux

In the work reported here the comprehensive physics‐based Integrated Hydrology Model (InHM) was employed to conduct both three‐ and two‐dimensional (3D and 2D) hydrologic‐response simulations for the small upland catchment known as C3 (located within the H. J. Andrews Experimental Forest in Oregon). Results from the 3D simulations for the steep unchannelled C3 (i) identify subsurface stormflow as the dominant hydrologic‐response mechanism and (ii) show the effect of the down‐gradient forest road on both the surface and subsurface flow systems. Comparison of the 3D results with the 2D results clearly illustrates the importance of convergent subsurface flow (e.g. greater pore‐water pressures in the hollow of the catchment for the 3D scenario). A simple infinite‐slope model, driven by subsurface pore‐water pressures generated from the 3D and 2D hydrologic‐response simulations, was employed to estimate slope stability along the long‐profile of the C3 hollow axis. As expected, the likelihood of slope failure is underestimated for the lower pore pressures from the 2D hydrologic‐response simulation compared, in a relative sense, to the higher pore pressures from the 3D hydrologic response simulation. The effort reported herein provides a firm quantitative foundation for generalizing the effects that forest roads can have on near‐surface hydrologic response and slope stability at the catchment scale. Copyright © 2006 John Wiley & Sons, Ltd.     

Effects of hillslope topography on hydrological responses in a weathered granite mountain,Japan: comparison of the runoff response between the valley‐head and the side slope

To evaluate the effects of hillslope topography on storm runoff in a weathered granite mountain, discharge rate, soil pore water pressures, and water chemistry were observed on two types of hillslope: a valley‐head (a concave hillslope) and a side slope (a planar hillslope). Hydrological responses on the valley‐head and side slope reflected their respective topographic characteristics and varied with the rainfall magnitude. During small rainfall events (<35 mm), runoff from the side slope occurred rapidly relative to the valley‐head. The valley‐head showed little response in storm runoff. As rainfall amounts increased (35–60 mm), the valley‐head yielded a higher flow relative to the side slope. For large rainfall events (>60 mm), runoff from both hillslopes increased with rainfall, although that from the valley‐head was larger than that from the side slope. The differences in the runoff responses were caused by differences in the roles of lower‐slope soils and the convergence of the hillslope. During small rainfall events, the side slope could store little water; in contrast, all rainwater could be stored in the soils at the valley‐head hollow. As the amount of rainfall increased, the subsurface saturated area of the valley‐head extended from the bottom to the upper portion of the slope, with the contributions of transient groundwater via lateral preferential flowpaths due to the high concentration of subsurface water. Conversely, saturated subsurface flow did not contribute to runoff responses, and the subsurface saturated area at the side slope did not extend to the upper slope for the same storm size. During large rainfall events, expansion of the subsurface saturated area was observed in both hillslopes. Thus, differences in the concentration of subsurface water, reflecting hillslope topography, may create differences in the extension of the subsurface saturated area, as well as variability in runoff responses. Copyright © 2007 John Wiley & Sons, Ltd.     

Modeling Variably Saturated Multispecies Reactive Groundwater Solute Transport with MODFLOW‐UZF and RT3D

A numerical model was developed that is capable of simulating multispecies reactive solute transport in variably saturated porous media. This model consists of a modified version of the reactive transport model RT3D (Reactive Transport in 3 Dimensions) that is linked to the Unsaturated‐Zone Flow (UZF1) package and MODFLOW. Referred to as UZF‐RT3D, the model is tested against published analytical benchmarks as well as other published contaminant transport models, including HYDRUS‐1D, VS2DT, and SUTRA, and the coupled flow and transport modeling system of CATHY and TRAN3D. Comparisons in one‐dimensional, two‐dimensional, and three‐dimensional variably saturated systems are explored. While several test cases are included to verify the correct implementation of variably saturated transport in UZF‐RT3D, other cases are included to demonstrate the usefulness of the code in terms of model run‐time and handling the reaction kinetics of multiple interacting species in variably saturated subsurface systems. As UZF1 relies on a kinematic‐wave approximation for unsaturated flow that neglects the diffusive terms in Richards equation, UZF‐RT3D can be used for large‐scale aquifer systems for which the UZF1 formulation is reasonable, that is, capillary‐pressure gradients can be neglected and soil parameters can be treated as homogeneous. Decreased model run‐time and the ability to include site‐specific chemical species and chemical reactions make UZF‐RT3D an attractive model for efficient simulation of multispecies reactive transport in variably saturated large‐scale subsurface systems.     

Modeling Variably Saturated Subsurface Solute Transport with MODFLOW‐UZF and MT3DMS

The MT3DMS groundwater solute transport model was modified to simulate solute transport in the unsaturated zone by incorporating the unsaturated‐zone flow (UZF1) package developed for MODFLOW. The modified MT3DMS code uses a volume‐averaged approach in which Lagrangian‐based UZF1 fluid fluxes and storage changes are mapped onto a fixed grid. Referred to as UZF‐MT3DMS, the linked model was tested against published benchmarks solved analytically as well as against other published codes, most frequently the U.S. Geological Survey's Variably‐Saturated Two‐Dimensional Flow and Transport Model. Results from a suite of test cases demonstrate that the modified code accurately simulates solute advection, dispersion, and reaction in the unsaturated zone. Two‐ and three‐dimensional simulations also were investigated to ensure unsaturated‐saturated zone interaction was simulated correctly. Because the UZF1 solution is analytical, large‐scale flow and transport investigations can be performed free from the computational and data burdens required by numerical solutions to Richards' equation. Results demonstrate that significant simulation runtime savings can be achieved with UZF‐MT3DMS, an important development when hundreds or thousands of model runs are required during parameter estimation and uncertainty analysis. Three‐dimensional variably saturated flow and transport simulations revealed UZF‐MT3DMS to have runtimes that are less than one tenth of the time required by models that rely on Richards' equation. Given its accuracy and efficiency, and the wide‐spread use of both MODFLOW and MT3DMS, the added capability of unsaturated‐zone transport in this familiar modeling framework stands to benefit a broad user‐ship.     

PHT3D‐UZF: A Reactive Transport Model for Variably‐Saturated Porous Media

A modified version of the MODFLOW/MT3DMS‐based reactive transport model PHT3D was developed to extend current reactive transport capabilities to the variably‐saturated component of the subsurface system and incorporate diffusive reactive transport of gaseous species. Referred to as PHT3D‐UZF, this code incorporates flux terms calculated by MODFLOW's unsaturated‐zone flow (UZF1) package. A volume‐averaged approach similar to the method used in UZF‐MT3DMS was adopted. The PHREEQC‐based computation of chemical processes within PHT3D‐UZF in combination with the analytical solution method of UZF1 allows for comprehensive reactive transport investigations (i.e., biogeochemical transformations) that jointly involve saturated and unsaturated zone processes. Intended for regional‐scale applications, UZF1 simulates downward‐only flux within the unsaturated zone. The model was tested by comparing simulation results with those of existing numerical models. The comparison was performed for several benchmark problems that cover a range of important hydrological and reactive transport processes. A 2D simulation scenario was defined to illustrate the geochemical evolution following dewatering in a sandy acid sulfate soil environment. Other potential applications include the simulation of biogeochemical processes in variably‐saturated systems that track the transport and fate of agricultural pollutants, nutrients, natural and xenobiotic organic compounds and micropollutants such as pharmaceuticals, as well as the evolution of isotope patterns.     

How hydrological factors initiate instability in a model sandy slope

Knowledge of the mechanisms of rain‐induced shallow landslides can improve the prediction of their occurrence and mitigate subsequent sediment disasters. Here, we examine an artificial slope's subsurface hydrology and propose a new slope stability analysis that includes seepage force and the down‐slope transfer of excess shear forces. We measured pore water pressure and volumetric water content immediately prior to a shallow landslide on an artificial sandy slope of 32°: The direction of the subsurface flow shifted from downward to parallel to the slope in the deepest part of the landslide mass, and this shift coincided with the start of soil displacement. A slope stability analysis that was restricted to individual segments of the landslide mass could not explain the initiation of the landslide; however, inclusion of the transfer of excess shear forces from up‐slope to down‐slope segments improved drastically the predictability. The improved stability analysis revealed that an unstable zone expanded down‐slope with an increase in soil water content, showing that the down‐slope soil initially supported the unstable up‐slope soil; destabilization of this down‐slope soil was the eventual trigger of total slope collapse. Initially, the effect of apparent soil cohesion was the most important factor promoting slope stability, but seepage force became the most important factor promoting slope instability closer to the landslide occurrence. These findings indicate that seepage forces, controlled by changes in direction and magnitude of saturated and unsaturated subsurface flows, may be the main cause of shallow landslides in sandy slopes. Copyright © 2013 John Wiley & Sons, Ltd.     

Two‐dimensional and three‐dimensional computational models in hydrodynamic and morphodynamic reconstructions of a river bend: sensitivity and functionality

This study assesses hydrodynamic and morphodynamic model sensitivity and functionality in a curved channel. The sensitivity of a depth‐averaged model to user‐defined parameters (grain size, roughness, transverse bed slope effect, transport relations and secondary flow) is tested. According to the sensitivity analysis, grain size, transverse bed slope effect and sediment transport relations are critical to simulated meander bend morphodynamics. The parametrization of grain size has the most remarkable effect: field‐based grain size parametrization is necessary in a successful morphodynamic reconstruction of a meander bend. The roughness parametrization method affects the distribution of flow velocities and therefore also morphodynamics. The combined effect of various parameters needs further research. Two‐dimensional (2D) and three‐dimensional (3D) reconstructions of a natural meander bend during a flood event are assessed against field measurements of acoustic Doppler current profiler and multi‐temporal mobile laser scanning data. The depth‐averaged velocities are simulated satisfactorily (differences from acoustic Doppler current profiler velocities 5–14%) in both 2D and 3D simulations, but the advantage of the 3D hydrodynamic model is unquestionable because of its ability to model vertical and near‐bed flows. The measured and modelled near‐bed flow, however, differed notably from each other's, the reason of which was left open for future research. It was challenging to model flow direction beyond the apex. The 3D flow features, which also affected the distribution of the bed shear stress, seem not to have much effect on the predicted morphodynamics: the 2D and 3D morphodynamic reconstructions over the point bar resembled each other closely. Although common features between the modelled and measured morphological changes were also found, some specific changes that occurred were not evident in the simulation results. Our results show that short‐term, sub‐bend scale morphodynamic processes of a natural meander bend are challenging to model, which implies that they are affected by factors that have been neglected in the simulations. The modelling of short‐term morphodynamics in natural curved channel is a challenge that requires further study. Copyright © 2014 John Wiley & Sons, Ltd.     

Formative conditions and sedimentary structures of sandy 3D antidunes: an application of the gravel step‐pool model to fine‐grained sand in an experimental flume

Antidunes and their sedimentary structures can be useful in reconstructing paleo‐hydraulic conditions, especially for large discharge events. However, three‐dimensional (3D) antidunes in sand‐sized sediments have not yet been studied extensively, as compared to either two‐dimensional (2D) antidunes or antidunes in gravel‐sized sediments. In this study, we estimated formative conditions of gravel step‐pool morphologies and applied them to the formation of 3D antidunes over a sand bed. Formative conditions are expressed in terms of a relationship between the water discharge per unit width and the bed slope. Flume experiments demonstrated that 3D mound‐like antidune configurations and their associated internal sedimentary structures could be preserved. Internal sedimentary structures were characterized by shallow lens‐like structures whose bases were erosional. Although gently‐dipping concave‐upward lamination was dominant, convex‐upward lamination was occasionally observed. The dimensions of lenticular lamina‐sets can be used to estimate antidune geometry. Thus if 3D antidunes can be interpreted in the stratigraphic record, it is possible to estimate the paleo‐hydraulic parameters such as water discharge and bed slope more precisely than previously. Copyright © 2010 John Wiley & Sons, Ltd.     

A low‐dimensional physically based model of hydrologic control of shallow landsliding on complex hillslopes

Hillslopes have complex three‐dimensional shapes that are characterized by their plan shape, profile curvature of surface and bedrock, and soil depth. To investigate the stability of complex hillslopes (with different slope curvatures and plan shapes), we combine the hillslope‐storage Boussinesq (HSB) model with the infinite slope stability method. The HSB model is based on the continuity and Darcy equations expressed in terms of storage along the hillslope. Solutions of the HSB equation account explicitly for plan shape by introducing the hillslope width function and for profile curvature through the bedrock slope angle and the hillslope soil depth function. The presented model is composed of three parts: a topography model conceptualizing three‐dimensional soil mantled landscapes, a dynamic hydrology model for shallow subsurface flow and water table depth (HSB model) and an infinite slope stability method based on the Mohr–Coulomb failure law. The resulting hillslope‐storage Boussinesq stability model (HSB‐SM) is able to simulate rain‐induced shallow landsliding on hillslopes with non‐constant bedrock slope and non‐parallel plan shape. We apply the model to nine characteristic hillslope types with three different profile curvatures (concave, straight, convex) and three different plan shapes (convergent, parallel, divergent). In the presented model, the unsaturated storage has been calculated based on the unit head gradient assumption. To relax this assumption and to investigate the effect of neglecting the variations of unsaturated storage on the assessment of slope stability in the transient case, we also combine a coupled model of saturated and unsaturated storage and the infinite slope stability method. The results show that the variations of the unsaturated zone storage do not play a critical role in hillslope stability. Therefore, it can be concluded that the presented dynamic slope stability model (HSB‐SM) can be used safely for slope stability analysis on complex hillslopes. Our results show that after a certain period of rainfall the convergent hillslopes with concave and straight profiles become unstable more quickly than others, whilst divergent convex hillslopes remain stable (even after intense rainfall). In addition, the relation between subsurface flow and hillslope stability has been investigated. Our analyses show that the minimum safety factor (FS) occurs when the rate of subsurface flow is a maximum. In fact, by increasing the subsurface flow, stability decreases for all hillslope shapes. Copyright © 2008 John Wiley & Sons, Ltd.     

Analysis of geomorphologic and hydrological characteristics in watershed saturated areas using topographic‐index threshold and geomorphology‐based runoff model

The objective of this paper is to investigate the variation of geomorphology and runoff characteristics in saturated areas under different partial contributing area (PCA) conditions. Geomorphologic information and hydrologic records from two mid‐size watersheds in northern Taiwan were selected for analysis. The PCA ratio in the watershed during a storm was assumed equal to the ratio of the surface‐flow volume to the direct runoff volume from measured hydrologic data. The extents of PCA regions were then determined by using a topographic‐index threshold. Consequently, the geomorphologic factors in saturated and unsaturated areas could be calculated using a digital elevation model, and these factors could then be linked to a geomorphology‐based IUH model for runoff simulation, which can consider both the surface‐ and subsurface‐flow processes in saturated and unsaturated areas, respectively. The results show that geomorphologic characteristics in the saturated areas vary significantly with different PCA ratios especially for higher order streams. A large PCA ratio results in a sharp hydrograph because the quick surface flow dominates the runoff process, whereas the hydrologic response in a low PCA case is dominated by the delayed subsurface flow. Copyright © 2007 John Wiley & Sons, Ltd.     

An estimation method for effective stress changes in a reservoir from 4D seismics data

Advances in seismics acquisition and processing and the widespread use of 4D seismics have made available reliable production‐induced subsurface deformation data in the form of overburden time‐shifts. Inversion of these data is now beginning to be used as an aid to the monitoring of a reservoir's effective stress. Past solutions to this inversion problem have relied upon analytic calculations for an unrealistically simplified subsurface, which can lead to uncertainties. To enhance the accuracy of this approach, a method based on transfer functions is proposed in which the function itself is calibrated using numerically generated overburden strain deformation calculated for a small select group of reference sources. This technique proves to be a good compromise between the faster but more accurate history match of the overburden strain using a geomechanical simulator and the slower, less accurate analytic method. Synthetic tests using a coupled geomechanical and fluid flow simulator for the South Arne field confirm the efficacy of the method. Application to measured time‐shifts from observed 4D seismics indicates compartmentalization in the Tor reservoir, more heterogeneity than is currently considered in the simulation model and moderate connectivity with the overlying Ekofisk formation.     

Evaluation of model complexity and space–time resolution on the prediction of long‐term soil salinity dynamics,western San Joaquin Valley,California

The numerical simulation of long‐term large‐scale (field to regional) variably saturated subsurface flow and transport remains a computational challenge, even with today's computing power. Therefore, it is appropriate to develop and use simplified models that focus on the main processes operating at the pertinent time and space scales, as long as the error introduced by the simpler model is small relative to the uncertainties associated with the spatial and temporal variation of boundary conditions and parameter values. This study investigates the effects of various model simplifications on the prediction of long‐term soil salinity and salt transport in irrigated soils. Average root‐zone salinity and cumulative annual drainage salt load were predicted for a 10‐year period using a one‐dimensional numerical flow and transport model (i.e. UNSATCHEM) that accounts for solute advection, dispersion and diffusion, and complex salt chemistry. The model uses daily values for rainfall, irrigation, and potential evapotranspiration rates. Model simulations consist of benchmark scenarios for different hypothetical cases that include shallow and deep water tables, different leaching fractions and soil gypsum content, and shallow groundwater salinity, with and without soil chemical reactions. These hypothetical benchmark simulations are compared with the results of various model simplifications that considered (i) annual average boundary conditions, (ii) coarser spatial discretization, and (iii) reducing the complexity of the salt‐soil reaction system. Based on the 10‐year simulation results, we conclude that salt transport modelling does not require daily boundary conditions, a fine spatial resolution, or complex salt chemistry. Instead, if the focus is on long‐term salinity, then a simplified modelling approach can be used, using annually averaged boundary conditions, a coarse spatial discretization, and inclusion of soil chemistry that only accounts for cation exchange and gypsum dissolution–precipitation. We also demonstrate that prediction errors due to these model simplifications may be small, when compared with effects of parameter uncertainty on model predictions. The proposed model simplifications lead to larger time steps and reduced computer simulation times by a factor of 1000. Copyright © 2006 John Wiley & Sons, Ltd.     

Nitrate subsurface transport and losses in response to its initial distributions in sloped soils: An experimental and modelling study

Transport and losses of nitrate from sloped soils are closely linked to nitrogen fertilizer management. Previous studies have always focused on different types of fertilizer applications and rarely analysed various initial nitrate distributions as a result of nitrogen fertilizer applications. Under certain conditions, both subsurface lateral saturated flow and vertical leaching dominate nitrate losses. Soil tank experiments and HYDRUS‐2D modelling were used to better understand the subsurface nitrate transport and losses through lateral saturated flow and vertical leaching under various initial nitrate distributions. Low (L: 180 mg L^?1), normal (N: 350 mg L^?1), and high (H: 500 mg L^?1) nitrate concentrations were used in five different distributions (NNNN, NLLN, LHHL, LNLN, and HNHN) along the slope of the tank. The first two treatments (NNNN and NLLN) were analysed both experimentally and numerically. Experiments were conducted under 12 rainfall events at intervals of 3 days. The HYDRUS‐2D model was calibrated and validated against the experimental data and demonstrated good model performance. The other three treatments (LHHL, LNLN, and HNHN) were investigated using the calibrated model. Nitrate concentrations in purple sloped soils declined exponentially with time under intermittent rainfalls, predominantly in the upper soil layers. Non‐uniform initial nitrate distributions contributed to larger differences between four locations along the slope in deeper soil layers. The non‐uniform nitrate distribution either enhanced or reduced decreases in nitrate concentrations in areas with higher or lower initial nitrate concentrations, respectively. Higher nitrate concentrations at the slope foot and along the slope were reduced mainly by lateral flow and vertical leaching, respectively. Increasing nitrogen application rates increased subsurface nitrate losses. Mean subsurface lateral nitrate fluxes were twice as large as mean vertical leaching nitrate fluxes. However, due to longer leaching durations, total nitrate losses due to vertical leaching were comparable with those due to lateral flow, which indicated comparable environmental risks to surface waters and groundwater.     

Joint interpretation of hydrological and geophysical data: electrical resistivity tomography results from a process hydrological research site in the Black Forest Mountains,Germany

The use of electrical resistivity tomography (ERT; non‐intrusive geophysical technique) was assessed to identify the hydrogeological conditions at a surface water/groundwater test site in the southern Black Forest, Germany. A total of 111 ERT transects were measured, which adopted electrode spacings from 0·5 to 5 m as well as using either Wenner or dipole‐dipole electrode arrays. The resulting two‐dimensional (2D) electrical resistivity distributions are related to the structure and water content of the subsurface. The images were interpreted with respect to previous classical hillslope hydrological investigations within the same research basin using both tracer methods and groundwater level observations. A raster‐grid survey provided a quasi 3D resistivity pattern of the floodplain. Strong structural heterogeneity of the subsurface could be demonstrated, and (non)connectivities between surface and subsurface bodies were mapped. Through the spatial identification of likely flow pathways and source areas of runoff, the deep groundwater within the steeper valley slope seems to be much more connected to runoff generation processes within the valley floodplain than commonly credited in such environmental circumstances. Further, there appears to be no direct link between subsurface water‐bodies adjacent to the stream channel. Deep groundwater sources are also able to contribute towards streamflow from exfiltration at the edge of the floodplain as well as through the saturated areas overlying the floodplain itself. Such exfiltrated water then moves towards the stream as channelized surface flow. These findings support previous tracer investigations which showed that groundwater largely dominates the storm hydrograph of the stream, but the source areas of this component were unclear without geophysical measurements. The work highlighted the importance of using information from previous, complementary hydrochemical and hydrometric research campaigns to better interpret the ERT measurements. On the other hand, the ERT can provide a better spatial understanding of existing hydrochemical and hydrometric data. Copyright © 2009 John Wiley & Sons, Ltd.     

A low‐dimensional hillslope‐based catchment model for layered groundwater flow

Despite the strong interaction between surface and subsurface waters, groundwater flow representation is often oversimplified in hydrological models. For instance, the interplay between local or shallow aquifers and deeper regional‐scale aquifers is typically neglected. In this work, a novel hillslope‐based catchment model for the simulation of combined shallow and deep groundwater flow is presented. The model consists of the hillslope‐storage Boussinesq (hsB) model representing shallow groundwater flow and an analytic element (AE) model representing deep regional groundwater flow. The component models are iteratively coupled via a leakage term based on Darcy's law, representing delayed recharge to the regional aquifer through a low conductivity layer. Simulations on synthetic single hillslopes and on a two‐hillslope open‐book catchment are presented, and the results are compared against a benchmark three‐dimensional Richards equation model. The impact of hydraulic conductivity, hillslope plan geometry (uniform, convergent, divergent), and hillslope inclination (0.2%, 5%, and 30%) under drainage and recharge conditions are examined. On the single hillslopes, good matches for heads, hydrographs, and exchange fluxes are generally obtained, with the most significant differences in outflows and heads observed for the 30% slope and for hillslopes with convergent geometry. On the open‐book catchment, cumulative outflows are overestimated by 1–4%. Heads in the confined and unconfined aquifers are adequately reproduced throughout the catchment, whereas exchange fluxes are found to be very sensitive to the hillslope drainable porosity. The new model is highly efficient computationally compared to the benchmark model. The coupled hsB/AE model represents an alternative to commonly used groundwater flow representations in hydrological models, of particular appeal when surface–subsurface exchanges, local aquifer–regional aquifer interactions, and low flows play a key role in a watershed's dynamics. Copyright © 2011 John Wiley & Sons, Ltd.     

Recent Landslides in Western Hokkaido,Japan

—Western Hokkaido, the northernmost island of Japan, is prone to landsliding due to geologic, geomorphologic and climatic change. From 1985 to 1997, many rapid large-scale landslides occurred in western Hokkaido, several of which are reviewed in this paper. The 1988 Kamaya Slide, the 1991 Tachimachi-misaki Slide (which was preceded by the 1985 Orito Slide), the 1993 Okushiri-Port Slide, and the 1994 Motochi Earthflow are described. Finally, two sea-cliff rockfalls are also described.¶The Okushiri-port Slide and the Toyohama Tunnel Rockfall claimed 29 and 20 lives, respectively. Except for the Okushiri-port Slide, which was induced by Hokkaido Nansei-oki Earthquake, most of the slides and rockfalls were probably related to geological structures, such as gentle-dipping strata interbedded with clayey tuffs and were triggered by long-duration and/or high precipitation, causing increased groundwater levels and/or high water pressures.     

Effects of the sea floor topography on the 1D inversion of time‐domain marine controlled source electromagnetic data

Time‐domain marine controlled source electromagnetic methods have been used successfully for the detection of resistive targets such as hydrocarbons, gas hydrate, or marine groundwater aquifers. As the application of time‐domain marine controlled source electromagnetic methods increases, surveys in areas with a strong seabed topography are inevitable. In these cases, an important question is whether bathymetry information should be included in the interpretation of the measured electromagnetic field or not. Since multi‐dimensional inversion is still not common in time‐domain marine controlled source electromagnetic methods, bathymetry effects on the 1D inversion of single‐offset and multi‐offset joint inversions of time‐domain controlled source electromagnetic methods data are investigated. We firstly used an adaptive finite element algorithm to calculate the time‐domain controlled source electromagnetic methods responses of 2D resistivity models with seafloor topography. Then, 1D inversions are applied on the synthetic data derived from marine resistivity models, including the topography in order to study the possible topography effects on the 1D interpretation. To evaluate the effects of topography with various steepness, the slope angle of the seabed topography is varied in the synthetic modelling studies for deep water (air interaction is absent or very weak) and shallow water (air interaction is dominant), respectively. Several different patterns of measuring configurations are considered, such as the systems adopting nodal receivers and the bottom‐towed system. According to the modelling results for deep water when air interaction is absent, the 2D topography can distort the measured electric field. The distortion of the data increases gradually with the enlarging of the topography's slope angle. In our test, depending on the configuration, the seabed topography does not affect the 1D interpretation significantly if the slope angle is less or around 10°. However, if the slope angle increases to 30° or more, it is possible that significant artificial layers occur in inversion results and lead to a wrong interpretation. In a shallow water environment with seabed topography, where the air interaction dominates, it is possible to uncover the true subsurface resistivity structure if the water depth for the 1D inversion is properly chosen. In our synthetic modelling, this scheme can always present a satisfactory data fit in the 1D inversion if only one offset is used in the inversion process. However, the determination of the optimal water depth for a multi‐offset joint inversion is challenging due to the various air interaction for different offsets.     

Estimates of flow resistance and eddy viscosity coefficients for 2D modelling on vegetated floodplains

The problem of quantifying the effects of flexible plants on flow resistance and eddy viscosity by vegetated floodplains is first addressed with a one‐dimensional (1D) approximation based upon the so‐called lateral distribution method. The estimates so obtained are then tested with two‐dimensional (2D) numerical simulations based on the full shallow water equations through the use of the computational code Telemac‐2D. Data obtained on a physical model of the Besòs River (Spain), whose floodplains were covered with plastic ornamental plants to mimic the effect of flexible vegetation, is used for the validation of the numerical results. Additionally, the values of flow resistance estimated numerically with the 1D and 2D simulations are compared with values obtained in a rectangular flume under flow conditions (slope, water depth and artificial lining) similar to those used on the reduced model. It is then established that as more physical mechanisms are included in the mathematical model used to study the problem, the ratio between the floodplain and the main channel flow resistance coefficient increases. The approach demonstrates that whenever enough flow data is available, the lateral distribution method delivers values of flow resistance and eddy viscosity which are highly consistent with 2D numerical modelling. This finding could mean considerable savings in the burdensome task of specifying flow resistance and turbulence dissipation values for 2D modelling of large compound channel systems. Copyright © 2004 John Wiley & Sons, Ltd.     

Three‐dimensional dynamic response of underground openings in saturated rock masses

In this paper a three‐dimensional dynamic analysis program for saturated porous rocks and soils is presented. The theoretical formulations incorporated in the proposed computer program are the extension of Biot's two‐phase theory to non‐linear region. Numerical study using the program is carried out to demonstrate the three‐dimensional features of the program and to show the effects of pore water on the dynamic response of underground openings in saturated rock masses. It is shown that underground openings in saturated porous media could be significantly more vulnerable to the potential damages associated with high motions and shear failure than those in dry media. Copyright © 2001 John Wiley & Sons, Ltd.