Late Quaternary kinematics of the Pallatanga strike-slip fault (Central Ecuador) from topographic measurements of displaced morphological features

The northeast-trending Pallatanga right-lateral strike-slip fault runs across the Western Cordillera connecting N50E-N70E trending normal faults in the Gulf of Guayaquil with N-S reverse faults in the Interandean Depression. Over most of its length, the fault trace has been partly obscured by erosional processes and can be inferred in the topography only at the large scale. Only the northern fault segment, which follows the upper Rio Pangor valley at elevations above 3600 m, is prominent in the morphology. Valleys and ridges cut and offset by the fault provide an outstanding record of right-lateral cumulative fault displacement. The fault geometry and kinematics of this particular fault segment can be determined from detailed topographic levellings. The fault strikes N30E and dips 75 to the NW. Depending on their size and nature, transverse morphological features such as tributaries of the Rio Pangor and intervening ridges, reveal right-lateral offsets which cluster around 27 ± 11m, 41.5 ± 4 m, 590 ± 65 m and 960 ± 70 m. The slip vector deduced from the short-term offsets shows a slight reverse component with a pitch of about 11.5 SW. The 41.5 ± 4 m displacements are assumed to be coeval with the last glacial termination, yielding a mean Holocene slip-rate of 2.9- 4.6 mm yr⁻¹. Assuming a uniform slip rate on the fault in the long term, the 27 m offset appears to correlate with an identified middle Holocene morphoclimatic event, and the long term offsets of 590 m and 960 m coincide with the glacial terminations at the beginning of the last two interglacial periods.     

Mapping landslide susceptibility from small datasets: A case study in the Pays de Herve (E Belgium)

A landslide susceptibility map is proposed for the Pays de Herve (E Belgium), where large landslides affect Cretaceous clay outcrop areas. Based on a Bayesian approach, this GIS-supported probabilistic map identifies the areas most susceptible to deep landslides. The database is comprised of the source areas of ten pre-existing landslides (i.e. a sample of 154 grid cells) and of six environmental data layers, namely lithology, proximity to active faults, slope angle and aspect, elevation and distance to the nearest valley-floor. A 30-m-resolution DEM from the Belgian National Geographical Institute is used for the analysis. Owing to the small size of the sample, a special cross-validation procedure of the susceptibility map is performed, which uses in an iterative way each of the landslides to test the predictive power of the map derived from the other landslides. Four different sets of variables are used to produce four susceptibility maps, whose prediction curves are compared. While the prediction rates associated with the models not involving the “proximity to active fault” criterion are comparable to those of the models considering this variable, strong weaknesses inherent in the fault data on which the latter rely suggest that the final susceptibility map should be based on a model that excludes any reference to fault. This highlights the difference between a triggering factor and determining factors, and in the same time broadens the scope of the produced map. A single reactivated slide is also used to test the possibility of predicting future reactivation of existing landslides in the area. Finally, the need for geomorphological control over the mathematical treatment is underlined in order to obtain realistic prediction maps.     

Reconciling pyroclastic flow and surge: the multiphase physics of pyroclastic density currents

Two end-member types of pyroclastic density current are commonly recognized: pyroclastic surges are dilute currents in which particles are carried in turbulent suspension and pyroclastic flows are highly concentrated flows. We provide scaling relations that unify these end-members and derive a segregation mechanism into basal concentrated flow and overriding dilute cloud based on the Stokes number (S_T), the stability factor (Σ_T) and the dense-dilute condition (D_D). We recognize five types of particle behaviors within a fluid eddy as a function of S_T and Σ_T: (1) particles sediment from the eddy, (2) particles are preferentially settled out during the downward motion of the eddy, but can be carried during its upward motion, (3) particles concentrate on the periphery of the eddy, (4) particles settling can be delayed or ‘fast-tracked’ as a function of the eddy spatial distribution, and (5) particles remain homogeneously distributed within the eddy. We extend these concepts to a fully turbulent flow by using a prototype of kinetic energy distribution within a full eddy spectrum and demonstrate that the presence of different particle sizes leads to the density stratification of the current. This stratification may favor particle interactions in the basal part of the flow and D_D determines whether the flow is dense or dilute. Using only intrinsic characteristics of the current, our model explains the discontinuous features between pyroclastic flows and surges while conserving the concept of a continuous spectrum of density currents.     

