EXPLORATION OF THE MTSAT2 SATELLITE CAPABILITIES FOR REAL TIME DETECTION AND CHARACTERIZATION OF VOLCANIC EMISSIONS

In this report, we attempt to define the capabilities of the infrared satellite remote sensor, Multifunctional Transport Satellite-2 (MTSAT-2) (i.e. a geosynchronous instrument), in characterizing volcanic eruptive behavior in the highly active region of Indonesia. Sulfur dioxide data from NASA's Ozone Monitoring Instrument (OMI) (i.e. a polar orbiting instrument) are presented here for validation of the processes interpreted using the thermal infrared datasets. Data provided from two case studies are analyzed specifically for eruptive products producing large thermal anomalies (i.e. lava flows, lava domes, etc.), volcanic ash and SO2 clouds; three distinctly characteristic and abundant volcanic emissions. Two primary methods used for detection of heat signatures are used and compared in this report including, single-channel thermal radiance (4-µm) and the normalized thermal index (NTI) algorithm. For automated purposes, fixed thresholds must be determined for these methods. A base minimum detection limit (MDL) for single-channel thermal radiance of 2.30E+05 Wm- 2sr-1m-1 and -0.925 for NTI generate false alarm rates of 35.78% and 34.16%, respectively. A spatial comparison method, developed here specifically for use in Indonesia and used as a second parameter for detection, is implemented to address the high false alarm rate. For the single-channel thermal radiance method, the utilization of the spatial comparison method eliminated 100% of the false alarms while maintaining every true anomaly. The NTI algorithm showed similar results with only 2 false alarms remaining. No definitive difference is observed between the two thermal detection methods for automated use; however, the single-channel thermal radiance method coupled with the SO2 mass abundance data can be used to interpret volcanic processes including the identification of lava dome activity at Sinabung as well as the mechanism for the dome emplacement (i.e. endogenous or exogenous). Only one technique, the brightness temperature difference (BTD) method, is used for the detection of ash. Trends of ash area, water/ice area, and their respective concentrations yield interpretations of increased ice formation, aggregation, and sedimentation processes that only a high-temporal resolution instrument like the MTSAT-2 can analyze. A conceptual model of a secondary zone of aggregation occurring in the migrating Kelut ash cloud, which decreases the distal fine-ash component and hazards to flight paths, is presented in this report. Unfortunately, SO2 data was unable to definitively reinforce the concept of a secondary zone of aggregation due to the lack of a sufficient temporal resolution. However, a detailed study of the Kelut SO2 cloud is used to determine that there was no climatic impacts generated from this eruption due to the atmospheric residence times and e-folding rate of ~14 days for the SO2. This report applies the complementary assets offered by utilizing a high-temporal and a high-spatial resolution satellite, and it demonstrates that these two instruments can provide unparalleled observations of dynamic volcanic processes.

Development and characterization of fluorescent pH sensors based on porous silica and hydrogel support matrices

