Search USGS⌕ Search

SEARCH · Search USGS

Results for “Science Trends”

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.

At least 595 records · Page 33Linked to original sources

Mercury source changes and food web shifts alter contamination signatures of predatory fish from Lake Michigan

To understand the impact reduced mercury (Hg) loading and invasive species have had on methylmercury bioaccumulation in predator fish of Lake Michigan, we reconstructed bioaccumulation trends from a fish archive (1978 to 2012). By measuring fish Hg stable isotope ratios, we related temporal changes in Hg concentrations to varying Hg sources. Additionally, dietary tracers were necessary to identify food web influences. Through combined Hg, C, and N stable isotopic analyses, we were able to differentiate between a shift in Hg sources to fish and periods when energetic transitions (from dreissenid mussels) led to the assimilation of contrasting Hg pools (2000 to present). In the late 1980s, lake trout δ 202 Hg increased (0.4‰) from regulatory reductions in regional Hg emissions. After 2000, C and N isotopes ratios revealed altered food web pathways, resulting in a benthic energetic shift and changes to Hg bioaccumulation. Continued increases in δ 202 Hg indicate fish are responding to several United States mercury emission mitigation strategies that were initiated circa 1990 and continued through the 2011 promulgation of the Mercury and Air Toxics Standards rule. Unlike archives of sediments, this fish archive tracks Hg sources susceptible to bioaccumulation in Great Lakes fisheries. Analysis reveals that trends in fish Hg concentrations can be substantially affected by shifts in trophic structure and dietary preferences initiated by invasive species in the Great Lakes. This does not diminish the benefits of declining emissions over this period, as fish Hg concentrations would have been higher without these actions.

Lake Michigan↗

Impact of precipitation and increasing temperatures on drought trends in eastern Africa

In eastern Africa droughts can cause crop failure and lead to food insecurity. With increasing temperatures, there is an a priori assumption that droughts are becoming more severe. However, the link between droughts and climate change is not sufficiently understood. Here we investigate trends in long-term agricultural drought and the influence of increasing temperatures and precipitation deficits. Using a combination of models and observational datasets, we studied trends, spanning the period from 1900 (to approximate pre-industrial conditions) to 2018, for six regions in eastern Africa in four drought-related annually averaged variables: soil moisture, precipitation, temperature, and evaporative demand ( E 0 ). In standardized soil moisture data, we found no discernible trends. The strongest influence on soil moisture variability was from precipitation, especially in the drier or water-limited study regions; temperature and E 0 did not demonstrate strong relations to soil moisture. However, the error margins on precipitation trend estimates are large and no clear trend is evident, whereas significant positive trends were observed in local temperatures. The trends in E 0 are predominantly positive, but we do not find strong relations between E 0 and soil moisture trends. Nevertheless, the E 0 trend results can still be of interest for irrigation purposes because it is E 0 that determines the maximum evaporation rate. We conclude that until now the impact of increasing local temperatures on agricultural drought in eastern Africa is limited and we recommend that any soil moisture analysis be supplemented by an analysis of precipitation deficit.

Earth Systems Science Dynamics↗

Reach-scale effects of riparian forest cover on urban stream ecosystems

We compared habitat and biota between paired open and forested reaches within five small streams (basin area 10?20 km2) in suburban catchments (9%?49% urban land cover) in the Piedmont of Georgia, USA. Stream reaches with open canopies were narrower than forested reaches (4.1 versus 5.0 m, respectively). There were no differences in habitat diversity (variation in velocity, depth, or bed particle size) between open and forested reaches. However, absence of local forest cover corresponded to decreased large wood and increased algal chlorophyll a standing crop biomass. These differences in basal food resources translated into higher densities of fishes in open (9.0 individuals?m?2) versus forested (4.9 individuals?m?2) reaches, primarily attributed to higher densities of the herbivore Campostoma oligolepis. Densities of terrestrial invertebrate inputs were higher in open reaches; however, trends suggested higher biomass of terrestrial inputs in forested reaches and a corresponding higher density of terrestrial prey consumed by water column feeding fishes. Reach-scale biotic integrity (macroinvertebrates, salamanders, and fishes) was largely unaffected by differences in canopy cover. In urbanizing areas where catchment land cover drives habitat and biotic quality, management practices that rely exclusively on forested riparian areas for stream protection are unlikely to be effective at maintaining ecosystem integrity.

