Nutrients and zooplankton composition and dynamics in relation to the hydrological pattern in a confined Mediterranean salt marsh (NE Iberian Peninsula)

The role of the hydrological regime in the nutrients and zooplankton composition and dynamics has been analysed in five lagoons of La Pletera salt marshes (NE Iberian Peninsula) during a complete hydrological cycle (2002–2003). Two of the lagoons have their origin in the old river mouths while the other three were recently created in the framework of a Life Restoration project. This fact has also allowed us to study the effect of the lagoon age on nutrient and zooplankton composition and dynamics. The salt marsh hydrology is determined by a prolonged period of confinement without water inputs, irregularly interrupted by sudden water inputs due to flooding events (sea storms or intense rainfalls). While the dynamics of oxidized nitrogen compounds in the lagoons depends on the water inputs variability within each hydrological cycle, the internal load of phosphorus, total nitrogen and organic matter is related more to the cumulative mechanisms during the confinement periods. Accumulation processes may be easily related to lagoon age, since old lagoons have higher content of nutrients and organic matter, suggesting that these lagoons progressively accumulate nutrients during the successive confinement events. This is the usual case for most Mediterranean salt marshes without an artificially manipulated water regime. The zooplankton community in La Pletera integrates the effects of both the hydrological regime and the lagoon age since the former determines the temporal pattern of the main zooplankton species and the latter explains differences in composition and structure between old and new lagoons.     

The electrical structure of the Slave craton

The Slave craton in northwestern Canada, a relatively small Archean craton (600×400 km), is ideal as a natural laboratory for investigating the formation and evolution of Mesoarchean and Neoarchean sub-continental lithospheric mantle (SCLM). Excellent outcrop and the discovery of economic diamondiferous kimberlite pipes in the centre of the craton during the early 1990s have led to an unparalleled amount of geoscientific information becoming available.

Over the last 5 years deep-probing electromagnetic surveys were conducted on the Slave, using the natural-source magnetotelluric (MT) technique, as part of a variety of programs to study the craton and determine its regional-scale electrical structure. Two of the four types of surveys involved novel MT data acquisition; one through frozen lakes along ice roads during winter, and the second using ocean-bottom MT instrumentation deployed from float planes.

The primary initial objective of the MT surveys was to determine the geometry of the topography of the lithosphere–asthenosphere boundary (LAB) across the Slave craton. However, the MT responses revealed, completely serendipitously, a remarkable anomaly in electrical conductivity in the SCLM of the central Slave craton. This Central Slave Mantle Conductor (CSMC) anomaly is modelled as a localized region of low resistivity (10–15 Ω m) beginning at depths of 80–120 km and striking NE–SW. Where precisely located, it is spatially coincident with the Eocene-aged kimberlite field in the central part of the craton (the so-called “Corridor of Hope”), and also with a geochemically defined ultra-depleted harzburgitic layer interpreted as oceanic or arc-related lithosphere emplaced during early tectonism. The CSMC lies wholly within the NE–SW striking central zone defined by Grütter et al. [Grütter, H.S., Apter, D.B., Kong, J., 1999. Crust–mantle coupling; evidence from mantle-derived xenocrystic garnets. Contributed paper at: The 7th International Kimberlite Conference Proceeding, J.B. Dawson Volume, 1, 307–313] on the basis of garnet geochemistry (G10 vs. G9) populations.

Deep-probing MT data from the lake bottom instruments infer that the conductor has a total depth-integrated conductivity (conductance) of the order of 2000 Siemens, which, given an internal resistivity of 10–15 Ω m, implies a thickness of 20–30 km. Below the CSMC the electrical resistivity of the lithosphere increases by a factor of 3–5 to values of around 50 Ω m. This change occurs at depths consistent with the graphite–diamond transition, which is taken as consistent with a carbon interpretation for the CSMC.

