Éditeur(s) : HAL CCSD Oxford University Press (OUP)
Résumé : International audience
Several studies, both theoretical and experimental, show that sedimentary rocks have a fractal pore-grain interface. In this paper a computer simulated 3-D sedimentary rock structure generated by the Relaxed Ballistic Bidisperse Deposition Model (RBBDM), is investigated to characterize the micro structure of its pores. The pore volume and the rock-pore interface show the same fractal dimension indicating that the pore volume is a fractal. The two point density correlation is computed for the pore space and the results compare favourably with the range reported from experiments. An array of 2-D X-ray tomography micrograph sections of a real sedimentary rock, an oolitic limestone (pure calcite) from the Mondeville formation of Middle Jurassic age (Paris Basin, France), was used to generate a 3-D bitmap. The 3-D real rock sample generated in this manner, was analysed for similar studies as the simulated structure. The results were compared with those obtained from simulation. The simulation results agree qualitatively with the real rock sample. Diffusion through the connected pore space of the simulated structure was studied using a random walk algorithm and the results compared with the similar simulation study done on the 3-D oolitic limestone specimen. In both cases diffusion was found to be anomalous indicating that the sedimentary rock has a fractal geometry. The favourable comparability of results between the simulated and real rock supports the usefulness of the model of sedimentary rock generation which can be applicable to transport phenomena.
ISSN: 0956-540X

hal-00823061
https://hal.archives-ouvertes.fr/hal-00823061
DOI : 10.1093/gji/ggs084

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
Several studies, both theoretical and experimental, show that sedimentary rocks have fractal pore-grain interface. The real rocks are 3-d structures with highly tortuous and often fractal pore spaces. Before attempting simulation of this daunting geometry, we present here, as a preliminary study, a simpler 2-d version. In this paper, a computer simulated 2-d sedimentary rock structure is generated by the relaxed bidisperse ballistic deposition model. Grains of two different sizes are dropped ballistically on a linear substrate. By changing the fraction of the two types of particles, the porosity of the rock structure can be tuned. The structure undergoes compaction through the relaxation of possible unstable overhangs. The micro structure of the pore space is investigated. The pore mass and the rock-pore interface show a fractal behaviour with the same fractal dimension indicating that the pore volume is a fractal. Our simulation results indicate that the process of compaction of grains during the deposition process seems to erase the dependency of the fractal dimension on the grain size distribution. The two point density correlation is measured for the pore space. It shows anisotropy which is an outcome of the growth rule. X-ray tomography of two-dimensional sections of real sedimentary rocks obtained from Mallorca Island is subjected to the same study and the results compared with those obtained from simulation. The simulation results agree qualitatively with the real rock sample. We also study diffusion on the pore space. Diffusion is found to be anomalous as is expected in fractal spaces. It also bears the signature of anisotropy of the structure. Diffusion studies on the real rock sample could not yield conclusive results as the system size is not large enough.
ISSN: 0926-9851

hal-00761818
https://hal.archives-ouvertes.fr/hal-00761818
DOI : 10.1016/j.jappgeo.2012.09.002

Éditeur(s) : HAL CCSD American Geophysical Union
Résumé : International audience
An analytical solution in the Laplace domain is derived for modeling anomalous pressure diffusion during pumping tests in aquifer displaying hierarchical fractal fracture networks. The proposed solution generalizes all of the analytical models for fractal flow published previously by combining multifractal diffusion and nested multiporosity with transient exchanges, interface skin effects, and well storage effect. Solutions are derived for fracture-delimited blocks with planar, cylindrical, and spherical shapes, as well as with any fractional dimensional shapes. Any combinations of these shapes can be defined in order to model a large range of situations. Within each permeability level, the fractal properties of the fracture network can be specified and the fractal dimension can be distinct from the shape dimension of the block.
ISSN: 0043-1397

hal-00411549
https://hal.archives-ouvertes.fr/hal-00411549
DOI : 10.1029/2008WR007125

