Tsunami simulation and detection using global navigation satellite system reflectometry

In this thesis, research for tsunami remote sensing using the Global Navigation Satellite System-Reflectometry (GNSS-R) delay-Doppler maps (DDMs) is presented. Firstly, a process for simulating GNSS-R DDMs of a tsunami-dominated sea sur- face is described. In this method, the bistatic scattering Zavorotny-Voronovich (Z-V) model, the sea surface mean square slope model of Cox and Munk, and the tsunami- induced wind perturbation model are employed. The feasibility of the Cox and Munk model under a tsunami scenario is examined by comparing the Cox and Munk model- based scattering coefficient with the Jason-1 measurement. A good consistency be- tween these two results is obtained with a correlation coefficient of 0.93. After con- firming the applicability of the Cox and Munk model for a tsunami-dominated sea, this work provides the simulations of the scattering coefficient distribution and the corresponding DDMs of a fixed region of interest before and during the tsunami. Fur- thermore, by subtracting the simulation results that are free of tsunami from those with presence of tsunami, the tsunami-induced variations in scattering coefficients and DDMs can be clearly observed. Secondly, a scheme to detect tsunamis and estimate tsunami parameters from such tsunami-dominant sea surface DDMs is developed. As a first step, a procedure to de- termine tsunami-induced sea surface height anomalies (SSHAs) from DDMs is demon- strated and a tsunami detection precept is proposed. Subsequently, the tsunami parameters (wave amplitude, direction and speed of propagation, wavelength, and the tsunami source location) are estimated based upon the detected tsunami-induced SSHAs. In application, the sea surface scattering coefficients are unambiguously re- trieved by employing the spatial integration approach (SIA) and the dual-antenna technique. Next, the effective wind speed distribution can be restored from the scat- tering coefficients. Assuming all DDMs are of a tsunami-dominated sea surface, the tsunami-induced SSHAs can be derived with the knowledge of background wind speed distribution. In addition, the SSHA distribution resulting from the tsunami-free DDM (which is supposed to be zero) is considered as an error map introduced during the overall retrieving stage and is utilized to mitigate such errors from influencing sub- sequent SSHA results. In particular, a tsunami detection procedure is conducted to judge the SSHAs to be truly tsunami-induced or not through a fitting process, which makes it possible to decrease the false alarm. After this step, tsunami parameter estimation is proceeded based upon the fitted results in the former tsunami detec- tion procedure. Moreover, an additional method is proposed for estimating tsunami propagation velocity and is believed to be more desirable in real-world scenarios. The above-mentioned tsunami-dominated sea surface DDM simulation, tsunami detection precept and parameter estimation have been tested with simulated data based on the 2004 Sumatra-Andaman tsunami event.

GPS linked photos of benthic cover transect surveys in Heron Reef

Underwater photo-transect surveys were conducted on September 23-27, 2007 at different sections of the reef flat, reef crest and reef slope in Heron Reef. This survey was done by swimming along pre-defined transect sites and taking a picture of the bottom substrate parallel to the bottom at constant vertical distance (30cm) every two to three metres. A total of 3,586 benthic photos were taken. A floating GPS setup connected to the swimmer/diver by a line enabled recording of coordinates of transect surveys. Approximation of the coordinates for each benthic photo was based on the photo timestamp and GPS coordinate time stamp, using GPS Photo Link Software. Coordinates of each photo were interpolated by finding the the gps coordinates that were logged at a set time before and after the photo was captured. The output of this process was an ArcMap point shapefile, a Google Earth KML file and a thumbnail of each benthic photo taken. The data in the ArcMap shapefile and in the Google Earth KML file consisted of the approximated coordinate of each benthic photo taken during the survey. Using the GPS Photo Link extension within the ArcMap environment, opening the ArcMap shapefile will enable thumbnail to be displayed on the associated benthic cover photo whenever hovering with the mouse over a point on the transect. By downloading the GPSPhotoLink software from the www.geospatialexperts.com, and installing it as a trial version the ArcMap exstension will be installed in the ArcMap environment.

Circumarctic wetland dataset based on ENVISAT ASAR Global Monitoring Mode with links to geotiff image and NetCDF files (2005-01 to 2011-12)

