Lava discharge rate estimates from thermal infrared satellite data at Pacaya Volcano, Guatemala

Time-averaged discharge rates (TADR) were calculated for five lava flows at Pacaya Volcano (Guatemala), using an adapted version of a previously developed satellite-based model. Imagery acquired during periods of effusive activity between the years 2000 and 2010 were obtained from two sensors of differing temporal and spatial resolutions; the Moderate Resolution Imaging Spectroradiometer (MODIS), and the Geostationary Operational Environmental Satellites (GOES) Imager. A total of 2873 MODIS and 2642 GOES images were searched manually for volcanic “hot spots”. It was found that MODIS imagery, with superior spatial resolution, produced better results than GOES imagery, so only MODIS data were used for quantitative analyses. Spectral radiances were transformed into TADR via two methods; first, by best-fitting some of the parameters (i.e. density, vesicularity, crystal content, temperature change) of the TADR estimation model to match flow volumes previously estimated from ground surveys and aerial photographs, and second by measuring those parameters from lava samples to make independent estimates. A relatively stable relationship was defined using the second method, which suggests the possibility of estimating lava discharge rates in near-real-time during future volcanic crises at Pacaya.

Investigations into the degassing and eruption mechanisms of Nyamuragira volcano, Democratic Republic of the Congo (Africa)

One of two active volcanoes in the western branch of the East African Rift, Nyamuragira (1.408ºS, 29.20ºE; 3058 m) is located in the D.R. Congo. Nyamuragira emits large amounts of SO2 (up to ~1 Mt/day) and erupts low-silica, alkalic lavas, which achieve flow rates of up to ~20 km/hr. The source of the large SO2 emissions and pre-eruptive magma conditions were unknown prior to this study, and 1994-2010 lava volumes were only recently mapped via satellite imagery, mainly due to the region’s political instability. In this study, new olivine-hosted melt inclusion volatile (H2O, CO2, S, Cl, F) and major element data from five historic Nyamuragira eruptions (1912, 1938, 1948, 1986, 2006) are presented. Melt compositions derived from the 1986 and 2006 tephra samples best represent pre-eruptive volatile compositions because these samples contain naturally glassy inclusions that underwent less post-entrapment modification than crystallized inclusions. The total amount of SO2 released from the 1986 (0.04 Mt) and 2006 (0.06 Mt) eruptions are derived using the petrologic method, whereby S contents in melt inclusions are scaled to erupted lava volumes. These amounts are significantly less than satellite-based SO2 emissions for the same eruptions (1986 = ~1 Mt; 2006 = ~2 Mt). Potential explanations for this observation are: 1) accumulation of a vapor phase within the magmatic system that is only released during eruptions, and/or 2) syn-eruptive gas release from unerupted magma. Post-1994 Nyamuragira lava volumes were not available at the beginning of this study. These flows (along with others since 1967) are mapped with Landsat MSS, TM, and ETM+, Hyperion, and ALI satellite data and combined with published flow thicknesses to derive volumes. Satellite remote sensing data was also used to evaluate Nyamuragira SO2 emissions. These results show that the most recent Nyamuragira eruptions injected SO2 into the atmosphere between 15 km (2006 eruption) and 5 km (2010 eruption). This suggests that past effusive basaltic eruptions (e.g., Laki 1783) are capable of similar plume heights that reached the upper troposphere or tropopause, allowing SO2 and resultant aerosols to remain longer in the atmosphere, travel farther around the globe, and affect global climates.

Evolution of the Southern Kenya Rift from Miocene to present with a focus on the Magadi area

