Search USGSSearch

SEARCH · Search USGS

Results for “Journal of African Earth Sciences”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,229 recordsLinked to original sources

Sedimentological and geochemical characterization of lacustrine deposits of the Babouri-Figuil basin, northern Cameroon: Implications for source rocks distribution and petroleum exploration

The West and Central African Rift System (WCARS) refers to the series of Cretaceous rift basins where commercial hydrocarbon accumulations have been discovered. Some of the WCARS frontier basins are currently being investigated to increase our understanding of these basins in light of new commercial discoveries. The present study was performed in the Babouri-Figuil Basin (BFB), which is genetically related to the WCARS and constitutes an area of interest in terms of petroleum prospecting, where the distribution of petroleum source rocks and potential targets for petroleum exploration across the entire basin is poorly understood. For the current study, an integrated facies analysis along with organic and inorganic geochemical techniques were applied to the basin's Cretaceous deposits with the aim of reconstructing the paleodepositional environment, assessing factors that triggered the input of organic matter, and providing a spatial overview of the organic matter accumulation in the basin based on outcrop samples. An alluvial fan-lacustrine-braided river system is inferred from the facies analysis of the stratigraphic sequence consisting of conglomerate, sandstone, siltstone, limestone, marlstone, and claystone. Bulk analysis of organic matter reveals that black shale and massive claystone are the main prospective petroleum source rocks in the basin. Inorganic geochemical analyses reveal the influence of anoxic conditions, moderate to high primary productivity, and low terrigenous inputs in organic matter enrichment. The formations rich in organic matter are predominantly concentrated in the western and eastern parts of the basin which may represent areas with depressions, characterized by high accommodation space. In terms of the regional context of the WCARS rift basins, typical hydrocarbon exploration in the BFB may target basal-conglomerate, sandstone beds situated directly above or/and interbedded with the Lower Cretaceous source rocks, and the Upper Cretaceous sandstone beds. Basement rocks (granite, granodiorite, and gneisses) and oil shale deposits may represent potential unconventional hydrocarbon exploration. The current integrated study provides an insight that should guide future hydrocarbon exploration campaigns in the basin.

Babouri-Figuil basin

Geophysical architecture and geochronology of the Neoarchean Mentor anorthosite intrusive complex, northwestern Minnesota: Largest anorthosite complex of the Superior Province?

The Neoarchean Mentor anorthosite intrusive complex (MAIC) lies within the Wawa subprovince in northwestern Minnesota, in a region where the Wawa, Quetico, and Wabigoon subprovinces are juxtaposed in close proximity. Archean rocks are entirely concealed, and interpretations are developed from aeromagnetic, gravity, and borehole samples. The MAIC includes both anorthosite (dense, weakly magnetized) and oxide-rich gabbro (strongly magnetized) lithologies. Anorthosite is coarse-grained to megacrystic, intensely altered, and locally brecciated. Pervasive epidote alteration enhanced the density of the anorthosite via introduction of higher-density mineral assemblages, explaining why the MAIC produces a significant gravity high. Oxide-rich gabbro forms a border phase of the MAIC and has potential for vanadium, chromium, and titanium mineralization, and produces a strong aeromagnetic high. The MAIC is interpreted to extend over an area of 640 km 2 , making it the largest known anorthosite complex of the Superior Province, as measured by preserved areal extent. Modeling indicates the MAIC extends more than 6 km into the subsurface. A new Pb–Pb zircon age of 2737.2 ± 4.5 Ma is interpretated as the crystallization age of anorthosite within the MAIC, showing that the MAIC formed well before the ca. 2690 Ma Shebandowanian orogeny, and raising new questions about correlations with other parts of the Wawa subprovince. A low-density batholith, here informally called the Fertile batholith, is interpreted to intrude the southern part of the MAIC. A new Pb–Pb zircon age of 2701.1 ± 6 Ma is interpreted as the magmatic age of the Fertile batholith.

Minnesota

New constraints on location and timing of the Great Lakes tectonic zone, central Upper Peninsula, Michigan, USA

