Quantitative integration of geophysical and hydraulic data : from the local towards the regional scale range

Des progrès significatifs ont été réalisés dans le domaine de l'intégration quantitative des données géophysique et hydrologique l'échelle locale. Cependant, l'extension à de plus grandes échelles des approches correspondantes constitue encore un défi majeur. Il est néanmoins extrêmement important de relever ce défi pour développer des modèles fiables de flux des eaux souterraines et de transport de contaminant. Pour résoudre ce problème, j'ai développé une technique d'intégration des données hydrogéophysiques basée sur une procédure bayésienne de simulation séquentielle en deux étapes. Cette procédure vise des problèmes à plus grande échelle. L'objectif est de simuler la distribution d'un paramètre hydraulique cible à partir, d'une part, de mesures d'un paramètre géophysique pertinent qui couvrent l'espace de manière exhaustive, mais avec une faible résolution (spatiale) et, d'autre part, de mesures locales de très haute résolution des mêmes paramètres géophysique et hydraulique. Pour cela, mon algorithme lie dans un premier temps les données géophysiques de faible et de haute résolution à travers une procédure de réduction déchelle. Les données géophysiques régionales réduites sont ensuite reliées au champ du paramètre hydraulique à haute résolution. J'illustre d'abord l'application de cette nouvelle approche dintégration des données à une base de données synthétiques réaliste. Celle-ci est constituée de mesures de conductivité hydraulique et électrique de haute résolution réalisées dans les mêmes forages ainsi que destimations des conductivités électriques obtenues à partir de mesures de tomographic de résistivité électrique (ERT) sur l'ensemble de l'espace. Ces dernières mesures ont une faible résolution spatiale. La viabilité globale de cette méthode est testée en effectuant les simulations de flux et de transport au travers du modèle original du champ de conductivité hydraulique ainsi que du modèle simulé. Les simulations sont alors comparées. Les résultats obtenus indiquent que la procédure dintégration des données proposée permet d'obtenir des estimations de la conductivité en adéquation avec la structure à grande échelle ainsi que des predictions fiables des caractéristiques de transports sur des distances de moyenne à grande échelle. Les résultats correspondant au scénario de terrain indiquent que l'approche d'intégration des données nouvellement mise au point est capable d'appréhender correctement les hétérogénéitées à petite échelle aussi bien que les tendances à gande échelle du champ hydraulique prévalent. Les résultats montrent également une flexibilté remarquable et une robustesse de cette nouvelle approche dintégration des données. De ce fait, elle est susceptible d'être appliquée à un large éventail de données géophysiques et hydrologiques, à toutes les gammes déchelles. Dans la deuxième partie de ma thèse, j'évalue en détail la viabilité du réechantillonnage geostatique séquentiel comme mécanisme de proposition pour les méthodes Markov Chain Monte Carlo (MCMC) appliquées à des probmes inverses géophysiques et hydrologiques de grande dimension . L'objectif est de permettre une quantification plus précise et plus réaliste des incertitudes associées aux modèles obtenus. En considérant une série dexemples de tomographic radar puits à puits, j'étudie deux classes de stratégies de rééchantillonnage spatial en considérant leur habilité à générer efficacement et précisément des réalisations de la distribution postérieure bayésienne. Les résultats obtenus montrent que, malgré sa popularité, le réechantillonnage séquentiel est plutôt inefficace à générer des échantillons postérieurs indépendants pour des études de cas synthétiques réalistes, notamment pour le cas assez communs et importants où il existe de fortes corrélations spatiales entre le modèle et les paramètres. Pour résoudre ce problème, j'ai développé un nouvelle approche de perturbation basée sur une déformation progressive. Cette approche est flexible en ce qui concerne le nombre de paramètres du modèle et lintensité de la perturbation. Par rapport au rééchantillonage séquentiel, cette nouvelle approche s'avère être très efficace pour diminuer le nombre requis d'itérations pour générer des échantillons indépendants à partir de la distribution postérieure bayésienne. - Significant progress has been made with regard to the quantitative integration of geophysical and hydrological data at the local scale. However, extending corresponding approaches beyond the local scale still represents a major challenge, yet is critically important for the development of reliable groundwater flow and contaminant transport models. To address this issue, I have developed a hydrogeophysical data integration technique based on a two-step Bayesian sequential simulation procedure that is specifically targeted towards larger-scale problems. The objective is to simulate the distribution of a target hydraulic parameter based on spatially exhaustive, but poorly resolved, measurements of a pertinent geophysical parameter and locally highly resolved, but spatially sparse, measurements of the considered geophysical and hydraulic parameters. To this end, my algorithm links the low- and high-resolution geophysical data via a downscaling procedure before relating the downscaled regional-scale geophysical data to the high-resolution hydraulic parameter field. I first illustrate the application of this novel data integration approach to a realistic synthetic database consisting of collocated high-resolution borehole measurements of the hydraulic and electrical conductivities and spatially exhaustive, low-resolution electrical conductivity estimates obtained from electrical resistivity tomography (ERT). The overall viability of this method is tested and verified by performing and comparing flow and transport simulations through the original and simulated hydraulic conductivity fields. The corresponding results indicate that the proposed data integration procedure does indeed allow for obtaining faithful estimates of the larger-scale hydraulic conductivity structure and reliable predictions of the transport characteristics over medium- to regional-scale distances. The approach is then applied to a corresponding field scenario consisting of collocated high- resolution measurements of the electrical conductivity, as measured using a cone penetrometer testing (CPT) system, and the hydraulic conductivity, as estimated from electromagnetic flowmeter and slug test measurements, in combination with spatially exhaustive low-resolution electrical conductivity estimates obtained from surface-based electrical resistivity tomography (ERT). The corresponding results indicate that the newly developed data integration approach is indeed capable of adequately capturing both the small-scale heterogeneity as well as the larger-scale trend of the prevailing hydraulic conductivity field. The results also indicate that this novel data integration approach is remarkably flexible and robust and hence can be expected to be applicable to a wide range of geophysical and hydrological data at all scale ranges. In the second part of my thesis, I evaluate in detail the viability of sequential geostatistical resampling as a proposal mechanism for Markov Chain Monte Carlo (MCMC) methods applied to high-dimensional geophysical and hydrological inverse problems in order to allow for a more accurate and realistic quantification of the uncertainty associated with the thus inferred models. Focusing on a series of pertinent crosshole georadar tomographic examples, I investigated two classes of geostatistical resampling strategies with regard to their ability to efficiently and accurately generate independent realizations from the Bayesian posterior distribution. The corresponding results indicate that, despite its popularity, sequential resampling is rather inefficient at drawing independent posterior samples for realistic synthetic case studies, notably for the practically common and important scenario of pronounced spatial correlation between model parameters. To address this issue, I have developed a new gradual-deformation-based perturbation approach, which is flexible with regard to the number of model parameters as well as the perturbation strength. Compared to sequential resampling, this newly proposed approach was proven to be highly effective in decreasing the number of iterations required for drawing independent samples from the Bayesian posterior distribution.

