Precise Orbit Determination of Low Earth Satellites at AIUB using GPS and SLR data

An ever increasing number of low Earth orbiting (LEO) satellites is, or will be, equipped with retro-reflectors for Satellite Laser Ranging (SLR) and on-board receivers to collect observations from Global Navigation Satellite Systems (GNSS) such as the Global Positioning Sys- tem (GPS) and the Russian GLONASS and the European Galileo systems in the future. At the Astronomical Insti- tute of the University of Bern (AIUB) LEO precise or- bit determination (POD) using either GPS or SLR data is performed for a wide range of applications for satellites at different altitudes. For this purpose the classical numeri- cal integration techniques, as also used for dynamic orbit determination of satellites at high altitudes, are extended by pseudo-stochastic orbit modeling techniques to effi- ciently cope with potential force model deficiencies for satellites at low altitudes. Accuracies of better than 2 cm may be achieved by pseudo-stochastic orbit modeling for satellites at very low altitudes such as for the GPS-based POD of the Gravity field and steady-state Ocean Circula- tion Explorer (GOCE).

Earth gravity field recovery using GPS, GLONASS, and SLR satellites

The time variable Earth’s gravity field provides the information about mass transport within the system Earth, i.e., the relationship of mass transport between atmosphere, oceans, and land hydrology. We recover the low-degree parameters of the time variable gravity field using microwave observations from GPS and GLONASS satellites and from SLR data to five geodetic satellites, namely LAGEOS-1/2, Starlette, Stella, and AJISAI. GPS satellites are particularly sensitive to specific coefficients of the Earth's gravity field, because of the deep 2:1 orbital resonance with Earth rotation (two revolutions of the GPS satellites per sidereal day). The resonant coefficients cause, among other, a “secular” drift (actually periodic variations of very long periods) of the semi-major axes of up to 5.3 m/day of GPS satellites. We processed 10 years of GPS and GLONASS data using the standard orbit models from the Center of Orbit Determination in Europe (CODE) with a simultaneous estimation of the Earth gravity field coefficients and other parameters, e.g., satellite orbit parameters, station coordinates, Earth rotation parameters, troposphere delays, etc. The weekly GNSS gravity solutions up to degree and order 4/4 are compared to the weekly SLR gravity field solutions. The SLR-derived geopotential coefficients are compared to monthly GRACE and CHAMP results.

Processing 20 years of SLR observations to GNSS satellites

We process 20 years of SLR observations to GPS and GLONASS satellites using the reprocessed 3-day and 1-day microwave orbits provided by the Center for Orbit Determination in Europe (CODE) for the period 1994-2013. We study the dependency of the SLR residuals on the type, size, and a number of corner cubes in satellite laser reflector arrays (LRA). We show that the mean SLR residuals and the RMS of residuals depend on the coating of LRA and the block or type of GNSS satellites. The SLR mean residuals are also a function of the equipment used at SLR stations including detector types and detecting modes.

Time variable Earth’s gravity field from SLR satellites

The time variable Earth’s gravity field contains information about the mass transport within the system Earth, i.e., the relationship between mass variations in the atmosphere, oceans, land hydrology, and ice sheets. For many years, satellite laser ranging (SLR) observations to geodetic satellites have provided valuable information of the low-degree coefficients of the Earth’s gravity field. Today, the Gravity Recovery and Climate Experiment (GRACE) mission is the major source of information for the time variable field of a high spatial resolution. We recover the low-degree coefficients of the time variable Earth’s gravity field using SLR observations up to nine geodetic satellites: LAGEOS-1, LAGEOS-2, Starlette, Stella, AJISAI, LARES, Larets, BLITS, and Beacon-C. We estimate monthly gravity field coefficients up to degree and order 10/10 for the time span 2003–2013 and we compare the results with the GRACE-derived gravity field coefficients. We show that not only degree-2 gravity field coefficients can be well determined from SLR, but also other coefficients up to degree 10 using the combination of short 1-day arcs for low orbiting satellites and 10-day arcs for LAGEOS-1/2. In this way, LAGEOS-1/2 allow recovering zonal terms, which are associated with long-term satellite orbit perturbations, whereas the tesseral and sectorial terms benefit most from low orbiting satellites, whose orbit modeling deficiencies are minimized due to short 1-day arcs. The amplitudes of the annual signal in the low-degree gravity field coefficients derived from SLR agree with GRACE K-band results at a level of 77 %. This implies that SLR has a great potential to fill the gap between the current GRACE and the future GRACE Follow-On mission for recovering of the seasonal variations and secular trends of the longest wavelengths in gravity field, which are associated with the large-scale mass transport in the system Earth.