The hydrogen ion activity (pH) is a very important parameter in environment monitoring, biomedical research and other applications. Optical pH sensors have several advantages over traditional potentiometric pH measurement, such as high sensitivity, no need of constant calibration, easy for miniaturization and possibility for remote sensing. Several pH indicators has been successfully immobilized in three different solid porous materials to use as pH sensing probes. The fluorescent pH indicator fluorescein-5-isothiocyanate (FITC) was covalently bound onto the internal surface of porous silica (pore size ~10 nm) and retained its pH sensitivity. The excited state pK* a of FITC in porous silica (5.58) was slightly smaller than in solution (5.68) due to the free silanol groups (Si-OH) on the silica surface. The pH sensitive range for this probe is pH 4.5 - 7.0 with an error less than 0.1 pH units. The probe response was reproducible and stable for at least four month, stored in DI water, but exhibit a long equilibrium of up to 100 minutes. Sol-gel based pH sensors were developed with immobilization of two fluorescent pH indicators fluorescein-5-(and-6)-sulfonic acid, trisodium salt (FS) and 8-hydroxypyrene- 1,3,6-trisulfonic acid (HPTS) through physical entrapment. Prior to immobilization, the indicators were ion-paired with a common surfactant hexadecyltrimethylammonium bromide (CTAB) in order to prevent leaching. The sol-gel films were synthesized through the hydrolysis of two different precursors, ethyltriethoxysilane (ETEOS) and 3- glycidoxypropyltrimethoxysilane (GPTMS) and deposited on a quartz slide through spin coating. The pK a of the indicators immobilized in sol-gel films was much smaller than in solutions due to silanol groups on the inner surface of the sol-gel films and ammonium groups from the surrounding surfactants. Unlike in solution, the apparent pK a of the indicators in sol-gel films increased with increasing ionic strength. The equilibrium time for these sensors was within 5 minutes (with film thickness of ~470 nm). Polyethylene glycol (PEG) hydrogel was of interest for optical pH sensor development because it is highly proton permeable, transparent and easy to synthesize. pH indicators can be immobilized in hydrogel through physical entrapment and copolymerization. FS and HPTS ion-pairs were physically entrapped in hydrogel matrix synthesized via free radical initiation. For covalent immobilization, three indicators, 6,8-dihydroxypyrene-1,3- disulfonic acid (DHPDS), 2,7-dihydroxynaphthalene-3,6-disulfonic acid (DHNDS) and cresol red were first reacted with methacrylic anhydride (MA) to form methacryloylanalogs for copolymerization. These hydrogels were synthesized in aqueous solution with a redox initiation system. The thickness of the hydrogel film is controlled as ~ 0.5 cm and the porosity can be adjusted with the percentage of polyethylene glycol in the precursor solutions. The pK a of the indicators immobilized in the hydrogel both physically and covalently were higher than in solution due to the medium effect. The sensors are stable and reproducible with a short equilibrium time (less than 4 minutes). In addition, the color change of cresol red immobilized hydrogel is vivid from yellow (acidic condition) to purple (basic condition). Due to covalently binding, cresol red was not leaching out from the hydrogel, making it a good candidate of reusable "pH paper".

Methane emissions from floodplains in the Amazon Basin: challenges in developing a process-based model for global applications

Tropical wetlands are estimated to represent about 50% of the natural wetland methane (CH4) emissions and explain a large fraction of the observed CH4 variability on timescales ranging from glacial–interglacial cycles to the currently observed year-to-year variability. Despite their importance, however, tropical wetlands are poorly represented in global models aiming to predict global CH4 emissions. This publication documents a first step in the development of a process-based model of CH4 emissions from tropical floodplains for global applications. For this purpose, the LPX-Bern Dynamic Global Vegetation Model (LPX hereafter) was slightly modified to represent floodplain hydrology, vegetation and associated CH4 emissions. The extent of tropical floodplains was prescribed using output from the spatially explicit hydrology model PCR-GLOBWB. We introduced new plant functional types (PFTs) that explicitly represent floodplain vegetation. The PFT parameterizations were evaluated against available remote-sensing data sets (GLC2000 land cover and MODIS Net Primary Productivity). Simulated CH4 flux densities were evaluated against field observations and regional flux inventories. Simulated CH4 emissions at Amazon Basin scale were compared to model simulations performed in the WETCHIMP intercomparison project. We found that LPX reproduces the average magnitude of observed net CH4 flux densities for the Amazon Basin. However, the model does not reproduce the variability between sites or between years within a site. Unfortunately, site information is too limited to attest or disprove some model features. At the Amazon Basin scale, our results underline the large uncertainty in the magnitude of wetland CH4 emissions. Sensitivity analyses gave insights into the main drivers of floodplain CH4 emission and their associated uncertainties. In particular, uncertainties in floodplain extent (i.e., difference between GLC2000 and PCR-GLOBWB output) modulate the simulated emissions by a factor of about 2. Our best estimates, using PCR-GLOBWB in combination with GLC2000, lead to simulated Amazon-integrated emissions of 44.4 ± 4.8 Tg yr−1. Additionally, the LPX emissions are highly sensitive to vegetation distribution. Two simulations with the same mean PFT cover, but different spatial distributions of grasslands within the basin, modulated emissions by about 20%. Correcting the LPX-simulated NPP using MODIS reduces the Amazon emissions by 11.3%. Finally, due to an intrinsic limitation of LPX to account for seasonality in floodplain extent, the model failed to reproduce the full dynamics in CH4 emissions but we proposed solutions to this issue. The interannual variability (IAV) of the emissions increases by 90% if the IAV in floodplain extent is accounted for, but still remains lower than in most of the WETCHIMP models. While our model includes more mechanisms specific to tropical floodplains, we were unable to reduce the uncertainty in the magnitude of wetland CH4 emissions of the Amazon Basin. Our results helped identify and prioritize directions towards more accurate estimates of tropical CH4 emissions, and they stress the need for more research to constrain floodplain CH4 emissions and their temporal variability, even before including other fundamental mechanisms such as floating macrophytes or lateral water fluxes.

