Accuracy and efficiency of time integration methods for 1D diffusive wave equation     

Sylvain Weill  Raphael di Chiara-Roupert  Philippe Ackerer 《Computational Geosciences》2014,18(5):697-709

The diffusive wave approximation of the Saint-Venant equations is commonly used in hydrological models to describe surface flow processes. Numerous numerical approaches can be used to solve this highly nonlinear equation. Nonlinear time integration schemes—also called methods of lines (MOL)—were proven very efficient to solve other nonlinear problems in geosciences but were never considered to deal with surface flow modeling with the diffusive wave equation. In this paper, we study the relative performance of different time and space integration schemes by comparing the results obtained with classical approaches and with nonlinear time integration approaches. The results show that (i) the integration method with a higher order in space shows high accuracy regarding an integrated indicator such as the global mass balance error but is less accurate regarding local indicators, and (ii) nonlinear time integration techniques perform better than classical ones. Overall, it seems that integration techniques combining nonlinear time integration and a low spatial order need to be considered when developing hydrological modeling tools owing to their simplicity of implementation and very good performance.  相似文献   

Machine Learning Feature Selection Methods for Landslide Susceptibility Mapping     

Natan Micheletti  Loris Foresti  Sylvain Robert  Michael Leuenberger  Andrea Pedrazzini  Michel Jaboyedoff  Mikhail Kanevski 《Mathematical Geosciences》2014,46(1):33-57

This paper explores the use of adaptive support vector machines, random forests and AdaBoost for landslide susceptibility mapping in three separated regions of Canton Vaud, Switzerland, based on a set of geological, hydrological and morphological features. The feature selection properties of the three algorithms are studied to analyze the relevance of features in controlling the spatial distribution of landslides. The elimination of irrelevant features gives simpler, lower dimensional models while keeping the classification performance high. An object-based sampling procedure is considered to reduce the spatial autocorrelation of data and to estimate more reliably generalization skills when applying the model to predict the occurrence of new unknown landslides. The accuracy of the models, the relevance of features and the quality of landslide susceptibility maps were found to be high in the regions characterized by shallow landslides and low in the ones with deep-seated landslides. Despite providing similar skill, random forests and AdaBoost were found to be more efficient in performing feature selection than adaptive support vector machines. The results of this study reveal the strengths of the classification algorithms, but evidence: (1) the need for relying on more than one method for the identification of relevant variables; (2) the weakness of the adaptive scaling algorithm when used with landslide data; and (3) the lack of additional features which characterize the spatial distribution of deep-seated landslides.  相似文献   

Geology and activity around volcanoes on Io from the analysis of NIMS spectral images     

Sylvain Douté  Rosaly Lopes  Robert Carlson  the Galileo NIMS Team 《Icarus》2004,169(1):175-196

The last two successful flybys of Io by Galileo in 2001 (orbits I31, I32) allowed the Near Infrared Mapping Spectrometer to enrich its collection of IR spectral image cubes of the satellite. These data cover hemispheric portions of Io, several volcanic centers as well as their surroundings with a spatial resolution ranging from 2 to 93 km pixel⁻¹. They map thermal emission from the hot-spots and the distribution of solid SO₂ in the 1.0-4.7 μm spectral range. We obtain maps of SO₂ abundance and granularity from the NIMS data using the method of Douté et al. (2002, Icarus 158, 460-482). The maps are correlated to distinguish four different physical units that indicate zones of SO₂ condensation, metamorphism and sublimation. We relate these information with visible images from Galileo's Solid State Imaging System and with detailed mapping of the thermal emission produced by Io's surface. Our principal goal is to understand the mechanisms controlling how lava, pyroclastics and gas are emitted by different types of volcanoes and how these products evolve. The 800 km diameter white ring of fallout created by a violent “Pillanian” eruption during summer of 2001 is at least partly composed of solid SO₂ and has enriched preexisting regional deposits. Orange materials have been recently or are currently emplaced 240 km south from the main eruption site, possibly as sulfur flows. A similar event may have taken place in the past at Ababinili Patera (12.5° N, 142° W). Carefull study of SO₂ maps covering the Emakong region also suggests that sulfur forms the bright channel-fed flow emerging from the south eastern side of the caldera. Within the main caldera of Tvashtar Catena completely cooled patches of crust exist. Elsewhere, the caldera is still cooling from previous episodes of flooding. We confirm that Amirani emits constantly large amount of SO₂ gas by interaction of fresh lava with the volatiles of the underlying plains. Nevertheless SO₂ frost is not the major component of the bright white ring seen in the SSI images. Over the whole Gish Bar region, SO₂ frost seems barely stable and is constantly regenerated. The stability increases along gray filamentary structures which could be faults filled with materials having peculiar thermal properties. Northwest of Gish Bar Patera, a localized bright deposit shows an unusual spectral signature potentially indicative of H₂O molecules forming ice crystals or being trapped in a nonidentified matrix. The Chaac region may present a thickened old crust reducing the geothermal flux to levels lower than 0.5 W m⁻² and thus creating a cold trap for SO₂. Looking at the abundance and degree of metamorphose of SO₂, we establish the relative age of different flows and ejecta for the Sobo Fluctus. Finally the assumption that the white patches in visible images indicate SO₂ rich deposits is once again challenged. In the Camaxtli region we identify a topographically controlled compact white deposit showing only moderate SO₂ abundance. In contrast, we detect two spots of quite pure SO₂ ice on the gray flanks of Emakong. Furthermore, the close association of fumarolic SO₂ and red S₂ already noted for several volcanic centers is observed at Tupan.  相似文献   

