PPARγ Agonist-Induced Fluid Retention Depends on αENaC Expression in Connecting Tubules.

BACKGROUND/AIMS: Thiazolidinediones (TZDs, like rosiglitazone (RGZ)) are peroxisome proliferator-activated receptor γ (PPARγ) agonists used to treat type 2 diabetes. Clinical limitations include TZD-induced fluid retention and body weight (BW) increase, which are inhibited by amiloride, an epithelial-sodium channel (ENaC) blocker. RGZ-induced fluid retention is maintained in mice with αENaC knockdown in the collecting duct (CD). Since ENaC in the connecting tubule (CNT) rather than in CD appears to be critical for normal NaCl retention, we aimed to further explore the role of ENaC in CNT in RGZ-induced fluid retention. METHODS: Mice with conditional inactivation of αENaC in both CNT and CD were used (αENaC lox/lox AQP2-Cre; 'αENaC-CNT/CD-KO') and compared with littermate controls (αENaC lox/lox mice; 'WT'). BW was monitored and total body water (TBW) and extracellular fluid volume (ECF) were determined by bioelectrical impedance spectroscopy (BIS) before and after RGZ (320 mg/kg diet for 10 days). RESULTS: On regular NaCl diet, αENaC-CNT/CD-KO had normal BW, TBW, ECF, hematocrit, and plasma Na(+), K(+), and creatinine, associated with an increase in plasma aldosterone compared with WT. Challenging αENaC-CNT/CD-KO with a low NaCl diet unmasked impaired NaCl and K homeostasis, consistent with effective knockdown of αENaC. In WT, RGZ increased BW (+6.1%), TBW (+8.4%) and ECF (+10%), consistent with fluid retention. These changes were significantly attenuated in αENaC-CNT/CD-KO (+3.4, 1.3, and 4.3%). CONCLUSION: Together with the previous studies, the current results are consistent with a role of αENaC in CNT in RGZ-induced fluid retention, which dovetails with the physiological relevance of ENaC in this segment. © 2014 S. Karger AG, Basel.

Creatine and guanidinoacetate transport at blood-brain and blood-cerebrospinal fluid barriers.

While it was thought that most of cerebral creatine is of peripheral origin, AGAT and GAMT are well expressed in CNS where brain cells synthesize creatine. While the creatine transporter SLC6A8 is expressed by microcapillary endothelial cells (MCEC) at blood-brain barrier (BBB), it is absent from their surrounding astrocytes. This raised the concept that BBB has a limited permeability for peripheral creatine, and that the brain supplies a part of its creatine by endogenous synthesis. This review brings together the latest data on creatine and guanidinoacetate transport through BBB and blood-CSF barrier (BCSFB) with the clinical evidence of AGAT-, GAMT- and SLC6A8-deficient patients, in order to delineate a clearer view on the roles of BBB and BCSFB in the transport of creatine and guanidinoacetate between periphery and CNS, and on brain synthesis and transport of creatine. It shows that in physiological conditions, creatine is taken up by CNS from periphery through SLC6A8 at BBB, but in limited amounts, and that CNS also needs its own creatine synthesis. No uptake of guanidinoacetate from periphery occurs at BBB except under GAMT deficiency, but a net exit of guanidinoacetate seems to occur from CSF to blood at BCSFB, predominantly through the taurine transporter TauT.

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.

Metallogenic model of the Trepca Pb-Zn-Ag skarn deposit, Kosovo: Evidence from fluid inclusions, Rare Earth Elements, and stable isotope data

