Search

Glacier displacement can in principle be measured at the large-scale by cross-correlation of satellite images. At weekly to monthly scales, the expected displacement is often of the same order as the noise for the commonly used satellite images, complicating the retrieval of accurate glacier velocity. Assessments of velocity changes on short time scales and over complex areas such as mountain ranges are therefore still lacking, but are essential to better understand how glacier dynamics are driven by internal and external factors. In this study, we take advantage of the wide availability and redundancy of satellite imagery over the Western Pamir to retrieve glacier velocity changes over 10 days for 7 years for a wide range of glacier geometry and dynamics. Our results reveal strong seasonal trends. In spring/summer, we observe velocity increases of up to 300% compared to a slow winter period. These accelerations clearly migrate upglacier throughout the melt-season, which we link to changes in subglacial hydrology efficiency. In autumn, we observe glacier accelerations that have rarely been observed before. These episodes are primarily confined to the upper ablation zone with a clear downglacier migration. We suggest that they result from glacier instabilities caused by sudden subglacial pressurization in response to (1) supraglacial pond drainage and/or (2) gradual closure of the hydrological system. Our 10-day resolved measurements allow us to characterize the short-term response of glacier to changing meteorological and climatic conditions.

Accurate measurements of ice flow are essential to predict future changes in glaciers and ice caps. Glacier displacement can in principle be measured at the large-scale by cross-correlation of satellite images. At weekly to monthly scales, the expected displacement is often of the same order noise for the commonly used satellite images, which limits the retrieval of accurate glacier velocity. Assessments of velocity changes on short time scales and over complex areas such as mountain ranges are therefore still lacking, but are essential to better understand how glacier dynamics are driven by internal and external factors. In this study, we take advantage of the wide availability and redundancy of satellite imagery over the Western Pamir to retrieve 10-day glacier velocity changes over 7 years for a wide range of glacier geometry and dynamics. Our results reveal strong seasonal trends. In spring/summer, we observe velocity increases of up to 300 % compared to a slow winter period. These accelerations clearly migrate upglacier throughout the melt-season, which we link to changes in subglacial hydrology efficiency. In autumn, we observe glacier accelerations that have rarely been observed before. These episodes are primarily confined to the upper ablation zone with a clear downglacier migration. We suggest that they result from glacier instabilities caused by sudden subglacial pressurization in response to (1) supraglacial pond drainage and/or (2) gradual closure of the hydrological system. Our 10-day resolved measurements allow us to characterize the short-term response of glacier to changing meteorological and climatic conditions.

Small ocean eddies are observed in fine-resolution satellite imagery within the coastal regions and ice margins of the world's oceans. The derivation of surface velocity fields of these features would enable geophysical estimates of the kinematic energy contained within the eddies, leading to improved understanding of the role these eddies play in ocean circulation. In this paper, we develop a method to derive surface velocity fields of small ocean eddies using fine resolution, multi-sensor imagery obtained over short time intervals. We consider the evolving image to be embedded in a deformable medium, and perform image matching using nonlinear viscous fluid registration model. We employ sum of squared differences and mutual information as matching functionals. The physical continuum equation is solved using an efficient multigrid full approximation scheme. We test our method using synthetic aperture radar (SAR) image pairs and show that it produces promising results.

Water flow patterns across coastal and deltaic wetlands affect biogeochemical cycling, denitrification, organic carbon burial, and coastal landscape evolution. Our understanding of such patterns across these important landscapes is incomplete, however, because of the inherent difficulty of conducting spatially and temporally dense ground- or boat-based surveys in shallow, vegetated terrain. We conducted an airborne L-band synthetic aperture radar (SAR) acquisition campaign on Wax Lake Delta, Louisiana, USA, in May 2015, to investigate whether water velocities and flow patterns over kilometer scales can be determined from remote sensing. Thirteen SAR flight lines over the delta region were acquired in 3 h with six different flight directions, concurrently with a small boat campaign. We show that SAR azimuth displacement due to Doppler shift can be used to estimate the surface water flow relative to the static and submerged vegetation interspersed on delta islands, using a simple Bragg wave scattering model and accounting for the Bragg wave's free velocity and wind drift. At Wax Lake Delta, we find that ~0.40 m/s water velocities within the main deltaic channels slow to 0.1–0.2 m/s as flow spreads laterally across, and converges within, the vegetated islands, coinciding with shallow (

