Characterization of complex groundwater flows in the environment of singular buildings by combining hydrogeological and non-destructive geophysical (ground-penetrating radar) techniques: Punta Begoña Galleries (Getxo,Spain)

Locating and quantifying groundwater flow in many built-up areas are a priority with regard to its complete restoration. In this work, a hydrogeological survey of the surroundings of the Punta Begoña Galleries (Getxo, Bizkaia), built on a coastal cliff, was completed by using ground penetrating radar (GPR) testing. Thus, the preliminary characterization of soils and rocks in accessible areas of the cliff was first improved by hydrogeological information gathered from a single survey borehole, including permeability measurements by low pressure injection tests (LPTs) and continuous water level monitoring. As a complementary method, the non-destructive GPR technique was performed during both dry and wet hydrological periods and in tandem with the injection tests, providing more complete spatial and temporal images of water flows. Specifically, GPR allows mapping of flow paths in soils and assessing the continuity of fractures in rock masses. Altogether, this complementary approach provides greater knowledge of complex underground flow dynamics in built environments, thus making it easier to make decisions for their management.     

A stress and displacement discontinuity element method for elastic transversely anisotropic rock

The paper presents closed‐form solutions for stress and displacement influence functions for stress discontinuity (SD) and displacement discontinuity (DD) elements, for a two‐dimensional plane‐strain elastic, transversely anisotropic medium. The solutions for SD elements are based on Kelvin's problem and for DD elements on the concept of dipoles. Stress and displacement influence functions are derived for the following elements: constant SD, linear SD, constant DD, linear DD, square root DD, parabolic DD, constant DD surface, and linear DD surface elements. The formulations are incorporated into FROCK, a hybridized boundary element method code, and are validated by providing comparisons between the results from FROCK and the finite element code ABAQUS. A limited parametric analysis shows the effects of slight anisotropy on the stress field around the tip of a crack and of the orientation of the crack with respect to the axes of elastic symmetry. Copyright © 2014 John Wiley & Sons, Ltd.     

Extracting drainage networks and their connectivity using LiDAR data

Policies, measures, and models geared towards flood prevention and managing surface waters benefit from high quality data on the presence and characteristics of drainage ditches. As a cost and labour effective alternative for acquiring such data through field surveys, we propose a method (a) to extract vector data representing ditch drainage networks based on local morphologic features derived from high resolution digital elevation models (DEM) and (b) to identify possible connections in the ditch network by calculating a probability of the connectivity using a logistic regression where the predictor variables are characteristics of the ditch centre lines or derived from the DEM. Using Light Detection and Ranging (LiDAR) derived DEMs with a 1 m resolution, the method was developed and tested for a mixed agricultural residential area in north‐eastern Belgium. The derived ditch segments had an error of omission of 8% and an error of commission of 5%. The original positional accuracy of the centre lines of the extracted ditches was 0.6 m and could be improved to 0.4 m by shifting each vertex to the position of the lowest LiDAR point located within a radius equal to the spatial resolution of the used DEM. About 69% of the false disconnections in the network were identified and corrected leading to a reduction of the unconnected parts of the ditch network by 71%. The extracted and connected network approximated the reference ditch network fairly well.     

Spatially distributed tracer‐aided modelling to explore water and isotope transport,storage and mixing in a pristine,humid tropical catchment