Data assimilation (4D-VAR) to forecast flood in shallow-waters with sediment erosion

In this paper, the four-dimensional variational data assimilation technique (4D-VAR) is presented as a tool to forecast floods. Our study is limited to purely hydrological flows and supposes that the weather, here a big rain, has been already forecasted by meteorological services. The technique consists in minimizing, in the sense of Lagrange, the cost function: a measure of the difference between calculated data and available observations, here the water level. This is done under constraints that are the equations of the physical model. In our case, we modified the shallow-water equations to include a simplified sediment transport model. The steepest descent algorithm is then used to find the minimum. This is made possible because we can compute analytically the gradient of the cost function by using the adjoint equations of the model. As an application of the 4D-VAR technique, the overflowing of the Chicoutimi River at the Chute-Garneau dam, during the 1996 flood, is investigated. It is found that the 4D-VAR method reduces the error in the water height forecast even when the erosion model is not activated. In terms of Lyapunov exponents, we estimate the predictability horizon of such an event to be about half-an-hour after a big rain. However, this limit of predictability can be increased by using more observations or by using a finer computational grid.     

Changes in stream water quality due to logging of the boreal forest in the Montmorency Forest,Québec

Summer stream water quality was monitored before and following the logging of 50% of the boreal forest within three small watersheds (<50 ha) nested in the ‘Ruisseau des Eaux‐Volées’ Experimental Watershed, Montmorency Forest (Québec, Canada). Logging was conducted in winter, on snow cover according to recommended best management practices (BMPs) to minimize soil disturbance and protect advance growth. A 20‐m forest buffer was maintained along perennial streams. In watershed 7·2, cut‐blocks were located near the stream network and logging was partially allowed within the riparian buffer zone. In watersheds 7·5 and 7·7, logging occurred farther away from the stream network. Observations were also made for watershed 7·3 that collected the runoff from watersheds 7·2 and 7·5, and watershed 7·6, the uproad portion of watershed 7·7. The control watershed 0·2 was contiguous to the impacted watersheds and remained undisturbed. Following clearcutting, changes in summer daily maximum and minimum stream temperatures remained within ± 1 °C while changes in diurnal variation did not decrease by more than 0·5 °C. Concentrations of NO₃^? greatly increased by up to 6000% and concentrations of K⁺ increased by up to 300% during the second summer after logging. Smaller increases were observed for Fe_total (up to 71%), specific conductance (up to 26%), and Mg²⁺ (up to 19%). Post‐logging pH decreased slightly by no more than 7% while PO₄^3? concentration remained relatively constant. Suspended sediment concentrations appeared to increase during post‐logging, but there was not enough pre‐logging data to statistically confirm this result. Logging of moderate intensity and respecting established BMPs may account for the limited changes of water quality parameters and the low exceedances of the criteria for the protection of aquatic life. The proximity of the cutover to the stream network and logging within the riparian zone did not appear to affect water quality. Copyright © 2008 John Wiley & Sons, Ltd.     

Wave groupiness and spectral bandwidth as relevant parameters for the performance assessment of wave energy converters

To date the estimation of long-term wave energy production at a given deployment site has commonly been limited to a consideration of the significant wave height H_s and mean energy period T_e. This paper addresses the sensitivity of power production from wave energy converters to the wave groupiness and spectral bandwidth of sea states. Linear and non-linear systems are implemented to simulate the response of converters equipped with realistic power take-off devices in real sea states. It is shown in particular that, when the converters are not much sensitive to wave directionality, the bandwidth characteristic is appropriate to complete the set of overall wave parameters describing the sea state for the purpose of estimating wave energy production.     

Integrating isolated and riparian wetland modules in the PHYSITEL/HYDROTEL modelling platform: model performance and diagnosis