The International Trans-Antarctic Scientific Expedition (ITASE): An Overview

From its original formulation in 1990 the International Trans-Antarctic Scientific Expedition (ITASE) has had as its primary aim the collection and interpretation of a continent-wide array of environmental parameters assembled through the coordinated efforts of scientists from several nations. ITASE offers the ground-based opportunities of traditional-style traverse travel coupled with the modern technology of CPS, crevasse detecting radar, satellite communications and multidisciplinary research. By operating predominantly in the mode of an oversnow traverse, ITASE offers scientists the opportunity to experience the dynamic range of the Antarctic environment. ITASE also offers an important interactive venue for research similar to that afforded by oceanographic research vessels and large polar field camps, without the cost of the former or the lack of mobility of the latter. More importantly, the combination of disciplines represented by ITASE provides a unique, multidimensional (space and time) view of the ice sheet and its history. ITASE has now collected > 20 000 km of snow radar, recovered more than 240 firn/ice cores (total length 7000m), remotely penetrated to similar to 4000m into the ice sheet, and sampled the atmosphere to heights of > 20 km.

Daytime Low Stratiform Cloud Detection on AVHRR Imagery

The near-real time retrieval of low stratiform cloud (LSC) coverage is of vital interest for such disciplines as meteorology, transport safety, economy and air quality. Within this scope, a novel methodology is proposed which provides the LSC occurrence probability estimates for a satellite scene. The algorithm is suited for the 1 × 1 km Advanced Very High Resolution Radiometer (AVHRR) data and was trained and validated against collocated SYNOP observations. Utilisation of these two combined data sources requires a formulation of constraints in order to discriminate cases where the LSC is overlaid by higher clouds. The LSC classification process is based on six features which are first converted to the integer form by step functions and combined by means of bitwise operations. Consequently, a set of values reflecting a unique combination of those features is derived which is further employed to extract the LSC occurrence probability estimates from the precomputed look-up vectors (LUV). Although the validation analyses confirmed good performance of the algorithm, some inevitable misclassification with other optically thick clouds were reported. Moreover, the comparison against Polar Platform System (PPS) cloud-type product revealed superior classification accuracy. From the temporal perspective, the acquired results reported a presence of diurnal and annual LSC probability cycles over Europe.

Sideband-Ratio and Harmonic Response of Sub-Millimeter Wave Mixers measured with a Martin Puplett Interferometer

Heterodyne receivers at millimeter and submillimeter wavelength are widely used for radiometric spectral line observations for atmospheric remote sensing or radio astronomy. The quantitative analysis of such observations requires an accurate knowledge of the mixers's sideband ratio. In addition, its potential sensitivity to spurious harmonics needs to be well understood. In this paper, we discuss a measurement technique for these receiver characteristics, which is based on a scanning Martin Puplett Interferometer used in conjunction with a wide band digital autocorrelation spectrometer for the analysis of the intermediate frequency band. We present measurement results of different double sideband and sideband separating mixers, which were developed for the proposed 340GHz multi-beam limb sounder STEAMR.

Lake surface water temperatures of European Alpine lakes (1989–2013) based on the Advanced Very High Resolution Radiometer (AVHRR) 1 km data set

