TanDEM-X High Resolution DEMs and Their Applications to Flow Modeling

Lava flow modeling can be a powerful tool in hazard assessments; however, the ability to produce accurate models is usually limited by a lack of high resolution, up-to-date Digital Elevation Models (DEMs). This is especially obvious in places such as Kilauea Volcano (Hawaii), where active lava flows frequently alter the terrain. In this study, we use a new technique to create high resolution DEMs on Kilauea using synthetic aperture radar (SAR) data from the TanDEM-X (TDX) satellite. We convert raw TDX SAR data into a geocoded DEM using GAMMA software [Werner et al., 2000]. This process can be completed in several hours and permits creation of updated DEMs as soon as new TDX data are available. To test the DEMs, we use the Harris and Rowland [2001] FLOWGO lava flow model combined with the Favalli et al. [2005] DOWNFLOW model to simulate the 3-15 August 2011 eruption on Kilauea's East Rift Zone. Results were compared with simulations using the older, lower resolution 2000 SRTM DEM of Hawaii. Effusion rates used in the model are derived from MODIS thermal infrared satellite imagery. FLOWGO simulations using the TDX DEM produced a single flow line that matched the August 2011 flow almost perfectly, but could not recreate the entire flow field due to the relatively high DEM noise level. The issues with short model flow lengths can be resolved by filtering noise from the DEM. Model simulations using the outdated SRTM DEM produced a flow field that followed a different trajectory to that observed. Numerous lava flows have been emplaced at Kilauea since the creation of the SRTM DEM, leading the model to project flow lines in areas that have since been covered by fresh lava flows. These results show that DEMs can quickly become outdated on active volcanoes, but our new technique offers the potential to produce accurate, updated DEMs for modeling lava flow hazards.

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.

Probabilistic approach to cloud and snow detection on Advanced Very High Resolution Radiometer (AVHRR) imagery

Derivation of probability estimates complementary to geophysical data sets has gained special attention over the last years. Information about a confidence level of provided physical quantities is required to construct an error budget of higher-level products and to correctly interpret final results of a particular analysis. Regarding the generation of products based on satellite data a common input consists of a cloud mask which allows discrimination between surface and cloud signals. Further the surface information is divided between snow and snow-free components. At any step of this discrimination process a misclassification in a cloud/snow mask propagates to higher-level products and may alter their usability. Within this scope a novel probabilistic cloud mask (PCM) algorithm suited for the 1 km × 1 km Advanced Very High Resolution Radiometer (AVHRR) data is proposed which provides three types of probability estimates between: cloudy/clear-sky, cloudy/snow and clear-sky/snow conditions. As opposed to the majority of available techniques which are usually based on the decision-tree approach in the PCM algorithm all spectral, angular and ancillary information is used in a single step to retrieve probability estimates from the precomputed look-up tables (LUTs). Moreover, the issue of derivation of a single threshold value for a spectral test was overcome by the concept of multidimensional information space which is divided into small bins by an extensive set of intervals. The discrimination between snow and ice clouds and detection of broken, thin clouds was enhanced by means of the invariant coordinate system (ICS) transformation. The study area covers a wide range of environmental conditions spanning from Iceland through central Europe to northern parts of Africa which exhibit diverse difficulties for cloud/snow masking algorithms. The retrieved PCM cloud classification was compared to the Polar Platform System (PPS) version 2012 and Moderate Resolution Imaging Spectroradiometer (MODIS) collection 6 cloud masks, SYNOP (surface synoptic observations) weather reports, Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) vertical feature mask version 3 and to MODIS collection 5 snow mask. The outcomes of conducted analyses proved fine detection skills of the PCM method with results comparable to or better than the reference PPS algorithm.

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.

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.

Lexicon manuale Hebraicum et Chaldaicum in Veteris Testamenti Libros

post ed. Germanicam 3. Latine elaboravit multisque modis retractavit et auxit Guil. Gesenius

Percent cropland coverage for sub-Saharan Africa