Estimation of the lateral correlation structure of subsurface water content from surface-based ground-penetrating radar reflection images

Over the past decade, significant interest has been expressed in relating the spatial statistics of surface-based reflection ground-penetrating radar (GPR) data to those of the imaged subsurface volume. A primary motivation for this work is that changes in the radar wave velocity, which largely control the character of the observed data, are expected to be related to corresponding changes in subsurface water content. Although previous work has indeed indicated that the spatial statistics of GPR images are linked to those of the water content distribution of the probed region, a viable method for quantitatively analyzing the GPR data and solving the corresponding inverse problem has not yet been presented. Here we address this issue by first deriving a relationship between the 2-D autocorrelation of a water content distribution and that of the corresponding GPR reflection image. We then show how a Bayesian inversion strategy based on Markov chain Monte Carlo sampling can be used to estimate the posterior distribution of subsurface correlation model parameters that are consistent with the GPR data. Our results indicate that if the underlying assumptions are valid and we possess adequate prior knowledge regarding the water content distribution, in particular its vertical variability, this methodology allows not only for the reliable recovery of lateral correlation model parameters but also for estimates of parameter uncertainties. In the case where prior knowledge regarding the vertical variability of water content is not available, the results show that the methodology still reliably recovers the aspect ratio of the heterogeneity.