Preliminary three-dimensional MT modelling supports the NE–SW striking geometry for the conductor, and also suggests a NW dip. This geometry is taken as implying that the tectonic processes that emplaced this geophysical–geochemical body are likely related to the subduction of a craton of unknown provenance from the SE (present-day coordinates) during 2630–2620 Ma. It suggests that the lithospheric stacking model of Helmstaedt and Schulze [Helmstaedt, H.H., Schulze, D.J., 1989. Southern African kimberlites and their mantle sample: implications for Archean tectonics and lithosphere evolution. In Ross, J. (Ed.), Kimberlites and Related Rocks, Vol. 1: Their Composition, Occurrence, Origin, and Emplacement. Geological Society of Australia Special Publication, vol. 14, 358–368] is likely correct for the formation of the Slave's current SCLM.     

Disjunctive Kriging with Hard and Imprecise Data

This paper presents a methodology for assessing local probability distributions by disjunctive kriging when the available data set contains some imprecise measurements, like noisy or soft information or interval constraints. The basic idea consists in replacing the set of imprecise data by a set of pseudohard data simulated from their posterior distribution; an iterative algorithm based on the Gibbs sampler is proposed to achieve such a simulation step. The whole procedure is repeated many times and the final result is the average of the disjunctive kriging estimates computed from each simulated data set. Being data-independent, the kriging weights need to be calculated only once, which enables fast computing. The simulation procedure requires encoding each datum as a pre-posterior distribution and assuming a Markov property to allow the updating of pre-posterior distributions into posterior ones. Although it suffers some imperfections, disjunctive kriging turns out to be a much more flexible approach than conditional expectation, because of the vast class of models that allows its computation, namely isofactorial models.     

Linking multi-level governance to local common-pool resource theory using fuzzy-set qualitative comparative analysis: Insights from twenty years of biodiversity conservation in Costa Rica

Understanding the relationship between multi-level institutional linkages and conditions influencing the likelihood of successful collective action has practical and theoretical relevance to sustainable local resource governance. This paper studies the relationship between multi-level linkages and local autonomy, a facilitating condition found to increase the likelihood of local successful collective action. A technique known as fuzzy-set qualitative comparative analysis (fsQCA) was applied to a longitudinal comparative data set. In the context of the decentralization of a protected area system in Costa Rica (1986–2006), it traced the emergence and endurance of autonomy among local institutions for biodiversity conservation. The technique illustrates which set of multi-level linkages combined in different ways, and at different points in time, to reach the same outcome (local autonomy). The findings show that a unique set of combinations of multi-level linkages led to the emergence of local autonomy among institutions for biodiversity conservation governance. In contrast, a more diverse set was associated with the endurance of local autonomy over time, suggesting that institutional diversity may play a more prominent role in the maintenance of institutional robustness than in processes of institutional emergence.     

Spatiotemporal patterns of burned areas,fire drivers,and fire probability across the equatorial Andes

Field-based fire studies in the equatorial Andes indicate that fires are strongly associated with biophysical and anthropogenic variables. However, fire controls and fire regimes at the regional scale remain undocumented. Therefore, this paper describes spatial and temporal burned-area patterns, identifies biophysical and anthropogenic fire drivers, and quantifies fire probability across 6° of latitude and 3° of longitude in the equatorial Andes. The spatial and temporal burned-area analysis was carried out based on 18 years(2001-2018) of the MCD64 A1 MODIS burned-area product. Climate, topography, vegetation, and anthropogenic variables were integrated in a logistic regression model to identify the significance of explanatory variables and determine fire occurrence probability. A total of 5779 fire events were registered during the 18 years of this study, located primarily along the western cordillera of the Andes and spreading from North to South. Eighty-eight percent of these fires took place within two fire hotspots located in the northwestern and southwestern corners of the study area. Ninety-nine percent occurred during the second part of the year, between June and December. The largest density of fires was primarily located on herbaceous vegetation and shrublands. Results show that mean monthly temperature, precipitation and NDVI during the prefire season, the location of land cover classes such as forest and agriculture, distance to roads and urban areas, slope, and aspect were the most important determinants of spatial and temporal fire distribution. The logistic regression model achieved a good accuracy in predicting fire probability(80%). Probability was higher in the southwestern and northern corners of the study area, and lower towards the north in the western and eastern piedmonts of the Andes. This analysis contributes to the understanding of fires in mountains within the tropics. The results here presented have the potential to contribute to fire management and control in the region.     

