Search USGSSearch

SEARCH · Search USGS

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 811 records · Page 45Linked to original sources

Separating sea and slow slip signals on the seafloor

Seafloor pressure measurements hold promise for estimating vertical displacements from transient slow slip events on submarine faults. We assess the accuracy of pressure offset estimates that evolve over days to weeks and the confidence with which they may be attributed to tectonic deformation or to the ocean water column. One common approach to resolve this ambiguity assumes water column pressures vary insignificantly over the study region and are represented by stable reference site pressures. Assessing the validity of this assumption requires independent evidence. Correlations between pressures and colocated temperatures collected during the Hikurangi Ocean Bottom Investigation of Tremor and Slow Slip experiment suggest temperatures might provide a useful independent proxy for water column pressures. We compared offsets estimated using several methods, with temperature and other proxies. The use of a temperature proxy was unsuccessful, because seafloor temperatures did not track the seasonal signal that contributes significantly to seafloor pressure changes over the slow slip event period. Regardless of the estimation method, offsets varied within a few cm around some uncertain reference level. Commonly used statistical measures are shown not to be reliable indicators of offset accuracy since offsets contribute minimally to the total variance. Offsets estimated using identical methods but with seafloor pressures simulated using a regional ocean model were larger than those derived from the data but had a similar pattern. Since the model simulates only water column processes, this suggests a significant fraction of the estimated pressure offsets are due to seasonal water column signal and are not of tectonic origin.

Journal of Geophysical Research- Solid Earth

The effect of brine on the electrical properties of methane hydrate

Gas hydrates possess lower electrical conductivity (inverse of resistivity) than either seawater or ice, but higher than clastic silts and sands, such that electromagnetic methods can be employed to help identify their natural formation in marine and permafrost environments. Controlled laboratory studies offer a means to isolate and quantify the effects of changing individual components within gas‐hydrate‐bearing systems, in turn yielding insight into the behavior of natural systems. Here we investigate the electrical properties of polycrystalline methane hydrate with ≥25% gas‐filled porosity and in mixture with brine. Initially, pure methane hydrate was synthesized from H 2 O ice and CH 4 gas while undergoing electrical impedance measurement, then partially dissociated to assess the effects of pure pore water accumulation on electrical conductivity. Methane hydrate + brine mixtures were then formed by either adding NaCl (0.25–2.5 wt %) to high‐purity ice or by using frozen seawater as a reactant. Conductivity was obtained from impedance measurements made in situ throughout synthesis while temperature cycled between +15 °C and −25 °C. Several possible conduction mechanisms were subsequently determined using equivalent circuit modeling. Samples with low NaCl concentration show a doping/impurity effect and a log linear conductivity response as a function of temperature. For higher salt content samples, conductivity increases exponentially with temperature and the log linear relationship no longer holds; instead, we observe phase changes within the samples that follow NaCl–H 2 O–CH 4 phase equilibrium predictions. Final samples were quenched in liquid nitrogen and imaged by cryogenic scanning electron microscopy (cryo‐SEM) to assess grain‐scale characteristics.

Journal of Geophysical Research

Coseismic slip and early afterslip of the M6.0 August 24, 2014 South Napa, California, earthquake

We employ strong motion seismograms and static offsets from the Global Positioning System, Interferometric Synthetic Aperture Radar, and other measurements in order to derive a coseismic slip and afterslip model of the M6.0 24 August 2014 South Napa earthquake. This earthquake ruptured an ∼13‐km‐long portion of the West Napa fault with predominantly right‐lateral strike slip. In the kinematic seismic slip inversions, we couple the coseismic slip and afterslip distributions by requiring both distributions to involve right‐lateral strike slip with positive amplitude, with the net static slip being the sum of the two. We consider several candidate fault geometries: a first involving two steeply east dipping fault planes that reach Earth's surface at the western surface trace (STW), where most surface rupture was observed, a second involving a steeply west dipping plane that also reaches Earth's surface at the STW, and a third involving a combination of two variably west dipping planes constrained to pass through the locus of postseismic seismicity located ∼1 km west of the STW. The data are best fit using the model of two east dipping fault planes, with coseismic slip up to ∼1.2 m on a dominant shallow asperity about 10 km north of the hypocenter and on deeper asperities on the southern part of the rupture. Afterslip up to 1 m is concentrated along the southern part of the rupture at depths 5 km, consistent with surface observations of afterslip. Seismic moments associated with coseismic slip and afterslip are 1.13×10 18 N m (Mw 6.00) and 3.64×10 17 N m, respectively.

