Petrology, geochemistry and geodynamic implications of Jurassic island arc magmatism as revealed by mafic volcanic rocks in the Mesozoic low-grade sequence, eastern Rhodope, Bulgaria

In the eastern Bulgarian Rhodope, mafic extrusive rocks and underlying greenschists are found in the Mesozoic low-grade unit, which represents the northern extension of similar sequences including the Evros ophiolites in Thrace (Greece). Both rock types define a suite of low-Ti tholeiitic basalts to transitional boninitic basaltic andesites and andesites and associated metapyroclastites (greenschists), intruded at its base by diorite dikes of a boninitic affinity. Mafic lavas and greenschists display large ion lithophile element (LILE) enrichment relative to high-field strength elements (HFSE), flat REE patterns of a slight light REE depletion, a strong island arc tholeiite (IAT) and weak MORB-like signature. All these rocks are characterized by negative Nb anomalies ascribed to arc lavas. They have positive epsilon Nd(i) values in the range of +4.87 to +6.09, approaching the lower limit of MORB-like source, and relatively high ((207)Pb/(204)Pb)(i) (15.57-15.663) at low ((206)Pb/(204)Pb)(i) (18.13-18.54) ratios. The Nd isotopic compositions coupled with trace element data imply a dominantly depleted MORB-like mantle source and a contribution of subduction modified LILE-enriched component derived from the mantle wedge. The diorite dike has a low eNdi value of -2.61 and is slightly more Pb radiogenic ((207)Pb/(204)Pb)(i) (15.64) and ((206)Pb/(204)Pb)(i) (18.56), respectively, reflecting crustal contamination. Petrologic and geochemical data indicate that the greenschists and mafic extrusive rocks represent a magmatic assemblage formed in an island arc setting. The magmatic suite is interpreted as representing an island arc-accretionary complex related to the southward subduction of the Meliata-Maliac ocean under the supra-subduction back-arc Vardar ocean/island arc system. Magmatic activity appears to have initiated in the north during the inception of the island arc system by the Early-Middle Jurassic time in the eastern Rhodope that most likely graded to back-arc spreading southwards as represented by the Late Jurassic MORB-type Samothraki Island ophiolites. This tectonic scenario is further constrained by paleotectonic reconstructions. The arc-trench system collided with the Rhodope in the Late Jurassic times. (c) 2007 Elsevier B.V. All rights reserved.

Organic geochemistry of Jurassic-Cretaceous source rocks and oil seeps from the profile across the Adriatic-Dinaric carbonate platform

Organic geochemical and stable isotope investigations were performed to provide an insight into the depositional environments, origin and maturity of the organic matter in Jurassic and Cretaceous formations of the External Dinarides. A correlation is made among various parameters acquired from Rock-Eval, gas chromatography-mass spectrometry data and isotope analysis of carbonates and kerogen. Three groups of samples were analysed. The first group includes source rocks derived from Lower Jurassic limestone and Upper Jurassic ``Leme'' beds, the second from Upper Cretaceous carbonates, while the third group comprises oil seeps genetically connected with Upper Cretaceous source rocks. The carbon and oxygen isotopic ratios of all the carbonates display marine isotopic composition. Rock-Eval data and maturity parameter values derived from biomarkers define the organic matter of the Upper Cretaceous carbonates as Type I-S and Type II-S kerogen at the low stage of maturity up to entering the oil-generating window. Lower and Upper Jurassic source rocks contain early mature Type III mixed with Type IV organic matter. All Jurassic and Cretaceous potential source rock extracts show similarity in triterpane and sterane distribution. The hopane and sterane distribution pattern of the studied oil seeps correspond to those from Cretaceous source rocks. The difference between Cretaceous oil seeps and potential source rock extracts was found in the intensity and distribution of n-alkanes, as well as in the abundance of asphaltenes which is connected to their biodegradation stage. In the Jurassic and Cretaceous potential source rock samples a mixture of aromatic hydrocarbons with their alkyl derivatives were indicated, whereas in the oil seep samples extracts only asphaltenes were observed.

Trace element distribution among rock-forming minerals from metamorphosed to partially molten basic igneous rocks in a contact aureole (Fuerteventura, Canaries)