Lake water temperature (LWT) is an important driver of lake ecosystems and it has been identified as an indicator of climate change. Consequently, the Global Climate Observing System (GCOS) lists LWT as an essential climate variable. Although for some European lakes long in situ time series of LWT do exist, many lakes are not observed or only on a non-regular basis making these observations insufficient for climate monitoring. Satellite data can provide the information needed. However, only few satellite sensors offer the possibility to analyse time series which cover 25 years or more. The Advanced Very High Resolution Radiometer (AVHRR) is among these and has been flown as a heritage instrument for almost 35 years. It will be carried on for at least ten more years, offering a unique opportunity for satellite-based climate studies. Herein we present a satellite-based lake surface water temperature (LSWT) data set for European water bodies in or near the Alps based on the extensive AVHRR 1 km data record (1989–2013) of the Remote Sensing Research Group at the University of Bern. It has been compiled out of AVHRR/2 (NOAA-07, -09, -11, -14) and AVHRR/3 (NOAA-16, -17, -18, -19 and MetOp-A) data. The high accuracy needed for climate related studies requires careful pre-processing and consideration of the atmospheric state. The LSWT retrieval is based on a simulation-based scheme making use of the Radiative Transfer for TOVS (RTTOV) Version 10 together with ERA-interim reanalysis data from the European Centre for Medium-range Weather Forecasts. The resulting LSWTs were extensively compared with in situ measurements from lakes with various sizes between 14 and 580 km2 and the resulting biases and RMSEs were found to be within the range of −0.5 to 0.6 K and 1.0 to 1.6 K, respectively. The upper limits of the reported errors could be rather attributed to uncertainties in the data comparison between in situ and satellite observations than inaccuracies of the satellite retrieval. An inter-comparison with the standard Moderate-resolution Imaging Spectroradiometer (MODIS) Land Surface Temperature product exhibits RMSEs and biases in the range of 0.6 to 0.9 and −0.5 to 0.2 K, respectively. The cross-platform consistency of the retrieval was found to be within ~ 0.3 K. For one lake, the satellite-derived trend was compared with the trend of in situ measurements and both were found to be similar. Thus, orbital drift is not causing artificial temperature trends in the data set. A comparison with LSWT derived through global sea surface temperature (SST) algorithms shows lower RMSEs and biases for the simulation-based approach. A running project will apply the developed method to retrieve LSWT for all of Europe to derive the climate signal of the last 30 years. The data are available at doi:10.1594/PANGAEA.831007.

Classifying and mapping the urban transition in Vietnam

The urban transition almost always involves wrenching social adjustment as small agricultural communities are forced to adjust rapidly to industrial ways of life. Large-scale in-migration of young people, usually from poor regions, creates enormous demand and expectations for community and social services. One immediate problem planners face in approaching this challenge is how to define, differentiate, and map what is rural, urban, and transitional (i.e., peri-urban). This project established an urban classification for Vietnam by using national census and remote sensing data to identify and map the smallest administrative units for which data are collected as rural, peri-urban, urban, or urban core. We used both natural and human factors in the quantitative model: income from agriculture, land under agriculture and forests, houses with modern sanitation, and the Normalized Difference Vegetation Index. Model results suggest that in 2006, 71% of Vietnam's 10,891 communes were rural, 18% peri-urban, 3% urban, and 4% urban core. Of the communes our model classified as peri-urban, 61% were classified by the Vietnamese government as rural. More than 7% of Vietnam's land area can be classified as peri-urban and approximately 13% of its population (more than 11 million people) lives in peri-urban areas. We identified and mapped three types of peri-urban places: communes in the periphery of large towns and cities; communes along highways; and communes associated with provincial administration or home to industrial, energy, or natural resources projects (e.g., mining). We validated this classification based on ground observations, analyses of multi-temporal night-time lights data, and an examination of road networks. The model provides a method for rapidly assessing the rural–urban nature of places to assist planners in identifying rural areas undergoing rapid change with accompanying needs for investments in building, sanitation, road infrastructure, and government institutions.

The rotation state of 67P/Churyumov-Gerasimenko from approach observations with the OSIRIS cameras on Rosetta

