USGS ScienceSearch

SEARCH · USGS Science

Results for “Journal of Volcanology and Geothermal 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,689 records · Page 35Linked to original sources

Hazard potential of compound flooding from rainfall, storm surge, and groundwater in coastal New York and Connecticut

Compound flood events, the co-occurrence of multiple flood drivers, can result in flood hazard potential exceeding that of any single driver alone. To evaluate compound flooding in a semi-urbanized coastal area, historical records dating back to 1970 are used to study the co-occurrences of high precipitation, storm surge, and shallow groundwater conditions along the coastlines of New York and Connecticut. Joint return periods for coincident precipitation-surge events were computed using statistical dependence models and compared to the assumption of independence as a ratio, referred to here as a return period adjustment. Results indicate distinct seasonality where compound events in the area disproportionately occur in the cold season between October and April. Return period adjustments range from a factor of 1 to almost 9, demonstrating the range in precipitation-storm surge dependence across the study area. Across all 24 station triad locations, groundwater levels were elevated during times of precipitation- surge co-occurrence, reflecting the tendency for coastal storms and shallow groundwater conditions to co-occur seasonally. The result is a pseudo-trivariate compound flood hazard score and corresponding hazard map that integrates dependence between daily precipitation-surge events and overall monthly groundwater levels (as a precondition) into a relative compound hazard score. The location with the highest compound flood hazard score is on the south shore of Long Island, as well as locations across coastal Connecticut where groundwater levels compound the co-occurrence of heavy precipitation and storm surge.

Connecticut, New York

Conceptual hydrogeologic framework and groundwater budget near the southeastern part of Puget Sound, Washington

More than 1 million people live within the active model area (AMA) in the southeastern part of the lowlands surrounding Puget Sound, or Puget Lowland, Washington, and groundwater is the source for approximately one-half of their public, domestic, and irrigation water demands. The 887-square-mile AMA, located in King and Pierce Counties, represents the area of analysis for the conceptual hydrogeologic framework and numerical groundwater-flow models within the study area and includes the Puyallup River and Chambers-Clover Creek watersheds. To assess the potential hydrologic and anthropogenic impacts to groundwater and the connected surface-water resources, conceptual and numerical groundwater-flow models of groundwater flow were developed by the U.S. Geological Survey Washington Water Science Center in close cooperation with 18 water-resource agencies and stakeholders. This report presents information used to characterize the groundwater-flow system and the development of a numerical model in the AMA. Included are descriptions of the geology and conceptual hydrogeologic framework, groundwater levels and flow directions, groundwater recharge and discharge, numerical groundwater-flow model construction and results, and model limitations. The study area encompasses the western part of Pierce County and the southwestern part of King County, Washington. The study area extends south to the Nisqually River, southwest to Tanwax Creek, northeast to the Green River, and north through the valley near Auburn and adjacent uplands. It is bounded on the east by foothills of the Cascade Range, and on the northwest by Puget Sound.

Washington

Analysis of minimum 7-day discharges and estimation of minimum 7-day, 2-year discharges for streamflow-gaging stations in the Brazos River basin, Texas

The minimum 7-day, 2-year discharge (7Q2) was estimated for each of 97 streamflow-gaging stations in the Brazos River Basin, Texas, with at least 10 years of unregulated or 10 years of regulated discharge data using the Pearson Type III distribution fitted by L-moment statistics of the annual minimum 7-day discharges. Forty-six stations have at least 10 years of unregulated discharge data only, 26 stations have at least 10 years of regulated discharge data only, and 25 stations have at least 10 years of unregulated and 10 years of regulated discharge data. Statistical tests were applied to the annual minimum 7-day discharges to determine if the unregulated and regulated data have trends; and to determine if the data before regulation are significantly different from the data after regulation. Results of the Mann-Kendall test indicate a significant trend at the 95-percent confidence level for 8 of the 46 unregulated-only stations, for 6 of the 26 regulated-only stations, and for 5 stations before regulation and for 10 stations after regulation of the 25 stations with both unregulated and regulated data. Results of the Wilcoxon rank-sum test indicate a significant difference between the unregulated and regulated annual minimum 7-day discharges at the 95-percent confidence level for 16 of the 25 stations with both types of data. The 7Q2 is zero for 31 of the 46 unregulated-only stations and for 9 of the 26 regulated-only stations. The 7Q2 also is zero for 9 stations before regulation and for 3 stations after regulation of the 25 stations with both unregulated and regulated data. The 7Q2 ranges from 33 to 631 cubic feet per second for the eight regulated-only stations on the mainstem Brazos River. For the seven stations on the mainstem Brazos River with at least 10 years of unregulated and 10 years of regulated discharge data, the 7Q2 ranges from 0 to 568 cubic feet per second before regulation and from 0.30 to 670 cubic feet per second after regulation. The 7Q2 increased during regulation for 19 of the 25 stations with both unregulated and regulated data. The effect of regulation for most stations generally was an increase in the magnitude of the 7Q2.

