Éditeur(s) : HAL CCSD American Association for the Advancement of Science
Résumé :  
Amazon forests are a key but poorly understood component of the global carbon cycle. If, as anticipated, they dry this century, they might accelerate climate change through carbon losses and changed surface energy balances. We used records from multiple long-term monitoring plots across Amazonia to assess forest responses to the intense 2005 drought, a possible analog of future events. Affected forest lost biomass, reversing a large long-term carbon sink, with the greatest impacts observed where the dry season was unusually intense. Relative to pre-2005 conditions, forest subjected to a 100-millimeter increase in water deficit lost 5.3 megagrams of aboveground biomass of carbon per hectare. The drought had a total biomass carbon impact of 1.2 to 1.6 petagrams (1.2 × 1015 to 1.6 × 1015 grams). Amazon forests therefore appear vulnerable to increasing moisture stress, with the potential for large carbon losses to exert feedback on climate change.
ISSN: 0036-8075

hal-01032111
https://hal.archives-ouvertes.fr/hal-01032111
DOI : 10.1126/science.1164033

Éditeur(s) : HAL CCSD Nature Publishing Group
Résumé : While Amazonian forests are extraordinarily diverse, the abundance of trees is skewed strongly towards relatively few 'hyperdominant' species. In addition to their diversity, Amazonian trees are a key component of the global carbon cycle, assimilating and storing more carbon than any other ecosystem on Earth. Here we ask, using a unique data set of 530 forest plots, if the functions of storing and producing woody carbon are concentrated in a small number of tree species, whether the most abundant species also dominate carbon cycling, and whether dominant species are characterized by specific functional traits. We find that dominance of forest function is even more concentrated in a few species than is dominance of tree abundance, with only ≈1% of Amazon tree species responsible for 50% of carbon storage and productivity. Although those species that contribute most to biomass and productivity are often abundant, species maximum size is also influential, while the identity and ranking of dominant species varies by function and by region.
ISSN: 2041-1723

hal-01204225
https://hal.archives-ouvertes.fr/hal-01204225
DOI : 10.1038/ncomms7857
PRODINRA : 305614
PUBMED : 25919449

Éditeur(s) : HAL CCSD Wiley
Résumé : International audience
Aim To examine variation in the phylogenetic diversity (PD) of tree communities across geographical and environmental gradients in Amazonia.Location Two hundred and eighty-three c.1ha forest inventory plots from across Amazonia.Methods We evaluated PD as the total phylogenetic branch length across species in each plot (PDss), the mean pairwise phylogenetic distance between species (MPD), the mean nearest taxon distance (MNTD) and their equivalents standardized for species richness (ses.PDss, ses.MPD, ses.MNTD). We compared PD of tree communities growing (1) on substrates of varying geological age; and (2) in environments with varying ecophysiological barriers to growth and survival.Results PDss is strongly positively correlated with species richness (SR), whereas MNTD has a negative correlation. Communities on geologically young- and intermediate-aged substrates (western and central Amazonia respectively) have the highest SR, and therefore the highest PDss and the lowest MNTD. We find that the youngest and oldest substrates (the latter on the Brazilian and Guiana Shields) have the highest ses.PDss and ses.MNTD. MPD and ses.MPD are strongly correlated with how evenly taxa are distributed among the three principal angiosperm clades and are both highest in western Amazonia. Meanwhile, seasonally dry tropical forest (SDTF) and forests on white sands have low PD, as evaluated by any metric.Main conclusions High ses.PDss and ses.MNTD reflect greater lineage diversity in communities. We suggest that high ses.PDss and ses.MNTD in western Amazonia results from its favourable, easy-to-colonize environment, whereas high values in the Brazilian and Guianan Shields may be due to accumulation of lineages over a longer period of time. White-sand forests and SDTF are dominated by close relatives from fewer lineages, perhaps reflecting ecophysiological barriers that are difficult to surmount evolutionarily. Because MPD and ses.MPD do not reflect lineage diversity per se, we suggest that PDss, ses.PDss and ses.MNTD may be the most useful diversity metrics for setting large-scale conservation priorities.
ISSN: 1366-9516

hal-01557462
https://hal.univ-lorraine.fr/hal-01557462
DOI : 10.1111/ddi.12357