An evaluation of new parsimony-based versus parametric inference methods in biogeography: a case study using the globally distributed plant family Sapindaceae

Aim  Recently developed parametric methods in historical biogeography allow researchers to integrate temporal and palaeogeographical information into the reconstruction of biogeographical scenarios, thus overcoming a known bias of parsimony-based approaches. Here, we compare a parametric method, dispersal-extinction-cladogenesis (DEC), against a parsimony-based method, dispersal-vicariance analysis (DIVA), which does not incorporate branch lengths but accounts for phylogenetic uncertainty through a Bayesian empirical approach (Bayes-DIVA). We analyse the benefits and limitations of each method using the cosmopolitan plant family Sapindaceae as a case study.Location  World-wide.Methods  Phylogenetic relationships were estimated by Bayesian inference on a large dataset representing generic diversity within Sapindaceae. Lineage divergence times were estimated by penalized likelihood over a sample of trees from the posterior distribution of the phylogeny to account for dating uncertainty in biogeographical reconstructions. We compared biogeographical scenarios between Bayes-DIVA and two different DEC models: one with no geological constraints and another that employed a stratified palaeogeographical model in which dispersal rates were scaled according to area connectivity across four time slices, reflecting the changing continental configuration over the last 110 million years.Results  Despite differences in the underlying biogeographical model, Bayes-DIVA and DEC inferred similar biogeographical scenarios. The main differences were: (1) in the timing of dispersal events - which in Bayes-DIVA sometimes conflicts with palaeogeographical information, and (2) in the lower frequency of terminal dispersal events inferred by DEC. Uncertainty in divergence time estimations influenced both the inference of ancestral ranges and the decisiveness with which an area can be assigned to a node.Main conclusions  By considering lineage divergence times, the DEC method gives more accurate reconstructions that are in agreement with palaeogeographical evidence. In contrast, Bayes-DIVA showed the highest decisiveness in unequivocally reconstructing ancestral ranges, probably reflecting its ability to integrate phylogenetic uncertainty. Care should be taken in defining the palaeogeographical model in DEC because of the possibility of overestimating the frequency of extinction events, or of inferring ancestral ranges that are outside the extant species ranges, owing to dispersal constraints enforced by the model. The wide-spanning spatial and temporal model proposed here could prove useful for testing large-scale biogeographical patterns in plants.

Focal dependent pleural thickening at MDCT: pleural lesion or functional abnormality?

PURPOSE: To describe the characteristics of reversible focal pleural thickenings (PTs) mimicking real plaques, that firstly suggest asbestos exposure or pleural metastasis; to propose an imaging strategy and propose an explanation for their mechanism of formation. PATIENTS AND METHODS: Retrospective review of data from 19 patients with PTs fitting the description of pleural plaques at chest computed tomography (CT) and presenting modifications (clearance or appearance) of at least one PT at an additional chest examination in prone position. RESULTS: A total of 152 PTs were recorded on the first chest CT examinations with a range of two to 19 pleural opacities per patient. All PTs had a posterior distribution in the lower lobes. On the additional acquisitions, 144 PTs disappeared. Seventeen patients presented complete regression of PTs and two patients presented persistence of eight PTs. CONCLUSION: Additional low dose acquisition in prone position should be performed in all patients presenting with focal PT in a dependent and basal location. This may allow to exclude a pleural plaque in case of asbestos exposure but also a pleural metastasis in oncologic patients. These reversible dependent PTs could be related to physiological focal accumulation of lymphatic fluid in subpleural area.

Pharmacokinetics and posterior segment biodistribution of ESBA105, an anti-TNF-alpha single-chain antibody, upon topical administration to the rabbit eye.

