Plutonium and americium concentrations in sea water and in sinking particles

Vertical fluxes of 239+240Pu and 241Am and temporal changes in their inventories in the northwestern Mediterranean Sea have been examined through high-resolution water column sampling coupled with direct measurements of the vertical flux of particle-bound transuranics using time-series sediment traps. Water column profiles of both radionuclides showed well-defined sub-surface maxima (2391240Pu between 100-400 m; 241Am at 100-200 m and 800 m), the depths of which are a result of the different biogeochemical scavenging behavior of the two radionuclides. Comparison of deep water column (0-2,000 m) transuranic inventories with those derived from earlier measurements demonstrate that the total 2391240Pu inventory had not substantially changed between 1976-1990 whereas 241Am had decreased by approximately 24%. Enhanced scavenging of 241Am and a resultant, more rapid removal from the water column relative to 239+240Pu was also supported by the observation of elevated Am/Pu activity ratios in sinking particles collected in sediment traps at depth. Direct measurements of the downward flux of particulate 239+240Pu and 241Am compared with transuranic removal rates derived from observed total water column inventory differences over time, show that particles sinking out of deep waters (1,000-2,000 m) could account for 26-72% of the computed total annual 239+240Pu loss and virtually all of the 241Am removal from the water column. Upper water column (0-200 m) residence times based on direct flux measurements ranged from 20-30 yr for 239+240Pu and 5-10 yr for 241Am. The observation that 241Am/239+240Pu activity ratios in unfiltered Mediterranean seawater are six times lower than those in the north Pacific suggests the existence of a specific mechanism for enhanced scavenging and removal of 241Am from the generally oligotrophic waters of the open Mediterranean. It is proposed that atmospheric inputs of aluminosilicate particles transported by Saharan dust events which frequently occur in the Mediterranean region could enhance the geochemical scavenging and resultant removal of 241Am to the sediments.

Glacier mass balances of Stubacher Sonnblickkees, Hohe Tauern Range, Eastern Alps, Austria, 1958/1959 to 2012/2013

