ALBERT

Description: Postdisaster search and rescue is an important application of ultrawideband (UWB) radar systems, which mainly detect trapped victims by their respiratory-motion response. The development of a respiratory-motion detection (RMD) algorithm that can eliminate nonstationary clutter and noise is a challenging task for the application. A new algorithm is proposed to deal with the task in this letter. It uses the multichannel singular spectrum analysis (MSSA) technique to reconstruct the respiratory-motion response detected by a UWB radar. During the reconstruction, the periodicity and range interrelation characteristics of the response are exploited to adaptively identify signal subspaces. The performance of the algorithm is verified both by simulated and real data. The results show its improved performance over the reference algorithms, e.g., a singular-value-decomposition-based algorithm. The adaptive-MSSA-based RMD algorithm has great promise not only in practical use but also for future research of UWB-radar-based human being remote sensing.

Description: Land–sea contamination observed in Soil Moisture and Ocean Salinity (SMOS) brightness temperature images is found to have two main contributions: the floor error inherent of image reconstruction and a multiplicative error either in the antenna temperature or in the visibility samples measured by the correlator. The origin of this last one is traced down to SMOS calibration parameters to yield a simple correction scheme, which is validated against several geophysical scenarios. Autoconsistency rules in interferometric synthesis together with redundant and complementary calibration procedures provide a robust SMOS calibration scheme.

Description: This letter proposes a novel algorithm, which is based on the generalized method of moments (GMM), for the estimation and correction of phase errors induced in synthetic aperture radar (SAR) imagery. The GMM algorithm is used to replace the original phase-estimation kernel in the basic structure of the phase-gradient-autofocus algorithm. Since this novel algorithm does not require the observed signal to be a certain distribution model, it is able to estimate arbitrary phase errors. The GMM algorithm has the ability of estimating range-dependent phase errors, which makes it an efficient estimator. As a result, higher accuracy of the estimated phase errors and a better focused image can be achieved. Excellent results have been obtained in autofocusing and imaging experiments on real SAR data.

Description: Strong clutter reflections of terrain and marine surfaces obscure the contrast between the target-of-interest and clutter (terrain and marine surface reflections) in synthetic aperture radar (SAR) images and consequently hinder the efficiency of image interpretation and analysis. To overcome this problem, this letter proposes an efficient clutter suppression method in SAR images, which is named shedding irrelevant patterns (SIP). The essence is to construct a regression function that can suppress clutter and preserve the target patterns concurrently. We assume that the clutter is irrelevant to the target-of-interest and distinguishable in patterns in terms of image-pixel distribution and intensity (spatial information). Experimental results show the efficiency of the proposed method in both clutter suppression and target pattern preservation.

Description: A concern in hyperspectral image classification is the high number of required training samples. When traditional classifiers are applied, feature reduction (FR) techniques are the most common approaches to deal with this problem. Subspace-based classifiers, which are developed based on high-dimensional space characteristics, are another way to handle the high dimension of hyperspectral images. In this letter, a novel subspace-based classification approach is proposed and compared with basic and improved subspace-based classifiers. The proposed classifier is also compared with traditional classifiers that are accompanied by an FR technique and the well-known support vector machine classifier. Experimental results prove the efficiency of the proposed method, especially when a limited number of training samples are available. Furthermore, the proposed method has a very high level of automation and simplicity, as it has no parameters to be set.

Description: To improve the spatial density of measurement points of persistent-scatterer interferometry, distributed scatterer (DS) should be considered and processed. An important procedure in DS interferometry is the phase triangulation (PT). This letter introduces two modified PT algorithms (i.e., equal-weighted PT and coherence-weighted PT) and analyzes the mathematical relations between different published PT methods (i.e., the maximum-likelihood phase estimator, least squares estimator, and eigendecomposition-based phase estimators). The analysis shows that the above five PT methods share very similar mathematical forms with different weight values in the estimation procedure.

Description: Sparse representation-based classifier and its variants have been widely adopted for hyperspectral image (HSI) classification recently. However, sparse representation is unstable so that similar features might obtain significantly different sparse codes. Despite the instability, we find that the sparse codes follow a class-dependent distribution under the structured dictionary consisting of training samples from all classes. Based on this observation, a novel discriminative feature, sparse code histogram (SCH), is developed for HSI classification. By counting the SCH of each sample from the sparse codes of its spatial neighbors, we can statistically obtain the distribution pattern of sparse codes of the class to which the sample belongs, and then treat the SCH as a new feature for classification. To reduce the possible outliers among the neighbors, a shape-adaptive neighborhood extractor is also employed to enhance the stability of the histogram feature. Experimental results demonstrate that SCH enjoys a strong discriminative power, which can achieve notably better performance than several state-of-the-art methods for HSI classification with limited training samples.

Description: Several detection statistics have been proposed for detecting fine ground disturbances between two synthetic aperture radar (SAR) images, such as vehicle tracks. The standard method involves estimating a local correlation coefficient between images. Other methods have been proposed using various statistical hypothesis tests. One of these alternative methods is a generalized likelihood ratio test (GLRT), which compares a full-correlation image model to a no-correlation image model. In this letter, we expand the GLRT to polarimetric SAR data and derive the appropriate GLRT detection statistics. Additionally, we explore relaxing the equal variance/equal polarimetric covariance assumptions used in previous results and find improved performance on macroscopic scene changes.

Description: For downward-looking linear array 3-D synthetic aperture radar, the resolution in cross-track direction is much lower than the ones in range and azimuth. Hence, superresolution reconstruction algorithms are desired. Since the cross-track signal to be reconstructed is sparse in the object domain, compressive sensing algorithm has been used. However, the imaging processing on the 3-D scene brings large computational loads, which renders challenges in both data acquisition and processing. To cover this shortage, truncated singular value decomposition is utilized to reconstruct a reduced-redundancy spatial measurement matrix. The proposed algorithm provides advantages in terms of computational time while maintaining the quality of the scene reconstructions. Moreover, our results on uniform linear array are generally applicable to sparse nonuniform linear array. Superresolution properties and reconstruction accuracies are demonstrated using simulations under the noise and clutter scenarios.

Description: This letter proposes a signal processing method of passive bistatic radar (PBR) exploiting an uncooperative radar as an illuminator. Compared with other opportunity illuminators, the transmitting signal of a radar usually has a better ambiguity function, which leads to a higher range resolution. Two channels are needed in PBR system. The reference channel is used to estimate radar signal parameters and reconstruct directly propagated signal. The surveillance channel is used to receive scattered wave. An array antenna and a simultaneous multibeam algorithm are necessary in the surveillance channel due to the flexible beam scanning of the uncooperative radar. The procedure of the proposed method is explained in detail, which is then followed by a field experiment. Preliminary results from the field experiment show that the proposed method can be applied to target angle and bistatic range measurement successfully.

Description: In order to achieve 3-D imaging with an airborne down-looking linear-array synthetic aperture radar (LASAR), a uniform virtual antenna array may be obtained by aperture synthesis of the cross-track sparse multiple-input–multiple-output array. However, the actual 3-D imaging quality is unavoidably degraded by errors in the virtual element position. In this letter, we investigate the effects of these errors on the forms and the degrees of image quality degradation by decomposing the error-related stochastic processes via an orthogonal transform based on discrete Legendre polynomials. It should be noted that these analyses are helpful for designing a LASAR system and providing a reference for specifying the requisite precision of measurement devices and calibration methods. Finally, we briefly consider the use of calibration methods to eliminate the effects of errors.

Description: In problems where labeled data are scarce, semisupervised learning (SSL) techniques are an attractive framework that can exploit both labeled and unlabeled data. These approaches typically rely on a smoothness assumption such that examples that are similar in input space should also be similar in label space. In many domains, such as remotely sensed hyperspectral image (HSI) classification, the data violate this assumption. In response, we propose a general method by which a neighborhood graph used in SSL is learned using binary classifiers that are trained to predict whether a pair of pixels shares the same label. Working within the framework of semisupervised neural networks (SSNNs), we show that our approach improves on the performance of the SSNN on two HSI data sets.

Description: In this letter, new models for the spatial correlation of sea clutter texture and intensity are proposed as improved versions of current power law models or exponential decay model. The models for texture have three unknown parameters, and thus can be called triparametric models. The structure of the models is a weighted sum of two components, which can describe the decaying process of the correlation coefficient with spatial lags, as well as the periodic behavior due to the existence of transient coherent structures in sea clutter. Unknown parameters are optimized by the nonlinear least square fit method. Models for sea clutter intensity can be obtained through a linear transform for uncorrelated speckle based on the compound-Gaussian representation of sea clutter. The proposed models are validated and compared with current models using S- and C-band measured sea clutter data. Analysis results indicate the effectiveness of the proposed models in that they can describe the behavior of spatial correlation coefficients with higher accuracy.

Description: Terrestrial laser scanning (TLS) has become a popular tool for acquiring source data points which can be used to construct digital elevation models (DEMs) for a wide number of applications. A TLS point cloud often has a very fine spatial resolution, which can represent well the spatial variation of a terrain surface. However, the uncertainty in DEMs created from this relatively new type of source data is not well understood, which forms the focus of this letter. TLS survey data representing four terrain surfaces of different characteristics were used to explore the effects of surface complexity and typical TLS data density (in terms of data point spacing) on DEM accuracy. The spatial variation in TLS data can be decomposed into parts corresponding to the signal of spatial variation (of terrain surfaces) and noise due to measurement error. We found a linear relation between the DEM error and the typical TLS data spacings considered (30–100 mm) which arises as a function of the interpolation error, and a constant contribution from the propagated data noise. This letter quantifies these components for each of the four surfaces considered and shows that, for the interpolation method considered here, higher density sampling would not be beneficial.

Description: A novel way to estimate the live fuel moisture content (LFMC) was explored from the ratio of canopy water content (CWC) and foliage dry biomass (FDB). The CWC was estimated using the PROSAIL (PROSPECT + SAIL) radiative transfer model from the Landsat 8 product. A weak constraint 4-D variational data assimilation method was employed to assimilate the temporally estimated leaf area index into a soil-water-atmosphere-plant (SWAP) model for optimizing the model control variables. Then, the SWAP model was reinitialized with this optimum set of control variables, and better prediction of FDB was obtained. Results showed that a high accuracy level was achieved for the estimated CWC ( $R^{2}=0.91$ , $mbox{RMSE}=84.74 mbox{g/m}^2$ ) and FDB ( $R^2=0.88$ , $mbox{RMSE}=48.54 mbox{g/m} ^2$ ) when compared with in situ measured values. However, the accuracy level of estimated LFMC was poor ( $R^2=0.59$ , $mbox{RMSE} =30.85%$ ) . Further analyses find that the estimated LFMC is reliable for low LFMC but challenged for high LFMC, which indicates that the presented method still makes sense to the assessment of wildfire risk since the wildfire generally occurs when the vegetation is in low LFMC condition.

Description: In this letter, we present the use of experimental human micro-Doppler signature data gathered by a multistatic radar system to discriminate between unarmed and potentially armed personnel walking along different trajectories. Different ways of extracting suitable features from the spectrograms of the micro-Doppler signatures are discussed, particularly empirical features such as Doppler bandwidth, periodicity, and others, and features extracted from singular value decomposition (SVD) vectors. High classification accuracy of armed versus unarmed personnel (between 90% and 97% depending on the walking trajectory of the people) can be achieved with a single SVD-based feature, in comparison with using four empirical features. The impact on classification performance of different aspect angles and the benefit of combining multistatic information is also evaluated in this letter.