PURPOSE: The purpose of this study was to characterize local distribution and systemic absorption of the tumor necrosis factor (TNF)-alpha inhibitory single-chain antibody fragment (scFv) ESBA105 following topical administration to the eye in vivo. METHODS: Rabbits received ESBA105 as topical eye drops in two dosing regimens. First, pharmacokinetics after the topical route of administration was compared to the intravenous (i.v.) route by means of applying the identical cumulative daily dose of ESBA105. In a second study rabbits received five eye drops daily for six consecutive days in a lower frequency topical dosing regimen. Kinetics and biodistribution of ESBA105 in ocular tissues and fluids as well as in sera were determined in all animals. RESULTS: After topical administration to the eye, ESBA105 quickly reaches therapeutic concentrations in all ocular compartments. Systemic exposure after topical administration is 25,000-fold lower than exposure after i.v. injection of the identical cumulative daily dose. ESBA105 levels in vitreous humor and neuroretina are significantly higher on topical administration than after i.v. injection. Absolute and relative intraocular biodistribution of ESBA105 is different with topical and systemic delivery routes. Compared to its terminal half-life in circulation (7 hours), the vitreal half-life of ESBA105 is significantly enhanced (16-24 hours). CONCLUSIONS: On topical administration, ESBA105 is efficiently absorbed and distributed to all compartments of the eye, whereby systemic drug exposure is very low. Based on its unique intraocular biodistribution and pharmacokinetics and the absolute intraocular levels reached, topical ESBA105 appears highly attractive for treatment of various ophthalmological disorders.

Changes in the transcriptional profile of transporters in the intestine along the anterior-posterior and crypt-villus axes.

BACKGROUND: The purpose of this work was to characterize the expression of drug and nutrient carriers along the anterior-posterior and crypt-villus axes of the intestinal epithelium and to study the validity of utilizing whole gut tissue rather than purified epithelial cells to examine regional variations in gene expression. RESULTS: We have characterized the mRNA expression profiles of 76 % of all currently known transporters along the anterior-posterior axis of the gut. This is the first study to describe the expression profiles of the majority of all known transporters in the intestine. The expression profiles of transporters, as defined according to the Gene Ontology consortium, were measured in whole tissue of the murine duodenum, jejunum, ileum and colon using high-density microarrays. For nine transporters (Abca1, Abcc1, Abcc3, Abcg8, Slc10a2, Slc28a2, Slc2a1, Slc34a2 and Slc5a8), the mRNA profiles were further measured by RT-PCR in laser micro-dissected crypt and villus epithelial cells corresponding to the aforementioned intestinal regions. With respect to differentially regulated transporters, the colon had a distinct expression profile from small intestinal segments. The majority (59 % for p cutoff < or = 0.05) of transporter mRNA levels were constant across the intestinal sections studied. For the transporter subclass "carrier activity", which contains the majority of known carriers for biologically active compounds, a significant change (p < or = 0.05) along the anterior-posterior axis was observed. CONCLUSION: All nine transporters examined in laser-dissected material demonstrated good replication of the region-specific profiles revealed by microarray. Furthermore, we suggest that the distribution characteristics of Slc5a8 along the intestinal tract render it a suitable candidate carrier for monocarboxylate drugs in the posterior portion of the intestine. Our findings also predict that there is a significant difference in the absorption of carrier-mediated compounds in the different intestinal segments. The most pronounced differences can be expected between the adjoining segments ileum and colon, but the differences between the other adjoining segments are not negligible. Finally, for the examined genes, profiles measured in whole intestinal tissue extracts are representative of epithelial cell-only gene expression.

Ocular and systemic bio-distribution of rhodamine-conjugated liposomes loaded with VIP injected into the vitreous of Lewis rats.