Texas

The geometry of fault reactivation and uplift along the central part of the Maacama fault zone, northern California Coast Ranges (USA)

Fault reactivation of bedrock structures in active fault zones influences stress state and earthquake rupture phenomena through the introduction of weak slip surfaces that impact fault zone geometry and width. Yet, geometric relationships between modern faults and older reactivated faults are difficult to quantify in rocks that have experienced multiple deformation episodes. We used new geologic mapping, geomorphic tools, and structural modeling to quantify rock uplift and subsurface fault geometry of the central part of the Maacama Fault Zone near Ukiah, California, USA, and the surrounding area. Results suggest that the northern Mayacamas Mountains are in a tectonically driven disequilibrium, with differential rock uplift focused on the western side of the range. Steeply east-dipping fault surfaces and splays characterize the geometry of the Maacama Fault Zone. We mapped two newly identified faults to the east of the main Maacama Fault, the Cow Mountain–Mill Creek Fault, and Willow Creek Fault, which align with a moderately east-dipping cluster of microseismicity between 4–10 km depth beneath the Mayacamas Mountains. Static stress modeling on the Maacama Fault Zone and newly identified faults to the east quantify slip tendency values of 0.5–0.4, which suggests that the faults are moderately to poorly suited for slip in the modern stress field and may be weak. We infer that modern uplift is driven by oblique reverse, up-to-the-east, dip-slip motion on the reactivated Cenozoic Cow Mountain–Mill Creek and Willow Creek Faults as material is advected through a restraining bend on the Maacama Fault. This study shows that reactivated bedrock faults increase the fault zone width and introduce fault surfaces that contribute a component of vertical deformation and uplift in major strike-slip fault zones. Deformation is accommodated on an interconnected network of new and reactivated faults that delineate a complex seismic hazard.

California

Noble and base metal distribution and processes affecting ore tenors in the disrupted lower stratigraphy of the Stillwater Complex, USA

Exploration continues for contact-style Ni-Cu sulfide and chromitite-associated PGE mineralization in ultramafic rocks of the Stillwater Complex. At the Iron and Chrome Mountain areas, massive sulfides occur along the complex’s footwall contact and anomalous concentrations of PGE+Au are associated with the three lowermost chromitite seams. Southeast of Chrome Mountain, magmatic layering is highly disrupted by the presence of faults, magmatic breccias, serpentinized discordant dunites, pyroxenite pegmatoids, and disaggregated chromitite seams. The bulk rock chemistry, sulfide chemistry, and noble metal mineralogy of samples from this area were examined to determine the deportment of PGE and processes that led to enrichments in PGE, Au, Cu, Co and Ni. Results show that a sulfide liquid was the principal collector of PGE. If sulfide liquid was initially deposited with chromite, it was disaggregated or redistributed by subsequent melt or fluid infiltration, which may have resulted in the offset of peak PGE(+Cu, Ni) from peak Cr 2 O 3 concentrations, and upgraded PGE tenors. Upon cooling, PGE exsolved from sulfides to form discrete bismuth tellurides, arsenides, arsenic sulfides, antimonides, and alloys, commonly along the margins of sulfide globules. Calculated metal tenors are highest in the disseminated sulfides southwest of Chrome Mountain, whereas massive and net-textured sulfides near the Iron Mountain-Camp zone represent monosulfide solid solution cumulates. At progressively shallower levels, higher metal tenors combined with lower S/Se ratios are consistent with increasing R-factors from 100 to 100,000. Serpentinization and talc-tremolite alteration resulted in S loss through partial replacement of sulfides by secondary silicate+carbonate+magnetite+sulfide assemblages, further upgrading Ni-Cu-PGE tenors. The present work shows that processes responsible for the disruption of magmatic layering and post-magmatic fluid alteration along the intrusion’s lower contact led to noble and base metal enrichments.

Montana

Divisions of geologic time—Major chronostratigraphic and geochronologic units