Rapidly transforming headwater catchments in the humid tropics provide important resources for drinking water, irrigation, hydropower, and ecosystem connectivity. However, such resources for downstream use remain unstudied. To improve understanding of the behaviour and influence of pristine rainforests on water and tracer fluxes, we adapted the relatively parsimonious, spatially distributed tracer‐aided rainfall–runoff (STARR) model using event‐based stable isotope data for the 3.2‐km² San Lorencito catchment in Costa Rica. STARR was used to simulate rainforest interception of water and stable isotopes, which showed a significant isotopic enrichment in throughfall compared with gross rainfall. Acceptable concurrent simulations of discharge (Kling–Gupta efficiency [KGE] ~0.8) and stable isotopes in stream water (KGE ~0.6) at high spatial (10 m) and temporal (hourly) resolution indicated a rapidly responding system. Around 90% of average annual streamflow (2,099 mm) was composed of quick, near‐surface runoff components, whereas only ~10% originated from groundwater in deeper layers. Simulated actual evapotranspiration (ET) from interception and soil storage were low (~420 mm/year) due to high relative humidity (average 96%) and cloud cover limiting radiation inputs. Modelling suggested a highly variable groundwater storage (~10 to 500 mm) in this steep, fractured volcanic catchment that sustains dry season baseflows. This groundwater is concentrated in riparian areas as an alluvial–colluvial aquifer connected to the stream. This was supported by rainfall–runoff isotope simulations, showing a “flashy” stream response to rainfall with only a moderate damping effect and a constant isotope signature from deeper groundwater (~400‐mm additional mixing volume) during baseflow. The work serves as a first attempt to apply a spatially distributed tracer‐aided model to a tropical rainforest environment exploring the hydrological functioning of a steep, fractured‐volcanic catchment. We also highlight limitations and propose a roadmap for future data collection and spatially distributed tracer‐aided model development in tropical headwater catchments.     

Late Quaternary (Marine Isotope Stage 3 to Recent) sedimentary evolution of the Basque shelf (southern Bay of Biscay)

Late Quaternary (MIS 3 to Recent) oceanographic evolution of the Basque shelf has been analysed for the first time based on the sedimentological analysis of three cores obtained from the middle and outer shelves. The cores are located in two interfluves separated by the San Sebastian canyon. The variability of the percentage of the planktonic foraminifera species Neogloboquadrina pachyderma sin. and of δ¹⁸O_bull allowed us to identify the influence of colder and warmer waters in the Basque shelf during the late Quaternary. From ～56 cal. ka BP to the end of the Last Glacial Maximum (19 cal. ka BP) the sedimentary record shows a decreasing trend in the mean grain size that correlates with the eustatic sea‐level fall. The last Deglaciation (19–11.5 cal. ka BP) is characterized by a sea‐level rise that produced an important hiatus in the western outer shelf. During the Holocene, the middle and outer shelves present different behaviours. From 11.5 to 6.7 cal. ka BP, in the outer shelf the sea‐level rise that started during the Deglaciation produced a hiatus, whereas in the middle shelf the sedimentary succession records the presence of warm to temperate waters. Between 6.7–4.9 cal. ka BP, the entrance of cold surface water‐masses that only affected the middle shelf has been identified, and temperate to warm waters occurred in the outer shelf. The cold surface water‐masses retreated during 4.9–4.3 cal. ka BP in the middle shelf. Finally, from 4.3 cal. ka BP to Recent, the middle shelf registers a hiatus due to sea‐level stabilization after a generalized transgression, synchronous to a decrease in the energy of the water‐masses in the outer shelf. In conclusion, the environmental changes detected in the Basque shelf are attributed to both regional and eustatic sea‐level changes.     

Stratigraphic architecture of Solander Basin records Southern Ocean currents and subduction initiation beneath southwest New Zealand

Solander Basin is characterized by subduction initiation at the Pacific‐Australia plate boundary, where high biological productivity is found at the northern edge of the Antarctic Circumpolar Current. Sedimentary architecture results from tectonic influences on accommodation space, sediment supply and ocean currents (via physiography); and climate influence on ocean currents and biological productivity. We present the first seismic‐stratigraphic analysis of Solander Basin based on high‐fold seismic‐reflection data (voyage MGL1803, SISIE). Solander Trough physiography formed by Eocene rifting, but basinal strata are mostly younger than ca. 17 Ma, when we infer Puysegur Ridge formed and sheltered Solander Basin from bottom currents, and mountain growth onshore increased sediment supply. Initial inversion on the Tauru Fault started at ca. 15 Ma, but reverse faulting from 12 to ca. 8 Ma on both the Tauru and Parara Faults was likely associated with reorganization and formation of the subduction thrust. The new seabed topography forced sediment pathways to become channelized at low points or antecedent gorges. Since 5 Ma, southern Puysegur Ridge and Fiordland mountains spread out towards the east and Solander Anticline grew in response to ongoing subduction and growth of a slab. Solander Basin had high sedimentation rates because (1) it is sheltered from bottom currents by Puysegur Ridge; and (2) it has a mountainous land area that supplies sediment to its northern end. Sedimentary architecture is asymmetric due to the Subtropical Front, which moves pelagic and hemi‐pelagic sediment, including dilute parts of gravity flows, eastward and accretes contourites to the shelf south of Stewart Island. Levees, scours, drifts and ridges of folded sediment characterize western Solander Basin, whereas hemi‐pelagic drape and secondary gravity flows are found east of the meandering axial Solander Channel. The high‐resolution record of climate and tectonics that Solander Basin contains may yield excellent sites for future scientific ocean drilling.     

