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

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

Estudio de la distribución espacial de la temperatura en el valle del río Neuquén con la utilización de imágenes satelitales

En el presente trabajo se estudia la distribución espacial de la temperatura de brillo y superficie, con la utilización de imágenes satelitales LANDSAT7 ETM+ y NOAA-AVHRR en el valle cultivado del río Neuquén. El estudio de la distribución espacial de la temperatura en una zona con un terreno algo complejo implica la utilización de una gran densidad de mediciones meteorológicas. Muchas veces no es posible concretar la densidad adecuada de la red agrometeorológica debido a los altos costos de instalación y mantenimiento. Los sensores remotos proporcionan un gran caudal de información en una variedad de resoluciones, con costos considerablemente menores. La determinación de zonas cálidas y frías en el valle permitiría hacer más eficientes los métodos de riego y de lucha contra heladas y proporcionaría herramientas para la mejor planificación productiva de la zona.

Estudio de la distribución espacial de la temperatura en el valle del río Neuquén con la utilización de imágenes satelitales

En el presente trabajo se estudia la distribución espacial de la temperatura de brillo y superficie, con la utilización de imágenes satelitales LANDSAT7 ETM+ y NOAA-AVHRR en el valle cultivado del río Neuquén. El estudio de la distribución espacial de la temperatura en una zona con un terreno algo complejo implica la utilización de una gran densidad de mediciones meteorológicas. Muchas veces no es posible concretar la densidad adecuada de la red agrometeorológica debido a los altos costos de instalación y mantenimiento. Los sensores remotos proporcionan un gran caudal de información en una variedad de resoluciones, con costos considerablemente menores. La determinación de zonas cálidas y frías en el valle permitiría hacer más eficientes los métodos de riego y de lucha contra heladas y proporcionaría herramientas para la mejor planificación productiva de la zona.

Estudio de la distribución espacial de la temperatura en el valle del río Neuquén con la utilización de imágenes satelitales

En el presente trabajo se estudia la distribución espacial de la temperatura de brillo y superficie, con la utilización de imágenes satelitales LANDSAT7 ETM+ y NOAA-AVHRR en el valle cultivado del río Neuquén. El estudio de la distribución espacial de la temperatura en una zona con un terreno algo complejo implica la utilización de una gran densidad de mediciones meteorológicas. Muchas veces no es posible concretar la densidad adecuada de la red agrometeorológica debido a los altos costos de instalación y mantenimiento. Los sensores remotos proporcionan un gran caudal de información en una variedad de resoluciones, con costos considerablemente menores. La determinación de zonas cálidas y frías en el valle permitiría hacer más eficientes los métodos de riego y de lucha contra heladas y proporcionaría herramientas para la mejor planificación productiva de la zona.

Comparisons of observed ice properties and Envisat Advanced Synthetic Aperture Radar (ASAR) data during Nathanial B. Palmer cruise NBP09-01 and Oden cruise OSO08/09

Envisat Advanced Synthetic Aperture Radar (ASAR) Wide Swath Mode (WSM) images are used to derive C-band HH-polarization normalized radar cross sections (NRCS). These are compared with ice-core analysis and visual ship-based observations of snow and ice properties observed according to the Antarctic Sea Ice Processes and Climate (ASPeCt) protocol during two International Polar Year summer cruises (Oden 2008 and Palmer 2009) in West Antarctica. Thick first-year (TFY) and multi-year (MY) ice were the dominant ice types. The NRCS value ranges between -16.3 ± 1.1 and -7.6 ± 1.0 dB for TFY ice, and is -12.6 ± 1.3 dB for MY ice; for TFY ice, NRCS values increase from ~-15 dB to -9 dB from December/January to mid-February. In situ and ASPeCt observations are not, however, detailed enough to interpret the observed NRCS change over time. Co-located Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E) vertically polarized 37 GHz brightness temperatures (TB37V), 7 day and 1 day averages as well as the TB37V difference between ascending and descending AMSR-E overpasses suggest the low NRCS values (-15 dB) are associated with snowmelt being still in progress, while the change towards higher NRCS values (-9dB) is caused by commencement of melt-refreeze cycles after about mid-January.

