A land-use-based modelling chain to assess the impacts of Natural Water Retention Measures on Europe’s Green Infrastructure

This article details the process of integrating models to answer a specific policy-driven question: ‘What could be the impact of proposed Natural Water Retention Measures (NWRMs) on Europe’s Green Infrastructure (GI)?’ It describes the new Land Use Modelling Platform (LUMP), now enabling a high spatial scale (100-m) and large coverage (pan-European), whereby several sector-specific models contribute to assessing the impact of regional-level policy on a given spatial topic of concern. The configuration (land claims and land allocations modules) and calibration (accessibility and biophysical suitability) of the LUMP are explained. Four NWRM scenarios (riparian areas, afforestation, grassland and baseline scenario) are configured to run the simulations. For the reference: year 2006, the spatial representation of GI is based on land-use features of a refined version of the CORINE Land Cover (CLC), and resumed as connected components made of nodes and links.

Mathematical morphological image processing and network graph theory model, available from the free software package GUIDOS (the Joint Research Center of the European Commission), enabled the measurement of the GI connectivity and identified most critical links. Results show that the competition for land claimed by different economic sectors, combined with policy-driven rule-sets for the implementation of different NWRMs, yields very different results for the 2030 land-use projections, and subsequently for the morphology of GI. Three indicators associated with the morphology of GI are computed in order to assess the model outputs for 2030. The indicators are computed to answer the following questions: (1) How is the quantity of GI affected by each of the NWRM, and what proportion of that GI is most valuable? (2) What is the location of the most critical nodes and connectors of GI, and what land-use conversions occur under these? (3) Are the average components getting larger or smaller?

Whereas the grassland measure results in the largest net increase of GI, the afforestation measure results in the overall largest number of hectares of key nodes and links within the network. Land conversions occur under the critical GI nodes and links, with a large increase in agricultural areas, especially for the riparian measure under critical nodes and the grassland measure under critical links. Also predominant is the swapping of land from pasture to forest under critical links with the afforestation measure. The riparian measure most increases the average size of GI components, and all three measures contribute to bridging two large components which were divided in the 2006 land-use map, thus increasing the size of the largest component by more than 50%.     

Spatial dimensions of New Zealand's environmental management

New Zealand's regional council boundaries drawn as part of the 1989 local government reforms privilege integrated catchment management. This privilege is now challenged by calls to both fragment and amalgamate councils and centralise decision‐making but with little analytic assessment of the spatial differentiation to support them. The spatial underpinnings of these calls were explored using a case study that assessed the spatial dimensions of environment management. Analysis of the Land Environments and River Environment Classification databases and councils' planning documents highlights challenges local governments face in applying integrated natural resource management within single purpose jurisdictions and implications for wider sustainable development policy.     

Uniwavelength lidar sensitivity to spherical aerosol microphysical properties for the interpretation of Lagrangian stratospheric observations

The determination of stratospheric particle microphysical properties from multiwavelength lidar, including Rayleigh and/or Raman detection, has been widely investigated. However, most lidar systems are uniwavelength operating at 532 nm. Although the information content of such lidar data is too limited to allow the retrieval of the full size distribution, the coupling of two or more uniwavelength lidar measurements probing the same moving air parcel may provide some meaningful size information. Within the ORACLE-O3 IPY project, the coordination of several ground-based lidars and the CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation) space-borne lidar is planned during measurement campaigns called MATCH-PSC (Polar Stratospheric Clouds). While probing the same moving air masses, the evolution of the measured backscatter coefficient (BC) should reflect the variation of particles microphysical properties. A sensitivity study of 532 nm lidar particle backscatter to variations of particles size distribution parameters is carried out. For simplicity, the particles are assumed to be spherical (liquid) particles and the size distribution is represented with a unimodal log-normal distribution. Each of the four microphysical parameters (i.e. log-normal size distribution parameters, refractive index) are analysed separately, while the three others are remained set to constant reference values. Overall, the BC behaviour is not affected by the initial values taken as references. The total concentration (N₀) is the parameter to which BC is least sensitive, whereas it is most sensitive to the refractive index (m). A 2% variation of m induces a 15% variation of the lidar BC, while the uncertainty on the BC retrieval can also reach 15%. This result underlines the importance of having both an accurate lidar inversion method and a good knowledge of the temperature for size distribution retrieval techniques. The standard deviation (σ) is the second parameter to which BC is most sensitive to. Yet, the impact of m and σ on BC variations is limited by the realistic range of their variations. The mean radius (r_m) of the size distribution is thus the key parameter for BC, as it can vary several-fold. BC is most sensitive to the presence of large particles. The sensitivity of BC to r_m and σ variations increases when the initial size distributions are characterized by low r_m and large σ. This makes lidar more suitable to detect particles growing on background aerosols than on volcanic aerosols.     

The relationship between zooplankton,conductivity and lake-water ionic composition in 111 lakes from the Interior Plateau of British Columbia,Canada

