USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Geophysical Research - Planets”

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 1,045 records · Page 58Linked to original sources

Apparent age dependence of the fault weakening distance in rock friction

During rock friction experiments at large displacement, room temperature and humidity, and following a hold test, the fracture energy increases approximately as the square of the logarithm of hold duration. While it's been long known that failure strength increases with log hold time, here the slip weakening distance, d h , also increases. The weakening distance increase is large, hundreds of percent change over a few thousand seconds. The initial bare surface and simulated fault gouge experiments were conducted in rotary shear at 25 MPa normal stress, 21 MPa confining stress and at displacements greater than 100 mm. In contrast, initially bare surface experiments at 5 MPa normal stress, unconfined at displacements less than 10 mm show effectively no change in d h . We attribute the difference to the presence of an appreciable shear zone that develops due to wear over significant displacements, confined at elevated normal stress. Prior published studies of sheared simulated fault gouge at short displacement show both acknowledged and unacknowledged increases in d h that may relate to our observations. Since natural faults have well-developed shear zones, the observations have more direct relevance to earthquake nucleation than prior laboratory studies that use short displacement data and focus on frictional strength recovery alone. However, the physics underlying this increase in weakening distance are not known; candidates are compaction (Nakatani, 1998) and delocalization (Sleep et al., 2000). Additional caveats are that these are room temperature and humidity experiments, at a single normal stress that have not yet been reproduced in other laboratories.

Journal of Geophysical Research↗

Three Mw ≥ 4.7 earthquakes within the Changning (China) shale gas field ruptured shallow faults intersecting with hydraulic fracturing wells

From 2017 to 2019, three destructive earthquakes (27 January 2017 Mw 4.7, 16 December 2018 Mw 5.2, and 3 January 2019 Mw 4.8) occurred in the Changning shale gas field in the southwest Sichuan Basin, China. Previous seismological studies attributed these events to hydraulic fracturing (HF), but were unable to identify the causative seismogenic faults and their slip behaviors. Here, we use Sentinel-1 synthetic aperture radar data to measure surface deformation triggered by the three events and conduct geodetic inversions to characterize their rupture models. The resulting coseismic interferograms show prominent surface deformation with the maximum line-of-sight displacements of up to 4 cm. The inversion results show that all three earthquakes mainly ruptured sedimentary formations above the shale gas bed, in the upper 3 km of the crust, with slip magnitudes ranging from 8.5 to 15 cm, and stress drops ranging from ∼1.8 to ∼3.3 MPa. Their source faults intersect with horizontal HF wells, but do not root in the crystalline basement. Combined with the reported difficulty of increasing HF operation pressures prior to the three events, we argue that they were most likely induced by direct injection of pressurized fluids into the fault zones. Crustal deformation patterns inferred from regional topography and GPS velocities highlight that the Changning field is located within a triple junction region near the southeastern margin of the Tibetan Plateau with large deformation gradients; such conditions are not only favorable to the development of critically stressed faults, but also facilitate the occurrence of at least moderate magnitude earthquakes.

Changning shale gas field, Sichuan Basin↗

Possible effects of multiphase methane evolution during a glacial cycle on underpressure development in sedimentary basins: An analysis with application to the northeast Michigan Basin

In low-permeability systems, groundwater may be accompanied by separate-phase fluids, and measured pore water pressures may deviate from those expected in steady-state, single-phase systems. These same systems may be of interest for storage of nuclear waste in Deep Geologic Repositories. Therefore, it is important to understand the relationship between the presence of a separate phase and anomalous pressure development. At the Bruce site in Southern Ontario, a significant underpressure was observed, and there is evidence for the presence of gas-phase methane in situ. This study used a one-dimensional (vertical) numerical model of the subsurface down to a depth of 844 m beneath the Bruce site to evaluate possible effects of hydromechanical coupling with multiphase flow on pressure evolution during glacial loading and unloading. The simulated pressure conditions were affected strongly by the amount of methane initially present in the system, and the maximum simulated underpressure varied nonmonotonically with increasing initial methane content. When the initial methane content was below the solubility limit, exsolution led to underpressures that briefly exceeded those that formed in the single-phase case. At intermediate initial methane contents (sufficient to produce an immobile gas phase), the gas phase dampened the hydromechanical effects of the glacial cycle. At large initial methane contents (when a mobile gas phase was present), gas migration caused a large decrease in relative liquid permeability, which further contributed to underpressure development in the pore water. Multiple scenarios that spanned a range of initial methane contents yielded underpressures like those observed at the Bruce site.