The dataset contains a cropland percent coverage map for Africa created through the combination of five existing land cover products: GLC-2000, MODIS Land Cover, GlobCover, MODIS Crop Likelihood and AfriCover. A synergy map was created in which the products are ranked by experts, which reflects the likelihood or probability that a given pixel is cropland. The cropland map was calibrated with national and sub-national crop statistics using a novel approach. Preliminary validation of the map was undertaken. The resulting cropland map has an accuracy of 83%, which is higher than the accuracy of any of the individual maps. The cropland percent coverage map for Africa is available for overlay on Google Earth or for download at http://agriculture.geo-wiki.org.

Particulate and phytoplankton absorption during POLARSTERN cruises ANT-XXVI/3 and ANT-XXVIII/3

The euphotic depth (Zeu) is a key parameter in modelling primary production (PP) using satellite ocean colour. However, evaluations of satellite Zeu products are scarce. The objective of this paper is to investigate existing approaches and sensors to estimate Zeu from satellite and to evaluate how different Zeu products might affect the estimation of PP in the Southern Ocean (SO). Euphotic depth was derived from MODIS and SeaWiFS products of (i) surface chlorophyll-a (Zeu-Chla) and (ii) inherent optical properties (Zeu-IOP). They were compared with in situ measurements of Zeu from different regions of the SO. Both approaches and sensors are robust to retrieve Zeu, although the best results were obtained using the IOP approach and SeaWiFS data, with an average percentage of error (E) of 25.43% and mean absolute error (MAE) of 0.10 m (log scale). Nevertheless, differences in the spatial distribution of Zeu-Chla and Zeu-IOP for both sensors were found as large as 30% over specific regions. These differences were also observed in PP. On average, PP based on Zeu-Chla was 8% higher than PP based on Zeu-IOP, but it was up to 30% higher south of 60°S. Satellite phytoplankton absorption coefficients (aph) derived by the Quasi-Analytical Algorithm at different wavelengths were also validated and the results showed that MODIS aph are generally more robust than SeaWiFS. Thus, MODIS aph should be preferred in PP models based on aph in the SO. Further, we reinforce the importance of investigating the spatial differences between satellite products, which might not be detected by the validation with in situ measurements due to the insufficient amount and uneven distribution of the data.

(Table 1) Normalized difference vegetation index in control and warming plots, Abisco Research Station

Extreme winter warming events in the sub-Arctic have caused considerable vegetation damage due to rapid changes in temperature and loss of snow cover. The frequency of extreme weather is expected to increase due to climate change thereby increasing the potential for recurring vegetation damage in Arctic regions. Here we present data on vegetation recovery from one such natural event and multiple experimental simulations in the sub-Arctic using remote sensing, handheld passive proximal sensors and ground surveys. Normalized difference vegetation index (NDVI) recovered fast (2 years), from the 26% decline following one natural extreme winter warming event. Recovery was associated with declines in dead Empetrum nigrum (dominant dwarf shrub) from ground surveys. However, E. nigrum healthy leaf NDVI was also reduced (16%) following this winter warming event in experimental plots (both control and treatments), suggesting that non-obvious plant damage (i.e., physiological stress) had occurred in addition to the dead E. nigrum shoots that was considered responsible for the regional 26% NDVI decline. Plot and leaf level NDVI provided useful additional information that could not be obtained from vegetation surveys and regional remote sensing (MODIS) alone. The major damage of an extreme winter warming event appears to be relatively transitory. However, potential knock-on effects on higher trophic levels (e.g., rodents, reindeer, and bear) could be unpredictable and large. Repeated warming events year after year, which can be expected under winter climate warming, could result in damage that may take much longer to recover.

Daily polynya area and ice production in the North Water polynya region

