A pseudo-spectral method for simulating of poro-elastic seismic wave propagation in complex borehole environments

We present a novel numerical approach for the comprehensive, flexible, and accurate simulation of poro-elastic wave propagation in cylindrical coordinates. An important application of this method is the modeling of complex seismic wave phenomena in fluid-filled boreholes, which represents a major, and as of yet largely unresolved, computational problem in exploration geophysics. In view of this, we consider a numerical mesh consisting of three concentric domains representing the borehole fluid in the center, the borehole casing and the surrounding porous formation. The spatial discretization is based on a Chebyshev expansion in the radial direction, Fourier expansions in the other directions, and a Runge-Kutta integration scheme for the time evolution. A domain decomposition method based on the method of characteristics is used to match the boundary conditions at the fluid/porous-solid and porous-solid/porous-solid interfaces. The viability and accuracy of the proposed method has been tested and verified in 2D polar coordinates through comparisons with analytical solutions as well as with the results obtained with a corresponding, previously published, and independently benchmarked solution for 2D Cartesian coordinates. The proposed numerical solution also satisfies the reciprocity theorem, which indicates that the inherent singularity associated with the origin of the polar coordinate system is handled adequately.

Time-lapse and probabilistic inversion strategies for plane-wave electromagnetic data.

L'utilisation efficace des systèmes géothermaux, la séquestration du CO2 pour limiter le changement climatique et la prévention de l'intrusion d'eau salée dans les aquifères costaux ne sont que quelques exemples qui démontrent notre besoin en technologies nouvelles pour suivre l'évolution des processus souterrains à partir de la surface. Un défi majeur est d'assurer la caractérisation et l'optimisation des performances de ces technologies à différentes échelles spatiales et temporelles. Les méthodes électromagnétiques (EM) d'ondes planes sont sensibles à la conductivité électrique du sous-sol et, par conséquent, à la conductivité électrique des fluides saturant la roche, à la présence de fractures connectées, à la température et aux matériaux géologiques. Ces méthodes sont régies par des équations valides sur de larges gammes de fréquences, permettant détudier de manières analogues des processus allant de quelques mètres sous la surface jusqu'à plusieurs kilomètres de profondeur. Néanmoins, ces méthodes sont soumises à une perte de résolution avec la profondeur à cause des propriétés diffusives du champ électromagnétique. Pour cette raison, l'estimation des modèles du sous-sol par ces méthodes doit prendre en compte des informations a priori afin de contraindre les modèles autant que possible et de permettre la quantification des incertitudes de ces modèles de façon appropriée. Dans la présente thèse, je développe des approches permettant la caractérisation statique et dynamique du sous-sol à l'aide d'ondes EM planes. Dans une première partie, je présente une approche déterministe permettant de réaliser des inversions répétées dans le temps (time-lapse) de données d'ondes EM planes en deux dimensions. Cette stratégie est basée sur l'incorporation dans l'algorithme d'informations a priori en fonction des changements du modèle de conductivité électrique attendus. Ceci est réalisé en intégrant une régularisation stochastique et des contraintes flexibles par rapport à la gamme des changements attendus en utilisant les multiplicateurs de Lagrange. J'utilise des normes différentes de la norme l2 pour contraindre la structure du modèle et obtenir des transitions abruptes entre les régions du model qui subissent des changements dans le temps et celles qui n'en subissent pas. Aussi, j'incorpore une stratégie afin d'éliminer les erreurs systématiques de données time-lapse. Ce travail a mis en évidence l'amélioration de la caractérisation des changements temporels par rapport aux approches classiques qui réalisent des inversions indépendantes à chaque pas de temps et comparent les modèles. Dans la seconde partie de cette thèse, j'adopte un formalisme bayésien et je teste la possibilité de quantifier les incertitudes sur les paramètres du modèle dans l'inversion d'ondes EM planes. Pour ce faire, je présente une stratégie d'inversion probabiliste basée sur des pixels à deux dimensions pour des inversions de données d'ondes EM planes et de tomographies de résistivité électrique (ERT) séparées et jointes. Je compare les incertitudes des paramètres du modèle en considérant différents types d'information a priori sur la structure du modèle et différentes fonctions de vraisemblance pour décrire les erreurs sur les données. Les résultats indiquent que la régularisation du modèle est nécessaire lorsqu'on a à faire à un large nombre de paramètres car cela permet d'accélérer la convergence des chaînes et d'obtenir des modèles plus réalistes. Cependent, ces contraintes mènent à des incertitudes d'estimations plus faibles, ce qui implique des distributions a posteriori qui ne contiennent pas le vrai modèledans les régions ou` la méthode présente une sensibilité limitée. Cette situation peut être améliorée en combinant des méthodes d'ondes EM planes avec d'autres méthodes complémentaires telles que l'ERT. De plus, je montre que le poids de régularisation des paramètres et l'écart-type des erreurs sur les données peuvent être retrouvés par une inversion probabiliste. Finalement, j'évalue la possibilité de caractériser une distribution tridimensionnelle d'un panache de traceur salin injecté dans le sous-sol en réalisant une inversion probabiliste time-lapse tridimensionnelle d'ondes EM planes. Etant donné que les inversions probabilistes sont très coûteuses en temps de calcul lorsque l'espace des paramètres présente une grande dimension, je propose une stratégie de réduction du modèle ou` les coefficients de décomposition des moments de Legendre du panache de traceur injecté ainsi que sa position sont estimés. Pour ce faire, un modèle de résistivité de base est nécessaire. Il peut être obtenu avant l'expérience time-lapse. Un test synthétique montre que la méthodologie marche bien quand le modèle de résistivité de base est caractérisé correctement. Cette méthodologie est aussi appliquée à un test de trac¸age par injection d'une solution saline et d'acides réalisé dans un système géothermal en Australie, puis comparée à une inversion time-lapse tridimensionnelle réalisée selon une approche déterministe. L'inversion probabiliste permet de mieux contraindre le panache du traceur salin gr^ace à la grande quantité d'informations a priori incluse dans l'algorithme. Néanmoins, les changements de conductivités nécessaires pour expliquer les changements observés dans les données sont plus grands que ce qu'expliquent notre connaissance actuelle des phénomenès physiques. Ce problème peut être lié à la qualité limitée du modèle de résistivité de base utilisé, indiquant ainsi que des efforts plus grands devront être fournis dans le futur pour obtenir des modèles de base de bonne qualité avant de réaliser des expériences dynamiques. Les études décrites dans cette thèse montrent que les méthodes d'ondes EM planes sont très utiles pour caractériser et suivre les variations temporelles du sous-sol sur de larges échelles. Les présentes approches améliorent l'évaluation des modèles obtenus, autant en termes d'incorporation d'informations a priori, qu'en termes de quantification d'incertitudes a posteriori. De plus, les stratégies développées peuvent être appliquées à d'autres méthodes géophysiques, et offrent une grande flexibilité pour l'incorporation d'informations additionnelles lorsqu'elles sont disponibles. -- The efficient use of geothermal systems, the sequestration of CO2 to mitigate climate change, and the prevention of seawater intrusion in coastal aquifers are only some examples that demonstrate the need for novel technologies to monitor subsurface processes from the surface. A main challenge is to assure optimal performance of such technologies at different temporal and spatial scales. Plane-wave electromagnetic (EM) methods are sensitive to subsurface electrical conductivity and consequently to fluid conductivity, fracture connectivity, temperature, and rock mineralogy. These methods have governing equations that are the same over a large range of frequencies, thus allowing to study in an analogous manner processes on scales ranging from few meters close to the surface down to several hundreds of kilometers depth. Unfortunately, they suffer from a significant resolution loss with depth due to the diffusive nature of the electromagnetic fields. Therefore, estimations of subsurface models that use these methods should incorporate a priori information to better constrain the models, and provide appropriate measures of model uncertainty. During my thesis, I have developed approaches to improve the static and dynamic characterization of the subsurface with plane-wave EM methods. In the first part of this thesis, I present a two-dimensional deterministic approach to perform time-lapse inversion of plane-wave EM data. The strategy is based on the incorporation of prior information into the inversion algorithm regarding the expected temporal changes in electrical conductivity. This is done by incorporating a flexible stochastic regularization and constraints regarding the expected ranges of the changes by using Lagrange multipliers. I use non-l2 norms to penalize the model update in order to obtain sharp transitions between regions that experience temporal changes and regions that do not. I also incorporate a time-lapse differencing strategy to remove systematic errors in the time-lapse inversion. This work presents improvements in the characterization of temporal changes with respect to the classical approach of performing separate inversions and computing differences between the models. In the second part of this thesis, I adopt a Bayesian framework and use Markov chain Monte Carlo (MCMC) simulations to quantify model parameter uncertainty in plane-wave EM inversion. For this purpose, I present a two-dimensional pixel-based probabilistic inversion strategy for separate and joint inversions of plane-wave EM and electrical resistivity tomography (ERT) data. I compare the uncertainties of the model parameters when considering different types of prior information on the model structure and different likelihood functions to describe the data errors. The results indicate that model regularization is necessary when dealing with a large number of model parameters because it helps to accelerate the convergence of the chains and leads to more realistic models. These constraints also lead to smaller uncertainty estimates, which imply posterior distributions that do not include the true underlying model in regions where the method has limited sensitivity. This situation can be improved by combining planewave EM methods with complimentary geophysical methods such as ERT. In addition, I show that an appropriate regularization weight and the standard deviation of the data errors can be retrieved by the MCMC inversion. Finally, I evaluate the possibility of characterizing the three-dimensional distribution of an injected water plume by performing three-dimensional time-lapse MCMC inversion of planewave EM data. Since MCMC inversion involves a significant computational burden in high parameter dimensions, I propose a model reduction strategy where the coefficients of a Legendre moment decomposition of the injected water plume and its location are estimated. For this purpose, a base resistivity model is needed which is obtained prior to the time-lapse experiment. A synthetic test shows that the methodology works well when the base resistivity model is correctly characterized. The methodology is also applied to an injection experiment performed in a geothermal system in Australia, and compared to a three-dimensional time-lapse inversion performed within a deterministic framework. The MCMC inversion better constrains the water plumes due to the larger amount of prior information that is included in the algorithm. The conductivity changes needed to explain the time-lapse data are much larger than what is physically possible based on present day understandings. This issue may be related to the base resistivity model used, therefore indicating that more efforts should be given to obtain high-quality base models prior to dynamic experiments. The studies described herein give clear evidence that plane-wave EM methods are useful to characterize and monitor the subsurface at a wide range of scales. The presented approaches contribute to an improved appraisal of the obtained models, both in terms of the incorporation of prior information in the algorithms and the posterior uncertainty quantification. In addition, the developed strategies can be applied to other geophysical methods, and offer great flexibility to incorporate additional information when available.

