Lifelong dynamics of human CD4+CD25+ regulatory T cells: insights from in vivo data and mathematical modeling.

Despite their limited proliferation capacity, regulatory T cells (T(regs)) constitute a population maintained over the entire lifetime of a human organism. The means by which T(regs) sustain a stable pool in vivo are controversial. Using a mathematical model, we address this issue by evaluating several biological scenarios of the origins and the proliferation capacity of two subsets of T(regs): precursor CD4(+)CD25(+)CD45RO(-) and mature CD4(+)CD25(+)CD45RO(+) cells. The lifelong dynamics of T(regs) are described by a set of ordinary differential equations, driven by a stochastic process representing the major immune reactions involving these cells. The model dynamics are validated using data from human donors of different ages. Analysis of the data led to the identification of two properties of the dynamics: (1) the equilibrium in the CD4(+)CD25(+)FoxP3(+)T(regs) population is maintained over both precursor and mature T(regs) pools together, and (2) the ratio between precursor and mature T(regs) is inverted in the early years of adulthood. Then, using the model, we identified three biologically relevant scenarios that have the above properties: (1) the unique source of mature T(regs) is the antigen-driven differentiation of precursors that acquire the mature profile in the periphery and the proliferation of T(regs) is essential for the development and the maintenance of the pool; there exist other sources of mature T(regs), such as (2) a homeostatic density-dependent regulation or (3) thymus- or effector-derived T(regs), and in both cases, antigen-induced proliferation is not necessary for the development of a stable pool of T(regs). This is the first time that a mathematical model built to describe the in vivo dynamics of regulatory T cells is validated using human data. The application of this model provides an invaluable tool in estimating the amount of regulatory T cells as a function of time in the blood of patients that received a solid organ transplant or are suffering from an autoimmune disease.

Inversion of crosshole seismic data in heterogenous environments : Comparison of waveform and ray-based approaches

High-resolution tomographic imaging of the shallow subsurface is becoming increasingly important for a wide range of environmental, hydrological and engineering applications. Because of their superior resolution power, their sensitivity to pertinent petrophysical parameters, and their far reaching complementarities, both seismic and georadar crosshole imaging are of particular importance. To date, corresponding approaches have largely relied on asymptotic, ray-based approaches, which only account for a very small part of the observed wavefields, inherently suffer from a limited resolution, and in complex environments may prove to be inadequate. These problems can potentially be alleviated through waveform inversion. We have developed an acoustic waveform inversion approach for crosshole seismic data whose kernel is based on a finite-difference time-domain (FDTD) solution of the 2-D acoustic wave equations. This algorithm is tested on and applied to synthetic data from seismic velocity models of increasing complexity and realism and the results are compared to those obtained using state-of-the-art ray-based traveltime tomography. Regardless of the heterogeneity of the underlying models, the waveform inversion approach has the potential of reliably resolving both the geometry and the acoustic properties of features of the size of less than half a dominant wavelength. Our results do, however, also indicate that, within their inherent resolution limits, ray-based approaches provide an effective and efficient means to obtain satisfactory tomographic reconstructions of the seismic velocity structure in the presence of mild to moderate heterogeneity and in absence of strong scattering. Conversely, the excess effort of waveform inversion provides the greatest benefits for the most heterogeneous, and arguably most realistic, environments where multiple scattering effects tend to be prevalent and ray-based methods lose most of their effectiveness.

Estimation of inelastic seismic material properties of a surgical sea-bed from multi-component marine seismic data

To date, state-of-the-art seismic material parameter estimates from multi-component sea-bed seismic data are based on the assumption that the sea-bed consists of a fully elastic half-space. In reality, however, the shallow sea-bed generally consists of soft, unconsolidated sediments that are characterized by strong to very strong seismic attenuation. To explore the potential implications, we apply a state-of-the-art elastic decomposition algorithm to synthetic data for a range of canonical sea-bed models consisting of a viscoelastic half-space of varying attenuation. We find that in the presence of strong seismic attenuation, as quantified by Q-values of 10 or less, significant errors arise in the conventional elastic estimation of seismic properties. Tests on synthetic data indicate that these errors can be largely avoided by accounting for the inherent attenuation of the seafloor when estimating the seismic parameters. This can be achieved by replacing the real-valued expressions for the elastic moduli in the governing equations in the parameter estimation by their complex-valued viscoelastic equivalents. The practical application of our parameter procedure yields realistic estimates of the elastic seismic material properties of the shallow sea-bed, while the corresponding Q-estimates seem to be biased towards too low values, particularly for S-waves. Given that the estimation of inelastic material parameters is notoriously difficult, particularly in the immediate vicinity of the sea-bed, this is expected to be of interest and importance for civil and ocean engineering purposes.