Journal of Geophysical Research - Solid Earth↗

Earthquake-derived seismic velocity changes during the 2018 caldera collapse of Kīlauea volcano

The 2018 Kīlauea caldera collapse produced extraordinary sequences of seismicity and deformation, with 62 episodic collapse events which significantly altered the landscape of the summit region. Despite decades of focused scientific studies at Kīlauea, detailed information about the internal structure of the volcano is limited. Recently developed techniques in seismic interferometry can be used to monitor the internal structure of an active volcano more directly by detecting subtle spatiotemporal changes in seismic wave velocity, but their utility relies on accurate interpretations of the underlying phenomena causing those velocity changes. Here, we retrospectively apply repeating-earthquake-based seismic interferometry to the 2018 Kīlauea eruption sequence. We find that seismic velocities changed over two distinct time scales: a sudden increase followed by a slower decrease in velocity in the hours following each collapse event, and a gradual, long-term decrease in velocity over several weeks that ceased approximately 1 month prior to the end of the eruption. Modeling suggests that short-term changes can be explained by magma reservoir pressurization which specifically closed vertical ring fractures. Long-term changes are related to subsidence of the caldera and likely include the influence of inelastic strain from the formation of new fractures. These observations provide new insights into the evolution of Kīlauea during its progressive collapse and will inform future interpretations for near-real-time monitoring at hazardous volcanoes around the world using similar techniques, especially where a dominant fracture orientation is present.

Hawaii↗

Detection of aseismic slip and poroelastic reservoir deformation at the North Brawley Geothermal Field from 2009 to 2019

The North Brawley Geothermal Field, located within the Brawley Seismic Zone of Southern California, presents a case study for understanding seismic hazards linked to fluid injection and geothermal energy extraction. An earthquake swarm near the geothermal field in 2012 included two earthquakes with magnitudes greater than 5 and was potentially preceded by a years-long aseismic slip transient. To better understand ground deformation around the geothermal field, including its evolution with time and its physical mechanisms, we analyze deformation before, during, and after the swarm using ground- and satellite-based geodetic techniques between 2009 and 2019. We integrate observations from GNSS, Sentinel-1, TerraSAR-X, UAVSAR, and leveling surveys into a single deformation history. Modeling of this new collection of observations at the North Brawley Geothermal Field provides evidence for 80% more pre-swarm aseismic slip than previously recognized from 2009 to 2012. During the 2012 Brawley swarm, our geodetic slip inversions closely match the results of seismic waveform inversions from the swarm events. After the 2012 swarm, surface deformation is dominated by poroelastic deformation of a shallow fluid reservoir at <1 km depth rather than fault slip. The deformation history and seismicity catalogs at North Brawley suggest a cessation of fault-related slip during the ∼7 years after the 2012 earthquake swarm.

California↗

Insights on multistage rock avalanche behavior from runout modeling constrained by seismic inversions

Inversion of low-frequency regional seismic records to solve for a time series of bulk forces exerted on the earth by a landslide (a force-time function) is increasingly being used to infer volumes and dynamics of large, highly energetic landslides, such as rock avalanches and flowslides, and to provide calibration information on event dynamics and volumes for numerical landslide runout models. Much of the work to date using landslide runout modeling constrained by seismic data has focused on using single-phase models with frictional or velocity-weakening rheologies. Awareness of multistage landslide initiations is increasing, with discrete failures separated in time contributing to the final impact of an event. Our work utilizes a method for incorporating seismic data as a calibration constraint for landslide runout models, considering variable rheologies and different initiation conditions. This study presents a systematic examination of multiple rheologies and initiation conditions, and shows how these factors affect the force-time function derived from the landslide runout model. Our work confirms that, while rheology and fragmenting or initially coherent initiations affect the force-time function, multiple collapses separated by tens of seconds have the greatest impact on the shape and amplitude. We apply this method to the analysis of three real rock avalanches to better constrain plausible initiation conditions and rheology parameters using both seismic and field data. This study provides insights on how assumptions about the initiation dynamics of the source zone and the runout model definition can aid in the interpretation of seismic inversions for multistage rock avalanches.