Description: Automatic urban area detection in remote sensing images is an important application in the field of earth observation. Most of the existing methods employ feature classifiers and thereby contain a data training process. Moreover, some methods cannot detect urban areas in complex scenes accurately. This letter proposes an automatic urban area detection method that uses multiple features that have different resolutions. First, a downsampled low-resolution image is used to segment the candidate area. After the corner points of the urban area are extracted, a weighted Gaussian voting matrix technique is employed to integrate the corner points into the candidate area. Then, the edge features and homogeneous region are extracted by using the original high-resolution image. Using these results as the input, the processes of guided filtering and contrast enhancement can finally detect accurately the urban areas. This method combines multiple features, such as corner, edge, and regional characteristics, to detect the urban areas. The experimental results show that the proposed method has better detection accuracy for urban areas than the existing algorithms.

Description: In marine sciences, time series are often nonlinear and nonstationary. Adequate and specific methods are needed to analyze such series. In this letter, an application of the empirical mode decomposition method (EMD) associated to the Hilbert spectral analysis (HSA) is presented. Furthermore, EMD-based time-dependent intrinsic correlation (TDIC) analysis is applied to consider the correlation between two nonstationary time series. Four temperature time series obtained from automatic measurements in nearshore waters of the Réunion island are considered, recorded every 10 min from July 2011 to January 2012. The application of the EMD on these series and the estimation of their power spectra using the HSA are illustrated. The authors identify low-frequency tidal waves and display the pattern of correlations at different scales and different locations. By TDIC analysis, it was concluded that the high-frequency modes have small correlation, whereas the trends are perfectly correlated.

Description: Detecting vehicles in aerial images provides important information for traffic management and urban planning. Detecting the cars in the images is challenging due to the relatively small size of the target objects and the complex background in man-made areas. It is particularly challenging if the goal is near-real-time detection, i.e., within few seconds, on large images without any additional information, e.g., road database and accurate target size. We present a method that can detect the vehicles on a 21-MPixel original frame image without accurate scale information within seconds on a laptop single threaded. In addition to the bounding box of the vehicles, we extract also orientation and type (car/truck) information. First, we apply a fast binary detector using integral channel features in a soft-cascade structure. In the next step, we apply a multiclass classifier on the output of the binary detector, which gives the orientation and type of the vehicles. We evaluate our method on a challenging data set of original aerial images over Munich and a data set captured from an unmanned aerial vehicle (UAV).

Description: We compare five slope correction methods developed by Walter et al. , Montes et al. , Schleppi et al. , España et al. , and Gonsamo et al. (referred to as WAL, MON, SCH, ESP, and GON, respectively) using artificial fisheye pictures simulated by graphics software and a lookup table (LUT) retrieval method. The LUT is built by simulating the directional gap fraction as a function of leaf area index (LAI) and average leaf inclination angle (ALIA) using the Poisson law. LAI and ALIA estimates correspond to the case of the LUT that provides the lowest root-mean-square error between the observed gap fractions after slope correction and the simulated ones. Three LAI values (1.5, 3.5, and 5.5), four ALIA values (26.8°, 45°, 57.5°, and 63.2°), and three slope angles (0°, 20°, and 50°) constituted 36 samples of random scenes. ESP is recommended because its results are accurate and independent on the leaf angle distribution (LAD), while GON only performs well for spherical LAD. The three other methods present less good performances with underestimation or overestimation of LAI and/or ALIA depending on the LAD, and the recommended order for them is MON, SCH, and WAL.

Description: In this letter, an improved phase correlation (PC) method based on 2-D plane fitting and the maximum kernel density estimator (MKDE) is proposed, which combines the idea of Stone's method and robust estimator MKDE. The proposed PC method first utilizes a vector filter to minimize the noise errors of the phase angle matrix and then unwraps the filtered phase angle matrix by the use of the minimum cost network flow unwrapping algorithm. Afterward, the unwrapped phase angle matrix is robustly fitted via MKDE, and the slope coefficients of the 2-D plane indicate the subpixel shifts between images. The experiments revealed that the improved method can effectively avoid the impact of outliers on the phase angle matrix during the plane fitting and is robust to aliasing and noise. The matching accuracy can reach 1/50th of a pixel using simulated data. The real image sequence tracking experiment was also undertaken to demonstrate the effectiveness of the proposed PC method with a registration accuracy of root-mean-square error better than 0.1 pixels.

Description: Accurately mapping forest carbon density by combining sample plots and remotely sensed images has become popular because this method provides spatially explicit estimates. However, mixed pixels often impede the improvement of the estimation. In this letter, regression modeling and spectral unmixing analysis were integrated to improve the estimation of forest carbon density for the You County of Hunan, China, using Landsat Thematic Mapper images. Linear spectral unmixing with and without a constraint (LSUWC and LSUWOC) and nonlinear spectral unmixing (NSU) were compared to derive the fractions of five endmembers, particularly forests. Stepwise regression, logistic regression, and polynomial regression (PR) with and without the forest fraction used as an independent variable and the product of the forest fraction image and the map from the best model without the forest fraction were compared. The models were developed using 56 sample plots, and their results were validated using 26 test plots. The decomposition of mixed pixels was assessed using higher spatial resolution SPOT images and a corresponding land cover map. The results showed that 1) LSUWC more accurately estimated the endmember fractions than LSUWOC and NSU, 2) PR had the greatest estimation accuracy of forest carbon, and 3) combining regression modeling and spectral unmixing increased the estimation accuracy by 31%–39%, and introducing the forest fraction into the regressions performed better than the product of forest fraction image and the results from PR without the fraction. This implied that the integrations provided great potential in reducing the impacts of mixed pixels in mapping forest carbon.

Description: This letter proposes a multiresolution technique to address the high computational cost in remote sensing image registration. The scale-invariant feature transform is applied to detect keypoints and descriptors, and then, global information combined with descriptors is utilized to establish keypoint mappings. Keypoints are first classified according to their octaves. Then, in the lowest resolution, the keypoints of the largest octave are mapped with descriptors and the global information, giving an initial affine transformation $T_0$ . In the next octave, the keypoints of the second largest octave are mapped by employing $T_0$ to narrow the space of matching keypoints. By this means, the process of establishing keypoint correspondences is conducted from one resolution (octave) to the next as the obtained transformation gets finer until we get to the highest resolution. Due to the high computational expense of computing global information, the proposed technique is important for aligning large-size remote sensing imagery. Experimental results show that the proposed method can achieve a comparable registration accuracy but with a less computational cost than directly building keypoint mappings on images of large size.

Description: In this letter, we present an efficient parallel implementation of composite kernels in support vector machines (SVMs) for hyperspectral image (HSI) classification. Our implementation makes effective use of commodity graphics processing units (GPUs). Specifically, we port the calculation of composite kernels to GPUs, perform intensive computations based on NVidia's compute unified device architecture, and execute the rest of the operations related with control and small data calculations in the CPU. Our experimental results, conducted using real hyperspectral data sets and NVidia GPU platforms, indicate significant improvements in terms of computational effectiveness, achieving near-real-time performance of spatial–spectral HSI classification for the first time in the literature.

Description: Compared with airborne laser scanning, terrestrial laser scanning (TLS) offers ground-based point cloud data of trees and provides greater potential to accurately estimate tree and stand parameters. However, there is a lack of effective methods to accurately identify locations of individual trees from TLS point cloud data. It is also unknown whether the estimation accuracy of the parameters, including tree height (H), diameter at breast height (DBH), and so on, using TLS can meet the requirement of forest management and planning. In this letter, a novel method to effectively process point cloud data and further determine the locations of individual trees in a stand based on the central coordinates of point cloud data on a defined grid according to the largest DBH was developed. Moreover, a point-cloud-data-based convex hull algorithm and the cylinder method were, respectively, used to estimate DBH and H of individual trees. This study was conducted in a pure Chinese fir plantation of 45 trees located in Huang-Feng-Qiao forest farm, You County of Hunan, China. The comparison of the estimated and observed values showed that the obtained tree locations had errors of less than 20 cm, and the relative root mean square errors for the estimates of both DBH and H were less than 5%. This implies that TLS is very promising for the retrieval of tree and stand parameters in forest stands. For the applications of these methods to mixed forests with a structure of multilayer canopies, further examination is needed.

Description: In this letter, a novel algorithm for attitude measurement based on a 3-D electromagnetic model (3-D em-model) is proposed. The 3-D em-model is established offline based on the geometric structure of the target, and it can be used to predict the scattering features at different target attitudes. In order to measure the attitude of the air target, we design a bistatic step frequency radar system. The directions of the two radars' lines of sight (LOSs) relative to the target are acquired by matching the high-resolution range profiles (HRRPs) from the target echoes to the HRRPs generated from the 3-D em-model. Since the directions of two radars' LOSs relative to the Earth are already known, the absolute attitude of the target can be acquired. The innovative contributions of this letter are as follows: 1) A comprehensive theoretical analysis of air target attitude measurement based on its own 3-D em-model is proposed; 2) the method can be applied to different kinds of air targets such as aircraft, satellite, missile, etc.; 3) the proposed attitude measurement method does not require target motion model in advance; and 4) the proposed algorithm can be applied to any kind of step frequency waveforms. Experiments using both data predicted by a high-frequency electromagnetic code and data measured in the chamber verify the validity of the method.

Description: The Soil Moisture and Ocean Salinity (SMOS) and Soil Moisture Active Passive (SMAP) missions provide Level-1 brightness temperature (Tb) observations that are used for global soil moisture estimation. However, the nature of these Tb data differs: the SMOS Tb observations contain atmospheric and select reflected extraterrestrial (“Sky”) radiation, whereas the SMAP Tb data are corrected for these contributions, using auxiliary near-surface information. Furthermore, the SMOS Tb observations are multiangular, whereas the SMAP Tb is measured at 40° incidence angle only. This letter discusses how SMOS Tb, SMAP Tb, and radiative transfer modeling components can be aligned in order to enable a seamless exchange of SMOS and SMAP Tb data in soil moisture retrieval and assimilation systems. The aggregated contribution of the atmospheric and reflected Sky radiation is, on average, about 1 K for horizontally polarized Tb and 0.5 K for vertically polarized Tb at 40° incidence angle, but local and short-term values regularly exceed 5 K.

Description: Numerical simulations of the Cordilleran ice sheet through the last glacial cycle The Cryosphere Discussions, 9, 4147-4203, 2015 Author(s): J. Seguinot, I. Rogozhina, A. P. Stroeven, M. Margold, and J. Kleman Despite more than a century of geological observations, the Cordilleran ice sheet of North America remains poorly understood in terms of its former extent, volume and dynamics. Although geomorphological evidence is abundant, its complexity is such that whole ice-sheet reconstructions of advance and retreat patterns are lacking. Here we use a numerical ice sheet model calibrated against field-based evidence to attempt a quantitative reconstruction of the Cordilleran ice sheet history through the last glacial cycle. A series of simulations is driven by time-dependent temperature offsets from six proxy records located around the globe. Although this approach reveals large variations in model response to evolving climate forcing, all simulations produce two major glaciations during marine oxygen isotope stages 4 (61.9–56.5 ka) and 2 (23.2–16.8 ka). The timing of glaciation is better reproduced using temperature reconstructions from Greenland and Antarctic ice cores than from regional oceanic sediment cores. During most of the last glacial cycle, the modelled ice cover is discontinuous and restricted to high mountain areas. However, widespread precipitation over the Skeena Mountains favours the persistence of a central ice dome throughout the glacial cycle. It acts as a nucleation centre before the Last Glacial Maximum and hosts the last remains of Cordilleran ice until the middle Holocene (6.6–6.2 ka).