Canadian Journal of Fisheries and Aquatic Sciences↗

Patterns, pace and processes of water-quality variability: Examples from a long-studied estuary

Environmental time series have rich information content that is invaluable for measuring and understanding changes over time and guiding policies to manage change. I extracted information from measurements of 10 water‐quality constituents in upper San Francisco Bay from 1975 to 2016, one of the longest observational records in a U.S. estuary. Changes were detected at every time scale captured by monthly sampling. Long‐term trends included increased ammonium (+53%), nitrate + nitrate (+50%), silicate (+14%), Secchi depth (+42%), and decreased chlorophyll a (Chl a ) (−74%) and suspended particulate matter (−45%). Changes at the decadal scale included abrupt shifts (Chl a , nitrate + nitrite) and oscillations between shorter trends of increase and decrease (Secchi depth, phosphate). Long‐term trends were not expressed equally across all seasons, and seasonal patterns of change varied across constituents. These examples illustrate key features of environmental variability at the land–sea interface: (1) water‐quality components change continually at time scales from months to decades; (2) patterns of seasonal, multiyear, and multidecadal change are complex and vary across constituents; (3) primary drivers of change are freshwater inflow, the master regulator of estuarine dynamics, and human activities such as river damming, water diversions, wastewater discharge, environmental policies, and species introductions; (4) extracting the full information content of time series requires multiple analyses, each revealing a different layer of insight into how changes develop over time; (5) water‐quality variability is nonstationary, so future changes cannot be forecast reliably; (6) repeated observation is an essential method of Earth system science with applications in the design and performance measures of environmental policies.

California↗

Long Term Resource Monitoring procedures—Aquatic vegetation monitoring

This standard operating procedure (SOP) manual describes the collection of standardized, long-term data for aquatic vegetation communities in selected study pools of the Upper Mississippi River System in the United States. The primary intent of the data collection is to assess the status and trends that aid in understanding the unique river ecosystem and to guide large-scale ecological restoration of the river and its biological communities, like aquatic plants and their dependent wildlife. This SOP is an update to the version published in 2000 and reflects modifications to sample sizes and additions of new data collection procedures. All long-term monitoring programs and their SOPs must be adapted to changing conditions and be improved through learning, and this SOP clarifies procedures and adds new elements since the initial SOP was written more than 25 years ago. The SOP is intended for multiple audiences, including vegetation specialists through the Upper Mississippi River Restoration Program, data analysts using the publicly available data generated through this SOP, and natural resource managers and restoration practitioners who need data and science to guide some decisions. This SOP may be transferable and adaptable to other ecosystems when the aquatic plant community is the focus.

Illinois, Iowa, Minnesota, Missouri, Wisconsin↗

Behavior of potentially toxic elements from stoker-boiler fly ash in Interior Alaska: Paired batch leaching and solid-phase characterization

