Automated extraction of hydrographically corrected contours for the conterminous United States: the US Geological Survey US Topo product

The US Topo is a new generation of digital topographic maps delivered by the US Geological Survey (USGS). These maps include contours in the traditional 7.5-min quadrangle format. The process for producing digital elevation contours has evolved over several years, but automated production of contours for the US Topo product began in 2010. This process, which is quite complex yet fairly elegant, automatically chooses the proper USGS quadrangle, captures the corresponding 1/3 as grid points from the national elevation data set (3D Elevation Program), and adjusts elevation data to better fit water features from the National Hydrography Dataset. After additional processing, such as identifying and tagging depressions, constructing proper contours across double-line streams, and omitting contours from water bodies, contours are automatically produced for the quadrangle. The resulting contours are then compared to the contours on the original (legacy) topographic map sheets, or to the 10-m contours from the original map sheets. Where the elevation surface used to generate the contours has been derived from the previously published contours for a quadrangle, the generated contours match the legacy contours quite well. Where newer elevation sources, such as lidar, originate the elevation surface, generated contours may vary significantly from the previous cartographically produced contours due to more accurate representations of the surface, and less reliance on cartographic interpretation.     

Crystallization and eruption ages of Breccia Museo (Campi Flegrei caldera,Italy) plutonic clasts and their relation to the Campanian ignimbrite

The Campi Flegrei volcanic district (Naples region, Italy) is a 12-km-wide, restless caldera system that has erupted at least six voluminous ignimbrites during the late Pleistocene, including the >300 km³ Campanian ignimbrite (CI) which originated from the largest known volcanic event of the Mediterranean region. The Breccia Museo (BM), a petrologically heterogeneous and stratigraphically complex volcanic deposit extending over 200 km² in close proximity to Campi Flegrei, has long remained contentious regarding its age and stratigraphic relation to the CI. Here, we present crystallization and eruption ages for BM plutonic ejecta clasts that were determined via uranium decay series and (U–Th)/He dating of zircon, respectively. Despite mineralogical and textural heterogeneity of these syenitic clasts, their U–Th zircon rim crystallization ages are indistinguishable with an average age of 49.7 ± 2.5 ka (2σ errors; mean square of weighted deviates MSWD = 1.2; n = 34). A subset of these crystals was used to obtain disequilibrium-corrected (U–Th)/He zircon ages which average 41.7 ± 1.8 ka (probability of fit P = 0.54; n = 15). This age closely overlaps with published CI ⁴⁰Ar/³⁹Ar eruption ages (40.6 ± 0.1 ka) after recalibration to recently revised flux monitor ages. Concordant eruption ages for BM and CI agree with previous chemostratigraphic and paleomagnetic correlations, suggesting their origin from the same eruption. However, they are at variance with recalibrated ⁴⁰Ar/³⁹Ar ages which have BM postdate CI by 3 ± 1 ka. BM syenites show similar geochemical and Sr–Nd isotopical features of pre-caldera rocks erupted between 58 and 46 ka, but are distinctive from subsequent caldera-forming magmas. Energy-constrained assimilation and fractional crystallization modeling of Nd–Sr isotopic data suggests that pre-caldera magmas formed a carapace of BM-type intrusions in a mid-crust magma chamber (≥8 km depth) shielding the younger CI magma from contamination by Hercynian basement wall rocks. An ~41–50 ka hiatus in crystallization ages implies rapid solidification of these pre-CI intrusions. This argues against protracted pre-eruptive storage of a large volume of CI magma at shallow crustal levels.     

Southern Ocean biogeochemical impact on the tropical ocean: Stable isotope records from the Pacific for the past 25,000 years

Modeling of oceanic nutrient fields indicates that the Southern Ocean may have a strong impact on the chemistry of the thermocline waters which upwell in the eastern tropical oceans and feed biological productivity there. The subantarctic is a primary source of equatorial undercurrent (EUC) waters. The Southern Ocean to equator connection has been shown through modeling to have a potential influence on atmospheric carbon dioxide content via an increase in the efficiency of the tropical biotic pump (silica leakage hypothesis). On the glacial–interglacial timescale the eastern equatorial Pacific (EEP) has a thermocline (EUC) carbon isotope record which is consistent with that idea and a stronger ice age biologic drawdown in the subantarctic. This carbon isotope record features glacial values more positive than those of the Holocene, which is the reverse of what is seen in the non-equatorial, stratified, ocean. We report planktonic carbon isotope records from the Pacific subantarctic in an effort to trace the unique EEP carbon isotope signature to its source. Our results are compatible with the subantarctic as a source of the tropical carbon isotope signature. Analysis of the glacial to Holocene isotopic pattern in terms of causative process indicates that an increased glacial subantarctic biotic pump accounts best for our observations. This supports the hypothesis of Southern Ocean drive on tropical biological production, and potential impact on the global carbon cycle.     

Crystallization and impact history of a meteoritic sample of early lunar crust(NWA 3163)refined by atom probe geochronology