During its ascent up Mount Sharp, the Mars Science Laboratory Curiosity rover traversed the Bagnold Dune Field. We model sand modal mineralogy and grain size at four locations near the rover traverse, using orbital shortwave infrared single-scattering albedo spectra and a Markov chain Monte Carlo implementation of Hapke's radiative transfer theory to fully constrain uncertainties and permitted solutions. These predictions, evaluated against in situ measurements at one site from the Curiosity rover, show that X-ray diffraction-measured mineralogy of the basaltic sands is within the 95% confidence interval of model predictions. However, predictions are relatively insensitive to grain size and are nonunique, especially when modeling the composition of minerals with solid solutions. We find an overall basaltic mineralogy and show subtle spatial variations in composition in and around the Bagnold Dunes, consistent with a mafic enrichment of sands with cumulative aeolian-transport distance by sorting of olivine, pyroxene, and plagioclase grains. Furthermore, the large variations in Fe and Mg abundances (~20 wt %) at the Bagnold Dunes suggest that compositional variability may be enhanced by local mixing of well-sorted sand with proximal sand sources. Our estimates demonstrate a method for orbital quantification of composition with rigorous uncertainty determination and provide key constraints for interpreting in situ measurements of compositional variability within Martian aeolian sandstones.

The transport and deposition of sand is the most prevalent agent of landscape modification on Mars today, with fluxes comparable to some sand dunes on Earth. Until now, the relationship between sand flux and dune field morphology has been poorly constrained. By tracking dune movement over ∼10 km-long dune fields in Herschel Crater and Nili Patera, representative of many dune fields on Mars, we find a downwind flux decrease that correlates with a sequence of changing morphology from barchans to barchanoids and seifs (longitudinal dunes) to isolated dome dunes and ending with sand sheets. We show empirical consistency with atmospheric Internal Boundary Layer (IBL) theory which can describe these broad flux and morphology changes in Martian dune fields. Deviations from IBL flux predictions are from wind streamline compressions up slopes, leading to a speedup effect. By establishing a dune field morphology type example and correlating it with measured and predicted flux changes, we provide an integrated morphology and flux model that can be applied to other areas of Mars and be used to infer paleo-environmental conditions from preserved sandstone.

Shallow coastal regions are among the fastest evolving landscapes but are notoriously difficult to measure with high spatiotemporal resolution. Using Uninhabited Aerial Vehicle Synthetic Aperture Radar (UAVSAR) data, we demonstrate that high signal‐to‐noise L band synthetic aperture radar (SAR) can reveal subaqueous channel networks at the distal ends of river deltas. Using 27 UAVSAR images collected between 2009 and 2015 from the Wax Lake Delta in coastal Louisiana, USA, we show that under normal tidal conditions, planform geometry of the distributary channel network is frequently resolved in the UAVSAR images, including ~700 m of seaward network extension over 5 years for one channel. UAVSAR also reveals regions of subaerial and subaqueous vegetation, streaklines of biogenic surfactants, and what appear to be small distributary channels aliased by the survey grid, all illustrating the value of fine resolution, low noise, L band SAR for mapping the nearshore subaqueous delta channel network.

This study explores validating and calibrating the wind regime predicted by Global Circulation Models (GCM) on Earth and other planets using optical remote sensing of dune dynamics. We use Spot-5 images to track the migration of 64 Barchan dunes in Qatar using the COSI-Corr technique. We estimate the volume of the dunes using a scaling law calibrated from one particular dune, which was surveyed in the field. Using volume and migration rate, we determine the sand flux from a single dune, Q_(Dunes), and scale this estimate to the whole dune field. We compare the measured sand flux with those derived from wind velocity measurements at a local meteorological station as well as with those predicted from ERA-Interim (a Global Circulation Model). The comparison revealed that the wind velocity predicted by ERA-Interim is inappropriate to calculate the sand flux. This is due to the 6-h sampling rate and to systematic bias revealed by a comparison with the local wind data. We describe a simple procedure to correct for these effects. With the proposed correction, similar sand flux are predicted using the local and ERA-Interim data, independently of the value of the value of the shear velocity threshold, u_(*t). The predicted sand flux is about 65% of Q_(Dunes). The agreement is best assuming the value u_(*t)=0.244 m/s, which is only slightly larger than the value of u_(*t)=0.2612 m/s estimated based in the sand granulometry measured from field samples. The influence of the dune topography on the wind velocity field could explain the underestimation. In any case, the study demonstrates the possibility of validating GCM model and calibrating aeolian sand transport laws using remote sensing measurements of dune dynamics and highlights the caveats associated to such an approach.

