SEARCH · USGS Science

Results for “Journal of Maps”

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

SourceResearch connectionFrom yearThrough year

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

At least 307 records · Page 17Linked to original sources

USGS2009

Mapping Curie temperature depth in the western United States with a fractal model for crustal magnetization

We have revisited the problem of mapping depth to the Curie temperature isotherm from magnetic anomalies in an attempt to provide a measure of crustal temperatures in the western United States. Such methods are based on the estimation of the depth to the bottom of magnetic sources, which is assumed to correspond to the temperature at which rocks lose their spontaneous magnetization. In this study, we test and apply a method based on the spectral analysis of magnetic anomalies. Early spectral analysis methods assumed that crustal magnetization is a completely uncorrelated function of position. Our method incorporates a more realistic representation where magnetization has a fractal distribution defined by three independent parameters: the depths to the top and bottom of magnetic sources and a fractal parameter related to the geology. The predictions of this model are compatible with radial power spectra obtained from aeromagnetic data in the western United States. Model parameters are mapped by estimating their value within a sliding window swept over the study area. The method works well on synthetic data sets when one of the three parameters is specified in advance. The application of this method to western United States magnetic compilations, assuming a constant fractal parameter, allowed us to detect robust long-wavelength variations in the depth to the bottom of magnetic sources. Depending on the geologic and geophysical context, these features may result from variations in depth to the Curie temperature isotherm, depth to the mantle, depth to the base of volcanic rocks, or geologic settings that affect the value of the fractal parameter. Depth to the bottom of magnetic sources shows several features correlated with prominent heat flow anomalies. It also shows some features absent in the map of heat flow. Independent geophysical and geologic data sets are examined to determine their origin, thereby providing new insights on the thermal and geologic crustal structure of the western United States.

Journal of Geophysical Research B: Solid Earth↗

USGS2011

The magnitude distribution of earthquakes near Southern California faults

We investigate seismicity near faults in the Southern California Earthquake Center Community Fault Model. We search for anomalously large events that might be signs of a characteristic earthquake distribution. We find that seismicity near major fault zones in Southern California is well modeled by a Gutenberg-Richter distribution, with no evidence of characteristic earthquakes within the resolution limits of the modern instrumental catalog. However, the b value of the locally observed magnitude distribution is found to depend on distance to the nearest mapped fault segment, which suggests that earthquakes nucleating near major faults are likely to have larger magnitudes relative to earthquakes nucleating far from major faults. Copyright 2011 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2019

Three-dimensional basin and fault structure from a detailed seismic velocity model of Coachella Valley, Southern California

The Coachella Valley in the northern Salton Trough is known to produce destructive earthquakes, making it a high seismic hazard area. Knowledge of the seismic velocity structure and geometry of the sedimentary basins and fault zones is required to improve earthquake hazard estimates in this region. We simultaneously inverted first P wave travel times from the Southern California Seismic Network (39,998 local earthquakes) and explosions (251 land/sea shots) from the 2011 Salton Seismic Imaging Project to obtain a 3-D seismic velocity model. Earthquakes with focal depths ≤10 km were selected to focus on the upper crustal structure. Strong lateral velocity contrasts in the top ~3 km correlate well with the surface geology, including the low-velocity (<5 km/s) sedimentary basin and the high-velocity crystalline basement rocks outside the valley. Sediment thickness is ~4 km in the southeastern valley near the Salton Sea and decreases to <2 km at the northwestern end of the valley. Eastward thickening of sediments toward the San Andreas fault within the valley defines Coachella Valley basin asymmetry. In the Peninsular Ranges, zones of relatively high seismic velocities (~6.4 km/s) between 2 to 4 km depth may be related to Late Cretaceous mylonite rocks or older inherited basement structures. Other high-velocity domains exist in the model down to 9 km depth and help define crustal heterogeneity. We identify a potential fault zone in Lost Horse Valley unassociated with mapped faults in Southern California from the combined interpretation of surface geology, seismicity, and lateral velocity changes in the model.

California↗

USGS2024

Framework for mapping liquefaction hazard–Targeted design ground motions