Éditeur(s) : HAL CCSD
Résumé : International audience
Mylonite bodies are exhumed extensions at depth of shallow structural level brittle faults. They are characterised by intense strain localisation involving large volumes of middle to lower continental crust. The analysis of shear zone geometry and microstructures is a useful tool to decipher deformation kinematics, though finite strain estimation is still relatively qualitative. Timing of the deformation has also been the focus of numerous studies using the complete spectrum of geochronometers. However, only rare examples provide geochronological data that can be unambiguously interpreted and cooling, recrystallisation and neocrystallisation are difficult to tell from each other. Thus, time-scales over which major mylonite zones develop and remain active under ductile conditions as well as their strain rate often remain poorly documented. We investigated two sections across the South-Armorican Shear Zone (SASZ), a crustal-scale pluri-km thick dextral shear zone active during a single tectonic event of late-Variscan age and associated with syntectonic granite emplacement.. Finite strain profiles on these cross-sections, quantified from average Quartz grain-size in 34 samples, show a strong deformation partitioning towards the center of the SASZ, with a decrease, for each section, from 150µm in slightly deformed zones to less than 7-5 µm in the core of the ultramylonite. Furthermore, toward the core of the SASZ, development of a penetrative shear band network progressively eliminates lenses of less deformed rocks and results in homogenisation of the Quartz grain sizes. In parallel, we performed in-situ UV-laser Ar/Ar dating on the fabric-forming white micas. White micas population can be decomposed into inherited muscovites carried by the folation and newly-formed, syn-tectonic substituted phengites located along the shear bands yielding 2 age groups for each sample. Along the profiles, the difference between these two ages evolves from 10-15 Ma in weakly deformed domains to zero toward the SASZ core. In addition, mylonite-forming phengites yielded the same age irrespectively of their position on the section or of the rock strain.. On the regional point of view, this study brings deformation-ages for the late-orogenic evolution of the South Armorican domain. Strain accumulated along the SASZ in a time-span extending to 300-298 Ma, i.e. contemporaneously with the activity of nearby major extensional shear zones. These results support a new model of crustal-scale shear zone activity, where the whole body of the shear zone is active throughout its history (defined as the largest age difference through the section), i.e. the actively deforming zone does not get progressively narrower. The strain distribution nevertheless evolves with time, with an increasingly large fraction partitioned in the ultramylonitic core.
EGU General Assembly
Vienne, Austria

insu-00839320
https://hal-insu.archives-ouvertes.fr/insu-00839320
BIBCODE : 2010EGUGA..12.9881A

Éditeur(s) : HAL CCSD Elsevier
Résumé : International audience
We develop a block model for the Aegean and surrounding areas, constrained by Global Positioning System (GPS), in order to investigate the degree of coupling on the Hellenic subduction interface (i.e., the fraction of the motion across the plate boundary accommodated by elastic strain accumulation). We use previously published models, and seismicity to define the geometry of the interface separating the down-going Nubian slab from the overriding Aegean. This model provides a good fit to the GPS observations; for the ∼200,000 km2 Aegean block the wrms of the residual velocities is 1.4 mm/yr for 80 GPS velocity estimates, approximately the 95% level of the GPS velocity uncertainties. We investigate the degree of coupling on the seismically active plate interface, the Hellenic trench splay fault (believed to be the source of the 365 AD Great Crete Earthquake and Tsunami), and the Kephalonia transform fault by comparing the modeled GPS residual velocity field for a range of coupling values. The GPS observations are almost insensitive to coupling on the Kephalonia transform fault, because of the vertical dip of the fault that creates interseismic deformation only close to the fault where few GPS sites exist. The absence of resolvable shortening of the leading edge of the Aegean Plate precludes coupling of more than 0.2 (20% of the full Nubia-Aegean convergence rate) on the modeled plate interface. Because of the shallow dip of the plate interface and trench splay fault, and high rate of convergence, if these boundaries were fully coupled, high elastic strain rates would be expected to extend well into the overriding Aegean plate. Based on our preferred value for the degree of coupling (0.1), and assuming characteristic earthquake behavior, we estimate a recurrence time for great earthquakes with slip similar to that for the 365 Crete event of 5700-8300 yr, consistent with the absence of subsequent great earthquakes on this segment of the subduction zone.
ISSN: 0012-821X

hal-01053486
https://hal.archives-ouvertes.fr/hal-01053486
DOI : 10.1016/j.epsl.2013.10.018

