A 16,000 C yr B.P. packrat midden series from the USA–Mexico Borderlands

A new packrat midden chronology from Playas Valley, southwestern New Mexico, is the first installment of an ongoing effort to reconstruct paleovegetation and paleoclimate in the U.S.A.–Mexico Borderlands. Playas Valley and neighboring basins supported pluvial lakes during full and/or late glacial times. Plant macrofossil and pollen assemblages from nine middens in the Playas Valley allow comparisons of two time intervals: 16,000–10,000 and 4000–0 ¹⁴C yr B.P. Vegetation along pluvial lake margins consisted of open pinyon–juniper communities dominated by Pinus edulis, Juniperus scopulorum, Juniperus cf. coahuilensis, and a rich understory of C₄ annuals and grasses. This summer-flowering understory is also characteristic of modern desert grassland in the Borderlands and indicates at least moderate summer precipitation. P. edulis and J. scopulorum disappeared or were rare in the midden record by 10,670 ¹⁴C yr B.P. The late Holocene is marked by the arrival of Chihuahuan desert scrub elements and few departures as the vegetation gradually became modern in character. Larrea tridentata appears as late as 2190 ¹⁴C yr B.P. based on macrofossils, but may have been present as early as 4095 ¹⁴C yr B.P. based on pollen. Fouquieria splendens, one of the dominant desert species present at the site today, makes its first appearance only in the last millennium. The midden pollen assemblages are difficult to interpret; they lack modern analogs in surface pollen assemblages from stock tanks at different elevations in the Borderlands.     

Biomorphodynamics of alternate bars in a channelized,regulated river: An integrated historical and modelling analysis

The development of alternate bars in channelized rivers can be explained theoretically as an instability of the riverbed when the active channel width to depth ratio exceeds a threshold. However, the development of a vegetation cover on the alternate bars of some channelized rivers and its interactions with bar morphology have not been investigated in detail. Our study focused on the co‐evolution of alternate bars and vegetation along a 33 km reach of the Isère River, France. We analysed historical information to investigate the development of alternate bars and their colonization by vegetation within a straightened, embanked river subject to flow regulation, sediment mining, and vegetation management. Over an 80 year period, bar density decreased, bar length increased, and bar mobility slowed. Vegetation encroachment across bar surfaces accompanied these temporal changes and, once established, vegetation cover persisted, shifting the overall system from an unvegetated to a vegetated dynamic equilibrium state. The unvegetated morphodynamics of the impressively regular sequence of alternate bars that developed in the Isère following channelization is consistent with previous theoretical morphodynamic work. However, the apparent triggering dynamics of vegetation colonization needs to be investigated, based on complex biophysical instability processes. If instability related to vegetation colonization is confirmed, further work needs to focus on the relevance of initial conditions for this instability, and on related feedback effects such as how the morphodynamics of bare‐sediment alternate bars may have affected vegetation development and, in turn, how vegetation has created a new dynamic equilibrium state. Copyright © 2018 John Wiley & Sons, Ltd.     

Trophic importance of diatoms in an intertidal Zostera noltii seagrass bed: Evidence from stable isotope and fatty acid analyses

A current predominant paradigm emphasizes the role of epiphytic algae for invertebrates in most seagrass food webs. However, in some intertidal Zostera noltii beds, epiphyte biomass is very low compared to microphytobenthos and seagrass biomasses. We assessed the role of microphytobenthos in a temperate intertidal Z. noltii bed by combining stable isotope and fatty acid (FA) analyses on primary producers, composite sources — suspended particulate organic matter (SPOM) and sediment surface organic matter (SSOM) — and the main macrofaunal consumers. Z. noltii showed high δ¹³C (−9.9‰) and high 18:2(n-6) and 18:3(n-3) contents. Microphytobenthos was slightly more ¹³C-depleted (−15.4‰) and had high levels of diatom markers: 14:0, 16:1(n-7)c, 20:5(n-3). Low mean δ¹³C (−22.0‰) and large amounts of diatom and bacteria (18:1(n-7)c) markers indicated that SPOM was mainly composed of a mixture of fresh and decayed pelagic diatoms. Higher mean δ¹³C (−17.9‰) and high amounts of diatom FAs were found in SSOM, showing that microphytobenthic diatoms dominate. Very low percentages of 18:2(n-6) and 18:3(n-3) in consumers indicated a low contribution of Z. noltii material to their diets. Grazers, deposit and suspension-deposit feeders had δ¹³C close to microphytobenthos and high levels of diatom FAs, confirming that microphytobenthos represented the main part of their diet. Lower δ¹³C and higher amounts of flagellate FAs – 22:6(n-3) and 16:4(n-3) – in suspension feeders indicated that their diet resulted from a mixture of SPOM and microphytobenthos. These results demonstrate that invertebrates do not consume high amounts of seagrass and highlight the main role of benthic diatoms in this intertidal seagrass bed.     

How basin model results enable the study of multi-layer aquifer response to pumping: the Paris Basin, France