Impact of Terrain Heterogeneity on Coherent Structure Properties: Numerical Approach   总被引：2，自引：2，他引：0  

Clement Fesquet  Sylvain Dupont  Philippe Drobinski  Thomas Dubos  Christian Barthlott 《Boundary-Layer Meteorology》2009,133(1):71-92

A three-dimensional large-eddy simulation (LES) model, which includes the effects of plant–atmosphere interactions, is used to study the effects of surface inhomogeneities on near-surface coherent structures over an open field and behind a forest canopy. These simulated conditions are representative of two wind sectors of the Site Instrumental de Recherche par Télédétection Atmosphérique (SIRTA) experimental site at the Institut Pierre Simon Laplace, Palaiseau, France. Coherent structure properties deduced from wavelet transforms of the simulated near-surface vertical velocity time series are not modified by upstream terrain heterogeneities, in agreement with site measurements. This feature is related to the nature of structures detected from the vertical velocity time series. The turbulence close to the surface seems composed of both local coherent structures and large coherent structures reflecting outer-layer properties, which depend on the overall surface heterogeneity or upstream heterogeneity. It is argued that the streamwise velocity is representative of these large outer-layer structures that impinge onto the ground through a top-down mechanism as identified through the space–time correlation of the wind velocity components. In contrast, the vertical velocity is more representative of small structures resulting from the impingement of the large outer-layer structures. These small structures represent locally-generated, active turbulence, which adjusts rapidly to local surface conditions, and consequently they are only weakly dependent on upstream heterogeneities.  相似文献   

Conflicting structural and geochronological data from the Ibituruna quartz-syenite (SE Brazil): Effect of protracted “hot” orogeny and slow cooling rate?     

Sylvain Petitgirard  Alain Vauchez  Marcos Egydio-Silva  Olivier Bruguier  Pierre Camps  Patrick Moni  Marly Babinski  Mathieu Mondou 《Tectonophysics》2009,477(3-4):174

The Ibituruna quartz-syenite was emplaced as a sill in the Ribeira-Araçuaí Neoproterozoic belt (Southeastern Brazil) during the last stages of the Gondwana supercontinent amalgamation. We have measured the Anisotropy of Magnetic Susceptibility (AMS) in samples from the Ibituruna sill to unravel its magnetic fabric that is regarded as a proxy for its magmatic fabric. A large magnetic anisotropy, dominantly due to magnetite, and a consistent magnetic fabric have been determined over the entire Ibituruna massif. The magmatic foliation and lineation are strikingly parallel to the solid-state mylonitic foliation and lineation measured in the country-rock. Altogether, these observations suggest that the Ibituruna sill was emplaced during the high temperature (~ 750 °C) regional deformation and was deformed before full solidification coherently with its country-rock. Unexpectedly, geochronological data suggest a rather different conclusion. LA-ICP-MS and SHRIMP ages of zircons from the Ibituruna quartz-syenite are in the range 530–535 Ma and LA-ICP-MS ages of zircons and monazites from synkinematic leucocratic veins in the country-rocks suggest a crystallization at ~ 570–580 Ma, i.e., an HT deformation > 35My older than the emplacement of the Ibituruna quartz-syenite. Conclusions from the structural and the geochronological studies are therefore conflicting. A possible explanation arises from ⁴⁰Ar–³⁹Ar thermochronology. We have dated amphiboles from the quartz-syenite, and amphiboles and biotites from the country-rock. Together with the ages of monazites and zircons in the country-rock, ⁴⁰Ar–³⁹Ar mineral ages suggest a very low cooling rate: < 3 °C/My between 570 and ~ 500 Ma and ~ 5 °C/My between 500 and 460 Ma. Assuming a protracted regional deformation consistent over tens of My, under such stable thermal conditions the fabric and microstructure of deformed rocks may remain almost unchanged even if they underwent and recorded strain pulses separated by long periods of time. This may be a characteristic of slow cooling “hot orogens” that rocks deformed at significantly different periods during the orogeny, but under roughly unchanged temperature conditions, may display almost indiscernible microstructure and fabric.  相似文献   