The state of seismic activity on Mars is currently unknown. On Earth, coseismic displacement has been observed using visible wavelength images and subpixel coregistration and correlation techniques. We apply this method to Mars with the COSI‐Corr (Co‐registration of Optically Sensed Images and Correlation) software package using High‐Resolution Imaging Science Experiment (HiRISE) images, focusing on part of the Cerberus Fossae fault system. We derive best practices for applying this method to the study of coseismic displacements on Mars. Using a time series of eight overlapping HiRISE images, we achieve pixel coregistration with a mean accuracy of about 1/50 of a pixel. We see no clear evidence for coseismic displacement in this region during a time period of over 8 years. One possible displacement signal (1–2 m of west‐east displacement over a length scale of about 50 m) that has similarities to terrestrial coseismic deformation was dismissed as the result of incomplete correction of steep topography during the coregistration stage. Ancillary observations of recurring slope lineae (RSL) activity in the surrounding fault system offer no supporting evidence for the occurrence of coseismic displacement but do seem to suggest RSL activity that does not fit into previous seasonal timescales. Although it is unlikely that we have observed coseismic displacement in our study area during this time period, the best practice method and the accuracy of our results offer encouragement for future studies. HiRISE, Context Camera, and Color and Stereo Surface Imaging System images can be used to complement and independently verify in situ seismic observations by the InSight (Interior Exploration using Seismic Investigations, Geodesy, and Heat Transport) lander or source location changes from the ExoMars Trace Gas Orbiter spacecraft.

Coseismic surface deformation in large earthquakes is typically measured using field mapping and with a range of geodetic methods (e.g., InSAR, lidar differencing, and GPS). Current methods, however, either fail to capture patterns of near-field coseismic surface deformation or lack preevent data. Consequently, the characteristics of off-fault deformation and the parameters that control it remain poorly understood. We develop a standardized method to fully measure the surface, near-field, coseismic deformation patterns at high resolution using the COSI-Corr program by correlating pairs of aerial photographs taken before and after the 1992 M_w 7.3 Landers earthquake. COSI-Corr offers the advantage of measuring displacement across the entire zone of surface deformation and over a wider aperture than that available to field geologists. For the Landers earthquake, our measured displacements are systematically larger than the field measurements, indicating the presence of off-fault deformation. We show that 46% of the total surface displacement occurred as off-fault deformation, over a mean deformation width of 154 m. The magnitude and width of off-fault deformation along the rupture is primarily controlled by the macroscopic structural complexity of the fault system, with a weak correlation with the type of near-surface materials through which the rupture propagated. Both the magnitude and width of distributed deformation are largest in stepovers, bends, and at the southern termination of the surface rupture. We find that slip along the surface rupture exhibits a consistent degree of variability at all observable length scales and that the slip distribution is self-affine fractal with dimension of 1.56.

We define a global matching framework based on energy pyramid, the Global Matching via Energy Pyramid (GM-EP) algorithm, which estimates the disparity map from a single stereo-pair by solving an energy minimization problem. We efficiently address this minimization by globally optimizing a coarse to fine sequence of sparse Conditional Random Fields (CRF) directly defined on the energy. This global discrete optimization approach guarantees that at each scale we obtain a near optimal solution, and we demonstrate its superiority over state of the art image pyramid approaches through application to real stereo-pairs. We conclude that multiscale approaches should be build on energy pyramids rather than on image pyramids.