The Trepca Pb-Zn-Ag skarn deposit (29 Mt of ore at 3.45% Pb, 2.30% Zn, and 80 g/t Ag) is located in the Kopaonik block of the western Vardar zone, Kosovo. The mineralization, hosted by recrystallized limestone of Upper Triassic age, was structurally and lithologically controlled. Ore deposition is spatially and temporally related with the postcollisional magmatism of Oligocene age (23-26 Ma). The deposit was formed during two distinct mineralization stages: an early prograde closed-system and a later retrograde open-system stage. The prograde mineralization consisting mainly of pyroxenes (Hd(54-100)Jo(0-45)Di(0-45)) resulted from the interaction of magmatic fluids associated with Oligocene (23-26 Ma) postcollisional magmatism. Whereas there is no direct contact between magmatic rocks and the mineralization, the deposit is classified as a distal Pb-Zn-Ag skarn. Abundant pyroxene reflects low oxygen fugacity (<10(-31) bar) and anhydrous environment. Fluid inclusion data and mineral assemblage limit the prograde stage within a temperature range between 390 degrees and 475 degrees C. Formation pressure is estimated below 900 bars. Isotopic composition of aqueous fluid, inclusions hosted by hedenbergite (delta D = -108 to -130 parts per thousand; delta O-18 = 7.5-8.0 parts per thousand), Mn-enriched mineralogy and high REE content of the host carbonates at the contact with the skarn mineralization suggest that a magmatic fluid was modified during its infiltration through the country rocks. The retrograde mineral assemblage comprises ilvaite, magnetite, arsenopyrite, pyrrhotite, marcasite, pyrite, quartz, and various carbonates. Increases in oxygen and sulfur fugacities, as well as a hydrous character of mineralization, require an open-system model. The opening of the system is related to phreatomagmatic explosion and formation of the breccia. Arsenopyrite geothermometer limits the retrograde stage within the temperature range between 350 degrees and 380 degrees C and sulfur fugacity between 10(-8.8) and 10(-7.2) bars. The principal ore minerals, galena, sphalerite, pyrite, and minor chalcopyrite, were deposited from a moderately saline Ca-Na chloride fluid at around 350 degrees C. According to the isotopic composition of fluid inclusions hosted by sphalerite (delta D = -55 to -74 parts per thousand; delta O-18 = -9.6 to -13.6 parts per thousand), the fluid responsible for ore deposition was dominantly meteoric in origin. The delta S-31 values of the sulfides spanning between -5.5 and +10 parts per thousand point to a magmatic origin of sulfur. Ore deposition appears to have been largely contemporaneous with the retrograde stage of the skarn development. Postore stage accompanied the precipitation of significant amount of carbonates including the travertine deposits at the deposit surface. Mineralogical composition of travertine varies from calcite to siderite and all carbonates contain significant amounts of Mn. Decreased formation temperature and depletion in the REE content point to an influence of pH-neutralized cold ground water and dying magmatic system.

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.

Intrathecal analgesia and restrictive perioperative fluid management within enhanced recovery pathway: hemodynamic implications.

BACKGROUND: Intrathecal analgesia and avoidance of perioperative fluid overload are key items within enhanced recovery pathways. Potential side effects include hypotension and renal dysfunction. STUDY DESIGN: From January 2010 until May 2010, all patients undergoing colorectal surgery within enhanced recovery pathways were included in this retrospective cohort study and were analyzed by intrathecal analgesia (IT) vs none (noIT). Primary outcomes measures were systolic and diastolic blood pressure, mean arterial pressure, and heart rate for 48 hours after surgery. Renal function was assessed by urine output and creatinine values. RESULTS: One hundred and sixty-three consecutive colorectal patients (127 IT and 36 noIT) were included in the analysis. Both patient groups showed low blood pressure values within the first 4 to 12 hours and a steady increase thereafter before return to baseline values after about 24 hours. Systolic and diastolic blood pressure and mean arterial pressure were significantly lower until 16 hours after surgery in patients having IT compared with the noIT group. Low urine output (<0.5 mL/kg/h) was reported in 11% vs 29% (IT vs noIT; p = 0.010) intraoperatively, 20% vs 11% (p = 0.387), 33% vs 22% (p = 0.304), and 31% vs 21% (p = 0.478) for postanesthesia care unit and postoperative days 1 and 2, respectively. Only 3 of 127 (2.4%) IT and 1 of 36 (2.8%) noIT patients had a transitory creatinine increase >50%; no patients required dialysis. CONCLUSIONS: Postoperative hypotension affects approximately 10% of patients within an enhanced recovery pathway and is slightly more pronounced in patients with IT. Hemodynamic depression persists for <20 hours after surgery; it has no measurable negative impact and therefore cannot justify detrimental postoperative fluid overload.

Role of fluid mixing in deep dissolution of carbonates

The presence of cavities filled with new minerals in carbonate rocks is a common feature in oil reservoirs and lead-zinc deposits. Since groundwater equilibrates rapidly with carbonates, the presence of dissolution cavities in deep carbonate host rocks is a paradox. Two alternative geochemical processes have been proposed to dissolve carbonates at depth: hydrogen sulfide oxidation to sulfuric acid, and metal sulfide precipitation. With the aid of geochemical modeling we show that mixing two warm solutions saturated with carbonate results in a new solution that dissolves limestone. Variations in the proportion of the end-member fluids can also form a supersaturated mixture and fill the cavity with a new generation of carbonate. Mixing is in general more effective in dissolving carbonates than the aforementioned processes. Moreover, mixing is consistent with the wide set of textures and mineral proportions observed in cavity infillings.

