Geology Reports⌕ Search

Geology topics

Randolph L. Kirk

Publications and source records attributed to Randolph L. Kirk.

At least 19 recordsLinked to original sources
USGS2026

Geodetic investigations of the Europa Clipper mission

The Europa Clipper mission will investigate the geophysical properties of Europa, one of Jupiter’s Galilean moons, to assess its habitability. Geodetic measurements will play a critical role in determining Europa’s internal structure, rotational state, and tidal deformation. The Geodesy Focus Group (GFG) coordinates cross-instrument efforts to measure Europa’s global shape, rotational parameters, gravity field, and degree-2 tidal Love numbers ( k 2 and h 2 ). The data provided by the Gravity/Radio Science (G/RS) investigation, Europa Imaging System (EIS), Radar for Europa Assessment and Sounding (REASON), and Europa Ultraviolet Spectrometer (UVS) will be used to infer geodetic constraints on the interior structure and construct a precise cartographic reference system. This combined dataset will provide new constraints on Europa’s tidal response, ice shell thickness, and the properties of its subsurface ocean. The resulting geodetic information will contribute to a deeper understanding of Europa’s internal dynamics and the potential habitability of its ocean.

Space Science Reviews↗
USGS2024

The Europa Imaging System (EIS) investigation

The Europa Imaging System (EIS) consists of a Narrow-Angle Camera (NAC) and a Wide-Angle Camera (WAC) that are designed to work together to address high-priority science objectives regarding Europa’s geology, composition, and the nature of its ice shell. EIS accommodates variable geometry and illumination during rapid, low-altitude flybys with both framing and pushbroom imaging capability using rapid-readout, 8-megapixel (4k × 2k) detectors. Color observations are acquired using pushbroom imaging with up to six broadband filters. The data processing units (DPUs) perform digital time delay integration (TDI) to enhance signal-to-noise ratios and use readout strategies to measure and correct spacecraft jitter. The NAC has a 2.3° × 1.2° field of view (FOV) with a 10-μrad instantaneous FOV (IFOV), thus achieving 0.5-m pixel scale over a swath that is 2 km wide and several km long from a range of 50 km. The NAC is mounted on a 2-axis gimbal, ±30° cross- and along-track, that enables independent targeting and near-global (≥90%) mapping of Europa at ≤100-m pixel scale (to date, only ∼15% of Europa has been imaged at ≤900 m/pixel), as well as stereo imaging from as close as 50-km altitude to generate digital terrain models (DTMs) with ≤4-m ground sample distance (GSD) and ≤0.5-m vertical precision. The NAC will also perform observations at long range to search for potential erupting plumes, achieving 10-km pixel scale at a distance of one million kilometers. The WAC has a 48° × 24° FOV with a 218-μrad IFOV, achieving 11-m pixel scale at the center of a 44-km-wide swath from a range of 50 km, and generating DTMs with 32-m GSD and ≤4-m vertical precision. The WAC is designed to acquire three-line pushbroom stereo and color swaths along flyby ground-tracks.

Space Science Reviews↗
USGS2024

Measuring erosional and depositional patterns across Comet 67P's imhotep region

Comet 67P/Churyumov-Gerasimenko displays a pronounced hemispherical dichotomy in surface morphology, where the southern hemisphere exhibits more erosional features than the northern hemisphere due to receiving much greater solar radiation. Consequently, it is generally assumed that particles are ejected from the southern hemisphere through sublimation and a significant fraction eventually descends as airfall, covering the northern terrains. To investigate this south-to-north material transfer during the comet's perihelion passage, we used photoclinometry to measure material redistribution within its most extensive smooth terrain deposit around the Imhotep region. However, our findings do not align with this expected trend. Instead, we show that local-scale processes substantially impact the erosion and accumulation of material, with one area experiencing net erosion while another nearby region, just a few dozen meters away, sees sediment buildup. Our analysis underscores the complex interplay of processes shaping Comet 67P's surface and likely comets more generally.

Journal of Geophysical Research↗
USGS2024

Current status of the community sensor model standard for the generation of planetary digital terrain models

The creation of accurate elevation models (topography) from stereo images are critical for a large variety of geospatial activities, including the production of digital orthomosaics, change detection, landing site analysis, geologic mapping, rover traverse planning, and spectral analysis. The United Stated Geological Survey, Astrogeology Science Center, continues to transition the supported planetary sensor models to the Community Sensor Model (CSM) standard. This paper describes the current state of use for this photogrammetric standard, supported sensor model types, and qualitatively compares derived topography between SOCET SET and SOCET GXP ( ® BAE Systems) using HiRISE stereo images of Mars. Our transition to the CSM standard will ensure an uninterrupted capability to make these valuable products for Mars and many other extraterrestrial planets and moons.

