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.

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.

Validation of a modified AVHRR aerosol optical depth retrieval algorithm over Central Europe

The Advanced Very High Resolution Radiometer (AVHRR) carried on board the National Oceanic and Atmospheric Administration (NOAA) and the Meteorological Operational Satellite (MetOp) polar orbiting satellites is the only instrument offering more than 25 years of satellite data to analyse aerosols on a daily basis. The present study assessed a modified AVHRR aerosol optical depth τa retrieval over land for Europe. The algorithm might also be applied to other parts of the world with similar surface characteristics like Europe, only the aerosol properties would have to be adapted to a new region. The initial approach used a relationship between Sun photometer measurements from the Aerosol Robotic Network (AERONET) and the satellite data to post-process the retrieved τa. Herein a quasi-stand-alone procedure, which is more suitable for the pre-AERONET era, is presented. In addition, the estimation of surface reflectance, the aerosol model, and other processing steps have been adapted. The method's cross-platform applicability was tested by validating τa from NOAA-17 and NOAA-18 AVHRR at 15 AERONET sites in Central Europe (40.5° N–50° N, 0° E–17° E) from August 2005 to December 2007. Furthermore, the accuracy of the AVHRR retrieval was related to products from two newer instruments, the Medium Resolution Imaging Spectrometer (MERIS) on board the Environmental Satellite (ENVISAT) and the Moderate Resolution Imaging Spectroradiometer (MODIS) on board Aqua/Terra. Considering the linear correlation coefficient R, the AVHRR results were similar to those of MERIS with even lower root mean square error RMSE. Not surprisingly, MODIS, with its high spectral coverage, gave the highest R and lowest RMSE. Regarding monthly averaged τa, the results were ambiguous. Focusing on small-scale structures, R was reduced for all sensors, whereas the RMSE solely for MERIS substantially increased. Regarding larger areas like Central Europe, the error statistics were similar to the individual match-ups. This was mainly explained with sampling issues. With the successful validation of AVHRR we are now able to concentrate on our large data archive dating back to 1985. This is a unique opportunity for both climate and air pollution studies over land surfaces.

A satellite-based snow cover climatology (1985–2011) for the European Alps derived from AVHRR data

Seasonal snow cover is of great environmental and socio-economic importance for the European Alps. Therefore a high priority has been assigned to quantifying its temporal and spatial variability. Complementary to land-based monitoring networks, optical satellite observations can be used to derive spatially comprehensive information on snow cover extent. For understanding long-term changes in alpine snow cover extent, the data acquired by the Advanced Very High Resolution Radiometer (AVHRR) sensors mounted onboard the National Oceanic and Atmospheric Association (NOAA) and Meteorological Operational satellite (MetOp) platforms offer a unique source of information. In this paper, we present the first space-borne 1 km snow extent climatology for the Alpine region derived from AVHRR data over the period 1985–2011. The objective of this study is twofold: first, to generate a new set of cloud-free satellite snow products using a specific cloud gap-filling technique and second, to examine the spatiotemporal distribution of snow cover in the European Alps over the last 27 yr from the satellite perspective. For this purpose, snow parameters such as snow onset day, snow cover duration (SCD), melt-out date and the snow cover area percentage (SCA) were employed to analyze spatiotemporal variability of snow cover over the course of three decades. On the regional scale, significant trends were found toward a shorter SCD at lower elevations in the south-east and south-west. However, our results do not show any significant trends in the monthly mean SCA over the last 27 yr. This is in agreement with other research findings and may indicate a deceleration of the decreasing snow trend in the Alpine region. Furthermore, such data may provide spatially and temporally homogeneous snow information for comprehensive use in related research fields (i.e., hydrologic and economic applications) or can serve as a reference for climate models.

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.

Lake ice phenology from AVHRR data for European lakes: An automated two-step extraction method

Several lake ice phenology studies from satellite data have been undertaken. However, the availability of long-term lake freeze-thaw-cycles, required to understand this proxy for climate variability and change, is scarce for European lakes. Long time series from space observations are limited to few satellite sensors. Data of the Advanced Very High Resolution Radiometer (AVHRR) are used in account of their unique potential as they offer each day global coverage from the early 1980s expectedly until 2022. An automatic two-step extraction was developed, which makes use of near-infrared reflectance values and thermal infrared derived lake surface water temperatures to extract lake ice phenology dates. In contrast to other studies utilizing thermal infrared, the thresholds are derived from the data itself, making it unnecessary to define arbitrary or lake specific thresholds. Two lakes in the Baltic region and a steppe lake on the Austrian–Hungarian border were selected. The later one was used to test the applicability of the approach to another climatic region for the time period 1990 to 2012. A comparison of the extracted event dates with in situ data provided good agreements of about 10 d mean absolute error. The two-step extraction was found to be applicable for European lakes in different climate regions and could fill existing data gaps in future applications. The extension of the time series to the full AVHRR record length (early 1980 until today) with adequate length for trend estimations would be of interest to assess climate variability and change. Furthermore, the two-step extraction itself is not sensor-specific and could be applied to other sensors with equivalent near- and thermal infrared spectral bands.

