In-Flight CCD Distortion Calibration for Pushbroom Satellites Based on Subpixel Correlation

We describe a method that allows for accurate inflight calibration of the interior orientation of any pushbroom camera and that in particular solves the problem of modeling the distortions induced by charge coupled device (CCD) misalignments. The distortion induced on the ground by each CCD is measured using subpixel correlation between the orthorectified image to be calibrated and an orthorectified reference image that is assumed distortion free. Distortions are modeled as camera defects, which are assumed constant over time. Our results show that in-flight interior orientation calibration reduces internal camera biases by one order of magnitude. In particular, we fully characterize and model the Satellite Pour l'Observation de la Terre (SPOT) 4-HRV1 sensor, and we conjecture that distortions mostly result from the mechanical strain produced when the satellite was launched rather than from effects of on-orbit thermal variations or aging. The derived calibration models have been integrated to the software package Coregistration of Optically Sensed Images and Correlation (COSI-Corr), freely available from the Caltech Tectonics Observatory website. Such calibration models are particularly useful in reducing biases in digital elevation models (DEMs) generated from stereo matching and in improving the accuracy of change detection algorithms

Co-Registration of Optically Sensed Images and Correlation (COSI-Corr): an Operational Methodology for Ground Deformation Measurements

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

The 2001 M_w 7.6 Bhuj earthquake, low fault friction, and the crustal support of plate driving forces in India

We present a source model for the 2001 M_w 7.6 Bhuj earthquake of northwest India. The slip distribution suggests a high stress drop (~35 MPa) and, together with the depth distribution of aftershocks, that the entire crust is seismogenic. We suggest that the active faults have an effective coefficient of friction of ~0.08, which is sufficient for the seismogenic crust to support the majority of the compressive force transmitted through the Indian lithosphere. This model is consistent with the midcrustal depth of the transition from extension to compression beneath the Ganges foreland basin where India underthrusts southern Tibet. If the coefficient of friction were the more traditional value of 0.6, the lithosphere would be required to support a net force roughly an order of magnitude higher than current estimates in order to match the observed depth of the neutral fiber

Monitoring Earth Surface Changes from Space

This report gives an overview of the activities which have been undertaken as part of the technical follow-on to the large study “Monitoring Earth Surface Changes from Space”. In addition to the support provided by the Keck Institute for Space Studies, these activities have been supported by matching funds from the Gordon and Betty Moore Foundation, from UAE and Kuwait, and from the MDAP NASA program. Activities were organized under five different themes, each lead by a different PI: 1- Optical Image Time-Series (PI: Sebastien Leprince). These activities aim at developing techniques to analyze optical images acquired by different imaging systems and at different times to look at general landscape evolution (evolutions due to tectonic activity, glacier flow, landslides, sand dunes migration, etc.). They also aim at building a framework for large scale processing to look at global changes. 2- SAR Time-Series Analysis (PI: Mark Simons). These activities aim at developing techniques to analyze radar image time series, in particular via interferrometric techniques. These activities involve close interactions with JPL via the ARIA project (PI: Susan Owen). 3- Seismic Waves Imaging (PI: Pablo Ampuero). These activities aim at developing techniques for seismic inversion with dense measurement in time and space, such as measurement that would be provided by a space seismometer. These activities involve close interactions with JPL, which received a matching R&TD funding to investigate the development of a space optical seismometer (PI: David Redding). 4- Sub-surface Imaging (PI: Essam Heggy). These activities involve close interactions at testing the possibility of an Earth orbiting Ground Penetrating Radar (GPR). Within the scope of this project, only airborne applications will be sought after, with study for space applications. 5- Science Applications (PI: Mike Lamb). These activities involve taking advantage of the techniques developed by the other groups. It also drives the technical developments and foresees the external visitor program. We detail below these activities. Each sub-section has software products, publications, and/or conference posters/talks as outcome. All publications and presentations in international meetings are listed again at the end of the report together with a few other publications produced by collaborators who have participated in the KISS study but did not receive funding from us. Regarding the ‘seismic waves imaging’ project, we have explored different designs and mission concepts for a 4 m-class Seismic Imager Geostationnary satellite system. We are currently working on estimating the cost and preparing a draft GSI Mission Whitepaper

Influence of camera distortions on satellite image registration and change detection applications

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

Deformation during the 1975–1984 Krafla rifting crisis, NE Iceland, measured from historical optical imagery

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

Inference by Learning: Speeding-up Graphical Model Optimization via a Coarse-to-Fine Cascade of Pruning Classifiers