Remote Sensing↗
USGS2023

Exploring the interior of Europa with the Europa Clipper

The Galileo mission to Jupiter revealed that Europa is an ocean world. The Galileo magnetometer experiment in particular provided strong evidence for a salty subsurface ocean beneath the ice shell, likely in contact with the rocky core. Within the ice shell and ocean, a number of tectonic and geodynamic processes may operate today or have operated at some point in the past, including solid ice convection, diapirism, subsumption, and interstitial lake formation. The science objectives of the Europa Clipper mission include the characterization of Europa’s interior; confirmation of the presence of a subsurface ocean; identification of constraints on the depth to this ocean, and on its salinity and thickness; and determination of processes of material exchange between the surface, ice shell, and ocean. Three broad categories of investigation are planned to interrogate different aspects of the subsurface structure and properties of the ice shell and ocean: magnetic induction, subsurface radar sounding, and tidal deformation. These investigations are supplemented by several auxiliary measurements. Alone, each of these investigations will reveal unique information. Together, the synergy between these investigations will expose the secrets of the Europan interior in unprecedented detail, an essential step in evaluating the habitability of this ocean world.

Space Science Reviews↗
USGS2023

The cartographic lunar reference frame

We describe how lunar cartographic coordinates are currently determined internationally, based on recommendations of the International Astronomical Union (IAU) Working Group on Cartographic Coordinates and Rotational Elements. Currently, a mean Earth/polar axis (ME) system-based frame and ephemeris (JPL DE 421 ME) is used for that purpose. Recently, there has been some discussion that a lunar principal axis system-based frame should be used instead. We describe how that would be a major change to how lunar data are currently acquired, and how data and mapping products are created and archived, searched for, accessed, and used, and our concerns that such a change would cause significant problems. Existing lunar data, data products, papers, and databases could need to be converted, or substantial confusion could result due to the ∼1 km difference between systems, perhaps resulting in navigation or safety issues. The publication and dissemination of any proposal(s) to change systems is encouraged, as are comments on a comparatively much smaller update to the current ME frame proposed here.

Conference Paper↗
USGS2022

Comparison of Digital Terrain Models from two photoclinometry methods

We evaluate the horizontal resolution and vertical precision for digital topographic models (DTMs) of the Moon derived from image radiance information, a process known as photoclinometry (PC) or shape-from-shading (SfS). We use the implementations in two available planetary image processing software systems, single image PC in the U.S. Geological Survey Integrated Software for Imagers and Spectrometers (ISIS) system, and multi-image SfS in the Ames Stereo Pipeline (ASP), and test results obtained with and without use of a starting solution from stereo, with single and multiple images, and for varying illumination conditions. To obtain the higher quality reference DTMs against which the products can be evaluated, we derived DTMs by stereoanalysis of Lunar Reconnaissance Orbiter Narrow-Angle Camera (LROC NAC) images at their native pixel spacing of ∼0.5 m, then produced a 16-m/post stereo DTM from images downsampled to 4 m/pixel and refined it with images at 16 m/pixel. When used with a single image, both algorithms improved resolution (by a factor of 1.4 for PC and 2.4 for SfS compared to stereo). An albedo map produced in ISIS by ratioing the image to a simulation based on the stereo DTM was well correlated with one output by SfS. The albedo correction was crucial for PC with ∼60° incidence but not at ∼80°. DTMs produced by PC and SfS without a starting stereo DTM had larger errors but good detail, and could be useful for many applications. In SfS, it was necessary to increase smoothing to get a usable DTM when the weighting on an a priori DTM was reduced. Multi-image SfS including modeling of spatially varying albedo reduced vertical errors by factors of 1.5 or more compared to single-image SfS.

International Archives of the Photogrammetry, Remo↗
USGS2022

Revealing active Mars with HiRISE digital terrain models

Many discoveries of active surface processes on Mars have been made due to the availability of repeat high-resolution images from the High Resolution Imaging Science Experiment (HiRISE) onboard the Mars Reconnaissance Orbiter. HiRISE stereo images are used to make digital terrain models (DTMs) and orthorectified images (orthoimages). HiRISE DTMs and orthoimage time series have been crucial for advancing the study of active processes such as recurring slope lineae, dune migration, gully activity, and polar processes. We describe the process of making HiRISE DTMs, orthoimage time series, DTM mosaics, and the difference of DTMs, specifically using the ISIS/SOCET Set workflow. HiRISE DTMs are produced at a 1 and 2 m ground sample distance, with a corresponding estimated vertical precision of tens of cm and ∼1 m, respectively. To date, more than 6000 stereo pairs have been acquired by HiRISE and, of these, more than 800 DTMs and 2700 orthoimages have been produced and made available to the public via the Planetary Data System. The intended audiences of this paper are producers, as well as users, of HiRISE DTMs and orthoimages. We discuss the factors that determine the effective resolution, as well as the quality, precision, and accuracy of HiRISE DTMs, and provide examples of their use in time series analyses of active surface processes on Mars.