Liquefaction-induced ground failure poses substantial challenges to geotechnical earthquake engineering design. Current approaches for designing against liquefaction hazards, as specified in most seismic provisions, focus on estimating a liquefaction factor of safety ( 𝐹⁢𝑆𝐿 ) and typically characterize earthquake loading using design parameters based on probabilistic or deterministic ground motion levels. Because 𝐹⁢𝑆𝐿 is estimated deterministically, this basis of design neglects considerable uncertainties for estimating liquefaction triggering and its consequences and results in a lack of liquefaction-specific design criteria, particularly as structural design has advanced toward risk-targeted performance objectives. This study presents a framework for developing liquefaction-targeted design criteria based on a minimum acceptable return period of liquefaction, informed by probabilistic liquefaction hazard analysis (PLHA). PLHA quantifies annualized rates of liquefaction by considering contributions from (1) the full ground-motion probability space, and (2) uncertainties in liquefaction triggering using probabilistic models. PLHA is used in this study to characterize the current, effective return periods of 𝐹⁢𝑆𝐿 ( 𝑇𝑅,𝐹⁢𝑆 ) obtained from conventional liquefaction hazard analysis (CLHA) using uniform-hazard ground motions. 𝑇𝑅,𝐹⁢𝑆 is evaluated in a parametric study of nearly 100 sites throughout the conterminous United States. The results indicate large geographic variations in acceptable liquefaction hazard levels, with implied 𝑇𝑅,𝐹⁢𝑆 ranging between approximately 1,000 to 3,000 years. To address these inconsistencies without the computational demands of full PLHA, a framework is proposed for developing a liquefaction-targeted design peak ground acceleration, 𝑃⁢𝐺⁢𝐴𝐿 , for use in liquefaction models that result in consistent liquefaction design levels across all geographic locations. The mapped 𝑃⁢𝐺⁢𝐴𝐿 is shown to be somewhat sensitive to site-specific properties, and adjustment factors are developed and presented. The proposed 𝑃⁢𝐺⁢𝐴𝐿 mapping procedure produces 𝐹⁢𝑆𝐿 estimates that are consistent with those obtained from full PLHA at a target 𝑇𝑅,𝐹⁢𝑆 , providing a promising roadmap to incorporating PLHA concepts into current liquefaction design methods.

Journal of Geotechnical and Geoenvironmental Engin↗

USGS2003

Mapping the seabed and habitats in National Marine Sanctuaries - Examples from the East, Gulf and West Coasts

The National Marine Sanctuary System requires seabed and habitat maps to serve as a basis for managing sanctuary resources and for conducting research. NOAA, the agency that manages the sanctuaries, and the USGS have conducted mapping projects in three sanctuaries (Stellwagen Bank NMS, Flower Garden Banks NMS, and Channel Islands NMS) with an emphasis on collaboration of geologists and biologists from the two agencies and from academic institutions. Mapping of seabed habitats is a developing field that requires the integration of geologic and biologic studies and the use of swath imaging techniques such as multibeam and sidescan sonar. Major products of swath mapping are shaded-relief topographic imagery which shows seabed features in great detail, and backscatter imagery which provides an indication of the types of materials that constitute the seabed. Sea floor images provide an excellent basis for conducting the groundtruthing studies (using video, photo, and sampling techniques) that are required to collect the data necessary for making meaningful interpretative maps of the seabed. The compilation of interpretive maps showing seabed environments and habitats also requires the development of a sea floor classification system that will be a basis for comparing, managing, and researching characteristic areas of the seabed. Seabed maps of the sanctuaries are proving useful for management and research decisions that address commercial and recreational fishing, habitat disturbance, engineering projects, tourism, and cultural resources.

Marine Technology Society Journal↗

USGS2015

Mapping of coal quality using stochastic simulation and isometric logratio transformation with an application to a Texas lignite