Diagnostic value of MRI for the detection of suspected placental invasion: Does it depend on observers' experience?

Purpose: To evaluate the diagnostic value of specific MR features for detection of suspected placental invasion according to observers' experience.Methods and Materials: Our study population included 25 pregnant women (mean age 35.16) investigated by prenatal MRI. In twelve out of them placental invasion was histopathologically proven, the 13 other women (52%) without placental invasion served as control group. Multiplanar T1- and T2-weighted sequences had been performed mostly without IV contrast injection (1.5 T). MR examinations of the two groups were rendered anonymous, mixed, then independently and retrospectively reviewed by two senior and two junior radiologists in view of 8 MR features indicating placentar invasion including the degree. Results were compared with surgical diagnosis (placenta normal/increta/accreta/percreta). Interobserver agrement between senior and junior readers were calculated. Stepwise logistic regression and receiver operating (ROC) curvers were performed.Results: Demographics between the two groups were not statistically different. Overall sensitivity and specificity for detecting placentar invasion was 90.9% and 75.0% for senior readers, and 81.8% and 61.8% for junior readers respectively. The most significant MR features indicating placentar invasion were T2 hypointense placental bands, followed by placenta praevia, focally interrupted myometrial border, posterior placental insertion, and heterogeneous placental signal. For each of the evaluated MR features the interobserver agreement kappa between the two senior readers was superior than that between the junior readers, ranging from bad (<0.4) to good (0.4-0.75).Conclusions: MRI can be a reliable and reproducible tool for detection of suspected placentar invasion, however very variable according to the observers' experience.

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.