Stubacher Sonnblickkees (SSK) is located in the Hohe Tauern Range (Eastern Alps) in the south of Salzburg Province (Austria) in the region of Oberpinzgau in the upper Stubach Valley. The glacier is situated at the main Alpine crest and faces east, starting at elevations close to 3050 m and in the 1980s terminated at 2500 m a.s.l. It had an area of 1.7 km² at that time, compared with 1 km² in 2013. The glacier type can be classified as a slope glacier, i.e. the relief is covered by a relatively thin ice sheet and there is no regular glacier tongue. The rough subglacial topography makes for a complex shape in the surface topography, with various concave and convex patterns. The main reason for selecting this glacier for mass balance observations (as early as 1963) was to verify on a complex glacier how the mass balance methods and the conclusions - derived during the more or less pioneer phase of glaciological investigations in the 1950s and 1960s - could be applied to the SSK glacier. The decision was influenced by the fact that close to the SSK there was the Rudolfshütte, a hostel of the Austrian Alpine Club (OeAV), newly constructed in the 1950s to replace the old hut dating from 1874. The new Alpenhotel Rudolfshütte, which was run by the Slupetzky family from 1958 to 1970, was the base station for the long-term observation; the cable car to Rudolfshütte, operated by the Austrian Federal Railways (ÖBB), was a logistic advantage. Another factor for choosing SSK as a glaciological research site was the availability of discharge records of the catchment area from the Austrian Federal Railways who had turned the nearby lake Weißsee ('White Lake') - a former natural lake - into a reservoir for their hydroelectric power plants. In terms of regional climatic differences between the Central Alps in Tyrol and those of the Hohe Tauern, the latter experienced significantly higher precipitation , so one could expect new insights in the different response of the two glaciers SSK and Hintereisferner (Ötztal Alps) - where a mass balance series went back to 1952. In 1966 another mass balance series with an additional focus on runoff recordings was initiated at Vernagtfener, near Hintereisferner, by the Commission of the Bavarian Academy of Sciences in Munich. The usual and necessary link to climate and climate change was given by a newly founded weather station (by Heinz and Werner Slupetzky) at the Rudolfshütte in 1961, which ran until 1967. Along with an extension and enlargement to the so-called Alpine Center Rudolfshütte of the OeAV, a climate observatory (suggested by Heinz Slupetzky) has been operating without interruption since 1980 under the responsibility of ZAMG and the Hydrological Service of Salzburg, providing long-term met observations. The weather station is supported by the Berghotel Rudolfshütte (in 2004 the OeAV sold the hotel to a private owner) with accommodation and facilities. Direct yearly mass balance measurements were started in 1963, first for 3 years as part of a thesis project. In 1965 the project was incorporated into the Austrian glacier measurement sites within the International Hydrological Decade (IHD) 1965 - 1974 and was afterwards extended via the International Hydrological Program (IHP) 1975 - 1981. During both periods the main financial support came from the Hydrological Survey of Austria. After 1981 funds were provided by the Hydrological Service of the Federal Government of Salzburg. The research was conducted from 1965 onwards by Heinz Slupetzky from the (former) Department of Geography of the University of Salzburg. These activities received better recognition when the High Alpine Research Station of the University of Salzburg was founded in 1982 and brought in additional funding from the University. With recent changes concerning Rudolfshütte, however, it became unfeasible to keep the research station going. Fortunately, at least the weather station at Rudolfshütte is still operating. In the pioneer years of the mass balance recordings at SSK, the main goal was to understand the influence of the complicated topography on the ablation and accumulation processes. With frequent strong southerly winds (foehn) on the one hand, and precipitation coming in with storms from the north to northwest, the snow drift is an important factor on the undulating glacier surface. This results in less snow cover in convex zones and in more or a maximum accumulation in concave or flat areas. As a consequence of the accentuated topography, certain characteristic ablation and accumulation patterns can be observed during the summer season every year, which have been regularly observed for many decades . The process of snow depletion (Ausaperung) runs through a series of stages (described by the AAR) every year. The sequence of stages until the end of the ablation season depends on the weather conditions in a balance year. One needs a strong negative mass balance year at the beginning of glacier measurements to find out the regularities; 1965, the second year of observation resulted in a very positive mass balance with very little ablation but heavy accumulation. To date it is the year with the absolute maximum positive balance in the entire mass balance series since 1959, probably since 1950. The highly complex ablation patterns required a high number of ablation stakes at the beginning of the research and it took several years to develop a clearer idea of the necessary density of measurement points to ensure high accuracy. A great number of snow pits and probing profiles (and additional measurements at crevasses) were necessary to map the accumulation area/patterns. Mapping the snow depletion, especially at the end of the ablation season, which coincides with the equilibrium line, is one of the main basic data for drawing contour lines of mass balance and to calculate the total mass balance (on a regular-shaped valley glacier there might be an equilibrium line following a contour line of elevation separating the accumulation area and the ablation area, but not at SSK). - An example: in 1969/70, 54 ablation stakes and 22 snow pits were used on the 1.77 km² glacier surface. In the course of the study the consistency of the accumulation and ablation patterns could be used to reduce the number of measurement points. - At the SSK the stratigraphic system, i.e. the natural balance year, is used instead the usual hydrological year. From 1964 to 1981, the yearly mass balance was calculated by direct measurements. Based on these records of 17 years, a regression analysis between the specific net mass balance and the ratio of ablation area to total area (AAR) has been used since then. The basic requirement was mapping the maximum snow depletion at the end of each balance year. There was the advantage of Heinz Slupetzky's detailed local and long-term experience, which ensured homogeneity of the series on individual influences of the mass balance calculations. Verifications took place as often as possible by means of independent geodetic methods, i.e. monoplotting , aerial and terrestrial photogrammetry, more recently also the application of PHOTOMODELLER and laser scans. The semi-direct mass balance determinations used at SSK were tentatively compared with data from periods of mass/volume change, resulting in promising first results on the reliability of the method. In recent years re-analyses of the mass balance series have been conducted by the World Glacier Monitoring Service and will be done at SSK too. - The methods developed at SSK also add to another objective, much discussed in the 1960s within the community, namely to achieve time- and labour-saving methods to ensure continuation of long-term mass balance series. The regression relations were used to extrapolate the mass balance series back to 1959, the maximum depletion could be reconstructed by means of photographs for those years. R. Günther (1982) calculated the mass balance series of SSK back to 1950 by analysing the correlation between meteorological data and the mass balance; he found a high statistical relation between measured and determined mass balance figures for SSK. In spite of the complex glacier topography, interesting empirical experiences were gained from the mass balance data sets, giving a better understanding of the characteristics of the glacier type, mass balance and mass exchange. It turned out that there are distinct relations between the specific net balance, net accumulation (defined as Bc/S) and net ablation (Ba/S) to the AAR, resulting in characteristic so-called 'turnover curves'. The diagram of SSK represents the type of a glacier without a glacier tongue. Between 1964 and 1966, a basic method was developed, starting from the idea that instead of measuring years to cover the range between extreme positive and extreme negative yearly balances one could record the AAR/snow depletion/Ausaperung during one or two summers. The new method was applied on Cathedral Massif Glacier, a cirque glacier with the same area as the Stubacher Sonnblickkees, in British Columbia, Canada. during the summers of 1977 and 1978. It returned exactly the expected relations, e.g. mass turnover curves, as found on SSK. The SSK was mapped several times on a scale of 1:5000 to 1:10000. Length variations have been measured since 1960 within the OeAV glacier length measurement programme. Between 1965 and 1981, there was a mass gain of 10 million cubic metres. With a time lag of 10 years, this resulted in an advance until the mid-1980s. Since 1982 there has been a distinct mass loss of 35 million cubic metres by 2013. In recent years, the glacier has disintegrated faster, forced by the formation of a periglacial lake at the glacier terminus and also by the outcrops of rocks (typical for the slope glacier type), which have accelerated the meltdown. The formation of this lake is well documented. The glacier has retreated by some 600 m since 1981. - Since August 2002, a runoff gauge installed by the Hydrographical Service of Salzburg has recorded the discharge of the main part of SSK at the outlet of the new Unterer Eisboden See. The annual reports - submitted from 1982 on as a contractual obligation to the Hydrological Service of Salzburg - document the ongoing processes on the one hand, and emphasize the mass balance of SSK and outline the climatological reasons, mainly based on the met-data of the observatory Rudolfshütte, on the other. There is an additional focus on estimating the annual water balance in the catchment area of the lake. There are certain preconditions for the water balance equation in the area. Runoff is recorded by the ÖBB power stations, the mass balance of the now approx. 20% glaciated area (mainly the Sonnblickkees) is measured andthe change of the snow and firn patches/the water content is estimated as well as possible. (Nowadays laserscanning and ground radar are available to measure the snow pack). There is a net of three precipitation gauges plus the recordings at Rudolfshütte. The evaporation is of minor importance. The long-term annual mean runoff depth in the catchment area is around 3.000 mm/year. The precipitation gauges have measured deficits between 10% and 35%, on average probably 25% to 30%. That means that the real precipitation in the catchment area Weißsee (at elevations between 2,250 and 3,000 m) is in an order of 3,200 to 3,400 mm a year. The mass balance record of SSK was the first one established in the Hohe Tauern region (and now since the Hohe Tauern National Park was founded in 1983 in Salzburg) and is one of the longest measurement series worldwide. Great efforts are under way to continue the series, to safeguard against interruption and to guarantee a long-term monitoring of the mass balance and volume change of SSK (until the glacier is completely gone, which seems to be realistic in the near future as a result of the ongoing global warming). Heinz Slupetzky, March 2014

