Geology ReportsSearch

SEARCH · Geology Reports

Results for “Journal of Coastal Research”

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,660 records · Page 11Linked to original sources

Improving crop-specific groundwater use estimation in the Mississippi Alluvial Plain: Implications for integrated remote sensing and machine learning approaches in data-scarce regions

Study region The Mississippi Alluvial Plain (MAP) in the United States (US). Study focus Understanding local-scale groundwater use, a critical component of the water budget, is necessary for implementing sustainable water management practices. The MAP is one of the most productive agricultural regions in the US and extracts more than 11 km 3 /year for irrigation activities. Consequently, groundwater-level declines in the MAP region pose a substantial challenge to water sustainability, and hence, we need reliable groundwater pumping monitoring solutions to manage this resource appropriately. New hydrological insights for the region We incorporate remote sensing datasets and machine learning to improve an existing lookup table-based model of groundwater use previously developed by the U.S. Geological Survey (USGS). Here, we employ Distributed Random Forests, an ensemble machine learning algorithm to predict annual and monthly groundwater use (2014–2020) throughout this region at 1-km resolution, using pumping data from existing flowmeters in the Mississippi Delta. Our model compares favorably with the existing USGS model, with higher R 2 (0.51 compared to 0.42 in the previous model), and lower root mean square error (RMSE) and mean absolute error (MAE)— 0.14 m and 0.09 m, respectively in our model, compared to 0.15 m and 0.1 m in the previous model. Therefore, this work advances our ability to predict groundwater use in regions with scarce or limited in-situ groundwater withdrawal data availability.

Journal of Hydrology Regional Studies

Modeling seawater intrusion along the Alabama coastline using physical and machine learning models to evaluate the effects of multiscale natural and anthropogenic stresses

Seawater intrusion threatens groundwater resources in coastal regions, including southern Baldwin County, Alabama, where the freshwater-saltwater interface dynamics remain poorly understood. To address this gap, this study uses combined physics-based and machine-learning models to quantify seawater intrusion caused by natural (storm surges) and anthropogenic (human activities) perturbations. The long short-term memory network and wavelet analysis were used to assess vertical aquifer vulnerabilities, revealing that the shallow part of the Coastal lowlands aquifer system (CL1) in the southern Baldwin County region is more susceptible to sea level rise and groundwater extraction than deeper aquifers. Based on these findings, a cross-sectional numerical model (physics approach) for the CL1 aquifer was developed to evaluate tidal and storm surge effects, using Tropical Storm Claudette (June 2021) as a case study. Results showed that tidal fluctuations had a minimal impact on the saltwater-freshwater interface location, whereas storm surges caused substantial inland movement, with effects lasting for nine months. The steady-state version of the three-dimensional (3D) physical model predicted seawater intrusion across the entire area, and convolutional neural network-based modeling further validated the model results. The 3D physical model was also applied to a smaller area to assess human impact on the saltwater interface due to two groundwater pumping scenarios (± 50% of the baseline pumping rate). Results revealed that a 50% increase in groundwater withdrawals caused seawater to advance ~ 320 m inland, whereas a 50% reduction led to a ~ 270-meter retreat. This study highlights the vulnerability of Alabama’s shallow coastal aquifers to seawater intrusion due to storm surges and human activities, and demonstrates that combining physics-based models with machine learning approaches can improve groundwater predictions, though its accuracy depends on the availability of site-specific data.

Alabama

Earthquake-hazard exposure of residents with potential access and functional needs in the United States