Éditeur(s) : HAL CCSD Springer Verlag (Germany) Springer Verlag
Résumé : Production rates in horizontal shale gas wells display declines that are controlled by the low permeability and the intrinsic heterogeneity of the shale matrix. We present an original multi-continuum approach that yields a physical model able to reproduce the complexity of the decreasing gas rates. The model describes the dynamics of gas rate as function of the physical reservoir parameters and geometry, while the shale matrix heterogeneity is accounted for by a stochastic description of transmissivity field. From the 3D (Dimensional) problem setting, including the heterogeneous shale matrix, the fractures generated by the hydrofracking operations, as well as the production well characteristics, we establish an effective upscaled 1D model for the gas pressures in fracture and matrix as well as the volumetric flux. We analyse the decline curves behaviour, and we identify the time scales that characterize the dynamics of the gas rate decline using explicit analytical Laplace space solutions of the upscaled process model. Asymptotically, the flux curves decrease exponentially, while in an intermediate regime we find a power-law behaviour, in which the flux scales with a power law in time as (t^{-beta }), where (beta ) reflects the medium heterogeneity. We use this solution to fit a set of real data displaying distinctly different decline trends and study the sensitivity of the model to the reservoir parameters in order to identify their respective controls at the different stages of the decline curve dynamics. Results indicate that the initial value of the gas rate is determined by the transmissivity of the fractures and the initial pressure of the gas in the shale matrix. The latter causes mainly a shift in the entire decline curve. The early time of decline curve shape is primarily controlled by the fracture properties (compressibility and transmissivity). During the main part of the economically valuable production times, i.e. before the production rate drops exponentially, the decline curve is strongly controlled by the properties of the shale rocks including their heterogeneity, which is modelled by two parameters describing the non-Fickian pressure diffusion effects in a stochastic framework.
ISSN: 0169-3913

hal-01217247
https://hal.archives-ouvertes.fr/hal-01217247
DOI : 10.1007/s11242-015-0504-y

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Mongolia represents the northernmost area affected by the India–Asia collision, and it is actively deformed along transpressive belts closely associated with large-scale strike-slip faults. The active and past mantle flow beneath this region is, however, poorly known. In order to investigate deep mantle deformation beneath central Mongolia and its relation with the surrounding major structures such as the Siberian craton, the Gobi–Altay belt and the Baikal rift, a NS-trending profile of broadband seismic stations has been deployed in the summer 2003 from the southern Siberian craton to the Gobi–Altay range, crossing the entire Hangay dome. Mantle flow is deduced from the splitting of teleseismic shear waves such as SKS phases. In eastern Mongolia, the permanent station ULN in Ulaanbaatar reveals the presence of two anisotropic layers, the upper one being oriented NE–SW, close to the trend of lithospheric structures and the lower one NW–SE, close to the trend of Eurasia absolute plate motion. Along the NS profile in central Mongolia, seismic anisotropy deduced from SKS splitting reveals a homogeneous NW–SE trending structure, fully consistent with the observations made in the Altay–Sayan in western Mongolia. The observed delay times of 1.5 to more than 2.0 s favor consistent mantle flow over large mantle thicknesses. Since the lithosphere is less than 100 km thick beneath central Mongolia and since the observed fast directions are parallel to the trend of the lithospheric structures but also close to the trend of the absolute plate motion, we propose that both the lithosphere and the asthenosphere may join their anisotropic effects beneath central Mongolia to explain the large delay times. Although GPS vectors represent the instantaneous displacement of the Earth's surface and SKS splitting the time and vertical integration of finite strain at depth, we use the opportunity of the dense geodetic measurements available in this region to discuss the anisotropy pattern in term of present-day deformation. In the Eurasia-fixed reference frame, GPS and SKS both depict a similar trend beneath central Mongolia, suggesting a lithospheric block “escaping” toward the east that could orient olivine a-axes in the upper mantle, within a transpressive tectonic regime. A different behaviour is observed in western Mongolia: the GPS vectors trend NS, close to the regional compression direction, whereas the fast SKS directions trend EW, suggesting a tectonic regime close to a mode of axial shortening, generating the development of an EW-trending foliation at depth. We therefore propose that Mongolia is a place where active and frozen lithospheric deformation may add their effects, together with the sublithospheric flow. In the three sources of anisotropy inferred, a fundamental role is played by the Siberian craton that acted as an undeformable core of the continent through time: the frozen Paleozoic and Mesozoic structures wrap around the craton, building up the fast anisotropic direction pattern; the present-day sublithospheric flow induced by the plate motion is likely deflecting around its deep roots; finally, the present-day tectonic regime appears to be controlled by the presence of the craton to the north.
ISSN: 0012-821X