Chlorophyll a concentration, and proportion of time in the water column obtained by instrumented bowhead whales

Sixty hours of direct measurements of fluorescence were collected from six bowhead whales (Balaena mysticetus) instrumented with fluorometers in Greenland in April 2005 and 2006. The data were used to (1) characterize the three-dimensional spatial pattern of chlorophyll-a (Chl-a) in the water column, (2) to examine the relationships between whale foraging areas and productive zones, and (3) to examine the correlation between whale-derived in situ values of Chl-a and those from concurrent satellite images using the NASA MODIS (Moderate Resolution Imaging Spectroradiometer) EOS-AQUA satellite (MOD21, SeaWifs analogue OC3M and SST MOD37). Bowhead whales traversed 1600 km**2, providing information on diving, Chl-a structure and temperature profiles to depths below 200 m. Feeding dives frequently passed through surface waters ( >50 m) and targeted depths close to the bottom, and whales did not always target patches of high concentrations of Chl-a in the upper 50 m. Five satellite images were available within the periods whales carried fluorometers. Whales traversed 91 pixels collecting on average 761 s (SD 826) of Chl-a samples per pixel (0-136 m). The depth of the Chl-a maximum ranged widely, from 1 to 66 m. Estimates of Chl-a made from the water-leaving radiance measurements using the OC3M algorithm were highly skewed with most samples estimated as <1 mg/m**3 Chl-a, while data collected from whales had a broad distribution with Chl-a reaching >9 mg/m**3. The correlation between the satellite-derived and whale-derived Chl-a maxima was poor, a linear fit explained only 10% of the variance.