The Kenya (a.k.a., Gregory) Rift is a geologically active area located within the eastern branch of the larger East African Rift System (EARS). The study area is located in the southern Kenya Rift between 1° South and the Kenya-Tanzania border (covering approximately 1.5 square degrees, semi-centered on Lake Magadi) and is predominantly filled with extrusive igneous rocks (mostly basalts, phonolites and trachytes) of Miocene age or younger. Sediments are thin, less than 1.5Ma, and are confined to small grabens. The EARS can serve both as an analogue for ancient continental rifting and as a modern laboratory to observe the geologic processes responsible for rifting. This study demonstrates that vintage (as in older, quality maps published by the Kenya Geological Survey, that may be outdated based on newer findings) quarter-degree maps can be successfully combined with recently published data, and used to interpret satellite (mainly Landsat 7) images to produce versatile, updated digital maps. The study area has been remapped using this procedure and although it covers a large area, the mapping retains a quadrangle level of detail. Additionally, all geologic mapping elements (formations, faults, etc.) have been correlated across older map boundaries so that geologic units don't end artificially at degree boundaries within the study area. These elements have also been saved as individual digital files to facilitate future analysis. A series of maps showing the evolution of the southern Kenya rift from the Miocene to the present was created by combining the updated geologic map with age dates for geologic formations and fault displacements. Over 200 age dates covering the entire length of the Kenya Rift have been compiled for this study, and 6 paleo-maps were constructed to demonstrate the evolution of the area, starting with the eruption of the Kishalduga and Lisudwa melanephelinites onto the metamorphic basement around 15Ma. These eruptions occurred before the initial rift faulting and were followed by a massive eruption of phonolites between 13-10 Ma that covered most of the Kenya dome. This was followed by a period of relative quiescence, until the initial faulting defined the western boundary of the rift around 7Ma. The resulting graben was asymmetrical until corresponding faults to the east developed around 3Ma. The rift valley was flooded by basalts and trachytes between 3Ma and 700ka, after which the volcanic activity slowed to a near halt. Since 700ka most of the deposition has been comprised of sediments, mainly from lakes occupying the various basins in the area. The main results of this study are, in addition to a detailed interpretation of the rift development, a new geologic map that correlates dozens of formations across old map boundaries and a compilation of over 300 age dates. Specific products include paleomaps, tables of fault timing and displacement, and volume estimates of volcanic formations. The study concludes with a generalization of the present environment at Magadi including discussions of lagoon chemistry, mantle gases in relation to the trona deposit, and biology of the hot springs. Several biologic samples were collected during the 2006 field season in an attempt to characterize the organisms that are commonly seen in the present Lake Magadi environment. Samples were selected to represent the different, distinctive forms that are found in the hotsprings. Each sample had it own distinctive growth habit, and analysis showed that each was formed by a different cyanobacterial. Actual algae was rare in the collected samples, and represented by a few scattered diatoms.

A comparison of pixel and object-based land cover classification. A case study of the Asmara region, Eritrea

In this paper we compare the performance of two image classification paradigms (object- and pixel-based) for creating a land cover map of Asmara, the capital of Eritrea and its surrounding areas using a Landsat ETM+ imagery acquired in January 2000. The image classification methods used were maximum likelihood for the pixel-based approach and Bhattacharyya distance for the object-oriented approach available in, respectively, ArcGIS and SPRING software packages. Advantages and limitations of both approaches are presented and discussed. Classifications outputs were assessed using overall accuracy and Kappa indices. Pixel- and object-based classification methods result in an overall accuracy of 78% and 85%, respectively. The Kappa coefficient for pixel- and object-based approaches was 0.74 and 0.82, respectively. Although pixel-based approach is the most commonly used method, assessment and visual interpretation of the results clearly reveal that the object-oriented approach has advantages for this specific case-study.

HELIOSPHERIC ENERGETIC NEUTRAL HYDROGEN MEASURED WITH ASPERA-3 AND ASPERA-4

We re-analyze the signal of non-planetary energetic neutral atoms (ENAs) in the 0.4-5.0 keV range measured with the Neutral Particle Detector (NPD) of the ASPERA-3 and ASPERA-4 experiments on board the Mars and Venus Express satellites. Due to improved knowledge of sensor characteristics and exclusion of data sets affected by instrument effects, the typical intensity of the ENA signal obtained by ASPERA-3 is an order of magnitude lower than in earlier reports. The ENA intensities measured with ASPERA-3 and ASPERA-4 now agree with each other. In the present analysis, we also correct the ENA signal for Compton-Getting and for ionization loss processes under the assumption of a heliospheric origin. We find spectral shapes and intensities consistent with those measured by the Interstellar Boundary Explorer (IBEX). The principal advantage of ASPERA with respect to the IBEX sensors is the two times better spectral resolution. In this study, we discuss the physical significance of the spectral shapes and their potential variation across the sky. At present, these observations are the only independent test of the heliospheric ENA signal measured with IBEX in this energy range. The ASPERA measurements also allow us to check for a temporal variation of the heliospheric signal as they were obtained between 2003 and 2007, whereas IBEX has been operational since the end of 2008.

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.