hal-00407853
https://hal.archives-ouvertes.fr/hal-00407853
https://hal.archives-ouvertes.fr/hal-00407853/document
https://hal.archives-ouvertes.fr/hal-00407853/file/Barruol_SKS_Mongolia_EPSL2008_hal.pdf
DOI : 10.1016/j.epsl.2008.07.027

Éditeur(s) : HAL CCSD Landes Bioscience
Résumé : International audience
In 2008 we published the first set of guidelines for standardizing research in autophagy. Since then, research on this topic has continued to accelerate, and many new scientists have entered the field. Our knowledge base and relevant new technologies have also been expanding. Accordingly, it is important to update these guidelines for monitoring autophagy in different organisms. Various reviews have described the range of assays that have been used for this purpose. Nevertheless, there continues to be confusion regarding acceptable methods to measure autophagy, especially in multicellular eukaryotes. A key point that needs to be emphasized is that there is a difference between measurements that monitor the numbers or volume of autophagic elements (e.g., autophagosomes or autolysosomes) at any stage of the autophagic process vs. those that measure flux through the autophagy pathway (i.e., the complete process); thus, a block in macroautophagy that results in autophagosome accumulation needs to be differentiated from stimuli that result in increased autophagic activity, defined as increased autophagy induction coupled with increased delivery to, and degradation within, lysosomes (in most higher eukaryotes and some protists such as Dictyostelium) or the vacuole (in plants and fungi). In other words, it is especially important that investigators new to the field understand that the appearance of more autophagosomes does not necessarily equate with more autophagy. In fact, in many cases, autophagosomes accumulate because of a block in trafficking to lysosomes without a concomitant change in autophagosome biogenesis, whereas an increase in autolysosomes may reflect a reduction in degradative activity. Here, we present a set of guidelines for the selection and interpretation of methods for use by investigators who aim to examine macroautophagy and related processes, as well as for reviewers who need to provide realistic and reasonable critiques of papers that are focused on these processes. These guidelines are not meant to be a formulaic set of rules, because the appropriate assays depend in part on the question being asked and the system being used. In addition, we emphasize that no individual assay is guaranteed to be the most appropriate one in every situation, and we strongly recommend the use of multiple assays to monitor autophagy. In these guidelines, we consider these various methods of assessing autophagy and what information can, or cannot, be obtained from them. Finally, by discussing the merits and limits of particular autophagy assays, we hope to encourage technical innovation in the field.
Autophagy

hal-00735751
https://hal.archives-ouvertes.fr/hal-00735751
PUBMED : 22966490

Éditeur(s) : HAL CCSD Elsevier
Résumé : Siberian cratonic mantle
International audience
The processes that control water distribution in nominally anhydrous minerals from peridotites are twofold. Melt depletion will remove water while metasomatism can potentially add water to these minerals. These processes can lead to a wide range of outcomes in water contents, which in turn could play a role in mantle rheology and long-term cratonic root stability. To examine these complexities, water concentrations in minerals from well-characterized peridotites from the Udachnaya kimberlite in the central Siberian craton were analyzed by FTIR. The peridotites span a complete top to bottom cross-section of typical cratonic lithospheric mantle (2-7 GPa and 700-1400 °C). Diffusion modeling of water content profiles across olivine grains shows that water loss during decompression is limited to the 100 μm rims of olivines; the cores preserved their mantle water contents. Water contents range from 6 to 323 ppm wt H2O in olivine, 28-301 ppm H2O in orthopyroxene (opx), 100-272 ppm H2O in clinopyroxene (cpx) and 0-23 ppm H2O in garnet. Melting modeling cannot reproduce the high water contents of cratonic mantle peridotites and any potential partial melting trend must have been erased by later events. The water contents of minerals, however, are correlated with modal abundances of clinopyroxene and garnet, bulk rock FeO, TiO2 and SiO2 as well as with light and middle rare earth elements in clinopyroxene and garnet. These relationships are best interpreted as interaction of residual, melt-depleted peridotites with silicate melt, which produced modal and cryptic metasomatism. Importantly, the water enrichment in the Siberian cratonic mantle took place prior to kimberlite magmatism and eruption. Water addition by metasomatism occurred from pressures >4 GPa all the way to the base of the cratonic root below central Siberia, but was limited to shallower levels (<5 GPa) in the Kaapvaal cratonic lithosphere. The difference in olivine water contents at the deepest levels of the Kaapvaal (<5 ppm H2O) and Siberian (6-323 ppm H2O) cratonic roots may be linked to oxygen fugacity and resulting fluid speciation or, alternatively, to reaction with different metasomatic agents. Calculated viscosities for the deepest Udachnaya samples are similar to those inferred for the asthenosphere. If these xenoliths are representative of the deep cratonic lithosphere, water is not as important a parameter as previously thought in the strength of cratonic lithosphere, otherwise the cratonic root beneath Udachnaya would have been delaminated. Alternatively, the metasomatic xenoliths may not be representative of the Siberian cratonic root and kimberlites preferentially sample cratonic mantle lithosphere material located near, and metasomatized by, melt conduits, which served as channels for upward migration of water-rich melts and fluids including kimberlites. In that case, the cratonic root overall still may have relatively low water contents, which in addition to its less metasomatized (more refractory) and thereby buoyant nature, still play a role in making it strong enough to resist delamination by the surrounding asthenosphere.
ISSN: 0016-7037

