Hydrologic landscape regionalisation using deductive classification and random forests

Landscape classification and hydrological regionalisation studies are being increasingly used in ecohydrology to aid in the management and research of aquatic resources. We present a methodology for classifying hydrologic landscapes based on spatial environmental variables by employing non-parametric statistics and hybrid image classification. Our approach differed from previous classifications which have required the use of an a priori spatial unit (e.g. a catchment) which necessarily results in the loss of variability that is known to exist within those units. The use of a simple statistical approach to identify an appropriate number of classes eliminated the need for large amounts of post-hoc testing with different number of groups, or the selection and justification of an arbitrary number. Using statistical clustering, we identified 23 distinct groups within our training dataset. The use of a hybrid classification employing random forests extended this statistical clustering to an area of approximately 228,000 km2 of south-eastern Australia without the need to rely on catchments, landscape units or stream sections. This extension resulted in a highly accurate regionalisation at both 30-m and 2.5-km resolution, and a less-accurate 10-km classification that would be more appropriate for use at a continental scale. A smaller case study, of an area covering 27,000 km2, demonstrated that the method preserved the intra- and inter-catchment variability that is known to exist in local hydrology, based on previous research. Preliminary analysis linking the regionalisation to streamflow indices is promising suggesting that the method could be used to predict streamflow behaviour in ungauged catchments. Our work therefore simplifies current classification frameworks that are becoming more popular in ecohydrology, while better retaining small-scale variability in hydrology, thus enabling future attempts to explain and visualise broad-scale hydrologic trends at the scale of catchments and continents.

PET imaging of tumours with a 64Cu labeled macrobicyclic cage amine ligand tethered to Tyr3-octreotate

The use of copper radioisotopes in cancer diagnosis and radionuclide therapy is possible using chelators that are capable of binding Cu(II) with sufficient stability in vivo to provide high tumour-to-background contrast. Here we report the design and synthesis of a new bifunctional chelator, 5-(8-methyl-3,6,10,13,16,19-hexaaza-bicyclo[6.6.6]icosan-1-ylamino)-5-oxopentanoic acid (MeCOSar), that forms copper complexes of exceptional stability by virtue of a cage amine (sarcophagine) ligand and a new conjugate referred to as SarTATE, obtained by the conjugation of MeCOSar to the tumour-targeting peptide Tyr(3)-octreotate. Radiolabeling of SarTATE with (64)Cu(II), a radioisotope suitable for positron emission tomography (PET), was fast (~20 min), easily performed at room temperature and consistently resulted in high radiochemical purity (>99%). In vitro and in vivo evaluation of (64)CuSarTATE demonstrated its high selectivity for tumour cells expressing somatostatin receptor 2 (sstr2). Biodistribution and PET imaging comparisons were made between (64)CuSarTATE and (64)Cu-labeled DOTA-Tyr(3)-octreotate ((64)CuDOTATATE). Both radiopharmaceuticals showed excellent uptake in sstr2-positive tumours at 2 h post-injection. While tumour uptake of (64)CuDOTATATE decreased significantly at 24 h, (64)CuSarTATE activity was retained, improving contrast at later time points. (64)CuSarTATE accumulated less than (64)CuDOTATATE in the non-target organs, liver and lungs. The uptake of (64)CuSarTATE in the kidneys was high at 2 h but showed significant clearance by 24 h. The new chemistry and pre-clinical evaluation presented here demonstrates that MeCOSar is a promising bifunctional chelator for Tyr(3)-octreotate that could be applied to a combined imaging and therapeutic regimen using a combination of (64)Cu- and (67)CuSarTATE complexes, owing to improved tumour-to-non-target organ ratios compared to (64)CuDOTATATE at longer time points.

NDVI With Artificial Neural Networks For SRTM Elevation Model Improvement – Hydrological Model Application