Low pressure partial melting of basanitic and ankaramitic dykes gave rise to unusual, zebra-like migmatites, in the contact aureole of a layered pyroxenite-gabbro intrusion, in the root zone of an ocean island (Basal Complex, Fuerteventura, Canary Islands). These migmatites are characterised by a dense network of closely spaced, millimetre-wide leucocratic segregations. Their mineralogy consists of plagioclase (An(32-36)), diopside, biotite, oxides (magnetite, ilmenite), +/-amphibole, dominated by plagioclase in the leucosome and diopside in the melanosome. The melanosome is almost completely recrystallised, with the preservation of large, relict igneous diopside phenocrysts in dyke centres. Comparison of whole-rock and mineral major- and trace-element data allowed us to assess the redistribution of elements between different mineral phases and generations during contact metamorphism and partial melting. Dykes within and outside the thermal aureole behaved like closed chemical systems. Nevertheless, Zr, Hf, Y and REEs were internally redistributed, as deduced by comparing the trace element contents of the various diopside generations. Neocrystallised diopside - in the melanosome, leucosome and as epitaxial phenocryst rims - from the migmatite zone, are all enriched in Zr, Hf, Y and REEs compared to relict phenocrysts. This has been assigned to the liberation of trace elements on the breakdown of enriched primary minerals, kaersutite and sphene, on entering the thermal aureole. Major and trace element compositions of minerals in migmatite melanosomes and leucosomes are almost identical, pointing to a syn- or post-solidus reequilibration on the cooling of the migmatite terrain i.e. mineral-melt equilibria were reset to mineral-mineral equilibria. (C) 2007 Elsevier B.V. All rights reserved.

Origin of cadmium enrichments in carbonate rocks deposited in the Alpine Tethys area during the Middle-Late Jurassic

ABSTRACTThe pollution of air, soil and water by heavy metals through anthropogenic activities is an object of numerous environmental studies since long times. A number of natural processes, such as volcanic activity, hydrothermal fluid circulation and weathering of metal-rich deposits may lead to an additional and potentially important input and accumulation of heavy metals in the environment. In the Swiss and French Jura Mountains, anomalous high cadmium (Cd) concentrations (up to 16 ppm) in certain soils are related to the presence of underlying Cd-enriched (up to 21 ppm) carbonate rocks of Middle to Late Jurassic age. The aim of this study is to understand the processes controlling Cd incorporation into carbonate rocks of Middle and Late Jurassic age and to reconstruct the sedimentary and environmental conditions, which have led to Cd enrichments in these sedimentary rocks.Cd concentrations in studied hemipelagic sections in France vary between 0.1 and 0.5 ppm (mean 0.15 ppm). Trace-element behavior and high Mn concentrations suggest that sediment accumulation occurred in a well-oxygenated environment. Increases in Cd contents in the bulk-rock carbonate sediments may be related to increases in surface-water productivity under oxic conditions and important remineralization of organic matter within the water column. In platform settings preserved in the Swiss Jura Mountains, no correlation is observed between Cd contents and evolution of environmental conditions. Cd concentrations in these platform sections are often below the detection limit, with isolated peaks of up to 21 ppm. These important Cd enrichments are associated with peaks in Zn concentrations and are present in carbonate rocks independently of facies and age. The high Cd contents in these shallow-water carbonate rocks are partly related to the presence of disseminated, Cd-rich (up to 1.8%), sphalerite (ZnS) mineralization. The basement rocks are considered to be the source of metals for sulfide mineralization in the overlying Jurassic strata, as the sphalerite Pb isotope pattern is comparable to that of granite rocks from the nearby southern Black Forest crystalline massif. The Rb-Sr ages of sphalerite samples indicate that a main phase of sphalerite formation occurred near the boundary between the late Middle and early Late Jurassic, at around 162 Ma, as a result of enhanced tectonic and hydrothermal activity in Europe, related to the opening of the Central Atlantic and to the tectonic/thermal subsidence during spreading of the Alpine Tethys. I therefore propose to use unusually high Cd concentrations in carbonates as a tracer of tectonic activity in the Jura Mountains area, especially in the case when important enrichments in Zn co-occur. Paleoproductivity reconstructions based on records of authigenic Cd may be compromised not only by post-depositional redistribution, but also by incorporation of additional Cd from hydrothermal solutions circulating in the rock.The circulation of metal-rich hydrothermal fluids through the sediment sequence, in addition to specific environmental conditions during sedimentation, contributes to the incorporation of Cd into the carbonate rocks. However, only hydrothermal activity has led to the unusually high concentrations of Cd in carbonate rocks of Bajocian-Oxfordian age, through the formation of sphalerite mineralization.

Jurassic and Early Cretaceous Radiolaria of the Lower Amurian Terrane: Khabarovsk region, far east of Russia