Remote Sensing↗
USGS2021

How well do we know Europa’s topography? An evaluation of the variability in digital terrain models of Europa.

Jupiter’s moon Europa harbors one of the most likely environments for extant extraterrestrial life. Determining whether Europa is truly habitable requires understanding the structure and thickness of its ice shell, including the existence of perched water or brines. Stereo-derived topography from images acquired by NASA Galileo’s Solid State Imager (SSI) of Europa are often used as a constraint on ice shell structure and heat flow, but the uncertainty in such topography has, to date, not been rigorously assessed. To evaluate the current uncertainty in Europa’s topography we generated and compared digital terrain models (DTMs) of Europa from SSI images using both the open-source Ames Stereo Pipeline (ASP) software and the commercial SOCET SET ® software. After first describing the criteria for assessing stereo quality in detail, we qualitatively and quantitatively describe both the horizontal resolution and vertical precision of the DTMs. We find that the horizontal resolution of the SOCET SET ® DTMs is typically 8–11× the root mean square (RMS) pixel scale of the images, whereas the resolution of the ASP DTMs is 9–13× the maximum pixel scale of the images. We calculate the RMS difference between the ASP and SOCET SET ® DTMs as a proxy for the expected vertical precision (EP), which is a function of the matching accuracy and stereo geometry. We consistently find that the matching accuracy is ~0.5 pixels, which is larger than well-established “rules of thumb” that state that the matching accuracy is 0.2–0.3 pixels. The true EP is therefore ~1.7× larger than might otherwise be assumed. In most cases, DTM errors are approximately normally distributed, and errors that are several times the derived EP occur as expected. However, in two DTMs, larger errors (differences) occur and correlate with real topography. These differences primarily result from manual editing of the SOCET SET ® DTMs. The product of the DTM error and the resolution is typically 4–8 pixel 2 if calculated using the RMS image scale for SOCET SET ® DTMs and the maximum images scale for the ASP DTMs, which is consistent with recent work using martian data sets and suggests that the relationship applies more broadly. We evaluate how ASP parameters affect DTM quality and find that using a smaller subpixel refinement kernel results in DTMs with smaller (better) resolution but, in some cases, larger gaps, which are sometimes reduced by increasing the size of the correlation kernel. We conclude that users of ASP should always systematically evaluate the choice of parameters for a given dataset.

Remote Sensing↗
USGS2021

Evaluating stereo digital terrain model quality at Mars Rover Landing Sites with HRSC, CTX, and HiRISE Images

We have used high-resolution digital terrain models (DTMs) of two rover landing sites based on mosaicked images from the High-Resolution Imaging Science Experiment (HiRISE) camera as a reference to evaluate DTMs based on High-Resolution Stereo Camera (HRSC) and Context Camera (CTX) images. The Next-Generation Automatic Terrain Extraction (NGATE) matcher in the SOCET SET and GXP® commercial photogrammetric systems produces DTMs with good (small) horizontal resolution but large vertical error. Somewhat surprisingly, results for NGATE are terrain dependent, with poorer resolution and smaller errors on smoother surfaces. Multiple approaches to smoothing the NGATE DTMs give similar tradeoffs between resolution and error; a 5 × 5 lowpass filter is near optimal in terms of both combined resolution-error performance and local slope estimation. Smoothing with an area-based matcher, the standard processing for U.S. Geological Survey planetary DTMs, yields similar errors to the 5 × 5 filter at slightly worse resolution. DTMs from the HRSC team processing pipeline fall within this same trade space but are less sensitive to terrain roughness. DTMs produced with the Ames Stereo Pipeline also fall in this space at resolutions intermediate between NGATE and the team pipeline. Considered individually, resolution and error each varied by approximately a factor of 2. Matching errors were 0.2–0.5 pixels but most results fell in the 0.2–0.3 pixel range that has been stated as a rule of thumb in multiple prior studies. Horizontal resolutions of 10–20 image pixels were found, consistently greater than the 3–5 pixel spacing generally used for stereo DTM production. Resolution and precision were inversely correlated; their product varied by ≤20% (4–5 pixels squared). Refinement of the stereo DTM by photoclinometry can yield quantitative improvement in resolution (more than a factor of 2), provided that albedo variations over distances smaller than the stereo DTM resolution are not too severe. We offer specific guidance for both producers and users of planetary stereo DTMs, based on our results.