Earthquake response plans and earthquake early warning (EEW) systems designed for general populations may not consider potential access and functional needs (AFN) of individuals with physical, sensory, cognitive, or social limitations. Previous efforts to map the distribution of these populations have focused on social-vulnerability indices that ignore or oversimply these limitations. The descriptive and exploratory analysis summarized in this United States (U.S.) case study addresses this gap by identifying and integrating spatially explicit data for AFN-related residential populations, earthquake hazards, and county and county equivalents for the conterminous U.S., Alaska, Puerto Rico, and Hawaii. We focus on 13 AFN-related attributes that relate to an individual's ability to access information contained in an EEW alert, to understand and process earthquake information or observed ground shaking, and to take self-protective actions based on this information and physical cues of an earthquake. Depending on the demographic attribute, there are millions to tens of millions of U.S. residents with AFN-related attributes in areas considered to have varying likelihoods (2%, 10%, and 50%) of exceedance of a damaging earthquake in the next 50 years. Although these amounts represent low percentages at the national level, the percentage of individuals with AFN-related attributes in many counties and county equivalents substantially exceeds national percentages. No one county, county equivalent, U.S. state, or U.S. territory has the highest percentage of individuals in all AFN-related attributes; therefore, future efforts to increase individual resilience to earthquakes may benefit from understanding the local context of individuals with potential access and functional needs.

International Journal of Disaster Risk Reduction

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

Top-down targeted network analysis of critical mineral commodities applied to international geochemistry database

The global demand for critical mineral commodities is rapidly increasing, making domestic production an important factor in supporting the economy and national security. Large scale, publicly available geochemical databases allow for the application of data informatics methods to interrogate critical mineral commodities data for correlations in deposit formation and distribution, particularly for identifying enrichment of multiple critical mineral commodities at the same deposit. In this study, we applied network analysis to the Critical Minerals Mapping Initiative (CMMI) ore geochemistry (Critical Minerals in Ores, CMiO) database to identify the high concentration (defined as 10× bulk crustal abundance) co-occurrence of different critical mineral commodities across a mineral system hierarchy from deposit environments to individual deposits. Identifying patterns or unique outliers in enrichment in network communities will allow for the location of secondary critical mineral commodity resources from under-utilized deposits. We find trends in the enrichment of critical mineral commodities in network-communities between the elements praseodymium (Pr), neodymium (Nd), terbium (Tb), and dysprosium (Dy) across multiple CMiO database deposit environments and groups down to specific deposit types and sites. A separate trend in network community deposition is observed as well between iridium (Ir) and platinum (Pt) in deposit environments, groups, types, and sites. Network analysis focused on critical minerals in magmatic-hydrothermal deposits identified multiple deposit sites from different deposit types within the CMiO database with concentrations of Dy, Nd, Tb, Pr, Ir, and Pt that are at least ten times greater than the crustal average. This approach can be applied to any target element(s) or deposit(s) of interest, allowing broad investigation of co-enriched critical mineral commodities.

Journal of Geochemical Exploration

Complex hydrology and variability of nitrogen sources in a karst watershed

Streams draining karst areas with rapid groundwater transit times may respond relatively quickly to nitrogen reduction strategies, but the complex hydrologic network of interconnected sinkholes and springs is challenging for determining the placement and effectiveness of management practices. This study aims to inform nitrogen reduction strategies in a representative agricultural karst setting of the Chesapeake Bay watershed (Fishing Creek watershed, Pennsylvania) with known elevated nitrate contamination and a previous documented groundwater residence time of less than a decade. During baseflow conditions, streamflow did not increase with drainage area. Headwaters and the main stem lost substantial flow to sinkholes until eventually discharging along large springs downstream. Seasonal hydrologic conditions shift the flow and nitrogen load spatially among losing and gaining stream sections. A compilation of nitrogen source inputs with the geochemistry and the pattern of enrichment of δ 15 N and δ 18 O suggest that the nitrogen in streams and springs during baseflow represents a mixture of manure, fertilizer, and wastewater sources with low potential for denitrification. The pH and calcite saturation index increased along generalized flow paths from headwaters to springs and indicate shorter groundwater residence times in baseflow during the spring versus summer. Given the substantial investment in management practices, fixed monitoring sites could incorporate synoptic water sampling to properly monitor long-term progress and help inform management actions in karst watersheds. Although karst watersheds have the potential to respond to nitrogen reduction strategies due to shorter groundwater residence times, high nitrogen inputs, effectiveness of conservation practices, and release of legacy nutrients within the karst cavities could confound progress of water quality goals.