The Great Lakes tectonic zone (GLTZ) forms the boundary between the Wawa–Abitibi and Minnesota River Valley subprovinces within the Archean Superior Province. The GLTZ is concealed for all of its 1100 km length, except for a segment in the central Upper Peninsula of Michigan. There, it is exposed as a northwest-striking mylonite zone along a 11 km segment, extending to the onlap of Paleozoic rocks to the east. Farther east, its location has been unknown. Here, we use aeromagnetic and gravity data to develop interpretations of the expression of the GLTZ and to define its extent under cover. Aeromagnetic gradients over the mylonite zone are interpreted to be produced by structurally juxtaposed rocks with varying magnetizations. Gravity data show a regional gradient along the GLTZ, produced by the juxtaposition of a dense greenstone belt on the north against lower-density gneisses and granites on the south. The GLTZ is interpreted to extend ∼55 km under cover to the east. The GLTZ is terminated on the east by the buried eastern arm of the ca. 1100 Ma Midcontinent Rift. An undeformed granitic dike that cuts the mylonitic foliation produces a U–Pb apatite age of 2523 ± 33 Ma, implying no major post-Archean shearing occurred, and is at odds with previous interpretations of major Proterozoic reactivation. A granite intrusion in the Minnesota River Valley subprovince produces a Pb–Pb zircon age of 2606.9 ± 3.6/7.4 Ma. This suggests that magmatism related to the Sacred Heart orogeny, previously known in Minnesota, extended to Michigan.

Michigan

Rare earth element-mineralized carbonatite in the Bear Lodge Alkaline Complex, USA—Ore genesis implications from fluid inclusion characterization

Rare earth element (REE) resources of the Bear Lodge Alkaline Complex, Wyoming, are hosted in variably leached carbonatite spatially related to diatreme breccia pipes. We investigated the genesis of REE and lesser-known gold resources through fluid inclusion analysis of carbonatite, fluorite breccia, and smoky quartz vein samples. Physicochemical characteristics of inclusion-trapped fluids were evaluated using petrography, microthermometry, Raman spectroscopy, decrepitate mound analysis, energy-dispersive spectroscopy, laser ablation inductively coupled plasma mass spectrometry, and noble gas isotope analysis. Microthermometry results reveal three fluid types that affected carbonatite dikes within deeper zones that escaped near-surface, ore-grade REE enrichment: (1) high-temperature (330–432°C) magmatic fluid captured in fine-grained calcite; (2) REE-enriched alkali bicarbonate-sulfate brine; and (3) low-temperature (117–182°C) diluted magmatic or meteoric water. Multiphase brine-like inclusions contain burbankite, nahcolite, strontianite, celestine and alkali sulfate daughter crystals, linking them to early burbankite mineralization. Peripheral smoky quartz and fluorite occurrences at Smith Ridge, 1.5 km from the central carbonatite dike swarm, contain primary inclusions that are Cl-poor and rich in Na-HCO 3 -SO 4 , similar to secondary and pseudosecondary inclusions in carbonatites. Helium isotopes reveal a MORB-like source for carbonatite samples and an older crust signature at Smith Ridge, consistent with the proximal ridge-top exposures of Archean granite. Results from this fluid inclusion study coupled with previous studies of carbonatite mineral paragenesis, show that light REEs (LREEs) were not mobilized great distances. Instead, burbankite crystallized within carbonatite from alkali bicarbonate fluids. With sodium retained in early burbankite, outward-emanating fluids enriched in potassium relative to sodium (higher K:Na) resulted in potassium–ferric iron metasomatism of silicate host rocks. This alkali fractionation was accompanied by fractionation of LREEs and heavy REEs (HREEs), with LREEs dominating the central carbonatite resources. In contrast, areas of peripheral REE mineralization at Bear Lodge are commonly characterized by higher HREE:LREE ratios. The K:Na ratio of associated fenites or alteration assemblages could be indicative of early crystallized burbankite in carbonatites and REE fractionation processes potentially leading to areas of concentrated HREEs with greater supply vulnerabilities.

South Dakota, Wyoming

Rare earth element potential in coal and coal ash in the U.S. Gulf Coast

