Identifying lichenometrically datable,glacierized terrains: a case study in the cascade range of western North America

Climatic interpretations of recent glacier fluctuations rely on ice-extent chronologies developed from lichenometric ages of Holocene landforms. However, lichenometry requires time- and resource-consuming field surveys, which limit our understanding of glacier chronologies, especially in remote locations. This study presents a rapid, coarse, a priori approach to predicting new field sites where lichenometry can be applied. Geologic, geographic, climatic, and landcover data were used in spatial and supervised classification analyses to identify areas in the Cascade Range of Washington and northern Oregon with similar environmental conditions to those where lichenometric dating techniques had previously been applied. These results focus the attention of researchers to only 1100 km², or 3%, of the broader Cascade Range study area. Though this study concentrates on the utility of lichenometry for dating recent glacier activity in the Cascade Range, the screening method presented is easily translatable to a variety of geomorphic and environmental applications.     

Monitoring the regional and temporal variability of winter snowfall in the arid Andes using digital NOAA/AVHRR data

Abstract

This paper deals with the spatial distribution and the temporal variability of snowfall in the most arid part of the Andes (18°‐ 28°S) during southern hemisphere winter (May‐September). As the official precipitation data is of poor quality, analyses were carried out by means of digital image processing techniques, using NOAA/AVHRR satellite‐data. Through analysis of 24 different snowfall events from six winters, a previously unknown spatial and temporal precipitation pattern in this remote and unexplored area was revealed. Snowfall is most abundant in the southernmost part of the research area and on the western side of the Andes, indicating the Pacific origin of the snowfall.

Nevertheless, the typical snowfall pattern is modified during different periods of the winter. Three typical time periods could be defined and distinguished from one another. Each of these three periods is characterized by typical weather conditions (cold fronts and “cut‐offs “) leading to a distinct snowfall pattern.

As this study is part of a broader paleoclimatic project, the results will serve as a basis for paleoclimatic reconstruction of past climate. Only by knowing the modern circulation and precipitation patterns is it possible to interpret paleoclimatic signals and archives found in the study area (e.g. paleosol, moraines) correctly.     

Recent literature in cartography and geographic information science

Tissot's Indicatrix and regular grids have been used for assessing map projection accuracies. Despite their broad applicability for accuracy assessment, they have limitations in quantifying resampling errors caused by map projections. This is due to the structural uncertainty with regard to the placement and pattern of grids. It is also difficult to calculate the absolute amount of resampling error in each projection. As an alternative to traditional testing methods, the use of random points was investigated. Specifically, random point generation, resampling with spherical block search algorithms, resampling accuracy with a perfect grid, and resampling accuracy with eight projections were investigated and are discussed here. Eight global referencing methods were tested: the equal-area cylindrical, sinusoidal, Mollweide, Eckert IV, Hammer-Aitoff, interrupted Goode homolosine, integerized sinusoidal projections, and the equal area global gridding with a fixed latitudinal metric distance. The resampling accuracy with a perfect grid is about 75 percent. Results showed the sinusoidal and the integerized sinusoidal projections and equal-area global gridding to achieve the highest accuracies.     

Graph theory for analyzing pair-wise data: application to geophysical model parameters estimated from interferometric synthetic aperture radar data at Okmok volcano,Alaska