A circumpolar representative and consistent wetland map is required for a range of applications ranging from upscaling of carbon fluxes and pools to climate modelling and wildlife habitat assessment. Currently available data sets lack sufficient accuracy and/or thematic detail in many regions of the Arctic. Synthetic aperture radar (SAR) data from satellites have already been shown to be suitable for wetland mapping. Envisat Advanced SAR (ASAR) provides global medium-resolution data which are examined with particular focus on spatial wetness patterns in this study. It was found that winter minimum backscatter values as well as their differences to summer minimum values reflect vegetation physiognomy units of certain wetness regimes. Low winter backscatter values are mostly found in areas vegetated by plant communities typically for wet regions in the tundra biome, due to low roughness and low volume scattering caused by the predominant vegetation. Summer to winter difference backscatter values, which in contrast to the winter values depend almost solely on soil moisture content, show expected higher values for wet regions. While the approach using difference values would seem more reasonable in order to delineate wetness patterns considering its direct link to soil moisture, it was found that a classification of winter minimum backscatter values is more applicable in tundra regions due to its better separability into wetness classes. Previous approaches for wetland detection have investigated the impact of liquid water in the soil on backscatter conditions. In this study the absence of liquid water is utilized. Owing to a lack of comparable regional to circumpolar data with respect to thematic detail, a potential wetland map cannot directly be validated; however, one might claim the validity of such a product by comparison with vegetation maps, which hold some information on the wetness status of certain classes. It was shown that the Envisat ASAR-derived classes are related to wetland classes of conventional vegetation maps, indicating its applicability; 30% of the land area north of the treeline was identified as wetland while conventional maps recorded 1-7%.

Impact of protected forests on rural household fuel tree planting and deforestation: Chiro District, Ethiopia

This study aims at exploring the potential impact of forest protection intervention on rural households’ private fuel tree planting in Chiro district of eastern Ethiopia. The study results revealed a robust and significant positive impact of the intervention on farmers’ decisions to produce private household energy by growing fuel trees on their farm. As participation in private fuel tree planting is not random, the study confronts a methodological issue in investigating the causal effect of forest protection intervention on rural farm households’ private fuel tree planting through non-parametric propensity score matching (PSM) method. The protection intervention on average has increased fuel tree planting by 503 (580.6%) compared to open access areas and indirectly contributed to slowing down the loss of biodiversity in the area. Land cover/use is a dynamic phenomenon that changes with time and space due to anthropogenic pressure and development. Forest cover and land use changes in Chiro District, Ethiopia over a period of 40 years was studied using remotely sensed data. Multi temporal satellite data of Landsat was used to map and monitor forest cover and land use changes occurred during three point of time of 1972,1986 and 2012. A pixel base supervised image classification was used to map land use land cover classes for maps of both time set. The result of change detection analysis revealed that the area has shown a remarkable land cover/land use changes in general and forest cover change in particular. Specifically, the dense forest cover land declined from 235 ha in 1972 to 51 ha in 1986. However, government interventions in forest protection in 1989 have slowed down the drastic change of dense forest cover loss around the protected area through reclaiming 1,300 hectares of deforested land through reforestation program up to 2012.

Permanent open water surfaces in central Siberia from ENVISAT ASAR wide swath data 2003/2004, link to shapefile

Permanent water bodies not only store dissolved CO2 but are essential for the maintenance of wetlands in their proximity. From the viewpoint of greenhouse gas (GHG) accounting wetland functions comprise sequestration of carbon under anaerobic conditions and methane release. The investigated area in central Siberia covers boreal and sub-arctic environments. Small inundated basins are abundant on the sub-arctic Taymir lowlands but also in parts of severe boreal climate where permafrost ice content is high and feature important freshwater ecosystems. Satellite radar imagery (ENVISAT ScanSAR), acquired in summer 2003 and 2004, has been used to derive open water surfaces with 150 m resolution, covering an area of approximately 3 Mkm**2. The open water surface maps were derived using a simple threshold-based classification method. The results were assessed with Russian forest inventory data, which includes detailed information about water bodies. The resulting classification has been further used to estimate the extent of tundra wetlands and to determine their importance for methane emissions. Tundra wetlands cover 7% (400,000 km**2) of the study region and methane emissions from hydromorphic soils are estimated to be 45,000 t/d for the Taymir peninsula.

Description of manganese deposits from the Todoroki mine, Hokkaido area, Japan

Manganese deposits are abundant in various places in the Oshima Peninsula southwest of Hokkaido. This is particular the case of Todoroki Mine situated about 25 kilometers to the south-east of the Ginzan railway station in Siribesi Province. It consists of manganese beds intermixed with a tertiary volcanic tuff complex overlaying granite.

Simple heterogeneity parametrization for sea surface temperature and chlorophyll