Remote Sensing↗
USGS2021

Further adventures in Mars DTM quality: Smoothing errors, sharpening details

We have used high-precision, high-resolution digital terrain models (DTMs) of the NASA Mars Science Laboratory (MSL) and Mars 2020 rover landing sites based on mosaicked images from the Mars Reconnaissance Orbiter High Resolution Imaging Science Experiment (MRO HiRISE) camera as a reference data set to evaluate DTMs based on Mars Express High Resolution Stereo Camera (MEX HRSC) images. The Next Generation Automatic Terrain Extraction (NGATE) matcher in the SOCET SET/GXP † commercial photogrammetric system produces DTMs with relatively good (small) horizontal resolution but high error, and results are terrain dependent, with poorer resolution and smaller errors on smoother surfaces. Multiple approaches to smoothing the NGATE DTMs give very similar tradeoffs between resolution and error. Smoothing the NGATE DTMs with a single pass of an area-based matcher, which has been the standard approach to generating planetary DTMs at the U.S. Geological Survey (USGS) to date is probably near-optimal in terms of both combined resolution-error performance and local slope estimation, but smoothing with a 5x5 lowpass filter performs as well or better. DTMs from the HRSC team processing pipeline fall within this same trade space but are less sensitive to terrain roughness. DTMs produced with the Ames Stereo Pipeline also fall in this space at resolutions intermediate between NGATE and the team pipeline. Although DTM resolution and error each vary by a factor of two, their product is much more consistent, varying by <20% across multiple image sets and matching algorithms. Refinement of the stereo DTM by photoclinometry can yield significant quantitative improvement in resolution and some improvement in error (improving their product by as much as a factor of two), provided that albedo variations over distances smaller than the stereo DTM resolution are not too severe.

Conference Paper↗
USGS2020

Coordination of planetary coordinate system recommendations by the IAU Working Group on Cartographic Coordinates and Rotational Elements–2020 status and future

Our goal is to request input from the lunar and planetary community regarding issues of planetary coordinate systems and cartography standards. We begin with an overview of the work of the International Astronomical Union Working Group on Cartographic Coordinates and Rotational Elements. We briefly describe the operations and membership of the Working Group, some of the various uses of the recommendations it makes, our most recent (2018) published report and the recommendations therein, and the outlook for our next such report. We then consider several issues and questions regarding the future of the Working Group and regarding planetary cartography and planetary data spatial infrastructure in general. This includes possible near-term projects, how we and others might collect and consider community input and includes some ideas regarding possible outcomes or future work that will need to be addressed by the Working Group or other organizations.

Conference Paper↗
USGS2020

Evaluating stereo DTM quality at Jezero Crater, Mars with HRSC, CTX, and HiRISE images

We have used a high-precision, high-resolution digital terrain model (DTM) of the NASA Mars 2020 rover Perseverance landing site in Jezero crater based on mosaicked images from the Mars Reconnaissance Orbiter High Resolution Imaging Science Experiment (MRO HiRISE) camera as a reference dataset to evaluate DTMs based on Mars Express High Resolution Stereo Camera (MEX HRSC) and MRO Context camera (CTX) images. Results are consistent with our earlier HRSC-HiRISE comparisons at the Mars Science Laboratory (MSL) Curiosity landing site in Gale crater, confirming that those results were not compromised by the small area compared and potential problems with spatial registration. Specifically, height errors are on the order of half a pixel and correspond to an image matching error of 0.2–0.3 pixel but estimates of horizontal resolution are 10–20 pixels. Products from the HRSC team pipeline at DLR are smoother but more precise vertically than those produced by using the commercial stereo package SOCET SET®. The DLR products are also homogenous in quality, whereas the SOCET products are less smoothed and have higher errors in rougher terrain. Despite this weak variation, our results are consistent with a rule of thumb of 0.2–0.3 pixel matching precision based on many prior studies. Horizontal resolution is significantly coarser than the DTM ground sample distance (GSD), which is typically 3–5 pixels.

Conference Paper↗
USGS2020

Community tools for cartographic and photogrammetric processing of Mars Express HRSC images