Spatial and temporal patterns of greenness on the Yamal Peninsula, Russia: interactions of ecological and social factors affecting the Arctic normalized difference vegetation index

The causes of a greening trend detected in the Arctic using the normalized difference vegetation index (NDVI) are still poorly understood. Changes in NDVI are a result of multiple ecological and social factors that affect tundra net primary productivity. Here we use a 25 year time series of AVHRR-derived NDVI data (AVHRR: advanced very high resolution radiometer), climate analysis, a global geographic information database and ground-based studies to examine the spatial and temporal patterns of vegetation greenness on the Yamal Peninsula, Russia. We assess the effects of climate change, gas-field development, reindeer grazing and permafrost degradation. In contrast to the case for Arctic North America, there has not been a significant trend in summer temperature or NDVI, and much of the pattern of NDVI in this region is due to disturbances. There has been a 37% change in early-summer coastal sea-ice concentration, a 4% increase in summer land temperatures and a 7% change in the average time-integrated NDVI over the length of the satellite observations. Gas-field infrastructure is not currently extensive enough to affect regional NDVI patterns. The effect of reindeer is difficult to quantitatively assess because of the lack of control areas where reindeer are excluded. Many of the greenest landscapes on the Yamal are associated with landslides and drainage networks that have resulted from ongoing rapid permafrost degradation. A warming climate and enhanced winter snow are likely to exacerbate positive feedbacks between climate and permafrost thawing. We present a diagram that summarizes the social and ecological factors that influence Arctic NDVI. The NDVI should be viewed as a powerful monitoring tool that integrates the cumulative effect of a multitude of factors affecting Arctic land-cover change.

3D documentation of footwear impressions and tyre tracks in snow with high resolution optical surface scanning

The three-dimensional documentation of footwear and tyre impressions in snow offers an opportunity to capture additional fine detail for the identification as present photographs. For this approach, up to now, different casting methods have been used. Casting of footwear impressions in snow has always been a difficult assignment. This work demonstrates that for the three-dimensional documentation of impressions in snow the non-destructive method of 3D optical surface scanning is suitable. The new method delivers more detailed results of higher accuracy than the conventional casting techniques. The results of this easy to use and mobile 3D optical surface scanner were very satisfactory in different meteorological and snow conditions. The method is also suitable for impressions in soil, sand or other materials. In addition to the side by side comparison, the automatic comparison of the 3D models and the computation of deviations and accuracy of the data simplify the examination and delivers objective and secure results. The results can be visualized efficiently. Data exchange between investigating authorities at a national or an international level can be achieved easily with electronic data carriers.

New data on the Lateglacial period of SW Europe: a high resolution multiproxy record from Laguna de la Roya (NW Iberia)

High-resolution multiproxy analyses were performed on a 128 cm section of organic sediments accumulated in a small mountain lake in NW Iberia (Laguna de la Roya, 1608 m asl). The pollen stratigraphy together with radiocarbon dating provided the basis for a chronology ranging from 15,600 to 10,500 cal yr BP. Chironomid-inferred July air temperatures suggest a temperature range from 7 to 13 degrees C, also evidencing two well-established cold periods which may be equivalent to the INTIMATE stages GS-2a and GS-1. Furthermore, a number of short cold events (with summer temperatures dropping about 0.5-1 degrees C) appear intercalated within the Lateglacial Interstadial (possibly equivalent to the INTIMATE cold events GI-1d, GI-1c2 and GI-1b) and the early Holocene (possibly equivalent to the 11.2 k event). The temperature variations predicted by our reconstruction allow explaining the changes in local conditions and productivity of the lake inferred from the biological record of the same sediment core. Furthermore, they also agree with the local and regional vegetation dynamics, and the main oscillations deduced for the vegetation belts. Based on its chronology our multiproxy record indicates a similar temperature development in NW Iberia as inferred by the Greenland delta O-18 record, the marine deep-sea records off the Atlantic Iberian Margin, and other chironomid-based Lateglacial temperature reconstructions from Europe. Nevertheless, the impact of most of the less intense Lateglacial/early Holocene cold events in NW Iberia was most probably limited to very sensitive sites that were very close to ecotonal situations. Particularly, our new pollen record indicates that they were represented as three minor environmental crises occurring during the Lateglacial Interestadial in this area. The Older Dryas event (in our usage corresponding to the Aegelsee Oscillation in Central Europe and event GI-1d in central Greenland) has previously been described in this region, but its age and duration (ca 14,250-14050 cal yr BP) is now better constrained. The two subsequent stages, La Roya I (ca 13,600-13,400 cal yr BP) and La Roya II (ca 13,300-12,900 cal yr BP) have been described for first time in NW Iberia. (C) 2013 Elsevier Ltd. All rights reserved.