Neodymium isotope composition of planktonic foraminifera from ODP Hole 121-758A

This study presents neodymium isotope and elemental data for cleaned planktonic foraminifera from ODP site 758 in the southernmost reaches of the Bay of Bengal in the north-east Indian Ocean. Cleaning experiments using oxidative-reductive techniques suggest that diagenetic Fe-Mn oxyhydroxide coatings can be effectively removed, and that the measured Nd isotope composition reflects the composition of seawater from which the foraminiferal calcium carbonate was precipitated. Modern core-top Pulleniatina obliquiloculata and Globorotalia menardii give epsilon-Nd values of 310.12 +/- 0.16 and 310.28 +/- 0.16, respectively, indistinguishable from recent direct measurements of surface seawater in this area. A high-resolution Nd isotope record obtained from G. menardii for the past 150 kyr shows systematic variations (Delta epsilon-Nd = 3) on glacial-interglacial timescales. The timing of those variations shows a remarkable correspondence with the global oxygen isotope record, which suggests a process controlling the Nd isotope composition that responds in phase with global climate cycles. Palaeoclimate reconstruction indicates that during the last glacial maximum changes in monsoon circulation resulted in a reduction in rainfall over the Indian subcontinent, and a decrease in the flux of river water delivered to the Bay of Bengal. Thus, changes in the riverine input of Nd, a change in either flux or composition, most likely caused the isotope variations, although changes in dust source or local ocean circulation may have also played a role. These results clearly establish a link between climate change and variations in radiogenic isotopes in the oceans, and illustrate the potential of Nd isotopes in foraminifera for highresolution palaeoceanographic reconstruction.

Production, oxygen uptake, and sinking velocity of copepod fecal pellets originated from the central North Sea

Production, oxygen uptake, and sinking velocity of copepod fecal pellets egested by Temora longicornis were measured using a nanoflagellate (Rhodomonas sp.), a diatom (Thalassiosira weissflogii), or a coccolithophorid (Emiliania huxleyi) as food sources. Fecal pellet production varied between 0.8 pellets ind**-1 h**-1 and 3.8 pellets ind**-1 h**-1 and was significantly higher with T. weissflogii than with the other food sources. Average pellet size varied between 2.2 x 10**5 µm**3 and 10.0 x 10**5 µm**3. Using an oxygen microsensor, small-scale oxygen fluxes and microbial respiration rates were measured directly with a spatial resolution of 2 µm at the interface of copepod fecal pellets and the surrounding water. Averaged volume-specific respiration rates were 4.12 fmol O2 µm**-3 d**-1, 2.86 fmol O2 µm**-3 d**-1, and 0.73 fmol O2 µm**-3 d**-1 in pellets produced on Rhodomonas sp., T. weissflogii, and E. huxleyi, respectively. The average carbon-specific respiration rate was 0.15 d**-1 independent on diet (range: 0.08-0.21 d**-1). Because of ballasting of opal and calcite, sinking velocities were significantly higher for pellets produced on T. weissflogii (322 +- 169 m d**-1) and E. huxleyi (200 +- 93 m d**-1) than on Rhodomonas sp. (35 +- 29 m d**-1). Preservation of carbon was estimated to be approximately 10-fold higher in fecal pellets produced when T. longicornis was fed E. huxleyi or T. weissflogii rather than Rhodomonas sp. Our study directly demonstrates that ballast increases the sinking rate of freshly produced copepod fecal pellets but does not protect them from decomposition.