The role of MMP-9 and its fibronectin-like domain in tumor growth & invasion : certainties, controversies & discrepancies

Tumors are often compared to wounds that do not heal, where the crosstalk between tumor cells and their surrounding stroma is crucial at all stages of development, from the initial primary growth to metastasis. Similar to wound healing, fibroblasts in the tumor stroma differentiate into myofibroblasts, also referred to as "cancer-associated fibroblasts" (CAFs), primarily, but not exclusively, in response to transforming growth factor-ß (TGF-ß). Myofibroblasts in turn enhance tumor progression by remodeling the stroma. Among molecules implicated in stroma remodeling, matrix metalloproteinases (MMPs), and MMP-g in particular, play a prominent role. However, the mechanisms that regulate MMP-g activation and function remain poorly understood. Recent evidence indicates that tumor cell surface association of MMP-g is an important event in its activation, and more generally in tumor growth and invasion. In the present work we address the potential association of MMP-g activity with cell-surface recruitment to human fibroblasts. We show for the first time that recruitment of MMP-g to the MRC-5 fibroblast cell surface occurs through the fibronectin-like (FN) domain, shared only by MMP-g and MMP-2 among all the MMPs. Functional assays suggest that both the pro- and active form of MMP-g trigger a-smooth muscle actin (aSMA) expression in resting fibroblasts that reflects myofibroblast differentiation, possibly through TGF-ß activation. Moreover, the FN domain of MMP-g inhibits both MMP-g-induced TGF-ß activation and aSMA expression by sequestering MMP-g. Xenograft experiments in NOD/SCID mice using HT1080 fibrosarcoma or MDA-MD231 breast adenocarcinoma cells stably expressing the FN domain of MMP-g revealed no changes in primary tumor growth. However, in the context of metastasis, expression of the FN domain by these same tumor cells dramatically increased their metastatic proclivity whereas expression of wt MMP-g either promoted no change or actually reduced the number of metastases. We observed a decrease of an active form of MMP-g in MDA-MB231 cells overexpressing the FN domain suggesting that the FN domain may inhibit MMP-g activity in Tumors are often compared to wounds that do not heal, where the crosstalk between tumor cells and their surrounding stroma is crucial at all stages of development, from the initial primary growth to metastasis. Similar to wound healing, fibroblasts in the tumor stroma differentiate into myofibroblasts, also referred to as "cancer-associated fibroblasts" (CAFs), primarily, but not exclusively, in response to transforming growth factor-ß (TGF-ß). Myofibroblasts in turn enhance tumor progression by remodeling the stroma. Among molecules implicated in stroma remodeling, matrix metalloproteinases (MMPs), and MMP-g in particular, play a prominent role. However, the mechanisms that regulate MMP-g activation and function remain poorly understood. Recent evidence indicates that tumor cell surface association of MMP-g is an important event in its activation, and more generally in tumor growth and invasion. In the present work we address the potential association of MMP-g activity with cell-surface recruitment to human fibroblasts. We show for the first time that recruitment of MMP-g to the MRC-5 fibroblast cell surface occurs through the fibronectin-like (FN) domain, shared only by MMP-g and MMP-2 among all the MMPs. Functional assays suggest that both the pro- and active form of MMP-g trigger a-smooth muscle actin (aSMA) expression in resting fibroblasts that reflects myofibroblast differentiation, possibly through TGF-ß activation. Moreover, the FN domain of MMP-g inhibits both MMP-g-induced TGF-ß activation and aSMA expression by sequestering MMP-g. Xenograft experiments in NOD/SCID mice using HT1080 fibrosarcoma or MDA-MD231 breast adenocarcinoma cells stably expressing the FN domain of MMP-9 revealed no changes in primary tumor growth. However, in the context of metastasis, expression of the FN domain by these same tumor cells dramatically increased their metastatic proclivity whereas expression of wt MMP-g either promoted no change or actually reduced the number of metastases. We observed a decrease of an active form of MMP-9 in MDA-MB231 cells overexpressing the FN domain suggesting that the FN domain may inhibit MMP-9 activity in those cells and therefore prevent MMP-9-induced activation of TGF-b, which results in increased invasion. Curiously, xenografts of SW480 colorectal adenocarcinoma cells stably expressing the FN domain of MMP-9 displayed reduced growth at both the primary (subcutaneous) injection site and the lungs of NOD/SCID mice, in experimental metastasis assays, whilst the same cells overexpressing wt MMP-9 showed enhanced growth and dissemination. Gelatin zymography of conditioned medium revealed that these effects may be due to the FN domain, which displaces MMP-9 from SW480 cell surface. These observations suggest a dual role of MMP-9 and its FN domain in primary tumor growth and metastasis, underscoring the notion that the effect of MMP-9 on tumor cells may depend on the cell type and highlighting possible protective effects of MMPs in tumor progression.

