A novel approach for estimating karst groundwater recharge in mountainous regions and its application in Switzerland

A pragmatic and simple approach for estimating the groundwater recharge of karst aquifers in mountainous regions by extrapolation of the hydrological regimes of gauged and well‐documented systems is presented. Specific discharge rates are derived using annual precipitation and spring measurements by taking into account catchment size and elevation, which are assumed to be the dominant factors. Reference sites with high data reliability are used for calibration and regional extrapolation. This is performed with normalized values employing spatial precipitation deviations and correlation with the elevation of the catchment areas. A tiered step procedure provides minimum and maximum normalized gradients for the relationship between recharge quantity and elevation for karst regions. The normalized recharge can therefore be obtained and extrapolated for any location using the spatial precipitation variability to provide an estimate of annual groundwater recharge. The approach was applied to Switzerland (approximately 7500 km² of karst terrain situated between 200 and over 4000 m a.s.l.) using annual precipitation data from meteorological stations for the years 2000 to 2011. Results show that the average recharge rates of different Swiss karst domains range from 20 to 46 L/km²s, which corresponds to an infiltration ratio between 0.6 and 0.9 of total precipitation. Despite uncertainties inherent in the approach, these results provide a benchmark for renewable karst groundwater resources in Switzerland of about 8.4 km³/year. The approach can be applied to any other mountainous karst region, that is, where a clear relationship between elevation, precipitation and recharge can be assumed. Copyright © 2015 John Wiley & Sons, Ltd.     

Sédimentation récifale rythmique dans le Jurassique superieur du Grand-Saleve (Haute-Savoie,France)

Résumé L'étude statistique des composants minéraux et organiques, selon les méthodes de l'auteur, a été effectuée sur plus de 600 coupes minces régulièrement distribuées dans une série récifale de 225 m d'épaisseur. Elle a révélé des soulèvements et des subsidences encadrant des récifs et se traduisant par une succession de microfaciès clastiques dont les assemblages faunistiques sont très bien délimités au point de vue bathymétrique. Une comparaison avec les formations récifales actuelles est tentée sur la base des faunes périrécifales.     

Possible climatic impacts of tropical deforestation

Large-scale conversion of tropical forests into pastures or annual crops will likely lead to changes in the local microclimate of those regions. Larger diurnal fluctuations of surface temperature and humidity deficit, increased surface runoff during rainy periods and decreased runoff during the dry season, and decreased soil moistrue are to be expected.It is likely that evapotranspiration will be reduced because of less available radiative energy at the canopy level since grass presents a higher albedo than forests, also because of the reduced availability of soil moisture at the rooting zone primarily during the dry season. Recent results from general circulation model (GCM) simulations of Amazonian deforestation seem to suggest that the equilibrium climate for a grassy vegetation in Amazonia would be one in which regional precipitation would be significantly reduced.Global climate changes probably will occur if there is a marked change in rainfall patterns in tropical forest regions as a result of deforestation. Besides that, biomass burning of tropical forests is likely adding CO₂ into the atmosphere, thus contributing to the enhanced greenhouse warming.     

Graphical methods of river profile analysis to unravel drainage area change,uplift and erodibility contrasts in the Central Range of Taiwan

To reveal river channel steepness patterns and variance in settings with significant variation in rock uplift rate, rock erodibility and moving water divides, we present a series of graphical methods to interpret channel profiles. To complement Cartesian χ plots, longitudinal profiles and mapping methods, we introduce a new method based on a radial coordinate system. We map each basin onto polar coordinates in which the radial dimension is χ and the azimuthal coordinate, ?, is calculated with an increment (Δ?) scaled to the distance to neighboring channel heads. The elevation is contoured on this mapping. Average channel steepness is estimated by fitting a conical surface to the elevation. The graph simplifies the comparison of χ and elevation between channels that share a divide, and helps identify spatial changes in drainage area and patterns of erodibility. We apply this approach to derive the uplift pattern in the eastern and southern Central Range of Taiwan, where the high tectonic convergence and uplift rates combined with sub‐tropical climate and frequent typhoons results in high exhumation rate, and well‐developed, detachment‐limited river networks. Additionally, the tectonic activity leads to drainage basin reorganization. We identify examples of divide migration, discrete river capture as well as anomalous steepness that we attribute to local variability in rock erodibility. Estimated basin‐average steepness values show the highest and a near constant value from Hsinwulu basin to Liwu basin in the center of the Island. To the north and south of this region, the values gradually decrease. These estimates show good correlation with the topography of the Central Range and erosion rates derived from in situ ¹⁰Be concentrations in river‐borne quartz. We conclude that the basin steepness reflects systematic differences in rock uplift rate and erosion rate. Copyright © 2016 John Wiley & Sons, Ltd.     

Improved algorithms for generalized finite element simulations of three-dimensional hydraulic fracture propagation

This paper reports improvements to algorithms for the simulation of 3-D hydraulic fracturing with the Generalized Finite Element Method (GFEM). Three optimizations are presented and analyzed. First, an improved initial guess based on solving a 3-D elastic problem with the pressure from the previous step is shown to decrease the number of Newton iterations and increase robustness. Second, an improved methodology to find the time step that leads to fracture propagation is proposed and shown to decrease significantly the number of iterations. Third, reduced computational cost is observed by properly recycling the linear part of the coupled stiffness matrix. Two representative examples are used to analyze these improvements. Additionally, a methodology to include the leak-off term is presented and verified against asymptotic analytical solutions. Conservation of mass is shown to be well satisfied in all examples.     

Improved Priority‐Flood method for depression filling by redundant calculation optimization in local micro‐relief areas