Introduction The reports and maps of our Nation’s geological surveys inform and benefit the public, private industry, government officials, and scientists. The use of clear and consistent nomenclature and classifications can improve communication of data and interpretations. Since 1899, the U.S. Geological Survey (USGS) Geologic Names Committee (GNC) has been responsible for defining standards that promote uniform geologic nomenclature and classifications among geoscientists. The GNC periodically publishes a geologic time scale, the “Divisions of Geologic Time,” that serves as the national standard for USGS publications (for example, refer to Orndorff and others, 2023). Authors may use other published geologic time scales, such as those of the Geological Society of America (GSA) or the International Commission on Stratigraphy (ICS), provided that they are clearly specified and referenced. Access to the USGS, GSA, and ICS geologic time scales is also available from the U.S. Geologic Names Lexicon (Geolex) website https://ngmdb.usgs.gov/Geolex/stratres/timescales. The geologic time scale serves a dual purpose by enabling authors to distinguish earth material units by position (chronostratigraphic) and time (geochronologic), as outlined in order of decreasing rank.

Fact Sheet

Arctic fold-and-thrust belts

The modern Arctic has been formed through a series of continent–continent collisions, accretion of terranes and phases of crustal extension. The Neoproterozoic Timanian, Paleozoic Caledonian and Uralian, and late Mesozoic Verkhoyansk–Kolyma, Chukotkan and Brookian orogenies formed several large fold-and-thrust belts (FTBs). The FTBs are exposed across vast areas of continents and continue offshore to form a complex tectonic basement for thick sedimentary basins, playing an important role in the history of accumulation and deformation of younger unmetamorphosed sedimentary successions that are the subject of this volume. Recognition of the importance of FTBs in the Arctic geological history and their role as a controlling factor of development of Arctic sedimentary basins resulted in this chapter, in which we review the current state of knowledge about Arctic FTBs and highlight questions that remain to be addressed. Enclosure D , a map showing boundaries of the FTB and their internal first-order structural fabric, is a part of the overview.

Book chapter

Bayesian ETAS modeling for the Pacific Northwest: Uncovering effects of tectonic regimes, regional differences, and swarms on aftershock parameters

The Pacific Northwest (PNW) of North America has high seismic hazard due to numerous earthquake sources under populated areas. It hosts several tectonic regimes and subregional seismic zones that are hypothesized to have different patterns of earthquake and aftershock occurrence. It is also predisposed to earthquake swarms, which can complicate the statistical modeling of these patterns. We present the first statistical seismicity model of the PNW catalog using the epidemic‐type aftershock sequence (ETAS) framework. We develop a Bayesian inference procedure that provides a stable estimation of both ETAS parameters and their uncertainties for different sets of PNW earthquakes, even those with very sparse catalogs. The Bayesian approach allows us to investigate how parameter estimates change between the intraslab and crustal tectonic regimes, the northern and southern PNW, and when swarms are included and excluded from the catalog. We also utilize our Bayesian framework to calculate parameter estimates under different prior beliefs about PNW seismicity, as well as to propagate catalog measurement errors into ETAS parameter estimates. We discuss the implications of parameter differences across the region for aftershock forecasting for the PNW.

Pacific Northwest

Models no not provide proof: An example of model ambiguity and application of isotopic data in a mine pit lake

Geochemical and hydrologic models of pit lakes are commonly used in environmental regulatory decisions to predict future water quality and hydrologic conditions and to understand existing pit lakes. Models may be used to quantify sulfide oxidation, predict thermal/chemical stratification and mixing, and better understand connections between pit lakes and aquifers. One concern related to the hydrologic character of pit lakes is if they are terminal (a groundwater sink with no outflow) or flowthrough (both receiving groundwater inflow and discharging to groundwater). This question was pertinent to the Liberty pit lake, a small acidic pit lake formed in a former Cu deposit in south-central Nevada where potentiometric and geochemical data potentially indicate pit-lake outflow. Potential discharge to groundwater from the pit lake was evaluated using a water-balance model, but uncertainty in hydraulic parameters led to ambiguity in the hydrologic character. Stable isotopes of water were then sampled from the pit lake and adjacent groundwater wells, which unambiguously indicated the lack of an evaporative signature in downgradient groundwater because the groundwater did not plot on a hypothetical mixing line between evaporated pit lake water and observed meteoric recharge. This methodology provided a more effective and more data-driven approach for understanding pit-lake hydrology. Although predictive models are required to quantify reasonable bounds on future conditions, many models contain substantial uncertainty and are not well suited in some environments. Datasets that provide more clear lines of evidence could be collected from existing pit lakes whenever possible to inform water-rock interaction, limnological behavior, and connectivity to adjacent groundwater.