Digital elevation model (DEM) plays a substantial role in hydrological study, from understanding the catchment characteristics, setting up a hydrological model to mapping the flood risk in the region. Depending on the nature of study and its objectives, high resolution and reliable DEM is often desired to set up a sound hydrological model. However, such source of good DEM is not always available and it is generally high-priced. Obtained through radar based remote sensing, Shuttle Radar Topography Mission (SRTM) is a publicly available DEM with resolution of 92m outside US. It is a great source of DEM where no surveyed DEM is available. However, apart from the coarse resolution, SRTM suffers from inaccuracy especially on area with dense vegetation coverage due to the limitation of radar signals not penetrating through canopy. This will lead to the improper setup of the model as well as the erroneous mapping of flood risk. This paper attempts on improving SRTM dataset, using Normalised Difference Vegetation Index (NDVI), derived from Visible Red and Near Infra-Red band obtained from Landsat with resolution of 30m, and Artificial Neural Networks (ANN). The assessment of the improvement and the applicability of this method in hydrology would be highlighted and discussed.

Performance Evaluation Of SMOS Soil Moisture Retrieval Parameters For Hydrological Application

Microwave remote sensing has high potential for soil moisture retrieval. However, the efficient retrieval of soil moisture depends on optimally choosing the soil moisture retrieval parameters. In this study first the initial evaluation of SMOS L2 product is performed and then four approaches regarding soil moisture retrieval from SMOS brightness temperature are reported. The radiative transfer equation based tau-omega rationale is used in this study for the soil moisture retrievals. The single channel algorithms (SCA) using H polarisation is implemented with modifications, which includes the effective temperatures simulated from ECMWF (downscaled using WRF-NOAH Land Surface Model (LSM)) and MODIS. The retrieved soil moisture is then utilized for soil moisture deficit (SMD) estimation using empirical relationships with Probability Distributed Model based SMD as a benchmark. The square of correlation during the calibration indicates a value of R2 =0.359 for approach 4 (WRF-NOAH LSM based LST with optimized roughness parameters) followed by the approach 2 (optimized roughness parameters and MODIS based LST) (R2 =0.293), approach 3 (WRF-NOAH LSM based LST with no optimization) (R2 =0.267) and approach 1(MODIS based LST with no optimization) (R2 =0.163). Similarly, during the validation a highest performance is reported by approach 4. The other approaches are also following a similar trend as calibration. All the performances are depicted through Taylor diagram which indicates that the H polarisation using ECMWF based LST is giving a better performance for SMD estimation than the original SMOS L2 products at a catchment scale.

Incorporação de incertezas associadas à predição de modelos hidrológicos aos instrumentos de gestão em recursos hídricos

Water management has in the watershed plans an important tool to plan the territory and adjust the activities develop over it to the natural resources availability. The incorporation of uncertainty analysis associated with hydrological modelling predictions is a manner to simulate scenarios and work with chances and probabilities that certain events happens inside these plans. Using stochastic methods is possible to consider uncertainty from estimations and even model it. Stochastic methods developed considerably during the last 30 years, but its applications to real-world problems have been limited, and did not turn into routine in hydrology. This paper brings an overview from eminent hydrologists about this subject and discuss the Brazilian and Paulista situation in the scope of groundwater monitoring.

Non-Markovian stochastic processes and their applications: from anomalous diffusion to time series analysis