California

Semiautomated estimates of directivity and related source properties of small to moderate southern California earthquakes using second seismic moments

We develop a semiautomated method for estimating with second seismic moments the directivity, rupture area, duration, and centroid velocity of earthquakes. The method is applied to 41 southern California earthquakes with magnitude in the range 3.5–5.2 and provides stable results for 28 events. Apparent source time functions (ASTFs) of P and S phases are derived using deconvolution with three stacked empirical Green's functions (seGf). The use of seGf suppresses nongeneric source effects, improves the focal mechanism correspondence to the analyzed earthquakes, and typically allows inclusion of 5 to 15 more ASTFs compared with analysis using a single eGf. Most analyzed earthquakes in the Trifurcation area of the San Jacinto Fault have directivities toward the northwest, while events around Cajon Pass and San Gabriel Mountain tend to propagate toward the southeast. These results are generally consistent with predictions for dynamic rupture on bimaterial interfaces associated with the imaged velocity contrasts in the area. The second moment inversions also provide constraints on the upper and lower bounds of rupture areas in our data set. Stress drops and uncertainties are estimated for elliptical ruptures using the derived characteristic rupture length and width. The semiautomated second moment method with seGfs can be used for routine application to moderate earthquakes in locations with good station coverage.

California

Offset channels may not accurately record strike-slip fault displacement: Evidence from landscape evolution models

Slip distribution, slip rate, and slip per event for strike‐slip faults are commonly determined by correlating offset stream channels—under the assumption that they record seismic slip—but offset channels are formed by the interplay of tectonic and geomorphic processes. To constrain offset channel development under known tectonic and geomorphic conditions, we use numerical landscape evolution simulations along a theoretical strike‐slip fault with uniform and steady uplift, erosion, and diffusion. We investigate the influence of four tectonic parameters (fault zone width, earthquake recurrence interval, variance of the recurrence interval, and total slip relative to channel spacing) on offset channel development through multiple earthquake cycles. Analysis of >3,000 automatically measured offsets from >135 simulations suggests ~30% variability in individual measurements, but modeled displacement is recovered by averaging multiple measurements. However, the average of multiple offset measurements systematically underestimates modeled slip except when the fault zone is less than ~5 m wide, total slip is less than channel spacing, and offsets are measured shortly after an earthquake. In these simulations, postearthquake landscape evolution widens the geomorphic expression of the fault zone and modifies apparent channel offsets. We distinguish this “geomorphic fault zone” from the tectonic fault zone (zone of coseismic distributed deformation). This study highlights the capability of landscape evolution models to explore a range of conditions not easily defined in natural examples and the importance of averaging multiple measurements. Our results verify that paleoseismic studies must consider how geomorphic change has modified offset markers and use caution interpreting slip histories with multiple earthquakes.

Journal of Geophysical Research Solid Earth

Are the stress drops of small earthquakes good predictors of the stress drops of moderate-to-large earthquakes?

The stress drops of small earthquakes often exhibit spatial patterns of variability. If moderate and large earthquakes follow the same spatial patterns, the stress drops of possible future damaging earthquakes could be better predicted by considering the stress drops of nearby small events. Better stress drop predictability could reduce ground-motion uncertainty in Probabilistic Seismic Hazard Assessment (PSHA) and Earthquake Early Warning (EEW). I find that for an internally consistent stress drop catalog of M1.8-3.1 events in southern California, the stress drops of the bigger earthquakes are predictable from the nearby smaller events. However, this catalog only weakly spatially correlates with another catalog of M3.0-5.8 earthquakes, and is spatially uncorrelated with five other stress drop catalogs of M≥3.4 earthquakes. For southern California events M5.5-7.5, stress drops compiled from the literature are weakly spatially correlated with the stress drops of the M1.8-3.1 events, although the correlations are not statistically significant. The lack of strong spatial correlation may be due to actual differences in the controlling factors of stress drop, for example dynamic weakening in moderate-to-large earthquakes. Alternatively, a stronger spatial correlation may exist that is obscured by methodological heterogeneity and large errors in the stress drop estimates. Either way, the stress drops of small earthquakes do not appear to be good predictors of the stress drops of nearby moderate-to-large earthquakes, at least for current techniques of stress drop estimation. If these results are representative, small-earthquake stress drops are not currently useful for substantially reducing uncertainty in PSHA and EEW.