Nevada

Collision structures of the Prince William terrane and Chugach terrane docking along the Shumagin and Unimak convergent margins, Alaska, USA

Western Alaska’s convergent margins are composed of tectonostratigraphic terranes. On land, terrane assembly is recognized along boundaries or sutures between neighboring geologic elements with distinctly different origins. In marine areas where rock outcrops are covered by sediment, recognizing terrane sutures is problematic. A fault in seismic dip line 5 of the ALEUT project has been interpreted as a terrane suture. It is imaged intermittently down to the 30+-km-deep plate interface. Processing of ALEUT strike line 7 revealed the suture at ~18 km depths extending 300 km along the margin. Upper structures in line 5 are like the structures of adjacent seismic transects where imaging is only 8−10 km deep. They were previously not recognized as the upper reaches of terrane sutures and show structural details obscured at greater depths. The composite data are the basis for a simple tectonic model of terrane docking.

Alaska

Effects of septic-tank effluent on ground-water quality in northern Williamson County and southern Davidson County, Tennessee

An investigation of the potential contamination of ground water from septic tank systems blasted in bedrock in Williamson and Davidson Counties, Tennessee, was conducted during 1988-89. Water samples were collected from domestic and observation wells, springs, and surface-water sites in a residential subdivision in the northern part of Williamson County near Nashville. The subdivision has a high density of septic-tank field lines installed into blasted bedrock Water samples also were collected from a well located in an area of Davidson County where field lines were installed in 5 feet of soil. Samples were analyzed for major inorganic constituents, nutrients, total organic carbon, optical brighteners, and bacteria. Although results of analyses of water samples from wells indicate no effect of septic-tank effluent on ground-water quality at these sites, water from two springs located downgradient from the subdivision had slightly larger concentrations of nitrite plus nitrate (2.2 and 2.7 milligrams per liter N), and much larger concentrations of fecal coliform and fecal streptococci bacteria (2,000 to 3,200 and 700 to 900 colonies per 100 milliliters of sample, respectively), than other wells and springs sampled during 1988. Water from one of these springs contained optical brighteners, which indicates that septic-tank effluent is affecting ground-water quality.

Tennessee

Non-destructive sampling of larval amphibians shows tail tissues reflect whole-body methylmercury concentrations and trophic-related differences in bioaccumulation

Mercury (Hg), and specifically methylmercury (MeHg), is a contaminant of global concern to humans and wildlife, but there is still limited understanding of its effects on many taxa. Amphibians are closely associated with water, where Hg is converted to MeHg, and high concentrations of MeHg can reduce survival of amphibians in the wild. However, limited non-lethal proxies exist for estimating MeHg bioaccumulation in amphibian larvae, which could be used to reduce destructive sampling and allow sampling of imperiled species. We evaluated whether tail clips were indicative of whole-body MeHg concentrations of larval amphibians, whether dragonflies were an effective bioindicator of Hg concentrations of larval amphibians, and compared MeHg bioaccumulation in predatory caudate (salamanders, newts) and non-predatory anuran (frogs, toads) larvae. Tail-clip and whole-body MeHg concentrations were strongly correlated, especially for caudates. There was high variability among sites and species, but caudates had higher MeHg than anurans, and MeHg concentrations of caudates and anurans were strongly correlated. Dragonflies were a weak indicator of MeHg for anurans but strong indicators for larval caudates, likely because caudates and dragonflies are carnivorous and occupy similar trophic positions. Our study revealed that sampling tails can be an effective non-destructive index of whole-body MeHg for larval amphibians and demonstrates differences in MeHg bioaccumulation related to trophic position in larval amphibians and dragonflies. All post-embryo life stages for a broad suite of amphibian species can now be sampled non-lethally for MeHg, allowing for a broader examination of population-level effects and focus on imperiled species that cannot be sampled destructively.

conterminous United States

Geologic map and structure sections along the southern part of the Bartlett Springs Fault Zone and adjacent areas from Cache Creek to Lake Berryessa, northern Coast Ranges, California