Pennsylvania

Impact of gas/liquid phase change of CO2 during injection for sequestration

CO 2 sequestration in deep saline formations is an effective and important process to control the rapid rise in CO 2 emissions. The process of injecting CO 2 requires reliable predictions of the stress in the formation and the fluid pressure distributions – particularly since monitoring of the CO 2 migration is difficult – to mitigate leakage, prevent induced seismicity, and analyze wellbore stability. A key aspect of CO 2 is the gas–liquid phase transition at the temperatures and pressures of relevance to leakage and sequestration, which has been recognized as being critical for accurate predictions but has been challenging to model without ad hoc empiricisms. This paper presents a robust multiphase thermodynamics-based poromechanics model to capture the complex phase transition behavior of CO 2 and predict the stress and pressure distribution under super- and sub- critical conditions during the injection process. A finite element implementation of the model is applied to analyze the behavior of a multiphase porous system with CO 2 as it displaces the fluid brine phase. We find that if CO 2 undergoes a phase transition in the geologic reservoir, the spatial variation of the density is significantly affected, and the migration mobility of CO 2 decreases in the reservoir. A key feature of our approach is that we do not a priori assume the location of the CO 2 gas/liquid interface – or even if it occurs at all – but rather, this is a prediction of the model, along with the spatial variation of the phase of CO 2 and the change of the saturation profile due to the phase change.

Journal of the Mechanics and Physics of Solids

Land application of drill waste: A scope analysis

Drilling fluid waste land application, a process where drilling wastes are spread and tilled into the land surface, has become common in some petroleum-producing states, however, the potential benefits and risks of this practice are not well studied. Drilling fluids can be water- or oil-based and can have high concentrations of total soluble salts and total petroleum hydrocarbons. Comprehensive chemical characterization of these fluids is not well documented in the literature, and the extent of land application is largely unknown. We hypothesized that the land application of drill waste would fluctuate over time due to economic factors. To begin to understand the extent of historical and potential future land application, we analyzed data from over 5,800 drilling fluid land application permits collected by the Oklahoma Corporation Commission for years 2000, 2005, 2010, and 2015–2020. During the years studied, drilling fluid wastes were applied to more than 250,000 acres in Oklahoma, with over 54,000 thousand barrels (Mbbl) of liquids and nearly 21,000 Mbbl of solids applied. Land application is widespread (occurring in 59/77 counties), however recent drilling activity, land availability, and the economics of transportation have created conditions favorable for land application specifically in the Anadarko Basin. Land application can co-occur with sensitive areas, such as important groundwater and surface-water drinking sources and agricultural fields used for subsistence or feed crop production. Our approach for quantifying the extent of land application, along with further chemical characterization studies, can aid operators and land managers who are considering this practice in assessing the associated benefits and risks.

Oklahoma

Concentration-discharge relations and transient metal loads reveal spatiotemporal variability in solute-generation mechanisms in a mine-affected watershed

Concentration-discharge (CQ) relations are commonly used to understand geochemical and hydrologic controls on the generation of solutes in watersheds. Despite the widespread application of CQ relations, this technique has been infrequently applied to acid mine drainage (AMD) sites, but the CQ framework may allow mechanistic understanding of remedial outcomes such as impoundment of water within underground mines. Results of CQ analyses and changes in metal loads in an AMD affected watershed in Colorado, USA indicate that dissolved loads increased at many individual locations following water impoundment within mine workings. Although increased loads were observed at most individual locations, these increases were offset by a large decrease in loading from the largest mine. A loading analysis that included data from an instream monitoring location showed a statistically significant decrease in Fe and Zn after bulkhead emplacement, indicating a net positive effect of bulkheads. Streams generally displayed dilution CQ patterns whereas mines and springs showed either flushing or chemostatic patterns prior to bulkheading, which transitioned to chemostatic patterns following bulkheading, indicating a transition from dynamic to equilibrium geochemical processes. Saturation indices for sulfide and secondary minerals indicated that mines and springs were near equilibrium for phases including schwertmannite, fluorite, and gypsum. Saturation indices vary through time for mines suggesting progressive leaching of sulfide minerals as the mass of available minerals in the mine workings decreases. Together, these diverse analyses provide an integrated understanding of the variability in solute generating processes in this watershed and may inform remediation plans for similarly affected sites by indicating the nature of mineralogic controls on water quality.