California

Geodetic measurements of slow slip events southeast of Parkfield, CA

Tremor and low-frequency earthquakes are presumed to be indicative of surrounding slow, aseismic slip that is often below geodetic detection thresholds. This study uses data from borehole seismometers and long-baseline laser strainmeters to observe both the seismic and geodetic signatures of episodic tremor and slip on the Parkfield region of the San Andreas Fault near Cholame, CA. The observed occurrence rates of both the tremors and co-located families of low-frequency earthquakes are not steady but instead exhibit quasiperiodic bursts of increased activity. We show that these periods of elevated seismic activity correlate with statistically significant stacked strain signals consisting of 44 slow-slip events. Modeled individual slow-slip events and their total summed moment, which are constrained by seismic signals and stacked strain, respectively, indicate that the individual moment magnitudes of these events range from 4.6–5.2. We find that the measured geodetic signal likely precedes the seismic signal by several hours, consistent with the aseismic slip preceding and driving the observed seismic tremor activity. We confirm that strike-slip faults, in addition to subduction zones, are capable of producing episodic tremor and slip.

California

Methods for rapidly estimating velocity precision from GNSS time series in the presence of temporal correlation: A new method and comparison of existing methods

Time series of position estimates from Global Navigational Satellite System (GNSS) are used to measure the velocities of points on the surface of the Earth. Along with the velocity estimates, a measure of the precision is needed to assess the quality of the velocity measurement. Here, I evaluate rate uncertainties provided by four different methods that have been applied to geodetic time series. The most rigorous approach uses a data covariance that incorporates a variety of noise processes relevant to geodetic time series but is computationally demanding. Two other approaches are efficient algorithms and are used widely, but both can provide less rigorous estimates of the rate uncertainty. I propose and evaluate a fourth method, which provides estimates of rate uncertainty closer to the rigorous approach but is significantly less computationally demanding. I have evaluated all three methods against the more rigorous method using both simulations and time series from 190 GNSS sites. For data best characterized as having a flicker type noise process, one of the widely used methods overestimates the uncertainty by up to a factor of 2, while the other widely used method underestimates the uncertainty by less than a factor of 2. For a random-walk process, both methods underestimate the rate uncertainty by a factor of 3 to 5.

Journal of Geophysical Research

Peak ground velocity spatial variability revealed by dense seismic array in southern California

Understanding and modeling variability of ground motion is essential for building accurate and precise ground motion prediction equations, which can net site‐specific characterization and reduced hazard levels. Here, we explore the spatial variability in peak ground velocity (PGV) at Sage Brush Flats along the San Jacinto Fault in southern California. We use data from a dense array (0.6 x 0.6 km 2 , 1,108 geophones, station spacings 10‐30 m) deployed in 2014 for ~1‐month. These data offer an opportunity to study small scale variability in this region. We examine 38 2≤M L ≤4.2 earthquakes within 200 km of the array. Fault strands and a small basin impact the ground motions, producing PGV variations up to 22% of the mean and a 40% reduction in P and S wave near‐surface velocities. We find along‐fault rupture directivity, source, and path effects can increase PGVs by 167%. Surface PGV measurements exceed the co‐located borehole station (depth 148 m) PGV by factors of 3‐10, confirming the impact on PGV from near surface fault structures, basins, topography, and amplifications from soft sediments. Consistently we find high PGVs within the basin structure. A pair of ~co‐located M L 2.6 events produce repeatable PGV values with similar spatial patterns. The average corner frequencies of these two events are 11‐16 Hz and viable measurements of stress drop can differ by 6.45MPa. Within this small array, the PGV values are variable implying spatial extrapolation of PGV to regions of known faults and basins, even across a small area, should be done with caution.

California

Hydrodynamic and morphologic response of a back-barrier estuary to an extratropical storm