PURPOSE: Local delivery of therapeutic molecules encapsulated within liposomes is a promising method to treat ocular inflammation. The purpose of the present study was to define the biodistribution of rhodamine-conjugated liposomes loaded with vasoactive intestinal peptide (VIP), an immunosuppressive neuropeptide, following their intravitreal (IVT) injection in normal rats. METHODS: Healthy seven- to eight-week-old Lewis male rats were injected into the vitreous with empty rhodamine-conjugated liposomes (Rh-Lip) or with VIP-loaded Rh-Lip (VIP-Rh-Lip; 50 mM of lipids with an encapsulation efficiency of 3.0+/-0.4 mmol VIP/mol lipids). Twenty-four h after IVT injection, the eyes, the cervical, mesenteric, and inguinal lymph nodes (LN), and spleen were collected. The phenotype and distribution of cells internalizing Rh-Lip and VIP-Rh-Lip were studied. Determination of VIP expression in ocular tissues and lymphoid organs and interactions with T cells in cervical LN was performed on whole mounted tissues and frozen tissue sections by immunofluorescence and confocal microscopy. RESULTS: In the eye, 24 h following IVT injection, fluorescent liposomes (Rh-Lip and VIP-Rh-Lip) were detected mainly in the posterior segment of the eye (vitreous, inner layer of the retina) and to a lesser extent at the level of the iris root and ciliary body. Liposomes were internalized by activated retinal Müller glial cells, ocular tissue resident macrophages, and rare infiltrating activated macrophages. In addition, fluorescent liposomes were found in the episclera and conjunctiva where free VIP expression was also detected. In lymphoid organs, Rh-Lip and VIP-Rh-Lip were distributed almost exclusively in the cervical lymph nodes (LN) with only a few Rh-Lip-positive cells detected in the spleen and mesenteric LN and none in the inguinal LN. In the cervical LN, Rh-Lip were internalized by resident ED3-positive macrophages adjacent to CD4 and CD8-positive T lymphocytes. Some of these T lymphocytes in close contact with macrophages containing VIP-Rh-Lip expressed VIP. CONCLUSIONS: Liposomes are specifically internalized by retinal Müller glial cells and resident macrophages in the eye. A limited passage of fluorescent liposomes from the vitreous to the spleen via the conventional outflow pathway and the venous circulation was detected. The majority of fluorescent liposomes deposited in the conjunctiva following IVT injection reached the subcapsular sinus of the cervical LN via conjuntival lymphatics. In the cervical LN, Rh-Lip were internalized by resident subcapsular sinus macrophages adjacent to T lymphocytes. Detection of VIP in both macrophages and T cells in cervical LN suggests that IVT injection of VIP-Rh-Lip may increase ocular immune privilege by modulating the loco-regional immune environment. In conclusion, our observations suggest that IVT injection of VIP-loaded liposomes is a promising therapeutic strategy to dampen ocular inflammation by modulating macrophage and T cell activation mainly in the loco-regional immune system.

Mass conservative three-dimensional water tracer distribution from Markov chain Monte Carlo inversion of time-lapse ground-penetrating radar data

Time-lapse geophysical measurements are widely used to monitor the movement of water and solutes through the subsurface. Yet commonly used deterministic least squares inversions typically suffer from relatively poor mass recovery, spread overestimation, and limited ability to appropriately estimate nonlinear model uncertainty. We describe herein a novel inversion methodology designed to reconstruct the three-dimensional distribution of a tracer anomaly from geophysical data and provide consistent uncertainty estimates using Markov chain Monte Carlo simulation. Posterior sampling is made tractable by using a lower-dimensional model space related both to the Legendre moments of the plume and to predefined morphological constraints. Benchmark results using cross-hole ground-penetrating radar travel times measurements during two synthetic water tracer application experiments involving increasingly complex plume geometries show that the proposed method not only conserves mass but also provides better estimates of plume morphology and posterior model uncertainty than deterministic inversion results.

Cellular distribution of Hsp70 expression in rat skeletal muscles. Effects of moderate exercise training and chronic hypoxia.

Rat hindlimb muscles constitutively express the inducible heat shock protein 72 (Hsp70), apparently in proportion to the slow myosin content. Since it remains controversial whether chronic Hsp70 expression reflects the overimposed stress, we investigated Hsp70 cellular distribution in fast muscles of the posterior rat hindlimb after (1) mild exercise training (up to 30 m/min treadmill run for 1 h/day), which induces a remodeling in fast fiber composition, or (2) prolonged exposure to normobaric hypoxia (10%O(2)), which does not affect fiber-type composition. Both conditions increased significantly protein Hsp70 levels in the skeletal muscle. Immunohistochemistry showed the labeling for Hsp70 in subsets of both slow/type 1 and fast/type 2A myofibers of control, sedentary, and normoxic rats. Endurance training increased about threefold the percentage of Hsp70-positive myofibers (P < 0.001), and changed the distribution of Hsp70 immunoreactivity, which involved a larger subset of both type 2A and intermediate type 2A/2X myofibers (P < 0.001) and vascular smooth muscle cells. Hypoxia induced Hsp70 immunoreactivity in smooth muscle cells of veins and did not increase the percentage of Hsp70-positive myofibers; however, sustained exposure to hypoxia affected the distribution of Hsp70 immunoreactivity, which appeared detectable in a very small subset of type 2A fibers, whereas it concentrated in type 1 myofibers (P < 0.05) together with the labeling for heme-oxygenase isoform 1, a marker of oxidative stress. Therefore, the chronic induction of Hsp70 expression in rat skeletal muscles is not obligatory related to the slow fiber phenotype but reveals the occurrence of a stress response.