Numerical analysis of wave-induced fluid flow effects on seismic data: Application to monitoring of CO2 storage at the Sleipner Field

In this work we analyze how patchy distributions of CO2 and brine within sand reservoirs may lead to significant attenuation and velocity dispersion effects, which in turn may have a profound impact on surface seismic data. The ultimate goal of this paper is to contribute to the understanding of these processes within the framework of the seismic monitoring of CO2 sequestration, a key strategy to mitigate global warming. We first carry out a Monte Carlo analysis to study the statistical behavior of attenuation and velocity dispersion of compressional waves traveling through rocks with properties similar to those at the Utsira Sand, Sleipner field, containing quasi-fractal patchy distributions of CO2 and brine. These results show that the mean patch size and CO2 saturation play key roles in the observed wave-induced fluid flow effects. The latter can be remarkably important when CO2 concentrations are low and mean patch sizes are relatively large. To analyze these effects on the corresponding surface seismic data, we perform numerical simulations of wave propagation considering reservoir models and CO2 accumulation patterns similar to the CO2 injection site in the Sleipner field. These numerical experiments suggest that wave-induced fluid flow effects may produce changes in the reservoir's seismic response, modifying significantly the main seismic attributes usually employed in the characterization of these environments. Consequently, the determination of the nature of the fluid distributions as well as the proper modeling of the seismic data constitute important aspects that should not be ignored in the seismic monitoring of CO2 sequestration problems.

Estimation of the correlation structure of crustal velocity heterogeneity from seismic reflection data

Numerous sources of evidence point to the fact that heterogeneity within the Earth's deep crystalline crust is complex and hence may be best described through stochastic rather than deterministic approaches. As seismic reflection imaging arguably offers the best means of sampling deep crustal rocks in situ, much interest has been expressed in using such data to characterize the stochastic nature of crustal heterogeneity. Previous work on this problem has shown that the spatial statistics of seismic reflection data are indeed related to those of the underlying heterogeneous seismic velocity distribution. As of yet, however, the nature of this relationship has remained elusive due to the fact that most of the work was either strictly empirical or based on incorrect methodological approaches. Here, we introduce a conceptual model, based on the assumption of weak scattering, that allows us to quantitatively link the second-order statistics of a 2-D seismic velocity distribution with those of the corresponding processed and depth-migrated seismic reflection image. We then perform a sensitivity study in order to investigate what information regarding the stochastic model parameters describing crustal velocity heterogeneity might potentially be recovered from the statistics of a seismic reflection image using this model. Finally, we present a Monte Carlo inversion strategy to estimate these parameters and we show examples of its application at two different source frequencies and using two different sets of prior information. Our results indicate that the inverse problem is inherently non-unique and that many different combinations of the vertical and lateral correlation lengths describing the velocity heterogeneity can yield seismic images with the same 2-D autocorrelation structure. The ratio of all of these possible combinations of vertical and lateral correlation lengths, however, remains roughly constant which indicates that, without additional prior information, the aspect ratio is the only parameter describing the stochastic seismic velocity structure that can be reliably recovered.

Joint inversion of crosshole GPR and seismic traveltime data