Colorado

Remote sensing evapotranspiration in ensemble-based framework to enhance cascade routing and re-infiltration concept in integrated hydrological model applied to support decision making

Integrated hydrological models (IHMs) help characterize the complexity of surface–groundwater interactions. The cascade routing and re-infiltration (CRR) concept, recently applied to a MODFLOW 6 IHM, improved conceptualization and simulation of overland flow processes. The CRR controls the transfer of rejected infiltration and groundwater exfiltration from upslope areas to adjacent downslope areas where that water can be evaporated, re-infiltrated back to subsurface, or discharged to streams as direct runoff. The partitioning between these three components is controlled by uncertain parameters that must be estimated. Thus, by quantifying and reducing those uncertainties, next to uncertainties of the other model parameters (e.g. hydraulic and storage parameters), the reliability of the CRR is improved and the IHM is better suited for decision support modelling, the two key objectives of this work. To this end, the remotely sensed MODIS-ET product was incorporated into the calibration process for complementing traditional hydraulic head and streamflow observations. A total of approximately 150,000 observations guided the calibration of a 13-year MODFLOW 6 IHM simulation of the Sardon catchment (Spain) with daily stress periods. The model input uncertainty was represented by grid-cell-scale parameterization, yielding approximately 500,000 unknown input parameters to be conditioned. The calibration was carried out through an iterative ensemble smoother. Incorporating the MODIS-ET data improved the CRR implementation, and reduced uncertainties associated with other model parameters. Additionally, it significantly reduced the uncertainty associated with net recharge, a critical flux for water management that cannot be directly measured and rather is commonly estimated by IHM simulations.

Sardon catchment

Meet the people where they are: Assessing user needs for aftershock forecast products in El Salvador, Mexico and the United States

Aftershock forecasts can help communities reduce their seismic risk by conveying how many aftershocks can be expected following a large earthquake, and how the expected number of aftershocks and their corresponding ground shaking evolves over time and space. Prior work finds that graphical forecast products may communicate such information better than only text or numbers. To identify which visual products can serve multiple user groups, we held workshops with members of several professions, including emergency managers, engineers, critical infrastructure operators, public health specialists, science communicators, and more. We conducted these workshops in El Salvador, Mexico and the United States to understand which forecast products may be effective across different countries. In these workshops, users performed small-group activities to elicit the types of aftershock forecast information that would support decisions in their respective roles and how this information would optimally be displayed. Maps of shaking hazards were frequently requested across all professions and countries, even for dissimilar forecast uses. The design of these maps, and other product needs, were differentiated by profession, country and other factors. Other forecast products, including those showing forecasts about the magnitudes and time periods of future aftershocks, also served a variety of users, but for different types of decisions. We found a greater variation in user needs by profession than by country, and that user needs also vary with time, communication channel and other contextual factors. We discuss practical implications for user-centered visual communication of operational aftershock forecasts.

International Journal for Disaster Risk Reduction

Distinguishing natural from mining-related metal sources by including streambank groundwater data in a stream mass loading study