Accurate calculations of intermolecular interaction energies using explicitly correlated coupled cluster wave functions and a dispersion-weighted MP2 method.

Explicitly correlated coupled-cluster calculations of intermolecular interaction energies for the S22 benchmark set of Jurecka, Sponer, Cerny, and Hobza (Chem. Phys. Phys. Chem. 2006, 8, 1985) are presented. Results obtained with the recently proposed CCSD(T)-F12a method and augmented double-zeta basis sets are found to be in very close agreement with basis set extrapolated conventional CCSD(T) results. Furthermore, we propose a dispersion-weighted MP2 (DW-MP2) approximation that combines the good accuracy of MP2 for complexes with predominately electrostatic bonding and SCS-MP2 for dispersion-dominated ones. The MP2-F12 and SCS-MP2-F12 correlation energies are weighted by a switching function that depends on the relative HF and correlation contributions to the interaction energy. For the S22 set, this yields a mean absolute deviation of 0.2 kcal/mol from the CCSD(T)-F12a results. The method, which allows obtaining accurate results at low cost, is also tested for a number of dimers that are not in the training set.

Invader immunology: invasion history alters immune system function in cane toads (Rhinella marina) in tropical Australia.

Because an individual's investment into the immune system may modify its dispersal rate, immune function may evolve rapidly in an invader. We collected cane toads (Rhinella marina) from sites spanning their 75-year invasion history in Australia, bred them, and raised their progeny in standard conditions. Evolved shifts in immune function should manifest as differences in immune responses among the progeny of parents collected in different locations. Parental location did not affect the offspring's cell-mediated immune response or stress response, but blood from the offspring of invasion-front toads had more neutrophils, and was more effective at phagocytosis and killing bacteria. These latter measures of immune function are negatively correlated with rate of dispersal in free-ranging toads. Our results suggest that the invasion of tropical Australia by cane toads has resulted in rapid genetically based compensatory shifts in the aspects of immune responses that are most compromised by the rigours of long-distance dispersal.