The North Water (NOW) Polynya is a regularly-forming area of open-water and thin-ice, located between northwestern Greenland and Ellesmere Island (Canada) at the northern tip of Baffin Bay. Due to its large spatial extent, it is of high importance for a variety of physical and biological processes, especially in wintertime. Here, we present a long-term remote sensing study for the winter seasons 1978/1979 to 2014/2015. Polynya characteristics are inferred from (1) sea ice concentrations and brightness temperatures from passive microwave satellite sensors (Advanced Microwave Scanning Radiometer (AMSR-E and AMSR2), Scanning Multichannel Microwave Radiometer (SMMR), Special Sensor Microwave Imager/Sounder (SSM/I-SSMIS)) and (2) thin-ice thickness distributions, which are calculated using MODIS ice-surface temperatures and European Center for Medium-Range Weather Forecasts (ECMWF) atmospheric reanalysis data in a 1D thermodynamic energy-balance model. Daily ice production rates are retrieved for each winter season from 2002/2003 to 2014/2015, assuming that all heat loss at the ice surface is balanced by ice growth. Two different cloud-cover correction schemes are applied on daily polynya area and ice production values to account for cloud gaps in the MODIS composites. Our results indicate that the NOW polynya experienced significant seasonal changes over the last three decades considering the overall frequency of polynya occurrences, as well as their spatial extent. In the 1980s, there were prolonged periods of a more or less closed ice cover in northern Baffin Bay in winter. This changed towards an average opening on more than 85% of the days between November and March during the last decade. Noticeably, the sea ice cover in the NOW polynya region shows signs of a later-appearing fall freeze-up, starting in the late 1990s. Different methods to obtain daily polynya area using passive microwave AMSR-E/AMSR2 data and SSM/I-SSMIS data were applied. A comparison with MODIS data (thin-ice thickness < 20 cm) shows that the wintertime polynya area estimates derived by MODIS are about 30 to 40% higher than those derived using the polynya signature simulation method (PSSM) with AMSR-E data. In turn, the difference in polynya area between PSSM and a sea ice concentration (SIC) threshold of 70% is fairly low (approximately 10%) when applied to AMSR-E data. For the coarse-resolution SSM/I-SSMIS data, this difference is much larger, particularly in November and December. Instead of a sea ice concentration threshold, the PSSM method should be used for SSM/I-SSMIS data. Depending on the type of cloud-cover correction, the calculated ice production based on MODIS data reaches an average value of 264.4 ± 65.1 km**3 to 275.7 ± 67.4 km**3 (2002/2003 to 2014/2015) and shows a high interannual variability. Our achieved long-term results underline the major importance of the NOW polynya considering its influence on Arctic ice production and associated atmosphere/ocean processes.

Pigment data set from the SeaWiFS Bio-optical Archive and Storage System (1961-1996)

The oceans play a critical role in the Earth's climate, but unfortunately, the extent of this role is only partially understood. One major obstacle is the difficulty associated with making high-quality, globally distributed observations, a feat that is nearly impossible using only ships and other ocean-based platforms. The data collected by satellite-borne ocean color instruments, however, provide environmental scientists a synoptic look at the productivity and variability of the Earth's oceans and atmosphere, respectively, on high-resolution temporal and spatial scales. Three such instruments, the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) onboard ORBIMAGE's OrbView-2 satellite, and two Moderate Resolution Imaging Spectroradiometers (MODIS) onboard the National Aeronautic and Space Administration's (NASA) Terra and Aqua satellites, have been in continuous operation since September 1997, February 2000, and June 2002, respectively. To facilitate the assembly of a suitably accurate data set for climate research, members of the NASA Sensor Intercomparison and Merger for Biological and Interdisciplinary Oceanic Studies (SIMBIOS) Project and SeaWiFS Project Offices devote significant attention to the calibration and validation of these and other ocean color instruments. This article briefly presents results from the SIMBIOS and SeaWiFS Project Office's (SSPO) satellite ocean color validation activities and describes the SeaWiFS Bio-optical Archive and Storage System (SeaBASS), a state-of-the-art system for archiving, cataloging, and distributing the in situ data used in these activities.

Secchi depth maps for the Baltic Sea, 2003-2012, with links to GeoTIFs