Aims. Approach observations with the Optical, Spectroscopic, and Infrared Remote Imaging System (OSIRIS) experiment onboard Rosetta are used to determine the rotation period, the direction of the spin axis, and the state of rotation of comet 67P’s nucleus. Methods. Photometric time series of 67P have been acquired by OSIRIS since the post wake-up commissioning of the payload in March 2014. Fourier analysis and convex shape inversion methods have been applied to the Rosetta data as well to the available ground-based observations. Results. Evidence is found that the rotation rate of 67P has significantly changed near the time of its 2009 perihelion passage, probably due to sublimation-induced torque. We find that the sidereal rotation periods P1 = 12.76129 ± 0.00005 h and P2 = 12.4043 ± 0.0007 h for the apparitions before and after the 2009 perihelion, respectively, provide the best fit to the observations. No signs of multiple periodicity are found in the light curves down to the noise level, which implies that the comet is presently in a simple rotation state around its axis of largest moment of inertia. We derive a prograde rotation model with spin vector J2000 ecliptic coordinates λ = 65° ± 15°, β = + 59° ± 15°, corresponding to equatorial coordinates RA = 22°, Dec = + 76°. However, we find that the mirror solution, also prograde, at λ = 275° ± 15°, β = + 50° ± 15° (or RA = 274°, Dec = + 27°), is also possible at the same confidence level, due to the intrinsic ambiguity of the photometric problem for observations performed close to the ecliptic plane.

Trajectory mapping of middle atmospheric water vapor by a mini network of NDACC instruments

The important task to observe the global coverage of middle atmospheric trace gases like water vapor or ozone usually is accomplished by satellites. Climate and atmospheric studies rely upon the knowledge of trace gas distributions throughout the stratosphere and mesosphere. Many of these gases are currently measured from satellites, but it is not clear whether this capability will be maintained in the future. This could lead to a significant knowledge gap of the state of the atmosphere. We explore the possibilities of mapping middle atmospheric water vapor in the Northern Hemisphere by using Lagrangian trajectory calculations and water vapor profile data from a small network of five ground-based microwave radiometers. Four of them are operated within the frame of NDACC (Network for the Detection of Atmospheric Composition Change). Keeping in mind that the instruments are based on different hardware and calibration setups, a height-dependent bias of the retrieved water vapor profiles has to be expected among the microwave radiometers. In order to correct and harmonize the different data sets, the Microwave Limb Sounder (MLS) on the Aura satellite is used to serve as a kind of traveling standard. A domain-averaging TM (trajectory mapping) method is applied which simplifies the subsequent validation of the quality of the trajectory-mapped water vapor distribution towards direct satellite observations. Trajectories are calculated forwards and backwards in time for up to 10 days using 6 hourly meteorological wind analysis fields. Overall, a total of four case studies of trajectory mapping in different meteorological regimes are discussed. One of the case studies takes place during a major sudden stratospheric warming (SSW) accompanied by the polar vortex breakdown; a second takes place after the reformation of stable circulation system. TM cases close to the fall equinox and June solstice event from the year 2012 complete the study, showing the high potential of a network of ground-based remote sensing instruments to synthesize hemispheric maps of water vapor.

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.

Trajectory mapping of middle atmospheric water vapor by a mini network of NDACC instruments

The important task to observe the global coverage of middle atmospheric trace gases like water vapor or ozone usually is accomplished by satellites. Climate and atmospheric studies rely upon the knowledge of trace gas distributions throughout the stratosphere and mesosphere. Many of these gases are currently measured from satellites, but it is not clear whether this capability will be maintained in the future. This could lead to a significant knowledge gap of the state of the atmosphere. We explore the possibilities of mapping middle atmospheric water vapor in the Northern Hemisphere by using Lagrangian trajectory calculations and water vapor profile data from a small network of five ground-based microwave radiometers. Four of them are operated within the frame of NDACC (Network for the Detection of Atmospheric Composition Change). Keeping in mind that the instruments are based on different hardware and calibration setups, a height-dependent bias of the retrieved water vapor profiles has to be expected among the microwave radiometers. In order to correct and harmonize the different data sets, the Microwave Limb Sounder (MLS) on the Aura satellite is used to serve as a kind of traveling standard. A domain-averaging TM (trajectory mapping) method is applied which simplifies the subsequent validation of the quality of the trajectory-mapped water vapor distribution towards direct satellite observations. Trajectories are calculated forwards and backwards in time for up to 10 days using 6 hourly meteorological wind analysis fields. Overall, a total of four case studies of trajectory mapping in different meteorological regimes are discussed. One of the case studies takes place during a major sudden stratospheric warming (SSW) accompanied by the polar vortex breakdown; a second takes place after the reformation of stable circulation system. TM cases close to the fall equinox and June solstice event from the year 2012 complete the study, showing the high potential of a network of ground-based remote sensing instruments to synthesize hemispheric maps of water vapor.