Confronto tra stime GPS e SAR di deformazioni del suolo e di ritardi troposferici

Il Global Positioning System (GPS) e l’Interferometric Synthetic Aperture Radar (InSAR) sono due tecniche osservative di grande importanza che utilizzano segnali nel campo delle microonde. Questa tesi intende contribuire a sviluppare una base di confronto tra i risultati derivati da queste due tecniche osservative. Una parte del lavoro riguarda uno studio delle deformazioni del suolo, in particolare, la stima dei movimenti verticali e di quelli che riguardano la componente Est della posizione delle stazioni. Un secondo ambito di ricerca è invece focalizzato alla determinazione del ritardo introdotto, nella propagazione dei segnali GPS e SAR, dal loro passaggio in atmosfera. In particolare, si è studiato l’effetto della componente umida della troposfera.

Multi-node TOA-DOA cooperative LOS-NLOS localization : enabling high accuracy and reliability

This dissertation investigates high performance cooperative localization in wireless environments based on multi-node time-of-arrival (TOA) and direction-of-arrival (DOA) estimations in line-of-sight (LOS) and non-LOS (NLOS) scenarios. Here, two categories of nodes are assumed: base nodes (BNs) and target nodes (TNs). BNs are equipped with antenna arrays and capable of estimating TOA (range) and DOA (angle). TNs are equipped with Omni-directional antennas and communicate with BNs to allow BNs to localize TNs; thus, the proposed localization is maintained by BNs and TNs cooperation. First, a LOS localization method is proposed, which is based on semi-distributed multi-node TOA-DOA fusion. The proposed technique is applicable to mobile ad-hoc networks (MANETs). We assume LOS is available between BNs and TNs. One BN is selected as the reference BN, and other nodes are localized in the coordinates of the reference BN. Each BN can localize TNs located in its coverage area independently. In addition, a TN might be localized by multiple BNs. High performance localization is attainable via multi-node TOA-DOA fusion. The complexity of the semi-distributed multi-node TOA-DOA fusion is low because the total computational load is distributed across all BNs. To evaluate the localization accuracy of the proposed method, we compare the proposed method with global positioning system (GPS) aided TOA (DOA) fusion, which are applicable to MANETs. The comparison criterion is the localization circular error probability (CEP). The results confirm that the proposed method is suitable for moderate scale MANETs, while GPS-aided TOA fusion is suitable for large scale MANETs. Usually, TOA and DOA of TNs are periodically estimated by BNs. Thus, Kalman filter (KF) is integrated with multi-node TOA-DOA fusion to further improve its performance. The integration of KF and multi-node TOA-DOA fusion is compared with extended-KF (EKF) when it is applied to multiple TOA-DOA estimations made by multiple BNs. The comparison depicts that it is stable (no divergence takes place) and its accuracy is slightly lower than that of the EKF, if the EKF converges. However, the EKF may diverge while the integration of KF and multi-node TOA-DOA fusion does not; thus, the reliability of the proposed method is higher. In addition, the computational complexity of the integration of KF and multi-node TOA-DOA fusion is much lower than that of EKF. In wireless environments, LOS might be obstructed. This degrades the localization reliability. Antenna arrays installed at each BN is incorporated to allow each BN to identify NLOS scenarios independently. Here, a single BN measures the phase difference across two antenna elements using a synchronized bi-receiver system, and maps it into wireless channel’s K-factor. The larger K is, the more likely the channel would be a LOS one. Next, the K-factor is incorporated to identify NLOS scenarios. The performance of this system is characterized in terms of probability of LOS and NLOS identification. The latency of the method is small. Finally, a multi-node NLOS identification and localization method is proposed to improve localization reliability. In this case, multiple BNs engage in the process of NLOS identification, shared reflectors determination and localization, and NLOS TN localization. In NLOS scenarios, when there are three or more shared reflectors, those reflectors are localized via DOA fusion, and then a TN is localized via TOA fusion based on the localization of shared reflectors.

Detecting offsets in GPS time series: first results from the Detection of Offsets in GPS Experiment