Chlorophyll and true-color maps of sediments from the intertidal zones in Winyah Bay, South Carolina, USA and List, Sylt, Germany

We used hyperspectral imaging to study short-term effects of bioturbation by lugworms (Arenicola marina) on the surficial biomass of microphytobenthos (MPB) in permeable marine sediments. Within days to weeks after the addition of a lugworm to a homogenized and recomposed sediment, the average surficial MPB biomass and its spatial heterogeneity were, respectively, 150 - 250% and 280% higher than in sediments without lugworms. The surficial sediment area impacted by a single medium-sized lugworm (~4 g wet weight) over this time-scale was at least 340 cm**2. While sediment reworking was the primary cause of the increased spatial heterogeneity, experiments with lugworm-mimics together with modeling showed that bioadvective porewater transport from depth to the sediment surface, as induced by the lugworm ventilating its burrow, was the main cause of the increased surficial MPB biomass. Although direct measurements of nutrient fluxes are lacking, our present data show that enhanced advective supply of nutrients from deeper sediment layers induced by faunal ventilation is an important mechanism that fuels high primary productivity at the surface of permeable sediments even though these systems are generally characterized by low standing stocks of nutrients and organic material.

Concentration of hydrogen sulfide in waters of the Black Sea in the layer of its coexistence with oxygen

About 100 parallel determinations of hydrogen sulfide by the volumetric and photometric methods were made in the layer of coexistence of oxygen with hydrogen sulfide (C layer). Thiosulfates were determined simultaneously. Regardless of locations of the stations, determinations by two methods coincided for the entire range of depths of occurrence of the C layer upper boundary. Within the C layer hydrogen sulfide readings obtained by these two independent methods agreed; thiosulfates were not found by direct measurements. Difference in the readings appears at the lower boundary of the C layer and below it, accompanied by appearance of thiosulfates. It is therefore concluded that it is correct to determine the upper boundary of the C layer by the iodometric method and to use concentration of hydrogen sulfide obtained by this method in the C layer to calculate rate of chemical oxidation of hydrogen sulfide in quasistationary processes.

Stable carbon isotope record of planktonic foraminifera and organic matter from sediment cores in the Gulf of Aqaba

The stable carbon isotopic composition of the planktonic foraminifera Globigerinoides sacculifer and G. ruber (white) and sedimentary organic matter from the northern Gulf of Aqaba have been investigated to estimate changes in delta13CDIC in surface waters during the last 1,000 years. The high sedimentation rates at the core sites (about 54 cm/Kyear) provide high temporal resolution (~10 years). Recent sediments at the top of the cores reflect conditions younger than 1950. The delta13C records of the planktonic foraminifera from three multicores display similar trends, showing a uniform and consistent pattern before the 1750s, and a gradual decrease of approximately 0.63? over the last two centuries. This decrease seems to track the decrease of delta13CDIC in surface waters, which is mainly caused by the increase of anthropogenic input of 13C-depleted CO2 into the atmosphere. Similarly, a trend towards lighter values of the carbon isotopic composition of sedimentary organic matter (delta13Corg) during the last 200 years supports the interpretation obtained from the planktonic foraminiferal delta13C. Furthermore, direct measurements of seawater show that delta13C of the dissolved inorganic carbon (DIC) in the northern Gulf of Aqaba has decreased by about 0.44 per mil during the period 1979-2000. The average annual decrease is 0.021 per mil, which is similar to that observed globally. The delta13C values of planktonic foraminifera combined with organic matter delta13C from marine sediments are good indicators for reconstructing past changes in atmospheric CO2 concentrations from the northern Gulf of Aqaba.