Submacular fluid after encircling buckle surgery for inferior macula-off retinal detachment in young patients.

PURPOSE: Characterization of persistent diffuse subretinal fluid using optical coherence tomography (OCT) after successful encircling buckle surgery for inferior macula-off retinal detachment in young patients. METHODS: Institutional retrospective review of six young patients (mean age 31 +/- 6 years; five female, one male) with spontaneous inferior rhegmatogenous macula-off retinal detachment. All patients were treated with encircling buckle surgery and five out of six underwent additional external drainage of subretinal fluid. Mean follow-up was 37 +/- 25 months (range 17-75 months) and included complete ophthalmic and OCT examination. RESULTS: At 6 months, 100% of patients showed persistence of subretinal fluid on OCT. Four patients had diffuse fluid accumulation, whereas two patients showed a 'bleb-like' accumulation of fluid. This fluid was present independent of whether or not patients had been treated with external fluid drainage. Subretinal fluid only started to disappear on OCT between 6 and more than 12 months after surgery. CONCLUSION: Young patients with inferior macula-off retinal detachments and a marginally liquefied vitreous may show persisting postoperative subclinical fluid under the macula for longer periods of time than described previously.

Local Renal Circadian Clocks Control Fluid-Electrolyte Homeostasis and BP.

The circadian timing system is critically involved in the maintenance of fluid and electrolyte balance and BP control. However, the role of peripheral circadian clocks in these homeostatic mechanisms remains unknown. We addressed this question in a mouse model carrying a conditional allele of the circadian clock gene Bmal1 and expressing Cre recombinase under the endogenous Renin promoter (Bmal1(lox/lox)/Ren1(d)Cre mice). Analysis of Bmal1(lox/lox)/Ren1(d)Cre mice showed that the floxed Bmal1 allele was excised in the kidney. In the kidney, BMAL1 protein expression was absent in the renin-secreting granular cells of the juxtaglomerular apparatus and the collecting duct. A partial reduction of BMAL1 expression was observed in the medullary thick ascending limb. Functional analyses showed that Bmal1(lox/lox)/Ren1(d)Cre mice exhibited multiple abnormalities, including increased urine volume, changes in the circadian rhythm of urinary sodium excretion, increased GFR, and significantly reduced plasma aldosterone levels. These changes were accompanied by a reduction in BP. These results show that local renal circadian clocks control body fluid and BP homeostasis.

BIO-INSPIRED COMPUTATIONAL TECHNIQUES APPLIED TO THE CLUSTERING AND VISUALIZATION OF SPATIO-TEMPORAL GEOSPATIAL DATA