Distinguishing stream metal loading caused by mine features from that caused by natural background sources remains challenging, yet this distinction is essential for making effective remedial decisions at many legacy mine sites. We combine a stream tracer injection and synoptic sampling study with data from shallow near-stream groundwater wells to estimate left-bank versus right-bank metal loading contributions at the 100-m spatial scale. The study was performed in the East Mancos River, a mountain headwater stream in Colorado, USA. The dominant source of elevated stream metal concentrations could be either groundwater infiltration through right-bank Doyle Mine waste piles or natural acid rock drainage from hydrothermally altered bedrock located mainly on the left bank. For the five metals of concern (Cu, Al, Zn, Cd, and Mn), we find that 15 % of the load contributed by diffuse groundwater inputs in the section potentially influenced by Doyle mine waste originates from the right bank. This right-bank potential mine contribution equates to only 3 % of the total watershed instream load for these metals. Furthermore, apparent 3 H/ 3 He groundwater ages in segments contributing most of the right-bank metal loading are sufficiently old (9–12 yr) to suggest that infiltration through the waste piles, located only 140–180 m from the stream, is unlikely. Estimated potential Doyle mine loading contributions can therefore be considered maximum values. Study results thus indicate that Doyle mine waste piles are a minor source of metal loading under low-flow conditions, and streambank groundwater data can provide valuable additional information in stream mass loading studies.

Colorado

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

A spatiotemporal deep learning approach for predicting daily air-water temperature signal coupling and identification of key watershed physical parameters in a montane watershed

Seasonal shifts from runoff to groundwater dominance influence daily headwater stream temperatures, especially where local groundwater input is strong. This input buffers temperature during hot periods, supporting cold-water habitats. Recent studies use air–water temperature signal metrics to identify zones of strong stream–groundwater connectivity. While Previous studies used air–water signal ratios as proxies for groundwater influence but were limited to specific sites and periods, without dynamic forecasting. This study is the first to forecast daily A r as a spatiotemporal signal using a Graph Convolutional Network–Long Short-Term Memory (GCN-LSTM) model. The model was trained using hydroclimate data (air temperature, precipitation, shortwave radiation, streamflow) and watershed physical features (e.g., sand content, slope). Results showed high predictive skill, achieving R 2 (NSE, RMSE) of 0.86 (0.73, 0.0004) for one-day-ahead to 0.52 (0.50, 0.0009) for seven-days ahead forecasts. Prior studies often have not explicitly incorporated spatial hydrogeologic drivers, but this model explicitly incorporates them to assess their impact on A r forecasting and stream-groundwater connectivity. Feature analysis identified mean sand, elevation, slope, clay, and TWI as key predictors of A r . Stronger groundwater signals appeared in hillslopes, elevations, and tributaries, highlighting watershed influence on streamflow. However, limitations include reliance on historical air–water temperature patterns for training and limited representation of extreme climate conditions. Despite these limitations, unlike previous studies relying on measured in-situ stream and air temperature, this study forecasts A r directly from climate and physiographic features after training, avoiding in-situ data requirements. Findings aiding predictions of stream ecosystem resilience.

New York

Detrital zircons and the magmatic history of Viti Levu, Fiji

We integrate the existing detrital zircon data from multiple modern river sediment samples on Viti Levu, Fiji, with the most current available geological and topographic mapping of the respective river drainage basins to compare detrital populations with potential bedrock sources. The temporal and spatial variations in zircon geochemistry supplement what is known from igneous rocks and confirm the petrological differences between plutonic and volcanic rocks from the Eocene to early Oligocene (Yavuna age, >30 Ma), middle Oligocene to middle Miocene (Wainimala age, 30–12.5 Ma), late Miocene (Colo age, 12.5–6.5 Ma) and latest Miocene (Namosi age, 6.5–5 Ma). The >30 Ma Yavuna-age zircons are restricted to areas that drain the previously mapped Yavuna Group. The 30–12.5 Ma zircons are found across central Viti Levu from west to east, and the 30–15 Ma zircons have distinctively low U/Yb and high Dy/Yb ratios. They are the best radiometric evidence of widespread early to middle Miocene arc magmatism in Fiji that was relatively U-poor. Peak deconvolution of the Colo age zircons from individual basins suggests the following ages for undated or poorly dated plutons from central Viti Levu. The large Mavuvu pluton is probably composed of multiple intrusions in the 12–10 Ma range, the Waiqa pluton is probably ca 10 Ma, and the Noikoro pluton is probably ca 9 Ma. There are zircons from unknown plutonic or volcanic sources between 8 and 7 Ma in western Viti Levu that have distinct Eu/Eu* ratios. We attribute the highest U/Yb ratios in some Colo age zircons to crustal anatexis. Namosi-age zircons are abundant in the Medrausucu Group and can be found scattered across Viti Levu.