hal-01053920
https://hal.archives-ouvertes.fr/hal-01053920
DOI : 10.1016/j.gca.2014.04.011

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
This paper presents a multicontinuum approach to model fractal temporal scaling of catchment response in hydrological systems. The temporal scaling of discharge is quantified in frequency domain by the transfer function HðxÞ, which is defined as the ratio between the spectra of catchment response and recharge time series. The transfer function may scale with frequency x as HðxÞ x2b. While the classical linear and Dupuit models predict exponents of b52 and b51, observations indicate scalings with noninteger exponents b. Such behaviors have been described by multifractal models, which, however, often lack a relation to the medium characteristics. We revisit and extend the classical linear Dupuit aquifer models and discuss their physical meanings in the light of the resulting aquifer dynamics. On the basis of these classical models, we derive a multicontinuum approach that provides physical recharge models which are able to explain fractal behaviors with exponents 1=2 < b < 2. Furthermore, this approach allows to link the fractal dynamics of the discharge process to the physical aquifer characteristics as reflected in the distribution of storage time scales. We systematically analyze the catchment responses in the proposed multicontinuum models, and identify and quantify the time scales which characterize the dynamics of the catchment response to recharge.
ISSN: 0043-1397

insu-00952537
https://hal-insu.archives-ouvertes.fr/insu-00952537
https://hal-insu.archives-ouvertes.fr/insu-00952537/document
https://hal-insu.archives-ouvertes.fr/insu-00952537/file/russian.pdf
DOI : 10.1002/2013WR014255

Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
Two types of olivine occur in kimberlites from Greenland, Canada and southern Africa. The first, xenocrystic olivine, displays several different forms. Most distinctive are 'nodules', a term we use to describe the large (1-5 mm), rounded, single crystals or polycrystalline aggregates that are a common constituent of many kimberlites. Olivine compositions are uniform within single nodules but vary widely from nodule to nodule, from Fo81 to 93. Within many nodules, sub- to euhedral asymmetric tablets have grown within larger anhedral olivine grains. Dislocation structures, particularly in the anhedral grains, demonstrate that this olivine was deformed before being incorporated into the kimberlite magma. Olivine grains in the kimberlite matrix between the nodules have morphologies similar to those of the tablets, suggesting that most matrix olivine grains are parts of disaggregated nodules. We propose that a sub- to euhedral form is not sufficient to identify phenocrysts in kimberlites and provide some criteria, based on morphology, internal deformation and composition, that distinguish phenocrysts from xenocrysts. The second type of olivine is restricted to rims surrounding xenocrystic grains. Only this olivine crystallized from the kimberlite magma. Major and trace element data for the rim olivine are used to calculate the composition of the parental kimberlite liquid, which is found to contain between about 20 and 30% MgO. The bulk compositions of many kimberlites contain higher MgO contents as a result of the presence of xenocrystic olivine. The monomineralic, dunitic, character of the nodules, and the wide range from Fo-rich to Fo-poor olivine compositions, provide major constraints on the origin of the nodules. Dunite is a relatively rare rock in the mantle and where present its olivine is persistently Fo-rich. The dunitic source of the nodules in kimberlites lacked minerals such as pyroxene and an aluminous phase, which make up about half of most mantle-derived rocks. It appears that these minerals were removed from the mantle peridotite that was to become the source of the nodules, and the Fo content of the retained olivine was modified during interaction with CO2-rich fluids whose arrival at the base of the lithosphere immediately preceded the passage of the kimberlite magmas. Fragments of the resultant dunite were entrained into the kimberlite, where they were retained both as intact nodules and as disaggregated grains in the matrix.
ISSN: 0022-3530