Anatomical distribution of areas of preserved cartilage in advanced femorotibial osteoarthritis using CT arthrography (Part 1).

OBJECTIVE: To determine subregions of normal and abnormal cartilage in advanced stages of femorotibial osteoarthritis (OA) by mapping the entire femorotibial joint in a cohort of pre-total knee replacement (TKR) OA knees. DESIGN: We defined an areal subdivision of the femorotibial articular cartilage surface on CT arthrography (CTA), allowing the division of the femorotibial articular surface into multiple (up to n = 204 per knee) subregions and the comparison of the same areas between different knees. Two readers independently classified each cartilage area as normal, abnormal or non-assessable in 41 consecutive pre-TKR OA knees. RESULTS: A total of 6447 cartilage areas (from 41 knees) were considered assessable by both readers. The average proportion of preserved cartilage was lower in the medial femorotibial joint than in the lateral femorotibial joint for both readers (32.0/69.8% and 33.9/68.5% (medial/lateral) for reader 1 and 2 respectively, all P < 0.001). High frequencies of normal cartilage were observed at the posterior aspect of the medial condyle (up to 89%), and the anterior aspect of the lateral femorotibial compartment (up to 100%). The posterior aspect of the medial condyle was the area that most frequently exhibited preserved cartilage in the medial femorotibial joint, contrasting with the high frequency of cartilage lesions in the rest of that compartment. CONCLUSIONS: Cartilage at the posterior aspect of the medial condyle, and at the anterior aspect of the lateral femorotibial compartment, may be frequently preserved in advanced grades of OA.

Ambulatory assessment of 3D ground reaction force using plantar pressure distribution.

This study aimed to use the plantar pressure insole for estimating the three-dimensional ground reaction force (GRF) as well as the frictional torque (T(F)) during walking. Eleven subjects, six healthy and five patients with ankle disease participated in the study while wearing pressure insoles during several walking trials on a force-plate. The plantar pressure distribution was analyzed and 10 principal components of 24 regional pressure values with the stance time percentage (STP) were considered for GRF and T(F) estimation. Both linear and non-linear approximators were used for estimating the GRF and T(F) based on two learning strategies using intra-subject and inter-subjects data. The RMS error and the correlation coefficient between the approximators and the actual patterns obtained from force-plate were calculated. Our results showed better performance for non-linear approximation especially when the STP was considered as input. The least errors were observed for vertical force (4%) and anterior-posterior force (7.3%), while the medial-lateral force (11.3%) and frictional torque (14.7%) had higher errors. The result obtained for the patients showed higher error; nevertheless, when the data of the same patient were used for learning, the results were improved and in general slight differences with healthy subjects were observed. In conclusion, this study showed that ambulatory pressure insole with data normalization, an optimal choice of inputs and a well-trained nonlinear mapping function can estimate efficiently the three-dimensional ground reaction force and frictional torque in consecutive gait cycle without requiring a force-plate.

Distribution of the human intracellular serpin protease inhibitor 8 in human tissues.

Ovalbumin-like serine protease inhibitors are mainly localized intracellularly and their in vivo functions are largely unknown. To elucidate their physiological role(s), we studied the expression of one of these inhibitors, protease inhibitor 8 (PI-8), in normal human tissues by immunohistochemistry using a PI-8-specific monoclonal antibody. PI-8 was strongly expressed in the nuclei of squamous epithelium of mouth, pharynx, esophagus, and epidermis, and by the epithelial layer of skin appendages, particularly by more differentiated epithelial cells. PI-8 was also expressed by monocytes and by neuroendocrine cells in the pituitary gland, pancreas, and digestive tract. Monocytes showed nuclear and cytoplasmic localization of PI-8, whereas neuroendocrine cells showed only cytoplasmic staining. In vitro nuclear localization of PI-8 was confirmed by confocal analysis using serpin-transfected HeLa cells. Furthermore, mutation of the P(1) residue did not affect the subcellular distribution pattern of PI-8, indicating that its nuclear localization is independent of the interaction with its target protease. We conclude that PI-8 has a unique distribution pattern in human tissues compared to the distribution patterns of other intracellular serpins. Additional studies must be performed to elucidate its physiological role.