Despite significant investigation of fly ash spills and mineralogical controls on the release of potentially toxic elements (PTEs) from fly ash, interactions with the surficial environment remain relatively poorly understood. We conducted 90-day batch leaching studies with paired analysis of supernatant and solid-phase mineralogy to assess the elemental release and transformation of fly ash upon reaction with aquatic media (18 MΩ cm −1 water and simulated rainwater). The fly ash in this study, collected from the University of Alaska Fairbanks stoker-boiler power plant, is high in unburned carbon (~20% LOI) and highly enriched in several PTEs relative to the upper continental crust. Supernatant concentrations of oxyanion-forming elements (e.g., As, Se, Mo, Sb) remained relatively low and constant, suggesting equilibrium with the solid phase, possibly ettringite [Ca 6 Al 2 (SO 4 ) 3 (OH) 12 •26H 2 O], which is known to incorporate and sorb oxyanion-forming PTEs and was identified by X-ray diffraction. Synthetic precipitation leaching procedure (SPLP) results failed to capture important temporal trends. Lead and Ba supernatant concentrations consistently exceeded drinking water standards, as well as others upon exposure to simulated physiological solutions. Seven-day experiments with dissolved organic matter-isolate solutions indicated that for certain elements, liberation was influenced by carbon concentration and/or the identity of the isolate. Overall, this paired approach can serve as a model for future studies, bridging existing gaps between batch leaching and single-element mineralogical, sorption, or speciation studies.

Alaska↗

Variability and trends in global drought

Monthly precipitation ( P ) and potential evapotranspiration (PET) from the CRUTS3.1 data set are used to compute monthly P minus PET (PMPE) for the land areas of the globe. The percent of the global land area with annual sums of PMPE less than zero are used as an index of global drought (%drought) for 1901 through 2009. Results indicate that for the past century %drought has not changed, even though global PET and temperature ( T ) have increased. Although annual global PET and T have increased, annual global P also has increased and has mitigated the effects of increased PET on %drought.

Earth and Space Science↗

Inference of pattern-based geological CO2 sequestration and oil recovery potential in a commingled main pay and residual oil zone CO2-EOR flood

Several detailed studies have shown that residual oil zones (ROZs) can present significant resources for additional hydrocarbon recovery as well as subsurface carbon dioxide (CO 2 ) sequestration via enhanced oil recovery by injecting CO 2 (CO 2 -EOR). Field development strategies included new wells drilled dedicated to main pay zones (MPZ) and ROZs, or existing wells in MPZs deepened to ROZs for commingled injection-production using different well patterns. The latter presented a challenge when discerning the injection and production from each of the zones, and for subsequent quantification of CO 2 sequestration and EOR potential from different patterns and from the field. In this paper, an innovative method for analyzing commingled injections and productions from MPZs and ROZs, with application to pattern-based data from four staggered line drive patterns in Wasson Field's Denver Unit, Texas, USA, was developed. Decline curve and ratio-trend methods were used as means of history-matching and forecasting. Cumulative production-time and cumulative production-rate data for oil, gas, and water, as well as water-oil ratio (WOR) and gas-oil ratio (GOR), were analyzed along with injection data for time intervals covering major injection events in MPZ, or MPZ and ROZ combined. A combined analysis enabled inference of allocation of fluids into different zones during WAG (water alternating gas) injection and thereby estimation of CO 2 storage, utilization, and retention in different zones as a function of total injection. Results show that ROZs generally present higher CO 2 sequestration potential compared to MPZs, and a comparable incremental oil recovery factor of ∼20%, on average. Results based on ratio analysis further show that while the WOR trend of the pattern production is mostly dominated and controlled by ROZ, GOR is controlled by both intervals. Although the method relying on decline curves and the approach used in zonal fluid allocations are subject to their limitations, this study presents a practical and innovative well-pattern-based method to infer and forecast CO 2 sequestration and oil recovery quantities and fluid ratios from MPZs and ROZs in commingled operations and highlight the added potential offered by ROZs.

Texas↗

Estimated variability of National Atmospheric Deposition Program/Mercury Deposition Network measurements using collocated samplers