Coal is a chemically complex commodity that often contains most of the natural elements in the periodic table. Coal constituents are conventionally grouped into four components (proximate analysis): fixed carbon, ash, inherent moisture, and volatile matter. These four parts, customarily measured as weight losses and expressed as percentages, share all properties and statistical challenges of compositional data. Consequently, adequate modeling should be done in terms of a logratio transformation, a requirement that is commonly overlooked by modelers. The transformation of choice is the isometric logratio transformation because of its geometrical and statistical advantages. The modeling is done through a series of realizations prepared by applying sequential simulation for the purpose of displaying the parts in maps incorporating uncertainty. The approach makes realistic assumptions and the results honor the data and basic considerations, such as percentages between 0 and 100, all four parts adding to 100% at any location in the study area, and a style of spatial fluctuation in the realizations equal to that of the data. The realizations are used to prepare different results, including probability distributions across a deposit, E-type maps displaying average properties, and probability maps summarizing joint fluctuations of several parts. Application of these maps to a lignite bed clearly delineates the deposit boundary, reveals a channel cutting across, and shows that the most favorable coal quality is to the north and deteriorates toward the southeast.

International Journal of Coal Geology↗

USGS2016

Calorific value and compositional ultimate analysis with a case study of a Texas lignite

Measurements to determine coal quality as fuel include proximate analysis, ultimate analysis and calorific value. The latter is an attribute taking non-negative real values, so a simple transformation is sufficient for its spatial modeling applying geostatistics. The analyses, however, involve proportions that follow the properties of compositional data, thus requiring special preprocessing for an adequate modeling already described in a previous publication for the case of proximate analysis data. 1 Here we model the results of calorific value and ultimate analysis. We propose to use two different binary partitions, one per analysis, map the corresponding isometric logratio transformations, and backtransform the results. The methodology is illustrated using the same coal bed in the previous paper modeling proximate analysis data. Results are summarized using probability maps that, in the case of this deposit, show a prominent channel crossing the deposit and separating the best quality coal from that of lower quality.

International Journal of Coal Geology↗

USGS2024

Hyperspectral (VNIR-SWIR) analysis of roll front uranium host rocks and industrial minerals from Karnes and Live Oak Counties, Texas Coastal Plain

VNIR-SWIR (400–2500 nm) reflectance measurements were made on the surfaces of various cores, cuttings and sample splits of sedimentary rocks from the Tertiary Jackson Group, and Catahoula, Oakville and Goliad Formations. These rocks vary in composition and texture from mudstone and claystone to sandstone and are known host rocks for roll front uranium occurrences in Karnes and Live Oak Counties, Texas. Spectral reflectance profiles, 569 in total, were reduced to 125 representative spectral signatures, which were analyzed using the U.S. Geological Survey's (USGS) Material Identification and Characterization Algorithm (MICA). MICA uses an automated continuum-removal procedure together with a least-squares linear regression to determine the fit of observed sample spectral absorption features to those of reference mineral standards in a spectral library. The reference minerals include various clay, mica, carbonate, ferric and ferrous iron minerals and their mixtures. In addition, absorption feature band-depth analysis was done to identify rock surfaces exhibiting absorption features related to uranium and zeolite minerals, which were not included in the command files used to execute MICA. Rocks from each of the four geologic units produced broadly similar spectral signatures as a result of comparable mineral compositions, but there were some notable differences. For example, Ca- and Na-montmorillonite was matched most frequently to the spectral absorption features in 2-μm (∼2000–2500 nm) wavelengths, while goethite occurred often at 1-μm (∼400–1000 nm) wavelengths. The latter is related to limonitic iron-staining in and around oxidized zones of the uranium roll front as described in previous papers. Rocks of the Jackson Group differed from those of the Catahoula, Oakville and Goliad units in that the former exhibited spectral features we interpret as being due to the presence of lignite-bearing mudstone layers. Goliad rocks exhibit spectral features related to dolomite, gypsum, anhydrite, and an unidentified green clay mineral that is possibly glauconite . Jackson Group rocks also exhibit weak but well-resolved absorption features at 964 and 1157 nm related to either or both zeolite minerals clinoptilolite and heulandite. These zeolite minerals and a few spectra exhibiting hydrous silica absorption features are indicative of alteration of volcanic glass in tuffaceous mudstone and claystone layers. A few sample spectra exhibited strong absorption features at around 1135 nm related to the uranium mineral coffinite. Both the 1135 nm coffinite and 1157 nm zeolite absorption features overlap somewhat, potentially making them difficult to distinguish without additional hyperspectral field, laboratory or remote sensing data. The results of this study were compared to mixtures of minerals described for ore, gangue and alteration minerals in deposit models for sandstone-hosted uranium, sedimentary bentonite and sedimentary zeolite. Use of these spectra can help facilitate mapping of both waste materials from the legacy mining of the above commodities, as well as future exploration and resource assessment activities.