Granulitic lunar meteorites offer rare insights into the timing and nature of igneous,metamorphic and impact processes in the lunar crust.Accurately dating the different events recorded by these materials is very challenging,however,due to low trace element abundances(e.g.Sm,Nd,Lu,Hf),rare micrometerscale U-Th-bearing accessory minerals,and disturbed Ar-Ar systematics following a multi-stage history of shock and thermal metamorphism.Here we report on micro-baddeleyite grains in granulitic mafic breccia NWA 3163 for the first time and show that targeted microstructural analysis(electron backscatter diffraction)and nanoscale geochronology(atom probe tomography)can overcome these barriers to lunar chronology.A twinned(～90°/401)baddeleyite domain yields a ~(232)Th/~(208)Pb age of 4328 ± 309 Ma,which overlaps with a robust secondary ion mass spectrometry(SIMS)~(207)Pb/~(206)Pb age of 4308 ± 18.6 Ma and is interpreted here as the crystallization age for the igneous protolith of NWA 3163.A second microstructural domain,2μm in width,contains patchy overprinting baddeleyite and yields a Th-Pb age of 2175 ±143 Ma,interpreted as dating the last substantial impact event to affect the sample.This finding demonstrates the potential of combining microstructural characterization with nanoscale geochronology when resolving complex P-T-t histories in planetary materials,here yielding the oldest measured crystallization age for components of lunar granulite NWA 3163 and placing further constraints on the formation and evolution of lunar crust.     

Keck spectroscopy of the faint dwarf elliptical galaxy population in the Perseus Cluster core: mixed stellar populations and a flat luminosity function

We present the result of a photometric and Keck low-resolution imaging spectrometer (LRIS) spectroscopic study of dwarf galaxies in the core of the Perseus Cluster, down to a magnitude of   M _B =−12.5  . Spectra were obtained for 23 dwarf-galaxy candidates, from which we measure radial velocities and stellar population characteristics from absorption line indices. From radial velocities obtained using these spectra, we confirm 12 systems as cluster members, with the remaining 11 as non-members. Using these newly confirmed cluster members, we are able to extend the confirmed colour–magnitude relation for the Perseus Cluster down to   M _B =−12.5  . We confirm an increase in the scatter about the colour–magnitude relationship below   M _B =−15.5  , but reject the hypothesis that very red dwarfs are cluster members. We measure the faint-end slope of the luminosity function between   M _B =−18  and −12.5, finding  α=−1.26 ± 0.06  , which is similar to that of the field. This implies that an overabundance of dwarf galaxies does not exist in the core of the Perseus Cluster. By comparing metal and Balmer absorption line indices with α-enhanced single stellar population models, we derive ages and metallicities for these newly confirmed cluster members. We find two distinct dwarf elliptical populations: an old, metal-poor population with ages ∼8 Gyr and metallicities  [Fe/H] < −0.33  , and a young, metal-rich population with ages <5 Gyr and metallicities  [Fe/H] > −0.33  . Dwarf galaxies in the Perseus Cluster are therefore not a simple homogeneous population, but rather exhibit a range in age and metallicity.     

Resistant rock layers amplify cosmogenically-determined erosion rates

Prior numerical modeling work has suggested that incision into sub-horizontal layered stratigraphy with variable erodibility induces non-uniform erosion rates even if base-level fall is steady and sustained. Erosion rates of cliff bands formed in the stronger rocks in a stratigraphic sequence can greatly exceed the rate of base-level fall. Where quartz in downstream sediment is sourced primarily from the stronger, cliff-forming units, erosion rates estimated from concentrations of cosmogenic beryllium-10 (¹⁰Be) in detrital sediment will reflect the locally high erosion rates in retreating cliff bands. We derive theoretical relationships for threshold hillslopes and channels described by the stream-power incision model as a quantitative guide to the potential magnitude of this amplification of ¹⁰Be-derived erosion rates above the rate of base-level fall. Our analyses predict that the degree of erosion rate amplification is a function of bedding dip and either the ratio of rock erodibility in alternating strong and weak layers in the channel network, or the ratio of cliff to intervening-slope gradient on threshold hillslopes. We test our predictions in the cliff-and-bench landscape of the Grand Staircase in southern Utah, USA. We show that detrital cosmogenic erosion rates in this landscape are significantly higher (median 300 m/Ma) than the base-level fall rate (~75 m/Ma) determined from the incision rate of a trunk stream into a ~0.6 Ma basalt flow emplaced along a 16 km reach of the channel. We infer a 3–6-fold range in rock strength from near-surface P-wave velocity measurements. The approximately four-fold difference between the median ¹⁰Be-derived erosion rate and the long-term rate of base-level fall is consistent with our model and the observation that the stronger, cliff-forming lithologies in this landscape are the primary source of quartz in detrital sediments. © 2020 John Wiley & Sons, Ltd.     

Interaction of wind and cold-season hydrologic processes on erosion from complex topography following wildfire in sagebrush steppe