New data on biostratigraphy, sedimentology and tectonics of the Russian Far Eastern region (Lower Amurian terrane) are presented. This study shows that sedimentary sequence of the terrane consists of interbedded Radiolaria-bearing siliceous and volcaniclastic sediments spanning an interval of over 90 million years. It is shown that accumulation of radiolarian deposits on an oceanic plate was associated with alkaline (intraplate) volcanism in the Jurassic, while the plate was drifting, and with some are volcanism during the Early Cretaceous. The younger siliceous rocks contain volcaniclastic material and indicate that the studied sequence approached the trench in the Early Cretaceous (Hauterivian-Barremian) and became accreted in the late Albian-early Cenomanian. We describe and illustrate radiolarian species extracted fi om 21 samples. A taxonomic list of 194 taxa and nine plates of Jurassic-Early Cretaceous Radiolaria are presented.

Late mesozoic to neogene radiolarian biostratigraphy and palaeoceanography in the Caribbean and East pacific region

^Raduolarians constitute a good tool for contributing to the biostratigraphy of accreted terranes and in deep-sea sediment sequences. The use of radiolarians is also proven to be valuable as a palaeoceanographic indicator. The present study evaluates radiolarians in three different geological settings, in order to better constrain the age of the sites and to try to understand their palaeoenvironmental situation at different periods, particularly in the Caribbean-Central America area. On the Jarabacoa Block, in Central Dominican Republic, a hundred meters of siliceous mudstones (Pedro Brand section in the Tireo Group) was dated as Turonian- Coniancian in age using radiolarians. A 40Ar-39Ar whole rock age of 75.1±1.1 Ma (Campanian), obtained in a basalt dyke crosscutting the radiolarian bearing rocks, a consistent minimum age for the pelagic-hemipelagic Pedro Brand section. The Jarabacoa Block is considered as the most complete outcrop section of Pacific ocean crust overlain by a first Aptian-Albian phase of Caribbean Large Igneous Province-type activity (CLIP), followed by the development of a Cenomanian-Santonian intraoceanic arc, which is in turn overlain by a late Campanian-Maastrichtian CLIP-phase. The Tireo Group records an episode of pelagic to hemi-pelagic and intermediate to acidic arc-derived sedimentation, previous to the youngest magmatic phase of the CLIP. Thus, the section of Pedro Brand has been interpreted in this study as being part of the intraoceanic arc. In northern Venezuela, a greenish radiolarite section from Siquisique Ophiolite (basalts, gabbros and some associated cherts) in Guaparo Creek has been studied. In previous studies, the Ophiolite unit (Petacas Creek section) has been dated as Bajocian-Bathonian, based on ammonites present in interpillow sediments from basalt blocks. New dating of the present study concluded in an Aptian?-Albian-Cenomanian age for the Guaparo creek section (middle Cretaceous), based on radiolarian assemblage associated to basalts-gabbros rocks of the unit. Previous plagioclase 40Ar-39Ar ages from the Siquisuique Ophiolite may be slightly younger (94-90 Ma.) and may, therefore, represent younger dykes that intruded onto a well-developed sheeted dyke complex of the Siquisique. The geochemistry of these rocks and the palaeotectonic reconstruction of the Caribbean area during this period suggest that these rocks were derived from a mid-ocean ridge with an influence of deep mantle plume. The Siquisique Ophiolite most probably represents a fragment of the proto-Caribbean basin. The Integrated Ocean Drilling Program Expedition 344 drilled a transect across the convergent margin off Costa Rica. Two sites of this expedition were chosen for radiolarian biostratigraphy and palaeoceanographic studies. Both sites (U1381C and U1414A) are located in the incoming Cocos plate, in the eastern Equatorial Pacific. The succession of U1381C yields a Middle Miocene to Pleistocene age, and presents an important hiatus of approximately 10 Ma. The core of U1414A exposes a continuous sequence that deposited during Late Miocene to Pleistocene (radiolarian zones RN6-RN16). The ages were assigned based on radiolarians and correlated with nannofossil zonation and tephra 40Ar-39Ar datation. With those results, and considering the northward movement of the Cocos plate motion (about 7 cm/year), deduction is made that the sites U1381C and U1414A were initially deposited during the Miocene, several hundreds of kilometres from the current location, slightly south of the Equator. This suggests that the faunas of these sites have been subjected to different currents, first influenced by the cold tongue of the South Equatorial Current and followed by the warm Equatorial Countercurrent. At last, coastal upwelling influenced faunas of the Pleistocene. -- Les radiolaires sont considérés comme un outil utile à la biostratigraphie des terrains accrétés et des sédiments profonds. Leur utilité est aussi prouvée comme étant remarquable au niveau des reconstructions paléocéanographiques. La présente étude évalue l'importance et la présence des radiolaires de trois localités géologiquement différentes d'Amérique Centrale-Caraïbes, dans le but d'améliorer les model d'âges et de mieux comprendre la situation paléoenvironnementale à travers le temps. Dans le Bloque de Jarabacoa, au centre de la République Dominicaine, une section de cent mètres (section de Pedro Brand, Groupe de Tireo) a été datée comme faisant partie du Turonien-Santonien, en utilisant les radiolaires. Une datation 40Ar-39Ar sur roche totale de 75±1.1 Ma (Campanien) a été obtenu pour vin dyke traversant les sédiments riches en radiolaires, en cohérence avec l'âge minimum accordé à la section de Pedro Brand. Aux Caraïbes, le Bloque de Jarabacoa est considéré comme l'affleurement le plus complet présentant une succession de croûte océanique d'origine Pacifique recouverte d'une première phase d'activité volcanique de type CLIP (Caribbean Large Igneous Province) d'âge Aptien- Albien, de dépôts d'arc volcanique intra-océanique d'âge Cénomanien-Santonien, puis d'une seconde phase de type CLIP d'âge Campanien-Maastrichtien. Le Groupe de Tireo enregistre un épisode de dépôt pélagiques-hémipélagiques et d'arc volcanique, antérieur à la plus jeune phase de type CLIP. Cette étude place donc la formation de la section de Pedro Brand au moment du développement de l'arc intra-océanique. A Guaparo Creek (nord du Vénézuela), une section de radiolarite verdâtre faisant partie des ophiolites de Siquisique (basaltes, gabbros, cherts) a été étudiée. Dans des études précédentes, sur la localité de Petacas Creek, l'unité ophiolitique a été daté d'âge Bajocien- Bathonien (Jurassique) sur la base d'ammonites trouvées dans des sédiments intercalés entre des laves en coussins. Les nouvelles datations de notre étude, basées sur des assemblages à radiolaires de l'unité à basaltes-gabbros, donnent un âge Aptien?-Albien-Cénomanien (Crétacé moyen). Les âges de l'Ophiolite de Siquisique, précédement calculés par la méthode sur plagioclases, pourraient être légèrement plus jeune (94-90 Ma) et donc représenter des intrusions plus récentes de dykes dans le complexe filonien déjà bien dévelopé. La géochimie de ces roches magmatiques, ainsi que les reconstructions paléotectoniques de la zone Caraïbes durant cette période, suggèrent que ces formations sont dérivées d'une ride médio-océanique associée à l'influence d'un panache mantellique. L'ophiolite de Siquisique représente très probablement un fragment du bassin de proto¬Caraïbe. L'expédition 344 du programme IODP (Integrated Ocean Drilling Program) a eu lieu dans l'optique de forer et dresser une coupe de la marge convergente au large du Costa Rica. Deux sites de cette expédition ont été choisis pour les besoins des études de biostratigraphie et de reconstruction paléocéanographique. Ces deux sites (U1381C et U1414A) sont situés sur la plaque subductante de Cocos, dans la zone Pacifique est-équatoriale. La carotte U1381C expose une séquence s'étalant du Miocène moyen au Pléistocène, et présente un important hiatus d'environ 10 Ma. La carotte U1414A expose une séquence continue s'étalant du Miocène tardif au Pléistocène (zone à radiolaires RN6-RN16). Les âges ont été assignés sur la base des radiolaires et corrélés avec les zones à nanofossiles et les datations 40Ar-39Ar sur téphras. Avec ces résultats, et en considérant le mouvement nord de la plaque de Cocos (environ 7 cm/an), déduction est faite que les deux sites étaient initialement situés, au cours du Miocène, à plusieurs centaines de kilomètres de leur location actuelle, au sud de l'équateur. Cela suggère que les faunes de ces sites ont été sujettes à différents courants; premièrement influencées par la langue froide du SEC (South Equatorial Current), puis par les eaux chaudes du ECC (Equatorial Countercurrent). Pour terminer, les remontées d'eau côtières ont influencées les faunes Pléistocène.