We used aerial photos taken before and after the 21 September 1999, M_w 7.6, Chi-Chi earthquake in central Taiwan to measure the near-field ground deformation. A total of 12 pairs of images were processed with Co-registration of Optically Sensed Images and Correlation to produce a horizontal displacement map of a 10 km × 10 km area near Tsaotun. Using pairs of images with different viewing angles, both the horizontal and vertical slip across the fault zone can be measured. Our measurements when resampled into lower resolution are consistent with lower resolution measurements of horizontal displacements obtained from SPOT images, as well as with vertical displacements obtained from repeated leveling measurements and field observations. Horizontal strain is strongly localized along the Chelungpu fault (CLPF) and along a secondary scarp that runs parallel to the CLPF about 2 km to the east, the Ailiao fold scarp (ALF). This pattern closely matches the surface ruptures mapped in the field. Horizontal strain across CLPF correlates remarkably well with the topographic features produced by long-term deformation. The cumulative horizontal shortening across the CLPF and ALF amounts to 4.9 ± 0.4 and 6.1 ± 0.6 m, respectively, and fault-parallel displacement is 3.4 ± 0.4 m. The pattern of surface strain is consistent with the interpretation of the ALF as a fold scarp formed over an active axial hinge zone. This study shows that, even in this compressional setting, most surface deformation is localized within narrow fault zones or active axial hinges.

Correlation of KH9 spy and SPOT5 satellite images, airphotos, digital elevation model differencing, electronic distance measurement, and leveling survey data is used to constrain the deformation resulting from the 1975–1984 Krafla rifting crisis. We find that diking typically extends to depths of 5 km, while the dike tops range from 0 km in the caldera region to 3 km at the northern end of the rift. Extension is accommodated by diking at depth and normal faulting in the shallowest crust. In the southern section of the Krafla rift, surface opening is 80% of the dike opening at depth. Over the 70–80 km length of the rift, the average dike opening was 4.3–5.4 m. From these estimates, we calculate the total geodetic moment released over the Krafla rift crisis, 4.4–9.0×10^(19) Nm, which is an order of magnitude higher than the seismic moment released over the same time period, ~5.8×10^(18) Nm. The total volume of magma added to the upper crust was 1.1–2.1×10^9m^3. This study highlights how optical image correlation using inexpensive declassified spy satellite and airphotos, combined with simple models of crustal deformation, can provide important constraints on the deformation resulting from past earthquake and volcanic events.

Sub-pixel correlation of pre- and post-event air photos reveal the complete near-field, horizontal surface deformation patterns of the 1992 M_w 7.3 Landers and 1999 M_w 7.1 Hector Mine ruptures. Total surface displacement values for both earthquakes are systematically larger than ‘on-fault’ displacements from geologic field surveys, indicating significant distributed, inelastic deformation occurred along these ruptures. Comparison of these two datasets show 46 ± 10% and 39 ± 22% of the total surface deformation was distributed over fault zones averaging 154 m and 121 m in width for the Landers and Hector Mine events, respectively. Spatial variations of distributed deformation along both ruptures show correlations with the type of near-surface lithology and degree of fault complexity; larger amounts of distributed shear occur where the rupture propagated through loose unconsolidated sediments, and areas of more complex fault structure. These results have basic implications for geologic-geodetic rate comparisons, and probabilistic seismic hazard analysis.

Geodetic slip inversions for three major (M_w > 7) strike-slip earthquakes (1992 Landers, 1999 Hector Mine and 2010 El Mayor–Cucapah) show a 15–60 per cent reduction in slip near the surface (depth < 2 km) relative to the slip at deeper depths (4–6 km). This significant difference between surface coseismic slip and slip at depth has been termed the shallow slip deficit (SSD). The large magnitude of this deficit has been an enigma since it cannot be explained by shallow creep during the interseismic period or by triggered slip from nearby earthquakes. One potential explanation for the SSD is that the previous geodetic inversions lack data coverage close to surface rupture such that the shallow portions of the slip models are poorly resolved and generally underestimated. In this study, we improve the static coseismic slip inversion for these three earthquakes, especially at shallow depths, by: (1) including data capturing the near-fault deformation from optical imagery and SAR azimuth offsets; (2) refining the interferometric synthetic aperture radar processing with non-boxcar phase filtering, model-dependent range corrections, more complete phase unwrapping by SNAPHU (Statistical Non-linear Approach for Phase Unwrapping) assuming a maximum discontinuity and an on-fault correlation mask; (3) using more detailed, geologically constrained fault geometries and (4) incorporating additional campaign global positioning system (GPS) data. The refined slip models result in much smaller SSDs of 3–19 per cent. We suspect that the remaining minor SSD for these earthquakes likely reflects a combination of our elastic model's inability to fully account for near-surface deformation, which will render our estimates of shallow slip minima, and potentially small amounts of interseismic fault creep or triggered slip, which could ‘make up’ a small percentages of the coseismic SSD during the interseismic period. Our results indicate that it is imperative that slip inversions include accurate measurements of near-fault surface deformation to reliably constrain spatial patterns of slip during major strike-slip earthquakes.