The National Atmospheric Deposition Program/Mercury Deposition Network (MDN) provides long-term, quality-assured records of mercury in wet deposition in the USA and Canada. Interpretation of spatial and temporal trends in the MDN data requires quantification of the variability of the MDN measurements. Variability is quantified for MDN data from collocated samplers at MDN sites in two states, one in Illinois and one in Washington. Median absolute differences in the collocated sampler data for total mercury concentration are approximately 11% of the median mercury concentration for all valid 1999-2004 MDN data. Median absolute differences are between 3.0% and 14% of the median MDN value for collector catch (sample volume) and between 6.0% and 15% of the median MDN value for mercury wet deposition. The overall measurement errors are sufficiently low to resolve between NADP/MDN measurements by ??2 ng??l-1 and ??2 ????m-2?? year-1, which are the contour intervals used to display the data on NADP isopleths maps for concentration and deposition, respectively. ?? Springer Science+Business Media B.V. 2007.

Environmental Monitoring and Assessment↗

Reconstruction of historic sea ice conditions in a sub-Arctic lagoon

Historical sea ice conditions were reconstructed for Izembek Lagoon, Bering Sea, Alaska. This lagoon is a crucial staging area during migration for numerous species of avian migrants and a major eelgrass ( Zostera marina ) area important to a variety of marine and terrestrial organisms, especially Pacific Flyway black brant geese ( Branta bernicla nigricans ). Ice cover is a common feature of the lagoon in winter, but appears to be declining, which has implications for eelgrass distribution and abundance, and its use by wildlife. We evaluated ice conditions from a model based on degree days, calibrated to satellite observations, to estimate distribution and long-term trends in ice conditions in Izembek Lagoon. Model results compared favorably with ground observations and 26 years of satellite data, allowing ice conditions to be reconstructed back to 1943. Specifically, periods of significant (limited access to eelgrass areas) and severe (almost complete ice coverage of the lagoon) ice conditions could be identified. The number of days of severe ice within a single season ranged from 0 (e.g., 2001) to ≥ 67 (e.g., 2000). We detected a slight long-term negative trend in ice conditions, superimposed on high inter-annual variability in seasonal aggregate ice conditions. Based on reconstructed ice conditions, the seasonally cumulative number of significant or severe ice days correlated linearly with mean air temperature from January until March. Further, air temperature at Izembek Lagoon was correlated with wind direction, suggesting that ice conditions in Izembek Lagoon were associated with synoptic-scale weather patterns. Methods employed in this analysis may be transferable to other coastal locations in the Arctic.

Alaska↗

Storage capacity and sedimentation trends of Lago Garzas, Puerto Rico, 1996-2007

Lago Garzas is located in west-central Puerto Rico, about 3.5 kilometers southwest of the town of Adjuntas, in the confluence of the Río Vacas and three other unnamed tributaries (fig. 1). The dam is owned and operated by the Puerto Rico Electric Power Authority (PREPA), and was constructed in 1943 for hydroelectric power generation and municipal water use along the southern coast. The dam is a semi-hydraulic earthfill embankment lined with boulders, and has a height of 61.57 meters, a top width of 9.14 meters, a base width of 365.76 meters, and a crest length of 227.37 meters; State Road PR-518 crosses the top of the dam. A morning-glory-type spillway is located near the west abutment of the dam at an elevation of 736.12 meters above mean sea level (Puerto Rico Water Resources Authority, 1969). Figure 2 shows an aerial photograph of the Lago Garzas earthfill dam and the morning-glory spillway section. Additional information and operational procedures are provided in Soler-López and others (1999). During July 17-18, 2007, the U.S. Geological Survey (USGS) Caribbean Water Science Center, in cooperation with the Puerto Rico Aqueduct and Sewer Authority, conducted a bathymetric survey of Lago Garzas to update the reservoir storage capacity and update the reservoir sedimentation rate by comparing the 2007 data with the previous 1996 bathymetric survey results. The purpose of this report is to describe and document the USGS sedimentation survey conducted at Lago Garzas during July 2007, including the methods used to update the reservoir storage capacity, sedimentation rates, and areas of substantial sediment accumulation since 1996.

Lago Garzas Basin Area↗

Hydrologic connectivity in floodplain systems: A multiscale review of concepts, metrics and management