Description: Utilisation of CryoSat-2 SAR altimeter in operational ice charting The Cryosphere Discussions, 9, 4117-4145, 2015 Author(s): E. Rinne and M. Similä We present methods to utilise Cryosat-2 (CS-2) Synthetic Aperture (SAR) mode data in operational ice charting. We compare CS-2 data qualitatively to Synthetic Aperture Radar (SAR) mosaics over Barents and Kara seas. Furthermore, we compare the CS-2 to archived operational ice charts. We present distributions of four CS-2 waveform parameters for different ice types as presented in the ice charts. We go on to present an automatic classification method for CS-2 data which, after training with operational ice charts, is capable of determining open water from ice with a hit rate of 〉 90 %. The training data is dynamically updated every five days using the most recent 15 days CS-2 data and operative ice charts. This helps the adaption of the classifier to the evolving ice/snow conditions throughout winter. The classifier is also capable of detecting three different ice classes (thin and thick first year ice as well as old ice) with success rates good enough for the output to be usable to support operational ice charting. Finally, we present a near real time CS-2 product just plotting the waveform characteristics and conclude that even such a simple product is usable for some of the needs of ice charting.

Description: Long-term coastal-polynya dynamics in the Southern Weddell Sea from MODIS thermal-infrared imagery The Cryosphere Discussions, 9, 3959-3993, 2015 Author(s): S. Paul, S. Willmes, and G. Heinemann Based upon high-resolution thermal-infrared Moderate-Resolution Imaging Spectroradiometer (MODIS) satellite imagery in combination with ERA-Interim atmospheric reanalysis data, we derived long-term polynya parameters such as polynya area, thin-ice thickness distribution and ice-production rates from daily cloud-cover corrected thin-ice thickness composites. Our study is based on a thirteen year investigation period (2002–2014) for the austral winter (1 April to 30 September) in the Antarctic Southern Weddell Sea. The focus lies on coastal polynyas which are important hot spots for new-ice formation, bottom-water formation and heat/moisture release into the atmosphere. MODIS has the capability to resolve even very narrow coastal polynyas. Its major disadvantage is the sensor limitation due to cloud cover. We make use of a newly developed and adapted spatial feature reconstruction scheme to account for cloud-covered areas. We find the sea-ice areas in front of Ronne and Brunt Ice Shelf to be the most active with an annual average polynya area of 3018 ± 1298 and 3516 ± 1420 km 2 as well as an accumulated volume ice production of 31 ± 13 and 31 ± 12 km 3 , respectively. For the remaining four regions, estimates amount to 421 ± 294 km 2 and 4 ± 3 km 3 (Antarctic Peninsula), 1148 ± 432 km 2 and 12 ± 5 km 3 (Iceberg A23A), 901 ± 703 km 2 and 10 ± 8 km 3 (Filchner Ice Shelf) as well as 499 ± 277 km 2 and 5 ± 2 km 3 (Coats Land). Our findings are discussed in comparison to recent studies based on coupled sea-ice/ocean models and passive-microwave satellite imagery, each investigating different parts of the Southern Weddell Sea.

Description: Climatic controls and climate proxy potential of Lewis Glacier, Mt Kenya The Cryosphere Discussions, 9, 3887-3924, 2015 Author(s): R. Prinz, L. I. Nicholson, T. Mölg, W. Gurgiser, and G. Kaser The Lewis Glacier on Mt Kenya is one of the best studied tropical glaciers and has experienced considerable retreat since a maximum extent in the late 19th century (L19). From distributed mass and energy balance modelling, this study evaluates the current sensitivity of the surface mass and energy balance to climatic drivers, explores climate conditions under which the L19 maximum extent might have sustained, and discusses the potential for using the glacier retreat to quantify climate change. Multiyear meteorological measurements at 4828 m provide data for input, optimization and evaluation of a spatially distributed glacier mass balance model to quantify the exchanges of energy and mass at the glacier–atmosphere interface. Currently the glacier loses mass due to the imbalance between insufficient accumulation and enhanced melt, because radiative energy gains cannot be compensated by turbulent energy sinks. Exchanging model input data with synthetic climate scenarios, which were sampled from the meteorological measurements and account for coupled climatic variable perturbations, reveal that the current mass balance is most sensitive to changes in atmospheric moisture (via its impact on solid precipitation, cloudiness and surface albedo). Positive mass balances result from scenarios with an increase of annual (seasonal) accumulation of 30 % (100 %), compared to values observed today, without significant changes in air temperature required. Scenarios with lower air temperatures are drier and associated with lower accumulation and increased net radiation due to reduced cloudiness and albedo. If the scenarios currently producing positive mass balances are applied to the L19 extent, negative mass balances are the result, meaning that the conditions required to sustain the glacier in its L19 extent are not reflected in today's observations. Alternatively, a balanced mass budget for the L19 extent can be explained by changing model parameters that imply a distinctly different coupling between the glacier's local surface-air layer and its surrounding boundary-layer. This result underlines the difficulty of deriving paleoclimates for larger glacier extents on the basis of modern measurements of small glaciers.

Description: The electrical self-potential method is a non-intrusive snow-hydrological sensor The Cryosphere Discussions, 9, 4437-4457, 2015 Author(s): S. S. Thompson, B. Kulessa, R. L. H. Essery, and M. P. Lüthi Our ability to measure, quantify and assimilate hydrological properties and processes of snow in operational models is disproportionally poor compared to the significance of seasonal snowmelt as a global water resource and major risk factor in flood and avalanche forecasting. Encouraged by recent theoretical, modelling and laboratory work, we show here that the diurnal evolution of aerially-distributed self-potential magnitudes closely track those of bulk meltwater fluxes in melting in-situ snowpacks at Rhone and Jungfraujoch glaciers, Switzerland. Numerical modelling infers temporally-evolving liquid water contents in the snowpacks on successive days in close agreement with snow-pit measurements. Muting previous concerns, the governing physical and chemical properties of snow and meltwater became temporally invariant for modelling purposes. Because measurement procedure is straightforward and readily automated for continuous monitoring over significant spatial scales, we conclude that the self-potential geophysical method is a highly-promising non-intrusive snow-hydrological sensor for measurement practice, modelling and operational snow forecasting.

Description: Tilt error in cryospheric surface radiation measurements at high latitudes: a model study The Cryosphere Discussions, 9, 4355-4376, 2015 Author(s): W. S. Bogren, J. F. Burkhart, and A. Kylling We have evaluated the magnitude and makeup of error in cryospheric radiation observations due to small sensor misalignment in in-situ measurements of solar irradiance. This error is examined through simulation of diffuse and direct irradiance arriving at a detector with a cosine-response foreoptic. Emphasis is placed on assessing total error over the solar shortwave spectrum from 250 to 4500 nm, as well as supporting investigation over other relevant shortwave spectral ranges. The total measurement error introduced by sensor tilt is dominated by the direct component. For a typical high latitude albedo measurement with a solar zenith angle of 60°, a sensor tilted by 1, 3, and 5° can respectively introduce up to 2.6, 7.7, and 12.8 % error into the measured irradiance and similar errors in the derived albedo. Depending on the daily range of solar azimuth and zenith angles, significant measurement error can persist also in integrated daily irradiance and albedo.

Description: Halogen-based reconstruction of Russian Arctic sea ice area from the Akademii Nauk ice core (Severnaya Zemlya) The Cryosphere Discussions, 9, 4407-4436, 2015 Author(s): A. Spolaor, T. Opel, J. R. McConnell, O. J. Maselli, G. Spreen, C. Varin, T. Kirchgeorg, D. Fritzsche, and P. Vallelonga The role of sea ice in the Earth climate system is still under debate, although it is known to influence albedo, ocean circulation, and atmosphere-ocean heat and gas exchange. Here we present a reconstruction of AD 1950 to 1998 sea ice in the Laptev Sea based on the Akademii Nauk ice core (Severnaya Zemlya, Russian Arctic). The halogens bromine (Br) and iodine (I) are strongly influenced by sea ice processes. Bromine reacts with the sea ice surface in auto-catalyzing "Bromine explosion" events causing an enrichment of the Br / Na ratio and the bromine excess (Br exc ) in snow compared to that in seawater. Iodine is emitted from algal communities growing under sea ice. The results suggest a connection between Br exc and spring sea ice area, as well as a connection between iodine concentration and summer sea ice area. These two halogens are therefore good candidates for extended reconstructions of past sea ice changes in the Arctic.

Description: How to produce the difference data of the two temporal images is a crucial factor in image change detection. In this letter, we propose multicontextual mutual information data (MMID) based on the bivariate Gaussian distribution (BGD) for synthetic aperture radar (SAR) image change detection and illustrate their superiorities over the classical difference data. MMID, which are an improved form of image spatial mutual information, are constructed based on the quadrilateral Markov random field (QMRF) and can be factored into the linear combination of the entropies. Then to adapt MMID to the change detection, we construct the 2-D entropies based on the BGD. In this way, MMID are able to capture the intertemporal statistical dependence of the two temporal images and thus can be taken as the feature-level difference data rather than the pixel-level data. The maximum-likelihood method, the automatic threshold method, and the Markov random field method are performed on the MMID of the real two temporal SAR images for the change detection. Experimental results demonstrate the superiorities of MMID over the traditional difference data.

Description: Spectral unmixing has been a popular technique for analyzing remotely sensed hyperspectral images. The goal of unmixing is to find a collection of pure spectral constituents (called endmembers ) that can explain each (possibly mixed) pixel of the scene as a combination of endmembers, weighted by their coverage fractions in the pixel or abundances . Over the last years, many algorithms have been presented to address the three main parts of the spectral unmixing chain: 1) estimation of the number of endmembers; 2) identification of the endmember signatures; and 3) estimation of the per-pixel fractional abundances. However, to date, there is no standardized tool that integrates these algorithms in a unified framework. In this letter, we present HyperMix, an open-source tool for spectral unmixing that integrates different approaches for spectral unmixing and allows building unmixing chains in graphical fashion, so that the end-user can define one or several spectral unmixing chains in fully configurable mode. HyperMix provides efficient implementations of most of the algorithms used for spectral unmixing, so that the tool automatically recognizes if the computer has a graphics processing unit (GPU) available and optimizes the execution of these algorithms in the GPU. This allows for the execution of spectral unmixing chains on large hyperspectral scenes in computationally efficient fashion. The tool is available online from http://hypercomphypermix.blogspot.com.es and has been validated with real hyperspectral scenes, providing state-of-the-art unmixing results.

Description: Seismic signals are nonlinear, and the seismic state-space model can be described as a nonlinear system. The particle filter (PF) method, as an effective method for estimating the state of a nonlinear system, can be applied to deal with seismic random noise attenuation. However, PF suffers from sample impoverishment caused by resampling, which results in serious loss of valid seismic information and leads to inaccurate representation of the reflected signal. To address the impoverishment issue and to further improve the particle quality, we propose a novel method to suppress seismic random noise—the adaptive fission particle filter (AFPF). In AFPF, all the particles undergo a fission process and produce “offspring” particles to maintain particle diversity. To implement the adaptation and to monitor the degree of fission, we apply a fission factor, which takes into account weights that indicate the quality of the particles. This leads to significant improvements in the particle quality, i.e., the proportion of highly weighted particles is increased. The effective seismic information provided by the resulting particles reproduces the true signal more reliably, reducing the bias of PF. In addition, we establish a dynamic state-space model suitable for seismic signals. Experimental results on synthetic records and field data illustrate the superior performance of AFPF in noise attenuation and reflected signal preservation compared with the PF.