Joint inversion of crosshole ground-penetrating radar and seismic data can improve model resolution and fidelity of the resultant individual models. Model coupling obtained by minimizing or penalizing some measure of structural dissimilarity between models appears to be the most versatile approach because only weak assumptions about petrophysical relationships are required. Nevertheless, experimental results and petrophysical arguments suggest that when porosity variations are weak in saturated unconsolidated environments, then radar wave speed is approximately linearly related to seismic wave speed. Under such circumstances, model coupling also can be achieved by incorporating cross-covariances in the model regularization. In two case studies, structural similarity is imposed by penalizing models for which the model cross-gradients are nonzero. A first case study demonstrates improvements in model resolution by comparing the resulting models with borehole information, whereas a second case study uses point-spread functions. Although radar seismic wavespeed crossplots are very similar for the two case studies, the models plot in different portions of the graph, suggesting variances in porosity. Both examples display a close, quasilinear relationship between radar seismic wave speed in unconsolidated environments that is described rather well by the corresponding lower Hashin-Shtrikman (HS) bounds. Combining crossplots of the joint inversion models with HS bounds can constrain porosity and pore structure better than individual inversion results can.

Enhancing the vertical resolution of surface georadar data

There are far-reaching conceptual similarities between bi-static surface georadar and post-stack, "zero-offset" seismic reflection data, which is expressed in largely identical processing flows. One important difference is, however, that standard deconvolution algorithms routinely used to enhance the vertical resolution of seismic data are notoriously problematic or even detrimental to the overall signal quality when applied to surface georadar data. We have explored various options for alleviating this problem and have tested them on a geologically well-constrained surface georadar dataset. Standard stochastic and direct deterministic deconvolution approaches proved to be largely unsatisfactory. While least-squares-type deterministic deconvolution showed some promise, the inherent uncertainties involved in estimating the source wavelet introduced some artificial "ringiness". In contrast, we found spectral balancing approaches to be effective, practical and robust means for enhancing the vertical resolution of surface georadar data, particularly, but not exclusively, in the uppermost part of the georadar section, which is notoriously plagued by the interference of the direct air- and groundwaves. For the data considered in this study, it can be argued that band-limited spectral blueing may provide somewhat better results than standard band-limited spectral whitening, particularly in the uppermost part of the section affected by the interference of the air- and groundwaves. Interestingly, this finding is consistent with the fact that the amplitude spectrum resulting from least-squares-type deterministic deconvolution is characterized by a systematic enhancement of higher frequencies at the expense of lower frequencies and hence is blue rather than white. It is also consistent with increasing evidence that spectral "blueness" is a seemingly universal, albeit enigmatic, property of the distribution of reflection coefficients in the Earth. Our results therefore indicate that spectral balancing techniques in general and spectral blueing in particular represent simple, yet effective means of enhancing the vertical resolution of surface georadar data and, in many cases, could turn out to be a preferable alternative to standard deconvolution approaches.

Segmentation and grid generation for numerical simulations of vad connections with patient-specific data

Objectives: We are interested in the numerical simulation of the anastomotic region comprised between outflow canula of LVAD and the aorta. Segmenta¬tion, geometry reconstruction and grid generation from patient-specific data remain an issue because of the variable quality of DICOM images, in particular CT-scan (e.g. metallic noise of the device, non-aortic contrast phase). We pro¬pose a general framework to overcome this problem and create suitable grids for numerical simulations.Methods: Preliminary treatment of images is performed by reducing the level window and enhancing the contrast of the greyscale image using contrast-limited adaptive histogram equalization. A gradient anisotropic diffusion filter is applied to reduce the noise. Then, watershed segmentation algorithms and mathematical morphology filters allow reconstructing the patient geometry. This is done using the InsightToolKit library (www.itk.org). Finally the Vascular Model¬ing ToolKit (www.vmtk.org) and gmsh (www.geuz.org/gmsh) are used to create the meshes for the fluid (blood) and structure (arterial wall, outflow canula) and to a priori identify the boundary layers. The method is tested on five different patients with left ventricular assistance and who underwent a CT-scan exam.Results: This method produced good results in four patients. The anastomosis area is recovered and the generated grids are suitable for numerical simulations. In one patient the method failed to produce a good segmentation because of the small dimension of the aortic arch with respect to the image resolution.Conclusions: The described framework allows the use of data that could not be otherwise segmented by standard automatic segmentation tools. In particular the computational grids that have been generated are suitable for simulations that take into account fluid-structure interactions. Finally the presented method features a good reproducibility and fast application.