We describe and test a procedure to accurately co-register and correlate multi-temporal aerial images. We show that this procedure can be used to measure surface deformation, and explore the performance and limitations of the technique. The algorithms were implemented in a software package, COSI-Corr (available from the Caltech Tectonics Observatory website). The technique is validated on several case examples of co-seismic deformation. First, we measure co-seismic ground deformation due to the 1992, Mw 7.3, Landers, California, earthquake from 1 m resolution aerial photography of the National Aerial Photography Program (United States Geological Survey). The fault ruptures are clearly detected, including small kilometric segments with fault slip as small as a few tens of centimeters. We also obtained similar performance from images of the fault ruptures produced by the 1999 Mw 7.1 Hector Mine, California, earthquake. The measurements are shown to be biased due to the inaccuracy of the Digital Elevation Model, film distortions, scanning artifacts, and ignorance of ground displacements at the location of the tie points used to co-register the multi-temporal images. We show that some of these artifacts can be identified and corrected.

We describe a procedure to accurately measure ground deformations from optical satellite images. Precise orthorectification is obtained owing to an optimized model of the imaging system, where look directions are linearly corrected to compensate for attitude drifts, and sensor orientation uncertainties are accounted for. We introduce a new computation of the inverse projection matrices for which a rigorous resampling is proposed. The irregular resampling problem is explicitly addressed to avoid introducing aliasing in the ortho-rectified images. Image registration and correlation is achieved with a new iterative unbiased processor that estimates the phase plane in the Fourier domain for subpixel shift detection. Without using supplementary data, raw images are wrapped onto the digital elevation model and coregistered with a 1/50 pixel accuracy. The procedure applies to images from any pushbroom imaging system. We analyze its performance using Satellite pour l'Observation de la Terre (SPOT) images in the case of a null test (no coseismic deformation) and in the case of large coseismic deformations due to the Mw 7.1 Hector Mine, California, earthquake of 1999. The proposed technique would also allow precise coregistration of images for the measurement of surface displacements due to ice-flow or geomorphic processes, or for any other change detection applications. A complete software package, the Coregistration of Optically Sensed Images and Correlation, is available for download from the Caltech Tectonics Observatory website.

We analyze the M_w 7.6 Kashmir earthquake of October 8, 2005, using sub-pixel correlation of ASTER images to measure ground deformation, and modeling seismic waveforms. The surface rupture is continuous over a distance of 75 km and cuts across the Hazara syntaxis reactivating the Tanda and the Muzaffarabad faults. North of Muzaffarabad the surface rupture coincides approximately with the MBT, on the southwestern flank of the syntaxis, although the two faults have opposite dip angles. The rupture terminates abruptly at the hairpin turn of the MBT showing a strong structural control. The fault offset is 4 m on average and peaks to 7 m northwest of Muzaffarabad. The rupture lasted about 25 s and propagated updip and bi-laterally ~2 km/s, with a rise time of 2–5 s. The shallowness and compactness of the rupture, both in time and space, provide an explanation for the intensity of destructions. This kind of analysis could be achieved as soon as a post-earthquake image is available, and would provide key information for early assessment of damages. The study sheds some light on seismic hazard in the Himalaya, and raises concern regarding the possibility of a repetition of the 1555 event which presumably ruptured the Himalayan front south of the Kashmir basin and may have reached a magnitude M_w > 8.

Coupling between surface winds and saltation is a fundamental factor governing geological activity and climate on Mars. Saltation of sand is crucial for both erosion of the surface and dust lifting into the atmosphere. Wind tunnel experiments along with measurements from surface meteorology stations and modelling of wind speeds suggest that winds should only rarely move sand on Mars. However, evidence for currently active dune migration has recently accumulated. Crucially, the frequency of sand-moving events and the implied threshold wind stresses for saltation have remained unknown. Here we present detailed measurements of Nili Patera dune field based on High Resolution Imaging Science Experiment images, demonstrating that sand motion occurs daily throughout much of the year and that the resulting sand flux is strongly seasonal. Analysis of the seasonal sand flux variation suggests an effective threshold for sand motion for application to large-scale model wind fields (1-100 km scale) of τ(s)=0.01±0.0015 N m(-2).