Description: The primary objective of the TanDEM-X mission is the generation of a global high-precision digital elevation model (DEM) by using synthetic aperture radar interferometry. This letter presents the developed strategy for estimating the relative height error of the TanDEM-X DEM on a global scale. The mosaicking process of the final DEM combines all acquisitions at full resolution and is expected to be finished by late 2016. On the other hand, global mosaics can be generated starting from quicklook images already available for each single input data take. These downsized mosaics are operationally used to analyze the performance improvement that can be achieved by combining multiple acquisitions over the same ground areas and are a powerful mean for optimizing further acquisition planning. This letter reports the expected global performance of the final TanDEM-X product in advance of the full-resolution DEM. Knowledge of the global status of the TanDEM-X DEM relative height error is fundamental for optimizing the acquisition strategy and, therefore, the final performance and represents a valuable input for the scientific community as well as for selecting suitable areas for further interferometric experiments on a global scale.

Description: Support vector machines (SVMs) have been applied to land cover classification, and a number of studies have demonstrated their ability to increase classification accuracy. The high correlation between the data set and SVM training model parameters indicates the high performance of the classification model. To improve the correlation, research has focused on the integration of SVMs and other algorithms for data set selection and SVM training model parameter estimation. This letter proposes a novel method, based on a particle filter (PF), of estimating SVM training model parameters according to an observation system. By treating the SVM training function as the observation system of the PF, the new method automatically updates the SVM training model parameters to values that are more appropriate for the data set and can provide a better classification model than can the original model, wherein the parameters are set by trial and error. Various experiments were conducted using Radarsat-2 synthetic aperture radar data from the 2011 Thailand flood. The proposed method provides superior performance and a more accurate analysis compared with the standard SVM.

Description: Maize is a widely planted crop in China and in other areas of the world and plays an important role in grain production. Monitoring the growth status of maize using remote sensing technology is an important component of precision agriculture and height, as a crucial growth indicator for maize, can be retrieved from light detection and ranging (LIDAR) data. However, height extraction for crops, such as maize using airborne laser scanning point clouds results in a great number of uncertainties and challenges. Here, airborne full-waveform LIDAR data were used to extract maize height. In the first step, a workflow was designed based on the Gold deconvolution algorithm combined with a basic data process technique. The method was then tested and was determined to be effective for capturing the portion of the waveform interacting with the tops of vegetation, characterized by lower amplitude stemming from the ground. Therefore, the number of second returns from point clouds was dramatically increased. During the experiment, the number of point clouds increased nearly 50% for three of the four maize plots, as compared with the original point clouds. Compared with the commonly used Gaussian fitting algorithm, the deconvolution algorithm had the advantage of extracting an accurate position for overlapping weak signals. The height percentiles indicated that the original and Gaussian decomposition derived point clouds data underestimated and deconvolution algorithm can accurately reflect the true height of maize, particularly for the 75% and 95% height percentiles.

Description: The Moderate Resolution Imaging Spectroradiometer (MODIS) is an important data set in global burned-area mapping. The MODIS global burned-area product has a coarse spatial resolution at approximately 500 m, which often introduces errors to the mapped burned areas. In this letter, a novel subpixel mapping (SPM) approach was proposed to produce burned-area maps at the fine spatial resolution similar to Landsat imagery, by exploring the spectral and spatial information provided by the second and fifth bands of MODIS. The proposed SPM approach aims to refine the estimate of burned areas, which have been detected by the MODIS global burned-area product. The performance of the proposed SPM approach was assessed with an experiment area containing six burned areas, by comparing with the MODIS burned-area product MCD45. The result shows that the average omission error decreased from 52.26% for MCD45 to 16.74% for SPM, and the average commission error decreased from 21.76% for MCD45 to 12.54% for SPM. The kappa value increased from 0.5583 for MCD45 to 0.8756 for SPM, indicating that the proposed SPM approach is effective in reducing the influence of the coarse spatial resolution of MODIS imagery in mapping a burned area and refining existing global burned-area products.

Description: Radiometric distortions caused by rugged terrain make the classification of forest types from satellite imagery a challenge. Various band-specific topographic normalization models are expected to eliminate or reduce these effects. The quality of these models also depends on the approach to estimate empirical parameters. Generally, a global estimation of these parameters from a whole satellite image is simple, but it may tend to overcorrection, particularly for larger areas. A land-cover-specific method usually performs better, but it requires obtaining a priori land classification, which presents another challenge in many cases. Empirical parameters can be directly estimated from local pixels in a given window. In this letter, we propose and evaluate a central-pixel-based parameter estimation method for topographic normalization using local window pixels. We tested the method with Landsat 8 imagery and the Advanced Spaceborne Thermal Emission and Reflection Radiometer Global Digital Elevation Model (ASTER GDEM) in very rough terrain with diverse forest types. Visual comparison and statistical analyses showed that the proposed method performed better at a range of window sizes compared with an uncorrected image or with a global parameter estimation approach. The intraclass spectral variability of each forest type has been reduced significantly, and it can yield higher accuracy of forest type classification. The proposed method does not require the a priori knowledge of land covers. Its simplicity and robustness suggest that this method has the potential to be a standard preprocessing approach for optical satellite imagery, particularly for rough terrain.

Description: A moving point approach to model shallow ice sheets: a study case with radially-symmetrical ice sheets The Cryosphere Discussions, 9, 4237-4270, 2015 Author(s): B. Bonan, M. J. Baines, N. K. Nichols, and D. Partridge Predicting the evolution of ice sheets requires numerical models able to accurately track the migration of ice sheet continental margins or grounding lines. We introduce a physically-based moving point approach for the flow of ice sheets based on the conservation of local masses. This allows the ice sheet margins to be tracked explicitly and the waiting time behaviours to be modelled efficiently. A finite difference moving point scheme is derived and applied in a simplified context (continental radially-symmetrical shallow ice approximation). The scheme, which is inexpensive, is validated by comparing the results with moving-margin exact solutions and steady states. In both cases the scheme is able to track the position of the ice sheet margin with high precision.

Description: Fram Strait spring ice export and September Arctic sea ice The Cryosphere Discussions, 9, 4205-4235, 2015 Author(s): M. H. Halvorsen, L. H. Smedsrud, R. Zhang, and K. Kloster The Arctic Basin exports ~ 10 % of the sea ice area southwards annually through Fram Strait. A larger than normal export decreases the remaining mean thickness and ice area. A new updated timeseries from 1979–2013 of Fram Strait sea ice area export shows an overall increase until today, and that more than 1 million km 2 has been exported annually in recent years. The new timeseries has been constructed from high resolution radar satellite imagery of sea ice drift across 79° N from 2004–2013, regressed on the observed cross-strait surface pressure difference, and shows an increasing trend of 7 % per decade. The trend is caused by higher southward ice drift speeds due to stronger southward geostrophic winds, largely explained by increasing surface pressure on Greenland. Spring and summer area export increases more (~ 14 % per decade) than in autumn and winter, and these export anomalies have a large influence on the following September mean ice extent.

Description: Many methods have been developed to detect damaged buildings due to earthquake. However, little attention has been paid to analyze slightly affected buildings. In this letter, an unsupervised method is presented to detect earthquake-triggered “ roof-holes ” on rural houses from unmanned aerial vehicle (UAV) images. First, both orthomosaic and gradient images are generated from a set of UAV images. Then, a modified Chinese restaurant franchise model is used to learn an unsupervised model of the geo-object classes in the area by fusing both oversegmented orthomosaic and gradient images. Finally, “roof-holes” on rural houses are detected using the learned model. The performance of the proposed method is evaluated in terms of both qualitative and quantitative indexes.

Description: Automatic target generation process (ATGP) has been used in a wide range of applications in hyperspectral image analysis. It performs a sequence of orthogonal subspace projections to extract potential targets of interest. This letter presents a recursive version of the ATGP, which is referred to as the recursive ATGP (RATGP) and has three advantages over the ATGP as follows: 1) there is no need of inverting a matrix as the ATGP does for finding each new target; 2) there is a significant reduction in the computational complexity in the hardware design due to its recursive structure; and 3) there is an automatic stopping rule that can be derived by the Neyman–Pearson detection theory to terminate the algorithm.

Description: In this letter, we propose a novel automatic algorithm for road extraction from remote sensing images. The algorithm includes low- and high-level processing. In the low-level processing, we determine a normalized second derivative map of road profiles of a generalized bar shape, which is width invariant and contrast proportional, and accordingly obtain initial road center pixels. In the high-level processing, using the map and initial center pixels, we initially determine road segments. The segments are then locally refined using their orientation randomness and length-to-width ratio and further refined via global graph-cut optimization. A final road network is thereby extracted in a robust manner. Experimental results demonstrate that the proposed algorithm provides noticeably more robust and higher road extraction performance in various images compared with the existing algorithms.

Description: Remote sensing images often need to be coded and/or transmitted with constrained computational resources. Among other features, such images commonly have high spatial, spectral, and bit-depth resolution, which may render difficult their handling. This letter introduces an embedded quantization scheme based on two-step scalar deadzone quantization (2SDQ) that enhances the quality of transmitted images when coded with a constrained number of bits. The proposed scheme is devised for use in JPEG2000. It is named cell-based 2SDQ since it uses cells, i.e., small sets of wavelet coefficients within the codeblocks defined by JPEG2000. Cells permit a finer discrimination of coefficients in which to apply the proposed quantizer. Experimental results indicate that the proposed scheme is especially beneficial for high bit-depth hyperspectral images.

Description: High suspended solid (SS) concentrations in coastal waters are aesthetically undesirable, and adversely affect fisheries and coastal ecosystems. Environmental agencies usually require frequent measurements of SS over coastal regions at a spatially detailed level for water quality assessment and control. To develop a method for SS estimation in the complex coastal waters of Hong Kong, an archive of 57 Landsat Thematic Mapper (TM), Enhanced Thematic Mapper Plus (ETM+), and HJ-1 A/B Charged Couple Device (CCD) images over a 13-year period from January 2000 to December 2012 was used. Atmospherically corrected Landsat TM/ETM+ and HJ-1 A/B CCD bands 1–4 along with 240 in situ field samples of SS concentration collected within 2 h of image acquisition, were used to develop and validate regression models over a wide range of SS concentrations from 0.5–56.0 mg/L. The best representation of actual SS concentrations was given by the log-transformed combination of Band 2 (Green, 0.52–0.60 $mumbox{m}$ ) and Band 3 (Red, 0.63–0.69 $mumbox{m} $ ), with correlation coefficient (R) of 0.85, root-mean-square error of 2.60 mg/L and mean absolute error of 2.04 mg/L. This is attributed to the sensitivity of SS to green and red wavelengths specific to the characteristic refractive index and grain size of SS found in Hong Kong waters. This letter is considered more robust than previous studies, due to the much larger number of images and in situ samples used for model development and validation, as well as the different times of year and wide range of SS concentrations investigated.

Description: The quantitative estimation of the fractional cover of carbonate rock (CR) is critical for natural resource management and ecological conservation in karst areas. Based on the analysis of spectral properties of CR together with other land cover types, we proposed two CR indices (CRIs) and established the model that represents the relationships between the CRIs and the fractional cover of CR. Then, the fractional cover of CR was estimated by using the developed model. Experimental results on Landsat-8 Operational Land Imager images acquired at Southwestern China demonstrated the effectiveness of the developed model. Compared with other indices, the proposed CRIs show the highest correlations with the fractional cover of CR.