(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.

(Table 3) Air temperatures, and snow and ice characteristics of fresh and brackish water lakes in the Northwest Territories, Canada

The algorithms designed to estimate snow water equivalent (SWE) using passive microwave measurements falter in lake-rich high-latitude environments due to the emission properties of ice covered lakes on low frequency measurements. Microwave emission models have been used to simulate brightness temperatures (Tbs) for snowpack characteristics in terrestrial environments but cannot be applied to snow on lakes because of the differing subsurface emissivities and scattering matrices present in ice. This paper examines the performance of a modified version of the Helsinki University of Technology (HUT) snow emission model that incorporates microwave emission from lake ice and sub-ice water. Inputs to the HUT model include measurements collected over brackish and freshwater lakes north of Inuvik, Northwest Territories, Canada in April 2008, consisting of snowpack (depth, density, and snow water equivalent) and lake ice (thickness and ice type). Coincident airborne radiometer measurements at a resolution of 80x100 m were used as ground-truth to evaluate the simulations. The results indicate that subsurface media are simulated best when utilizing a modeled effective grain size and a 1 mm RMS surface roughness at the ice/water interface compared to using measured grain size and a flat Fresnel reflective surface as input. Simulations at 37 GHz (vertical polarization) produce the best results compared to airborne Tbs, with a Root Mean Square Error (RMSE) of 6.2 K and 7.9 K, as well as Mean Bias Errors (MBEs) of -8.4 K and -8.8 K for brackish and freshwater sites respectively. Freshwater simulations at 6.9 and 19 GHz H exhibited low RMSE (10.53 and 6.15 K respectively) and MBE (-5.37 and 8.36 K respectively) but did not accurately simulate Tb variability (R= -0.15 and 0.01 respectively). Over brackish water, 6.9 GHz simulations had poor agreement with airborne Tbs, while 19 GHz V exhibited a low RMSE (6.15 K), MBE (-4.52 K) and improved relative agreement to airborne measurements (R = 0.47). Salinity considerations reduced 6.9 GHz errors substantially, with a drop in RMSE from 51.48 K and 57.18 K for H and V polarizations respectively, to 26.2 K and 31.6 K, although Tb variability was not well simulated. With best results at 37 GHz, HUT simulations exhibit the potential to track Tb evolution, and therefore SWE through the winter season.

High resolution color reflectance and hematite and goethite record of ODP Hole 184-1143A in the South China Sea

Precipitation has a larger variability than temperature in tropical monsoon regions, thus it is an important climate variable. However, reconstructions of long-term rainfall histories are scarce because of the lack of reliable proxies. Here we document that iron oxide minerals, specifically the ratio of hematite to goethite (Hm/Gt), is a reasonable precipitation proxy. Using diffuse reflectance spectrophotometry, we measured samples from Ocean Drilling Program (ODP) 1143 drilling site (9°21.72'N, 113°17.11'E, 2777 m water depth) for hematite and goethite, whose formation processes are favored by opposing climate conditions. In order to determine the content of hematite and goethite we produced a set of calibration samples by removing the iron oxides to generate the natural matrix to which hematite and goethite in known percentages were added. From these calibration samples we developed a transfer function for determining hematite and goethite concentration from a sample's spectral reflectance. Applying this method to ODP 1143 sediments (top 34 m of a 510 m core with sampling interval of 10 cm) we were able to reconstruct a continuous precipitation history for SE Asia of the past 600 kyr using the Hm/Gt ratio as a proxy of the precipitation variability of Asian monsoon. The reliability of this Hm/Gt proxy is corroborated by its consistency with the stalagmite delta18O data from South China. Comparing long-term Hm/Gt records with the surface temperature gradient of equatorial Pacific Ocean, we found that monsoon precipitation and El Niño are correlated for the last 600 kyr. The development of El Niño-like conditions decreased SE Asia precipitation, whereas precipitation increases in response to La Niña intensification

(Table 1) Mean temperatures in january and july, and total annual snowfall at Great Bear Lake and Great Slave Lake, Canada

The sensitivity of brightness temperature (T(B)) at 6.9, 10.7, and 18.7 GHz from Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E) observations is investigated over five winter seasons (2002-2007) on Great Bear Lake and Great Slave Lake, Northwest Territories, Canada. The T(B) measurements are compared to ice thicknesses obtained with a previously validated thermodynamic lake ice model. Lake ice thickness is found to explain much of the increase of T(B) at 10.7 and 18.7 GHz. T(B) acquired at 18.7 GHz (V-pol) and 10.7 GHz (H-pol) shows the strongest relation with simulated lake ice thickness over the period of study (R**2 > 0.90). A comparison of the seasonal evolution of T(B) for a cold winter (2003-2004) and a warm winter (2005-2006) reveals that the relationship between T(B) and ice growth is stronger in the cold winter (2003-2004). Overall, this letter shows the high sensitivity of T(B) to ice growth and, thus, the potential of AMSR-E mid-frequency channels to estimate ice thickness on large northern lakes.

Temperature affects the morphology and calcification of Emiliania huxleyi strains