The evolutionary host switches of Polychromophilus: a multi-gene phylogeny of the bat malaria genus suggests a second invasion of mammals by a haemosporidian parasite.

BACKGROUND: The majority of Haemosporida species infect birds or reptiles, but many important genera, including Plasmodium, infect mammals. Dipteran vectors shared by avian, reptilian and mammalian Haemosporida, suggest multiple invasions of Mammalia during haemosporidian evolution; yet, phylogenetic analyses have detected only a single invasion event. Until now, several important mammal-infecting genera have been absent in these analyses. This study focuses on the evolutionary origin of Polychromophilus, a unique malaria genus that only infects bats (Microchiroptera) and is transmitted by bat flies (Nycteribiidae). METHODS: Two species of Polychromophilus were obtained from wild bats caught in Switzerland. These were molecularly characterized using four genes (asl, clpc, coI, cytb) from the three different genomes (nucleus, apicoplast, mitochondrion). These data were then combined with data of 60 taxa of Haemosporida available in GenBank. Bayesian inference, maximum likelihood and a range of rooting methods were used to test specific hypotheses concerning the phylogenetic relationships between Polychromophilus and the other haemosporidian genera. RESULTS: The Polychromophilus melanipherus and Polychromophilus murinus samples show genetically distinct patterns and group according to species. The Bayesian tree topology suggests that the monophyletic clade of Polychromophilus falls within the avian/saurian clade of Plasmodium and directed hypothesis testing confirms the Plasmodium origin. CONCLUSION: Polychromophilus' ancestor was most likely a bird- or reptile-infecting Plasmodium before it switched to bats. The invasion of mammals as hosts has, therefore, not been a unique event in the evolutionary history of Haemosporida, despite the suspected costs of adapting to a new host. This was, moreover, accompanied by a switch in dipteran host.