The Albian aquifer of the Paris Basin (France) has been exploited since 1841 and shows drastic drawdown. A three-dimensional (3D) groundwater flow model is used to study the hydrodynamic response of the multi-layered aquifers to pumping activity in the Albian, at basin scale over 167 years. This 3D flow model uses geometry and hydrodynamic parameter distributions that are inherited from a genetic approach through basin modelling, the basin model creating a geometric pattern of hydrodynamic properties constrained by geological history. The paper aims to promote the use of the basin model approach (long time scale, 248 Ma) for the study of deep-aquifer response to anthropogenic perturbation (short time scale, 167 years) in situations for which hydrodynamic data are scarce but geological data are numerous. The results show that parameter distribution is insufficient to reproduce the Albian aquifer behaviour, notably highlighting a different meaning of the specific storage coefficient between basin modelling and groundwater-flow modelling. Dividing the storage coefficient by 100 and including available transmissivity data significantly improved the model/data comparison. The potential impact on a deep aquitard is then discussed. This study sheds light on the advantages and limitations of the basin model approach for groundwater-flow modelling in 3D.     

A watershed-scale study of climate change impacts on groundwater recharge (Annapolis Valley,Nova Scotia,Canada)

Abstract

The potential impacts of future climate change on the evolution of groundwater recharge are examined at a local scale for a 546-km² watershed in eastern Canada. Recharge is estimated using the infiltration model Hydrologic Evaluation of Landfill Performance (HELP), with inputs derived from five climate runs generated by a regional climate model in combination with the A2 greenhouse gas emissions scenario. The model runs project an increase in annual recharge over the 2041–2070 period. On a seasonal basis, however, a marked decrease in recharge during the summer and a marked increase during the winter are observed. The results suggest that increased evapotranspiration resulting from higher temperatures does not offset the large increase in winter infiltration. In terms of individual water budget components, clear differences are obtained for the different climate change scenarios. Monthly recharge values are also found to be quite variable, even for a given climate scenario. These findings are compared with results from two regional-scale studies.

Editor D. Koutsoyiannis; Associate editor M. Besbes     

Tsunami Mapping Related to Local Earthquakes on the French–Italian Riviera (Western Mediterranean)

The Ligurian coast, located at the French–Italian border, is densely populated as well as a touristic area. It is also a location where earthquakes and underwater landslides are recurrent. The nature of the local tsunamigenesis is therefore a legitimate question, because no tsunami warning system can resolve tsunami arrival times of a few minutes, which is the case for the area. As far as the seismicity of the area is concerned, the frequent recurrent earthquakes are generally of moderate magnitude: most of them are lower than M _w 5. However, the relatively large M _w 6.9 earthquake (Larroque et al., in Geophys J Int, 2012. doi:10.1111/j.1365-246X.2012.05498.x) that occurred on the February 23, 1887, offshore of Imperia (Italian Riviera) is quite emblematic. This unusual event for the region merits a complete study: the quantification of its rupture mechanism is essential (1) to understand the regional active deformation, but also (2) to evaluate its tsunamigenesis potential by deriving relevant rupture scenarios obtained from our knowledge of the event; for that purpose the event is extensively described here. The first point has been the subject of quite a few studies based on the seismotectonics of the area. The last documented approach has been completed by Larroque et al. (Geophys J Int, 2012. doi:10.1111/j.1365-246X.2012.05498.x) who proposed a rupture scenario involving a reverse faulting along a north dipping fault and favoring a M _w 6.9 magnitude. In the present paper (1) we study the accuracy of their solutions in relation to the computational grid spacing and the dispersive/nondispersive parameterization, (2) based on an uncertainty on the recorded wave amplitude of the Genoa tide gauge they used, we propose a M _w 6.7 earthquake magnitude solution for the event (the kinematics is unchanged), co-existing with the M _w 6.9, (3) we evaluate the tsunami coastal impact of the 1887 event, and (4) we test a range of possible ruptures that local faults may undergo in order to propose a synoptic mapping of the tsunami threat in the area. The spatial distribution of the maximum wave height (MWH) is provided with a tentative identification of the processes that are responsible for it. This latter issue is imperative in order to make our mapping as generic as possible in the framework of our deterministic approach (based on realistic scenarios and not on ensemble statistics). The predictions suggest that the wave impact is mostly local, considering the relatively moderate size of the rupture planes. Although the present-day seismicity in this region is moderate, stronger earthquakes (M > 6.5) have occurred in the past. The studied scenarios show that for such events specific localities along the French–Italian Riviera may experience very significant MWH related to the shallow focal depth tested for such scenarios. We may reasonably conclude that the tsunami threat is relatively significant and uniform at the Italian side of the Riviera (from Ventimiglia to Imperia), while it is more localized (sporadic) at the French side from Antibes to Menton with, however, higher local level of inundation, e.g., Nice city center.     

A kernel spatial density estimation allowing for the analysis of spatial clustering. Application to Monsoon Asia Drought Atlas data

A nonparametric density estimate that incorporates spatial dependency has not been studied in the literature. In this article, we propose a new spatial density estimator that depends on two kernels: one controls the distance between observations while the other controls the spatial dependence structure. The uniform almost sure convergence of the density estimate is established with the rate of convergence. The consistency of the mode of this kernel density is also studied. Then a spatial hierarchical unsupervised clustering algorithm based on the mode estimate is presented. Some simulations as well as an application to the Monsoon Asia Drought Atlas data illustrate the efficiency of our algorithm, and a comparison of the spatial structures of these data detected by the density estimate and clustering algorithm are done.