The coverage and volume of geo-referenced datasets are extensive and incessantly¦growing. The systematic capture of geo-referenced information generates large volumes¦of spatio-temporal data to be analyzed. Clustering and visualization play a key¦role in the exploratory data analysis and the extraction of knowledge embedded in¦these data. However, new challenges in visualization and clustering are posed when¦dealing with the special characteristics of this data. For instance, its complex structures,¦large quantity of samples, variables involved in a temporal context, high dimensionality¦and large variability in cluster shapes.¦The central aim of my thesis is to propose new algorithms and methodologies for¦clustering and visualization, in order to assist the knowledge extraction from spatiotemporal¦geo-referenced data, thus improving making decision processes.¦I present two original algorithms, one for clustering: the Fuzzy Growing Hierarchical¦Self-Organizing Networks (FGHSON), and the second for exploratory visual data analysis:¦the Tree-structured Self-organizing Maps Component Planes. In addition, I present¦methodologies that combined with FGHSON and the Tree-structured SOM Component¦Planes allow the integration of space and time seamlessly and simultaneously in¦order to extract knowledge embedded in a temporal context.¦The originality of the FGHSON lies in its capability to reflect the underlying structure¦of a dataset in a hierarchical fuzzy way. A hierarchical fuzzy representation of¦clusters is crucial when data include complex structures with large variability of cluster¦shapes, variances, densities and number of clusters. The most important characteristics¦of the FGHSON include: (1) It does not require an a-priori setup of the number¦of clusters. (2) The algorithm executes several self-organizing processes in parallel.¦Hence, when dealing with large datasets the processes can be distributed reducing the¦computational cost. (3) Only three parameters are necessary to set up the algorithm.¦In the case of the Tree-structured SOM Component Planes, the novelty of this algorithm¦lies in its ability to create a structure that allows the visual exploratory data analysis¦of large high-dimensional datasets. This algorithm creates a hierarchical structure¦of Self-Organizing Map Component Planes, arranging similar variables' projections in¦the same branches of the tree. Hence, similarities on variables' behavior can be easily¦detected (e.g. local correlations, maximal and minimal values and outliers).¦Both FGHSON and the Tree-structured SOM Component Planes were applied in¦several agroecological problems proving to be very efficient in the exploratory analysis¦and clustering of spatio-temporal datasets.¦In this thesis I also tested three soft competitive learning algorithms. Two of them¦well-known non supervised soft competitive algorithms, namely the Self-Organizing¦Maps (SOMs) and the Growing Hierarchical Self-Organizing Maps (GHSOMs); and the¦third was our original contribution, the FGHSON. Although the algorithms presented¦here have been used in several areas, to my knowledge there is not any work applying¦and comparing the performance of those techniques when dealing with spatiotemporal¦geospatial data, as it is presented in this thesis.¦I propose original methodologies to explore spatio-temporal geo-referenced datasets¦through time. Our approach uses time windows to capture temporal similarities and¦variations by using the FGHSON clustering algorithm. The developed methodologies¦are used in two case studies. In the first, the objective was to find similar agroecozones¦through time and in the second one it was to find similar environmental patterns¦shifted in time.¦Several results presented in this thesis have led to new contributions to agroecological¦knowledge, for instance, in sugar cane, and blackberry production.¦Finally, in the framework of this thesis we developed several software tools: (1)¦a Matlab toolbox that implements the FGHSON algorithm, and (2) a program called¦BIS (Bio-inspired Identification of Similar agroecozones) an interactive graphical user¦interface tool which integrates the FGHSON algorithm with Google Earth in order to¦show zones with similar agroecological characteristics.

Measurements of seismic attenuation and transient fluid pressure in partially saturated Berea sandstone: Evidence of fluid flow on the mesoscopic scale

A novel laboratory technique is proposed to investigate wave-induced fluid flow on the mesoscopic scale as a mechanism for seismic attenuation in partially saturated rocks. This technique combines measurements of seismic attenuation in the frequency range from 1 to 100?Hz with measurements of transient fluid pressure as a response of a step stress applied on top of the sample. We used a Berea sandstone sample partially saturated with water. The laboratory results suggest that wave-induced fluid flow on the mesoscopic scale is dominant in partially saturated samples. A 3-D numerical model representing the sample was used to verify the experimental results. Biot's equations of consolidation were solved with the finite-element method. Wave-induced fluid flow on the mesoscopic scale was the only attenuation mechanism accounted for in the numerical solution. The numerically calculated transient fluid pressure reproduced the laboratory data. Moreover, the numerically calculated attenuation, superposed to the frequency-independent matrix anelasticity, reproduced the attenuation measured in the laboratory in the partially saturated sample. This experimental?numerical fit demonstrates that wave-induced fluid flow on the mesoscopic scale and matrix anelasticity are the dominant mechanisms for seismic attenuation in partially saturated Berea sandstone.

Do seismic waves and fluid flow sense the same permeability?

We study the discrepancy between the effective flow permeability and the effective seismic permeability, that is, the effective permeability controlling seismic attenuation due to wave-induced fluid flow, in 2D rock samples having mesoscopic heterogeneities and in the presence of strong permeability fluctuations. In order to do so, we employ a numerical oscillatory compressibility test to determine attenuation and velocity dispersion due to wave-induced fluid flow in these kinds of media and compare the responses with those obtained by replacing the heterogeneous permeability field by constant values, including the average permeability as well as the effective flow permeability of the sample. The latter is estimated in a separate upscaling procedure by solving the steady-state flow equation in the rock sample under study. Numerical experiments let us verify that attenuation levels are less significant and the attenuation peak gets broader in the presence of such strong permeability fluctuations. Moreover, we observe that for very low frequencies the effective seismic permeability is similar to the effective flow permeability, while for very high frequencies it approaches the arithmetic average of the permeability field.