Introduction Located in the Coast Ranges of northern California, the Bartlett Springs Fault Zone is the easternmost fault in the San Andreas Fault system in northern California. The fault is a right-lateral, strike-slip structure considered capable of producing an earthquake of moment magnitude 7. The purpose of this mapping is to better characterize the geology and earthquake hazards associated with the southern part of the Bartlett Springs Fault Zone and to help identify any evidence of active uplift on the faults bounding the Coast Ranges. Although the area immediately surrounding the Bartlett Springs Fault Zone is sparsely populated, its southern segment presents a potential seismic hazard to northern California communities as far away as the San Francisco Bay region and Sacramento. There are also nearby water resources, mineral resources, and public lands used for public recreation. The Coast Ranges of northern California are a series of northwest-southeast-oriented mountain ranges and valleys located north of the San Francisco Bay region, between the Pacific Ocean to the west and the Sacramento Valley to the east. The region has rugged terrain, high mountain peaks that reach more than 2,400 meters above sea level, isolated and narrow valley bottoms on which most human settlements are located, and large drainage systems that tend to follow the northwest-southeast-oriented topographic grain. The physiographic character of the region is shaped by its bedrock geology, deformational history, and active faulting. The basement rocks of the northern Coast Ranges consist of the Franciscan Complex and the Great Valley complex, the latter of which consists of two informal units, the Coast Range ophiolite and the Great Valley sequence. The Franciscan Complex and the Great Valley complex are in structural contact along the Coast Range Fault, a regional-scale structure and fundamental crustal boundary. The Franciscan Complex and the Great Valley complex are superposed by active, northwest-southeast-striking strike-slip faults that are associated with seismicity swarms. These active strike-slip faults can produce moderate to large earthquakes that have moment magnitudes of 7–8. In places, these active structures bound large ranges and valleys, suggesting that much of the modern topographic expression is the result of active deformation processes. This report contains new 1:24,000-scale geologic mapping along the southern part of the Bartlett Springs Fault Zone between Clear Lake and Lake Berryessa. The map area spans 738 square kilometers in northern Napa County, southern Lake County, and parts of Yolo and Colusa Counties. The south and east borders of the map are 90 kilometers north of San Francisco and 70 kilometers west of Sacramento, respectively. The map area is within the Knoxville mining district, which has a history of mercury and gold mining dating back to the mid-19th century. The two main towns in the region, Lower Lake and Clearlake, California, are west-northwest of the map area. Approximately 71,000 people live in the cities and rural communities located within a 40-kilometer radius of the center of the map area. The bedrock geology, cross sections, and structural data presented herein are critical for evaluating the long-term evolution of the Bartlett Springs Fault Zone. This work will supplement studies on local seismic hazards, liquefaction potential, landslide hazards, earthquake geology, natural resources, groundwater resources, engineering geology, and tectonic history by providing the background information for site-specific investigations on these subjects.

California

Factors contributing to pesticide contamination in riverine systems: The role of wastewater and landscape sources

Wastewater treatment plant (WWTP) discharges can be a source of organic contaminants, including pesticides, to rivers. An integrated model was developed for the Potomac River watershed (PRW) to determine the amount of accumulated wastewater percentage of streamflow (ACCWW) and calculate predicted environmental concentrations (PECs) for 14 pesticides in non-tidal National Hydrography Dataset Plus Version 2.1 stream segments. Predicted environmental concentrations were compared to measured environmental concentrations (MECs) from 32 stream sites that represented a range of ACCWW and land use to evaluate model performance and to assess possible non-WWTP loading sources. Statistical agreement between PECs and MECs was strongest for insecticides, followed by fungicides and herbicides. Principal component analysis utilizing optical fluorescence and ancillary water quality data identified wastewater and urban runoff sources. Pesticides that indicated relatively larger sources from WWTPs included dinotefuran, fipronil, carbendazim, thiabendazole, and prometon whereas imidacloprid, azoxystrobin, propiconazole, tebuconazole, and diuron were more related to urban runoff. In addition, PECs generally comprised a low proportion of MECs, which indicates possible dominant loading sources beyond WWTP discharges. Cumulative potential toxicity was higher for sites with greater ACCWW and/or located in developed areas. Imidacloprid, fipronil, and carbendazim accounted for the largest portion of predicted potential toxicity across sites. The chronic aquatic life toxicity benchmarks for freshwater invertebrates were exceeded for 82 % of the imidacloprid detections ( n = 28) and 47 % of the fipronil detections ( n = 19). These results highlight the ecological implications of pesticide contamination from WWTP discharges and also the potential legacy effects from accumulated soil and groundwater sources. Pesticide management strategies that mitigate both current and historical impacts may improve the health of aquatic ecosystems.