SLR-derived terrestrial reference frame using the observations to LAGEOS-1/2, Starlette, Stella, and AJISAI

Currently, the contributions of Starlette, Stella, and AJISAI are not taken into account when defining the International Terrestrial Reference Frame (ITRF), despite the large amount of data collected in a long time-span. Consequently, the SLR-derived parameters and the SLR part of the ITRF are almost exclusively defined by LAGEOS-1 and LAGEOS-2. We investigate the potential of combining the observations to several SLR satellites with different orbital characteristics. Ten years of SLR data are homogeneously processed using the development version 5.3 of the Bernese GNSS Software. Special emphasis is put on orbit parameterization and the impact of LEO data on the estimation of the geocenter coordinates, Earth rotation parameters, Earth gravity field coefficients, and the station coordinates in one common adjustment procedure. We find that the parameters derived from the multi-satellite solutions are of better quality than those obtained in single satellite solutions or solutions based on the two LAGEOS satellites. A spectral analysis of the SLR network scale w.r.t. SLRF2008 shows that artifacts related to orbit perturbations in the LAGEOS-1/2 solutions, i.e., periods related to the draconitic years of the LAGEOS satellites, are greatly reduced in the combined solutions.

Micrometer-sized ice particles for planetary-science experiments – II. Bidirectional reflectance

We have measured the bidirectional reflectance of spherical micrometer-sized water-ice particles in the visible spectral range over a wide range of incidence and emission angles. The small ice spheres were produced by spraying fine water droplets directly into liquid nitrogen. The resulting mean particle radii are 1.47 + 0.96 - 0.58 μm. Such a material shares many properties with ice in comets and at the surface of icy satellites. Measurements show that the fresh sample material is highly backscattering, contrasting with natural terrestrial snow and frost. The formation of agglomerates of particles during the sample production results in a noticeable variability of the photometric properties of the samples in their initial state. We have also observed significant temporal evolutions of the scattering behavior of the samples, shifting towards more forward scattering within some tens of hours, resulting most likely from sintering processes. All reflectance data are fitted by the Hapke photometric model (1993 and 2002 formulation) with a one/two/three-parameter Henyey-Greenstein phase function and the resulting Hapke parameters are provided. These parameters can be used to compare laboratory results with the observed photometric behaviors of astronomical objects. We show, in particular, that the optical properties of the fresh micrometer-sized ice samples can be used to reproduce the predominant backscattering in the phase curves of Enceladus and Europa.

Calving Giant Icebergs: Old Principles, New Applications

Earth-orbiting satellites can now monitor calving of large icebergs from ice shelves bordering the marine West Antarctic Ice Sheet, and recent calving events have stimulated interest in calving mechanisms. To advance this interest pioneering work in brittle and ductile fracture mechanics is reviewed, leading to a new application to calving of giant icebergs from Antarctic ice shelves. The aim is to view iceberg calving as more than terminal events for Antarctic ice when glaciologists lose interest. Instead calving launches Antarctic ice into the larger dynamic system of Earth's climate machine. This encourages a holistic approach to glaciology.

Model-Based Estimates of Calanus Finmarchicus Abundance in the Gulf of Maine

