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.     

Reconstruction of Horizontal Plasma Motions at the Photosphere from Intensitygrams: A Comparison Between DeepVel,LCT, FLCT,and CST

Direct measurements of plasma motions in the photosphere are limited to the line-of-sight component of the velocity. Several algorithms have therefore been developed to reconstruct the transverse components from observed continuum images or magnetograms. We compare the space and time averages of horizontal velocity fields in the photosphere inferred from pairs of consecutive intensitygrams by the LCT, FLCT, and CST methods and the DeepVel neural network in order to identify the method that is best suited for generating synthetic observations to be used for data assimilation. The Stein and Nordlund (Astrophys. J. Lett.753, L13, 2012) magnetoconvection simulation is used to generate synthetic SDO/HMI intensitygrams and reference flows to train DeepVel. Inferred velocity fields show that DeepVel performs best at subgranular and granular scales and is second only to FLCT at mesogranular and supergranular scales.     

Numerical analysis of frost effects in porous media. Benefits and limits of the finite element poroelasticity formulation

The aim of this paper is to analyse the performance of a finite element formulation usable for predicting the mechanical consequences of frost effects on porous media. It considers the characteristics of porous media and how the frost action can be assessed. The problem is then separated into two parts: thermal and poromechanical calculations. The constitutive equations developed in the framework of poromechanics are presented and the implementation in a usual finite element poroelasticity formulation based on Zuber's method is adopted. An analysis of the time‐step influence on the convergence rate is given and leads us to propose a simple method in order to obtain objectivity of the finite element response and avoid over‐long calculations. Frost effect simulations are carried out on real porous media (two fired clays) as a case study. Although the experimental behaviour of the porous media subjected to frost action is in accordance with some observations, the calculated strains appear to be overestimated compared with measurements. The problem could be largely attributable to the difficulty of assessing permeability evolution during frost development. Copyright © 2011 John Wiley & Sons, Ltd.     

Finite element analysis of the effect of corrosion on the behavior of reinforced earth walls

A simple model for the corrosion‐induced loss of stiffness and strength of the steel strips of earth‐reinforced walls was introduced in a finite element simulation of the long‐term behavior of the wall, in which the backfill‐strips interactions are taken into account by means of a generalized homogenization procedure (called a multiphase model). The results show an initial phase of slow displacements induced by the loss of stiffness, followed after a few decades by a steep acceleration of the displacements, leading to wall failure. The influences of the parameter controlling corrosion, the backfill cohesion and the heterogeneity of the corrosion process are discussed. Results are used to discuss a strategy for reinforced earth wall surveillance. Copyright © 2011 John Wiley & Sons, Ltd.     

Evaluation of antifouling properties of non-toxic marine paints

The most common way to prevent biofouling of immersed surfaces is to coat the exposed surfaces with a paint containing a biocide acting on marine organisms by diffusion into the surrounding water. However diffusion of toxins into the marine environment causes pollution and serious economic problems. This paper outlines the anti-microfouling properties of some insolubilized quaternary ammonium salts (grafted onto a vinyl copolymer by means of a covalent non-hydrolyzable bond). After immersion times in sea-water up to 4 months, the microbial cover was limited to bacterial forms, without microalgae or cyanobacteria as observed on the untreated, or tin-salt painted surfaces. Bacteria are mostly on unicell forms and numerous areas appear free of any microorganisms. As microbiofouling is an obligate step to macrofouling, this unleachable antifouling treatment could result in a prolonged protection.     

Dynamic faulting under rate-dependent friction

We discuss the effects of rate-dependent friction on the propagation of seismic rupture on active faults. Several physicists using Burridge and Knopoff's box and spring model of faulting have proposed that fault complexity may arise from the spontaneous development of a self-similar stress distribution on the fault plane. If this model proves to be correct, it has important consequences for the origin of the complexity of seismic sources. In order to test these ideas on a more realistic earthquake model, we developed a new boundary integral equation method for studying rupture propagation along an antiplane fault in the presence of nonlinear rate-dependent friction. We study rupture dynamics of models with single and twin asperities. In our models, asperities are places on the fault with a higher value of prestress. Othewise all fault parameters are homogeneous. We show that for models with such asperities, a slip velocity weakening friction leads to the propagation of supersonic healing phases and to the spontaneous arrest of fracture if the prestress outside the asperities is low enough. For models with asperities, we can also observe narrow slip velocity pulses, qualitatively similar to the so-called Heaton pulses observed in some earthquake accelerograms. We also observe a complex distribution of stress after the rupture that depends on details of the initial distribution of asperities and on the details of the friction law.     

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.     

Twin structure in coesite studied by high resolution electron microscopy

A synthetic crystal of coesite, SiO₂, has been found to contain microtwins. The (100) reflection twin has been made evident on incoherent (010) planes. The twin structure is directly imaged by high resolution electron microscopy at 2,5 Å resolution, allowing an atomic model to be proposed. This twin could be incoherent on large surfaces and the dislocation density is predicted. The X-ray crystallographic data are analyzed critically in view of the presence of such defects in single crystals. The mutual orientation between a silicon matrix and coesite precipitates in Czochralski-silicon annealed specimens is proposed to be non unique as a consequence of the easy twinning of coesite.     

Detailed gravity study of the offshore structure of Piton de la Fournaise volcano,Reunion Island

We present an interpretation of gravity data acquired in 1984 by the French R/V Jean Charcot on the submarine part of the eastern flank of Piton de la Fournaise volcano. We comment on the Bouguer anomaly map and give a quantitative interpretation of three gravity profiles. The main results are that a gravity high over Grand Brûlé, the lower subaerial part of the eastern flank, does not extend far offshore and that an anomalous topographic feature, discovered in 1982 on the submarine eastern flank, is characterized by a large negative anomaly. We propose three hypotheses to explain the origin of this anomaly, i. e., it marks the site of a new volcano, or it is a consequence of lateral volcanism from a volcano older than Piton de la Fournaise, or more probably, it represents a great landslide deposit.