This work provides a forward step in the study and comprehension of the relationships between stochastic processes and a certain class of integral-partial differential equation, which can be used in order to model anomalous diffusion and transport in statistical physics. In the first part, we brought the reader through the fundamental notions of probability and stochastic processes, stochastic integration and stochastic differential equations as well. In particular, within the study of H-sssi processes, we focused on fractional Brownian motion (fBm) and its discrete-time increment process, the fractional Gaussian noise (fGn), which provide examples of non-Markovian Gaussian processes. The fGn, together with stationary FARIMA processes, is widely used in the modeling and estimation of long-memory, or long-range dependence (LRD). Time series manifesting long-range dependence, are often observed in nature especially in physics, meteorology, climatology, but also in hydrology, geophysics, economy and many others. We deepely studied LRD, giving many real data examples, providing statistical analysis and introducing parametric methods of estimation. Then, we introduced the theory of fractional integrals and derivatives, which indeed turns out to be very appropriate for studying and modeling systems with long-memory properties. After having introduced the basics concepts, we provided many examples and applications. For instance, we investigated the relaxation equation with distributed order time-fractional derivatives, which describes models characterized by a strong memory component and can be used to model relaxation in complex systems, which deviates from the classical exponential Debye pattern. Then, we focused in the study of generalizations of the standard diffusion equation, by passing through the preliminary study of the fractional forward drift equation. Such generalizations have been obtained by using fractional integrals and derivatives of distributed orders. In order to find a connection between the anomalous diffusion described by these equations and the long-range dependence, we introduced and studied the generalized grey Brownian motion (ggBm), which is actually a parametric class of H-sssi processes, which have indeed marginal probability density function evolving in time according to a partial integro-differential equation of fractional type. The ggBm is of course Non-Markovian. All around the work, we have remarked many times that, starting from a master equation of a probability density function f(x,t), it is always possible to define an equivalence class of stochastic processes with the same marginal density function f(x,t). All these processes provide suitable stochastic models for the starting equation. Studying the ggBm, we just focused on a subclass made up of processes with stationary increments. The ggBm has been defined canonically in the so called grey noise space. However, we have been able to provide a characterization notwithstanding the underline probability space. We also pointed out that that the generalized grey Brownian motion is a direct generalization of a Gaussian process and in particular it generalizes Brownain motion and fractional Brownain motion as well. Finally, we introduced and analyzed a more general class of diffusion type equations related to certain non-Markovian stochastic processes. We started from the forward drift equation, which have been made non-local in time by the introduction of a suitable chosen memory kernel K(t). The resulting non-Markovian equation has been interpreted in a natural way as the evolution equation of the marginal density function of a random time process l(t). We then consider the subordinated process Y(t)=X(l(t)) where X(t) is a Markovian diffusion. The corresponding time-evolution of the marginal density function of Y(t) is governed by a non-Markovian Fokker-Planck equation which involves the same memory kernel K(t). We developed several applications and derived the exact solutions. Moreover, we considered different stochastic models for the given equations, providing path simulations.

Towards tailored radiopeptide therapy.

PURPOSE Somatostatin receptor-targeted radiopeptide therapy is commonly performed using single radioisotopes. We evaluated the benefits and harms of combining radioisotopes in radiopeptide therapy in patients with neuroendocrine tumor. METHODS Using multivariable-adjusted survival analyses and competing risk analyses we evaluated outcomes in patients with neuroendocrine tumor receiving (90)Y-DOTATOC, (177)Lu-DOTATOC or their combination. RESULTS (90)Y-DOTATOC plus (177)Lu-DOTATOC treatment was associated with longer survival than (90)Y-DOTATOC (66.1 vs. 47.5 months; n = 1,358; p < 0.001) or (177)Lu-DOTATOC alone (66.1 vs. 45.5 months; n = 390; p < 0.001). (177)Lu-DOTATOC was associated with longer survival than (90)Y-DOTATOC in patients with solitary lesions (HR 0.3, range 0.1 - 0.7; n = 153; p = 0.005), extrahepatic metastases (HR 0.5, range 0.3 - 0.9; n = 256; p = 0.029) and metastases with low uptake (HR 0.1, range 0.05 - 0.4; n = 113; p = 0.001). (90)Y-DOTATOC induced higher hematotoxicity rates than combined treatment (9.5% vs. 4.0%, p = 0.005) or (177)Lu-DOTATOC (9.5 vs. 1.4%, p = 0.002). Renal toxicity was similar among the treatments. CONCLUSIONS Using (90)Y and (177)Lu might facilitate tailoring radiopeptide therapy and improve survival in patients with neuroendocrine tumors.

(Table 1) Characteristics of snow pits on De Geerfonna and Vestfonna during the IPY 2007-2010