Zooplankton collected from vertical net tows were related to the environmental variables from 98 lakes from the Interior Plateau of British Columbia. Canonical correspondence analysis showed that both salinity and ionic composition (pH and Mg) of the lake-water made major and significant contributions to the first two ordination axes (=0.42 and 0.11 respectively,P<0.05). BothArtemia franciscana andMoina hutchinsoni had their highest relative abundance in meso-hypersaline waters. However,Artemia franciscana preferred waters that were higher in Mg and Ca, whileMoina hutchinsoni was found in waters that were lower in Mg and Ca. Similarly, at intermediate salinities,Daphnia pulex and the calanoid copepods preferred waters slightly lower in Mg and Ca, whereasCeriodaphnia laticaudata andSimocephalus spp. were relatively more common in waters higher in Mg and Ca. Because the freshest lakes studied varied much less in ionic composition, the zooplankton in these lakes did not show a preference to ionic composition. As expected, multi-generic groups, such as the calanoid copepods, cyclopoid copepods and nauplii, had wider tolerances to conductivity than groups identified to lower taxonomic levels. Significant weighted-averaging regression and calibration models of conductivity were developed based on zooplankton species composition from the study lakes (r ²=0.56,P<0.05). Samples composed largely of multi-generic taxa yielded the worst estimates of salinity in the reconstruction model. This study suggests that zooplankton community composition may be developed into a useful proxy for paleosalinity reconstruction.     

Rapid and cyclic aeolian deposition during the Last Glacial in European loess: a high-resolution record from Nussloch, Germany

This paper reports the results of an investigation of the Weichselian Upper Pleniglacial loess sequences of Nussloch (Rhine Valley, Germany) based on stratigraphy, palaeopedology, sedimentology, palynology, malacology and geochemistry (δ¹³C), supported by radiocarbon, TL and OSL dating. Grain-size and magnetic susceptibility records are taken at 5 cm intervals from the Upper Pleniglacial (UPG) loess. The data indicate cyclic variations in loess deposition between ca 34 and 17 ka, when the sedimentation rate is especially high (1.0–1.2 m per ka for more than 10 m). The grain-size index (GSI: ratio of coarse silt versus fine silt and clay) shows variations, which are assumed to be an indirect measurement of wind intensity. The sedimentation rate, interpreted from the profiles, indicates high values in loess (Loess events LE-1 to LE-7) and low or negligible values in tundra gley horizons G1 to G8. OSL ages from the loess and ¹⁴C dates from organic matter in the loess show that loess deposition was rapid but was interrupted by shorter periods of reduced aeolian sedimentation. Comparison between the data from Nussloch and other European sequences demonstrates a progressive coarsening of the loess deposits between ca 30 and 22 ka. This coarsening trend ends with a short but major grain-size decrease and is followed by an increase to a new maximum at 20 ± 2 ka (“W” shape). Correlation between the loess GSI and the Greenland ice-core dust records, suggests a global connection between North Atlantic and Western European global atmospheric circulation and wind regimes. In addition, the typical Upper Pleniglacial loess deposition begins at ca 30–31 ka, close to Heinrich event (HE) 3, and the main period of loess sedimentation at about 25 ± 2 ka is coeval to HE 2. Correlation of magnetic susceptibility and grain-size records shows that the periods, characterised by high GSI, coincide with an increase in the amount of ferromagnetic minerals reworked from the Rhine alluvial plain. They suggest enhancement in the frequency of the storms from N–NW. These results are integrated within a palaeogeographical model of dust transport and deposition in Western Europe for the Weichselian Upper Pleniglacial (or Late Pleniglacial).     

New insights into the P- and S-wave velocity structure of the D" discontinuity beneath the Cocos plate

Broad-band P - and S -waves from earthquakes in South America recorded at Californian network stations are analysed to image lateral variations of the D"-discontinuity beneath the Cocos plate. We apply two array processing methods to the data set: a simplified migration method to the P -wave data set and a double-array method to both the P - and S -wave data sets, allowing us to compare results from the two methods. The double-array method images a dipping reflector at a depth range from 2650 to 2700 km in the southern part of the study area. We observe a step-like topography of 100 km to a shallower reflector at about 2600 km depth to the north, as well as evidence for a second (deeper) reflector at a depth range from 2700 to 2750 km in the north. Results from the simplified migration agree well with those from the double-array method, similarly locating a large step in reflector depth in a similar location (about 2650 km depth in the south and about 2550 km in the north) as well as the additional deeper reflector at the depth of about 2750 km in the north. Waveform modelling of the reflected waves from both methods suggests a positive velocity contrast for S waves, but a negative velocity contrast for P waves for the upper reflector in agreement with predictions from mineral physical calculations for a post-perovskite phase transition. The data also show some evidence for the existence of another deeper reflector that could indicate a double intersection of the geotherm with the post-perovskite stability field, that is, the back-transformation of post-perovskite to perovskite close to the core–mantle boundary.     

Slowness-backazimuth weighted migration: A new array approach to a high-resolution image