The accuracy of Global Positioning System (GPS) time series is degraded by the presence of offsets. To assess the effectiveness of methods that detect and remove these offsets, we designed and managed the Detection of Offsets in GPS Experiment. We simulated time series that mimicked realistic GPS data consisting of a velocity component, offsets, white and flicker noises (1/f spectrum noises) composed in an additive model. The data set was made available to the GPS analysis community without revealing the offsets, and several groups conducted blind tests with a range of detection approaches. The results show that, at present, manual methods (where offsets are hand picked) almost always give better results than automated or semi‒automated methods (two automated methods give quite similar velocity bias as the best manual solutions). For instance, the fifth percentile range (5% to 95%) in velocity bias for automated approaches is equal to 4.2 mm/year (most commonly ±0.4 mm/yr from the truth), whereas it is equal to 1.8 mm/yr for the manual solutions (most commonly 0.2 mm/yr from the truth). The magnitude of offsets detectable by manual solutions is smaller than for automated solutions, with the smallest detectable offset for the best manual and automatic solutions equal to 5 mm and 8 mm, respectively. Assuming the simulated time series noise levels are representative of real GPS time series, robust geophysical interpretation of individual site velocities lower than 0.2–0.4 mm/yr is therefore certainly not robust, although a limit of nearer 1 mm/yr would be a more conservative choice. Further work to improve offset detection in GPS coordinates time series is required before we can routinely interpret sub‒mm/yr velocities for single GPS stations.

Geocenter coordinates estimated from GNSS data as viewed by perturbation theory

Time series of geocenter coordinates were determined with data of two global navigation satellite systems (GNSSs), namely the U.S. GPS (Global Positioning System) and the Russian GLONASS (Global’naya Nawigatsionnaya Sputnikowaya Sistema). The data was recorded in the years 2008–2011 by a global network of 92 permanently observing GPS/GLONASS receivers. Two types of daily solutions were generated independently for each GNSS, one including the estimation of geocenter coordinates and one without these parameters. A fair agreement for GPS and GLONASS was found in the geocenter x- and y-coordinate series. Our tests, however, clearly reveal artifacts in the z-component determined with the GLONASS data. Large periodic excursions in the GLONASS geocenter z-coordinates of about 40 cm peak-to-peak are related to the maximum elevation angles of the Sun above/below the orbital planes of the satellite system and thus have a period of about 4 months (third of a year). A detailed analysis revealed that the artifacts are almost uniquely governed by the differences of the estimates of direct solar radiation pressure (SRP) in the two solution series (with and without geocenter estimation). A simple formula is derived, describing the relation between the geocenter z-coordinate and the corresponding parameter of the SRP. The effect can be explained by first-order perturbation theory of celestial mechanics. The theory also predicts a heavy impact on the GNSS-derived geocenter if once-per-revolution SRP parameters are estimated in the direction of the satellite’s solar panel axis. Specific experiments using GPS observations revealed that this is indeed the case. Although the main focus of this article is on GNSS, the theory developed is applicable to all satellite observing techniques. We applied the theory to satellite laser ranging (SLR) solutions using LAGEOS. It turns out that the correlation between geocenter and SRP parameters is not a critical issue for the SLR solutions. The reasons are threefold: The direct SRP is about a factor of 30–40 smaller for typical geodetic SLR satellites than for GNSS satellites, allowing it in most cases to not solve for SRP parameters (ruling out the correlation between these parameters and the geocenter coordinates); the orbital arc length of 7 days (which is typically used in SLR analysis) contains more than 50 revolutions of the LAGEOS satellites as compared to about two revolutions of GNSS satellites for the daily arcs used in GNSS analysis; the orbit geometry is not as critical for LAGEOS as for GNSS satellites, because the elevation angle of the Sun w.r.t. the orbital plane is usually significantly changing over 7 days.

Spatial Patterns in Mass Balance of the Siple Coast and Amundsen Sea Sectors, West Antarctica

Local rates of change in ice-sheet thickness were calculated at IS sites in West Antarctica using the submergence velocity technique. This method entails a comparison of the vertical velocity of the ice sheet, measured using repeat global positioning system surveys of markers, and local long-term rates of snow accumulation obtained using firn-core stratigraphy. Any significant difference between these two quantities represents a thickness change with time. Measurements were conducted at sites located similar to 100-200 km apart along US ITASE traverse routes, and at several isolated locations. All but one of the sites are distributed in the Siple Coast and the Amundsen Sea basin along contours of constant elevation, along flowlines, across ice divides and close to regions of enhanced flow. Calculated rates of thickness change are different from site to site. Most of the large rates of change in ice thickness (similar to 10 cm a(-1) or larger) are observed in or close to regions of rapid flow, and are probably related to ice-dynamics effects. Near-steady-state conditions are calculated mostly at sites in the slow-moving ice-sheet interior and near the main West Antarctic ice divide. These results are consistent with regional estimates of ice-sheet change derived from remote-sensing measurements at similar locations in West Antarctica.