Viti Levu

Identification of representative earthquakes for probabilistic tsunami hazard analysis (PTHA) using earthquake rupture forecasts and machine learning

As probabilistic tsunami hazard analysis (PTHA) focuses more on assessments for localized, populous regions, techniques are needed to identify a subsample of representative earthquake ruptures to make the computational requirements for producing high-resolution hazard maps tractable. Moreover, the greatest epistemic uncertainty in seismic PTHA is related to source characterization, which is often poorly defined and subjective. We address these two salient issues by applying streamlined earthquake rupture forecasts (ERFs), based on combinatorial optimization methods, to an unsupervised machine learning workflow for identifying representative ruptures. ERFs determine the optimal distribution of a millennia-scale sample of earthquakes by inverting the observed slip rate on major faults. We use two previously developed combinatorial optimization ERFs, integer programming and greedy sequential, to produce the optimal location of ruptures with seismic moments sampled from a regional Gutenberg–Richter magnitude–frequency distribution. These ruptures in turn are used to calculate peak nearshore tsunami amplitude, using computationally efficient tsunami Green's functions. An unsupervised machine learning workflow is then used to identify a small subsample of the earthquakes input to ERFs for onshore PTHA analysis. We eliminate epistemic uncertainty related to source distribution under traditional PTHA analysis; in its place, a quantifiable, less subjective and generally smaller uncertainty related to the input to ERFs is included. The Nankai subduction zone is used as a test case, where previous ERFs have been conducted. Results indicate that the locations of representative earthquakes are sensitive to choice of magnitude–area relation and to whether a minimum cumulative stress objective is imposed on the fault. In general, incorporating ERFs into PTHA provide a physically self-consistent method to incorporate fault slip information in determining representative earthquakes for onshore PTHA, eliminating a major source of epistemic uncertainty.

Nankai subduction zone

Upper Triassic igneous rocks of the southern Kenai Peninsula, Alaska – Prelude to Early Jurassic subduction along the western Wrangellia composite terrane margin

New U–Pb zircon geochronology identifies a latest Triassic (ca 214–201 Ma) igneous suite of tuff, hypabyssal dikes, and a pluton on the southern Kenai Peninsula, Alaska. The igneous suite was emplaced within Upper Triassic sedimentary rocks along the southern margin of Western Wrangellia, the western-most fragment of the Wrangellia composite terrane. The igneous rocks range from mafic (50.6% SiO 2 ) to felsic (78.3% SiO 2 ), characteristically have less than 1.55% K 2 O, and generally have low trace element abundances. The tonalitic and trondhjemitic magmas were largely sourced in mafic-rich lower crust and incompletely assimilated quartz and other mineral xenocrysts are common. Fractionation involving plagioclase and amphibole is indicated for some magmas and composite intrusions and igneous xenoliths indicate magma mixing was possible. Paleozoic and Precambrian inherited zircons and initial 87 Sr/ 86 Sr (0.704103–0.705609) and 143 Nd/ 144 Nd (0.512396–0.512777) ratios indicate that the Western Wrangellia crustal sources are heterogeneous and contain sialic components. The latest Triassic magmatism reflects processes that preceded Early Jurassic subduction along the Wrangellia composite terrane and Pacific Ocean plate boundary. These processes involved heating and melting of mantle lithosphere and lower crust as mantle instabilities accompanied the breaking of the plate boundary linkages. The Late Triassic transition to subduction along the Wrangellia composite terrane margin coincided with the transition to subduction cessation in the Late Triassic arcs of the western Intermontane terranes of Canada. The shift to subduction along the outboard Wrangellia composite terrane margin marks the beginning of the Pacific Ocean–Cordillera plate interactions that came to dominate the tectonic evolution of the northern Cordillera from the Early Jurassic to today.

Alaska