United States heavy reliance on imports of critical minerals (CMs), including rare earth elements (REEs), underscores the importance of development of domestic sources. The study objective was to quantify CM and REE concentrations in coal and coal ash in the US Gulf Coast region. CM and REE concentrations were measured for 118 samples from outcrops and 14 mines in the Gulf Coast. Results show that total REE + Yttrium (REY) concentrations (dry coal basis) are comparable to those of the upper continental crust (UCC) with localized hot spots, such as the Texas Gibbons Creek mine (REY ≤ ~ 2860 ppm). When normalized to UCC REY concentration (169 ppm, dry coal basis), REY to UCC ratios for Gulf Coast coal samples range from 0.1 to 17 (median ratio 0.6). REE extractability from lignites is high (median: 63%–93%) using environmentally benign weak acid. In addition to raw coal, coal ash from power plants could also serve as an REE source with a median ratio of REY in ash relative to coal of 4; however, extractability from coal ash is generally much lower (≤ 5% using the same weak acid as in coal). The median basket price for extracted REY as oxides from coal, assuming 70% extractability, is $3.2 per tonne of coal and $186 billion based on 58 billion metric tonnes of dry coal in the Gulf Coast. REEs important for magnets (Pr + Nd + Tb + Dy) account for ~ 80% of the total value. The corresponding median basket price for extracted REY as oxides from coal ash, assuming ~ 30% extractability, is ~$4.4 per tonne of ash and $1.2 billion based on 258 million tonnes of ash. REE production from coal would likely require co-products, such as activated carbon or humic acids, to attain economic viability. Production of REEs from coal ash could offset remediation costs related to potential water contamination. This reconnaissance study shows the potential for REE production from coal and coal ash in the Gulf Coast; however, carbon coproducts and/or societal benefits would likely be required for socioeconomic viability.

Alabama, Florida, Louisiana, Mississippi, Texas

VIPER site analysis

We needed to evaluate available orbital data of NASA’s Volatiles Investigating Polar Exploration Rover (VIPER) mission area in order to derive a variety of maps to help the science team identify scientifically interesting places for the rover to visit and to provide scientific context for our mission. Some of these maps also fulfilled engineering and mission design needs to enable safe and efficient landing and roving. We incorporated data from the Lunar Reconnaissance Orbiter Camera, the Lunar Orbital Laser Altimeter, the Mini-RF instrument, the Chandrayaan-2 Orbital High Resolution Camera, the Korean Pathfinder Lunar Orbiter’s Shadowcam, the Kaguya Spectral Profiler and Multiband Imager, and the Chandrayaan-1 Moon Mineralogy Mapper. We used a variety of techniques to build these maps, including stereogrammetry, shape-from-shading, ice stability depth and surface temperature calculations, and the horizon method for solar illumination and direct-to-Earth communications maps. Altogether, these maps allowed us to survey for boulders, evaluate features in permanently shadowed regions that VIPER might explore, provide mineralogic context for what VIPER’s instruments may learn, estimate the ages and radar properties of craters in the VIPER mission area, and evaluate the potential for gravity traverses with the rover. These data and techniques provided a rich set of information from which both the VIPER science team and engineering teams were able to draw in order to plan a safe landing and to plan a VIPER surface mission that will be both scientifically valuable and robust from an operational perspective.

The Planetary Science Journal

Review: The hydrogeology of critical mineral resources relevant to the energy transition

Attaining the goals of the international treaty on climate change (the Paris Agreement) will greatly increase the demand for the critical minerals required to implement clean-energy technologies. This poses both challenges and opportunities to the hydrogeologic community from several perspectives. Here, important insights that the hydrogeological sciences have to offer for mineral exploration, mineral production, and addressing environmental issues related to mining and mine decommissioning are summarized. This study focuses on copper, cobalt, lithium, and rare earths to represent the broad spectrum of critical minerals and illustrate their relevance by referring to projected demands and production rates. The current understanding of the hydrogeologic processes that form major deposits of these minerals are then summarized. Ore is defined as the naturally occurring material from which minerals of economic value can be extracted, where most ore deposits are the products of complex hydrogeologic couplings between fluid flow, heat transport, solute transport, chemical reactions, and mechanical deformation. Exploration models for the discovery of deeper, hidden deposits are potentially informed by hydrogeologic theory and hydrogeochemical processes. Hydrogeologic understanding and methods are also essential to production and recovery. Longstanding challenges are mine dewatering and (conversely) mine water supply, as well as mineral-extraction practices such as spoil heap leaching and in situ mining. New challenges arise from element extraction from subsurface brines. Finally, the quantity of water use and potential environmental impacts of mining on water quality are at the core of ‘social license’: the approval and acceptance of society to mining activities.