Journal of Geophysical Research: Solid Earth↗

Cohesive sediment modeling in a shallow estuary: Model and environmental implications of sediment parameter variation

Numerical models of sediment transport in estuarine systems rely on parameter values that are often poorly constrained and can vary on timescales relevant to model processes. The selection of parameter values can affect the accuracy of model predictions, while environmental variation of these parameters can impact the temporal and spatial ranges of sediment fluxes, erosion, and deposition in the real world. We implemented a numerical model of San Pablo Bay, an embayment within San Francisco Bay, California, for November–December 2014, and compared model outputs to observations of water level, velocity, wave parameters, salinity, and suspended sediment concentration (SSC) in the shallow regions. Idealized model runs show that wind timing relative to the phase of the tides is the strongest control on sediment fluxes and bed erosion. We varied sediment erodibility in the outflow of the Petaluma River; while this causes erosion and deposition to vary strongly through the shallows system, total export from the shallows does not change. Model runs with realistic winds show that wind likely resuspends faster settling particles or allows for more particle flocculation; particle settling velocity controls system-wide sediment accumulation. At the margins of the system, the magnitude of SSC is closely tied to wind direction when winds occur during flood tide, but sediment deposition is less connected: Both bed evolution and SSC need to be considered in the prediction of marsh fate. Spatial patterns of light attenuation due to SSC is strongly tied to assumed settling velocity.

California↗

Distinguishing between regression model fits to global mean sea level reconstructions

Global mean sea level (GMSL) has been rising since the last century, posing a serious challenge for the coastal areas. A variety of regression models have been utilized for determining GMSL rise over the past one hundred years, resulting in a large spread of sea level rise rates and multidecadal variations. In this study, we develop a new nonparametric noise model that is data-dependent and considers overfitting due to regression. The noise model is used to determine whether one regression model has significantly better skill than others over the period 1900–2010. The choices of background noise and GMSL reconstruction influence whether two sea level models can be statistically distinguished. With our new nonparametric noise spectra, the differences of model skills in explaining sea level variance are significant only in 34% of model comparisons. However, stepwise trends with three inflection points are significantly more skillful than the linear, quadratic, or exponential trend for most GMSL reconstructions, suggesting the importance of multidecadal variability of sea level rise in the twentieth century. Nevertheless, stepwise trend models cannot be distinguished from models with a long-term harmonic oscillation, indicating that the shape of multidecadal variability is not conclusive. The multidecadal variability is also significant in the steric and barystatic sea level contributions and is related to both natural and anthropogenic forcings. GMSL predictions based on regression fits in the twentieth century underestimate the sea level rise rate over the period 2011–2020 because the sea level acceleration in the recent decade (2011–2020) is not well represented.

Journal of Geophysical Research- Oceans↗

Hydrate formation on marine seep bubbles and the implications for water column methane dissolution

Methane released from seafloor seeps contributes to a number of benthic, water column, and atmospheric processes. At seafloor seeps within the methane hydrate stability zone, crystalline gas hydrate shells can form on methane bubbles while the bubbles are still in contact with the seafloor or as the bubbles begin ascending through the water column. These shells reduce methane dissolution rates, allowing hydrate-coated bubbles to deliver methane to shallower depths in the water column than hydrate-free bubbles. Here, we analyze seafloor videos from six deepwater seep sites associated with a diverse range of bubble-release processes involving hydrate formation. Bubbles that grow rapidly are often hydrate-free when released from the seafloor. As bubble growth slows and seafloor residence time increases, a hydrate coating can form on the bubble's gas-water interface, fully coating most bubbles within ∼10 s of the onset of hydrate formation at the seafloor. This finding agrees with water-column observations that most bubbles become hydrate-coated after their initial ∼150 cm of rise, which takes about 10 s. Whether a bubble is coated or not at the seafloor affects how much methane a bubble contains and how quickly that methane dissolves during the bubble's rise through the water column. A simplified model shows that, after rising 150 cm above the seafloor, a bubble that grew a hydrate shell before releasing from the seafloor will have ∼5% more methane than a bubble of initial equal volume that did not grow a hydrate shell after it traveled to the same height.

Journal of Geophysical Research - Oceans↗

Postwildfire soil‐hydraulic recovery and the persistence of debris flow hazards