Extensive glaciological field measurements were carried out on the ice cap Vestfonna as well as on the minor ice body De Geerfonna (Nordaustlandet, Svalbard) within the framework of IPY Kinnvika. Field campaigns were conducted during the period 2007-2010 in spring (April/May) and summer (August). In this study we compile and present snow cover information obtained from 22 snow pits that were dug on Vestfonna during this period. Locations are along two transects on the northwestern, land terminating slope of the ice cap, on its central summit, Ahlmann Summit, and at a set of several other locations in the eastern and northern part of the ice cap. Snow-cover information acquired from four snow pits on adjacent De Geerfonna is also incorporated in this study. Field data are analysed regarding snow stratigraphy, snow density, snow hardness and snow temperature. Results reveal mean snow densities of around 400 kg/m**3 for the snowpack of Vestfonna with no apparent spatial or interannual variability. A distinctly higher value of more than 450 kg/m**3 was obtained for De Geerfonna. A spatial comparison of snow water equivalents above the previous end-of-summer surface serves for obtaining insights into the spatial distribution of snow accumulation across Vestfonna. Altitude was found to be the only significant spatial parameter for controlling snow accumulation across the ice cap.

Aspectos metodológicos para la evaluación de sistemas de tomografía por emisión de positrones empleando técnicas Montecarlo, protocolos estandarizados y diferentes trazadores

Esta tesis analiza los elementos que afectan a la evaluación del rendimiento dentro de la técnica de radiodiagnóstico mediante tomografía por emisión de positrones (PET), centrándose en escáneres preclínicos. Se exploran las posibilidades de los protocolos estándar de evaluación sobre los siguientes aspectos: su uso como herramienta para validar programas de simulación Montecarlo, como método para la comparación de escáneres y su validez en el estudio del efecto sobre la calidad de imagen al utilizar radioisótopos alternativos. Inicialmente se estudian los métodos de evaluación orientados a la validación de simulaciones PET, para ello se presenta el programa GAMOS como entorno de simulación y se muestran los resultados de su validación basada en el estándar NEMA NU 4-2008 para escáneres preclínicos. Esta validación se ha realizado mediante la comparación de los resultados simulados frente a adquisiciones reales en el equipo ClearPET, describiendo la metodología de evaluación y selección de los parámetros NEMA. En este apartado también se mencionan las aportaciones desarrolladas en GAMOS para aplicaciones PET, como la inclusión de herramientas para la reconstrucción de imágenes. Por otro lado, la evaluación NEMA del ClearPET es utilizada para comparar su rendimiento frente a otro escáner preclínico: el sistema rPET-1. Esto supone la primera caracterización NEMA NU 4 completa de ambos equipos; al mismo tiempo que se analiza cómo afectan las importantes diferencias de diseño entre ellos, especialmente el tamaño axial del campo de visión y la configuración de los detectores. El 68Ga es uno de los radioisótopos no convencionales en imagen PET que está experimentando un mayor desarrollo, sin embargo, presenta la desventaja del amplio rango o distancia recorrida por el positrón emitido. Además del rango del positrón, otra propiedad física característica de los radioisótopos PET que puede afectar a la imagen es la emisión de fotones gamma adicionales, tal como le ocurre al isótopo 48V. En esta tesis se evalúan dichos efectos mediante estudios de resolución espacial y calidad de imagen NEMA. Finalmente, se analiza el alcance del protocolo NEMA NU 4-2008 cuando se utiliza para este propósito, adaptándolo a tal fin y proponiendo posibles modificaciones. Abstract This thesis analyzes the factors affecting the performance evaluation in positron emission tomography (PET) imaging, focusing on preclinical scanners. It explores the possibilities of standard protocols of assessment on the following aspects: their use as tools to validate Monte Carlo simulation programs, their usefulness as a method for comparing scanners and their validity in the study of the effect of alternative radioisotopes on image quality. Initially we study the methods of performance evaluation oriented to validate PET simulations. For this we present the GAMOS program as a simulation framework and show the results of its validation based on the standard NEMA NU 4-2008 for preclinical PET scanners. This has been accomplished by comparing simulated results against experimental acquisitions in the ClearPET scanner, describing the methodology for the evaluation and selection of NEMA parameters. This section also mentions the contributions developed in GAMOS for PET applications, such as the inclusion of tools for image reconstruction. Furthermore, the evaluation of the ClearPET scanner is used to compare its performance against another preclinical scanner, specifically the rPET-1 system. This is the first complete NEMA NU 4 based characterization study of both systems. At the same time we analyze how do the significant design differences of these two systems, especially the size of the axial field of view and the detectors configuration affect their performance characteristics. 68Ga is one of the unconventional radioisotopes in PET imaging the use of which is currently significantly increasing; however, it presents the disadvantage of the long positron range (distance traveled by the emitted positron before annihilating with an electron). Besides the positron range, additional gamma photon emission is another physical property characteristic of PET radioisotopes that can affect the reconstructed image quality, as it happens to the isotope 48V. In this thesis we assess these effects through studies of spatial resolution and image quality. Finally, we analyze the scope of the NEMA NU 4-2008 to carry out such studies, adapting it and proposing possible modifications.