International audience; Great earthquakes rarely occur within active accretionary prisms, despite the intense long‐term deformation associated with the formation of these geologic structures. This paucity of earthquakes is often attributed to partitioning of deformation across multiple structures as well as aseismic deformation within and at the base of the prism (Davis et al., 1983). We use teleseismic data and satellite optical and radar imaging of the 2013 Mw 7.7 earthquake that occurred on the southeastern edge of the Makran plate boundary zone to study this unexpected earthquake. We first compute a multiple point‐source solution from W‐phase waveforms to estimate fault geometry and rupture duration and timing. We then derive the distribution of subsurface fault slip from geodetic coseismic offsets. We sample for the slip posterior probability density function using a Bayesian approach, including a full description of the data covariance and accounting for errors in the elastic structure of the crust. The rupture nucleated on a subvertical segment, branching out of the Chaman fault system, and grew into a major earthquake along a 50° north‐dipping thrust fault with significant along‐strike curvature. Fault slip propagated at an average speed of 3.0 km/s for about 180 km and is concentrated in the top 10 km with no displacement on the underlying décollement. This earthquake does not exhibit significant slip deficit near the surface, nor is there significant segmentation of the rupture. We propose that complex interaction between the subduction accommodating the Arabia-Eurasia convergence to the south and the Ornach Nal fault plate boundary between India and Eurasia resulted in the significant strain gradient observed prior to this earthquake. Convergence in this region is accommodated both along the subduction megathrust and as internal deformation of the accretionary wedge.

International audience; We analyse the Mw 7.7 Balochistan earthquake of 09/24/2013 based on ground surface deformation measured from sub-pixel correlation of Landsat-8 images, combined with back-projection and finite source modeling of teleseismic waveforms. The earthquake nucleated south of the Chaman strike-slip fault and propagated southwestward along the Hoshab fault at the front of the Kech Band. The rupture was mostly unilateral, propagated at 3 km/s on average and produced a 200 km surface fault trace with purely strike-slip displacement peaking to 10 m and averaging around 6 m. The finite source model shows that slip was maximum near the surface. Although the Hoshab fault is dipping by 45° to the North, in accordance with its origin as a thrust fault within the Makran accretionary prism, slip was nearly purely strike-slip during that earthquake. Large seismic slip on such a non-optimally oriented fault was enhanced possibly due to the influence of the free surface on dynamic stresses or to particular properties of the fault zone allowing for strong dynamic weakening. Strike-slip faulting on thrust fault within the eastern Makran is interpreted as due to eastward extrusion of the accretionary prism as it bulges out over the Indian plate. Portions of the Makran megathrust, some thrust faults in the Kirthar range and strike-slip faults within the Chaman fault system have been brought closer to failure by this earthquake. Aftershocks cluster within the Chaman fault system north of the epicenter, opposite to the direction of rupture propagation. By contrast, few aftershocks were detected in the area of maximum moment release. In this example, aftershocks cannot be used to infer earthquake characteristics.

We measure the displacement field resulting from the 1975–1984 Krafla rifting crisis, NE Iceland, using optical image correlation. Images are processed using the COSI-Corr software package. Surface extension is accommodated on normal faults and fissures which bound the rift zone, in response to dike injection at depth. Correlation of declassified KH-9 spy and SPOT5 satellite images reveals extension between 1977–2002 (2.5 m average opening over 80 km), while correlation of aerial photos between 1957–1990 provide measurements of the total extension (average 4.3 m opening over 80 km). Our results show ∼8 m of opening immediately north of Krafla caldera, decreasing to 3–4 m at the northern end of the rift. Correlation of aerial photos from 1957–1976 reveal a bi-modal pattern of opening along the rift during the early crisis, which may indicate either two different magma sources located at either end of the rift zone (a similar pattern of opening was observed in the 2005 Afar rift crisis in East Africa), or variations in rock strength along the rift. Our results provide new information on how past dike injection events accommodate long-term plate spreading, as well as providing more details on the Krafla rift crisis. This study also highlights the potential of optical image correlation using inexpensive declassified spy satellite and aerial photos to measure deformation of the Earth's surface going back many decades, thus providing a new tool for measuring Earth surface dynamics, e.g. glaciers, landsliding, coastal erosion, volcano monitoring and earthquake studies, when InSAR and GPS data are not available.