Ocean observing systems and satellites routinely collect a wealth of information on physical conditions in the ocean. With few exceptions, such as chlorophyll concentrations, information on biological properties is harder to measure autonomously. Here, we present a system to produce estimates of the distribution and abundance of the copepod Calanus finmarchicus in the Gulf of Maine. Our system uses satellite-based measurements of sea surface temperature and chlorophyll concentration to determine the developmental and reproductive rates of C. finmarchicus. The rate information then drives a population dynamics model of C. finmarchicus that is embedded in a 2-dimensional circulation field. The first generation of this system produces realistic information on interannual variability in C. finmarchicus distribution and abundance during the winter and spring. The model can also be used to identify key drivers of interannual variability in C. finmarchicus. Experiments with the model suggest that changes in initial conditions are overwhelmed by variability in growth rates after approximately 50 d. Temperature has the largest effect on growth rate. Elevated chlorophyll during the late winter can lead to increased C. finmarchicus abundance during the spring, but the effect of variations in chlorophyll concentrations is secondary to the other inputs. Our system could be used to provide real-time estimates or even forecasts of C. finmarchicus distribution. These estimates could then be used to support management of copepod predators such as herring and right whales.

Homogeneous reprocessing of GPS, GLONASS and SLR observations

The International GNSS Service (IGS) provides operational products for the GPS and GLONASS constellation. Homogeneously processed time series of parameters from the IGS are only available for GPS. Reprocessed GLONASS series are provided only by individual Analysis Centers (i. e. CODE and ESA), making it difficult to fully include the GLONASS system into a rigorous GNSS analysis. In view of the increasing number of active GLONASS satellites and a steadily growing number of GPS+GLONASS-tracking stations available over the past few years, Technische Universität Dresden, Technische Universität München, Universität Bern and Eidgenössische Technische Hochschule Zürich performed a combined reprocessing of GPS and GLONASS observations. Also, SLR observations to GPS and GLONASS are included in this reprocessing effort. Here, we show only SLR results from a GNSS orbit validation. In total, 18 years of data (1994–2011) have been processed from altogether 340 GNSS and 70 SLR stations. The use of GLONASS observations in addition to GPS has no impact on the estimated linear terrestrial reference frame parameters. However, daily station positions show an RMS reduction of 0.3 mm on average for the height component when additional GLONASS observations can be used for the time series determination. Analyzing satellite orbit overlaps, the rigorous combination of GPS and GLONASS neither improves nor degrades the GPS orbit precision. For GLONASS, however, the quality of the microwave-derived GLONASS orbits improves due to the combination. These findings are confirmed using independent SLR observations for a GNSS orbit validation. In comparison to previous studies, mean SLR biases for satellites GPS-35 and GPS-36 could be reduced in magnitude from −35 and −38 mm to −12 and −13 mm, respectively. Our results show that remaining SLR biases depend on the satellite type and the use of coated or uncoated retro-reflectors. For Earth rotation parameters, the increasing number of GLONASS satellites and tracking stations over the past few years leads to differences between GPS-only and GPS+GLONASS combined solutions which are most pronounced in the pole rate estimates with maximum 0.2 mas/day in magnitude. At the same time, the difference between GLONASS-only and combined solutions decreases. Derived GNSS orbits are used to estimate combined GPS+GLONASS satellite clocks, with first results presented in this paper. Phase observation residuals from a precise point positioning are at the level of 2 mm and particularly reveal poorly modeled yaw maneuver periods.

The new CMOS Tracking Camera used at the Zimmerwald Observatory

During the last years the use of tracking cameras for SLR observations became less important due to the high accuracy of the predicted orbits. Upcoming new targets like satellites in eccentric orbits and space debris objects, however, require tracking cameras again. In 2013 the interline CCD camera was replaced at the Zimmerwald Observatory with a so called scientific CMOS camera. This technology promises a better performance for this application than all kinds of CCD cameras. After the comparison of the different technologies the focus will be on the integration in the Zimmerwald SLR system.

Contribution of Starlette, Stella, and AJISAI to the SLR-derived global reference frame