Role of Fracturing in the Organization of the Karst Features of Azrou Plateau (Middle Atlas,Morroco) Studied by Remote Sensing Imagery

Karst formation geometry can be controlled by fractures and faults, and by other lithologies. Here we show the organization of kastic collapse features related to structures and to extensive basaltic lava flows in the Middle Atlas of Morocco. A lineament map of major faults and fractures has been created for the Middle Atlas using Landsat 7 ETM⁺ satellite images. This shows a dominant NE–SW regional direction and less prominent NNW–SSE and ENE–WSE directions. All these directions coincide with the alignments of karstic depressions that have formed in the Liassic limestones. The basaltic flows covering these formations on the Middle Atlas limestone plateau, have allowed the generation of cryptokarst, geometrically organized a long these major lineament directions. Karst landforms probabaly existed before the eruption of the lavas, but there were partly invaded by intrusions and volcanism. The extensive basaltic flows allowed for increased infiltration, and subsurface water flow, increasing the rate of kast formation after eruptions. Some basins show evidence of increased subsidence after lava emplacement (Aguelmam Sidi Ali Lake) and some maar-like craters also have subsided after eruption, by karts formation. We lay out the structural and lithological controls on Karstic formation in an intraplate volcanic field based on limestones and evaporites.     

Evaluating the Ability of CRCM5 to Simulate Mixed Precipitation

Precipitation episodes in the form of freezing rain and ice pellets represent natural hazards affecting eastern Canada during the cold season. These types of precipitation mainly occur in the St. Lawrence River valley and the Atlantic provinces of Canada. This study aims to evaluate the ability of the fifth-generation Canadian Regional Climate Model (CRCM5), using a 0.11° horizontal grid mesh, to hindcast mixed precipitation when driven by reanalyses produced by the European Centre for Medium-range Weather Forecasts (ERA-Interim) for a 35-year period. In general, the CRCM5 simulation slightly overestimates the occurrence of freezing rain, but the geographical distribution is well reproduced. The duration of freezing rain events and accompanying surface winds in the Montréal region are reproduced by CRCM5. A case study is performed for an especially catastrophic freezing-rain event in January 1998; the model succeeds in simulating the intensity and duration of the episode, as well as the propitious meteorological environment. Overall, the model is also able to reproduce the climatology and a specific event of freezing rain and ice pellets.     

Using low-density discrete Airborne Laser Scanning data to assess the potential carbon dioxide emission in case of a fire event in a Mediterranean pine forest

The aim of study is to map the carbon dioxide (CO₂) emission of the aboveground tree biomass (AGB) in case of a fire event. The suitability of low point density, discrete, multiple-return, Airborne Laser Scanning (ALS) data and the influence of several characteristics of these data and the study area on the results obtained have been evaluated. A sample of 45 circular plots representative of Pinus halepensis Miller stands were used to fit and validate the model of AGB. The ALS point clouds were processed to obtain the independent variables and a multivariate linear regression analysis between field data and ALS-derived variables allowed estimation of AGB. Then, the influence of several characteristics on the residuals of the model was analyzed. Finally, conversion factors were applied to obtain the CO₂ values. The AGB model presented a R² value of 0.84 with a relative root-mean-square error of 27.35%. This model included ALS variables related to vegetation height variability and to canopy density. Terrain slope, aspect, canopy cover, scan angle and the number of laser returns did not influence AGB estimations at plot level.