Depression filling is a critical step in distributed hydrological modeling using digital elevation models (DEMs). The traditional Priority‐Flood (PF) approach is widely used due to its relatively high efficiency when dealing with a small‐sized DEM. However, it seems inadequate and inefficient when dealing with large high‐resolution DEMs. In this work, we examined the relationship between the PF algorithm calculation process and the topographical characteristics of depressions, and found significant redundant calculations in the local micro‐relief areas in the conventional PF algorithm. As such calculations require more time when dealing with large DEMs, we thus propose a new variant of the PF algorithm, wherein redundant points and calculations are recognized and eliminated based on the local micro‐relief water‐flow characteristics of the depression‐filling process. In addition, depressions and flatlands were optimally processed by a quick queue to improve the efficiency of the process. The proposed method was applied and validated in eight case areas using the Shuttle Radar Topography Mission digital elevation model (SRTM‐DEM) with 1 arc‐second resolution. These selected areas have different data sizes. A comparative analysis among the proposed method, the Wang and Liu‐based PF, the improved Barnes‐based PF, the improved Zhou‐based PF, and the Planchon and Darboux (P&D) algorithms was conducted to evaluate the accuracy and efficiency of the proposed algorithm. The results showed that the proposed algorithm is 43.2% (maximum) faster than Wang and Liu's variant of the PF method, with an average of 31.8%. In addition, the proposed algorithm achieved similar performance to the improved Zhou‐based PF algorithm, though our algorithm has the advantage of being simpler. The optimal strategies using the proposed algorithm can be employed in various landforms with high efficiency. The proposed method can also achieve good depression filling, even with large amounts of DEM data.     

Stream channel erosion in a rapidly urbanizing region of the US–Mexico border: documenting the importance of channel hardpoints with Structure‐from‐Motion photogrammetry

Urbanization can lead to accelerated stream channel erosion, especially in areas experiencing rapid population growth, unregulated urban development on erodible soils, and variable enforcement of environmental regulations. A combination of field surveys and Structure‐from‐Motion (SfM) photogrammetry techniques was used to document spatial patterns in stream channel geometry in a rapidly urbanizing watershed, Los Laureles Canyon (LLCW), in Tijuana, Mexico. Ground‐based SfM photogrammetry was used to map channel dimensions with 1 to 2 cm vertical mean error for four stream reaches (100–300 m long) that were highly variable and difficult to survey with a differential GPS. Regional channel geometry curves for LLCW had statistically larger slopes and intercepts compared with regional curves developed for comparable, undisturbed reference channels. Cross‐sectional areas of channels downstream of hardpoints, such as concrete reaches or culverts, were up to 64 times greater than reference channels, with enlargement persisting, in some cases, up to 230 m downstream. Percentage impervious cover was not a good predictor of channel enlargement. Proximity to upstream hardpoint, and lack of riparian and bank vegetation paired with highly erodible bed and bank materials may account for the instability of the highly enlarged and unstable cross‐sections. Channel erosion due to urbanization accounts for approximately 25–40% of the total sediment budget for the watershed, and channel erosion downstream of hardpoints accounts for one‐third of all channel erosion. Channels downstream of hardpoints should be stabilized to prevent increased inputs of sediment to the Tijuana Estuary and local hazards near the structures, especially in areas with urban settlements near the stream channel. Copyright © 2017 John Wiley & Sons, Ltd.     

Accuracy assessment of the Mexican National Forest Inventory map: A study in four ecogeographical areas

The accuracy of the Mexican National Forest Inventory (NFI) map is derived in four distinct ecogeographical areas, using an assessment design tailored for the project. A main achievement of the design was to integrate the high diversity of classes encompassed at the most detailed subcommunity level of the classification scheme within a cost‐controlled statistically sound assessment. A hybrid double sampling strategy was applied to the 2.5 million‐ha study area. A total of 5955 reference sites were verified against their NFI map label. The availability of detailed quasi‐synchronous reference data for the 2000 Landsat‐derived NFI and the high diversity of mapped classes allowed a careful thematic analysis on the selected regions, relevant for national extrapolation. Global accuracy estimates of 64–78 per cent were registered among the four ecogeographical areas (two with mainly temperate climate and the other two with mainly tropical climate), with the lower accuracy levels found in areas more densely covered with forests. According to the estimates, the NFI map tends to underestimate the presence of temperate forest (especially oak) and overestimate the presence of tropical forest in the areas investigated. The analysis of confusions reveals difficulties in unambiguously interpreting or labelling forests with secondary vegetation, herbaceous and/or shrub‐like vegetation as well as distinguishing between aquatic vegetation types. The design proved useful from the perspective of accuracy assessments of regional maps in biodiverse regions.     

Analysis of the Devonian (Frasnian) platform from Belgium: a multi‐faceted approach for basin evolution reconstruction

This study undertakes a multi‐disciplinary approach (sedimentology, carbon isotopes, magnetic susceptibility and thickness distribution) to improve understanding of a major Palaeozoic carbonate platform, the Frasnian platform of Belgium. These combined techniques are used to reconstruct the platform history, which evolved in two main steps. During the first phase, the basin was strongly influenced by faulting, producing notable thickness and facies variations, with open ocean conditions, with good water circulation and no/or limited barrier reef. The second phase of platform development was less influenced by differential subsidence, as indicated by homogeneous facies distribution. However, this platform developed under restricted waters, with low circulation which is likely related to the development of a barrier reef.     

ISA-Hypoplasticity accounting for cyclic mobility effects for liquefaction analysis

Acta Geotechnica - The hypoplastic model for sands proposed by Wolffersdorff (Mech Cohes Frict Mater 1: 251–271, 1996) combined with the intergranular strain anisotropy by Fuentes and...