Search for satellites near comet 67P/Churyumov-Gerasimenko using Rosetta/OSIRIS images

Context. The European Space Agency Rosetta mission reached and started escorting its main target, the Jupiter-family comet 67P/Churyumov-Gerasimenko, at the beginning of August 2014. Within the context of solar system small bodies, satellite searches from approaching spacecraft were extensively used in the past to study the nature of the visited bodies and their collisional environment. Aims. During the approaching phase to the comet in July 2014, the OSIRIS instrument onboard Rosetta performed a campaign aimed at detecting objects in the vicinity of the comet nucleus and at measuring these objects' possible bound orbits. In addition to the scientific purpose, the search also focused on spacecraft security to avoid hazardous material in the comet's environment. Methods. Images in the red spectral domain were acquired with the OSIRIS Narrow Angle Camera, when the spacecraft was at a distance between 5785 km and 5463 km to the comet, following an observational strategy tailored to maximize the scientific outcome. From the acquired images, sources were extracted and displayed to search for plausible displacements of all sources from image to image. After stars were identified, the remaining sources were thoroughly analyzed. To place constraints on the expected displacements of a potential satellite, we performed Monte Carlo simulations on the apparent motion of potential satellites within the Hill sphere. Results. We found no unambiguous detections of objects larger than similar to 6 m within similar to 20 km and larger than similar to 1 m between similar to 20 km and similar to 110 km from the nucleus, using images with an exposure time of 0.14 s and 1.36 s, respectively. Our conclusions are consistent with independent works on dust grains in the comet coma and on boulders counting on the nucleus surface. Moreover, our analysis shows that the comet outburst detected at the end of April 2014 was not strong enough to eject large objects and to place them into a stable orbit around the nucleus. Our findings underline that it is highly unlikely that large objects survive for a long time around cometary nuclei.

Robust trends of landscape dynamics in the Arctic Lena Delta with temporally dense Landsat time-series stacks, with links to GeoTIFFs

Arctic permafrost landscapes are among the most vulnerable and dynamic landscapes globally, but due to their extent and remoteness most of the landscape changes remain unnoticed. In order to detect disturbances in these areas we developed an automated processing chain for the calculation and analysis of robust trends of key land surface indicators based on the full record of available Landsat TM, ETM +, and OLI data. The methodology was applied to the ~ 29,000 km**2 Lena Delta in Northeast Siberia, where robust trend parameters (slope, confidence intervals of the slope, and intercept) were calculated for Tasseled Cap Greenness, Wetness and Brightness, NDVI, and NDWI, and NDMI based on 204 Landsat scenes for the observation period between 1999 and 2014. The resulting datasets revealed regional greening trends within the Lena Delta with several localized hot-spots of change, particularly in the vicinity of the main river channels. With a 30-m spatial resolution various permafrost-thaw related processes and disturbances, such as thermokarst lake expansion and drainage, fluvial erosion, and coastal changes were detected within the Lena Delta region, many of which have not been noticed or described before. Such hotspots of permafrost change exhibit significantly different trend parameters compared to non-disturbed areas. The processed dataset, which is made freely available through the data archive PANGAEA, will be a useful resource for further process specific analysis by researchers and land managers. With the high level of automation and the use of the freely available Landsat archive data, the workflow is scalable and transferrable to other regions, which should enable the comparison of land surface changes in different permafrost affected regions and help to understand and quantify permafrost landscape dynamics.

Mapping decadal land cover changes in the woodlands of north eastern Namibia using the Landsat satellite archive (1975-2014)

Woodland savannahs provide essential ecosystem functions and services to communities. On the African continent, they are widely utilized and converted to intensive land uses. This study investigates the land cover changes of 108,038 km**2 in NE Namibia using multi-temporal, multi-sensor Landsat imagery, at decadal intervals from 1975 to 2014, with a post-classification change detection method and supervised Regression Tree classifiers. We discuss likely impacts of land tenure and reforms over the past four decades on changes in land use and land cover. These changes included losses, gains and exchanges between predominant land cover classes. Exchanges comprised logical conversions between woodland and agricultural classes, implying woodland clearing for arable farming, cropland abandonment and vegetation succession. The most dominant change was a reduction in the area of the woodland class due to the expansion of the agricultural class, specifically, small-scale cereal and pastoral production. Woodland area decreased from 90% of the study area in 1975 to 83% in 2014, while cleared land increased from 9% to 14%. We found that the main land cover changes are conversion from woodland to agricultural and urban land uses, driven by urban expansion and woodland clearing for subsistence-based agriculture and pastoralism.

Table 1: Glacier measurements from air photos 1969 and landsat images 1972 and 1973