Quantification of coronary artery stenosis with high-resolution CT in comparison with histopathology in an ex vivo study

PURPOSE To investigate the ex vivo performance of high-resolution computed tomography (CT) for quantitative assessment of percentage diameter stenosis in coronary arteries compared to histopathology. MATERIALS AND METHODS High-resolution CT was performed in 26 human heart specimens after the injection of iodinated contrast media into the coronary arteries. Coronary artery plaques were visually identified on CT images and the grade of stenosis for each plaque was measured with electronic calipers. All coronary plaques were characterized by histopathology according to the Stary classification, and the percentage of stenosis was measured. RESULTS CT depicted 84% (274/326) of all coronary plaques identified by histology. Missed plaques by CT were of Stary type I (n=31), type II (n=16), and type III (n=5). The stenosis degree significantly correlated between CT and histology (r=0.81, p<0.001). CT systematically overestimated the stenosis of calcified plaques (mean difference - 11.0 ± 9.5%, p<0.01) and systematically underestimated the stenosis of non-calcified plaques (mean difference -6.8 ± 10.4%, p<0.05), while there was no significant difference for mixed-type plaques (mean difference -0.4 ± 11.7%, p=0.85). There was a significant underestimation of stenosis degree as measured by CT for Stary II plaques (mean difference -14 ± 9%, p<0.01) and a significant overestimation for Stary VII plaques (mean difference 9 ± 10%, p<0.05), but there was no significant difference in stenosis degree between both modalities for other plaque types. CONCLUSIONS High-resolution CT reliably depicts advanced stage coronary plaques with an overall good correlation of stenosis degree compared to histology, however, the degree of stenosis is systematically overestimated in calcified and underestimated in non-calcified plaques.

Self-supporting CVD diamond charge state conversion surfaces for high resolution imaging of low-energy neutral atoms in space plasmas

Two polycrystalline diamond surfaces, manufactured by chemical vapour deposition (CVD) technique, are investigated regarding their applicability as charge state conversion surfaces (CS) for use in a low energy neutral atom imaging instrument in space research. The capability of the surfaces for converting neutral atoms into negative ions via surface ionisation processes was measured for hydrogen and oxygen with particle energies in the range from 100 eV to 1 keV and for angles of incidence between 6 deg and 15 deg. We observed surface charging during the surface ionisation processes for one of the CVD samples due to low electrical conductivity of the material. Measurements on the other CVD diamond sample resulted in ionisation efficiencies of ~2 % for H and up to 12 % for O. Analysis of the angular scattering revealed very narrow and almost circular scattering distributions. Comparison of the results with the data of the CS of the IBEX-Lo sensor shows that CVD diamond has great potential as CS material for future space missions.

High Resolution Transmission Spectroscopy as a Diagnostic for Jovian Exoplanet Atmospheres: Constraints from Theoretical Models

We present high resolution transmission spectra of giant planet atmospheres from a coupled 3-D atmospheric dynamics and transmission spectrum model that includes Doppler shifts which arise from winds and planetary motion. We model jovian planets covering more than two orders of magnitude in incident flux, corresponding to planets with 0.9 to 55 day orbital periods around solar-type stars. The results of our 3-D dynamical models reveal certain aspects of high resolution transmission spectra that are not present in simple 1-D models. We find that the hottest planets experience strong substellar to anti-stellar (SSAS) winds, resulting in transmission spectra with net blue shifts of up to 3 km s−1, whereas less irradiated planets show almost no net Doppler shifts. Compared to 1-D models, peak line strengths are significantly reduced for the hottest atmospheres owing to Doppler broadening from a combination of rotation (which is faster for close-in planets under the assumption of tidal locking) and atmospheric winds. Finally, high resolution transmission spectra may be useful in studying the atmospheres of exoplanets with optically thick clouds since line cores for very strong transitions should remain optically thick to very high altitude. High resolution transmission spectra are an excellent observational test for the validity of 3-D atmospheric dynamics models, because they provide a direct probe of wind structures and heat circulation. Ground-based exoplanet spectroscopy is currently on the verge of being able to verify some of our modeling predictions, most notably the dependence of SSAS winds on insolation. We caution that interpretation of high resolution transmission spectra based on 1-D atmospheric models may be inadequate, as 3-D atmospheric motions can produce a noticeable effect on the absorption signatures.

Intraocular pressure during a very high altitude climb

Reports on intraocular pressure (IOP) changes at high altitudes have provided inconsistent and even conflicting

Detection of degenerative cartilage disease: comparison of high-resolution morphological MR and quantitative T2 mapping at 3.0 Tesla

The aim of the study was to investigate the association of T2 relaxation times of the knee with early degenerative cartilage changes. Furthermore the impact of unloading the knee on T2 values was evaluated.