Species distribution models reveal apparent competitive and facilitative effects of a dominant species on the distribution of tundra plants

Abiotic factors are considered strong drivers of species distribution and assemblages. Yet these spatial patterns are also influenced by biotic interactions. Accounting for competitors or facilitators may improve both the fit and the predictive power of species distribution models (SDMs). We investigated the influence of a dominant species, Empetrum nigrum ssp. hermaphroditum, on the distribution of 34 subordinate species in the tundra of northern Norway. We related SDM parameters of those subordinate species to their functional traits and their co-occurrence patterns with E. hermaphroditum across three spatial scales. By combining both approaches, we sought to understand whether these species may be limited by competitive interactions and/or benefit from habitat conditions created by the dominant species. The model fit and predictive power increased for most species when the frequency of occurrence of E. hermaphroditum was included in the SDMs as a predictor. The largest increase was found for species that 1) co-occur most of the time with E. hermaphroditum, both at large (i.e. 750 m) and small spatial scale (i.e. 2 m) or co-occur with E. hermaphroditum at large scale but not at small scale and 2) have particularly low or high leaf dry matter content (LDMC). Species that do not co-occur with E. hermaphroditum at the smallest scale are generally palatable herbaceous species with low LDMC, thus showing a weak ability to tolerate resource depletion that is directly or indirectly induced by E. hermaphroditum. Species with high LDMC, showing a better aptitude to face resource depletion and grazing, are often found in the proximity of E. hermaphroditum. Our results are consistent with previous findings that both competition and facilitation structure plant distribution and assemblages in the Arctic tundra. The functional and co-occurrence approaches used were complementary and provided a deeper understanding of the observed patterns by refinement of the pool of potential direct and indirect ecological effects of E. hermaphroditum on the distribution of subordinate species. Our correlative study would benefit being complemented by experimental approaches.

Low-Dose Vascular Photodynamic Therapy Decreases Tumor Interstitial Fluid Pressure, which Promotes Liposomal Doxorubicin Distribution in a Murine Sarcoma Metastasis Model.

INTRODUCTION: Solid tumors are known to have an abnormal vasculature that limits the distribution of chemotherapy. We have recently shown that tumor vessel modulation by low-dose photodynamic therapy (L-PDT) could improve the uptake of macromolecular chemotherapeutic agents such as liposomal doxorubicin (Liporubicin) administered subsequently. However, how this occurs is unknown. Convection, the main mechanism for drug transport between the intravascular and extravascular spaces, is mostly related to interstitial fluid pressure (IFP) and tumor blood flow (TBF). Here, we determined the changes of tumor and surrounding lung IFP and TBF before, during, and after vascular L-PDT. We also evaluated the effect of these changes on the distribution of Liporubicin administered intravenously (IV) in a lung sarcoma metastasis model. MATERIALS AND METHODS: A syngeneic methylcholanthrene-induced sarcoma cell line was implanted subpleurally in the lung of Fischer rats. Tumor/surrounding lung IFP and TBF changes induced by L-PDT were determined using the wick-in-needle technique and laser Doppler flowmetry, respectively. The spatial distribution of Liporubicin in tumor and lung tissues following IV drug administration was then assessed in L-PDT-pretreated animals and controls (no L-PDT) by epifluorescence microscopy. RESULTS: L-PDT significantly decreased tumor but not lung IFP compared to controls (no L-PDT) without affecting TBF. These conditions were associated with a significant improvement in Liporubicin distribution in tumor tissues compared to controls (P < .05). DISCUSSION: L-PDT specifically enhanced convection in blood vessels of tumor but not of normal lung tissue, which was associated with a significant improvement of Liporubicin distribution in tumors compared to controls.