Deadly and destructive debris flows often follow wildfire, but understanding of changes in the hazard potential with time since fire is poor. We develop a simulation‐based framework to quantify changes in the hydrologic triggering conditions for debris flows as postwildfire infiltration properties evolve through time. Our approach produces time‐varying rainfall intensity‐duration thresholds for runoff‐ and infiltration‐generated debris flows with physics‐based hydrologic simulations that are parameterized with widely available hydroclimatic, vegetation reflectance, and soil texture data. When we apply our thresholding protocol to a test case in the San Gabriel Mountains (California, USA), the results are consistent with existing regional empirical thresholds and rainstorms that caused runoff‐ and infiltration‐generated debris flows soon after and three years following a wildfire, respectively. We find that the hydrologic triggering mechanisms for the two observed debris flow types are coupled with the effects of fire on the soil saturated hydraulic conductivity. Specifically, the rainfall intensity needed to generate debris flows via runoff increases with time following wildfire while the rainfall duration needed to produce debris flows via subsurface pore‐water pressures decreases. We also find that variations in soil moisture, rainfall climatology, median grain size, and root reinforcement could impact the median annual probability of postwildfire debris flows. We conclude that a simulation‐based method for calculating rainfall thresholds is a tractable approach to improve situational awareness of debris flow hazard in the years following wildfire. Further development of our framework will be important to quantify postwildfire hazard levels in variable climates, vegetation types, and fire regimes.

Journal of Geophysical Research: Earth Surface↗

Multi-model comparison of computed debris flow runout for the 9 January 2018 Montecito, California post-wildfire event

Hazard assessment for post-wildfire debris flows, which are common in the steep terrain of the western United States, has focused on the susceptibility of upstream basins to generate debris flows. However, reducing public exposure to this hazard also requires an assessment of hazards in downstream areas that might be inundated during debris flow runout. Debris flow runout models are widely available, but their application to hazard assessment for post-wildfire debris flows has not been extensively tested. Necessary inputs to these models include the total volume of the mobilized flow, flow properties (either inherent material properties or calibration coefficients), and site topography. Estimates of volume are possible in post-event (“back calculation”) studies, yet before an event, volume is an uncertain quantity. We simulated debris flow runout for the well-constrained 9 January 2018 Montecito event using three models (RAMMS, FLO2D, and D-Claw) to determine the relative importance of volume and flow properties. We broke the impacted area into three domains, and for each model-domain combination, we performed a numerical sampling study in which volume and flow properties varied within a wide, but plausible range. We assessed model performance based on inundation patterns and peak flow depths. We found all models could simulate the event with comparable results. Simulation performance was most sensitive to flow volume and less sensitive to flow properties. Our results emphasize the importance of reducing uncertainty in pre-event estimates of flow volume for hazard assessment.

California↗

Vadose zone thickness limits pore-fluid pressures and acceleration in a large, slow-moving landslide

The rate and timing of hydrologically forced landslides is a complex function of precipitation patterns, material properties, topography, and groundwater hydrology. In the simplest form, however, slopes fail when subsurface pore pressure grows large enough to exceed the Mohr-Coulomb failure criterion. The capacity for pore pressure rise in a landslide is determined in part by the thickness of the unsaturated zone above the water table, which itself is set by weathering patterns that should have predictable patterns across different lithologies. To investigate how this structure affects landslide behavior, we exploit a multi-year record of precipitation, pore pressure, and velocity from Oak Ridge earthflow, a slow-moving landslide set in Franciscan mélange, northern California, USA. In conjunction with electrical resistivity tomography and hydraulic conductivity measurements, these data show that Oak Ridge has a thin weathered profile that is comparable in thickness to other mélange landslides in California. We propose that due to the inherently thin vadose zone, mélange landscapes experience an unusually high water table that frequently brings them close to movement; however, the capacity to increase stress is limited by the small amount of dynamic storage available. Instead, excess pore pressure is shed via springs and saturation overland flow once the water table reaches the surface. Linkages between weathering patterns, hydrology, and deformation can explain behavior patterns exhibited by Franciscan mélange earthflows across a large precipitation gradient.

California↗

Multi-decadal simulation of marsh topography evolution under sea level rise and episodic sediment loads