Secchi depth is a measure of water transparency. In the Baltic Sea region, Secchi depth maps are used to assess eutrophication and as input for habitat models. Due to their spatial and temporal coverage, satellite data would be the most suitable data source for such maps. But the Baltic Sea's optical properties are so different from the open ocean that globally calibrated standard models suffer from large errors. Regional predictive models that take the Baltic Sea's special optical properties into account are thus needed. This paper tests how accurately generalized linear models (GLMs) and generalized additive models (GAMs) with MODIS/Aqua and auxiliary data as inputs can predict Secchi depth at a regional scale. It uses cross-validation to test the prediction accuracy of hundreds of GAMs and GLMs with up to 5 input variables. A GAM with 3 input variables (chlorophyll a, remote sensing reflectance at 678 nm, and long-term mean salinity) made the most accurate predictions. Tested against field observations not used for model selection and calibration, the best model's mean absolute error (MAE) for daily predictions was 1.07 m (22%), more than 50% lower than for other publicly available Baltic Sea Secchi depth maps. The MAE for predicting monthly averages was 0.86 m (15%). Thus, the proposed model selection process was able to find a regional model with good prediction accuracy. It could be useful to find predictive models for environmental variables other than Secchi depth, using data from other satellite sensors, and for other regions where non-standard remote sensing models are needed for prediction and mapping. Annual and monthly mean Secchi depth maps for 2003-2012 come with this paper as Supplementary materials.

Radiation and meteorological observations at Neumay Station (2002-2014)

Abstract has to submitted by author

Accumulation, velocities, fluxes and mass balance of glacial tributaries to the Lambert Glacier-Amery Ice Shelf system, Antarctica

By incorporating recently available remote sensing data, we investigated the mass balance for all individual tributary glacial basins of the Lambert Glacier-Amery Ice Shelf system, East Antarctica. On the basis of the ice flow information derived from SAR interferometry and ICESat laser altimetry, we have determined the spatial configuration of eight tributary drainage basins of the Lambert-Amery glacial system. By combining the coherence information from SAR interferometry and the texture information from SAR and MODIS images, we have interpreted and refined the grounding line position. We calculated ice volume flux of each tributary glacial basin based on the ice velocity field derived from Radarsat three-pass interferometry together with ice thickness data interpolated from Australian and Russian airborne radio echo sounding (RES) surveys and inferred from ICESat laser altimetry data. Our analysis reveals that three tributary basins have a significant net positive imbalance, while five other subbasins are slightly positive or close to zero balance. Overall, in contrast to previous studies, we find that the grounded ice in Lambert Glacier-Amery Ice Shelf system has a positive mass imbalance of 22.9 ± 4.4 Gt/a. The net basal melting for the entire Amery Ice Shelf is estimated to be 27.0 ± 7.0 Gt/a. The melting rate decreases rapidly from the grounding zone to the ice shelf front. Significant basal refreezing is detected in the downstream section of the ice shelf. The mass balance estimates for both the grounded ice sheet and the ice shelf mass differ substantially from other recent estimates.

(Table 1) Ice thickness at lake Nam Co, Tibet, during 2008-2009

Lake ice change is one of the sensitive indicators of regional and global climate change. Different sources of data are used in monitoring lake ice phenology nowadays. Visible and Near Infrared bands of imagery (VNIR) are well suited for the observation of freshwater ice change, for example data from AVHRR and MODIS. Active and passive microwave data are also used for the observation of lake ice, e.g., from satellite altimetry and radiometry, backscattering coefficient from QuickSCAT, brightness temperature (Tb) from SSM/I, SMMR, and AMSR-E. Most of the studies are about lake ice cover phenology, while few studies focus on lake ice thickness. For example, Hall et al. using 5 GHz (6 cm) radiometer data showed a good relationship between Tb and ice thickness. Kang et al. found the seasonal evolution of Tb at 10.65 GHz and 18.7 GHz from AMSR-E to be strongly influenced by ice thickness. Many studies on lake ice phenology have been carried out since the 1970s in cold regions, especially in Canada, the USA, Europe, the Arctic, and Antarctica. However, on the Tibetan Plateau, very little research has focused on lake ice-cover change; only a small number of published papers on Qinghai Lake ice observations. The main goal of this study is to investigate the change in lake ice phenology at Nam Co on the Tibetan Plateau using MODIS and AMSR-E data (monitoring the date of freeze onset, the formation of stable ice cover, first appearance of water, and the complete disappearance of ice) during the period 2000-2009.