Flood frequency analysis by a bivariate model based on copulas

La adecuada estimación de avenidas de diseño asociadas a altos periodos de retorno es necesaria para el diseño y gestión de estructuras hidráulicas como presas. En la práctica, la estimación de estos cuantiles se realiza normalmente a través de análisis de frecuencia univariados, basados en su mayoría en el estudio de caudales punta. Sin embargo, la naturaleza de las avenidas es multivariada, siendo esencial tener en cuenta características representativas de las avenidas, tales como caudal punta, volumen y duración del hidrograma, con el fin de llevar a cabo un análisis apropiado; especialmente cuando el caudal de entrada se transforma en un caudal de salida diferente durante el proceso de laminación en un embalse o llanura de inundación. Los análisis de frecuencia de avenidas multivariados han sido tradicionalmente llevados a cabo mediante el uso de distribuciones bivariadas estándar con el fin de modelar variables correlacionadas. Sin embargo, su uso conlleva limitaciones como la necesidad de usar el mismo tipo de distribuciones marginales para todas las variables y la existencia de una relación de dependencia lineal entre ellas. Recientemente, el uso de cópulas se ha extendido en hidrología debido a sus beneficios en relación al contexto multivariado, permitiendo superar los inconvenientes de las técnicas tradicionales. Una copula es una función que representa la estructura de dependencia de las variables de estudio, y permite obtener la distribución de frecuencia multivariada de dichas variables mediante sus distribuciones marginales, sin importar el tipo de distribución marginal utilizada. La estimación de periodos de retorno multivariados, y por lo tanto, de cuantiles multivariados, también se facilita debido a la manera en la que las cópulas están formuladas. La presente tesis doctoral busca proporcionar metodologías que mejoren las técnicas tradicionales usadas por profesionales para estimar cuantiles de avenida más adecuados para el diseño y la gestión de presas, así como para la evaluación del riesgo de avenida, mediante análisis de frecuencia de avenidas bivariados basados en cópulas. Las variables consideradas para ello son el caudal punta y el volumen del hidrograma. Con el objetivo de llevar a cabo un estudio completo, la presente investigación abarca: (i) el análisis de frecuencia de avenidas local bivariado centrado en examinar y comparar los periodos de retorno teóricos basados en la probabilidad natural de ocurrencia de una avenida, con el periodo de retorno asociado al riesgo de sobrevertido de la presa bajo análisis, con el fin de proporcionar cuantiles en una estación de aforo determinada; (ii) la extensión del enfoque local al regional, proporcionando un procedimiento completo para llevar a cabo un análisis de frecuencia de avenidas regional bivariado para proporcionar cuantiles en estaciones sin aforar o para mejorar la estimación de dichos cuantiles en estaciones aforadas; (iii) el uso de cópulas para investigar tendencias bivariadas en avenidas debido al aumento de los niveles de urbanización en una cuenca; y (iv) la extensión de series de avenida observadas mediante la combinación de los beneficios de un modelo basado en cópulas y de un modelo hidrometeorológico. Accurate design flood estimates associated with high return periods are necessary to design and manage hydraulic structures such as dams. In practice, the estimate of such quantiles is usually done via univariate flood frequency analyses, mostly based on the study of peak flows. Nevertheless, the nature of floods is multivariate, being essential to consider representative flood characteristics, such as flood peak, hydrograph volume and hydrograph duration to carry out an appropriate analysis; especially when the inflow peak is transformed into a different outflow peak during the routing process in a reservoir or floodplain. Multivariate flood frequency analyses have been traditionally performed by using standard bivariate distributions to model correlated variables, yet they entail some shortcomings such as the need of using the same kind of marginal distribution for all variables and the assumption of a linear dependence relation between them. Recently, the use of copulas has been extended in hydrology because of their benefits regarding dealing with the multivariate context, as they overcome the drawbacks of the traditional approach. A copula is a function that represents the dependence structure of the studied variables, and allows obtaining the multivariate frequency distribution of them by using their marginal distributions, regardless of the kind of marginal distributions considered. The estimate of multivariate return periods, and therefore multivariate quantiles, is also facilitated by the way in which copulas are formulated. The present doctoral thesis seeks to provide methodologies that improve traditional techniques used by practitioners, in order to estimate more appropriate flood quantiles for dam design, dam management and flood risk assessment, through bivariate flood frequency analyses based on the copula approach. The flood variables considered for that goal are peak flow and hydrograph volume. In order to accomplish a complete study, the present research addresses: (i) a bivariate local flood frequency analysis focused on examining and comparing theoretical return periods based on the natural probability of occurrence of a flood, with the return period associated with the risk of dam overtopping, to estimate quantiles at a given gauged site; (ii) the extension of the local to the regional approach, supplying a complete procedure for performing a bivariate regional flood frequency analysis to either estimate quantiles at ungauged sites or improve at-site estimates at gauged sites; (iii) the use of copulas to investigate bivariate flood trends due to increasing urbanisation levels in a catchment; and (iv) the extension of observed flood series by combining the benefits of a copula-based model and a hydro-meteorological model.