Spurious diapycnal mixing in terrain-following coordinate models: The problem and a solution

In this paper, we identify a crucial numerical problem in sigma coordinate models, leading to unacceptable spurious diapycnal mixing. This error is a by-product of recent advances in numerical methods, namely the implementation of high-order diffusive advection schemes. In the case of ROMS, spurious mixing is produced by its third-order upwind advection scheme, but our analysis suggests that all diffusive advection schemes would behave similarly in all sigma models. We show that the common idea that spurious mixing decreases with resolution is generally false. In a coarse-resolution regime, spurious mixing increases as resolution is refined, and may reach its peak value when eddy-driven lateral mixing becomes explicitly resolved. At finer resolution, diffusivities are expected to decrease but with values that only become acceptable at resolutions finer than the kilometer. The solution to this problem requires a specifically designed advection scheme. We propose and validate the RSUP3 scheme, where diffusion is split from advection and is represented by a rotated biharmonic diffusion scheme with flow-dependent hyperdiffusivity satisfying the Peclet constraint. The rotated diffusion operator is designed for numerical stability, which includes improvements of linear stability limits and a clipping method adapted to the sigma-coordinate. Realistic model experiments in a southwest Pacific configuration show that RSUP3 is able to preserve low dispersion and diffusion capabilities of the original third-order upwind scheme, while preserving water mass characteristics. There are residual errors from the rotated diffusion operator, but they remain acceptable. The use of a constant diffusivity rather than the Peclet hyperdiffusivity tends to increase these residual errors which become unacceptable with Laplacian diffusion. Finally, we have left some options open concerning the use of time filters as an alternative to spatial diffusion. A temporal discretization approach to the present problem (including implicit discretization) will be reported in a following paper.     

Change of Support for Estimating Local Block Grade Distributions

An important aspect in mineral resource evaluation is the reduction of variance when post-processing the grade distributions defined on the support (volume) of the available data into distributions defined on the support of the proposed selective mining units. Although the volume-variance relationship is well understood for the estimation of global grade distributions, it is still an unsolved issue for local estimation studies based on non-parametric geostatistical methods, such as indicator kriging, for which the support correction is not inherent to the method. To clarify this relationship, the local change of support problem is examined in the scope of two parametric models (multi-Gaussian and discrete Gaussian models). It is shown that the variance reduction factor between point and block-support local distributions depends on the block being considered and is less than the global variance reduction factor. As a consequence, post-processing the local point-support grade distributions on the basis of the latter systematically understates the importance of the change of support at the local scale and makes selective mining appear more economically attractive than it really is. In the light of these results, a methodology is proposed to post-process the local point-support distributions obtained via non-parametric (indicator) methods into block-support distributions. An application to simulated data indicates that this methodology provides an accurate estimation at the block support when dealing with diffusion-type random fields.     

Estimates of groundwater discharge to a coastal wetland using multiple techniques: Taylor Slough, Everglades National Park, USA

Quantifying water exchange between a coastal wetland and the underlying groundwater is important for closing water, energy and chemical budgets. The coastal wetlands of the Florida Everglades (USA) are at the forefront of a large hydrologic restoration project, and understanding of groundwater/surface-water interactions is needed to comprehend the effects of the project. Four independent techniques were used to identify water exchange at varying spatial and temporal scales in Taylor Slough, Everglades National Park. The techniques included a water-budget study and measurements of hydraulic head gradients, geochemical tracers, and temperature. During the 18-month study, the four methods converged as to the timing of groundwater discharge, typically between June and September, contemporaneous with the wet season and increasing surface-water levels. These results were unexpected, as groundwater discharge was predicted to be greatest when surface-water levels were low, typically during the dry season. Either a time lag of 1?C5?months in the response of groundwater discharge to low surface-water levels or precipitation-induced groundwater discharge may explain the results. Groundwater discharge was a significant contributor (27?%) to the surface water in Taylor Slough with greater rates of discharge observed towards the coastline in response to seawater intrusion.