Hydrologic connectivity (HC), particularly in floodplain systems, is pivotal in regulating ecosystem services by facilitating the movement of nutrients, sediments, chemicals, and biota. However, human interventions such as dam construction, levee installation, water management practices, and alterations in vegetation have significantly disrupted natural HC patterns globally. To provide a structured entry into the growing body of HC research, we conducted a systematic literature review of 1920 studies, analysing diverse definitions, influencing factors, quantification approaches, spatial and temporal scales, and management strategies. In addition to traditional review methods, our approach integrates keyword and cluster analysis to elucidate dominant research themes and trends across the literature. Our review reveals that the literature is heavily skewed towards research in North America and Europe (accounting for 72% of studies) and predominantly utilises field investigations, simulation modelling, and remote sensing integrated with geographic information systems. Although these methodologies have advanced our understanding, most studies focus on restricted spatial scales such as individual hillslopes, catchments, or stream networks and short temporal intervals, including single precipitation events or seasonal cycles. A narrow focus becomes a limitation when such studies do not contribute to broader efforts aimed at scaling insights across larger domains. These limitations highlight the potential benefits of innovative conceptual frameworks and quantification methods to better capture HC across broader environments and extended temporal scales. We conclude by discussing challenges in defining and quantifying floodplain HC and outlining potential future research directions to advance connectivity science and management, particularly in floodplain systems characterised by frequent hydrologic fluctuations, such as seasonal inundation and changing flow paths.

Hydrological Processes↗

Three-dimensional imaging, change detection, and stability assessment during the centerline trench levee seepage experiment using terrestrial light detection and ranging technology, Twitchell Island, California, 2012

A full scale field seepage test was conducted on a north-south trending levee segment of a now bypassed old meander belt on Twitchell Island, California, to understand the effects of live and decaying root systems on levee seepage and slope stability. The field test in May 2012 was centered on a north-south trench with two segments: a shorter control segment and a longer seepage test segment. The complete length of the trench area measured 40.4 meters (m) near the levee centerline with mature trees located on the waterside and landside of the levee flanks. The levee was instrumented with piezometers and tensiometers to measure positive and negative porewater pressures across the levee after the trench was flooded with water and held at a constant hydraulic head during the seepage test—the results from this component of the experiment are not discussed in this report. We collected more than one billion three-dimensional light detection and ranging (lidar) data points before, during, and after the centerline seepage test to assess centimeter-scale stability of the two trees and the levee crown. During the seepage test, the waterside tree toppled (rotated 20.7 degrees) into the water. The landside tree rotated away from the levee by 5 centimeters (cm) at a height of 2 m on the tree. The paved surface of the levee crown had three regions that showed subsidence on the waterside of the trench—discussed as the northern, central, and southern features. The northern feature is an elongate region that subsided 2.1 cm over an area with an average width of 1.35 m that extends 15.8 m parallel to the trench from the northern end of the trench to just north of the trench midpoint, and is associated with a crack 1 cm in height that formed during the seepage test on the trench wall. The central subsidence feature is a semicircular region on the waterside of the trench that subsided by as much as 6.2 cm over an area 3.4 m wide and 11.2 m long. The southern feature is an elongate region that has a maximum subsidence of 3.5 cm over an area 0.75 m wide and 8.1 m long and is associated with a number of small fractures in the pavement that are predominately north-south-trending and parallel to the trench. We determined that there was no significant motion of the levee flank during the last week of the seepage test. We also determined biomorphic parameters for the landside tree, such as the 3D positioning on the levee, tree height, levee parallel/perpendicular cross sectional area, and canopy centroid. These biomorphic parameters were requested to support a University of California Berkeley team studying seepage and stability on the levee. A gridded, 2-cm bare-earth digital elevation model of the levee crown and the landside levee flank from the final terrestrial lidar (T-Lidar) survey provided detailed topographic data for future assessment. Because the T-Lidar was not integrated into the project design, other than an initial courtesy dataset to help characterize the levee surface, our ability to contribute to the overall science goals of the seepage test was limited. Therefore, our analysis focused on developing data collection and processing methodology necessary to align ultra high-resolution T-Lidar data (with an average spot spacing 2–3 millimeters on the levee crown) from several instrument setup locations to detect, measure, and characterize dynamic centimeter-scale deformation and surface changes during the seepage test.