We investigated the hydrodynamic and morphologic response of Barnegat Bay-Little Egg Harbor, New Jersey, USA to Hurricane Sandy. We implemented a three-dimensional, coupled ocean-wave-sediment transport model of the estuary and explored the role of offshore water levels, offshore waves, local winds and waves by systematically removing forcings from a series of simulations. Offshore water levels had the largest impact on water levels in the bay while waves and local wind forcing created substantial spatial variation along the longitudinal axis of the bay. The shape of the bay and its orientation relative to the storm track influenced the response to winds and restricted the maximum water levels in the northern bay and reduced the maximum volume of surge. Basin-average hydrodynamic residence time was reduced by 40%, though its typical spatial distribution remained during the storm. Wave and current-induced bed shear stress resuspended fine sediment resulting in net erosion from the shoals with ensuing net deposition over fringing low-lying land. The net sediment exchange between the bay and the ocean was several times smaller than the exchange at the peak of the storm resulting in negligible net change in the bay volume. Overall, our results suggest that water level responses are highly sensitive to the specific orientation of storm winds relative to the estuary, thereby limiting the utility of simple inundation models. The sediment transport patterns indicate that storms are an important mechanism for redistributing sediment from shoals to fringing wetlands, while net change to sediment budget can be negligible.

New Jersey

Seasonal variation in sediment delivery across the bay-marsh interface of an estuarine salt marsh

Sediment transport across bay–marsh interfaces depends on wave energy, vegetation, and marsh-edge morphology, and varies over a range of timescales. We investigated these dynamics in a tidal salt marsh with a gently-sloped, vegetated edge adjacent to northern San Francisco Bay. Spartina foliosa (cordgrass) inhabits the lower marsh and Salicornia pacifica (pickleweed) predominates on the marsh plain. We measured suspended-sediment concentration (SSC) and hydrodynamics in bay shallows and along a 100-m cross-shore transect in the marsh, during winter and summer. Four-year averaged accretion measured with marker-horizon plots was twice as great along the marsh transect as adjacent to a tidal creek, 50 m from the bay. We estimated deposition and trapping efficiency from the time-series data to assess its variation with season and wave energy. At high tide the transition zone (between cordgrass and pickleweed) was usually erosional, the pickleweed zone was depositional, and both erosion and deposition increased with wave energy, as did the landward position of maximum deposition. Erosion from the transition zone accounted for approximately one-third of the sediment flux into the pickleweed. In the pickleweed zone, SSC, the difference between flood- and ebb-tide SSC and trapping efficiency were greater in summer than winter for comparable wave conditions, which we attribute to increased sediment trapping by dense summer cordgrass. Moderate waves in summer (46%) accounted for more annual accretion in the pickleweed zone than larger waves in winter (28%), although the contribution of winter storms was diminished by the dry winter during the study.

California

Alkalinity in tidal tributaries of the Chesapeake Bay

Despite the important role of alkalinity in estuarine carbon cycling, the seasonal and decadal variability of alkalinity, particularly within multiple tidal tributaries of the same estuary, is poorly understood. Here we analyze more than 26,000 alkalinity measurements, mostly from the 1980s and 1990s, in the major tidal tributaries of the Chesapeake Bay, a large, coastal-plain estuary of eastern North America. The long-term means of alkalinity in tidal-fresh waters vary by a factor of 6 among seven tidal tributaries, reflecting the alkalinity of non-tidal rivers draining to these estuaries. At 25 stations, mostly in the Potomac River Estuary, we find significant long-term increasing trends that exceed the trends in the non-tidal rivers upstream of those stations. Box model calculations in the Potomac River Estuary indicate that the main cause of the estuarine trends is a declining alkalinity sink. The magnitude of this sink is consistent with a simple model of calcification by the invasive bivalve Corbicula fluminea. More generally, in tidal tributaries fed by high-alkalinity non-tidal rivers, alkalinity is consumed, with sinks ranging from 8 to 27% of the upstream input. In contrast, tidal tributaries that are fed by low-alkalinity non-tidal rivers have sources of alkalinity amounting to 34 to 171% of the upstream input. For a single estuarine system, the Chesapeake Bay has diverse alkalinity dynamics and can thus serve as a laboratory for studying the numerous processes influencing alkalinity among the world’s estuaries.

Maryland, Virginia

Surface methane concentrations along the mid-Atlantic bight driven by aerobic subsurface production rather than seafloor gas seeps