Éditeur(s) : HAL CCSD
Résumé : International audience
Mylonite-bearing crustal shear zones are exhumed extensions at depth of (potentially seismogenic) brittle faultsat shallow structural level. By concentrating deformation, shear zones play a major role in the rheology and mechanicalbehavior of the continental crust. Analysis of shear zone geometry and microstructures is a useful toolto decipher deformation kinematics, though finite strain estimation is still relatively qualitative. Timing of thedeformation has also been the focus of numerous studies using the complete spectrum of geochronometers. However,only rare examples provide geochronological data that can be unambiguously linked with deformation. Thus,time-scales over which major mylonite zones develop and remain active under ductile conditions as well as theirstrain rates often remain poorly documented. Several sections across the South Armorican Shear Zone (SASZ),a crustal-scale, several km-thick dextral shear zone were investigated by structural, petrological and Ar/Ar radiochronologicalmethods. Finite strain profiles on these sections showed a strong deformation partitioning towardsthe ultramylonite core of the SASZ, with a decrease, of both the C’/S mean angle (from ca. 45 to less than5) and the Quartz grain-size (from ca. 150 to 5 m). In parallel, conventional and in-situ Ar/Ar datings on thefabric-forming white micas were performed on compositionally complex white-micas. Inherited magmatic muscovitescarried by the foliation and newly-formed, syn-tectonic substituted phengites located along the shear bandsyielded 2 age groups for the less deformed samples. Along the shear zone cross-section, the difference betweenthese two ages evolves from 10-15 Ma in weakly deformed domains to zero toward the SASZ core. In addition, thevery fine-grained syn-tectonic phengites yielded the same ca. 300-298 Ma age irrespectively of their distance to thecore of the shear zone, even in slightly deformed rocks distant by several km. These results support a new modelof crustal-scale shear zone activity, where the whole body of the shear zone is active throughout its history, i.e.the actively deforming zone does not get progressively narrower. The strain distribution nevertheless evolves withtime, with an increasingly larger fraction partitioned in the ultramylonitic core. Besides, combination of the finitestrain estimations and the Ar/Ar ages enables us to assess the range of strain rates associated with mylonitization.
EGU General Assembly 2014
Vienne, Austria

insu-01565057
https://hal-insu.archives-ouvertes.fr/insu-01565057

Éditeur(s) : HAL CCSD
Résumé : International audience
Volumes of magma intruded in and emitted by volcanoes through time can be estimated by reconstruction of vol-cano's morphology and time sequence. Classical approaches for quantifying magma volumes on active volcanoes are based on the difference between pre-and post-eruption digital elevation models (DEM), but this kind of approach needs the pre-eruptive surfaces to be available. For old and eroded volcanoes these surfaces are poorly constrained. However, because the geometrical form of many volcanic edifices exhibits a remarkable symmetry we propose, here, a new approach using primary volcanic landforms in order to estimate the amount of the both erupted and eroded material and to locate eruptive centers. A large fraction of composite volcanoes have near constant slope on their flanks and a form that is concave upwards near their summits. But many phenomena can lead to non-symetrical edifices and complex morphologies can result, for example from parasitic centers of volcanism on the flanks, from alternation of short effusive and explosive construction phases, from flank or caldera collapses, or from glacial and other types of erosion. In this study we propose that, on the first order approximation, complex morphologies can be modeled by piling regular cones. In this model, cones centers and slopes are derived by fitting primary volcanic landform with a linear function :elevation=f(distance from center). Such an approach allows to estimate both errors on location of the eruptive center and on the volume of the resulting cones. This model can then be used for quantifying volume of erupted and eroded material, and for quantifying catastrophic events as giant landslides or flank collapse. This approach is tested on four different active volcanoes : Mount Mayon (Philippines), Mount Fuji (Japan), Mount Etna (Sicily) and Mount Teide (Canary Island) to estimate errors in volume between modeled and actual edifices. It is then used on volcanoes of La Réunion hotspot to reconstruct the Piton des Neiges and Piton de la Fournaise volcanoes at its different stages of growing.
EGU 2014
Vienna, Austria

hal-01386609
http://hal.univ-reunion.fr/hal-01386609