Texas↗

USGS2008

Surface drift prediction in the Adriatic Sea using hyper-ensemble statistics on atmospheric, ocean and wave models: Uncertainties and probability distribution areas

Despite numerous and regular improvements in underlying models, surface drift prediction in the ocean remains a challenging task because of our yet limited understanding of all processes involved. Hence, deterministic approaches to the problem are often limited by empirical assumptions on underlying physics. Multi-model hyper-ensemble forecasts, which exploit the power of an optimal local combination of available information including ocean, atmospheric and wave models, may show superior forecasting skills when compared to individual models because they allow for local correction and/or bias removal. In this work, we explore in greater detail the potential and limitations of the hyper-ensemble method in the Adriatic Sea, using a comprehensive surface drifter database. The performance of the hyper-ensembles and the individual models are discussed by analyzing associated uncertainties and probability distribution maps. Results suggest that the stochastic method may reduce position errors significantly for 12 to 72??h forecasts and hence compete with pure deterministic approaches. ?? 2007 NATO Undersea Research Centre (NURC).

Journal of Marine Systems↗

USGS2010

What you should know about land-cover data

Wildlife biologists are using land-characteristics data sets for a variety of applications. Many kinds of landscape variables have been characterized and the resultant data sets or maps are readily accessible. Often, too little consideration is given to the accuracy or traits of these data sets, most likely because biologists do not know how such data are compiled and rendered, or the potential pitfalls that can be encountered when applying these data. To increase understanding of the nature of land-characteristics data sets, I introduce aspects of source information and data-handling methodology that include the following: ambiguity of land characteristics; temporal considerations and the dynamic nature of the landscape; type of source data versus landscape features of interest; data resolution, scale, and geographic extent; data entry and positional problems; rare landscape features; and interpreter variation. I also include guidance for determining the quality of land-characteristics data sets through metadata or published documentation, visual clues, and independent information. The quality or suitability of the data sets for wildlife applications may be improved with thematic or spatial generalization, avoidance of transitional areas on maps, and merging of multiple data sources. Knowledge of the underlying challenges in compiling such data sets will help wildlife biologists to better assess the strengths and limitations and determine how best to use these data.

Journal of Wildlife Management↗

USGS2016

Integrating remote sensing with species distribution models; Mapping tamarisk invasions using the Software for Assisted Habitat Modeling (SAHM)

Early detection of invasive plant species is vital for the management of natural resources and protection of ecosystem processes. The use of satellite remote sensing for mapping the distribution of invasive plants is becoming more common, however conventional imaging software and classification methods have been shown to be unreliable. In this study, we test and evaluate the use of five species distribution model techniques fit with satellite remote sensing data to map invasive tamarisk ( Tamarix spp.) along the Arkansas River in Southeastern Colorado. The models tested included boosted regression trees (BRT), Random Forest (RF), multivariate adaptive regression splines (MARS), generalized linear model (GLM), and Maxent. These analyses were conducted using a newly developed software package called the Software for Assisted Habitat Modeling (SAHM). All models were trained with 499 presence points, 10,000 pseudo-absence points, and predictor variables acquired from the Landsat 5 Thematic Mapper (TM) sensor over an eight-month period to distinguish tamarisk from native riparian vegetation using detection of phenological differences. From the Landsat scenes, we used individual bands and calculated Normalized Difference Vegetation Index (NDVI), Soil-Adjusted Vegetation Index (SAVI), and tasseled capped transformations. All five models identified current tamarisk distribution on the landscape successfully based on threshold independent and threshold dependent evaluation metrics with independent location data. To account for model specific differences, we produced an ensemble of all five models with map output highlighting areas of agreement and areas of uncertainty. Our results demonstrate the usefulness of species distribution models in analyzing remotely sensed data and the utility of ensemble mapping, and showcase the capability of SAHM in pre-processing and executing multiple complex models.

Journal of Visualized Experiments↗

USGS2024

Estimates of soil taxonomic change due to near-surface permafrost loss in Alaska