In this chapter we describe the software we have developed for photogrammetric processing of images from the Mars Express High Resolution Stereo Camera (MEX HRSC) to produce digital topographic models (DTMs) and orthoimages, as well as testing we have performed. HRSC has returned images, including stereo and color coverage of most of Mars at decameter scales. The instrument team has developed an extremely powerful processing pipeline and delivered a large number of high-level data products, but our independent software is nevertheless of interest because it provides a check on the standard products, sheds light on the capabilities of software elements we use for multiple missions besides HRSC, and is publicly available, giving users the opportunity to make products that may not (yet) be released by the team and custom products such as local mosaics. We have tested our software on images of three areas: Candor Chasma and Nanedi Valles (both the subject of past DTM comparisons reported by Heipke et al., 2007) and Gale crater, which was extensively mapped at pixel scales 50 times finer than HRSC before its selection as the landing site of the Curiosity rover. We find the vertical precision and mean deviation from the altimetry data used as a control reference for our DTMs to be comparable to the nadir image pixel size. The horizontal resolution of the DTMs appears to be an order of magnitude coarser than the lower limit of 3–5 image pixels that is commonly stated.

Book chapter↗
USGS2019

Overview of spirit microscopic imager results

This paper provides an overview of Mars Exploration Rover Spirit Microscopic Imager (MI) operations and the calibration, processing, and analysis of MI data. The focus of this overview is on the last five Earth years (2005–2010) of Spirit 's mission in Gusev crater, supplementing the previous overview of the first 450 sols of the Spirit MI investigation. Updates to radiometric calibration using in‐flight data and improvements in high‐level processing are summarized. Released data products are described, and a table of MI observations, including target/feature names and associated data sets, is appended. The MI observed natural and disturbed exposures of rocks and soils as well as magnets and other rover hardware. These hand‐lens‐scale observations have provided key constraints on interpretations of the formation and geologic history of features, rocks, and soils examined by Spirit . MI images complement observations by other Spirit instruments, and together show that impact and volcanic processes have dominated the origin and evolution of the rocks in Gusev crater, with aqueous activity indicated by the presence of silica‐rich rocks and sulfate‐rich soils. The textures of some of the silica‐rich rocks are similar to terrestrial hot spring deposits, and observations of subsurface cemented layers indicate recent aqueous mobilization of sulfates in places. Wind action has recently modified soils and abraded many of the rocks imaged by the MI, as observed at other Mars landing sites.

Journal of Geophysical Research E: Planets↗
USGS2018

Degradation of 100‐m‐scale rocky ejecta craters at the InSight Landing Site on Mars and implications for surface processes and erosion rates in the hesperian and amazonian

Rocky ejecta craters (RECs) at the Interior Exploration using Seismic Investigations, Geodesy and Heat Transport (InSight) landing site on Elysium Planitia, Mars, provide constraints on crater modification and rates for the Hesperian and Amazonian. The RECs are between 10 m and 1.2 km in diameter and exhibit five classes of preservation. Class 1 represents pristine craters with sharp rims and abundant ejected rocks. From Classes 2 to 5, rims become more subdued, craters are infilled, and the ejecta become discontinuously distributed. High‐Resolution Imaging Science Experiment digital elevation models indicate a maximum depth to diameter ratio of ~0.15, which is lower than pristine models for craters of similar size. The low ratio is related to the presence of a loosely consolidated regolith and early‐stage eolian infill. Rim heights have an average height to diameter ratio of ~0.03 for the most pristine class. The size‐frequency distribution of RECs, plotted using cumulative and differential methods, indicates that crater classes within the diameter range of 200 m to 1.2 km are separated by ~100 to 200 Myr. Smaller craters degrade faster, with classes separated by <100 Myr. Rim erosion can be entirely modeled by nonlinear diffusional processes using the calculated timescales and a constant diffusivity of 8 × 10 −7 m 2 /year for craters 200 to 500 m in diameter. Diffusion models only partly capture depth‐related degradation, which requires eolian infill. Depth degradation and rim erosion rates are 10 −2 to 10 −3 m/Myr, respectively. The rates are consistent with relatively slow modification that is typical of the last two epochs of Martian history.

Journal of Geophysical Research - Planets↗
USGS2018

A novel technique for precision geometric correction of jitter distortion for the Europa Imaging System and other rolling shutter cameras

We use simulated images to demonstrate a novel technique for mitigating geometric distortions caused by platform motion (“jitter”) as two-dimensional image sensors are exposed and read out line by line (“rolling shutter”). The results indicate that the Europa Imaging System (EIS) on NASA’s Europa Clipper can likely meet its scientific goals requiring 0.1-pixel precision. We are therefore adapting the software used to demonstrate and test rolling shutter jitter correction to become part of the standard processing pipeline for EIS. The correction method will also apply to other rolling-shutter cameras, provided they have the operational flexibility to read out selected “check lines” at chosen times during the systematic readout of the frame area

The International Archives of the Photogrammetry, ↗