We have developed a new array method combining conventional migration with a slowness-backazimuth deviation weighting scheme. All seismic traces are shifted based on the theoretical traveltime of the scattered wave from specific gridpoints in a 3-D volume. Observed slowness and backazimuth are calculated for each raypath and compared with theoretical values in order to estimate slowness and backazimuth deviations. Subsequently, stacked energy calculated by a conventional migration method is weighted by the slowness and backazimuth deviations to suppress any arrival energy whose slowness and backazimuth are inconsistent with the expected theoretical values. This new method was applied to two P- wave data sets which comprise (1) underside reflections at the 410 and 660 km mantle discontinuities and (2) D" reflections as well as their corresponding synthetic data sets. The results show that the weighting scheme dramatically increases the resolution of the migrated images and enables us to obtain well-constrained, focused images, making upper-mantle discontinuities and D" reflections more distinct by reducing their surrounding energy.     

Localization and characterization of an active fault in an urbanized area in central Guatemala by means of geoelectrical imaging

The Polochic and Motagua faults define the active plate boundary between the North American and Caribbean plates in central Guatemala. A splay of the Polochic Fault traverses the rapidly growing city of San Miguel Uspantán that is periodically affected by destructive earthquakes. This fault splay was located using a 2D electrical resistivity tomography (ERT) survey that also characterized the fault damage zone and evaluated the thickness and nature of recent deposits upon which most of the city is built. ERT images show the fault as a ~50 m wide, near-vertical low-resistivity anomaly, bounded within a few meters by high resistivity anomalies. Forward modeling reproduces the key aspects of the observed electrical resistivity data with remarkable fidelity thus defining the overall location, geometry, and internal structure of the fault zone as well as the affected lithologies. Our results indicate that the city is constructed on a ~20 m thick surficial layer consisting of poorly consolidated, highly porous, water-logged pumice. This soft layer is likely to amplify seismic waves and to liquefy upon moderate to strong ground shaking. The electrical conductivity as well as the major element chemistry of the groundwater provides evidence to suggest that the local aquifer might, at least in part, be fed by water rising along the fault. Therefore, the potential threat posed by this fault splay may not be limited to its seismic activity per se, but could be compounded its potential propensity to enhance seismic site effects by injecting water into the soft surficial sediments. The results of this study provide the basis for a rigorous analysis of seismic hazard and sustainable development of San Miguel Uspantán and illustrate the potential of ERT surveying for paleoseismic studies.     

Optimal water depth management on river-fed National Wildlife Refuges in a changing climate

The prairie pothole region (PPR) in the north-central United States and south-central Canada constitutes the most important waterfowl breeding area in North America. Projected long-term changes in precipitation and temperature may alter the drivers of waterfowl abundance: wetland availability and emergent vegetation cover. Previous studies have focused on isolated wetland dynamics, but the implications of changing precipitation on managed, river-fed wetlands have not been addressed. Using a structured decision making (SDM) approach, we derived optimal water management actions for 20 years at four river-fed National Wildlife Refuges (NWRs) in North and South Dakota under contrasting increasing/decreasing (+/?0.4 %/year) inflow scenarios derived from empirical trends. Refuge pool depth is manipulated by control structures. Optimal management involves setting control structure heights that have the highest probability of providing a desired mix of waterfowl habitat, given refuge capacities and inflows. We found optimal seasonal control structure heights for each refuge were essentially the same under increasing and decreasing inflow trends of 0.4 %/year over the next 20 years. Results suggest managed pools in the NWRs receive large inflows relative to their capacities. Hence, water availability does not constrain management; pool bathymetry and management tactics can be greater constraints on attaining management objectives than climate-mediated inflow. We present time-dependent optimal seasonal control structure heights for each refuge, which are resilient to the non-stationary precipitation scenarios we examined. Managers can use this information to provide a desired mixture of wildlife habitats, and to re-assess management objectives in reserves where pool bathymetry prevents attaining the currently stated objectives.     

Using Statistical Downscaling to Quantify the GCM-Related Uncertainty in Regional Climate Change Scenarios： A Case Study of Swedish Precipitation

There are a number of sources of uncertainty in regional climate change scenarios. When statistical downscaling is used to obtain regional climate change scenarios, the uncertainty may originate from the uncertainties in the global climate models used, the skill of the statistical model, and the forcing scenarios applied to the global climate model. The uncertainty associated with global climate models can be evaluated by examining the differences in the predictors and in the downscaled climate change scenarios based on a set of different global climate models. When standardized global climate model simulations such as the second phase of the Coupled Model Intercomparison Project （CMIP2） are used, the difference in the downscaled variables mainly reflects differences in the climate models and the natural variability in the simulated climates. It is proposed that the spread of the estimates can be taken as a measure of the uncertainty associated with global climate models. The proposed method is applied to the estimation of global-climate-model-related uncertainty in regional precipitation change scenarios in Sweden. Results from statistical downscaling based on 17 global climate models show that there is an overall increase in annual precipitation all over Sweden although a considerable spread of the changes in the precipitation exists. The general increase can be attributed to the increased large-scale precipitation and the enhanced westerly wind. The estimated uncertainty is nearly independent of region. However, there is a seasonal dependence. The estimates for winter show the highest level of confidence, while the estimates for summer show the least.