Models of crustal flow in the India–Asia collision zone

Surface velocities in parts of the India–Asia collision zone are compared to velocities calculated from equations describing fluid flow driven by topographically produced pressure gradients. A good agreement is found if the viscosity of the crust is ∼10²⁰ Pa s in southern Tibet and ∼10²² Pa s in the area between the Eastern Syntaxis and the Szechwan Basin. The lower boundary condition of the flow changes between these two areas, with a stress-free lower boundary in the area between the Szechwan basin and the Eastern Syntaxis, and a horizontally rigid but vertically deformable boundary where strong Indian lithospheric material underlies southern Tibet. Deformation maps for olivine, diopside and anorthite show our findings to be consistent with laboratory measurements of the rheology of minerals. Gravitationally driven flow is also suggested to be taking place in the Indo–Burman Ranges, with a viscosity of ∼10¹⁹–10²⁰ Pa s. Flow in both southern Tibet and the Indo–Burman Ranges provides an explanation for the formation of the geometry of the Eastern Himalayan Syntaxis. The majority of the normal faulting earthquakes in the Tibetan Plateau occur in the area of southern Tibet which we model as gravitationally spreading over the Indian shield.     

SUMER spectral observations of post-flare supra-arcade inflows

On 21 April 2002 a large eruptive flare on the west limb of the Sun developed a bright, very dynamic, post-flare arcade. In TRACE 195 Å images, a series of dark, sunward moving flows were seen against the bright extreme ultraviolet (EUV) arcade. SUMER obtained a series of spectra of the dark EUV flows in the lines Cii, Fexii, and Fexxi at a fixed position above the limb. These spectra give spatially resolved line-of-sight velocities and emission measures for the arcade plasma over a temperature range 2×10⁴ to 10⁷ K. The flows are dark in all SUMER lines. The UV continuum longward (1350 Å) and shortward (675 Å) of the hydrogen Lyman limit is used to determine whether the dark 195 Å inflows are due to regions of low plasma density (plasma voids) or cold absorbing material. There is some evidence of absorption near the front of one of the inflows; however, along most of the dark channels there is no change in continuum ratio and we therefore conclude, as originally suggested by McKenzie and Hudson (1999), that they are plasma voids.     

Mapping Paleo‐Fire Boundaries from Binary Point Data: Comparing Interpolation Methods

Fire history studies have traditionally emphasized temporal rather than spatial properties of paleo‐fire regimes. In this study we compare four methods of mapping paleo‐fires in central Washington from binary point data: indicator kriging, inverse distance weighting, Thiessen polygons, and an expert approach. We evaluate the results of each mapping method using a test (validation) dataset and receiver operating characteristic plots. Interpolation methods perform well, but results vary with fire size and spatial pattern of points. Though all methods involve some subjectivity, automated interpolation methods perform well, are replicable, and can be applied across varying landscapes.     

The East Anatolian Fault: A major structure in Eastern Turkey

Satellite photographs of Eastern Turkey show a large strike slip fault about 550 km long extending from the Gulf of Iskenderun to the North Anatolian Fault with a strike of 60°E. This fault takes up most of the motion between the Turkish and Arabian plates.     

Does hillslope trenching enhance groundwater recharge and baseflow in the Peruvian Andes?