The global warming debate has sparked an unprecedented interest in temperature effects on coccolithophores. The calcification response to temperature changes reported in the literature, however, is ambiguous. The two main sources of this ambiguity are putatively differences in experimental setup and strain-specificity. In this study we therefore compare three strains isolated in the North Pacific under identical experimental conditions. Three strains of Emiliania huxleyi type A were grown under non-limiting nutrient and light conditions, at 10, 15, 20 and 25 ºC. All three strains displayed similar growth rate versus temperature relationships, with an optimum at 20-25 ºC. Elemental production (particulate inorganic carbon (PIC), particulate organic carbon (POC), total particulate nitrogen (TPN)), coccolith mass, coccolith size, and width of the tube elements cycle were positively correlated with temperature over the sub-optimum to optimum temperature range. The correlation between PIC production and coccolith mass/size supports the notion that coccolith mass can be used as a proxy for PIC production in sediment samples. Increasing PIC production was significantly positively correlated with the percentage of incomplete coccoliths in one strain only. Generally, coccoliths were heavier when PIC production was higher. This shows that incompleteness of coccoliths is not due to time shortage at high PIC production. Sub-optimal growth temperatures lead to an increase in the percentage of malformed coccoliths in a strain-specific fashion. Since in total only six strains have been tested thus far, it is presently difficult to say whether sub-optimal temperature is an important factor causing malformations in the field. The most important parameter in biogeochemical terms, the PIC:POC, shows a minimum at optimum growth temperature in all investigated strains. This clarifies the ambiguous picture featuring in the literature, i.e. discrepancies between PIC:POC-temperature relationships reported in different studies using different strains and different experimental setups. In summary, global warming might cause a decline in coccolithophore's PIC contribution to the rain ratio, as well as improved fitness in some genotypes due to less coccolith malformations.

OSIRIS observations of meter-sized exposures of H2O ice at the surface of 67P/Churyumov-Gerasimenko and interpretation using laboratory experiments

Since OSIRIS started acquiring high-resolution observations of the surface of the nucleus of comet 67P/Churyumov-Gerasimenko, over one hundred meter-sized bright spots have been identified in numerous types of geomorphologic regions, but mostly located in areas receiving low insolation. The bright spots are either clustered, in debris fields close to decameter-high cliffs, or isolated without structural relation to the surrounding terrain. They can be up to ten times brighter than the average surface of the comet at visible wavelengths and display a significantly bluer spectrum. They do not exhibit significant changes over a period of a few weeks. All these observations are consistent with exposure of water ice at the surface of boulders produced by dislocation of the weakly consolidated layers that cover large areas of the nucleus. Laboratory experiments show that under simulated comet surface conditions, analog samples acquire a vertical stratification with an uppermost porous mantle of refractory dust overlaying a layer of hard ice formed by recondensation or sintering under the insulating dust mantle. The evolution of the visible spectrophotometric properties of samples during sublimation is consistent with the contrasts of brightness and color seen at the surface of the nucleus. Clustered bright spots are formed by the collapse of overhangs that is triggered by mass wasting of deeper layers. Isolated spots might be the result of the emission of boulders at low velocity that are redepositioned in other regions.

The VLT/NaCo large program to probe the occurrence of exoplanets and brown dwarfs at wide orbits II. Survey description, results, and performances

Context. Young, nearby stars are ideal targets for direct imaging searches for giant planets and brown dwarf companions. After the first-imaged planet discoveries, vast efforts have been devoted to the statistical analysis of the occurence and orbital distributions of giant planets and brown dwarf companions at wide (>= 5-6 AU) orbits. Aims. In anticipation of the VLT/SPHERE planet-imager, guaranteed-time programs, we have conducted a preparatory survey of 86 stars between 2009 and 2013 to identify new faint comoving companions to ultimately analyze the occurence of giant planets and brown dwarf companions at wide (10-2000 AU) orbits around young, solar-type stars. Methods. We used NaCo at VLT to explore the occurrence rate of giant planets and brown dwarfs between typically 0.1 and 8 ''. Diffraction-limited observations in H-band combined with angular differential imaging enabled us to reach primary star-companion brightness ratios as small as 10(-6) at 1.5 ''. Repeated observations at several epochs enabled us to discriminate comoving companions from background objects. Results. During our survey, twelve systems were resolved as new binaries, including the discovery of a new white dwarf companion to the star HD8049. Around 34 stars, at least one companion candidate was detected in the observed field of view. More than 400 faint sources were detected; 90% of them were in four crowded fields. With the exception of HD8049 B, we did not identify any new comoving companions. The survey also led to spatially resolved images of the thin debris disk around HD61005 that have been published earlier. Finally, considering the survey detection limits, we derive a preliminary upper limit on the frequency of giant planets for the semi-major axes of [10, 2000] AU: typically less than 15% between 100 and 500 AU and less than 10% between 50 and 500 AU for exoplanets that are more massive than 5 M-Jup and 10 M-Jup respectively, if we consider a uniform input distribution and a confidence level of 95%. Conclusions. The results from this survey agree with earlier programs emphasizing that massive, gas giant companions on wide orbits around solar-type stars are rare. These results will be part of a broader analysis of a total of similar to 210 young, solar-type stars to bring further statistical constraints for theoretical models of planetary formation and evolution.