Bases metodológicas para un sistema de apoyo a la selección de especies en la restauración de la vegetación

RESUMEN El apoyo a la selección de especies a la restauración de la vegetación en España en los últimos 40 años se ha basado fundamentalmente en modelos de distribución de especies, también llamados modelos de nicho ecológico, que estiman la probabilidad de presencia de las especies en función de las condiciones del medio físico (clima, suelo, etc.). Con esta tesis se ha intentado contribuir a la mejora de la capacidad predictiva de los modelos introduciendo algunas propuestas metodológicas adaptadas a los datos disponibles actualmente en España y enfocadas al uso de los modelos en la selección de especies. No siempre se dispone de datos a una resolución espacial adecuada para la escala de los proyectos de restauración de la vegetación. Sin embrago es habitual contar con datos de baja resolución espacial para casi todas las especies vegetales presentes en España. Se propone un método de recalibración que actualiza un modelo de regresión logística de baja resolución espacial con una nueva muestra de alta resolución espacial. El método permite obtener predicciones de calidad aceptable con muestras relativamente pequeñas (25 presencias de la especie) frente a las muestras mucho mayores (más de 100 presencias) que requería una estrategia de modelización convencional que no usara el modelo previo. La selección del método estadístico puede influir decisivamente en la capacidad predictiva de los modelos y por esa razón la comparación de métodos ha recibido mucha atención en la última década. Los estudios previos consideraban a la regresión logística como un método inferior a técnicas más modernas como las de máxima entropía. Los resultados de la tesis demuestran que esa diferencia observada se debe a que los modelos de máxima entropía incluyen técnicas de regularización y la versión de la regresión logística usada en las comparaciones no. Una vez incorporada la regularización a la regresión logística usando penalización, las diferencias en cuanto a capacidad predictiva desaparecen. La regresión logística penalizada es, por tanto, una alternativa más para el ajuste de modelos de distribución de especies y está a la altura de los métodos modernos con mejor capacidad predictiva como los de máxima entropía. A menudo, los modelos de distribución de especies no incluyen variables relativas al suelo debido a que no es habitual que se disponga de mediciones directas de sus propiedades físicas o químicas. La incorporación de datos de baja resolución espacial proveniente de mapas de suelo nacionales o continentales podría ser una alternativa. Los resultados de esta tesis sugieren que los modelos de distribución de especies de alta resolución espacial mejoran de forma ligera pero estadísticamente significativa su capacidad predictiva cuando se incorporan variables relativas al suelo procedente de mapas de baja resolución espacial. La validación es una de las etapas fundamentales del desarrollo de cualquier modelo empírico como los modelos de distribución de especies. Lo habitual es validar los modelos evaluando su capacidad predictiva especie a especie, es decir, comparando en un conjunto de localidades la presencia o ausencia observada de la especie con las predicciones del modelo. Este tipo de evaluación no responde a una cuestión clave en la restauración de la vegetación ¿cuales son las n especies más idóneas para el lugar a restaurar? Se ha propuesto un método de evaluación de modelos adaptado a esta cuestión que consiste en estimar la capacidad de un conjunto de modelos para discriminar entre las especies presentes y ausentes de un lugar concreto. El método se ha aplicado con éxito a la validación de 188 modelos de distribución de especies leñosas orientados a la selección de especies para la restauración de la vegetación en España. Las mejoras metodológicas propuestas permiten mejorar la capacidad predictiva de los modelos de distribución de especies aplicados a la selección de especies en la restauración de la vegetación y también permiten ampliar el número de especies para las que se puede contar con un modelo que apoye la toma de decisiones. SUMMARY During the last 40 years, decision support tools for plant species selection in ecological restoration in Spain have been based on species distribution models (also called ecological niche models), that estimate the probability of occurrence of the species as a function of environmental predictors (e.g., climate, soil). In this Thesis some methodological improvements are proposed to contribute to a better predictive performance of such models, given the current data available in Spain and focusing in the application of the models to selection of species for ecological restoration. Fine grained species distribution data are required to train models to be used at the scale of the ecological restoration projects, but this kind of data are not always available for every species. On the other hand, coarse grained data are available for almost every species in Spain. A recalibration method is proposed that updates a coarse grained logistic regression model using a new fine grained updating sample. The method allows obtaining acceptable predictive performance with reasonably small updating sample (25 occurrences of the species), in contrast with the much larger samples (more than 100 occurrences) required for a conventional modeling approach that discards the coarse grained data. The choice of the statistical method may have a dramatic effect on model performance, therefore comparisons of methods have received much interest in the last decade. Previous studies have shown a poorer performance of the logistic regression compared to novel methods like maximum entropy models. The results of this Thesis show that the observed difference is caused by the fact that maximum entropy models include regularization techniques and the versions of logistic regression compared do not. Once regularization has been added to the logistic regression using a penalization procedure, the differences in model performance disappear. Therefore, penalized logistic regression may be considered one of the best performing methods to model species distributions. Usually, species distribution models do not consider soil related predictors because direct measurements of the chemical or physical properties are often lacking. The inclusion of coarse grained soil data from national or continental soil maps could be a reasonable alternative. The results of this Thesis suggest that the performance of the models slightly increase after including soil predictors form coarse grained soil maps. Model validation is a key stage of the development of empirical models, such as species distribution models. The usual way of validating is based on the evaluation of model performance for each species separately, i.e., comparing observed species presences or absence to predicted probabilities in a set of sites. This kind of evaluation is not informative for a common question in ecological restoration projects: which n species are the most suitable for the environment of the site to be restored? A method has been proposed to address this question that estimates the ability of a set of models to discriminate among present and absent species in a evaluation site. The method has been successfully applied to the validation of 188 species distribution models used to support decisions on species selection for ecological restoration in Spain. The proposed methodological approaches improve the predictive performance of the predictive models applied to species selection in ecological restoration and increase the number of species for which a model that supports decisions can be fitted.