California↗

Seismic-reflection profiles of the central part of the Clarendon-Linden fault system of western New York in relation to regional seismicity

Geological and geophysical research in upstate New York, with few exceptions, has not definitively associated seismicity with specific Proterozoic basement or Paleozoic bedrock structures. The central part of the Clarendon-Linden fault system (CLFS) between Batavia and Dale, NY is one of those exceptions where seismicity has been studied and has been spatially associated with structure. The CLFS is either a complex system of long faults with associated shorter branches and parallel segments, or a region of many short faults aligned north-south from the Lake Ontario shore southward to Allegany County, NY. Interpretation of 38 km of Vibroseis and approximately 56 km of conventional seismic-reflection data along 13 lines suggests that the CLFS is a broad zone of small faults with small displacements in the lower Paleozoic bedrock section that is at least 77 km long and 7-17 km wide and spatially coincident with a northtrending geophysical (combined aeromagnetic and gravity) lineament within the basement. The relative offset across the faults of the system is more than 91 m near Attica, NY. The CLFS is the expression of tectonic crustal adjustments within the Paleozoic rock above the boundary of two basement megablocks of differing petrologic provinces and differing earthquake characteristics that forms the eastern side of the Elzevir-Frontenac boundary zone. Deep seismic-reflection profiles display concave-eastward listric faults that probably merge at depth near the mid-crustal boundary layer. An interpretive vertical section provides the setting for refined definitions of the CLFS, its extensions at depth and its relation to seismicity. Most modern seismicity in western New York and the Niagara Peninsula of Ontario occurs in apparent patterns of randomly dispersed activity. The sole exception is a line of seven epicenters of small earthquakes that trend east from Attica, NY into the Rochester basement megablock. Earthquakes may be triggered at the intersections of north- and east-trending brittle faults within the Niagara basement megablock. Current interpretations of the mechanisms for earthquake generation in western New York and the Niagara Peninsula of Ontario require conservative estimates of seismic hazards that assume that an earthquake the size of the 1929 Attica, NY, event (Mb = 5.2) or larger could occur anywhere in the Eastern Great Lakes Basin (EGLB). The broad zone of small-displacement faults that marks the CLFS in the lower Paleozoic sedimentary section and the uppermost basement may not provide the structural environment for generation of earthquakes in western New York. If this interpretation is correct, most seismicity is generated within the Niagara basement megablock beneath or west of the CLFS. Consequently, we may have to look to the deeper tectonic regime of basement megablocks to understand the distribution of modern seismicity in the EGLB. ?? 2002 Elsevier Science B.V. All rights reserved.

Tectonophysics↗

Environmental policy analysis, peer reviewed: Reservoir sediment cores show US lead declines

As a result of the Clean Air Act, lead (Pb) emissions to the atmosphere have been greatly reduced since the mid-1970s. As part of its National Water Quality Assessment, the U.S. Geological Survey has been using paleolimnological techniques to assess past trends in hydrophobic contaminants. In urban-suburban environments, reservoir sediment cores show prominent peaks in Pb distributions that correlate well with the rise and fall of leaded gasoline. However, Pb concentrations in sediments are approximately double those of baseline values prior to the 1950s and 1960s. It is apparent that significant concentrations of anthropogenic Pb still exist in soils and aquatic sediments and that it will take many years to reduce these concentrations to prepollution values, even if there are no new sources of Pb pollution.

Environmental Science & Technology↗

Nimbus IV view of the major structural features of Alaska