Description: Tomography-based observation of sublimation and snow metamorphism under temperature gradient and advective flow The Cryosphere Discussions, 9, 4845-4864, 2015 Author(s): P. P. Ebner, M. Schneebeli, and A. Steinfeld Snow at or close to the surface commonly undergoes temperature gradient metamorphism under advective flow, which alters its microstructure and physical properties. Time-lapse X-ray micro-tomography is applied to investigate the structural dynamics of temperature gradient snow metamorphism exposed to an advective airflow in controlled laboratory conditions. The sublimation of water vapor for saturated air flowing across the snow sample was experimentally determined via variations of the porous ice structure. The results showed that the exothermic gas-to-solid phase change is favorable vis-a-vis the endothermic solid-to-gas phase change, thus leading to more ice deposition than ice sublimation. Sublimation has a marked effect on the structural change of the ice matrix but diffusion of water vapor in the direction of the temperature gradient counteracted the mass transport of advection. Therefore, the total net ice change was negligible leading to a constant porosity profile. However, the strong reposition process of water molecules on the ice grains is relevant for atmospheric chemistry.

Description: Glacier dynamics over the last quarter of a century at Jakobshavn Isbræ The Cryosphere Discussions, 9, 4865-4892, 2015 Author(s): I. S. Muresan, S. A. Khan, A. Aschwanden, C. Khroulev, T. Van Dam, J. Bamber, M. R. van den Broeke, B. Wouters, P. Kuipers Munneke, and K. H. Kjær Observations over the past two decades show substantial ice loss associated with the speedup of marine terminating glaciers in Greenland. Here we use a regional 3-D outlet glacier model to simulate the behaviour of Jakobshavn Isbræ (JI) located in west Greenland. Using atmospheric and oceanic forcing we tune our model to reproduce the observed frontal changes of JI during 1990–2014. We identify two major accelerations. The first occurs in 1998, and is triggered by moderate thinning prior to 1998. The second acceleration, which starts in 2003 and peaks in summer 2004, is triggered by the final breakup of the floating tongue, which generates a reduction in buttressing at the JI terminus. This results in further thinning, and as the slope steepens inland, sustained high velocities have been observed at JI over the last decade. As opposed to other regions on the Greenland Ice Sheet (GrIS), where dynamically induced mass loss has slowed down over recent years, both modelled and observed results for JI suggest a continuation of the acceleration in mass loss. Further, we find that our model is not able to capture the 2012 peak in the observed velocities. Our analysis suggests that the 2012 acceleration of JI is likely the result of an exceptionally long melt season dominated by extreme melt events. Considering that such extreme surface melt events are expected to intensify in the future, our findings suggest that the 21st century projections of the GrIS mass loss and the future sea level rise may be larger than predicted by existing modelling results.

Description: About the consistency between Envisat and CryoSat-2 radar freeboard retrieval over Antarctic sea ice The Cryosphere Discussions, 9, 4893-4923, 2015 Author(s): S. Schwegmann, E. Rinne, R. Ricker, S. Hendricks, and V. Helm Knowledge about Antarctic sea-ice volume and its changes over the past decades has been sparse due to the lack of systematic sea-ice thickness measurements in this remote area. Recently, first attempts have been made to develop a sea-ice thickness product over the Southern Ocean from space-borne radar altimetry and results look promising. Today, more than 20 years of radar altimeter data are potentially available for such products. However, data come from different sources, and the characteristics of individual sensors differ. Hence, it is important to study the consistency between single sensors in order to develop long and consistent time series over the potentially available measurement period. Here, the consistency between freeboard measurements of the Radar Altimeter 2 on-board Envisat and freeboard measurements from the Synthetic-Aperture Interferometric Radar Altimeter on-board CryoSat-2 is tested for their overlap period in 2011. Results indicate that mean and modal values are comparable over the sea-ice growth season (May–October) and partly also beyond. In general, Envisat data shows higher freeboards in the seasonal ice zone while CryoSat-2 freeboards are higher in the perennial ice zone and near the coasts. This has consequences for the agreement in individual sectors of the Southern Ocean, where one or the other ice class may dominate. Nevertheless, over the growth season, mean freeboard for the entire (regional separated) Southern Ocean differs generally by not more than 2 cm (5 cm, except for the Amundsen/Bellingshausen Sea) between Envisat and CryoSat-2, and the differences between modal freeboard lie generally within ±10 cm and often even below.

Description: Brief Communication: Upper air relaxation in RACMO2 significantly improves modelled interannual SMB variability in Antarctica The Cryosphere Discussions, 9, 4981-4995, 2015 Author(s): W. J. van de Berg and B. Medley The regional climate model RACMO2 has been a powerful tool for improving SMB estimates from GCMs or reanalyses. However, new yearly SMB observations for West Antarctica show that the modelled interannual variability in SMB is poorly simulated by RACMO2, in contrast to ERA-Interim, which resolves this variability well. In an attempt to remedy RACMO2 performance, we included additional upper air relaxation (UAR) in RACMO2. With UAR, the correlation to observations is similar for RACMO2 and ERA-Interim. The spatial SMB patterns and ice sheet integrated SMB modelled using UAR remain very similar to the estimates of RACMO2 without UAR. We only observe an upstream smoothing of precipitation in regions with very steep topography like the Antarctic Peninsula. We conclude that UAR is a useful improvement for RCM simulations, although results in regions with steep topography should be treated with care.

Description: Modeling Antarctic subglacial lake filling and drainage cycles The Cryosphere Discussions, 9, 6545-6579, 2015 Author(s): C. F. Dow, M. A. Werder, S. Nowicki, and R. T. Walker The growth and drainage of active subglacial lakes in Antarctica has previously been inferred from analysis of ice surface altimetry data. We use a subglacial hydrology model applied to a synthetic Antarctic ice stream to determine internal controls on the filling and drainage of subglacial lakes and their impact on ice stream dynamics. Our model outputs suggest that the highly constricted subglacial environment of the ice stream, combined with relatively high rates of water flow funneled from large catchments, can combine to create a system exhibiting slow-moving pressure waves. Over a period of years, the accumulation of water in the ice stream onset region results in a buildup of pressure creating temporary channels, which then evacuate the excess water. This increased flux of water through the ice stream drives lake growth. As the water body builds up, it too steepens the hydraulic gradient and allows greater flux out of the overdeepened lake basin. Eventually this flux is large enough to create channels that cause the lake to drain. Due to the presence of the channels, the drainage of the lake causes high water pressures around 50 km downstream of the lake rather than immediately in the vicinity of the overdeepening. Following lake drainage, channels again shut down. Lake drainage depends on the internal hydrological development in the wider system and therefore does not directly correspond to a particular water volume or depth. This creates a highly temporally and spatially variable system, which is of interest for assessing the importance of subglacial lakes in ice stream hydrology and dynamics.

Description: Assessment of error in satellite derived lead fraction in Arctic The Cryosphere Discussions, 9, 6315-6344, 2015 Author(s): N. Ivanova, P. Rampal, and S. Bouillon Leads within consolidated sea ice control heat exchange between the ocean and the atmosphere during winter thus constituting an important climate parameter. These narrow elongated features occur when sea ice is fracturing under the action of wind and currents, reducing the local mechanical strength of the ice cover, which in turn impact the sea ice drift patterns. This makes a high quality lead fraction (LF) dataset to be in demand for sea ice model evaluation, initialization and for assimilation of such data in regional models. In this context, the available LF dataset retrieved from satellite passive microwave observations (Advanced Microwave Scanning Radiometer – Earth Observing System, AMSR-E) is of great value, providing pan-Arctic light- and cloud-independent daily coverage since 2002. Here we quantify errors in this dataset using accurate LF estimates retrieved from Synthetic Aperture Radar (SAR) images employing a threshold technique, also introduced in this work. We find a consistent overestimation by a factor of 2–4 of the LF estimates in the AMSR-E LF product. We show for a data sample from the AMSR-E LF dataset that a simple adjustment of the tie points used in the method to estimate the LF can reduce the pixel-wise error by a factor of 2 on average. Applying such adjustment to the full dataset may thus significantly increase the quality and value of the original dataset.

Description: Acquisition of isotopic composition for surface snow in East Antarctica and the links to climatic parameters The Cryosphere Discussions, 9, 6275-6313, 2015 Author(s): A. Touzeau, A. Landais, B. Stenni, R. Uemura, K. Fukui, S. Fujita, S. Guilbaud, A. Ekaykin, M. Casado, E. Barkan, B. Luz, O. Magand, G. Teste, E. Le Meur, M. Baroni, J. Savarino, I. Bourgeois, and C. Risi The isotopic composition of oxygen and hydrogen in ice cores are invaluable tools for the reconstruction of past climate variations. Used alone, they give insights into the variations of the local temperature, whereas taken together they can provide information on the climatic conditions at the point of origin of the moisture. However, recent analyses of snow from shallow pits indicate that the climatic signal can become erased in very low accumulation regions, due to local processes of snow reworking. The signal to noise ratio decreases and the climatic signal can then only be retrieved using stacks of several snow pits. Obviously, the signal is not completely lost at this stage, otherwise it would be impossible to extract valuable climate information from ice cores as has been done, for instance, for the last glaciation. To better understand how the climatic signal is passed from the precipitation to the snow, we present here results from varied snow samples from East Antarctica. First, we look at the relationship between isotopes and temperature from a geographical point of view, using results from three traverses across Antarctica, to see how the relationship is built up through the distillation process. We also take advantage of these measures to see how second order parameters (d-excess and 17 O-excess) are related to δ 18 O and how they are controlled. d-excess increases in the interior of the continent (i.e. when δ 18 O decreases), due to the distillation process, whereas 17 O-excess decreases in remote areas, due to kinetic fractionation at low temperature. In both cases, these changes are associated with the loss of original information regarding the source. Then, we look at the same relationships in precipitation samples collected over one year at Dome C and Vostok, as well as in surface snow at Dome C. We note that the slope of the δ 18 O / T relationship decreases in these samples compared to those from the traverses, and thus advocate caution when using spatial slopes for past climate reconstruction. The second-order parameters behave in the same way in the precipitation as in the surface snow from traverses, indicating that similar processes are active. Finally we check if the same relationships between δ 18 O and second-order parameters are also found in the snow from four snow pits. While the d-excess remains opposed to δ 18 O in most snow pits, the 17 O-excess is no longer positively correlated to δ 18 O and even shows anti-correlation to δ 18 O at Vostok. This may be due to a stratospheric influence at this site and/or to post-deposition processes.

Description: Provides a listing of the editors, board members, and current staff for this issue of the publication.