Mass Balance and Accumulation Rate Across Siple Dome, West Antarctica

Snow-accumulation rates and rates of ice-thickness change (mass balance) are studied at several sites on Siple Dome, West Antarctica. Accumulation rates are derived from analyses of gross beta radioactivity in shallow firn cores located along a 60 km transect spanning both flanks and the crest of the dome. There is a north-south gradient in snow-accumulation rate across the dome that is consistent with earlier radar mapping of internal stratigraphy. Orographic processes probably control this distribution. Mass balance is inferred from the difference between global positioning system (GPS)-derived vertical velocities and snow-accumulation rates for sites close to the firn-core locations. Results indicate that there is virtually no net thickness change at four of the five sites. The exception is at the northernmost site where a small amount of thinning is detected, that appears to be inconsistent with other studies. A possible cause of this anomalous thinning is recent retreat of the grounding line of Ice Stream D.

Local Rates of Ice-Sheet Thickness Change in Greenland

The rate of ice-sheet thickness change is calculated for 10 sites in Greenland by comparing measured values of ice vertical velocity and snow-accumulation rate. Vertical velocities are derived from repeat surveys of markers using precision global positioning system techniques, and accumulation rates are determined from stratigraphic analysis of firn cores. The results apply to time-scales covered by the firn-core records, which in most cases are a few decades. A spectrum of thickness-change rates is obtained, ranging from substantial thinning to slow thickening. The sites where ice-sheet thinning is indicated are located near the ice-sheet margin or in outlet glacier catchments. Interior and high-elevation sites are predominantly in balance or thickening slowly. Uncertainties in the rates of thickness change are dominated by errors in the determination of accumulation rates. The results of this work are broadly comparable with regional estimates of mass balance obtained from the analysis of catchment input vs discharge.

Variability in Accumulation Rates from GPR Profiling on the West Antarctic Plateau

lsochronal layers in firn detected with ground-penetrating radar (GPR) and dated using results from ice-core analyses are used to calculate accumulation rates along a 100 km across-flow profile in West Antarctica. Accumulation rates are shown to be highly variable over short distances. Elevation measurements from global positioning system surveys show that accumulation rates derived from shallow horizons correlate well with surface undulations, which implies that wind redistribution of snow is the leading cause of this variability. Temporal changes in accumulation rate over 25-185 year intervals are smoothed to along-track length scales comparable to surface undulations in order to identify trends in accumulation that are likely related to changes in climate. Results show that accumulation rates along this profile have decreased in recent decades, which is consistent with core-derived time series of annual accumulation rates measured at the two ends of the radar profile. These results suggest that temporal variability observed in accumulation-rate records from ice cores and GPR profiles can be obscured by spatial influences, although it is possible to resolve temporal signals if the effects of local topography and ice flow are quantified and removed.

Spatial and Temporal Melt Variability at Helheim Glacier, East Greenland, and Its Effect on Ice Dynamics

Understanding the behavior of large outlet glaciers draining the Greenland Ice Sheet is critical for assessing the impact of climate change on sea level rise. The flow of marine-terminating outlet glaciers is partly governed by calving-related processes taking place at the terminus but is also influenced by the drainage of surface runoff to the bed through moulins, cracks, and other pathways. To investigate the extent of the latter effect, we develop a distributed surface-energy-balance model for Helheim Glacier, East Greenland, to calculate surface melt and thereby estimate runoff. The model is driven by data from an automatic weather station operated on the glacier during the summers of 2007 and 2008, and calibrated with independent measurements of ablation. Modeled melt varies over the deployment period by as much as 68% relative to the mean, with melt rates approximately 77% higher on the lower reaches of the glacier trunk than on the upper glacier. We compare melt variations during the summer season to estimates of surface velocity derived from global positioning system surveys. Near the front of the glacier, there is a significant correlation (on >95% levels) between variations in runoff (estimated from surface melt) and variations in velocity, with a 1 day delay in velocity relative to melt. Although the velocity changes are small compared to accelerations previously observed following some calving events, our findings suggest that the flow speed of Helheim Glacier is sensitive to changes in runoff. The response is most significant in the heavily crevassed, fast-moving region near the calving front. The delay in the peak of the cross-correlation function implies a transit time of 12-36 h for surface runoff to reach the bed.