Seismic characterization of heterogeneous sedimentary and fractured rocks based on Biot's theory of poroelasticity

Understanding and quantifying seismic energy dissipation, which manifests itself in terms of velocity dispersion and attenuation, in fluid-saturated porous rocks is of considerable interest, since it offers the perspective of extracting information with regard to the elastic and hydraulic rock properties. There is increasing evidence to suggest that wave-induced fluid flow, or simply WIFF, is the dominant underlying physical mechanism governing these phenomena throughout the seismic, sonic, and ultrasonic frequency ranges. This mechanism, which can prevail at the microscopic, mesoscopic, and macroscopic scale ranges, operates through viscous energy dissipation in response to fluid pressure gradients and inertial effects induced by the passing wavefield. In the first part of this thesis, we present an analysis of broad-band multi-frequency sonic log data from a borehole penetrating water-saturated unconsolidated glacio-fluvial sediments. An inherent complication arising in the interpretation of the observed P-wave attenuation and velocity dispersion is, however, that the relative importance of WIFF at the various scales is unknown and difficult to unravel. An important generic result of our work is that the levels of attenuation and velocity dispersion due to the presence of mesoscopic heterogeneities in water-saturated unconsolidated clastic sediments are expected to be largely negligible. Conversely, WIFF at the macroscopic scale allows for explaining most of the considered data while refinements provided by including WIFF at the microscopic scale in the analysis are locally meaningful. Using a Monte-Carlo-type inversion approach, we compare the capability of the different models describing WIFF at the macroscopic and microscopic scales with regard to their ability to constrain the dry frame elastic moduli and the permeability as well as their local probability distribution. In the second part of this thesis, we explore the issue of determining the size of a representative elementary volume (REV) arising in the numerical upscaling procedures of effective seismic velocity dispersion and attenuation of heterogeneous media. To this end, we focus on a set of idealized synthetic rock samples characterized by the presence of layers, fractures or patchy saturation in the mesocopic scale range. These scenarios are highly pertinent because they tend to be associated with very high levels of velocity dispersion and attenuation caused by WIFF in the mesoscopic scale range. The problem of determining the REV size for generic heterogeneous rocks is extremely complex and entirely unexplored in the given context. In this pilot study, we have therefore focused on periodic media, which assures the inherent self- similarity of the considered samples regardless of their size and thus simplifies the problem to a systematic analysis of the dependence of the REV size on the applied boundary conditions in the numerical simulations. Our results demonstrate that boundary condition effects are absent for layered media and negligible in the presence of patchy saturation, thus resulting in minimum REV sizes. Conversely, strong boundary condition effects arise in the presence of a periodic distribution of finite-length fractures, thus leading to large REV sizes. In the third part of the thesis, we propose a novel effective poroelastic model for periodic media characterized by mesoscopic layering, which accounts for WIFF at both the macroscopic and mesoscopic scales as well as for the anisotropy associated with the layering. Correspondingly, this model correctly predicts the existence of the fast and slow P-waves as well as quasi and pure S-waves for any direction of wave propagation as long as the corresponding wavelengths are much larger than the layer thicknesses. The primary motivation for this work is that, for formations of intermediate to high permeability, such as, for example, unconsolidated sediments, clean sandstones, or fractured rocks, these two WIFF mechanisms may prevail at similar frequencies. This scenario, which can be expected rather common, cannot be accounted for by existing models for layered porous media. Comparisons of analytical solutions of the P- and S-wave phase velocities and inverse quality factors for wave propagation perpendicular to the layering with those obtained from numerical simulations based on a ID finite-element solution of the poroelastic equations of motion show very good agreement as long as the assumption of long wavelengths remains valid. A limitation of the proposed model is its inability to account for inertial effects in mesoscopic WIFF when both WIFF mechanisms prevail at similar frequencies. Our results do, however, also indicate that the associated error is likely to be relatively small, as, even at frequencies at which both inertial and scattering effects are expected to be at play, the proposed model provides a solution that is remarkably close to its numerical benchmark. -- Comprendre et pouvoir quantifier la dissipation d'énergie sismique qui se traduit par la dispersion et l'atténuation des vitesses dans les roches poreuses et saturées en fluide est un intérêt primordial pour obtenir des informations à propos des propriétés élastique et hydraulique des roches en question. De plus en plus d'études montrent que le déplacement relatif du fluide par rapport au solide induit par le passage de l'onde (wave induced fluid flow en anglais, dont on gardera ici l'abréviation largement utilisée, WIFF), représente le principal mécanisme physique qui régit ces phénomènes, pour la gamme des fréquences sismiques, sonique et jusqu'à l'ultrasonique. Ce mécanisme, qui prédomine aux échelles microscopique, mésoscopique et macroscopique, est lié à la dissipation d'énergie visqueuse résultant des gradients de pression de fluide et des effets inertiels induits par le passage du champ d'onde. Dans la première partie de cette thèse, nous présentons une analyse de données de diagraphie acoustique à large bande et multifréquences, issues d'un forage réalisé dans des sédiments glaciaux-fluviaux, non-consolidés et saturés en eau. La difficulté inhérente à l'interprétation de l'atténuation et de la dispersion des vitesses des ondes P observées, est que l'importance des WIFF aux différentes échelles est inconnue et difficile à quantifier. Notre étude montre que l'on peut négliger le taux d'atténuation et de dispersion des vitesses dû à la présence d'hétérogénéités à l'échelle mésoscopique dans des sédiments clastiques, non- consolidés et saturés en eau. A l'inverse, les WIFF à l'échelle macroscopique expliquent la plupart des données, tandis que les précisions apportées par les WIFF à l'échelle microscopique sont localement significatives. En utilisant une méthode d'inversion du type Monte-Carlo, nous avons comparé, pour les deux modèles WIFF aux échelles macroscopique et microscopique, leur capacité à contraindre les modules élastiques de la matrice sèche et la perméabilité ainsi que leur distribution de probabilité locale. Dans une seconde partie de cette thèse, nous cherchons une solution pour déterminer la dimension d'un volume élémentaire représentatif (noté VER). Cette problématique se pose dans les procédures numériques de changement d'échelle pour déterminer l'atténuation effective et la dispersion effective de la vitesse sismique dans un milieu hétérogène. Pour ce faire, nous nous concentrons sur un ensemble d'échantillons de roches synthétiques idéalisés incluant des strates, des fissures, ou une saturation partielle à l'échelle mésoscopique. Ces scénarios sont hautement pertinents, car ils sont associés à un taux très élevé d'atténuation et de dispersion des vitesses causé par les WIFF à l'échelle mésoscopique. L'enjeu de déterminer la dimension d'un VER pour une roche hétérogène est très complexe et encore inexploré dans le contexte actuel. Dans cette étude-pilote, nous nous focalisons sur des milieux périodiques, qui assurent l'autosimilarité des échantillons considérés indépendamment de leur taille. Ainsi, nous simplifions le problème à une analyse systématique de la dépendance de la dimension des VER aux conditions aux limites appliquées. Nos résultats indiquent que les effets des conditions aux limites sont absents pour un milieu stratifié, et négligeables pour un milieu à saturation partielle : cela résultant à des dimensions petites des VER. Au contraire, de forts effets des conditions aux limites apparaissent dans les milieux présentant une distribution périodique de fissures de taille finie : cela conduisant à de grandes dimensions des VER. Dans la troisième partie de cette thèse, nous proposons un nouveau modèle poro- élastique effectif, pour les milieux périodiques caractérisés par une stratification mésoscopique, qui prendra en compte les WIFF à la fois aux échelles mésoscopique et macroscopique, ainsi que l'anisotropie associée à ces strates. Ce modèle prédit alors avec exactitude l'existence des ondes P rapides et lentes ainsi que les quasis et pures ondes S, pour toutes les directions de propagation de l'onde, tant que la longueur d'onde correspondante est bien plus grande que l'épaisseur de la strate. L'intérêt principal de ce travail est que, pour les formations à perméabilité moyenne à élevée, comme, par exemple, les sédiments non- consolidés, les grès ou encore les roches fissurées, ces deux mécanismes d'WIFF peuvent avoir lieu à des fréquences similaires. Or, ce scénario, qui est assez commun, n'est pas décrit par les modèles existants pour les milieux poreux stratifiés. Les comparaisons des solutions analytiques des vitesses des ondes P et S et de l'atténuation de la propagation des ondes perpendiculaires à la stratification, avec les solutions obtenues à partir de simulations numériques en éléments finis, fondées sur une solution obtenue en 1D des équations poro- élastiques, montrent un très bon accord, tant que l'hypothèse des grandes longueurs d'onde reste valable. Il y a cependant une limitation de ce modèle qui est liée à son incapacité à prendre en compte les effets inertiels dans les WIFF mésoscopiques quand les deux mécanismes d'WIFF prédominent à des fréquences similaires. Néanmoins, nos résultats montrent aussi que l'erreur associée est relativement faible, même à des fréquences à laquelle sont attendus les deux effets d'inertie et de diffusion, indiquant que le modèle proposé fournit une solution qui est remarquablement proche de sa référence numérique.