insu-00543321
https://hal-insu.archives-ouvertes.fr/insu-00543321
BIBCODE : 10.1093/petrology/egp080
DOI : 10.1093/petrology/egp080

Éditeur(s) : HAL CCSD Springer Verlag
Résumé : International audience
Random walk (RW) or Continuous Time Random Walk (CTRW) are recurring Monte Carlo methods used to model convective and diffusive transport in complex heterogeneous media. Many applications can be found, including fluid mechanic, hydrology, and chemical reactors modeling. These methods are easy to implement, very versatile and flexible enough to become appealing for many applications because they generally overlook of deeply simplify the building of explicit complex meshes required by deterministic methods. RW and CTRW provide a good physical understanding of the interactions between the space scales of het-erogeneities and the transport phenomena under consideration. In addition, they can result in efficient up-scaling methods, especially in context of flow and transport in fractured media. In the present study, we review the applications of RW or CTRW for several situations coping with various spatial scales, and different insights into up-scaling applications. RW and CTRW advantages and downsides are also discussed, thus providing a few avenues for further works and applications.
ISSN: 0169-3913

hal-01449131
https://hal.archives-ouvertes.fr/hal-01449131
https://hal.archives-ouvertes.fr/hal-01449131/document
https://hal.archives-ouvertes.fr/hal-01449131/file/0025972-02.pdf
DOI : 10.1007/s11242-016-0693-z

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
Both flow field heterogeneity and mass transfer between mobile and immobile domains have been studied separately for explaining observed anomalous transport. Here we investigate non-Fickian transport using high-resolution 3-D X-ray microtomographic images of Berea sandstone containing microporous cement with pore size below the setup resolution. Transport is computed for a set of representative elementary volumes and results from advection and diffusion in the resolved macroporosity (mobile domain) and diffusion in the microporous phase (immobile domain) where the effective diffusion coefficient is calculated from the measured local porosity using a phenomenological model that includes a porosity threshold ( inline image) below which diffusion is null and the exponent n that characterizes tortuosity-porosity power-law relationship. We show that both flow field heterogeneity and microporosity trigger anomalous transport. Breakthrough curve (BTC) tailing is positively correlated to microporosity volume and mobile-immobile interface area. The sensitivity analysis showed that the BTC tailing increases with the value of inline image, due to the increase of the diffusion path tortuosity until the volume of the microporosity becomes negligible. Furthermore, increasing the value of n leads to an increase in the standard deviation of the distribution of effective diffusion coefficients, which in turn results in an increase of the BTC tailing. Finally, we propose a continuous time random walk upscaled model where the transition time is the sum of independently distributed random variables characterized by specific distributions. It allows modeling a 1-D equivalent macroscopic transport honoring both the control of the flow field heterogeneity and the multirate mass transfer between mobile and immobile domains.
ISSN: 0043-1397

hal-01278285
https://hal.archives-ouvertes.fr/hal-01278285
DOI : 10.1002/2015WR017645

Éditeur(s) : HAL CCSD Springer Verlag
Résumé : International audience
Random walk (RW) methods are recurring Monte Carlo methods used to modelconvective and diffusive transport in complex heterogeneous media. Many applications canbe found, including fluid mechanic, hydrology and chemical reactors modeling. These methodsare easy to implement, very versatile and flexible enough to become appealing for manyapplications because they generally overlook or deeply simplify the building of explicit complexmeshes required by deterministic methods. RWprovides a good physical understandingof the interactions between the space scales of heterogeneities and the transport phenomenaunder consideration. In addition, they can result in efficient upscaling methods, especiallyin the context of flow and transport in fractured media. In the present study, we review theapplications of RW to several situations that cope with diverse spatial scales and differentinsights into upscaling problems. The advantages and downsides of RW are also discussed,thus providing a few avenues for further works and applications.
ISSN: 0169-3913