Description: Characterizing Arctic sea ice topography using high-resolution IceBridge data The Cryosphere Discussions, 9, 6495-6543, 2015 Author(s): A. A. Petty, M. C. Tsamados, N. T. Kurtz, S. L. Farrell, T. Newman, J. P. Harbeck, D. L. Feltham, and J. A. Richter-Menge We present an analysis of Arctic sea ice topography using high resolution, three-dimensional, surface elevation data from the Airborne Topographic Mapper, flown as part of NASA's Operation IceBridge mission. Surface features in the sea ice cover are detected using a newly developed surface feature picking algorithm. We derive information regarding the height, volume and geometry of surface features from 2009–2014 within the Beaufort/Chukchi and Central Arctic regions. The results are delineated by ice type to estimate the topographic variability across first-year and multi-year ice regimes. The results demonstrate that Arctic sea ice topography exhibits significant spatial variability, mainly driven by the increased surface feature height and volume (per unit area) of the multi-year ice that dominates the Central Arctic region. The multi-year ice topography exhibits greater interannual variability compared to the first-year ice regimes, which dominates the total ice topography variability across both regions. The ice topography also shows a clear coastal dependency, with the feature height and volume increasing as a function of proximity to the nearest coastline, especially north of Greenland and the Canadian Archipelago. A strong correlation between ice topography and ice thickness (from the IceBridge sea ice product) is found, using a square-root relationship. The results allude to the importance of ice deformation variability in the total sea ice mass balance, and provide crucial information regarding the tail of the ice thickness distribution across the western Arctic. Future research priorities associated with this new dataset are presented and discussed, especially in relation to calculations of atmospheric form drag.

Description: Provides a listing of the editors, board members, and current staff for this issue of the publication.

Description: Committed near-future retreat of Smith, Pope, and Kohler Glaciers inferred by transient model calibration The Cryosphere Discussions, 9, 4459-4498, 2015 Author(s): D. Goldberg, P. Heimbach, I. Joughin, and B. Smith A glacial flow model of Smith, Pope and Kohler Glaciers has been calibrated by means of inverse methods against time-varying, annualy resolved observations of ice height and velocities, covering the period 2002 to 2011. The inversion – termed "transient calibration" – produces an optimal set of time-mean, spatially varying parameters together with a time-evolving state that accounts for the transient nature of observations and the model dynamics. Serving as an optimal initial condition, the estimated state for 2011 is used, with no additional forcing, for predicting grounded ice volume loss and grounding line retreat over the ensuing 30 years. The transiently calibrated model predicts a near-steady loss of grounded ice volume of approximately 21 km 3 a −1 over this period, as well as loss of 33 km 2 a −1 grounded area. We contrast this prediction with one obtained following a commonly used "snapshot" or steady-state inversion, which does not consider time dependence and assumes all observations to be contemporaneous. Transient calibration is shown to achieve a better fit with observations of thinning and grounding line retreat histories, and yields a quantitatively different projection with respect to ice volume loss and ungrounding. Sensitivity studies suggest large near-future levels of unforced, i.e. committed sea level contribution from these ice streams under reasonable assumptions regarding uncertainties of the unknown parameters.

Description: Photopolarimetric retrievals of snow properties The Cryosphere Discussions, 9, 3055-3074, 2015 Author(s): M. Ottaviani, B. van Diedenhoven, and B. Cairns Polarimetric observations of snow surfaces, obtained in the 410–2264 nm range with the Research Scanning Polarimeter onboard the NASA ER-2 high-altitude aircraft, are analyzed and presented. These novel measurements are of interest to the remote sensing community because the overwhelming brightness of snow plagues aerosol and cloud retrievals based on air- and space-borne total reflection measurements. The spectral signatures of the polarized reflectance of snow are therefore worthwhile investigating in order to provide guidance for the adaptation of algorithms currently employed for the retrieval of aerosol properties over soil and vegetated surfaces. At the same time, the increased information content of polarimetric measurements allows for a meaningful characterization of the snow medium. In our case, the grains are modeled as hexagonal prisms of variable aspect ratios and microscale roughness, yielding retrievals of the grains' scattering asymmetry parameter, shape and size. The results agree with our previous findings based on a more limited dataset, with the majority of retrievals leading to moderately rough crystals of extreme aspect ratios, for each scene corresponding to a single value of the asymmetry parameter.

Description: A prognostic model of the sea ice floe size and thickness distribution The Cryosphere Discussions, 9, 2955-2997, 2015 Author(s): C. Horvat and E. Tziperman Sea ice exhibits considerable seasonal and longer-term variations in extent, concentration, thickness and age, and is characterized by a complex and continuously changing distribution of floe sizes and thicknesses. Models of sea ice used in current climate models keep track of its concentration and of the distribution of ice thicknesses, but do not account for the floe size distribution and its potential effects on air–sea exchange and sea-ice evolution. Accurately capturing sea-ice variability in climate models may require a better understanding and representation of the distribution of floe sizes and thicknesses. We develop and demonstrate a model for the evolution of the joint sea-ice floe size and thickness distribution that depends on atmospheric and oceanic forcing fields. The model accounts for effects due to multiple processes that are active in the marginal and seasonal ice zones: freezing and melting along the lateral side and base of floes, mechanical interactions due to floe collisions (ridging and rafting) and sea-ice fracture due to swell propagation into the ice pack. The model is then examined and demonstrated in a series of idealized test cases.

Description: ENSO influence on surface energy and mass balance at Shallap Glacier, Cordillera Blanca, Peru The Cryosphere Discussions, 9, 2999-3053, 2015 Author(s): F. Maussion, W. Gurgiser, M. Großhauser, G. Kaser, and B. Marzeion The El Niño/Southern Oscillation (ENSO) is a major driver of climate variability in the tropical Andes, where recent Niño and Niña events left an observable footprint on glacier mass balance. The nature and strength of the relationship between ENSO and glacier mass balance, however, varies between regions and time periods, leaving several unanswered questions about its exact mechanisms. The starting point of this study is a four-year long time series of distributed surface energy and mass balance (SEB/SMB) calculated using a process-based model driven by observations at Shallap Glacier (Cordillera Blanca, Peru). These data are used to calibrate a regression-based downscaling model that links the local SEB/SMB fluxes to atmospheric reanalysis variables on a monthly basis, allowing an unprecedented quantification of the ENSO influence on the SEB/SMB at climatological time scales (1980–2013, ERA-Interim period). We find a stronger and steadier anti-correlation between pacific sea surface temperature (SST) and glacier mass balance than previously reported. This relationship is most pronounced during the wet season (December–May) and at low altitudes where Niño (Niña) events are accompanied with a snowfall deficit (excess) and a higher (lower) radiation energy input. We detect a weaker but significant ENSO anti-correlation with total precipitation (Niño dry signal) and positive correlation with the sensible heat flux, but find no ENSO influence on sublimation. Sensitivity analyses comparing several downscaling methods and reanalysis datasets resulted in stable mass balance correlations with pacific SST but also revealed large uncertainties in computing the mass balance trend of the last decades. The newly introduced open-source downscaling tool can be applied easily to other glaciers in the tropics, opening new research possibilities on even longer time scales.

Description: Changing surface–atmosphere energy exchange and refreezing capacity of the lower accumulation area, west Greenland The Cryosphere Discussions, 9, 2867-2913, 2015 Author(s): C. Charalampidis, D. van As, J. E. Box, M. R. van den Broeke, W. T. Colgan, M. MacFerrin, H. Machguth, and P. P. Smeets We present five years (2009–2013) of automatic weather station measurements from the lower accumulation area (1840 m a.s.l.) of the ice sheet in the Kangerlussuaq region, western Greenland. Here, the summers of 2010 and 2012 were both exceptionally warm, but only 2012 resulted in a strongly negative surface mass budget (SMB) and surface meltwater runoff. The observed runoff was due to a large ice fraction in the upper 10 m of firn that prevented meltwater from percolating to available pore volume below. Analysis reveals a relatively low 2012 summer albedo of ~0.7 as meltwater was present at the surface. Consequently, during the 2012 melt season the surface absorbed 29% (213 MJ m -2 ) more solar radiation than the average of all other years. A surface energy balance model is used to evaluate the seasonal and interannual variability of all surface energy fluxes. The model reproduces the observed melt rates as well as the SMB for each season. A sensitivity test reveals that 71% of the additional solar radiation in 2012 was used for melt, corresponding to 36% (0.64 m) of the 2012 surface lowering. The remaining 1.14 m was primarily due to the high atmospheric temperatures up to +2.6 °C daily average, indicating that 2012 would have been a negative SMB year at this site even without the melt-albedo feedback. Longer time series of SMB, regional temperature and remotely sensed albedo (MODIS) show that 2012 was the first strongly negative SMB year with the lowest albedo at this elevation on record. The warm conditions of the last years resulted in enhanced melt and reduction of the refreezing capacity at the lower accumulation area. If high temperatures continue the current lower accumulation area will turn into a region with superimposed ice in coming years.

Description: The relative contributions of calving and surface ablation to ice loss at a lake-terminating glacier The Cryosphere Discussions, 9, 2915-2953, 2015 Author(s): M. Chernos, M. Koppes, and R. D. Moore Bridge Glacier is a lake-terminating glacier in the Coast Mountains of British Columbia and has retreated over 3.55 km since 1972, with the majority of the retreat having occurred since 1991. This retreat is out of proportion to surface melt inferred from regional climate indices, suggesting that it has been driven primarily by calving as the glacier retreated across an over-deepened basin. In order to better understand the primary drivers of mass balance, the relative importance of surface melt and calving is investigated during the 2013 melt season using a distributed energy balance model and time-lapse imagery. Calving is responsible for 23% of the mass loss during the 2013 melt season, and is limited by modest flow speeds and a small terminus cross-section. Calving and summer balance estimates over the last 30 years suggest that calving is consistently a smaller contributor of mass loss relative to surface melt. Although calving is estimated to be responsible for up to 49% of ice loss for individual seasons, averaged over multiple summers it typically accounts for 10 to 25%. Calving has been driven primarily by buoyancy and water depths, and fluxes were greatest between 2005 and 2010 as the glacier retreated over the deepest part of Bridge Lake. These losses are part of a transient stage in the glacier's retreat, and are expected to diminish as the terminus recedes into shallower water. Surface melt is the primary driver of ice loss at Bridge Glacier, and future mass loss and retreat is dependent on governing climatic conditions.

Description: Coastal dynamics and submarine permafrost in shallow water of the central Laptev Sea, East Siberia Pier Paul Overduin, Sebastian Wetterich, Frank Günther, Mikhail N. Grigoriev, Guido Grosse, Lutz Schirrmeister, Hans-Wolfgang Hubberten, and Aleksandr Makarov The Cryosphere, 10, 1449-1462, doi:10.5194/tc-10-1449-2016, 2016 How fast does permafrost warm up and thaw after it is covered by the sea? Ice-rich permafrost in the Laptev Sea, Siberia, is rapidly eroded by warm air and waves. We used a floating electrical technique to measure the depth of permafrost thaw below the sea, and compared it to 60 years of coastline retreat and permafrost depths from drilling 30 years ago. Thaw is rapid right after flooding of the land and slows over time. The depth of permafrost is related to how fast the coast retreats.

Description: Age of the Mt. Ortles ice cores, the Tyrolean Iceman and glaciation of the highest summit of South Tyrol since the Northern Hemisphere Climatic Optimum Paolo Gabrielli, Carlo Barbante, Giuliano Bertagna, Michele Bertó, Daniel Binder, Alberto Carton, Luca Carturan, Federico Cazorzi, Giulio Cozzi, Giancarlo Dalla Fontana, Mary Davis, Fabrizio De Blasi, Roberto Dinale, Gianfranco Dragà, Giuliano Dreossi, Daniela Festi, Massimo Frezzotti, Jacopo Gabrieli, Stephan Galos, Patrick Ginot, Petra Heidenwolf, Theo M. Jenk, Natalie Kehrwald, Donald Kenny, Olivier Magand, Volkmar Mair, Vladimir Mikhalenko, Ping Nan Lin, Klaus Oeggl, Gianni Piffer, Mirko Rinaldi, Ulrich Schotterer, Margit Schwikowski, Roberto Seppi, Andrea Spolaor, Barbara Stenni, Davdi Tonidandel, CChiara Uglietti, Victor Zagorodnov, Thomas Zanoner, and Piero Zennaro The Cryosphere Discuss., doi:10.5194/tc-2016-159,2016 Manuscript under review for TC (discussion: open, 0 comments) New ice cores were extracted from Alto dell'Ortles, the highest glacier of South Tyrol in the Italian Alps, to check whether prehistoric ice, that is coeval to the famous 5200 years old Tyrolean Iceman, is still preserved in this region. Dating of the ice cores confirms the hypothesis and indicates that the drilling site was glaciated since the end of the Northern Hemisphere Climatic Optimum (7000 BP). We also infer that an unprecedented acceleration of the glacier flow has just recently begun.