We propose a general and versatile framework that significantly speeds-up graphical model optimization while maintaining an excellent solution accuracy. The proposed approach, refereed as Inference by Learning or in short as IbyL, relies on a multi-scale pruning scheme that progressively reduces the solution space by use of a coarse-to-fine cascade of learnt classifiers. We thoroughly experiment with classic computer vision related MRF problems, where our novel framework constantly yields a significant time speed-up (with respect to the most efficient inference methods) and obtains a more accurate solution than directly optimizing the MRF. We make our code available on-line [4]

Revisiting Past Earthquakes and Seismo-Volcanic Crises Using Declassified Optical Satellite Imagery

Recent development of the user-friendly software package “Co-registration of Optically Sensed Images and Correlation” (COSI-Corr), which allows for automatic and precise ortho-rectification, co-registration, and sub-pixel correlation of pushbroom satellite and aerial images, has enabled Earth’s surface dynamics to be accurately monitored using optical imagery [1]. This technique compares two images of the Earth’s surface that were acquired at different times, and estimates any potential pixel shifts between them with an accuracy typically better than 1/10 of the pixel size. Correlation of both satellite and aerial images has been successfully used to identify coseismic ground ruptures and quantify fault offsets during large earthquakes [2]–[4], as well as monitoring sand dune migration, landsliding, ice flow [5] [6], and volcanic activity [7] [8]. In this study, we demonstrate that recently declassified US spy satellite images can be used to measure ground deformation resulting from seismotectonic and volcanic events using optical sub-pixel correlation. KH-9 Hexagon satellite images, with a swath size of 250×125 km, were acquired by the US government between 1971 and 1980, and are available for purchase from the United States Geological Survey (USGS) at small cost ($30 per image). During this period, around 29,000 images were acquired globally [9], providing a comprehensive record of the Earth’s surface at 6–9m resolution

Revisiting Past Earthquakes and Seismo-Volcanic Crises Using Declassified Optical Satellite Imagery

Recent development of the user-friendly software package “Co-registration of Optically Sensed Images and Correlation” (COSI-Corr), which allows for automatic and precise ortho-rectification, co-registration, and sub-pixel correlation of pushbroom satellite and aerial images, has enabled Earth’s surface dynamics to be accurately monitored using optical imagery [1]. This technique compares two images of the Earth’s surface that were acquired at different times, and estimates any potential pixel shifts between them with an accuracy typically better than 1/10 of the pixel size. Correlation of both satellite and aerial images has been successfully used to identify coseismic ground ruptures and quantify fault offsets during large earthquakes [2]–[4], as well as monitoring sand dune migration, landsliding, ice flow [5] [6], and volcanic activity [7] [8]. In this study, we demonstrate that recently declassified US spy satellite images can be used to measure ground deformation resulting from seismotectonic and volcanic events using optical sub-pixel correlation. KH-9 Hexagon satellite images, with a swath size of 250×125 km, were acquired by the US government between 1971 and 1980, and are available for purchase from the United States Geological Survey (USGS) at small cost ($30 per image). During this period, around 29,000 images were acquired globally [9], providing a comprehensive record of the Earth’s surface at 6–9m resolution

Rupture Process of the 1999 M_w 7.1 Duzce Earthquake from Joint Analysis of SPOT, GPS, InSAR, Strong-Motion, and Teleseismic Data: A Supershear Rupture with Variable Rupture Velocity

We analyze the rupture process of the 1999 M_w 7.1 Duzce earthquake using seismological, remote sensing, and geodetic data. Ground deformation measured from the subpixel cross correlation of Satellite Pour l'Observation de la Terre (SPOT) images reveals a 55 km long fault trace and smooth surface-slip distribution peaking at 3.5–4 m. The westernmost segment overlaps for over 10 km with ruptures from the M_w 7.4 Izmit earthquake. The 15 km long easternmost segment, which cuts across mountainous topography, had not been reported previously. We determine a well-constrained source model using a four-segment fault geometry using constraints on surface fault slip and inverting Global Positioning System and Interferometric Synthetic Aperture Radar data along with strong-motion records. Our results show that some variability of the rupture velocity and an eastward supershear velocity are required to fit the strong-motion data. The rise time, up to 6 sec, correlates with cumulative slip, suggesting a sliding velocity of about 1 m/sec. The source model predicts teleseismic waveforms well, although early by 2 sec. This time shift is probably due to the weak beginning of the earthquake that is not observable at teleseismic distances. Strong-motion records are relatively well predicted from a source model derived from the teleseismic data using the fault geometry derived from the satellite images. This study demonstrates the benefit of using accurate fault geometries to determine finite-fault source models