Invasion of lesion territory by regenerating fibers after spinal cord injury in adult macaque monkeys.

In adult macaque monkeys subjected to an incomplete spinal cord injury (SCI), corticospinal (CS) fibers are rarely observed to grow in the lesion territory. This situation is little affected by the application of an anti-Nogo-A antibody which otherwise fosters the growth of CS fibers rostrally and caudally to the lesion. However, when using the Sternberger monoclonal-incorporated antibody 32 (SMI-32), a marker detecting a non-phosphorylated neurofilament epitope, numerous SMI-32-positive (+) fibers were observed in the spinal lesion territory of 18 adult macaque monkeys; eight of these animals had received a control antibody infusion intrathecally for 1month after the injury, five animals an anti-Nogo-A antibody, and five animals received an anti-Nogo-A antibody together with brain-derived neurotrophic factor (BDNF). These fibers occupied the whole dorso-ventral axis of the lesion site with a tendency to accumulate on the ventral side, and their trajectories were erratic. Most of these fibers (about 87%) were larger than 1.3μm and densely SMI-32 (+) stained. In the undamaged spinal tissue, motoneurons form the only large population of SMI-32 (+) neurons which are densely stained and have large diameter axons. These data therefore suggest that a sizeable proportion of the fibers seen in the lesion territory originate from motoneurons, although fibers of other origins could also contribute. Neither the presence of the antibody neutralizing Nogo-A alone, nor the presence of the antibody neutralizing Nogo-A combined with BDNF influenced the number or the length of the SMI-32 (+) fibers in the spinal lesion area. In summary, our data show that after a spinal cord lesion in adult monkeys, the lesion site is colonized by fibers, a large portion of which presumably originate from motoneurons.

Evidence of climatic niche shift during biological invasion

Niche-based models calibrated in the native range by relating species observations to climatic variables are commonly used to predict the potential spatial extent of species' invasion. This climate matching approach relies on the assumption that invasive species conserve their climatic niche in the invaded ranges. We test this assumption by analysing the climatic niche spaces of Spotted Knapweed in western North America and Europe. We show with robust cross-continental data that a shift of the observed climatic niche occurred between native and non-native ranges, providing the first empirical evidence that an invasive species can occupy climatically distinct niche spaces following its introduction into a new area. The models fail to predict the current invaded distribution, but correctly predict areas of introduction. Climate matching is thus a useful approach to identify areas at risk of introduction and establishment of newly or not-yet-introduced neophytes, but may not predict the full extent of invasions.

Wave-vector dependence of spin and density multipole excitations in quantum dots

We have employed time-dependent local-spin density-functional theory to analyze the multipole spin and charge density excitations in GaAs-AlxGa1-xAs quantum dots. The on-plane transferred momentum degree of freedom has been taken into account, and the wave-vector dependence of the excitations is discussed. In agreement with previous experiments, we have found that the energies of these modes do not depend on the transferred wave vector, although their intensities do. Comparison with a recent resonant Raman scattering experiment [C. Schüller et al., Phys. Rev. Lett. 80, 2673 (1998)] is made. This allows us to identify the angular momentum of several of the observed modes as well as to reproduce their energies