Hydrogeology Journal

Remote compositional analyses of space-weathered lunar maria

Visible-to-shortwave infrared (VSWIR) reflectance spectroscopy has revolutionized our understanding of planetary surface compositions. However, space-weathering processes on airless bodies complicate quantitative compositional analyses. Here, we present a framework to isolate the signatures of space weathering in VSWIR spectra of lunar maria by leveraging radiative transfer modeling under the assumptions that (i) a space-weathered target can be expressed as a mixture of fresh and fully space-weathered components and (ii) remaining signatures can be modeled by including agglutinates as an end-member component. We first validate this approach against laboratory spectra of space-weathered Apollo mare soils of known mineral compositions using a probabilistic Markov Chain Monte Carlo implementation of the Hapke radiative transfer model. Second, we illustrate how this approach can be applied to orbital Moon Mineralogy Mapper data. The proposed space-weathering correction workflow for lunar maria could be expanded to other lunar lithologies and applied to existing and future data sets.

Planetary Science Journal

Bayesian mapping of regionally grouped, sparse, univariate earth science data

Some earth science data are naturally grouped by region, and it is often desirable to map these data by region. However, if there are only a few samples within each region, then the map should be smoothed in an appropriate way to mitigate the problems that arise from having only a few samples. A smoothing algorithm based on a Bayesian hierarchical model is developed and presented in this report. This algorithm has several features that make it especially suitable for mapping earth science data: it can account for measurements that are censored, it can process multiple datasets with different measurement errors and different censoring thresholds, and it can calculate the uncertainty in any statistic that is mapped. The algorithm is demonstrated by mapping gold concentrations that are measured in streambed sediments in the Taylor Mountains quadrangle in southwestern Alaska.

Alaska

Prospectivity modeling of the NASA VIPER landing site at Mons Mouton near the Lunar South Pole

We use a high-resolution digital elevation model and a numerical thermal model to produce a variety of inputs for a water-ice prospectivity model for the Volatiles Investigating Polar Exploration Rover (VIPER) landing site. These input data are maps of topography, surface slope, surface aspect, surface curvature, maximum temperature, depth to ice stability, permanently shadowed regions (PSRs), distance to PSRs, and PSR density. This model predicts where water ice is most likely within the top meter of regolith, assuming plausible relationships between ice concentration and the various inputs. The model is designed to be adjusted in near-real time as data are collected during the VIPER mission. As such, it is a tool for both analyzing data from the mission as well as planning operations. Since the current model, at this point, relies only on orbital remote sensing, the final version will also be a tool to extrapolate the VIPER mission results across the lunar poles.

Planetary Science Journal

Keeping pace: Migrating sea lampreys adjust swim speed in response to shifting water velocity to maintain a constant ground speed

Upstream migrating fish must balance energetic costs with environmental and ecological constraints. Theoretical models predict that fish can conserve energy by adjusting swim speed in response to changing water velocity, but detailed field tests of these predictions remain rare. We used fine-scale acoustic telemetry and hydrodynamic modeling to track adult sea lamprey (Petromyzon marinus) migrating through a river and evaluated how movement speeds varied with environmental conditions and individual traits. Swim speed increased with faster water velocities, while ground speed remained constant, supporting the hypothesis that lampreys regulate swim effort to maintain a fixed migration pace. The average ground speed aligned with a predicted optimum of ∼1 body length per second. Both swim and ground speeds declined slightly in deeper water, suggesting a trade-off between movement efficiency and predator avoidance. These findings provide field-based support for a theoretical model of energy optimization and underscore how migratory fish adaptively adjust behavior to local hydrodynamic conditions. Our results highlight the value of using high-resolution tracking and environmental modeling to reveal behavioral tactics that underlie efficient migration strategies.

Michigan

Identifying overlap between native fish movements and a possible seasonal Grass Carp deterrent in the Sandusky River, Ohio, USA

