An accurate method to correct atmospheric phase delay for InSAR with the ERA5 global atmospheric model

Differential SAR Interferometry (DInSAR) has proven its unprecedented ability and merits of monitoring ground deformation on a large scale with centimeter to millimeter accuracy. However, atmospheric artifacts due to spatial and temporal variations of the atmospheric state often affect the reliability and accuracy of its results. The commonly-known Atmospheric Phase Screen (APS) appears in the interferograms as ghost fringes not related to either topography or deformation. Atmospheric artifact mitigation remains one of the biggest challenges to be addressed within the DInSAR community. State-of-the-art research works have revealed that atmospheric artifacts can be partially compensated with empirical models, point-wise GPS zenith path delay, and numerical weather prediction models. In this study, we implement an accurate and realistic computing strategy using atmospheric reanalysis ERA5 data to estimate atmospheric artifacts. With this approach, the Line-of-Sight (LOS) path along the satellite trajectory and the monitored points is considered, rather than estimating it from the zenith path delay. Compared with the zenith delay-based method, the key advantage is that it can avoid errors caused by any anisotropic atmospheric phenomena. The accurate method is validated with Sentinel-1 data in three different test sites: Tenerife island (Spain), Almería (Spain), and Crete island (Greece). The effectiveness and performance of the method to remove APS from interferograms is evaluated in the three test sites showing a great improvement with respect to the zenith-based approach.Peer ReviewedPostprint (published version

Geomorphic flood hazard mapping: from floodplain delineation to flood hazard characterization

Recent studies show how geomorphic descriptors, retrieved from digital elevation models (DEMs), can be used for flood hazard mapping. As they strictly depend on the accuracy of the input DEMs and reference flood hazard maps used for training, DEM-based flood hazard models may display severe inconsistencies. Our study shows the application of two advanced DEM-based models to a large study area, and presents two main innovative points. First, the delicate tasks of appropriately selecting the input DEM and flood hazard map are specifically addressed with newly defined methods. Second, the ability of DEM-based models to exploit their natural features to enhance flood hazard mapping over the study region is investigated. Our results show (a) the benefits of considering multiple geomorphic descriptors, (b) the potential of DEM-based models for completing the information of imperfect reference flood hazard maps, and (c) the advantages of continuous representation of hazard over binary flood maps

Radar satellite imagery for humanitarian response. Bridging the gap between technology and application

This work deals with radar satellite imagery and its potential to assist of humanitarian operations. As the number of displaced people annually increases, both hosting countries and relief organizations face new challenges which are often related to unclear situations and lack of information on the number and location of people in need, as well as their environments. It was demonstrated in numerous studies that methods of earth observation can deliver this important information for the management of crises, the organization of refugee camps, and the mapping of environmental resources and natural hazards. However, most of these studies make use of -high-resolution optical imagery, while the role of radar satellites is widely neglected. At the same time, radar sensors have characteristics which make them highly suitable for humanitarian response, their potential to capture images through cloud cover and at night in the first place. Consequently, they potentially allow quicker response in cases of emergencies than optical imagery. This work demonstrates the currently unused potential of radar imagery for the assistance of humanitarian operations by case studies which cover the information needs of specific emergency situations. They are thematically grouped into topics related to population, natural hazards and the environment. Furthermore, the case studies address different levels of scientific objectives: The main intention is the development of innovative techniques of digital image processing and geospatial analysis as an answer on the identified existing research gaps. For this reason, novel approaches are presented on the mapping of refugee camps and urban areas, the allocation of biomass and environmental impact assessment. Secondly, existing methods developed for radar imagery are applied, refined, or adapted to specifically demonstrate their benefit in a humanitarian context. This is done for the monitoring of camp growth, the assessment of damages in cities affected by civil war, and the derivation of areas vulnerable to flooding or sea-surface changes. Lastly, to foster the integration of radar images into existing operational workflows of humanitarian data analysis, technically simple and easily-adaptable approaches are suggested for the mapping of rural areas for vaccination campaigns, the identification of changes within and around refugee camps, and the assessment of suitable locations for groundwater drillings. While the studies provide different levels of technical complexity and novelty, they all show that radar imagery can largely contribute to the provision of a variety of information which is required to make solid decisions and to effectively provide help in humanitarian operations. This work furthermore demonstrates that radar images are more than just an alternative image source for areas heavily affected by cloud cover. In fact, what makes them valuable is their information content regarding the characteristics of surfaces, such as shape, orientation, roughness, size, height, moisture, or conductivity. All these give decisive insights about man-made and natural environments in emergency situations and cannot be provided by optical images Finally, the findings of the case studies are put into a larger context, discussing the observed potential and limitations of the presented approaches. The major challenges are summarized which need be addressed to make radar imagery more useful in humanitarian operations in the context of upcoming technical developments. New radar satellites and technological progress in the fields of machine learning and cloud computing will bring new opportunities. At the same time, this work demonstrated the large need for further research, as well as for the collaboration and transfer of knowledge and experiences between scientists, users and relief workers in the field. It is the first extensive scientific compilation of this topic and the first step for a sustainable integration of radar imagery into operational frameworks to assist humanitarian work and to contribute to a more efficient provision of help to those in need.Die vorliegende Arbeit beschäftigt sich mit bildgebenden Radarsatelliten und ihrem potenziellen Beitrag zur Unterstützung humanitärer Einsätze. Die jährlich zunehmende Zahl an vertriebenen oder geflüchteten Menschen stellt sowohl Aufnahmeländer als auch humanitäre Organisationen vor große Herausforderungen, da sie oft mit unübersichtlichen Verhältnissen konfrontiert sind. Effektives Krisenmanagement, die Planung und Versorgung von Flüchtlingslagern, sowie der Schutz der betroffenen Menschen erfordern jedoch verlässliche Angaben über Anzahl und Aufenthaltsort der Geflüchteten und ihrer natürlichen Umwelt. Die Bereitstellung dieser Informationen durch Satellitenbilder wurde bereits in zahlreichen Studien aufgezeigt. Sie beruhen in der Regel auf hochaufgelösten optischen Aufnahmen, während bildgebende Radarsatelliten bisher kaum Anwendung finden. Dabei verfügen gerade Radarsatelliten über Eigenschaften, die hilfreich für humanitäre Einsätze sein können, allen voran ihre Unabhängigkeit von Bewölkung oder Tageslicht. Dadurch ermöglichen sie in Krisenfällen verglichen mit optischen Satelliten eine schnellere Reaktion. Diese Arbeit zeigt das derzeit noch ungenutzte Potenzial von Radardaten zur Unterstützung humanitärer Arbeit anhand von Fallstudien auf, in denen konkrete Informationen für ausgewählte Krisensituationen bereitgestellt werden. Sie sind in die Themenbereiche Bevölkerung, Naturgefahren und Ressourcen aufgeteilt, adressieren jedoch unterschiedliche wissenschaftliche Ansprüche: Der Hauptfokus der Arbeit liegt auf der Entwicklung von innovativen Methoden zur Verarbeitung von Radarbildern und räumlichen Daten als Antwort auf den identifizierten Forschungsbedarf in diesem Gebiet. Dies wird anhand der Kartierung von Flüchtlingslagern zur Abschätzung ihrer Bevölkerung, zur Bestimmung von Biomasse, sowie zur Ermittlung des Umwelteinflusses von Flüchtlingslagern aufgezeigt. Darüber hinaus werden existierende oder erprobte Ansätze für die Anwendung im humanitären Kontext angepasst oder weiterentwickelt. Dies erfolgt im Rahmen von Fallstudien zur Dynamik von Flüchtlingslagern, zur Ermittlung von Schäden an Gebäuden in Kriegsgebieten, sowie zur Erkennung von Risiken durch Überflutung. Zuletzt soll die Integration von Radardaten in bereits existierende Abläufe oder Arbeitsroutinen in der humanitären Hilfe anhand technisch vergleichsweise einfacher Ansätze vorgestellt und angeregt werden. Als Beispiele dienen hier die radargestützte Kartierung von entlegenen Gebieten zur Unterstützung von Impfkampagnen, die Identifizierung von Veränderungen in Flüchtlingslagern, sowie die Auswahl geeigneter Standorte zur Grundwasserentnahme. Obwohl sich die Fallstudien hinsichtlich ihres Innovations- und Komplexitätsgrads unterscheiden, zeigen sie alle den Mehrwert von Radardaten für die Bereitstellung von Informationen, um schnelle und fundierte Planungsentscheidungen zu unterstützen. Darüber hinaus wird in dieser Arbeit deutlich, dass Radardaten für humanitäre Zwecke mehr als nur eine Alternative in stark bewölkten Gebieten sind. Durch ihren Informationsgehalt zur Beschaffenheit von Oberflächen, beispielsweise hinsichtlich ihrer Rauigkeit, Feuchte, Form, Größe oder Höhe, sind sie optischen Daten überlegen und daher für viele Anwendungsbereiche im Kontext humanitärer Arbeit besonders. Die in den Fallstudien gewonnenen Erkenntnisse werden abschließend vor dem Hintergrund von Vor- und Nachteilen von Radardaten, sowie hinsichtlich zukünftiger Entwicklungen und Herausforderungen diskutiert. So versprechen neue Radarsatelliten und technologische Fortschritte im Bereich der Datenverarbeitung großes Potenzial. Gleichzeitig unterstreicht die Arbeit einen großen Bedarf an weiterer Forschung, sowie an Austausch und Zusammenarbeit zwischen Wissenschaftlern, Anwendern und Einsatzkräften vor Ort. Die vorliegende Arbeit ist die erste umfassende Darstellung und wissenschaftliche Aufarbeitung dieses Themenkomplexes. Sie soll als Grundstein für eine langfristige Integration von Radardaten in operationelle Abläufe dienen, um humanitäre Arbeit zu unterstützen und eine wirksame Hilfe für Menschen in Not ermöglichen

Modelling the geographical distribution of soil-transmitted helminth infections in Bolivia

The prevalence of infection with the three common soil-transmitted helminths (i.e. Ascaris lumbricoides, Trichuris trichiura, and hookworm) in Bolivia is among the highest in Latin America. However, the spatial distribution and burden of soil-transmitted helminthiasis are poorly documented.; We analysed historical survey data using Bayesian geostatistical models to identify determinants of the distribution of soil-transmitted helminth infections, predict the geographical distribution of infection risk, and assess treatment needs and costs in the frame of preventive chemotherapy. Rigorous geostatistical variable selection identified the most important predictors of A. lumbricoides, T. trichiura, and hookworm transmission.; Results show that precipitation during the wettest quarter above 400 mm favours the distribution of A. lumbricoides. Altitude has a negative effect on T. trichiura. Hookworm is sensitive to temperature during the coldest month. We estimate that 38.0%, 19.3%, and 11.4% of the Bolivian population is infected with A. lumbricoides, T. trichiura, and hookworm, respectively. Assuming independence of the three infections, 48.4% of the population is infected with any soil-transmitted helminth. Empirical-based estimates, according to treatment recommendations by the World Health Organization, suggest a total of 2.9 million annualised treatments for the control of soil-transmitted helminthiasis in Bolivia.; We provide estimates of soil-transmitted helminth infections in Bolivia based on high-resolution spatial prediction and an innovative variable selection approach. However, the scarcity of the data suggests that a national survey is required for more accurate mapping that will govern spatial targeting of soil-transmitted helminthiasis control

Application of random forest classification and remotely sensed data in geological mapping on the Jebel Meloussi area (Tunisia)

Remotely sensed data such as satellite photos and radar images can be used to produce geological maps on arid regions, where the vegetation coverage does not have a significant effect. In central Tunisia, the Jebel Meloussi area has unique geological features and characteristic morphology (i.e. flat areas with dune fields in contrast with hills of folded and eroded stratigraphic sequences), which makes it an ideal area for testing new methods of automatic terrain classification. For this, data from the Sentinel 2 satellite sensor and the SRTM-based MERIT DEM (digital elevation model) were used in the present study. Using R scripts and the random forest classification method, modelling was performed on four lithological variables-derived from the different bands of the Sentinel 2 images-and two morphometric parameters for the area of the 1:50,000 geological map sheet no. 103. The four lithological variables were chosen to highlight the iron-bearing minerals since the spectral parameters of the Sentinel 2 sensors are especially useful for this purpose. The training areas of the classification were selected on the geological map. The results of the modelling identified Eocene and Cretaceous evaporite-bearing sedimentary series (such as the Jebs and the Bouhedma Formations) with the highest producer accuracy (> 60% of the predicted pixels match with the map). The pyritic argillites of the Sidi Khalif Formation were also recognized with the same accuracy, and the Quaternary sebhkas and dunes were also well predicted. The study concludes that the classification-based geological map is useful for field geologist prior to field surveys

Innovative techniques for the hydraulic and hydrological variables assessment

This Thesis focues on two main research topics: (1) the use of innovative techniques for the evaluation of main hydraulic variables of natural rivers (e.g. river bathymetry, discharge, water level) and (2) the monitoring and hydrological modelling of Monate Lake (Varese, Italy)

Long-term monitoring of geodynamic surface deformation using SAR interferometry

Thesis (Ph.D.) University of Alaska Fairbanks, 2014Synthetic Aperture Radar Interferometry (InSAR) is a powerful tool to measure surface deformation and is well suited for surveying active volcanoes using historical and existing satellites. However, the value and applicability of InSAR for geodynamic monitoring problems is limited by the influence of temporal decorrelation and electromagnetic path delay variations in the atmosphere, both of which reduce the sensitivity and accuracy of the technique. The aim of this PhD thesis research is: how to optimize the quantity and quality of deformation signals extracted from InSAR stacks that contain only a low number of images in order to facilitate volcano monitoring and the study of their geophysical signatures. In particular, the focus is on methods of mitigating atmospheric artifacts in interferograms by combining time-series InSAR techniques and external atmospheric delay maps derived by Numerical Weather Prediction (NWP) models. In the first chapter of the thesis, the potential of the NWP Weather Research & Forecasting (WRF) model for InSAR data correction has been studied extensively. Forecasted atmospheric delays derived from operational High Resolution Rapid Refresh for the Alaska region (HRRRAK) products have been compared to radiosonding measurements in the first chapter. The result suggests that the HRRR-AK operational products are a good data source for correcting atmospheric delays in spaceborne geodetic radar observations, if the geophysical signal to be observed is larger than 20 mm. In the second chapter, an advanced method for integrating NWP products into the time series InSAR workflow is developed. The efficiency of the algorithm is tested via simulated data experiments, which demonstrate the method outperforms other more conventional methods. In Chapter 3, a geophysical case study is performed by applying the developed algorithm to the active volcanoes of Unimak Island Alaska (Westdahl, Fisher and Shishaldin) for long term volcano deformation monitoring. The volcano source location at Westdahl is determined to be approx. 7 km below sea level and approx. 3.5 km north of the Westdahl peak. This study demonstrates that Fisher caldera has had continuous subsidence over more than 10 years and there is no evident deformation signal around Shishaldin peak.Chapter 1. Performance of the High Resolution Atmospheric Model HRRR-AK for Correcting Geodetic Observations from Spaceborne Radars -- Chapter 2. Robust atmospheric filtering of InSAR data based on numerical weather prediction models -- Chapter 3. Subtle motion long term monitoring of Unimak Island from 2003 to 2010 by advanced time series SAR interferometry -- Chapter 4. Conclusion and future work

Development of novel radiotracers as tools for imaging the human brain

Introduction: Brain imaging using single photon emission computed tomography (SPECT) or positron emission tomography (PET) can be used to study the processes underlying neurological and psychiatric disorders. In addition, in vivo brain imaging using SPECT or PET may provide new approaches for drug target identification, pre-clinical testing and occupancy studies, and therefore improve drug discovery. The utility of in vivo brain imaging using SPECT or PET relies on the ability of different radiotracers (typically organic compounds labelled with radionuclides) to bind to a wide variety of targets, including receptors, transporters and enzymes. Therefore the development of novel radiotracers for in vivo brain imaging using SPECT of PET is of vital importance. This thesis is focused on the process of developing novel radiotracers as tools for imaging the human brain, where the radiotracer discovery and development pipeline is discussed and each step prior to clinical trials investigated. Radiotracer discovery: Previously, discovery of novel brain radiotracers has largely relied on simplistic screening tools. Improved selection methods at the early stages of radiotracer discovery and an increased understanding of the relationships between in vitro physicochemical and in vivo radiotracer properties are needed. This thesis investigated if high performance liquid chromatography (HPLC) methodologies could provide criteria for lead candidate selection by comparing HPLC measurements with radiotracer properties in humans. In this study, ten molecules, previously used as radiotracers in humans, were analysed to obtain the following measures: partition coefficient (Log P); permeability (Pm); percentage of plasma protein binding (%PPB); and membrane partition coefficient (Km). Relationships between brain entry measurements (Log P, Pm and %PPB) and in vivo brain percentage injected dose (%ID); and Km and specific binding in vivo (BPND) were investigated. Results showed that HPLC measurements of Pm, %PPB and Km were potentially useful in predicting in vivo performance and hence allow evaluation and ranking of compound libraries for the selection of lead radiotracer candidates at early stages of radiotracer discovery. The HPLC tool developed provides information on in vivo non-specific binding and binding potential that is not possible using conventional screening methods. Another important finding reported in this thesis is that Log P should not be relied on as a predictor of brain entry. The HPLC tool developed, together with competition binding assays, was used to characterise a newly synthesised library of compounds for imaging of the translocator protein (TSPO) in brain using SPECT. Results showed that compound LS 1 was the most likely to succeed within the library investigated, but the high %PPB observed for LS 1 suggested novel compounds with improved %PPB were needed. Thus, a novel library of compounds for imaging of TSPO in brain using SPECT is currently been developed for future testing using the HPLC tool developed here and competition binding assays. Pre-clinical research: radiotracers for imaging the noradrenaline transporter (NAT) in brain using SPECT. In this thesis, NKJ64, a novel iodinated analogue of reboxetine, was successfully radiolabelled via electrophilic iododestannylation and evaluated as a potential SPECT radiotracer for imaging the NAT in brain using rodents and non-human primates. Biological evaluation of the novel radiotracer, 123/125I-NKJ64, in rodents included: in vitro ligand binding assays; in vitro and ex vivo autoradiography; in vivo biodistribution studies and ex vivo pharmacological blocking studies. In rats, 123/125I-NKJ64 displayed saturable binding with nanomolar affinity for the NAT in cortical homogenates, regional distribution consistent with the known density of NAT in the rodent brain and high maximum brain uptake of around 2.93 % of the injected dose. The specific: non-specific ratio (locus coeruleus:caudate putamen) of 123I-NKJ64 uptake was 2.8 at 30 minutes post intravenous injection and prior administration of reboxetine significantly reduced the accumulation of 123I-NKJ64 in the locus coeruleus (> 50% reduction). Data obtained using rodents indicated that further evaluation of 123I-NKJ64 in non-human primates was needed to determine its utility as a SPECT radiotracer for imaging of NAT in brain. Consequently, in vivo kinetic modelling studies using SPECT imaging with 123I-NKJ64 and two baboons were carried out to determine 123I-NKJ64 brain binding kinetics, brain distribution and plasma metabolism in non-human primates. Even though a high brain uptake of around 3.0% of the injected dose was determined, the high non-specific binding observed throughout the brain, a low binding potential (BPND<2) in NAT rich regions and a brain distribution that was inconsistent with the known NAT distribution in non-human primate brain precludes the translation of 123I-NKJ64 into humans. Another NAT radiotracer, 123I-INER, developed by Tamagnan and colleagues at Yale University and Institute for Degenerative Disorders, New Haven, USA, was also investigated as part of this thesis. Kinetic modelling analysis of 123I-INER in baboon brain was investigated for different models, namely invasive and reference tissue models. Bolus plus constant infusion experiments with displacement at equilibrium using six different doses of atomoxetine and four different doses of reboxetine were carried out in several baboons to obtain occupancy measurements as a function of injected dose (mg/kg) for the two NAT selective drugs. Results showed that reference tissue models were able to determine BPND values of 123I-INER in different brain regions. In addition the volume of distribution could be determined by dividing concentration in tissue by the concentration in venous blood at 3 hours post-injection. After administration of atomoxetine or reboxetine, dose-dependent occupancy was observed in brain regions known to contain high densities of NATs. Results supported the translation of 123I-INER into humans studies, despite the slow kinetics determined over the imaging period. Pharmacokinetic properties of 123I-INER described in this thesis may be used to simplify future data acquisition and image processing. Conclusion In conclusion, this thesis reported: (1) the development of novel radiotracers for brain imaging, namely NAT and TSPO; and (2) the development of a new methodology for aiding lead molecule identification at early stages of radiotracer discovery (i.e. prior to radiolabelling). In addition, an overview of radiotracer discovery and development process is provided in a single document, with a focus on brain radiotracers

The use of remotely sensed data for forest biomass monitoring : a case of forest sites in north-eastern Armenia

Dissertation submitted in partial fulfilment of the requirements for the Degree of Master of Science in Geospatial TechnologiesIn recent years there has been an increasing interest in the use of synthetic aperture radar (SAR) data and geospatial technologies for environmental monitoring․ Particularly, forest biomass evaluation was of high importance, as forests have a crucial role in global carbon emission. Within this study we evaluate the use of Sentinel 1 C-band multitemporal SAR data with combination of Alos Palsar L-band SAR and Sentinel 2 multispectral remote sensing (RS) data for mapping forest aboveground biomass (AGB) of dry subtropical forests in mountainous areas. Field observation from National Forest Inventory was used as a ground truth data. As the SAR data suffers greatly by the complex topography, a simple approach of aspect and slope information as forestry ancillary data was implemented directly in the regression model for the first time to mitigate the topography effect on radar backscattering value․ Dense time-series analysis allowed us to overcome the SAR saturation by the forest phenology and select the optimal C-band scene. Image texture measures of SAR data has been strongly related to the biomass distribution and has robustly contributed to the prediction․ Multilinear Stepwise Regression allowed to select and evaluate the most relevant variables for AGB. The prediction model combining RS with ancillary data explained the 62 % of variance with root-mean-square error of 56.6 t ha¯¹. The study also reveals that C-band SAR data on forest biomass prediction is limited due to their short wavelength. Further, the mountainous condition is a major constraint for AGB estimation. Additionally, this research demonstrates a positive outcome in forest AGB prediction with freely accessible RS data