Un método de diferencias finitas generalizadas para simular la actividad eléctrica de un tejido cardiaco

Una evolución del método de diferencias finitas ha sido el desarrollo del método de diferencias finitas generalizadas (MDFG) que se puede aplicar a mallas irregulares o nubes de puntos. En este método se emplea una expansión en serie de Taylor junto con una aproximación por mínimos cuadrados móviles (MCM). De ese modo, las fórmulas explícitas de diferencias para nubes irregulares de puntos se pueden obtener fácilmente usando el método de Cholesky. El MDFG-MCM es un método sin malla que emplea únicamente puntos. Una contribución de esta Tesis es la aplicación del MDFG-MCM al caso de la modelización de problemas anisótropos elípticos de conductividad eléctrica incluyendo el caso de tejidos reales cuando la dirección de las fibras no es fija, sino que varía a lo largo del tejido. En esta Tesis también se muestra la extensión del método de diferencias finitas generalizadas a la solución explícita de ecuaciones parabólicas anisótropas. El método explícito incluye la formulación de un límite de estabilidad para el caso de nubes irregulares de nodos que es fácilmente calculable. Además se presenta una nueva solución analítica para una ecuación parabólica anisótropa y el MDFG-MCM explícito se aplica al caso de problemas parabólicos anisótropos de conductividad eléctrica. La evidente dificultad de realizar mediciones directas en electrocardiología ha motivado un gran interés en la simulación numérica de modelos cardiacos. La contribución más importante de esta Tesis es la aplicación de un esquema explícito con el MDFG-MCM al caso de la modelización monodominio de problemas de conductividad eléctrica. En esta Tesis presentamos un algoritmo altamente eficiente, exacto y condicionalmente estable para resolver el modelo monodominio, que describe la actividad eléctrica del corazón. El modelo consiste en una ecuación en derivadas parciales parabólica anisótropa (EDP) que está acoplada con un sistema de ecuaciones diferenciales ordinarias (EDOs) que describen las reacciones electroquímicas en las células cardiacas. El sistema resultante es difícil de resolver numéricamente debido a su complejidad. Proponemos un método basado en una separación de operadores y un método sin malla para resolver la EDP junto a un método de Runge-Kutta para resolver el sistema de EDOs de la membrana y las corrientes iónicas. ABSTRACT An evolution of the method of finite differences has been the development of generalized finite difference (GFD) method that can be applied to irregular grids or clouds of points. In this method a Taylor series expansion is used together with a moving least squares (MLS) approximation. Then, the explicit difference formulae for irregular clouds of points can be easily obtained using a simple Cholesky method. The MLS-GFD is a mesh-free method using only points. A contribution of this Thesis is the application of the MLS-GFDM to the case of modelling elliptic anisotropic electrical conductivity problems including the case of real tissues when the fiber direction is not fixed, but varies throughout the tissue. In this Thesis the extension of the generalized finite difference method to the explicit solution of parabolic anisotropic equations is also given. The explicit method includes a stability limit formulated for the case of irregular clouds of nodes that can be easily calculated. Also a new analytical solution for homogeneous parabolic anisotropic equation has been presented and an explicit MLS- GFDM has been applied to the case of parabolic anisotropic electrical conductivity problems. The obvious difficulty of performing direct measurements in electrocardiology has motivated wide interest in the numerical simulation of cardiac models. The main contribution of this Thesis is the application of an explicit scheme based in the MLS-GFDM to the case of modelling monodomain electrical conductivity problems using operator splitting including the case of anisotropic real tissues. In this Thesis we present a highly efficient, accurate and conditionally stable algorithm to solve a monodomain model, which describes the electrical activity in the heart. The model consists of a parabolic anisotropic partial differential equation (PDE), which is coupled to systems of ordinary differential equations (ODEs) describing electrochemical reactions in the cardiac cells. The resulting system is challenging to solve numerically, because of its complexity. We propose a method based on operator splitting and a meshless method for solving the PDE together with a Runge-Kutta method for solving the system of ODE’s for the membrane and ionic currents.