Riverine fishes are vulnerable to habitat fragmentation caused by interrupted access to vital habitats. Fragmentation using artificial structures can be used to limit the spread of invasive fishes by inhibiting movement between critical habitats. Thus, tradeoffs exist between maintaining connectivity for native species and restricting movement of invasive fishes. A nonphysical deterrent has been proposed to prevent invasive Grass Carp ( Ctenopharyngodon idella ) from reaching spawning habitats in the Sandusky River, a tributary to Lake Erie, during late spring–summer. However, this timing overlaps with spawning periods of a suite of native fishes. We used acoustic telemetry and ichthyoplankton surveys to examine movement and spawning activity of multiple native species and to assess potential impacts of deterrent operation on reproduction. All species used the proposed deterrent location during the Grass Carp spawning period, and their movements corresponded with larval fish collections, indicating upstream spawning. These findings suggest the deterrent could unintentionally restrict native species’ access to essential habitats. Balancing habitat connectivity and invasive species control will likely depend on species-specific responses and careful timing of deterrent operation.

Ohio

Subsurface water ice mapping on Mars: A probabilistic approach

Subsurface water ice deposits on Mars are an important resource for potential future human exploration. They are also an indicator of the planet’s past climate. However, the distribution of subsurface water ice in Mars’s midlatitudes is uncertain because spacecraft imagery cannot directly observe subsurface ice in most cases. Various spacecraft remote sensing instruments are sensitive to subsurface water ice, including thermal imaging spectrometers, radar sounders, and neutron spectrometers. Geomorphic analyses of images can also implicate subsurface ice. Building upon the data products from the Mars Subsurface Water Ice Mapping project, we provide a probabilistic framework to jointly interpret existing data and estimate the likelihood of subsurface water ice in the Martian midlatitudes between 60 ∘ S and 60 ∘ N with uncertainty. Broadly, we find that near-surface ice is likely present poleward of ∼45 ∘ in both the northern and southern hemispheres. However, closer to the equator, existing remote sensing data cannot uniquely constrain the presence of subsurface water ice. Our probabilistic results provide a framework for quantifying the abundance of ice on Mars, and our uncertainty estimates allow future analysis and exploration to target regions of high uncertainty.

Planetary Science Journal

Integrating marine historical ecology into management of Alaska’s Pacific cod fishery for climate readiness

The Pacific cod ( Gadus macrocephalus ) fishery was closed in 2020 after a rapid decline in biomass caused by the marine heat waves of 2014–2019. Pacific cod are exceptionally thermally sensitive and management of this fishery is now challenged by increasingly unpredictable climate conditions. Fisheries monitoring is critical for climate readiness, but short-term monitoring data may be inadequate for recognizing and anticipating change under rapid climate changes. We propose an interdisciplinary, marine historical ecology framework that looks to long-term records (local and traditional knowledge, history, archaeology, and paleoclimatology) to capture a long range of ecological variability and provide historical context for management. In order to connect to contemporary fisheries management, this framework must be built on a common vocabulary and an understanding of the key metrics used in fisheries stock assessments. Here, we propose metrics derived from Pacific cod stock assessment and synthesize information relevant to understanding the effects of past warming periods on cod populations across the Gulf of Alaska and Bering Sea. This case study provides a framework for thinking about how to use these historical records in the context of fisheries management under rapidly changing climate conditions.

Alaska

Cascading land surface hazards as a nexus in the Earth system

Earth’s surface is sculpted by numerous processes that move sediment, ranging from gradual and benign to abrupt and catastrophic. Although infrequent, high-magnitude sediment mobilization events can be hazardous to people and infrastructure, leaving topographic imprints on the landscape and remarkable narratives in the historical record. Hazardous events such as fires, storms, and earthquakes accelerate erosion and sediment transport, increasing landscape sensitivity to subsequent perturbations, thus forming a cascading hazard. Although the redistribution of sediment across Earth’s landscape can result in higher risks to vulnerable populations, cascading processes are commonly unaccounted for in hazard assessments. Cascading hazards can occur almost immediately after triggering events, such as coseismic landslides, or over months, years, or even decades after an initial perturbation, such as debris flows after wildfires or flooding in channels alluviated by volcanic debris. Sediment cascades span Earth’s surface, from mountaintops to river valleys, where erosion, deposition, and aggradation can lead to a myriad of hazardous processes, including decreased river conveyance capacity, which increases the likelihood of downstream flooding. An improved understanding of the magnitude, frequency, and persistence of cascading hazards is critical given the rapid changes in the frequency and severity of storms, fires, sea-level change, and cryospheric melting, as well as the expansion of high-population-density urban footprints in regions susceptible to solid Earth hazards. Understanding the full consequences and underlying physics of Earth’s cascading land surface hazards can help minimize future human and economic losses.