Relatively minor amounts of methane, a potent greenhouse gas, are currently emitted from the oceans to the atmosphere, but such methane emissions have been hypothesized to increase as oceans warm. Here, we investigate the source, distribution, and fate of methane released from the upper continental slope of the U.S. Mid-Atlantic Bight, where hundreds of gas seeps have been discovered between the shelf-break and ~1600 m water depth. Using physical, chemical, and isotopic analyses, we identify two main sources of methane in the water column: seafloor gas seeps and in situ aerobic methanogenesis which primarily occurs at 100 – 200 m depth in the water column. Stable isotopic analyses reveal that water samples collected at all depths were significantly impacted by aerobic methane oxidation, the dominant methane sink in this region, with more than 50% of the methane being oxidized, on average. Due to methane oxidation in the deeper water column, below 200 m depth, surface concentrations of methane are influenced more by methane sources found near the surface (0 – 10 m depth) and in the subsurface (10 - 200 m depth), rather than seafloor emissions at greater depths.

Atlantic margin

Assessing the seasonal dynamics of nitrate and sulfate aerosols at the South Pole utilizing stable isotopes

Atmospheric nitrate (NO 3 − = particulate NO 3 − + gas‐phase nitric acid [HNO 3 ]) and sulfate (SO 4 2− ) are key molecules that play important roles in numerous atmospheric processes. Here, the seasonal cycles of NO 3 − and total suspended particulate sulfate (SO 4 2− (TSP) ) were evaluated at the South Pole from aerosol samples collected weekly for approximately 10 months (26 January to 25 October) in 2002 and analyzed for their concentration and isotopic compositions. Aerosol NO 3 − was largely affected by snowpack emissions in which [NO 3 − ] and δ 15 N(NO 3 − ) were highest (49.3 ± 21.4 ng/m 3 , n = 8) and lowest (−47.0 ± 11.7‰, n = 5), respectively, during periods of sunlight in the interior of Antarctica. The seasonal cycle of Δ 17 O(NO 3 − ) reflected tropospheric chemistry year‐round with lower values observed during sunlight periods and higher values observed during dark periods, reflecting shifts from HO x ‐ to O 3 ‐dominated oxidation chemistry. SO 4 2− (TSP) concentrations were highest during austral summer and fall (86.7 ± 73.7 ng/m 3 , n = 18) and are indicated to be derived from dimethyl sulfide (DMS) emissions, as δ 34 S(SO 4 2− ) (TSP) values (18.5 ± 1.0‰, n = 10) were similar to literature δ 34 S(DMS) values. The seasonal cycle of Δ 17 O(SO 4 2− ) (TSP) exhibited minima during austral summer (0.9 ± 0.1‰, n = 5) and maxima during austral fall (1.3 ± 0.3‰, n = 6) and austral spring (1.6 ± 0.1‰, n = 5), indicating a shift from HO x ‐ to O 3 ‐dominated chemistry in the atmospheric derived SO 4 2− component. Overall, the budgets of NO 3 − and SO 4 2− (TSP) at the South Pole were complex functions of transport, localized chemistry, biological activity, and meteorological conditions, and these results will be important for interpretations of oxyanions in ice core records in the interior of Antarctica.

Journal of Geophysical Research D: Atmospheres

Combining numerical and statistical models to predict storm-induced dune erosion

Dune erosion is an important aspect to consider when assessing coastal flood risk, as dune elevation loss makes the protected areas more susceptible to flooding. However, most advanced dune erosion numerical models are computationally expensive, which hinders their application in early-warning systems. Based on a combination of probabilistic and process-based numerical modeling, we develop an efficient statistical tool to predict dune erosion during storms. The analysis focuses on Dauphin Island, AL in the northern Gulf of Mexico, where we combine synthetic sea storms with a calibrated and validated XBeach model to develop and test a range of different surrogate models for their ability to predict barrier-island geometric parameters under storm conditions. Surrogate models are developed by combining the oceanographic forcing from 100 optimally sampled sea-storm events covering the entire multivariate parameter space (used as XBeach input) and associated changes in the dune system (XBeach output). We test four surrogate models using a k-fold approach for validation. All models perform well in predicting changes in dune elevation, barrier-island area, and width, but are less accurate in predicting alterations in the cross-shore locations of dune morphological features. Multivariate adaptive regression splines are identified as the best surrogate model based on their fast development and good performance, attaining a modified Mielke index of 0.81 for dune crest height. As demonstrated at Dauphin Island, our approach shows potential to be used in an operational framework to predict dune response (in particular crest elevation change) when water level and wave forecasts are available.