The contribution of Starlette, Stella, and AJI-SAI is currently neglected when defining the International Terrestrial Reference Frame, despite a long time series of precise SLR observations and a huge amount of available data. The inferior accuracy of the orbits of low orbiting geodetic satellites is the main reason for this neglect. The Analysis Centers of the International Laser Ranging Service (ILRS ACs) do, however, consider including low orbiting geodetic satellites for deriving the standard ILRS products based on LAGEOS and Etalon satellites, instead of the sparsely observed, and thus, virtually negligible Etalons. We process ten years of SLR observations to Starlette, Stella, AJISAI, and LAGEOS and we assess the impact of these Low Earth Orbiting (LEO) SLR satellites on the SLR-derived parameters. We study different orbit parameterizations, in particular different arc lengths and the impact of pseudo-stochastic pulses and dynamical orbit parameters on the quality of the solutions. We found that the repeatability of the East and North components of station coordinates, the quality of polar coordinates, and the scale estimates of the reference are improved when combining LAGEOS with low orbiting SLR satellites. In the multi-SLR solutions, the scale and the Z component of geocenter coordinates are less affected by deficiencies in solar radiation pressure modeling than in the LAGEOS-1/2 solutions, due to substantially reduced correlations between the Z geocenter coordinate and empirical orbit parameters. Eventually, we found that the standard values of Center-of-mass corrections (CoM) for geodetic LEO satellites are not valid for the currently operating SLR systems. The variations of station-dependent differential range biases reach 52 and 25 mm for AJISAI and Starlette/Stella, respectively, which is why estimating station dependent range biases or using station-dependent CoM, instead of one value for all SLR stations, is strongly recommended.This clearly indicates that the ILRS effort to produce CoM corrections for each satellite, which are site-specific and depend on the system characteristics at the time of tracking,is very important and needs to be implemented in the SLR data analysis.

Impact of Earth radiation pressure on LAGEOS orbits and on the global scale

The indirect solar radiation pressure caused by reflected or re-emitted radiation by the Earth’s surface is an important non-gravitational force perturbing the orbits of geodetic satellites (Rubincam and Weiss, 1986; Martin and Rubincam, 1996). In the case of LAGEOS this acceleration is of the order of 15% of the direct solar radiation pressure. Therefore, Earth radiation pressure has a non-negligible impact not only on LAGEOS orbits, but also on the SLR-derived terrestrial reference frame. We investigate the impact of the Earth radiation pressure on LAGEOS orbits and on the SLR-derived parameters. Earth radiation pressure has a remarkable impact on the semi-major axes of the LAGEOS satellites, causing a systematic reduction of 1.5 mm. The infrared Earth radiation causes a reduction of about 1.0 mm and the Earth’s reflectivity of 0.5 mm of the LAGEOS’ semi-major axes. The global scale defined by the SLR network is changed by 0.07 ppb, when applying Earth radiation pressure. The resulting station heights differ by 0.5-0.6 mm in the solution with and without Earth radiation pressure. However, when range biases are estimated, the height differences are absorbed by the range biases, and thus, the station heights are not shifted.

A system design for space-based space surveillance

This paper presents the capabilities of a Space-Based Space Surveillance (SBSS) demonstration mission for Space Surveillance and Tracking (SST) based on a micro- satellite platform. The results have been produced in the frame of ESA’s "As sessment Study for Space Based Space Surveillance Demonstration Mission (Phase A) " performed by the Airbus DS consortium. Space Surveillance and Tracking is part of Space Situational Awareness (SSA) and covers the detection, tracking and cataloguing of spa ce debris and satellites. Derived SST services comprise a catalogue of these man-made objects, collision warning, detection and characterisation of in-orbit fragmentations, sub-catalogue debris characterisation, etc. The assessment of SBSS in an SST system architecture has shown that both an operational SBSS and also already a well - designed space-based demonstrator can provide substantial performance in terms of surveillance and tracking of beyond - LEO objects. Especially the early deployment of a demonstrator, possible by using standard equipment, could boost initial operating capability and create a self-maintained object catalogue. Unlike classical technology demonstration missions, the primary goal is the demonstration and optimisation of the functional elements in a complex end-to-end chain (mission planning, observation strategies, data acquisition, processing and fusion, etc.) until the final products can be offered to the users. The presented SBSS system concept takes the ESA SST System Requirements (derived within the ESA SSA Preparatory Program) into account and aims at fulfilling some of the SST core requirements in a stand-alone manner. The evaluation of the concept has shown that an according solution can be implemented with low technological effort and risk. The paper presents details of the system concept, candidate micro - satellite platforms, the observation strategy and the results of performance simulations for GEO coverage and cataloguing accuracy