Graph theory is useful for analyzing time-dependent model parameters estimated from interferometric synthetic aperture radar (InSAR) data in the temporal domain. Plotting acquisition dates (epochs) as vertices and pair-wise interferometric combinations as edges defines an incidence graph. The edge-vertex incidence matrix and the normalized edge Laplacian matrix are factors in the covariance matrix for the pair-wise data. Using empirical measures of residual scatter in the pair-wise observations, we estimate the relative variance at each epoch by inverting the covariance of the pair-wise data. We evaluate the rank deficiency of the corresponding least-squares problem via the edge-vertex incidence matrix. We implement our method in a MATLAB software package called GraphTreeTA available on GitHub (https://github.com/feigl/gipht). We apply temporal adjustment to the data set described in Lu et al. (Geophys Res Solid Earth 110, 2005) at Okmok volcano, Alaska, which erupted most recently in 1997 and 2008. The data set contains 44 differential volumetric changes and uncertainties estimated from interferograms between 1997 and 2004. Estimates show that approximately half of the magma volume lost during the 1997 eruption was recovered by the summer of 2003. Between June 2002 and September 2003, the estimated rate of volumetric increase is \((6.2 \, \pm \, 0.6) \times 10^6~\mathrm{m}^3/\mathrm{year} \). Our preferred model provides a reasonable fit that is compatible with viscoelastic relaxation in the five years following the 1997 eruption. Although we demonstrate the approach using volumetric rates of change, our formulation in terms of incidence graphs applies to any quantity derived from pair-wise differences, such as range change, range gradient, or atmospheric delay.     

Adaptive filtering of GOCE-derived gravity gradients of the disturbing potential in the context of the space-wise approach

Filtering and signal processing techniques have been widely used in the processing of satellite gravity observations to reduce measurement noise and correlation errors. The parameters and types of filters used depend on the statistical and spectral properties of the signal under investigation. Filtering is usually applied in a non-real-time environment. The present work focuses on the implementation of an adaptive filtering technique to process satellite gravity gradiometry data for gravity field modeling. Adaptive filtering algorithms are commonly used in communication systems, noise and echo cancellation, and biomedical applications. Two independent studies have been performed to introduce adaptive signal processing techniques and test the performance of the least mean-squared (LMS) adaptive algorithm for filtering satellite measurements obtained by the gravity field and steady-state ocean circulation explorer (GOCE) mission. In the first study, a Monte Carlo simulation is performed in order to gain insights about the implementation of the LMS algorithm on data with spectral behavior close to that of real GOCE data. In the second study, the LMS algorithm is implemented on real GOCE data. Experiments are also performed to determine suitable filtering parameters. Only the four accurate components of the full GOCE gravity gradient tensor of the disturbing potential are used. The characteristics of the filtered gravity gradients are examined in the time and spectral domain. The obtained filtered GOCE gravity gradients show an agreement of 63–84 mEötvös (depending on the gravity gradient component), in terms of RMS error, when compared to the gravity gradients derived from the EGM2008 geopotential model. Spectral-domain analysis of the filtered gradients shows that the adaptive filters slightly suppress frequencies in the bandwidth of approximately 10–30 mHz. The limitations of the adaptive LMS algorithm are also discussed. The tested filtering algorithm can be connected to and employed in the first computational steps of the space-wise approach, where a time-wise Wiener filter is applied at the first stage of GOCE gravity gradient filtering. The results of this work can be extended to using other adaptive filtering algorithms, such as the recursive least-squares and recursive least-squares lattice filters.     

Searching for ancient gold mines in Filippoi area,Macedonia, Greece,using Worldview-2 satellite imagery

The study area is located near the town of Filippoi, north of the city of Kavala in northern Greece, known from ancient times for its rich gold mines, situated inside hydrothermal alteration zones (Fe–Mn oxide minerals). A Very High-Resolution (0.5 m pixel size) image of Worldview-2 satellite was digitally enhanced, yielding target areas of potential ore existence and lineaments. Ground-truth that followed digital image processing, revealed abandoned ancient mines, slags and ore occurrences. Also, a number of lineaments delineated on the satellite image were verified as faults.     

基于高光谱遥感数据提取森林结构参数的研究

以位于三峡库区的龙门河森林自然保护区为研究区,综合利用线性光谱混合模型和几何光学模型,基于高光谱遥感数据提取森林结构参数是本文研究的重点。在研究区地面调查数据的基础上,通过高光谱数据和混合光谱分解法,获得反演几何光学模型所需的四分量参数,根据背景光照分量与森林植被冠层各参数间的关系,反演得到森林冠层郁闭度及平均冠幅的定量分布图,并利用37个野外实测样本进行结果验证。     

GPS tomography in the polar cap: comparison with ionosondes and in situ spacecraft data

Tomographic 4D reconstructions of ionospheric anomalies appearing in the high-latitude polar cap region are compared with plasma density measurements by digital ionosonde located near the north magnetic pole at Eureka station and with in situ plasma measurements on-board DMSP spacecraft. The moderate magnetic storm of 14–17 October 2002 is taken as an example of a geomagnetic disturbance which generates large-scale ionospheric plasma anomalies at mid-latitudes and in the polar cap region. Comparison of the GPS tomographic reconstructions over Eureka station with the ionosonde measurements of the F layer peak densities indicates that the GPS tomography correctly predicts the time of arrival and passage of the ionospheric tongue of ionization over the magnetic pole area, although the tomographic technique appears to under-estimate the value of F peak plasma density. Comparison with the in situ plasma measurements by the DMSP SSIES instruments shows that the GPS tomography correctly reproduces the large-scale spatial structure of ionospheric anomalies over a wide range of latitudes from mid-latitudes to the high-latitude polar cap region, though the tomographic reconstructions tend to over-estimate the density of the topside ionosphere at 840 km DMSP orbit. This study is essential for understanding the quality and limitations of the tomographic reconstruction techniques, particularly in high-latitude regions where GPS TEC measurements and other ionospheric data sources are limited.     

Synoptic charts of solar 9.1 cm and coronal hole data

Synoptic charts for Carrington rotations 1601–1605 (May–August, 1973) were prepared using the central meridian column of the daily 9.1 cm Stanford solar radio maps. These charts were especially contoured to emphasize temperatures near the quiet solar disk level. Synoptic charts of coronal holes from the ATM-Skylab were superimposed on the radio data to investigate the ability of the radio charts to show coronal holes. This brief period is unfortunately the only interval for which both sets of data are available. The conclusion reached is that in spite of certain problems due to active regions, side-lobe effects and a rather large beamwidth, the 9.1 cm synoptic charts can be of substantial value in identifying large coronal holes, especially during periods of low solar activity. Such synoptic charts, therefore, for the years 1962–1973 that Stanford data are available, could enhance significantly the meagre data pool for coronal holes prior to the Skylab mission.