California

Evaluating the effects of a multi-modal deterrent on upstream passage by invasive carp and native species

Invasive carp, including Silver Carp ( Hypophthalmichthys molitrix ) and Grass Carp ( Ctenopharyngodon idella ), are increasing their range within the Mississippi River Basin, threatening native fishes and ecosystem function. We conducted a 3-year field test of a multi-modal deterrent that uses sound, bubbles, and light to test its ability to impede upstream passage by invasive carp. The test was performed at a lock and dam on the Cumberland River, Kentucky, USA, using acoustic telemetry to track the movements of Silver Carp, Grass Carp, and several native fish species. When the deterrent was operating, the estimated probability of upstream lock passage was 53% lower (95% CI: 41%–63%) for Silver Carp but 74% higher (95% CI: 3%–190%) for Grass Carp compared to times when the deterrent was not operating. Paddlefish ( Polyodon spathula ) passage was reduced by 43% (95% CI: 71% reduction to 9% increase), although this was not statistically significant. We detected no significant effects on Freshwater Drum ( Aplodinotus grunniens ) or Smallmouth Buffalo ( Ictiobus bubalus ). However, sample sizes were small for these native species, limiting our ability to identify effects. Water temperature, tailwater elevation, and vessel lockages affected passage rates, although this varied among species.

Kentucky

3D viscoelastic models of slip-deficit rate along the Cascadia subduction zone

Interseismic deformation in the Pacific Northwest is constrained by the horizontal crustal velocity field derived from the Global Positioning System (GPS) in addition to vertical rates derived from GPS, leveling, and tide gauge measurements. Such measurements were folded in to deformation models of fault slip rates as part of the 2023 National Seismic Hazard Model (NSHM) update. Here I build upon one of the contributing models, the viscoelastic earthquake-cycle model of Pollitz [2022]. This model permits inclusion of effects of time-dependent viscoelastic relaxation within earthquake cycles (i.e., ‘ghost transients’) and laterally variable elastic and/or ductile material properties. I lever-age these capabilities to incorporate the Cascadia megathrust into Western U.S.-wide deformation models in which crustal fault slip rates are estimated simultaneously with slip deficit rates along the interplate boundary between the descending Juan de Fuca plate and North American plate. This effort includes construction of a margin-wide model of viscoelastic structure founded on the Slab 2.0 model and probes different models of the ductile properties of the surrounding oceanic asthenosphere, continental lower crust, and mantle asthenosphere. This results in new estimates of the distribution of slip deficit rate along the ∼ 1000 km long margin, highlights the importance of correcting for glacial-isostatic adjustment effects, and permits assessment of sensitivity of results to assumed ductile properties.

California, Oregon, Washington

Thermal detectability of subsurface water ice on Mars: A comparative analysis for the Subsurface Water Ice Mapping (SWIM) Project

We have developed a new global map of shallow ground ice distribution, SWIM23, based on Mars Global Surveyor Thermal Emission Spectrometer data and made systematic comparisons between this new map and two similar, previously developed data products. We have explored the origin of differences between the three ice maps by detailing technical and procedural differences in their development, by making global pixel-by-pixel comparisons, and by carrying out a series of one-dimensional thermal simulations to explore fundamental physical limitations of thermal ice-detection techniques. These efforts and the production of a composite thermal ice-consistency map supported integration of multiple geophysical data products relevant to ice detection in the upper meter of the Martian regolith into the larger Mars Subsurface Water Ice Mapping project. Our work also highlights fundamental physical limitations to thermal ice detection as a technique, particularly the rapid fall-off in ice detection sensitivity at depths >30 cm, which produces maximum uncertainty in the presence and depth of ice within regions preferred for potential human landing sites. A future Mars orbiter mission designed to detect ice and support crewed landing site selection in the midlatitude region should give payload priority to an instrument capable of probing the 1–5 m depth range (i.e., a high-frequency radar), over a next-generation thermal spectrometer, which is unlikely to offer clarity on ice table depths or lateral continuity of the ice table in the locations of highest interest.

Planetary Science Journal