The extent of snow cover at the end of the ablation season on glaciers in the Tyrolean Alps in 1972 and 1973 was determined from Landsat-1 Multispectral Scanner (MSS) images. For snovv mapping the MSS-images with a ground resolution of 80 meters were enlarged to a scale of 1: 100.000 by photographic methods. Different appearance of snow cover in the 4 MSS-channels is discussed in connection with ground truth control. The accuracy of snow and ice mapping from Landsat images was checked on 15 glaciers with an area from 1 to 10 km2 by aerial photography and/or ground truth control. These comparisons imply the usefulness of Landsat images for snow mapping on glaciers of a few square kilometers. The altitude of the equilibrium line was determined from Landsat images for 53 glaciers in the Tyrolean Alps. The regional differences in the equilibrium line altitude correspond to the regional precipitation patterns. The equilibrium line was identical with the snow line at the end of the budget year 1971/1972; therefore it was possible to determine the equilibrium line from satellite images. For 1968/69 the equilibrium line was mapped from aerial photographs for several glaciers. In 1972/73 mass balance was strongly negative and the equilibrimn line was within the firn area of the glaciers. Therefore it was not possible to distinguish between accumulation and ablation areas from the Landsat images of September 1973; however, snow and ice areas could be olearly differentiated. The ratios of accumulation area 01' snow area to the total area of the glaciers were determineel from satellite images and aerial photography separately for aelvancing anel for retreating glaciers and were relateel to the mass balance. In the budget years 1968/69 and 1972/73 with negative mass balance the accumulation area ratios of the advancing glacien; were olearly different from the ratios of the retreating glaciers, in 1971/72 with positive 01' balanced mass budget the differences between advancing and retreating glaciers were not significant.

Landsat derived percent vegetation cover of Israel (as of 2003)

Wildfires are part of the Mediterranean ecosystem, however, in Israel all wildfires are human caused, either intentionally or un-intentionally. In this study we aimed to develop and test a new method for mapping fire scars from MODIS imagery, to examine the temporal and spatial patterns of wildfires in Israel in the 2000s and to examine the factors controlling Israel's wildfire regime. To map the fires we used two 'off-the-shelf' MODIS fire products as our basis-the 1 km MODIS Collection 5 fire hotspots, the 500 m MCD45A1 burnt areas-and we created a new set of fire scar maps from the 250 m MOD13Q1 product. We carried out a cross comparison of the three MODIS based wildfire scar maps and evaluated them independently against the wild fire scars mapped from 30 m Landsat TM imagery. To examine the factors controlling wildfires we used GIS layers of rainfall, land use, and a Landsat-based national vegetation map. Wildfires occurred in areas where annual rainfall was above 250 mm, mostly in areas with herbaceous vegetation. Wildfire frequency was especially high in the Golan Heights and in the foothills of the Judean mountains, and a high correspondence was found between military training zones and the spatial distribution of fire scars. The use of MODIS satellite images enabled us to map wildfires at a national scale due to the high temporal resolution of the sensor. Our MOD13Q1 based mapping of fire scars adequately mapped large (>1 km**2) fires with accuracies above 80%. Such large fires account for a large proportion of all fires, and pose the greatest threats. This database can aid managers in determining wildfire risks in space and in time.

Classified Adélie penguin colonies from Landsat data

Breeding distribution of the Adelie penguin, Pygoscelis adeliae, was surveyed with Landsat-7 Enhanced Thematic Mapper Plus (ETM+) data along the coastline of Antarctica, an area covering approximately 330° of longitude. An algorithm was designed to minimize the radiometric contribution from exogenous sources and to retrieve Adelie penguin colony location and spatial extent from the ETM+ data. In all, 9143 individual pixels were classified as belonging to an Adelie penguin colony class out of the entire dataset of 195 ETM+ scenes, where the dimension of each pixel is 30 m by 30 m, and each scene is approximately 180 km by 180 km. Pixel clustering identified a total of 187 individual Adelie penguin colonies, ranging in size from a single pixel (900 m**2) to a maximum of 875 pixels (0.788 km**2). Colony retrievals have a very low error of commission, on the order of 1 percent or less, and the error of omission was estimated to be 2.9 percent by population based on comparisons with direct observations from surveys across east Antarctica. Thus, the Landsat retrievals can successfully locate Adelie penguin colonies that account for ~97 percent of a regional population. Geographic coordinates and the spatial extent of each colony retrieved from the Landsat data are available publically. Regional analysis found several areas where the Landsat retrievals suggest populations that are significantly larger than published estimates. Six Adelie penguin colonies were found that are believed to be unreported in the literature.