As Andean glaciers rapidly retreat due to climate change, the balance of groundwater and glacial meltwater contributions to stream discharge in tropical, proglacial watersheds will change, potentially increasing vulnerability of water resources. The Shullcas River Watershed, near Huancayo, Peru, is fed only partly by the rapidly receding Huaytapallana glaciers (<20% of dry season flow). To potentially increase recharge and therefore increase groundwater derived baseflow, the government and not‐for‐profit organizations have installed trenches along large swaths of hillslope in the Shullcas Watershed. Our study focuses on a nonglacierized subcatchment of the Shullcas River Watershed and has 2 objectives: (a) create a model of the Shullcas groundwater system and assess the controls on stream discharge and (b) investigate the impact of the infiltration trenches on recharge and baseflow. We first collected hydrologic data from the field including a year‐long hydrograph (2015–2016), meteorological data (2011–2016), and infiltration measurements. We use a recharge model to evaluate the impact of trenched hillslopes on infiltration and runoff processes. Finally, we use a 3‐dimensional groundwater model, calibrated to the measured dry season baseflow, to determine the impact of trenching on the catchment. Simulations show that trenched hillslopes receive approximately 3.5% more recharge, relative to precipitation, compared with unaltered hillslopes. The groundwater model indicates that because the groundwater flow system is fast and shallow, incorporating trenched hillslopes (~2% of study subcatchment area) only slightly increases baseflow in the dry season. Furthermore, the location of trenching is an important consideration: Trenching higher in the catchment (further from the river) and in flatter terrain provides more baseflow during the dry season. The results of this study may have important implications for Andean landscape management and water resources.     

Thermal Imagery of Groundwater Seeps: Possibilities and Limitations

Quantifying groundwater flow at seepage faces is crucial because seepage faces influence the hydroecology and water budgets of watersheds, lakes, rivers and oceans, and because measuring groundwater fluxes directly in aquifers is extremely difficult. Seepage faces provide a direct and measurable groundwater flux but there is no existing method to quantitatively image groundwater processes at this boundary. Our objective is to determine the possibilities and limitations of thermal imagery in quantifying groundwater discharge from discrete seeps. We developed a conceptual model of temperature below discrete seeps, observed 20 seeps spectacularly exposed in three dimensions at an unused limestone quarry and conducted field experiments to examine the role of diurnal changes and rock face heterogeneity on thermal imagery. The conceptual model suggests that convective air‐water heat exchange driven by temperature differences is the dominant heat transfer mechanism. Thermal imagery is effective at locating and characterizing the flux of groundwater seeps. Areas of active groundwater flow and ice growth can be identified from thermal images in the winter, and seepage rates can be differentiated in the summer. However, the application of thermal imagery is limited by diverse factors including technical issues of image acquisition, diurnal changes in radiation and temperature, and rock face heterogeneity. Groundwater discharge rates could not be directly quantified from thermal imagery using our observations but our conceptual model and experiments suggest that thermal imagery could quantify groundwater discharge when there are large temperature differences, simple cliff faces, non‐freezing conditions, and no solar radiation.     

Regional variability in dust‐on‐snow processes and impacts in the Upper Colorado River Basin

Dust deposition onto mountain snow cover in the Upper Colorado River Basin frequently occurs in the spring when wind speeds and dust emission peaks on the nearby Colorado Plateau. Dust loading has increased since the intensive settlement in the western USA in the mid 1880s. The effects of dust‐on‐snow have been well studied at Senator Beck Basin Study Area (SBBSA) in the San Juan Mountains, CO, the first high‐altitude area of contact for predominantly southwesterly winds transporting dust from the southern Colorado Plateau. To capture variability in dust transport from the broader Colorado Plateau and dust deposition across a larger area of the Colorado River water sources, an additional study plot was established in 2009 on Grand Mesa, 150 km to the north of SBBSA in west central, CO. Here, we compare the 4‐year (2010–2013) dust source, deposition, and radiative forcing records at Grand Mesa Study Plot (GMSP) and Swamp Angel Study Plot (SASP), SBBSA's subalpine study plot. The study plots have similar site elevations/environments and differ mainly in the amount of dust deposited and ensuing impacts. At SASP, end of year dust concentrations ranged from 0.83 mg g^?1 to 4.80 mg g^?1, and daily mean spring dust radiative forcing ranged from 50–65 W m^?2, advancing melt by 24–49 days. At GMSP, which received 1.0 mg g^?1 less dust per season on average, spring radiative forcings of 32–50 W m^?2 advanced melt by 15–30 days. Remote sensing imagery showed that observed dust events were frequently associated with dust emission from the southern Colorado Plateau. Dust from these sources generally passed south of GMSP, and back trajectory footprints modelled for observed dust events were commonly more westerly and northerly for GMSP relative to SASP. These factors suggest that although the southern Colorado Plateau contains important dust sources, dust contributions from other dust sources contribute to dust loading in this region, and likely account for the majority of dust loading at GMSP. Copyright © 2015 John Wiley & Sons, Ltd.     