Tectonometamorphic evolution of the Atbashi high‐P units (Kyrgyz CAOB,Tien Shan): Implications for the closure of the Turkestan Ocean and continental subduction–exhumation of the South Kazakh continental margin          下载免费PDF全文

Chloé Loury  Yann Rolland  Stéphane Guillot  Pierre Lanari  Clément Ganino  Raphael Melis  Anthony Jourdon  Carole Petit  Olivier Beyssac  Sylvain Gallet  Patrick Monié 《Journal of Metamorphic Geology》2018,36(8):959-985

The South Tien Shan (STS) belt results from the last collision event in the western Central Asian Orogenic Belt (CAOB). Understanding its formation is of prime importance in the general framework of the CAOB. The Atbashi Range preserves high‐P (HP) rocks along the STS suture, but still, its global metamorphic evolution remains poorly constrained. Several HP units have been identified: (a) a HP tectonic mélange including boudins of mafic eclogites in a sedimentary matrix, (b) a large (>100 km long) high‐P metasedimentary unit (HPMU) and (c) a lower blueschist facies accretionary prism. Raman Spectroscopy on carbonaceous material combined with phengite and chlorite multiequilibria and isochemical phase diagram modelling indicates that the HPMU recorded homogeneous P–T conditions of 23–25 kbar and 560–570°C along the whole unit. ⁴⁰Ar/³⁹Ar dating on phengite from the HPMU ranges between 328 and 319 Ma at regional scale. These ages are interpreted as (re‐) crystallization ages of phengite during T_max conditions at a pressure range of 20–25 kbar. Thermobarometry on samples from the HP tectonic mélange provides similar metamorphic peak conditions. Thermobarometry on the blueschist to lower greenschist facies accretionary prism indicates that it underwent P–T conditions of 5–6 kbar and 290–340°C, highlighting a 17–20 kbar pressure gap between the HPMU‐tectonic mélange units and the accretionary prism. Comparison with available geochronological data suggests a very short time span between the prograde path (340 Ma), HP metamorphic peak (330 Ma), the T_max (328–319 Ma) and the final exhumation of the HPMU (303–295 Ma). Extrusion of the HPMU, accommodated by a basal thrust and an upper detachment, was driven by buoyant forces from 70–75 km up to 60 km depth, which directly followed continental subduction and detachment of the HPMU. At crustal depths, extrusion was controlled by collisional tectonics up to shallow levels. Lithological homogeneity of the HPMU and its continental‐derived character from the North Tien Shan suggest this unit corresponds to the hyper‐extended continental margin of the Kazakh continent, subducted southward below the north continental active margin of the Tarim craton. Integration of the available geological data allows us to propose a general geodynamic scenario for Tien Shan during the Carboniferous with a combination of (a) N‐dipping subduction below the Kazakh margin of Middle Tien Shan until 390–340 Ma and (b) S‐dipping subduction of remaining Turkestan marginal basins between 340 and 320 Ma.  相似文献   

Short-term impact of reservoir cleaning on the microhabitat use of three non-salmonid fishes in a piedmont river in south west France     

Sylvain Mastrorillo  Francis Dauba 《Aquatic Sciences - Research Across Boundaries》1999,61(4):323-336

As with many rivers of the Pyrenees Mountains, reservoirs on the River Ariège (SW Francc) are used for hydroelectricity generation, irrigation and to maintain summer flows in the River Garonne. Two river reservoirs on the River Ariège were simultaneously emptied in April 1991 to remove accumulated sediment. Fish and ten environmental variables were sampled using Point Abundance Sampling by electrofishing before (February 1991) and one year later (February 1992) at three study sites downstream of the dams. We examined the short-term impact of reservoir cleaning on the relative abundance, habitat availability and habitat use of three species of small fish: minnow Phoxinus phoxinus (L.), gudgeon Gobio gobio (L.) and stone loach Barbatula barbatula (L.).¶After the reservoir cleaning, areas covered by gravel, sand and mud increased significantly as the result of the sediment release. The three species did not react in the same manner to these environmental changes. Stone loach appeared to be the most sensitive species, with its relative abundance decreasing at all sites. Conversely, minnow seemed to be the most resistant species and as such is a leading pioneer species for colonizing new or impacted stretches of rivers channels. Finally, the number of gudgeon decreased particularly at the most downstream site where large deposits of mud were recorded.  相似文献   

