Search

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.

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

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

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.