Classic tie-point detection algorithms such as the Scale Invariant Feature Transform (SIFT) show their limitations when the images contain drastic changes or repetitive patterns. This is especially evident when considering multi-temporal series of images for change detection. In order to overcome this limitation we propose a new algorithm, the Affine Parameters Estimation by Random Sampling (APERS), which detects the outliers in a given set of matched points. This is accomplished by estimating the global affine transform defined by the largest subset of points and by detecting the points which are not coherent (outliers) with the transform. Comparisons with state-of-the-art methods such as GroupSAC or ORSA demonstrate the higher performance of the proposed method. In particular, when the proportion of outliers varies between 60% and 90% APERS is able to reject all the outliers while the others fail. Examples with real images and a shaded Digital Elevation Model are provided.

Strong and sustained winds on Mars have been considered rare, on the basis of surface meteorology measurements and global circulation models, raising the question of whether the abundant dunes and evidence for wind erosion seen on the planet are a current process. Recent studies showed sand activity, but could not determine whether entire dunes were moving--implying large sand fluxes--or whether more localized and surficial changes had occurred. Here we present measurements of the migration rate of sand ripples and dune lee fronts at the Nili Patera dune field. We show that the dunes are near steady state, with their entire volumes composed of mobile sand. The dunes have unexpectedly high sand fluxes, similar, for example, to those in Victoria Valley, Antarctica, implying that rates of landscape modification on Mars and Earth are similar.

Applications such as change detection and digital elevation model extraction from optical images require a rigorous modeling of the acquisition geometry. We show that the unrecorded satellite jitter during image acquisition, and the uncertainties on the CCD arrays geometry are the current major limiting factors for applications requiring high accuracy. These artifacts are identified and quantified on several optical satellites, i.e., SPOT, ASTER, QuickBird, and HiRISE.

The increasing availability of high-quality optical satellite images should allow, in principle, continuous monitoring of Earth's surface changes due to geologic processes, climate change, or anthropic activity. For instance, sequential optical images have been used to measure displacements at Earth's surface due to coseismic ground deformation [e.g., Van Puymbroeck et al., 2000], ice flow [Scambos et al., 1992; Berthier et al., 2005], sand dune migration [Crippen, 1992], and landslides [Kääb, 2002; Delacourt et al., 2004]. Surface changes related to agriculture, deforestation, urbanization, and erosion—which do not involve ground displacement—might also be monitored, provided that the images can be registered with sufficient accuracy. Although the approach is simple in principle, its use is still limited, mainly because of geometric distortion of the images induced by the imaging system, biased correlation techniques, and implementation difficulties.

Recent methodological progress, Co-Registration of Optically Sensed Images and Correlation, outlined here, makes it possible to measure horizontal ground deformation from optical images on an operational basis, using the COSI-Corr software package. In particular, its sub-pixel capabilities allow for accurate mapping of surface ruptures and measurement of co-seismic offsets. We retrieved the fault rupture of the 2005 Mw 7.6 Kashmir earthquake from ASTER images, and we also present a dense mapping of the 1992 Mw 7.3 Landers earthquake of California, from the mosaicking of 30 pairs of aerial images.

Accurate measurements of ice flow are essential to predict future changes in glaciers and ice caps. Glacier displacement can in principle be measured on the large scale by cross-correlation of satellite images. At weekly to monthly scales, the expected displacement is often of the same order as the noise for the commonly used satellite images, complicating the retrieval of accurate glacier velocity. Assessments of velocity changes on short timescales and over complex areas such as mountain ranges are therefore still lacking but are essential to better understand how glacier dynamics are driven by internal and external factors. In this study, we take advantage of the wide availability and redundancy of satellite imagery over the western Pamirs to retrieve glacier velocity changes over 10 d intervals for 7 years and for a wide range of glacier geometry and dynamics. Our results reveal strong seasonal trends. In spring/summer, we observe velocity increases of up to 300 % compared to a slow winter period. These accelerations clearly migrate upglacier throughout the melt season, which we link to changes in subglacial hydrology efficiency. In autumn, we observe glacier accelerations that have rarely been observed before. These episodes are primarily confined to the upper ablation zone with a clear downglacier migration. We suggest that they result from glacier instabilities caused by sudden subglacial pressurization in response to (1) supraglacial pond drainage and/or (2) gradual closure of the hydrological system. Our 10 d resolved measurements allow us to characterize the short-term response of glaciers to changing meteorological and climatic conditions.