Slow dimer dissociation of the TATA binding protein dictates the kinetics of DNA binding

The association of the TATA binding protein (TBP) to eukaryotic promoters is a possible rate-limiting step in gene expression. Slow promoter binding might be related to TBP’s ability to occlude its DNA binding domain through dimerization. Using a “pull-down” based assay, we find that TBP dimers dissociate slowly (t½ = 6–10 min), and thus present a formidable kinetic barrier to TATA binding. At 10 nM, TBP appears to exist as a mixed population of monomers and dimers. In this state, TATA binding displays burst kinetics that appears to reflect rapid binding of monomers and slow dissociation of dimers. The kinetics of the slow phase is in excellent agreement with direct measurements of the kinetics of dimer dissociation.

Submucosal gland secretions in airways from cystic fibrosis patients have normal [Na+] and pH but elevated viscosity

Fluid and macromolecule secretion by submucosal glands in mammalian airways is believed to be important in normal airway physiology and in the pathophysiology of cystic fibrosis (CF). An in situ fluorescence method was applied to measure the ionic composition and viscosity of freshly secreted fluid from airway glands. Fragments of human large airways obtained at the time of lung transplantation were mounted in a humidified perfusion chamber and the mucosal surface was covered by a thin layer of oil. Individual droplets of secreted fluid were microinjected with fluorescent indicators for measurement of [Na+], [Cl−], and pH by ratio imaging fluorescence microscopy and viscosity by fluorescence recovery after photobleaching. After carbachol stimulation, 0.1–0.5 μl of fluid accumulated in spherical droplets at gland orifices in ≈3–5 min. In gland fluid from normal human airways, [Na+] was 94 ± 8 mM, [Cl−] was 92 ± 12 mM, and pH was 6.97 ± 0.06 (SE, n = 7 humans, more than five glands studied per sample). Apparent fluid viscosity was 2.7 ± 0.3-fold greater than that of saline. Neither [Na+] nor pH differed in gland fluid from CF airways, but viscosity was significantly elevated by ≈2-fold compared to normal airways. These results represent the first direct measurements of ionic composition and viscosity in uncontaminated human gland secretions and indicate similar [Na+], [Cl−], and pH to that in the airway surface liquid. The elevated gland fluid viscosity in CF may be an important factor promoting bacterial colonization and airway disease.