Stable carbon isotope ratios of alkane of IODP Hole 302-M0004A

The Palaeocene/Eocene thermal maximum represents a period of rapid, extreme global warming approx ~55 million years ago, superimposed on an already warm world (Zachos et al., 2003, doi:10.1126/science.1090110; Bowen et al., 2004, doi:10.1038/nature03115; Thomas et al., 2002, doi:10.1130/0091-7613(2002)030<1067:WTFFTF>2.0.CO;2). This warming is associated with a severe shoaling of the ocean calcite compensation depth **4 and a >2.5 per mil negative carbon isotope excursion in marine and soil carbonates (Zachos et al., 2003, doi:10.1126/science.1090110; Bowen et al., 2004, doi:10.1038/nature03115; Thomas et al., 2002, doi:10.1130/0091-7613(2002)030<1067:WTFFTF>2.0.CO;2; Zachos et al., doi:10.1126/science.1109004). Together these observations indicate a massive release of 13C-depleted carbon (Zachos et al., doi:10.1126/science.1109004) and greenhouse-gas-induced warming. Recently, sediments were recovered from the central Arctic Ocean (Backman et al., 2006, doi:10.2204/iodp.proc.302.2006), providing the first opportunity to evaluate the environmental response at the North Pole at this time. Here we present stable hydrogen and carbon isotope measurements of terrestrial-plant- and aquatic-derived n-alkanes that record changes in hydrology, including surface water salinity and precipitation, and the global carbon cycle. Hydrogen isotope records are interpreted as documenting decreased rainout during moisture transport from lower latitudes and increased moisture delivery to the Arctic at the onset of the Palaeocene/Eocene thermal maximum, consistent with predictions of poleward storm track migrations during global warming (Backman et al., 2006, doi:10.2204/iodp.proc.302.2006). The terrestrial-plant carbon isotope excursion (about ~4.5 to ~6 per mil) is substantially larger than those of marine carbonates. Previously, this offset was explained by the physiological response of plants to increases in surface humidity (Bowen et al., 2004, doi:10.1038/nature03115). But this mechanism is not an effective explanation in this wet Arctic setting, leading us to hypothesize that the true magnitude of the excursion - and associated carbon input - was greater than originally surmised. Greater carbon release and strong hydrological cycle feedbacks may help explain the maintenance of this unprecedented warmth.of this unprecedented warmth.

Insect control by radiation and radiosotopes : a literature search /

This bibliography contains one hundred and one selected references on the uses of radiation and radioisotopes in entomology. These references were selected from the scientific literature published since 1961. An index to the insects studied is included.

Development of additional hazard assessment models /

"March 1977."

A system for using radar to record wave direction /

"September 1979."

Thyroid uptake calibration.

"Contract No. AT(40-1)-Gen-33."