Coastal marsh within Mediterranean climate zones is exposed to episodic watershed runoff and sediment loads that occur during storm events. Simulating future marsh accretion under sea level rise calls for attention to: (a) physical processes acting over the time scale of storm events and (b) biophysical processes acting over time scales longer than storm events. Using the upper Newport Bay in Southern California as a case study, we examine the influence of event-scale processes on simulated change in marsh topography by comparing: (a) a biophysical model that integrates with an annual time step and neglects event-scale processes (BP-Annual), (b) a physical model that resolves event-scale processes but neglects biophysical interactions (P-Event), and (c) a biophysical model that resolves event-scale physical processes and biophysical processes at annual and longer time scales (BP-Event). A calibrated BP-Event model shows that large (>20-year return period) episodic storm events are major drivers of marsh accretion, depositing up to 30 cm of sediment in one event. Greater deposition is predicted near fluvial sources and tidal channels and less on marshes further from fluvial sources and tidal channels. In contrast, the BP-Annual model poorly resolves spatial structure in marsh accretion as a consequence of neglecting event-scale processes. Furthermore, the P-Event model significantly overestimates marsh accretion as a consequence of neglecting marsh surface compaction driven by annual scale biophysical processes. Differences between BP-Event and BP-Annual models translate up to 20 cm per century in marsh surface elevation.

California↗

Storm-scale and seasonal dynamics of carbon export from a nested subarctic watershed underlain by permafrost

Subarctic catchments underlain by permafrost sequester a major stock of frozen organic carbon (C), which may be mobilized as the Arctic warms. Warming can impact C export from thawing soils by altering the depth and timing of runoff related to changing storm and fire regimes and altered soil thaw depths. We investigated C export in a first order headwater stream (West Twin Creek) and its receiving third order river (Beaver Creek) in interior Alaska using discrete sampling of dissolved organic and inorganic C (DOC and DIC) and 15-min collection of specific conductance (SC), fluorescent dissolved organic matter (fDOM) and water discharge (Q). Storm SC-Q relationships displayed negative slopes, indicating solute limitation and limited influence of seasonal soil thaw on storm runoff chemistry. Concurrently, fDOM-Q displayed positive slopes that decreased over the summer, indicating flushing of a limited fDOM pool. Baseflow DIC increased over the season concurrent with soil thaw, with higher DIC at the larger scale indicating greater influence of deeper, mineral-rich flow paths. Storm and seasonal trends were generally similar at both scales. The biggest difference was in fDOM, which displayed higher concentrations and slower depletion in the first order stream. Improved process understanding from this study can be used to better predict carbon export and cycling by stream networks as northern forests and arctic regions continue to warm.

Alaska↗

Event scale relationships of DOC and TDN fluxes in throughfall and stemflow diverge from stream exports in a forested catchment

Aquatic fluxes of carbon and nutrients link terrestrial and aquatic ecosystems. Within forests, storm events drive both the delivery of carbon and nitrogen to the forest floor and the export of these solutes from the land via streams. To increase understanding of the relationships between hydrologic event character and the relative fluxes of carbon and nitrogen in throughfall, stemflow and streams, we measured dissolved organic carbon (DOC) and total dissolved nitrogen (TDN) concentrations in each flow path for 23 events in a forested watershed in Vermont, USA. DOC and TDN concentrations increased with streamflow, indicating their export was limited by water transport of catchment stores. DOC and TDN concentrations in throughfall and stemflow decreased exponentially with increasing precipitation, suggesting that precipitation removed a portion of available sources from tree surfaces during the events. DOC and TDN fluxes were estimated for 76 events across a 2-year period. For most events, throughfall and stemflow fluxes greatly exceeded stream fluxes, but the imbalance narrowed for larger storms (>30 mm). The largest 10 stream events exported 40% of all stream event DOC whereas those same 10 events contributed 14% of all throughfall export. Approximately 2–5 times more DOC and TDN was exported from trees during rain events than left the catchment via streams annually. The diverging influence of event size on tree versus stream fluxes has important implications for forested ecosystems as hydrological events increase in intensity and frequency due to climate change.

Vermont↗

Aquatic-terrestrial linkages control metabolism and carbon dynamics in a mid-sized, urban stream influenced by snowmelt