Land Surface Temperature continuous time series from SSM/I passive microwave radiometers, over northern latitudes (Latitudes > 45°N), during summer snow free period (2000-2011), with links to NetCDF files

Remote sensing instruments are key players to map land surface temperature (LST) at large temporal and spatial scales. In this paper, we present how we combine passive microwave and thermal infrared data to estimate LST during summer snow-free periods over northern high latitudes. The methodology is based on the SSM/I-SSMIS 37 GHz measurements at both vertical and horizontal polarizations on a 25 km × 25 km grid size. LST is retrieved from brightness temperatures introducing an empirical linear relationship between emissivities at both polarizations as described in Royer and Poirier (2010). This relationship is calibrated at pixel scale, using cloud-free independent LST data from MODIS instruments. The SSM/I-SSMIS and MODIS data are synchronized by fitting a diurnal cycle model built on skin temperature reanalysis provided by the European Centre for Medium-Range Weather Forecasts (ECMWF). The resulting temperature dataset is provided at 25 km scale and at an hourly time step during the ten-year analysis period (2000-2011). This new product was locally evaluated at five experimental sites of the EU-PAGE21 project against air temperature measurements and meteorological model reanalysis, and compared to the MODIS LST product at both local and circumpolar scale. The results giving a mean RMSE of the order of 2.2 K demonstrate the usefulness of the microwave product, which is unaffected by clouds as opposed to thermal infrared products and offers a better resolution compared to model reanalysis.

Efficient Model-based 3D Tracking by Using Direct Image Registration

This thesis deals with the problem of efficiently tracking 3D objects in sequences of images. We tackle the efficient 3D tracking problem by using direct image registration. This problem is posed as an iterative optimization procedure that minimizes a brightness error norm. We review the most popular iterative methods for image registration in the literature, turning our attention to those algorithms that use efficient optimization techniques. Two forms of efficient registration algorithms are investigated. The first type comprises the additive registration algorithms: these algorithms incrementally compute the motion parameters by linearly approximating the brightness error function. We centre our attention on Hager and Belhumeur’s factorization-based algorithm for image registration. We propose a fundamental requirement that factorization-based algorithms must satisfy to guarantee good convergence, and introduce a systematic procedure that automatically computes the factorization. Finally, we also bring out two warp functions to register rigid and nonrigid 3D targets that satisfy the requirement. The second type comprises the compositional registration algorithms, where the brightness function error is written by using function composition. We study the current approaches to compositional image alignment, and we emphasize the importance of the Inverse Compositional method, which is known to be the most efficient image registration algorithm. We introduce a new algorithm, the Efficient Forward Compositional image registration: this algorithm avoids the necessity of inverting the warping function, and provides a new interpretation of the working mechanisms of the inverse compositional alignment. By using this information, we propose two fundamental requirements that guarantee the convergence of compositional image registration methods. Finally, we support our claims by using extensive experimental testing with synthetic and real-world data. We propose a distinction between image registration and tracking when using efficient algorithms. We show that, depending whether the fundamental requirements are hold, some efficient algorithms are eligible for image registration but not for tracking.

Tiled Projector Displays Correction for Dark Scenes in Railway Simulators

Tiled projector displays are a common choice for training simulators, where a high resolution output image is required. They are cheap for the resolution that they can reach and can be configured in many different ways. Nevertheless, such kinds of displays require geometric and color correction so that the composite image looks seamless. Display correction is an even bigger challenge when the projected images include dark scenes combined with brighter scenes. This is usually a problem for railway simulators when the train is positioned inside a tunnel and the black offset effect becomes noticeable. In this paper, a method for fast photometric and geometric correction of tiled display systems where dark and bright scenes are combined is presented. The image correction is carried out in two steps. First, geometric alignment and overlapping areas attenuation for brighter scenes is applied. Second, in the event of being inside a tunnel, the brightness of the scene is increased in certain areas using light sources in order to create the impression of darkness but minimizing the effect of the black offset