Wildfire is a natural component of sagebrush (Artemisia spp.) steppe rangelands that induces temporal shifts in plant community physiognomy, ground surface conditions, and erosion rates. Fire alteration of the vegetation structure and ground cover in these ecosystems commonly amplifies soil losses by wind- and water-driven erosion. Much of the fire-related erosion research for sagebrush steppe has focused on either erosion by wind over gentle terrain or water-driven erosion under high-intensity rainfall on complex topography. However, many sagebrush rangelands are geographically positioned in snow-dominated uplands with complex terrain in which runoff and sediment delivery occur primarily in winter months associated with cold-season hydrology. Current understanding is limited regarding fire effects on the interaction of wind- and cold-season hydrologic-driven erosion processes for these ecosystems. In this study, we evaluated fire impacts on vegetation, ground cover, soils, and erosion across spatial scales at a snow-dominated mountainous sagebrush site over a 2-year period post-fire. Vegetation, ground cover, and soil conditions were assessed at various plot scales (8 m² to 3.42 ha) through standard field measures. Erosion was quantified through a network of silt fences (n = 24) spanning hillslope and side channel or swale areas, ranging from 0.003 to 3.42 ha in size. Sediment delivery at the watershed scale (129 ha) was assessed by suspended sediment samples of streamflow through a drop-box v-notch weir. Wildfire consumed nearly all above-ground live vegetation at the site and resulted in more than 60% bare ground (bare soil, ash, and rock) in the immediate post-fire period. Widespread wind-driven sediment loading of swales was observed over the first month post-fire and extensive snow drifts were formed in these swales each winter season during the study. In the first year, sediment yields from north- and south-facing aspects averaged 0.99–8.62 t ha⁻¹ at the short-hillslope scale (~0.004 ha), 0.02–1.65 t ha⁻¹ at the long-hillslope scale (0.02–0.46 ha), and 0.24–0.71 t ha⁻¹ at the swale scale (0.65–3.42 ha), and watershed scale sediment yield was 2.47 t ha⁻¹. By the second year post fire, foliar cover exceeded 120% across the site, but bare ground remained more than 60%. Sediment yield in the second year was greatly reduced across short- to long-hillslope scales (0.02–0.04 t ha⁻¹), but was similar to first-year measures for swale plots (0.24–0.61 t ha⁻¹) and at the watershed scale (3.05 t ha⁻¹). Nearly all the sediment collected across all spatial scales was delivered during runoff events associated with cold-season hydrologic processes, including rain-on-snow, rain-on-frozen soils, and snowmelt runoff. Approximately 85–99% of annual sediment collected across all silt fence plots each year was from swales. The high levels of sediment delivered across hillslope to watershed scales in this study are attributed to observed preferential loading of fine sediments into swale channels by aeolian processes in the immediate post-fire period and subsequent flushing of these sediments by runoff from cold-season hydrologic processes. Our results suggest that the interaction of aeolian and cold-season hydrologic-driven erosion processes is an important component for consideration in post-fire erosion assessment and prediction and can have profound implications for soil loss from these ecosystems. © 2019 John Wiley & Sons, Ltd.     

GeoNat v1.0: A dataset for natural feature mapping with artificial intelligence and supervised learning

Machine learning allows “the machine” to deduce the complex and sometimes unrecognized rules governing spatial systems, particularly topographic mapping, by exposing it to the end product. Often, the obstacle to this approach is the acquisition of many good and labeled training examples of the desired result. Such is the case with most types of natural features. To address such limitations, this research introduces GeoNat v1.0, a natural feature dataset, used to support artificial intelligence‐based mapping and automated detection of natural features under a supervised learning paradigm. The dataset was created by randomly selecting points from the U.S. Geological Survey’s Geographic Names Information System and includes approximately 200 examples each of 10 classes of natural features. Resulting data were tested in an object‐detection problem using a region‐based convolutional neural network. The object‐detection tests resulted in a 62% mean average precision as baseline results. Major challenges in developing training data in the geospatial domain, such as scale and geographical representativeness, are addressed in this article. We hope that the resulting dataset will be useful for a variety of applications and shed light on training data collection and labeling in the geospatial artificial intelligence domain.     

Formation of low‐magnesium calcite at cold seeps in an aragonite sea

This study investigates the conditions of occurrence and petrographic characteristics of low‐Mg calcite (LMC) from cold seeps of the Gulf of Mexico at a water depth of 2340 m. Such LMC mineral phases should precipitate in calcite seas rather than today's aragonite sea. The ¹³C‐depleted carbonates formed as a consequence of anaerobic oxidation of hydrocarbons in shallow subsurface cold seep environments. The occurrence of LMC may result from brine fluid flows. Brines are relatively Ca²⁺‐enriched and Mg²⁺‐depleted (Mg/Ca mole ratio <0.7) relative to seawater, where the Mg/Ca mole ratio is ~5, which drives high‐Mg calcite and aragonite precipitation. The dissolution of aragonitic mollusk shells, grains and cements was observed. Aerobic oxidation of hydrocarbons and H₂S is the most likely mechanism to explain carbonate dissolution. These findings have important implications for understanding the occurrence of LMC in deep water marine settings and consequently their counterparts in the geological record.