Gelisols (permafrost-affected soils in U.S. Soil Taxonomy) are extensive in Alaska, currently occurring on approximately 45% of the land area of the state. Gelisol taxonomic criteria rely on the presence of near-surface (less than 2m deep) permafrost, but ongoing climatic and environmental change has the potential to affect the presence of near-surface permafrost across much of Alaska throughout the 21st century. In this study, we utilized scenarios of near-surface permafrost loss and active layer deepening through the 21st century under low (SRES B1, RCP 4.5), mid- (SRES A1B) and high (SRES A2, RCP8.5) emissions scenarios, in conjunction with the statewide STATSGO soil map, to generate spatially-explicit predictions of the susceptibility of Gelisols and Gelisol suborders to taxonomic change in Alaska. We find that 15-53% of Alaskan Gelisols are susceptible to taxonomic change by mid-century and 41-69% of Alaskan Gelisols are susceptible to taxonomic change by the end of the century. The extent of potential change varies between suborders and geographic regions, with Gelisols in Northern Alaska being the most resilient to taxonomic change and Western and Interior Alaskan Gelisols most susceptible to taxonomic change. The Orthel suborder is likely to be highly restricted by the late 21st century, while Histels and Tubels are more likely to be of greater extent. These results should be taken into consideration when designing initial survey and re-mapping efforts in Alaska, and suggest that Alaskan Gelisol taxa should be considered threatened soil taxa due to the proportional extent of likely loss.

Alaska↗

USGS2023

Providing a framework for seagrass mapping in United States coastal ecosystems using high spatial resolution satellite imagery

Seagrasses have been widely recognized for their ecosystem services, but traditional seagrass monitoring approaches emphasizing ground and aerial observations are costly, time-consuming, and lack standardization across datasets. This study leveraged satellite imagery from Maxar's WorldView-2 and WorldView-3 high spatial resolution, commercial satellite platforms to provide a consistent classification approach for monitoring seagrass at eleven study areas across the continental United States, representing geographically, ecologically, and climatically diverse regions. A single satellite image was selected at each of the eleven study areas to correspond temporally to reference data representing seagrass coverage and was classified into four general classes: land, seagrass, no seagrass, and no data. Satellite-derived seagrass coverage was then compared to reference data using either balanced agreement, the Mann-Whitney U test, or the Kruskal-Wallis test, depending on the format of the reference data used for comparison. Balanced agreement ranged from 58% to 86%, with better agreement between reference- and satellite-indicated seagrass absence (specificity ranged from 88% to 100%) than between reference- and satellite-indicated seagrass presence (sensitivity ranged from 17% to 73%). Results of the Mann-Whitney U and Kruskal-Wallis tests demonstrated that satellite-indicated seagrass percentage cover had moderate to large correlations with reference-indicated seagrass percentage cover, indicative of moderate to strong agreement between datasets. Satellite classification performed best in areas of dense, continuous seagrass compared to areas of sparse, discontinuous seagrass and provided a suitable spatial representation of seagrass distribution within each study area. This study demonstrates that the same methods can be applied across scenes spanning varying seagrass bioregions, atmospheric conditions, and optical water types, which is a significant step toward developing a consistent, operational approach for mapping seagrass coverage at the national and global scales. Accompanying this manuscript are instructional videos describing the processing workflow, including data acquisition, data processing, and satellite image classification. These instructional videos may serve as a management tool to complement field- and aerial-based mapping efforts for monitoring seagrass ecosystems.

Alaska, California, Delaware, Florida, Maine, Mary↗

USGS2007

Northward extension of Carolina slate belt stratigraphy and structure, South-Central Virginia: Results from geologic mapping