Source wavelet estimation for waveform inversion of crosshole georadar data

AbstractFor a wide range of environmental, hydrological, and engineering applications there is a fast growing need for high-resolution imaging. In this context, waveform tomographic imaging of crosshole georadar data is a powerful method able to provide images of pertinent electrical properties in near-surface environments with unprecedented spatial resolution. In contrast, conventional ray-based tomographic methods, which consider only a very limited part of the recorded signal (first-arrival traveltimes and maximum first-cycle amplitudes), suffer from inherent limitations in resolution and may prove to be inadequate in complex environments. For a typical crosshole georadar survey the potential improvement in resolution when using waveform-based approaches instead of ray-based approaches is in the range of one order-of- magnitude. Moreover, the spatial resolution of waveform-based inversions is comparable to that of common logging methods. While in exploration seismology waveform tomographic imaging has become well established over the past two decades, it is comparably still underdeveloped in the georadar domain despite corresponding needs. Recently, different groups have presented finite-difference time-domain waveform inversion schemes for crosshole georadar data, which are adaptations and extensions of Tarantola's seminal nonlinear generalized least-squares approach developed for the seismic case. First applications of these new crosshole georadar waveform inversion schemes on synthetic and field data have shown promising results. However, there is little known about the limits and performance of such schemes in complex environments. To this end, the general motivation of my thesis is the evaluation of the robustness and limitations of waveform inversion algorithms for crosshole georadar data in order to apply such schemes to a wide range of real world problems.One crucial issue to making applicable and effective any waveform scheme to real-world crosshole georadar problems is the accurate estimation of the source wavelet, which is unknown in reality. Waveform inversion schemes for crosshole georadar data require forward simulations of the wavefield in order to iteratively solve the inverse problem. Therefore, accurate knowledge of the source wavelet is critically important for successful application of such schemes. Relatively small differences in the estimated source wavelet shape can lead to large differences in the resulting tomograms. In the first part of my thesis, I explore the viability and robustness of a relatively simple iterative deconvolution technique that incorporates the estimation of the source wavelet into the waveform inversion procedure rather than adding additional model parameters into the inversion problem. Extensive tests indicate that this source wavelet estimation technique is simple yet effective, and is able to provide remarkably accurate and robust estimates of the source wavelet in the presence of strong heterogeneity in both the dielectric permittivity and electrical conductivity as well as significant ambient noise in the recorded data. Furthermore, our tests also indicate that the approach is insensitive to the phase characteristics of the starting wavelet, which is not the case when directly incorporating the wavelet estimation into the inverse problem.Another critical issue with crosshole georadar waveform inversion schemes which clearly needs to be investigated is the consequence of the common assumption of frequency- independent electromagnetic constitutive parameters. This is crucial since in reality, these parameters are known to be frequency-dependent and complex and thus recorded georadar data may show significant dispersive behaviour. In particular, in the presence of water, there is a wide body of evidence showing that the dielectric permittivity can be significantly frequency dependent over the GPR frequency range, due to a variety of relaxation processes. The second part of my thesis is therefore dedicated to the evaluation of the reconstruction limits of a non-dispersive crosshole georadar waveform inversion scheme in the presence of varying degrees of dielectric dispersion. I show that the inversion algorithm, combined with the iterative deconvolution-based source wavelet estimation procedure that is partially able to account for the frequency-dependent effects through an "effective" wavelet, performs remarkably well in weakly to moderately dispersive environments and has the ability to provide adequate tomographic reconstructions.

The Exotic Palaeo-Thethys Terrane in Central Iran: New geological data from Anarak, Jandaq and Posht-e-Badam areas