Freshwater streams can exchange nutrients and carbon with the surrounding terrestrial environment through various mechanisms including physical erosion, flooding, leaf drop, and snowmelt. These aquatic-terrestrial interactions are crucial in carbon mobilization, transformation, ecosystem productivity, and have important implications for the role of freshwater ecosystems in the global carbon budget. We utilized high-frequency oxygen, temperature, and carbon dioxide (CO 2 ) data to infer watershed connectivity in Boulder Creek, a mid-sized (1160 km 2 ) watershed located in Colorado, USA. Daily modeled gross primary production (GPP), ecosystem respiration (ER), net ecosystem production (NEP), and reaeration coefficients ( K 600 ) were paired with high-frequency, in-situ dissolved CO 2 data to characterize changes in metabolic regime and carbon flux on a stream influenced by seasonal snowmelt. GPP and ER were correlated ( ρ = −0.72, p ≪ 0.001) during the non-snowmelt period and NEP was frequently negative. Mean F CO2 during the non-snowmelt period was approximately 302 (±171) mmol C m −2 d −1 and was primarily supported by watershed CO 2 inputs. During snowmelt, GPP and ER were not significantly correlated ( ρ = −0.22, p = 0.05), and mean NEP was significantly more negative than during non-snowmelt. Watershed connectivity was higher during snowmelt, as evidenced by significantly higher F CO2 (843 ± 338 mmol C m −2 d −1 ) and greater allochthonous CO 2 inputs than during non-snowmelt periods, emphasizing the effects of seasonal differences in aquatic-terrestrial linkages in this stream. We suggest that our understanding of watershed carbon budgets is subject to temporal dynamics which control the degree of connectivity between terrestrial and aquatic ecosystems.

Colorado↗

Spatiotemporal methane emission from global reservoirs

Inland aquatic systems, such as reservoirs, contribute substantially to global methane (CH 4 ) emissions; yet are among the most uncertain components of the total CH 4 budget. Reservoirs have received recent attention as they may generate high CH 4 fluxes. Improved quantification of these CH 4 fluxes, particularly their spatiotemporal distribution, is key to realistically incorporating them in CH 4 modeling and budget studies. Here we report on a new global, gridded (0.25° lat × 0.25° lon) study of reservoir CH 4 emissions, accounting for new knowledge regarding reservoir areal extent and distribution, and spatiotemporal emission patterns influenced by diurnal variability, temperature-dependent seasonality, satellite-derived freeze-thaw dynamics, and eco-climatic zone. The results of this new data set comprise daily CH 4 emissions throughout the full annual cycle and show that reservoirs cover 297 × 10 3 km 2 globally and emit 10.1 Tg CH 4 yr −1 (1σ uncertainty range of 7.2–12.9 Tg CH 4 yr −1 ) from diffusive (1.2 Tg CH 4 yr −1 ) and ebullitive (8.9 Tg CH 4 yr −1 ) emission pathways. This analysis of reservoir CH 4 emission addresses multiple gaps and uncertainties in previous studies and represents an important contribution to studies of the global CH 4 budget. The new data sets and methodologies from this study provide a framework to better understand and model the current and future role of reservoirs in the global CH 4 budget and to guide efforts to mitigate reservoir-related CH 4 emissions.

Journal of Geophysical Research: Biogeosciences↗

Loss of ice cover, shifting phenology, and more extreme events in Northern Hemisphere lakes

Long-term lake ice phenological records from around the Northern Hemisphere provide unique sensitive indicators of climatic variations, even prior to the existence of physical meteorological measurement stations. Here, we updated ice phenology records for 60 lakes with time-series ranging from 107–204 years to provide the first re-assessment of Northern Hemispheric ice trends since 2004 by adding 15 additional years of ice phenology records and 40 lakes to our study. We found that, on average, ice-on was 11.0 days later, ice-off was 6.8 days earlier, and ice duration was 17.0 days shorter per century over the entire record for each lake. Trends in ice-on and ice duration were six times faster in the last 25-year period (1992–2016) than previous quarter centuries. More extreme events in recent decades, including late ice-on, early ice-off, shorter periods of ice cover, or no ice cover at all, contribute to the increasing rate of lake ice loss. Reductions in greenhouse gas emissions could limit increases in air temperature and abate losses in lake ice cover that would subsequently limit ecological, cultural, and socioeconomic consequences, such as increased evaporation rates, warmer water temperatures, degraded water quality, and the formation of toxic algal blooms.

Journal of Geophysical Research: Biogeosciences↗