Alabama

Widespread initiation, reactivation, and acceleration of landslides in the northern California Coast Ranges due to extreme rainfall

Episodically to continuously active slow-moving landslides are driven by precipitation. Climate change, which is altering both the frequency and magnitude of precipitation world21 wide, is therefore predicted to have a major impact on landslides. Here we examine the behavior of hundreds of slow-moving landslides in northern California in response to large changes in annual precipitation that occurred between 2016 and 2018. We quantify the landslide displacement using repeat-pass radar interferometry and pixel offset tracking techniques on a novel dataset from the airborne NASA/JPL Uninhabited Aerial Vehicle Synthetic Aperture Radar. We found that 312 landslides were moving due to extreme rainfall during 2017, compared to 119 during 2016, which was the final year of a historic multi-year drought. However, with a return to below average rainfall in 2018, only 146 landslides remained in motion. The increased number of landslides during 2017 was primarily accommodated by landslides that were smaller than the landslides that remained active between 2016 and 2018. Furthermore, by examining a subset of 51 landslides, we found that 49 had increased velocities during 2017 when compared to 2016. Our results show that slow-moving landslides are sensitive to large changes in annual precipitation, particularly the smaller and thinner landslides that likely experience larger basal pore water pressure changes. Based on climate model predictions for the next century in California, which include increases in annual precipitation and increases in the frequency of dry-to-wet extremes, we hypothesize that there will be an overall increase in landslide activity.

California

Causes of variability in suspended‐sand concentration evaluated using measurements in the Colorado River in Grand Canyon

Rivers commonly exhibit substantial variability in suspended‐sand concentration, even at constant water discharge. Here we derive an approach for evaluating how much of this variability arises from mean bed‐sand grain size. We apply this approach to the Colorado River in Grand Canyon, where discharge‐independent concentration of suspended sand varies by more than a factor of 23 ( N = 1.4 × 10 6 ). Theory predicts that where concentration is controlled by bed‐sand grain size, concentration and grain size in suspension will be inversely correlated (i.e., coarsening of the bed causes suspended sand to become coarser in grain size and lower in concentration). Although the observed correlation is negative, riverbed grain size accounts for only 40% of the variability in concentration. The residuals vary by an order of magnitude; they arise from other processes, such as changes in topography or distribution of sand that cause shear stress to change at constant discharge, changes in the fine tail of bed‐sand grain sizes or changing bedforms. Both bed sand and the other factors influence concentration for durations from less than 1 day to several years. Predictions of concentration based on bed‐sand grain size ( N = 4 × 10 4 ) are less accurate than predictions based on suspended‐sand grain size, probably because suspended sand is a natural integrator of sand‐transporting processes, giving more weight to those areas of the bed that exchange more sand with the flow. Although the causes of variability vary from one river to another, the approach illustrated here is applicable to any river in which concentration varies at constant water discharge.

Arizona

Water tracks enhance water flow above permafrost in upland Arctic Alaska hillslopes

Upland permafrost regions occupy approximately one third of the Arctic landscape. In upland regions, hydrologic fluxes are influenced by water tracks, curvilinear features on hillslopes that preferentially fill with and route water in response to snowmelt and rainfall when the soil above continuous permafrost thaws in the summer. As continued warming of the Arctic may alter hydrologic cycling leading to increased frequency of extreme hydrologic events like drought and flooding as well as modification of biogeochemical cycling, it is imperative to untangle the interplay between precipitation, runoff, and subsurface flow as water is routed from upland Arctic regions to the Arctic Ocean. This study quantifies how ground surface temperatures affect groundwater discharge from hillslopes with water tracks in the upland Arctic by employing a three-dimensional, physically based subsurface flow model with variable saturation and freeze and thaw capabilities that is calibrated to field measurements from the Upper Kuparuk River watershed on the North Slope of Alaska, USA. Model analysis indicates that higher ground surface temperatures along water track hillslopes promote increases in groundwater discharge where water tracks act as conduits for large recharge events and continue to discharge groundwater into the autumn after the adjacent hillslope has frozen. Simulating the conditions that distinguish water tracks from their hillslope watersheds changes subsurface water storage and ground thermal responses but does not alter the total magnitude of groundwater discharge outside of parameter uncertainty. These findings suggest that water tracks play a complex and critical role in hydrologic cycles of the upland Arctic.

Alaska