Résumé Le « terrane » d'Anarak-Jandak occupe une position géologique clé au nord-ouest du Microcontinent Centre-East Iranien (CE1M), connecté avec le Bloc du Grand Kavir et la ceinture métamorphique de Sanandaj-Sirjan. Nous discutons ici l'origine de ces différentes unités, reliées jusqu'à présent à des épisodes orogéniques d'âge Précambrien à Paléozoïque inférieur, pour conclure finalement de leur affinité paléotéthysienne. Leur histoire commence par un épisode de rifting d'âge Ordovicien supérieur-Dévonien inférieur, pour se terminer au Trias par la collision des blocs Cimmériens dérivé du Gondwana avec le Bloc du Turan d'affinité asiatique (événement Eocimmérien). La plus importante unité métamorphique affleurant au sud-ouest de la région de Jandak-Anarak-Kaboudan est une épaisse séquence silicoclastique à grains fins contenant des blocs ophiolitiques (marginal-sea-type), et des associations basalte-gabbro à signatures géochimiques de type supra-subduction. Dans la région de Nakhlak, nous avons daté ces gabbros par la méthode U-Pb à 387f0.11 Ma ; les roches métamorphiques pélitiques ont donné des âges de refroidissement Ar-Ar pour la muscovite de 320 à 333 Ma. Ce complexe d'accrétion "varisque" a été métamorphisé dans le faciès schiste vert-amphibolite au cours de l'accrétion de la ceinture granitique d'Airekan, d'âge Cambrien inférieur (549±15 Ma par la méthode U/Pb), qui affleure aujourd'hui à l'extrémité nord-ouest du terrane d'Anarak-Jandak . La subduction vers le nord de l'océan Paléotéthys depuis le Paléazoïque supérieur jusqu'au Trias, a permis l'accumulation de grandes quantités de matériel océanique dans la zone de subduction. Par exemple, une succession de guyots (Anarak, Kaboudan, et Meraji Seamounts) et de hauts sous-marins, entrés en collision oblique avec le prisme d'accrétion, est à l'origine d'un léger métamorphisme de type HP qui affecte ces séries {âges Ar-Ar de 280 à 230 Ma). De plus, le magmatisme bimodal de Chah Gorbeh est caractérisé d'une part par des roches de type trondjémite-gabbros (262 Ma), d'autre part par des laves en coussin de type basaltes alcalins-rhyolites; ces roches magmatiques ont recoupé l'ophiolite d'Anarak lors de la mise en place de cette dernière dans la fosse interne de subduction. Quant au prisme d'accrétion de Doshakh, d'âge essentiellement Permien supérieur, i1 a été accrété le long de la marge continentale et métamorphisé dans le faciès schiste vert. La fermeture de la Paléotéthys s'enregistre finalement par la sédimentation dans le bassin d'avant pays du flysch de Bayazeh, d'âge probable Triasique. Le matériel issu de l'arc magmatique de la Paléotéthys est très bien préservé dans les dépôts infra-arc Dévonien supérieur-Carbonifère de Godar-e-Siah, ainsi que dans la succession d'avant-arc de Nakhlak. Pendant l'intervalle Paléozoïque supérieur-Trias, la région de Jandak a été soumise à un régime extensif de type bassin d'arrière-arc, dont un témoin pourrait être la ceinture ophiolitique d'Arusan, elle-même comparable aux écailles ophiolitiques d'Aghdarband au nord-est de l'Iran. Cet ensemble métamorphique est recoupé par des granites d'arc à collisionnel datés à 215±15 Ma. Dans la région de Yazd, témoin de la marge passive Cimmérienne, la sédimentation syn-rift Silurienne à Dévonienne inférieure a été interrompue pendant l'intervalle Trias moyen-Trias supérieur; il en a été de même pour les dépôts de plate-forme Paléozoïque supérieur. L'érosion, qui dans ce dernier cas a atteint le Permien, pourrait être liée au bombement flexural de la marge passive. La collision finale n'a pas induit de déformations trop importantes, et se caractérise par la mise en place de nappes sur la marge passive. Cet événement est scellé par des dépôts molassique du Lias. D'un point de vue régional, la zone s'étendant actuellement de la Mer Noire au Pamir a été soumise à six épisodes d'extension-compression du Jurassique inférieur (début du l'ouverture en position arrière-arc de la Néotéthys) à l'Eocène moyen. Par exemple, le terrane d'AnarakJandak, probablement situé entre le Kopeh Dagh et la plate-forme nord Afghane, s'est complètement détaché de sa patrie d'origine au début du Crétacé supérieur. Des preuves de cet événement se retrouvent dans les séries de plate-forme de Khur (préservation de séries syn-rift puis de marge passive). Les ophiolites de Nain et de Sabzevar sont de plus interprétée comme un témoin de l'existence de ce bassin d'arrière-arc. Dans l'intervalle Eocène-Oligocène, l'indentation par la plaque indienne de l'Eurasie a été contemporaine de la rotation horaire de fragments de l'ancien microcontinent Iranien et de la formation du CEIM. Cette rotation est responsable du transport du terrane d'Anarak-Jandak vers sa position actuelle en Iran Central, et de la dislocation de Terranes de moindre importance, comme le bloc de Posht-e Badam. Depuis le Miocène supérieur, et à la suite de la collision entre l'Arabie et l'Iran, le ternane d'Anarak-Jandak a subi des déformations liées à l'activité d'une zone de cisaillement dextre parallèle à la suture du Zagros, à l'arrière de l'arc magmatique d'Uromieh-Dokhtar. Résumé large public Le Microcontinent Centre-Est Iranien occupe une position géologique clé au centre de l'Iran. Les différentes unités qui le composent, reliées jusqu'à présent à des épisodes orogéniques d'âge Précambrien à Paléozoïque inférieur, sont maintenant rajeunies et liés à la fermeture de l'océean Paléotéthys. Leur histoire commence par un épisode de rifting d'âge Ordovicien supérieur à Dévonien inférieur, pour se terminer au Trias par la collision des- blocs Cimmériens, dérivés du Gondwana, avec le Bloc du Turan d'affinité asiatique. Dans la marge active asiatique de la Paléotéthys, nous avons daté les restes d'un océan marginal à 387±0.11 Ma. Ce complexe d'accrétion a été métamorphisé au cours de la réaccrétion de la ceinture granitique d'Airekan, d'âge Cambrien inférieur (549±15 Ma), qui affleure aujourd'hui à l'extrémité nord-ouest du « terrane » d'Anarak-Jandak correspondant à la plus grande partie de la région étudiée. Le matériel issu de l'arc magmatique de la Paléotéthys est très bien préservé et daté du Dévonien supérieur-Carbonifère. Pendant l'intervalle Paléozoïque supérieur-Trias, la région a été soumise à un régime extensif de type bassin d'arrière-arc, dont un témoin pourrait être la ceinture ophiolitique d'Arusan, comparable aux écailles ophiolitiques d'Aghdarband au nord-est de l'Iran. Cet ensemble métamorphique est recoupé par des granites datés à 215±15 Ma. La subduction vers le nord de l'océan Paléotéthys depuis le Paléozoïque supérieur jusqu'au Trias, a permis l'accumulation de grandes quantités de matériel océanique dans la zone de subduction. Par exemple, une succession de volcans sous-marins, entrés en collision avec le prisme d'accrétion, est à l'origine d'un léger métamorphisme de type HP qui affecte ces séries (280 à 230 Ma). Quant au prisme d'accrétion de Doshakh, d'âge essentiellement Permien supérieur, il a été mis en place le long de la marge continentale et métamorphisé dans le faciès schiste vert. La fermeture de la Paléotéthys s'enregistre finalement par la sédimentation dans le bassin d'avant pays du flysch de Bayazeh, d'âge Triasique. Dans la région de Yazd, on trouve les témoins de la marge passive Cimmérienne, la sédimentation syn-rift Silurienne à Dévonienne inférieure a été interrompue pendant l'intervalle Trias moyen-Trias supérieur, marqué par la flexuration de la marge passive lorsqu'elle rentra en collision avec la marge active asiatique. Cet événement est scellé par des dépôts molassique à charbon du Lias. Le «terrane» d'Anarak-Jandak, probablement situé à l'origine entre le Kopeh Dagh et la plate-forme nord Afghane, s'est complètement détaché de cette région au début du Crétacé supérieur lors de l'ouverture d'un bassin d'arrière-arc, engendré, cette fois, par la subduction de l'océan Néotéthys situé au sud des blocs cimmériens. Des preuves de cet événement se retrouvent dans les séries syn-rift, puis de marge passive de Khour. Les ophiolites de Nain et de Sabzevar sont interprétées comme un témoin de l'existence de ce bassin d'arrière-arc. Dans l'intervalle Eocène-Oligocène, l'indentation de l'Eurasie par la plaque indienne a été contemporaine de la rotation horaire de fragments de l'ancien microcontinent centre-Iranien. Cette rotation de près de 90° est responsable du transport du « terrane » d'Anarak-Jandak vers sa position actuelle. Abstract The Anarak-Jandaq terrane occupies a strategic geological situation at the north-western part of the Central-East Iranian Microcontinent (CEIM) and in connection with the Great Kavir Block and Sanandaj-Sirjan metamorphic belt. Our recent findings redefine the origin of these mentioned areas so far attributed to the Precambrian-Early Palaeozoic orogenic episodes, to be now directly related to the tectonic evolution of the Palaeo-Tethys Ocean, commenced by Late Ordovician-Early Devonian rifting events and terminated in the Triassic by the Eocimmerian tectonic event due to the collision of the Cimmerian blocks with the Asiatic Turan block. The most distributed metamorphic unit that is exposed from the south-west of Jandaq to the Anarak and Kaboudan areas is a thick and fine grain siliciclastic sequence accompanied by marginal-sea-basin ophiolitic blocks including basalt-gabbro association with supra-subduction-geochemical signature. These gabbros in the Nakhlak area were dated by U/Pb method at 387.6 ± 0.11 Ma and the metamorphic pelitic rocks yielded a range of 320 to 333 Ma muscovite-cooling ages based on 40Ar/39 Ar method. This "Variscan" accretionary complex was metamorphosed in greenschist-amphibolite facies during accretion to the Lower Cambrian Airekan granitic belt (549 ± 15 Ma by U/Pb method) that crops out at the northwestern edge of the Anarak-Jandaq terrane. Continued northward subduction of the Palaeo-Tethys Ocean during the entire Late Palaeozoic-Middle Triassic brought huge amount of oceanic material to the subduction zone. One chain of Carboniferous-Triassic oceanic rises and seamounts (the Anarak, Kaboudan, and Meraji Seamounts) obliquely collided with the accretionary wedge and created a mild HP metamorphic event (280-230 Ma based on 40Ar/39Ar results). Bimodal magmatism of the Chah Gorbeh area is characterized by a 262 Ma trondjemite-gabbro as well as pillow alkalibasalts-rhyolites which intruded the Anarak ophiolite when it was being emplaced within the inner-wall trench. The mainly Late Permian-Triassic Doshakh wedge was accreted along the continent and metamorphosed under lower greenschist facies and the probable Triassic Bayazeh flysch filled the foreland basin during the final closure. The Palaeo-Tethys magmatic arc products have been well preserved in the Late Devonian-Carboniferous Godar-e-Siah intra-arc deposits and the Triassic Nakhlak fore-arc succession. During the Late Palaeozoic-Triassic times, the Jandaq area has been affected by back-arc extension and probably the Arusan ophiolitic belt is the remnant of this narrow basin comparable to the Aqdarband ophiolitic remnant in north-east Iran. This metamorphic belt was intruded by 215 ± 15 Ma arc to collisional granites. In the passive margin of the Cimmerian block, on the Yazd region, the Silurian-Early Devonian syn-rift succession as well as the nearly continuous Upper Palaeozoic platform-type deposition was interrupted during the Middle to Late Triassic time, local erosion down to Devonian levels may be related to flexural bulge erosion. The collision event was not so strong to generate intensive deformation but was accompanied by some nappe thrusting onto the passive margin. It is finally unconformably covered by Liassic continental molassic deposits. Related to the onset of Neo-Tethyan back-arc opening in Early Jurassic to Mid-Eocene times, six periods of extensional-compressional events have differently influenced an elongated area, extending from the West Black Sea to Pamir. The Anarak-Jandaq terrane which was situated somewhere in this affected area, probably between the Kopeh Dagh and North Afghan platform, was completely detached from its source at the beginning of the Late Cretaceous