insu-01308060
https://hal-insu.archives-ouvertes.fr/insu-01308060
DOI : 10.1007/s11242-016-0693-z

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
We study diffusion in a heterogeneous medium that is characterized by spatially varying diffusion properties from a random walk point of view. We show that an inhomogeneous continuous time random walk (CTRW) with a spatially variable exponential transition time distribution solves the spatially discretized heterogeneous diffusion equation. This demonstrates the equivalence of the widely used time-domain random walk (TDRW) scheme and spatially inhomogeneous CTRW and at the same time provides a demonstration of the formal equivalence of the TDRW particle formulation and the heterogeneous diffusion equation. Based on this equivalence, we develop a TDRW method for heterogeneous diffusion under spatially variable multirate mass transfer properties. We discuss the implementation of these schemes and study the diffusion behavior in the presence of traps that are characterized by a truncated power-law trapping time distribution.
ISSN: 0309-1708

hal-00767566
https://hal.archives-ouvertes.fr/hal-00767566
DOI : 10.1016/j.advwatres.2012.07.015

Éditeur(s) : HAL CCSD AGU and the Geochemical Society
Résumé : International audience
Hydrothermal circulation at ultramafic-hosted sites supports a large variety of high-and low-temperature hydrothermal vents and associated ecosystems. The discovery of abundant fossil vesicomyid and thyasirid shell accumulations at the ridge crest, approximately 2.5 km east of the active Rainbow vent field on the Mid-Atlantic Ridge (MAR, 36 degrees 13'N), increased our knowledge regarding the diversity of vent communities at slow spreading ridges. Bivalve molluscs of the family Vesicomyidae were represented by the genus Phreagena. Here we present the first record of this genus in the Atlantic Ocean. This second vesicomyid species known from the MAR, Phreagena sp., was found to be associated with a Thyasira species that is affiliated with T. southwardae (at the Logatchev vent field on the MAR) and with T. vulcolutre (in the Gulf of Cadiz). These two clams have close relationships with seep taxa along the continental margin, and were likely associated with sedimented vent fields. delta O-18 and delta C-13 analyses of the shells suggested that the burrowing bivalve Thyasira could incorporate isotopically light carbon, derived from the oxidation of methane in the sediment, while the signature of Phreagena sp. shells denoted a different carbonate source. C-14 dating of the shells denoted that the hydrothermal activity in the Rainbow area began at least similar to 25.5 kyr BP, which is similar to the model of the hydrothermal vent field distribution that was proposed for the Logatchev hydrothermal site. The results provide new insight regarding the diversity of chemosynthetic fauna on the MAR over geologic time. Ultramafic-hosted, on-axis sedimented vent fields extend the range of habitats for chemosynthetic communities, underlying the need to further explore the geology of these types of environments on slow-spreading ridges and to determine their role in the ecology of deep-sea vent communities.
ISSN: 1525-2027

hal-00523467
https://hal.archives-ouvertes.fr/hal-00523467
https://hal.archives-ouvertes.fr/hal-00523467/document
https://hal.archives-ouvertes.fr/hal-00523467/file/G3-Lartaud-_-al.pdf
DOI : 10.1029/2010GC003079

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
We derive a general formulation of the time domain random walk (TDRW) approach to model the hydrodynamic transport of inert solutes in complex geometries and heterogeneous media. We demonstrate its formal equivalence with the discretized advection-dispersion equation and show that the TDRW is equivalent to a continuous time random walk (CTRW) characterized by space-dependent transition times and transition probabilities. The transition times are exponentially distributed. We discuss the implementation of different concentration boundary conditions and initial conditions as well as the occurrence of numerical dispersion. Furthermore, we propose an extension of the TDRW scheme to account for mobile-immobile multirate mass transfer. Finally, the proposed TDRW scheme is validated by comparison to analytical solutions for spatially homogeneous and heterogeneous transport scenarios.
ISSN: 0043-1397

hal-01356319
https://hal.archives-ouvertes.fr/hal-01356319
DOI : 10.1002/2015WR018511

Éditeur(s) : HAL CCSD Wiley
Résumé : International audience
ISSN: 0028-646X

hal-01092529
https://hal-agroparistech.archives-ouvertes.fr/hal-01092529