Geologic mapping in south-central Virginia demonstrates that the stratigraphy and structure of the Carolina slate belt extend northward across a steep thermal gradient into upper amphibolite-facies correlative gneiss and schist. The Neoproterozoic greenschist-facies Hyco, Aaron, and Virgilina Formations were traced northward from their type localities near Virgilina, Virginia, along a simple, upright, northeast-trending isoclinal syncline. This syncline is called the Dryburg syncline and is a northern extension of the more complex Virgilina synclinorium. Progressively higher-grade equivalents of the Hyco and Aaron Formations were mapped northward along the axial trace of the refolded and westwardly-overturned Dryburg syncline through the Keysville and Green Bay 7.5-minute quadrangles, and across the northern end of the Carolina slate belt as interpreted on previous geologic maps. Hyco rocks, including felsic metatuff, metawacke, and amphibolite, become gneisses upgrade with areas of local anatexis and the segregation of granitic melt into leucosomes with biotite selvages. Phyllite of the Aaron Formation becomes garnet-bearing mica schist. Aaron Formation rocks disconformably overlie the primarily felsic volcanic and volcaniclastic rocks of the Hyco Formation as evidenced by repeated truncation of internal contacts within the Hyco on both limbs of the Dryburg syncline at the Aaron-Hyco contact. East-northeast-trending isograds, defined successively by the first appearance of garnet, then kyanite ?? staurolite in sufficiently aluminous rocks, are superposed on the stratigraphic units and synclinal structure at moderate to high angles to strike. The textural distinction between gneisses and identifiable sedimentary structures occurs near the kyanite ?? staurolite-in isograd. Development of the steep thermal gradient and regional penetrative fabric is interpreted to result from emplacement of the Goochland terrane adjacent to the northern end of the slate belt during Alleghanian orogenesis. This mapping study indicates that the Carolina slate belt does not terminate on the north against through-going faults or rest on higher-grade basement as previously suggested.

American Journal of Science↗

USGS2005

Geologic framework, evolution, and sediment resources for restoration of the Louisiana Coastal Zone

The Louisiana Coastal Zone along the north-central Gulf of Mexico represents one of America's most important coastal ecosystems in terms of natural resources, human infrastructure, and cultural heritage. This zone also has the highest rates of coastal erosion and wetland loss in the nation because of a complex combination of natural processes and anthropogenic activities during the past century. In response to the dramatic land loss, regional-scale restoration plans are being developed through a partnership of federal and state agencies. One objective is to maintain the barrier island and tidal inlet systems, thereby reducing the impact of storm surge and interior wetland loss. Proposed shore line restoration work relies primarily upon the use of large volumes of sand-rich sediment for shoreline stabilization and the implementation of the shoreline projects. Although sand-rich sediment is required for the Louisiana restoration projects, it is of limited availability within the generally clay to silt-rich, shallow strata of the Louisiana Coastal Zone. Locating volumetrically significant quantities of sand-rich sediment presents a challenge and requires detailed field investigations using direct sampling and geophysical sensing methods. Consequently, there is a fundamental need to thoroughly understand and map the distribution and textural character {e.g., sandiness) of sediment resources within the Coastal Zone for the most cost-effective design and completion of restoration projects.

Gulf of Mexico↗

USGS2016

Thermal mapping of a pahoehoe lava flow, Kilauea Volcano

Pāhoehoe lava flows are a major component of Hawaiian eruptive activity, and an important part of basaltic volcanism worldwide. In recent years, pāhoehoe lava has destroyed homes and threatened parts of Hawai‘i with inundation and disruption. In this study, we use oblique helicopter-borne thermal images to create high spatial resolution (~1 m) georeferenced thermal maps of the active pāhoehoe flow on Kīlauea Volcano’s East Rift Zone. Thermal maps were created on 27 days during 2014–2016 in the course of operational monitoring, encompassing a phase of activity that threatened the town of Pāhoa. Our results illustrate and reinforce how pāhoehoe flows are multicomponent systems consisting of the vent, master tube, distributary tubes and surface breakouts. The thermal maps accurately depict the distribution and character of pāhoehoe breakouts through time, and also delineate the subsurface lava tube. Surface breakouts were distributed widely across the pāhoehoe flow, with significant portions concurrently active well upslope of the flow front, often concentrated in clusters of activity that evolved through time. Gradual changes to surface breakout distribution and migration relate to intrinsic processes in the flow, including the slow evolution of the distributary tube system. Abrupt disruptions to this system, and the creation of new breakouts (and associated hazards), were triggered by extrinsic forcing—namely fluctuations in lava supply rate at the vent which disrupted the master lava tube. Although the total area of a pāhoehoe flow has been suggested to relate to effusion rate, our results show that changes in the proportion of expansion vs. overplating can complicate this relationship. By modifying existing techniques, we estimate time-averaged discharge rates for the flow during 2014–2016 generally in the range of 1–2 m3 s-1 (mean: 1.3±0.4 m3 s-1) – less than half of Kīlauea’s typical eruption rate on the East Rift Zone and suggestive of a weak eruptive regime during 2014–2016. We caution, however, that this discharge rate approach requires further independent corroboration. The thermal maps provide the first synoptic characterization of pāhoehoe flow activity at high spatial resolution, essential both for operational hazard assessment and fundamental understanding of pāhoehoe behavior.