Description: Combined diurnal variations of discharge and hydrochemistry of the Isunnguata Sermia outlet of the Greenland Ice Sheet give in sight on sub glacial conditions Joseph Graly, Joel Harrington, and Neil Humphrey The Cryosphere Discuss., doi:10.5194/tc-2016-137,2016 Manuscript under review for TC (discussion: open, 0 comments) At a major outlet of the Greenland Ice Sheet in West Greenland, we find that the chemical solutes is the emerging subglacial waters are out of phase with water discharge and can spike in concentration during waning flow. This suggests that the subglacial waters are spreading out across a large area of the glacial bed throughout the day, stimulating chemical weathering beyond the major water distribution channels.

Description: Inversion of geothermal heat flux in a thermomechanically coupled nonlinear Stokes ice sheet model Hongyu Zhu, Noemi Petra, Georg Stadler, Tobin Isaac, Thomas J. R. Hughes, and Omar Ghattas The Cryosphere, 10, 1477-1494, doi:10.5194/tc-10-1477-2016, 2016 We study how well the basal geothermal heat flux can be inferred from surface velocity observations using a thermomechanically coupled nonlinear Stokes ice sheet model. The prospects and limitations of this inversion is studied in two and three dimensional model problems. We also argue that a one-way coupled approach for the adjoint equations motivated by staggered solvers for forward multiphysics problems can lead to an incorrect gradient and premature termination of the optimization iteration.

Description: Increased nitrate and decreased δ 15 N–NO 3 − in the Greenland Arctic after 1940 attributed to North American oil burning Nathan J. Chellman, Meredith G. Hastings, and Joseph R. McConnell The Cryosphere Discuss., doi:10.5194/tc-2016-163,2016 Manuscript under review for TC (discussion: open, 0 comments) This manuscript analyzes the changing sources of nitrate deposition to Greenland since 1760 CE using a dataset consisting of sub-seasonally resolved nitrogen isotopes of nitrate and source tracers. Correlations amongst ion concentration, source tracers, and the δ 15 N–NO 3 − provide evidence of the impact of biomass burning and fossil fuel combustion emissions of nitrogen oxides and suggest that oil combustion is the likely driver of increased nitrate concentration in Greenland ice since 1940 CE.

Description: Slight glacier reduction over the northwestern Tibetan Plateau despite significant recent warming Yetang Wang, Shugui Hou, Wenling An, Hongxi Pang, and Yaping Liu The Cryosphere Discuss., doi:10.5194/tc-2016-165,2016 Manuscript under review for TC (discussion: open, 0 comments) This study further confirms "Pamir–Karakoram–Western-Kunlun-Mountain (northwestern Tibetan Plateau) Glacier Anomaly". Slight glacier reduction over the northwestern Tibetan Plateau may result from more accumulation from increased precipitation in winter which to great extent protects it from mass reductions under climate warming during 1961–2000. Warming slowdown since 2000 happening at this region may further mitigate glacier mass reduction.

Description: Response of seasonal soil freeze depth to climate change across China Xiaoqing Peng, Oliver W. Frauenfeld, Tingjun Zhang, Kang Wang, Bin Cao, Xinyue Zhong, Hang Su, and Cuicui Mu The Cryosphere Discuss., doi:10.5194/tc-2016-129,2016 Manuscript under review for TC (discussion: open, 0 comments) Previous researches paid more attention on the permafrost, e.g. active layer thickness, soil temperature, permafrost area extent, and associated with permafrost degradation leading other changes. However, seasonally frozen ground, vast area extent, did not focus by so much attention. Here, we combined more than 800 observation station data and gridded data to investigate soil freeze depth across China. The results indicate that soil freeze depth increase with climate warming.

Description: Characteristics of an avalanche-feeding and partially debris-covered glacier and its response to atmospheric warming in Mt. Tomor, Tian Shan, China Puyu Wang, Zhongqin Li, and Huilin Li The Cryosphere Discuss., doi:10.5194/tc-2016-138,2016 Manuscript under review for TC (discussion: open, 0 comments) A cirque-valley glacier with complex topography and partially debris-covered area was investigated in the Mt. Tomor, Tian Shan. The glacier is analogous to temperate one on movement and temperature regimes. The strongest ablation and most significant terminus retreat and area reduction occurred at the end of last century and the beginning of this century. Since inhibition of debris cover to melting, it is expected to keep shrinkage in the coming decades, but the terminus retreat is to be slower.

Description: An ice-sheet-wide framework for englacial attenuation from ice-penetrating radar data T. M. Jordan, J. L. Bamber, C. N. Williams, J. D. Paden, M. J. Siegert, P. Huybrechts, O. Gagliardini, and F. Gillet-Chaulet The Cryosphere, 10, 1547-1570, doi:10.5194/tc-10-1547-2016, 2016 Ice penetrating radar enables determination of the basal properties of ice sheets. Existing algorithms assume stationarity in the attenuation rate, which is not justifiable at an ice sheet scale. We introduce the first ice-sheet-wide algorithm for radar attenuation that incorporates spatial variability, using the temperature field from a numerical model as an initial guess. The study is a step toward ice-sheet-wide data products for basal properties and evaluation of model temperature fields.

Description: The inverse problem of synthetic aperture imaging radiometers (SAIRs) has been demonstrated to be not well posed. The regularization methods are crucial for providing unique and stable solutions in the reconstruction of radiometric brightness temperature (BT) maps. Different to deterministic ones, a new approach is presented by referring to the rule of Bayesian inference, providing a probability model of regularized constraints to combat the ill-posedness of finite-dimensional discrete inverse problems. In addition, the SAIR inverse problem can be converted into the probability estimation of the reconstructed BT. Furthermore, in application to both uniformly and nonuniformly spaced arrays, our method can obtain the optimal solution adaptively and avoid the dilemma of choosing the optimal regularization parameter. Finally, simulation results illustrating the effectiveness and performance of the proposed method are provided.

Description: A range-cell-focusing algorithm is proposed in order to improve the quality of the target image. In a synthetic aperture radar (SAR) imaging system, the range resolution depends on the frequency bandwidth and determines the ability to distinguish between targets that are very close to each other. In cases where the resolution and the SNR from the environment are not adequate, targets cannot be accurately visualized. In order to successively classify targets that are close, we are combining an enhanced-multiple-signal-classification spectrum as a weighting function to reproduce the raw data. The proposed algorithm improves classification and separation for close targets while suppressing artifacts in the final images. The targets of interest are stationary point scatterers. The results are obtained from both simulated and experimental data to demonstrate that the proposed algorithm provides better performance than a conventional SAR imaging algorithm, the range migration algorithm.

Description: This letter presents a multiscale edge detection method for multilook polarimetric synthetic aperture radar (PolSAR) images based on the nonsubsampled contourlet transform (NSCT). The NSCT can provide flexible multiscale and directional decomposition. In the multiscale decomposition, the coefficients of the nonsubsampled pyramid in the NSCT are calculated via maximizing the polarimetric contrast between the adjacent subband levels, instead of using the difference of the adjacent subbands as used in the additive noise model. By this way, we make the NSCT applicable to PolSAR data and multiband data. Then, the edges are detected in the NSCT domain based on a fusion of the directional subband coefficients at different scales. Experimental results with both simulated and real PolSAR data show that the present approach is robust to noise and the extracted edges are complete and continuous.

Description: A method for defining the spatial resolution of a Global Navigation Satellite System reflectometry delay–Doppler map (DDM) and of any derived geophysical product is proposed. An effective spatial resolution is derived as a function of measurement geometry and delay–Doppler (DD) interval, and as a more appropriate representation of resolution than the geometric resolution previously used in the literature. The definition more accurately accounts for variations in the scattered power across different pixels of the DDM and more accurately includes the power spreading effect caused by the Woodward ambiguity function. The dependence of the effective resolution on incidence angle, receiver altitude, and DD interval is analyzed and compared with the dependence of the geometric resolution with similar parameters.

Description: Laser scanner-captured 3-D point cloud data analysis is becoming more commonly used for remote sensing and plant science applications. Because of nonrigidity and complexity, reconstructing a 3-D model of a plant is extremely challenging. Existing algorithms often fail to find correct correspondences for plantlike thin structures. We address the problem of finding 3-D junction points in plant point cloud data as a first step of this correspondence matching process. Temporarily, we transform the 3-D problem into 2-D by performing appropriate coordinate transformations to the neighborhood of each 3-D point. Our proposed method has two steps. First, a statistical dip test of multimodality is performed to detect the nonlinearity of the local 2D structure. Then, each branch is approximated by sequential random-sample-consensus line fitting and a Euclidean clustering technique. The straight line parameters of each branch are extracted using total-least-squares estimation. Finally, the straight line equations are solved to determine if they intersect in the local neighborhood. Such junction points are good candidates for subsequent correspondence algorithms. Using these detected junction points, we formulate a correspondence algorithm as a subgraph matching problem and show that, without using traditional descriptor similarity-based matching, good correspondences can be obtained by simply considering geodesic distances among graph nodes. Experiments on synthetic and real ( Arabidopsis plant) data show that the proposed method outperforms the state of the art.

Description: Repeat-track analysis is commonly utilized to generate elevation change time series from satellite radar altimetry over ice sheets. It requires surface gradient (SG) correction due primarily to orbital drifts and radar-related empirical corrections caused by radar scatters from ice surface and potential subsurface. In this letter, two approaches, namely, the use of a digital elevation model (DEM) and the modified repeat-track analysis, which uses the accumulated Envisat altimetry profiles, are applied to correct the SG over both Greenland ice sheet (GrIS) and Antarctic ice sheet (AIS). By comparing the root mean square (rms) of elevation change time series after SG correction, the percentage of data (rms< 1 m) obtained by using modified repeat-track analysis is found to be 85% and 88% for the GrIS and AIS, respectively, as opposed to 45% and 44% if the DEM method is used. Furthermore, three cases are studied to assess empirical corrections for elevation retrieved from both ice-1 and ice-2 algorithms over the AIS. We conclude that the modified repeat-track analysis is more effective to remove topographic induced error. For the ice-2 algorithm, waveform shape parameters are needed in addition to applying corrections from changes in backscatter coefficients. The trend of elevation changes from the ice-1 algorithm with only backscatter analysis agrees with that from the ice-2 algorithm with corrections from backscatter coefficient changes and waveform shape parameters. This study could provide a potential data processing recipe for generating improved satellite radar altimetry elevation time series over ice sheets.