Using satellite maps this paper offers a complex analysis of chlorophyll & SST heterogeneity in the shelf seas around the southwest of the UK. The heterogeneity scaling follows a simple power law and is consequently pa- rametrized by two parameters. It is shown that in most cases these two parameters vary only relatively little with time. The paper offers a detailed comparison of field heterogeneity between different regions. How much hetero- geneity is in each region preserved in the annual median data is also determined. The paper explicitly demon- strates how one can use these results to calculate representative measurement area for in situ networks.

Simple heterogeneity parametrization for sea surface temperature and chlorophyll

Using satellite maps this paper offers a complex analysis of chlorophyll & SST heterogeneity in the shelf seas around the southwest of the UK. The heterogeneity scaling follows a simple power law and is consequently pa- rametrized by two parameters. It is shown that in most cases these two parameters vary only relatively little with time. The paper offers a detailed comparison of field heterogeneity between different regions. How much hetero- geneity is in each region preserved in the annual median data is also determined. The paper explicitly demon- strates how one can use these results to calculate representative measurement area for in situ networks.

Pulling Back the Veil: The Characterization and Habitability of Enshrouded Worlds

Thesis (Ph.D.)--University of Washington, 2016-08

Mapping of Submerged Aquatic Vegetation in Rivers From Very High Resolution Image Data, Using Object Based Image Analysis Combined with Expert Knowledge

The use of remote sensing for monitoring of submerged aquatic vegetation (SAV) in fluvial environments has been limited by the spatial and spectral resolution of available image data. The absorption of light in water also complicates the use of common image analysis methods. This paper presents the results of a study that uses very high resolution (VHR) image data, collected with a Near Infrared sensitive DSLR camera, to map the distribution of SAV species for three sites along the Desselse Nete, a lowland river in Flanders, Belgium. Plant species, including Ranunculus aquatilis L., Callitriche obtusangula Le Gall, Potamogeton natans L., Sparganium emersum L. and Potamogeton crispus L., were classified from the data using Object-Based Image Analysis (OBIA) and expert knowledge. A classification rule set based on a combination of both spectral and structural image variation (e.g. texture and shape) was developed for images from two sites. A comparison of the classifications with manually delineated ground truth maps resulted for both sites in 61% overall accuracy. Application of the rule set to a third validation image, resulted in 53% overall accuracy. These consistent results show promise for species level mapping in such biodiverse environments, but also prompt a discussion on assessment of classification accuracy.

A fast parallel hyperspectral coded aperture algorithm for compressive sensing using OpenCL

In this paper, we develop a fast implementation of an hyperspectral coded aperture (HYCA) algorithm on different platforms using OpenCL, an open standard for parallel programing on heterogeneous systems, which includes a wide variety of devices, from dense multicore systems from major manufactures such as Intel or ARM to new accelerators such as graphics processing units (GPUs), field programmable gate arrays (FPGAs), the Intel Xeon Phi and other custom devices. Our proposed implementation of HYCA significantly reduces its computational cost. Our experiments have been conducted using simulated data and reveal considerable acceleration factors. This kind of implementations with the same descriptive language on different architectures are very important in order to really calibrate the possibility of using heterogeneous platforms for efficient hyperspectral imaging processing in real remote sensing missions.

Assessing the influence of abiotic factors and leaf-level properties on the stability of growing-season canopy greenness in a deciduous forest

Maps depicting spatial pattern in the stability of summer greenness could advance understanding of how forest ecosystems will respond to global changes such as a longer growing season. Declining summer greenness, or “greendown”, is spectrally related to declining near-infrared reflectance and is observed in most remote sensing time series to begin shortly after peak greenness at the end of spring and extend until the beginning of leaf coloration in autumn,. Understanding spatial patterns in the strength of greendown has recently become possible with the advancement of Landsat phenology products, which show that greendown patterns vary at scales appropriate for linking these patterns to proposed environmental forcing factors. This study tested two non-mutually exclusive hypotheses for how leaf measurements and environmental factors correlate with greendown and decreasing NIR reflectance across sites. At the landscape scale, we used linear regression to test the effects of maximum greenness, elevation, slope, aspect, solar irradiance and canopy rugosity on greendown. Secondly, we used leaf chemical traits and reflectance observations to test the effect of nitrogen availability and intrinsic water use efficiency on leaf-level greendown, and landscape-level greendown measured from Landsat. The study was conducted using Quercus alba canopies across 21 sites of an eastern deciduous forest in North America between June and August 2014. Our linear model explained greendown variance with an R2=0.47 with maximum greenness as the greatest model effect. Subsequent models excluding one model effect revealed elevation and aspect were the two topographic factors that explained the greatest amount of greendown variance. Regression results also demonstrated important interactions between all three variables, with the greatest interaction showing that aspect had greater influence on greendown at sites with steeper slopes. Leaf-level reflectance was correlated with foliar δ13C (proxy for intrinsic water use efficiency), but foliar δ13C did not translate into correlations with landscape-level variation in greendown from Landsat. Therefore, we conclude that Landsat greendown is primarily indicative of landscape position, with a small effect of canopy structure, and no measureable effect of leaf reflectance. With this understanding of Landsat greendown we can better explain the effects of landscape factors on vegetation reflectance and perhaps on phenology, which would be very useful for studying phenology in the context of global climate change