The importance of incorporating diurnally fluctuating stream discharge in stream temperature energy balance models

Although stream temperature energy balance models are useful to predict temperature through time and space, a major unresolved question is whether fluctuations in stream discharge reduce model accuracy when not exactly represented. However, high‐frequency (e.g., subdaily) discharge observations are often unavailable for such simulations, and therefore, diurnal streamflow fluctuations are not typically represented in energy balance models. These fluctuations are common due to evapotranspiration, snow pack or glacial melt, tidal influences within estuaries, and regulated river flows. In this work, we show when to account for diurnally fluctuating streamflow. To investigate how diurnal streamflow fluctuations affect predicted stream temperatures, we used a deterministic stream temperature model to simulate stream temperature along a reach in the Quilcayhuanca Valley, Peru, where discharge varies diurnally due to glacial melt. Diurnally fluctuating streamflow was varied alongside groundwater contributions via a series of computational experiments to assess how uncertainty in reach hydrology may impact simulated stream temperature. Results indicated that stream temperatures were more sensitive to the rate of groundwater inflow to the reach compared with the timing and amplitude of diurnal fluctuations in streamflow. Although incorporating observed diurnal fluctuations in discharge resulted in a small improvement in model RMSE, we also assessed other diurnal discharge signals and found that high amplitude signals were more influential on modelled stream temperatures when the discharge peaked at specific times. Results also showed that regardless of the diurnal discharge signal, the estimated groundwater flux to the reach only varied from 1.7% to 11.7% of the upstream discharge. However, diurnal discharge fluctuations likely have a stronger influence over longer reaches and in streams where the daily range in discharge is larger, indicating that diurnal fluctuations in stream discharge should be considered in certain settings.     

Lacustrine oxygen isotope record of 20th-century climate change in central EuropeL evaluation of climatic controls on oxygen isotopes in precipitation

We report oxygen isotope data from a 108-yr (1885–1993) sequence with annual laminae of bio-induced authigenic calcite in a frozen core from Baldeggersee, a small lake in Central Switzerland. These isotope results provide proxy data on the isotopic composition of past precipitation in the Baldeggersee catchment region and are quantitatively compared with instrumental climate data (i.e. mean annual air temperature and atmospheric circulation pattern indices) to evaluate climatic controls on oxygen isotopes in precipitation.Monitoring the isotope hydrology of Baldeggersee demonstrates that the oxygen isotopic composition of the lake water is controlled by the isotopic composition of local atmospheric precipitation (¹⁸O_p) and that the isotopic signal of precipitation is preserved, albeit damped, in the lake calcite oxygen isotope record (¹⁸O_c). Comparison of the calcite oxygen isotope proxy for ¹⁸O_p in the catchment with historical mean annual air temperature measurements from Bern, Switzerland confirms that authigenic calcite reliably records past annual air temperature in the region. This ¹⁸O_c/temperature relationship is calculated to be 0.39/°C for the period 1900–1960, based on an isotope mass-balance model for Baldeggersee. An exception is a 0.8 anomalous negative shift in calcite ¹⁸O values since the 1960s. Possible explanations for this recent ¹⁸O_c shift, as it is not related to mean annual air temperature, include changes in ¹⁸O_p due to synoptic circulation patterns. In particular, the 0.8 negative shift coincides with a trend towards a more dominant North Atlantic Oscillation (NAO) index. This circulation pattern would tend to bring more isotopically more negative winter precipitation to the region and could account for the 0.8 offset in ¹⁸O_c data.