Description: Extracting ships from complex backgrounds is the bottleneck of ship detection in high-resolution optical satellite images. In this letter, we propose a nearly closed-form ship rotated bounding box space used for ship detection and design a method to generate a small number of highly potential candidates based on this space. We first analyze the possibility of accurately covering all ships by labeling rotated bounding boxes. Moreover, to reduce search space, we construct a nearly closed-form ship rotated bounding box space. Then, by scoring for each latent candidate in the space using a two-cascaded linear model followed by binary linear programming, we select a small number of highly potential candidates. Moreover, we also propose a fast version of our method. Experiments on our data set validate the effectiveness of our method and the efficiency of its fast version, which achieves a close detection rate in near real time.

Description: In this letter, a two-stage method for airport detection on remote sensing images is proposed. In the first stage, a new algorithm composed of several line-based processing steps is used for extraction of candidate airport regions. In the second stage, the scale-invariant feature transformation and Fisher vector coding are used for efficient representation of the airport and nonairport regions and support vector machines employed for classification. In order to evaluate the performance of the proposed method, extensive experiments are conducted on airports around the world with different layouts. The measures used in the evaluation are accuracy, sensitivity, and specificity. The proposed method achieved an accuracy of 94.6%, which was benchmarked with two previous methods to prove its superiority.

Description: High-resolution wide-swath synthetic aperture radar (SAR) systems are very attractive for the observation of dynamic processes on the Earth's surface, but they require the downlink of a huge volume of data. In order to comply with azimuth ambiguity requirements, in fact, a pulse repetition frequency much higher than the required processed Doppler bandwidth is often desirable. The volume of downlinked data, however, can be drastically reduced by performing Doppler filtering and decimation on board. A finite-impulse-response filter with a relatively small number of taps suffices to suppress the additional ambiguous components and to recover the original impulse response. This strategy is of special relevance for staggered SAR systems, which are typically characterized by a high oversampling factor. The proposed data reduction technique is also baseline for Tandem-L, where onboard Doppler filtering, resampling, and decimation will be jointly implemented.

Description: In this letter, we generate anisotropic bicontinuous media with different vertical and horizontal correlation functions. With the computer-generated bicontinuous medium, we then use numerical solutions of Maxwell equations in 3-dimensions (NMM3D) to calculate the anisotropic effective permittivities and the effective propagation constants of V and H polarizations. The copolarization phase difference (CPD) of VV and HH is then derived. The CPDs have recently been applied to the retrieval of snow water equivalent, snow depth, and anisotropy. The NMM3D simulation results are also compared with the results of the strong permittivity fluctuations in the low frequency limit and compared against the Maxwell–Garnett mixing formula.

Description: The IEEE GRS Society is grateful for the support given by the organizations listed and invites applications for Institutional Listings from other firms interested in the field of geoscience and remote sensing.

Description: These instructions give guidelines for preparing papers for this publication. Presents information for authors publishing in this journal.

Description: Circumpolar polynya regions and ice production in the Arctic: Results from MODIS thermal infrared imagery for 2002/2003 to 2014/2015 with a regional focus on the Laptev Sea Andreas Preußer, Günther Heinemann, Sascha Willmes, and Stephan Paul The Cryosphere Discuss., doi:10.5194/tc-2016-133,2016 Manuscript under review for TC (discussion: open, 0 comments) Our study presents spatial and temporal characteristics of 16 prominent polynya regions over the entire Arctic basin. By using an energy balance model, we are able to derive daily thin-ice thickness distributions from satellite and atmospheric reanalysis data. All polynyas combined yield an average ice production of about 1444 km 3 per winter. Interestingly, we find distinct regional differences in calculated trends over the last 13 years. Finally, we set a special focus on the Laptev Sea region.

Description: Statistical indicators of Arctic sea-ice stability – prospects and limitations Sebastian Bathiany, Bregje van der Bolt, Mark S. Williamson, Timothy M. Lenton, Marten Scheffer, Egbert H. van Nes, and Dirk Notz The Cryosphere, 10, 1631-1645, doi:10.5194/tc-10-1631-2016, 2016 We examine if a potential "tipping point" in Arctic sea ice, causing abrupt and irreversible sea-ice loss, could be foreseen with statistical early warning signals. We assess this idea by using several models of different complexity. We find robust and consistent trends in variability that are not specific to the existence of a tipping point. While this makes an early warning impossible, it allows to estimate sea-ice variability from only short observational records or reconstructions.

Description: Signature of Arctic first-year ice melt pond fraction in X-band SAR imagery Ane S. Fors, Dmitry V. Divine, Anthony P. Doulgeris, Angelika H. H. Renner, and Sebastian Gerland The Cryosphere Discuss., doi:10.5194/tc-2016-125,2016 Manuscript under review for TC (discussion: open, 0 comments) This paper investigates the signature of melt ponds in satellite borne synthetic aperture radar (SAR) imagery. The study shows that polarimetric X-band SAR can be used to extract melt pond fraction from drifting first-year sea ice, and discusses the influence of factors like wind speed, SAR incidence angle and instrument noise floor. Melt ponds strongly influence the Arctic sea ice energy budget, and the results imply prospective opportunities for expanded monitoring of melt ponds from space.

Description: Greenland annual accumulation along the EGIG line, 1959–2004, from ASIRAS airborne radar and neutron-probe density measurements Thomas B. Overly, Robert L. Hawley, Veit Helm, Elizabeth M. Morris, and Rohan N. Chaudhary The Cryosphere, 10, 1679-1694, doi:10.5194/tc-10-1679-2016, 2016 We demonstrate that snow accumulation rates across the Greenland Ice Sheet, determined from RADAR layers and modeled snow density profiles, are identical to ground-based measurements of snow accumulation. Three regional climate models underestimate snow accumulation compared to RADAR layer estimates. Using RADAR increases spatial coverage and improves accuracy of snow accumulation estimates. Incorporating our results into climate models may reduce uncertainty of sea-level rise estimates.

Description: Thermal impacts of engineering activities and vegetation layer on permafrost in different alpine ecosystems of the Qinghai–Tibet Plateau, China Qingbai Wu, Zhongqiong Zhang, Siru Gao, and Wei Ma The Cryosphere, 10, 1695-1706, doi:10.5194/tc-10-1695-2016, 2016 Climate warming and engineering activities have various impacts on the thermal regime of permafrost in alpine ecosystems of the Qinghai–Tibet Plateau. Using recent observations of permafrost thermal regimes along the Qinghai–Tibet highway and railway, the changes of such regimes beneath embankments constructed in alpine meadows and steppes are studied. The results show that alpine meadows on the Qinghai–Tibet Plateau can have a controlling role among engineering construction effects on permafrost beneath embankments. As before railway construction, the artificial permafrost table (APT) beneath embankments is not only affected by climate change and engineering activities but is also controlled by alpine ecosystems. However, the change rate of APT is not dependent on ecosystem type, which is predominantly affected by climate change and engineering activities. Instead, the rate is mainly related to cooling effects of railway ballast and heat absorption effects of asphalt pavement. No large difference between alpine and steppe can be identified regarding the variation of soil temperature beneath embankments, but this difference is readily identified in the variation of mean annual soil temperature with depth. The vegetation layer in alpine meadows has an insulation role among engineering activity effects on permafrost beneath embankments, but this insulation gradually disappears because the layer decays and compresses over time. On the whole, this layer is advantageous for alleviating permafrost temperature rise in the short term, but its effect gradually weakens in the long term.

Description: Impact of Icebergs on Net Primary Productivity in the Southern Ocean Shuang-Ye Wu and Shugui Hou The Cryosphere Discuss., doi:10.5194/tc-2016-166,2016 Manuscript under review for TC (discussion: open, 0 comments) The primary productivity in the Southern Ocean (SO) is limited by the amount of iron available for biological activities. Recent studies show that icebergs could be a main source of iron to the SO. Based on remote sensing data, our study shows that iceberg presence has a small, yet statistically significant, positive impact on productivity. As iceberg quantity increases, their positive influence on productivity also increases. This impact could serve as a negative feedback to the climate system.

Description: Possible groundwater dominance in the subglacial hydrology of ice sheet interiors: example at Dome C, East Antarctica Brad T. Gooch, Sasha P. Carter, Omar Ghattas, Duncan A. Young, and Donald D. Blankenship The Cryosphere Discuss., doi:10.5194/tc-2016-141,2016 Manuscript under review for TC (discussion: open, 0 comments) Our work investigates the potential significance of groundwater flow underneath the interior of East Antarctica where the ice doesn't rapidly melt. We attempt to describe the relationship between two hydrologic systems (water under the ice and in the ground) and how they might interact along a flow path between lakes under the ice. We find that groundwater is significant in regional water transport for melt water under the ice in areas of low melting in East Antarctica.

Description: Design of a scanning laser meter for monitoring the spatio-temporal evolution of snow depth and its application in the Alps and in Antarctica Ghislain Picard, Laurent Arnaud, Jean-Michel Panel, and Samuel Morin The Cryosphere, 10, 1495-1511, doi:10.5194/tc-10-1495-2016, 2016 A cost-effective automatic laser scan has been built to measure snow depth spatio-temporal variations. Deployed in the Alps and in Dome C (Antarctica), two devices acquired daily scans covering a surface area of 100–150 m 2 . The precision and long-term stability of the measurements are about 1 cm and the accuracy is better than 5 cm. These high performances are particularly suited at Dome C, where it was possible to reveal that most of the accumulation in the year 2015 stems from a single event.

Description: Arctic sea-ice diffusion from observed and simulated Lagrangian trajectories Pierre Rampal, Sylvain Bouillon, Jon Bergh, and Einar Ólason The Cryosphere, 10, 1513-1527, doi:10.5194/tc-10-1513-2016, 2016 Due to the increasing activity in Arctic, sea-ice–ocean models are now frequently used to produce operational forecasts, for oil spill trajectory modelling and to assist in offshore operations planning. In this study we evaluate the performance of two models with respect to their capability to reproduce observed sea ice diffusion properties by using metrics based on Lagrangian statistics. This paper presents a new and useful evaluation metric for current coupled sea ice–ocean models.

Description: Retrieval of the thickness of undeformed sea ice from simulated C-band compact polarimetric SAR images Xi Zhang, Wolfgang Dierking, Jie Zhang, Junmin Meng, and Haitao Lang The Cryosphere, 10, 1529-1545, doi:10.5194/tc-10-1529-2016, 2016 In this work, we introduced a parameter ("CP ratio") for the retrieval of the thickness of undeformed first-year sea ice that is specifically adapted to compact polarimetric SAR images. Based on a validation using other compact polarimetric SAR images from the Labrador Sea, we found a root mean square error of 8 cm and a maximum correlation coefficient of 0.94 for the retrieval procedure when applying it to level ice between 0.1 m and 0.8 m thick.

Description: High resolution boundary conditions of an old ice target near Dome C, Antarctica Duncan A. Young, Jason L. Roberts, Catherine Ritz, Massimo Frezzotti, Enrica Quartini, Marie G. P. Cavitte, Carly R. Tozer, Daniel Steinhage, Stefano Urbini, Hugh F. J. Corr, Tas Van Ommen, and Donald D. Blankenship The Cryosphere Discuss., doi:10.5194/tc-2016-169,2016 Manuscript under review for TC (discussion: open, 0 comments) Glacial ice can preserve records of Earth's ancient atmosphere (including greenhouse gases) that can be compared to ocean sediment records of ice sheet size to better understand climate feedbacks. To find records of the greenhouse gases found in key periods of climate transition, we need to find sites of unmelted old ice near the base of the ice sheet. We performed a high resolution survey of such a site near Europe's Concordia station in East Antarctica, using a multi instrument aircraft.