Notwithstanding the relatively low degree of ground resolution, many of the major structural features of Alaska can be identified on the Nimbus IV IDCS image, exposed at an altitude of 600 nautical miles (1100 km). In addition, linears of regional extent that may be structurally controlled can be seen, many of which have not yet been recognized in surface mapping. The synoptic view provided by the image brings into focus an orthogonal set of fractures trending north-northeast and east-southeast and not heretofore apparent in regional maps of Alaska. This orthogonal fracture set may reflect a conjugate set of fractures within the crust, which has exerted significant control over the geologic history of the state. Increased resolution in other images from space platforms, such as the resolution of 200 to 650 feet (60 to 200 m) planned for the satellite television cameras of the ERTS program, will permit the discernment of finer detail and a greater accuracy in identifying and locating geologic features.

Alaska↗

Performance of temperature and dissolved oxygen criteria to predict habitat use by lake trout (Salvelinus namaycush)

We compared theoretical habitat volumes, determined from traditional combinations of temperature and dissolved oxygen concentration (DO) boundaries, with in situ habitat use by acoustically tagged lake trout ( Salvelinus namaycush ). The widely used criteria of 8–12 °C underestimated lake trout habitat use by 68%–80%. Instead, combined temperature (<12 or 15 °C) and DO (>4 or 6 mg·L –1 ) criteria most closely matched lake trout habitat use, had a similar seasonal trend as the tagged fish, suggested modest reductions (5% of total lake volume) in habitat during a warmer year, and performed best when the constraints of temperature and DO were most limiting. All data were collected in a small boreal shield lake (27 ha, z max = 21 m) at the Experimental Lakes Area in northwestern Ontario, Canada, during two contrasting periods of thermal stratification (2003: warmer and longer; 2004: cooler and shorter), providing an assessment of observed and theoretical habitat volumes over current environmental extremes.

Canadian Journal of Fisheries and Aquatic Sciences↗

Holocene thermokarst lake dynamics in northern Interior Alaska: The interplay of climate, fire, and subsurface hydrology

The current state of permafrost in Alaska and meaningful expectations for its future evolution are informed by long-term perspectives of previous permafrost degradation. Thermokarst processes in permafrost landscapes often lead to widespread lake formation and the spatial and temporal evolution of thermokarst lake landscapes reflects the combined effects of climate, ground conditions, vegetation, and fire. This study provides detailed analyses of thermokarst lake sediments of Holocene age from the southern loess uplands of the Yukon Flats; including bathymetry and sediment core analyses across a water depth transect. The sediment core results, dated by radiocarbon and 210Pb, indicate the onset of finely laminated lacustrine sedimentation between ~10,000 and 9,000 cal yr BP following basin development through inferred thermokarst processes. Thermokarst expansion to modern shoreline configurations continued until ~5000 cal yr BP, which may have been influenced by increased fire. Between ~5000 and 2000 cal yr BP, the preservation of fine laminations at intermediate and deep-water depths indicate higher lake levels than present. At that time, the lake likely overflowed into an over-deepened gully system that is no longer occupied by perennial streams. By ~2000 cal yr BP, massive sedimentation at intermediate water depths indicates that lake levels lowered, which is interpreted to reflect a response to drier conditions based on correspondence with Yukon Flats regional fire and local paleoclimate reconstructions. Consideration of additional contributing mechanisms include the possible influence of catastrophic lake drainages on downgradient base flow levels that may have enhanced subsurface water loss, although this mechanism is untested. The overall consistency between the millennial lake level trends documented here with regional paleoclimate trends indicates that after lakes formed, their size and depth has likely been affected directly by North Pacific atmospheric circulation changes and indirectly through evolution of permafrost, ground ice and sub-surface hydrology. As the first detailed study of Holocene thermokarst basin expansion, stabilization and subsequent climate-driven lake level variations in a loess upland, results provide a framework for future investigations of paleoclimatic signals from similar lake systems that characterize large regions of Alaska and Siberia.

Alaska↗