Potomac River watershed

Pre-Acadian tectonics of the eastern Orange-Milford Belt, south-central Connecticut

This excursion presents a reinterpretation of mapping and new analytical data from the eastern Orange-Milford belt (OMB) in south-central Connecticut. The OMB is a fault-bound terrane of argillites and mafic rocks of anomalously low metamorphic grade—and of poorly constrained ages and tectonic affinity—wedged between kyanite/sillimanite-grade peri- Laurentian rocks to the west and anatectic peri-Gondwanan rocks to the east (Fig. 1A). Our data demonstrate that Ordovician(?) igneous and sedimentary rocks of the OMB were variably metamorphosed in the Ordovician and Silurian but escaped regionally pervasive, high-grade Devonian and later metamorphism. Previous interpretations (Fritts 1963a, 1965a, 1965b; Burger, 1967; Burger and others, 1968; Rodgers, 1985) described these rocks as a conformable sequence of low-grade, Ordovician to Devonian metasediments and metavolcanics. Our results reveal that the “metavolcanics” are not extrusive rocks but rather slivers of lower oceanic crust with complicated high- and low-grade metamorphic fabrics, intruded by a swarm of Silurian sheeted basalt dikes, and in fault contact with the surrounding metasediments. The purpose of this trip is to show evidence of early Paleozoic (Taconic) deformation and metamorphism preserved in rocks of the eastern OMB. These rocks remained shallow, cool, and sufficiently dry during the regionally dominant Acadian and Alleghanian orogenies to have avoided significant overprinting. As such, these rocks serve as windows into a geologic history otherwise unavailable between anatectic rocks of the peri-Gondwanan Bronson Hill, Avalon, and Gander terranes east of the Hartford basin and sillimanite-grade rocks of the peri-Laurentian Hartland and gneiss dome belts west of the OMB. We present major and trace element geochemistry including rare-earth element patterns for all mafic units in the eastern OMB as well as 40 Ar/ 39 Ar age spectra of amphibole, muscovite, and K-feldspar from rocks of the Maltby Lakes complex (of Deasy and others, 2017), Savin Schist, and Wepawaug Schist. Our evidence demonstrates that the units of the OMB have been assembled by faulting or intrusion, and that no stratigraphic relationships exist between the argillaceous schists and the metaigneous rocks.

Connecticut

Reconstruction of Holocene and Last Interglacial vegetation dynamics and wildfire activity in Southern Siberia

Wildfires are a rapidly increasing threat to boreal forests. While our understanding of the drivers behind wildfires and their environmental impact is growing, it is mostly limited to the observational period. Here we focus on the boreal forests of southern Siberia and exploit a U–Th-dated stalagmite from Botovskaya Cave, located in the upper Lena region of southern Siberia, to document wildfire activity and vegetation dynamics during parts of two warm periods: the Last Interglacial (LIG; specifically part of the Last Interglacial maximum between 124.1 and 118.8 ka) and the Holocene (10–0 ka). Our record is based on levoglucosan (Lev), a biomarker sensitive to biomass burning, and on lignin oxidation products (LOPs) that discriminate between open and closed forest and hard- or softwood vegetation. In addition, we used carbonate carbon stable isotope ratios ( δ 13 C ), which reflect a dominant control of the host rock, to evaluate soil respiration and local infiltration changes. Our LOP data suggest that, during the Last Interglacial, the region around Botovskaya Cave was characterised by open forest, which by ca. 121.5 ka underwent a transition from fire-resistant hardwood to fire-prone softwood. The Lev record indicates that fire activity was high and increased towards the end of Last Interglacial just before 119 ka. In contrast, the Holocene was characterised by a closed-forest environment with mixed hard- and softwood vegetation. Holocene fire activity varied but at a much lower level than during the Last Interglacial. We attribute the changes in wildfire activity during the intervals of interest to the interplay between vegetation and climate. The open forests of the Last Interglacial were more likely to ignite than their closed Holocene equivalents, and their flammability was aided by warmer and drier summers and a stronger seasonal temperature contrast due to the increase in seasonal insolation difference compared to the Holocene. Our comparison of the last two interglacial intervals suggests that, with increasing global temperatures, the boreal forest of southern Siberia may become progressively more vulnerable to higher wildfire activity.

Botovskaya Cave, Siberia