Surface Deformations as Indicators of Deep Ebullition Fluxes in a Large Northern Peatland

Peatlands deform elastically during precipitation cycles by small (+/- 3 cm) oscillations in surface elevation. In contrast, we used a Global Positioning System network to measure larger oscillations that exceeded 20 cm over periods of 4 - 12 hours during two seasonal droughts at a bog and fen site in northern Minnesota. The second summer drought also triggered 19 depressuring cycles in an overpressured stratum under the bog site. The synchronicity between the largest surface deformations and the depressuring cycles indicates that both phenomena are produced by the episodic release of large volumes of gas from deep semi-elastic compartments confined by dense wood layers. We calculate that the three largest surface deformations were associated with the release of 136 g CH4 m(-2), which exceeds by an order of magnitude the annual average chamber fluxes measured at this site. Ebullition of gas from the deep peat may therefore be a large and previously unrecognized source of radiocarbon depleted methane emissions from northern peatlands.

Passively Track WiFi Users with an Enhanced Particle Filter using Power-based Ranging

Passive positioning systems produce user location information for third-party providers of positioning services. Since the tracked wireless devices do not participate in the positioning process, passive positioning can only rely on simple, measurable radio signal parameters, such as timing or power information. In this work, we provide a passive tracking system for WiFi signals with an enhanced particle filter using fine-grained power-based ranging. Our proposed particle filter provides an improved likelihood function on observation parameters and is equipped with a modified coordinated turn model to address the challenges in a passive positioning system. The anchor nodes for WiFi signal sniffing and target positioning use software defined radio techniques to extract channel state information to mitigate multipath effects. By combining the enhanced particle filter and a set of enhanced ranging methods, our system can track mobile targets with an accuracy of 1.5m for 50% and 2.3m for 90% in a complex indoor environment. Our proposed particle filter significantly outperforms the typical bootstrap particle filter, extended Kalman filter and trilateration algorithms.

Introduction to the Maïdo Lidar Calibration Campaign dedicated to the validation of upper air meteorological parameters

The first operations at the new High-altitude Maïdo Observatory at La Réunion began in 2013. The Maïdo Lidar Calibration Campaign (MALICCA) was organized there in April 2013 and has focused on the validation of the thermodynamic parameters (temperature, water vapor, and wind) measured with many instruments including the new very large lidar for water vapor and temperature profiles. The aim of this publication consists of providing an overview of the different instruments deployed during this campaign and their status, some of the targeted scientific questions and associated instrumental issues. Some specific detailed studies for some individual techniques were addressed elsewhere. This study shows that temperature profiles were obtained from the ground to the mesopause (80 km) thanks to the lidar and regular meteorological balloon-borne sondes with an overlap range showing good agreement. Water vapor is also monitored from the ground to the mesopause by using the Raman lidar and microwave techniques. Both techniques need to be pushed to their limit to reduce the missing range in the lower stratosphere. Total columns obtained from global positioning system or spectrometers are valuable for checking the calibration and ensuring vertical continuity. The lidar can also provide the vertical cloud structure that is a valuable complementary piece of information when investigating the water vapor cycle. Finally, wind vertical profiles, which were obtained from sondes, are now also retrieved at Maïdo from the newly implemented microwave technique and the lidar. Stable calibrations as well as a small-scale dynamical structure are required to monitor the thermodynamic state of the middle atmosphere, ensure validation of satellite sensors, study the transport of water vapor in the vicinity of the tropical tropopause and study their link with cirrus clouds and cyclones and the impact of small-scale dynamics (gravity waves) and their link with the mean state of the mesosphere.

Community Trails in Community Context: GIS Analysis of Associations between Trail and Neighborhood Characteristics