A low-dimensional subsurface model for saturated and unsaturated flow processes: ability to address heterogeneity     

Sylvain Weill  Frederick Delay  Yi Pan  Philippe Ackerer 《Computational Geosciences》2017,21(2):301-314

A low-dimensional model that describes both saturated and unsaturated flow processes in a single equation is presented. Subsurface flow processes in the groundwater, the vadose zone, and the capillary fringe are accounted for through the computation of aggregated hydrodynamic parameters that result from the integration of the governing flow equations from the bedrock to the land surface. The three-dimensional subsurface flow dynamics are thus described by a two-dimensional equation, allowing for a drastic reduction of model unknowns and simplification of the model parameterizations. This approach is compared with a full resolution of the Richards equation in different synthetic test cases. Because the model reduction stems from the vertical integration of the flow equations, the test cases all use different configurations of heterogeneity for vertical cross-sections of a soil-aquifer system. The low-dimensional flow model shows strong consistency with results from a complete resolution of the Richards equation for both the water table and fluxes. The proposed approach is therefore well suited to the accurate reproduction of complex subsurface flow processes.  相似文献   

Sediment yield of the lower Tana River,Kenya, is insensitive to dam construction: sediment mobilization processes in a semi‐arid tropical river system          下载免费PDF全文

Naomi Geeraert  Fred Ochieng Omengo  Fredrick Tamooh  Paolo Paron  Steven Bouillon  Gerard Govers 《地球表面变化过程与地形》2015,40(13):1827-1838

Dam construction in the 1960s to 1980s significantly modified sediment supply from the Kenyan uplands to the lower Tana River. To assess the effect on suspended sediment fluxes of the Tana River, we monitored the sediment load at high temporal resolution for 1 year and complemented our data with historical information. The relationship between sediment concentration and water discharge was complex: at the onset of the wet season, discharge peaks resulted in high sediment concentrations and counterclockwise hysteresis, while towards the end of the wet season, a sediment exhaustion effect led to low concentrations despite the high discharge. The total sediment flux at Garissa (c. 250 km downstream of the lowermost dam) between June 2012 and June 2013 was 8.8 Mt yr^‐1. Comparison of current with historical fluxes indicated that dam construction had not greatly affected the annual sediment flux. We suggest that autogenic processes, namely river bed dynamics and bank erosion, mobilized large quantities of sediments stored in the alluvial plain downstream of the dams. Observations supporting the importance of autogenic processes included the absence of measurable activities of the fall‐out radionuclides ⁷Be and ¹³⁷Cs in the suspended sediment, the rapid lateral migration of the river course, and the seasonal changes in river cross‐section. Given the large stock of sediment in the alluvial valley of the Tana River, it may take centuries before the effect of damming shows up as a quantitative reduction in the sediment flux at Garissa. Many models relate the sediment load of rivers to catchment characteristics, thereby implicitly assuming that alterations in the catchment induce changes in the sediment load. Our research confirms that the response of an alluvial river to external disturbances such as land use or climate change is often indirect or non‐existent as autogenic processes overwhelm the changes in the input signal. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

Compensatory mitigation in marine ecosystems: Which indicators for assessing the “no net loss” goal of ecosystem services and ecological functions?   总被引：1，自引：0，他引：1  

Harold Levrel  Sylvain Pioch  Richard Spieler 《Marine Policy》2012

Recent years have seen increasing interest in the concepts of compensation and ecosystem services. Regulation systems in the United States dealing with environmental protection (Superfund Act, Oil Pollution Act, National Environment Policy Act, Clean Water Act, Endangered Species Act, etc.) require those responsible for damage to ecosystem services to compensate for it “physically” and restore these services for the benefit of the entire population. This article, using simple indicators of compensation identified in the literature, attempts to analyze what types of ecological compensation are adopted, how performance is assessed, how standards on ecological equivalencies are adopted, and what are the costs of this compensation. To perform this analysis, compensatory measures carried out during the last ten years in the case of coastal and marine ecosystems in Florida have been addressed. The results show that: analysis criteria for the equivalencies between ecosystem services lost due to damage and ecosystem services gained due to compensatory measures are questionable; most compensation monitoring is for a relatively brief period of time and the data obtained during this period may be insufficient for assessing the net effect of the compensatory measure; the weaknesses regarding criteria for the equivalencies and the uncertainty about the relevant time-scale can be counter-balanced by increasing the area of compensation, a problematic solution at best.  相似文献