Hawaii↗

USGS2022

Can lava flow like water? Assessing applications of critical flow theory to channelized basaltic lava flows

Flowing lava and water have dramatically different physical properties but can form similar hydraulic structures, including undular hydraulic jumps, or standing wave trains. In water flows, undular hydraulic jumps are evidence of critical flow (Froude number ∼1) and open-channel hydraulic theory provides a powerful tool for estimating flow depth and velocity. Monitoring these parameters in an active lava channel is inherently challenging, but essential for calculating lava discharge (effusion rate), a primary control on the rate of flow front advance and ultimate flow runout distance. We analyze undular hydraulic jumps in both water and lava flows to assess the conditions under which they form and, by extension, the potential use of critical flow theory to estimate, in real time, lava flow velocity, depth, and discharge. Experimental data for water flows show that these structures mark the transition from supercritical to subcritical flow. Undular hydraulic jumps in the near-vent lava channel of the 2018 lower East Rift Zone eruption of Kīlauea, Hawaiʻi also reflect critical flow conditions; their wavelengths scale with flow depth and velocity, consistent with hydraulic theory. Calculated lava effusion rates are similar to estimates made using more traditional approaches (Jeffreys', 1925, https://doi.org/10.1080/14786442508634662 , equation based on lava viscosity, density, and channel slope) and with lava volumes derived from topographic-change mapping. From this we conclude that critical flow phenomena show great potential to track flow dynamics and inform hazard assessment for a wide range of geophysical fluids.

Journal of Geophysical Research - Earth Surface↗

USGS2012

Global petrologic variations on the moon: A ternary-diagram approach.

A ternary-diagram approach for determination of global petrologic variations on the lunar surface is presented that incorporates valuable improvements in our previous method of using geochemical variation diagrams. This new approach uses a ternary diagram that is subdivided into equally spaced segments along each of its three sides and has a triangular area in its center. The segments are assigned distinct colors that cover the range of the visible spectrum. The apexes are assigned the three primary colors, and the center triangle, which represents equal proportions of the three apexes, is assigned the color gray. The apexes are assigned the average Fe and (Th/Ti) c (ratio normalized to chondrites) compositions of KREEP/Mg-suite rocks, mare basalts, and ferroan anorthosites. The composition of each picture element (pixel) within these Apollo orbital gamma ray databases is used to determine its position within the ternary diagram. The color corresponding to this ternary position is then placed on a classification map at that pixel's position within the orbital databases. Error databases produced for the Fe and (Th/Ti) c data were used to exclude pixels having high errors. The resultant classification map shows spatial transitions among petrologic units, allows direct determination of the relative proportions of each end-member composition in a pixel, and increases the geologic interpretability of these data over that of the individual elemental databases. The classification units correspond remarkably well to observed geologic units, when we consider the spatial resolution (100 km) of the gamma ray spectrometer. Our results are as follows: (1) The highlands contain large areas of relatively pure ferroan anorthosite. (2) The average composition of the upper lunar crust is represented by an “anorthositic gabbro” composition, which supports the “magma-ocean” hypothesis for lunar crustal genesis. (3) KREEP/Mg-suite rocks are a minor fraction of the upper lunar crust. (4) Within the farside highlands, areas of KREEP/Mg-suite rocks coincide mostly with areas of crustal thinning, which are probably areas of KREEP basalt extrusions or localized excavations of Mg-suite rocks or KREEP-rich rocks. (5) Portions of the east limb and farside highlands have considerable amounts of a mafic, chondritic Th/Ti component (like mare basalt) whose occurrences coincide with mapped concentrations of light plains that display dark-halo craters; the presence of this component supports the hypothesis that mare volcanism had occurred within this highland region before the end of final heavy bombardment.

Journal of Geophysical Research Solid Earth↗