Mathematical modelling is a well‐accepted framework to evaluate the effects of wetlands on stream flow and watershed hydrology in general. Although the integration of wetland modules into a distributed hydrological model represents a cost‐effective way to make this assessment, the added value brought by landscape‐specific modules to a model's ability to replicate basic hydrograph characteristics remains unclear. The objectives of this paper were the following: (i) to present the adaptation of PHYSITEL (a geographic information system) to parameterize isolated and riparian wetlands; (ii) to describe the integration of specific isolated wetland and riparian wetland modules into HYDROTEL, a distributed hydrological model; and (iii) to evaluate the performance of the updated modelling platform with respect to the capacity of replicating various hydrograph characteristics. To achieve this, two sets of simulations were performed (with and without wetland modules), and the added value was assessed at three river segments of the Becancour River watershed, Quebec, Canada, using six general goodness‐of‐fit indicators and 14 water flow criteria. A sensitivity analysis of the wetland module parameters was performed to characterize their impact on stream flows of the modelled watershed. Results of this study indicate the following: (i) integration of specific wetland modules can slightly increase the capacity of HYDROTEL to replicate basic hydrograph characteristics; and (ii) the updated modelling platform allows for the explicit assessment of the impact of wetlands (e.g. typology and location) on watershed hydrology. Copyright © 2015 John Wiley & Sons, Ltd.     

Cratering and modification of wet‐target craters: Projectile impact experiments and field observations of the Lockne marine‐target crater (Sweden)

Abstract— Marine impacts are one category of crater formation in volatile targets. At target water depths exceeding the diameter of the impactor, the zones of vaporization, melting, and excavation of the standard land‐target cratering model develop partially or entirely in the water column. The part of the crater that has a potential of being preserved (seafloor crater) may to a great extent be formed by material emplacement and excavation processes that are very different from land‐target craters. These processes include a high‐energy, water‐jet‐driven excavation flow. At greater water depths, the difference in strength of the target layers causes a concentric crater to evolve. The crater consists of a wide water cavity with a shallow excavation flow along the seabed surrounding a nested, deeper crater in the basement. The modification of the crater is likewise influenced by the water through its forceful resurge to fill the cavity in the water mass and the seafloor. The resurge flow is strongly erosive and incorporates both ejecta and rip‐up material from the seabed surrounding the excavated crater. A combination of field observations and impact experiments has helped us analyze the processes affecting the zone between the basement crater and the maximum extent of the water cavity. The resurge erosion is facilitated by fragmentation of the upper parts of the solid target caused by a) spallation and b) vibrations from the shallow excavation flow and, subsequently, c) the vertical collapse of the water cavity rim wall. In addition, poorly consolidated and saturated sediments may collapse extensively, possibly aided by a violent expansion of the pore water volume when it turns into a spray during passage of the rarefaction wave. This process may also occur at impacts into water‐saturated targets without an upper layer of seawater present. Our results have implications for impacts on both Earth and Mars, and possibly anywhere in the solar system where volatiles exist/have existed in the upper part of the target.     

Estimating River Conductance from Prior Information to Improve Surface‐Subsurface Model Calibration

Most groundwater models simulate stream‐aquifer interactions with a head‐dependent flux boundary condition based on a river conductance (CRIV). CRIV is usually calibrated with other parameters by history matching. However, the inverse problem of groundwater models is often ill‐posed and individual model parameters are likely to be poorly constrained. Ill‐posedness can be addressed by Tikhonov regularization with prior knowledge on parameter values. The difficulty with a lumped parameter like CRIV, which cannot be measured in the field, is to find suitable initial and regularization values. Several formulations have been proposed for the estimation of CRIV from physical parameters. However, these methods are either too simple to provide a reliable estimate of CRIV, or too complex to be easily implemented by groundwater modelers. This paper addresses the issue with a flexible and operational tool based on a 2D numerical model in a local vertical cross section, where the river conductance is computed from selected geometric and hydrodynamic parameters. Contrary to other approaches, the grid size of the regional model and the anisotropy of the aquifer hydraulic conductivity are also taken into account. A global sensitivity analysis indicates the strong sensitivity of CRIV to these parameters. This enhancement for the prior estimation of CRIV is a step forward for the calibration and uncertainty analysis of surface‐subsurface models. It is especially useful for modeling objectives that require CRIV to be well known such as conjunctive surface water‐groundwater use.