Efeitos diretos e indiretos de parâmetros dos semivariogramas referentes aos atributos químicos do solo sobre a produtividade de grãos de soja

In order to add value to soybens crops, and hence the marketing, medium and large producers have been using precision agriculture techniques (PA), as the Remote Sensing, Geographic Information Systems (GIS) and positioning satellite, to assist the management of crops. Thus, given the economic relevance of that culture to the southwest of Paraná State and Brazil, scientific studies to increase their productivity and profitability are of main importance. The objective of this study was to evaluate the correlation between the chemical soil properties and soybean yield for each estimated parameter of semivariogram (range, nugget and level effect), and the deployment of these correlations in direct and indirect effects, aiming to improve the mapping process of spatial variability of soil chemical properties for use in PA. The hypothesis is that not all attributes of soil used to estimate the semivariogram parameters has a direct effect on productivity, and that even in groups of plants within a larger area it is possible to estimate the parameters of the semivariograms. The experiment was conducted in a commercial area of 19.7 ha, located in the city of Pato Branco - PR, central geographic coordinates 26º 11 '35 "South, 52 43' 05" West longitude, and average altitude of 780 m. The area is planted with soybeans for over 30 years, currently being adopted to cultivate Brasmax Target RR - Don Mario 5.9i, with row spacing of 0.50 m and 13 plants m-1, totaling 260,000 plants ha-1. For georeferencing of the area of study and sampling points was used a couple of topographic ProMarkTM3 receptors, making a relative positioning to obtain the georeferenced coordinates. To collect data (chemical analyzes of soil and crop yield) were sampled 10 blocks in the experimental area, each with an area of 20 m2 (20 meters long x 1 meter wide) containing two spaced adjacent rows of 0.5 m. Each block was divided into 20 portions of 1 m2, and from each were collected four subsamples at a distance of 0.5 m in relation to the lines of blocks, making up a sample depth for 0-10 cm a sample to 10-20 cm for each plot, totaling 200 samples for each depth. The soybean crop was performed on the blocks depending on maturity, and in each block was considered a bundle at each meter. In the data analysis, it was performed a diagnosis of multicollinearity, and subsequently a path analysis of the main variables according to the explanatory variables (range of chemical attributes: pH, K, P, Ca, etc.). The results obtained by the path analysis of the parameters of the semivariogram of soil chemical properties, indicated that only the Fe, Mg, Mn, organic matter (OM), P and Saturation by bases (SB) exerted direct and indirect effects on soybean productivity, although they have not presented spatial variability, indicating that the distribution of blocks in the area was unable to identify the spatial dependence of these elements, making it impossible to draw up maps of the chemical attributes for use in PA.

DopSCAT: A mission concept for simultaneous measurements of marine winds and surface currents

A radar scatterometer operates by transmitting a pulse of microwave energy toward the ocean's surface and measuring the normalized (per-unit-surface) radar backscatter coefficient (σ°). The primary application of scatterometry is the measurement of near-surface ocean winds. By combining σ° measurements from different azimuth angles, the 10 m vector wind can be determined through a Geophysical Model Function (GMF), which relates wind and backscatter. This paper proposes a mission concept for the measurement of both oceanic winds and surface currents, which makes full use of earlier C-band radar remote sensing experience. For the determination of ocean currents, in particular, the novel idea of using two chirps of opposite slope is introduced. The fundamental processing steps required to retrieve surface currents are given together with their associated accuracies. A detailed description of the mission proposal and comparisons between real and retrieved surface currents are presented. The proposed ocean Doppler scatterometer can be used to generate global surface ocean current maps with accuracies better than 0.2 m/s at a spatial resolution better than 25 km (i.e., 12.5 km spatial sampling) on a daily basis. These maps will allow gaining some insights on the upper ocean mesoscale dynamics. The work lies at a frontier, given that the present inability to measure ocean currents from space in a consistent and synoptic manner represents one of the greatest weaknesses in ocean remote sensing.