This research examines the site and situation characteristics of community trails as landscapes promoting physical activity. Trail segment and neighborhood characteristics for six trails in urban, suburban, and exurban towns in northeastern Massachusetts were assessed from primary Global Positioning System (GPS) data and from secondary Census and land use data integrated in a geographic information system (GIS). Correlations between neighborhood street and housing density, land use mix, and sociodemographic characteristics and trail segment characteristics and amenities measure the degree to which trail segment attributes are associated with the surrounding neighborhood characteristics.

An evaluation of a Lake Houston tributary, Cypress Creek, for contamination and water quality

Introduction. Lake Houston serves as a reservoir for both recreational and drinking water for residents of Houston, Texas, and the metropolitan area. The Texas Commission on Environmental Quality (TCEQ) expressed concerns about the water quality and increasing amounts of pathogenic bacteria in Lake Houston (3). The objective of this investigation is to evaluate water quality for the presence of bacteria, nitrates, nitrites, carbon, phosphorus, dissolved oxygen, pH, turbidity, suspended solids, dissolved solids, and chlorine in Cypress Creek. The aims of this project are to analyze samples of water from Cypress Creek and to render a quantitative and graphical representation of the results. The collected information will allow for a better understanding of the aqueous environment in Cypress Creek.^ Methods. Water samples were collected in August 2009 and analyzed in the field and at UTSPH laboratory by spectrophotometry and other methods. Mapping software was utilized to develop novel maps of the sample sites using coordinates attained with the Global Positioning System (GPS). Sample sites and concentrations were mapped using Geographic Information System (GIS) software and correlated with permitted outfalls and other land use characteristic.^ Results. All areas sampled were positive for the presence of total coliform and Escherichia coli (E. coli). The presences of other water contaminants varied at each location in Cypress Creek but were under the maximum allowable limits designated by the Texas Commission on Environmental Quality. However, dissolved oxygen concentrations were elevated above the TCEQ limit of 5.0 mg/L at majority of the sites. One site had near-limit concentration of nitrates at 9.8 mg/L. Land use above this site included farm land, agricultural land, golf course, parks, residential neighborhoods, and nine permitted TCEQ effluent discharge sites within 0.5 miles upstream.^ Significance. Lake Houston and its tributary, Cypress Creek, are used as recreational waters where individuals may become exposed to microbial contamination. Lake Houston also is the source of drinking water for much of Houston/Harris and Galveston Counties. This research identified the presence of microbial contaminates in Cypress Creek above TCEQ regulatory requirements. Other water quality variables measured were in line with TCEQ regulations except for near-limit for nitrate at sample site #10, at Jarvis and Timberlake in Cypress Texas.^

Investigating the epidemic of heroin use in Dar es Salaam, Tanzania: Recognizing the issues of heroin smokers

Since heroin was introduced to East Africa during the 1980s, heroin use practices have changed rapidly in response to various internal and external pressures. The aim of this study was to identify and describe the population of heroin users and locations of heroin use in Dar es Salaam, Tanzania, in order to understand recent contexts of heroin use. The study took place between June 30 and August 19, 2011, in all three districts (Kinondoni, Ilala, and Temeke) of Dar es Salaam. We mapped sites using a Global Positioning System device, counted numbers of heroin users, and conducted informal interviews with heroin users. The mixed-methods analyses of the data included quantifying the basic demographic and aggregate information about the sites and heroin users, as well as qualitative analysis and coding of fieldnotes from observations and responses to interviews which was used to identify themes and characteristics of heroin users. ^ We identified a total of 150 sites and counted a total of 1046 male and 46 female non-injecting drug users and 78 male and 9 female injecting drug users (IDUs) of heroin. We found that social organization existed at some of the sites, with 31% (n=47) of sites reporting having a leader and 44% (n=66) of sites reporting mutual aid between users frequenting the site. We had difficulty locating IDUs and female drug users, and the majority of users we encountered were heroin smokers of kokteli, a mixture of heroin, cannabis, and/or tobacco which is smoked like a cigarette. ^ This research highlighted heroin smokers’ desire for access to drug treatment services. The current methadone-based medication assisted treatment (MAT) program is funded and operates as an HIV prevention program for IDUs to reduce HIV infection in this population and slow or stop the spread of a second wave of HIV infection in the general population. However, smokers perceived MAT to be primarily a drug use prevention or cessation program and felt unjustly neglected from the intervention, leading to a tense relationship with IDUs. From a public health standpoint, future interventions should include heroin smokers to prevent HIV transmission. ^