Uncertainty assessment of a process-based integrated catchment model of phosphorus

Despite the many models developed for phosphorus concentration prediction at differing spatial and temporal scales, there has been little effort to quantify uncertainty in their predictions. Model prediction uncertainty quantification is desirable, for informed decision-making in river-systems management. An uncertainty analysis of the process-based model, integrated catchment model of phosphorus (INCA-P), within the generalised likelihood uncertainty estimation (GLUE) framework is presented. The framework is applied to the Lugg catchment (1,077 km2), a River Wye tributary, on the England–Wales border. Daily discharge and monthly phosphorus (total reactive and total), for a limited number of reaches, are used to initially assess uncertainty and sensitivity of 44 model parameters, identified as being most important for discharge and phosphorus predictions. This study demonstrates that parameter homogeneity assumptions (spatial heterogeneity is treated as land use type fractional areas) can achieve higher model fits, than a previous expertly calibrated parameter set. The model is capable of reproducing the hydrology, but a threshold Nash-Sutcliffe co-efficient of determination (E or R 2) of 0.3 is not achieved when simulating observed total phosphorus (TP) data in the upland reaches or total reactive phosphorus (TRP) in any reach. Despite this, the model reproduces the general dynamics of TP and TRP, in point source dominated lower reaches. This paper discusses why this application of INCA-P fails to find any parameter sets, which simultaneously describe all observed data acceptably. The discussion focuses on uncertainty of readily available input data, and whether such process-based models should be used when there isn’t sufficient data to support the many parameters.

Insights into redox cycling of sulfur and iron in peatlands using high-resolution diffusive equilibrium thin film (DET) gel probe sampling

High spatial resolution vertical profiles of pore-water chemistry have been obtained for a peatland using diffusive equilibrium in thin films (DET) gel probes. Comparison of DET pore-water data with more traditional depth-specific sampling shows good agreement and the DET profiling method is less invasive and less likely to induce mixing of pore-waters. Chloride mass balances as water tables fell in the early summer indicate that evaporative concentration dominates and there is negligible lateral flow in the peat. Lack of lateral flow allows element budgets for the same site at different times to be compared. The high spatial resolution of sampling also enables gradients to be observed that permit calculations of vertical fluxes. Sulfate concentrations fall at two sites with net rates of 1.5 and 5.0nmol cm− 3 day− 1, likely due to a dominance of bacterial sulfate reduction, while a third site showed a net gain in sulfate due to oxidation of sulfur over the study period at an average rate of 3.4nmol cm− 3 day− 1. Behaviour of iron is closely coupled to that of sulfur; there is net removal of iron at the two sites where sulfate reduction dominates and addition of iron where oxidation dominates. The profiles demonstrate that, in addition to strong vertical redox related chemical changes, there is significant spatial heterogeneity. Whilst overall there is evidence for net reduction of sulfate within the peatland pore-waters, this can be reversed, at least temporarily, during periods of drought when sulfide oxidation with resulting acid production predominates.

A data assimilation method for using low-resolution Earth observation data in heterogeneous ecosystems

We present an approach for dealing with coarse-resolution Earth observations (EO) in terrestrial ecosystem data assimilation schemes. The use of coarse-scale observations in ecological data assimilation schemes is complicated by spatial heterogeneity and nonlinear processes in natural ecosystems. If these complications are not appropriately dealt with, then the data assimilation will produce biased results. The “disaggregation” approach that we describe in this paper combines frequent coarse-resolution observations with temporally sparse fine-resolution measurements. We demonstrate the approach using a demonstration data set based on measurements of an Arctic ecosystem. In this example, normalized difference vegetation index observations are assimilated into a “zero-order” model of leaf area index and carbon uptake. The disaggregation approach conserves key ecosystem characteristics regardless of the observation resolution and estimates the carbon uptake to within 1% of the demonstration data set “truth.” Assimilating the same data in the normal manner, but without the disaggregation approach, results in carbon uptake being underestimated by 58% at an observation resolution of 250 m. The disaggregation method allows the combination of multiresolution EO and improves in spatial resolution if observations are located on a grid that shifts from one observation time to the next. Additionally, the approach is not tied to a particular data assimilation scheme, model, or EO product and can cope with complex observation distributions, as it makes no implicit assumptions of normality.

Cloud bands over southern Africa: seasonality, contribution to rainfall variability, and modulation by the MJO

Tropical-extratropical cloud band systems over southern Africa, known as tropical temperate troughs (TTTs), are known to contribute substantially to South African summer rainfall. This study performs a comprehensive assessment of the seasonal cycle and rainfall contribution of TTTs by using a novel object-based strategy that explicitly tracks these systems for their full life cycle. The methodology incorporates a simple assignment of station rainfall data to each event, thereby creating a database containing detailed rainfall characteristics for each TTT. This is used to explore the importance of TTTs for rain days and climatological rainfall totals in October–March. Average contributions range from 30 to 60 % with substantial spatial heterogeneity observed. TTT rainfall contributions over the Highveld and eastern escarpment are lower than expected. A short analysis of TTT rainfall variability indicates TTTs provide substantial, but not dominant, intraseasonal and interannual variability in station rainfall totals. TTTs are however responsible for a high proportion of heavy rainfall days. Of 52 extreme rainfall events in the 1979–1999 period, 30 are associated with these tropical-extratropical interactions. Cut-off lows were included in the evolution of 6 of these TTTs. The study concludes with an analysis of the question: does the Madden-Julian Oscillation influence the intensity of TTT rainfall over South Africa? Results suggest a weak but significant suppression (enhancement) of intensity during phase 1(6).

How will organic carbon stocks in mineral soils evolve under future climate? Global projections using RothC for a range of climate change scenarios

We use a soil carbon (C) model (RothC), driven by a range of climate models for a range of climate scenarios to examine the impacts of future climate on global soil organic carbon (SOC) stocks. The results suggest an overall global increase in SOC stocks by 2100 under all scenarios, but with a different extent of increase among the climate model and emissions scenarios. The impacts of projected land use changes are also simulated, but have relatively minor impacts at the global scale. Whether soils gain or lose SOC depends upon the balance between C inputs and decomposition. Changes in net primary production (NPP) change C inputs to the soil, whilst decomposition usually increases under warmer temperatures, but can also be slowed by decreased soil moisture. Underlying the global trend of increasing SOC under future climate is a complex pattern of regional SOC change. SOC losses are projected to occur in northern latitudes where higher SOC decomposition rates due to higher temperatures are not balanced by increased NPP, whereas in tropical regions, NPP increases override losses due to higher SOC decomposition. The spatial heterogeneity in the response of SOC to changing climate shows how delicately balanced the competing gain and loss processes are, with subtle changes in temperature, moisture, soil type and land use, interacting to determine whether SOC increases or decreases in the future. Our results suggest that we should stop looking for a single answer regarding whether SOC stocks will increase or decrease under future climate, since there is no single answer. Instead, we should focus on improving our prediction of the factors that determine the size and direction of change, and the land management practices that can be implemented to protect and enhance SOC stocks.

Diagnosing atmosphere–land feedbacks in CMIP5 climate models

Human-made transformations to the environment, and in particular the land surface, are having a large impact on the distribution (in both time and space) of rainfall, upon which all life is reliant. Focusing on precipitation, soil moisture and near-surface temperature, we compare data from Phase 5 of the Climate Modelling Intercomparison Project (CMIP5), as well as blended observational–satellite data, to see how the interaction between rainfall and the land surface differs (or agrees) between the models and reality, at daily timescales. As expected, the results suggest a strong positive relationship between precipitation and soil moisture when precipitation leads and is concurrent with soil moisture estimates, for the tropics as a whole. Conversely a negative relationship is shown when soil moisture leads rainfall by a day or more. A weak positive relationship between precipitation and temperature is shown when either leads by one day, whereas a weak negative relationship is shown over the same time period between soil moisture and temperature. Temporally, in terms of lag and lead relationships, the models appear to be in agreement on the overall patterns of correlation between rainfall and soil moisture. However, in terms of spatial patterns, a comparison of these relationships across all available models reveals considerable variability in the ability of the models to reproduce the correlations between precipitation and soil moisture. There is also a difference in the timings of the correlations, with some models showing the highest positive correlations when precipitation leads soil moisture by one day. Finally, the results suggest that there are 'hotspots' of high linear gradients between precipitation and soil moisture, corresponding to regions experiencing heavy rainfall. These results point to an inability of the CMIP5 models to simulate a positive feedback between soil moisture and precipitation at daily timescales. Longer timescale comparisons and experiments at higher spatial resolutions, where the impact of the spatial heterogeneity of rainfall on the initiation of convection and supply of moisture is included, would be expected to improve process understanding further.

Statistical–dynamical downscaling of present day and future precipitation regimes in the Aksu river catchment in Central Asia

A statistical–dynamical downscaling (SDD) approach is applied to determine present day and future high-resolution rainfall distributions in the catchment of the river Aksu at the southern slopes of the Tienshan Mountains, Central Asia. First, a circulation weather type (CWT) classification is employed to define typical lower atmospheric flow regimes from ERA-40 reanalysis data. Selected representatives of each CWT are dynamically downscaled with the regional climate model COSMO-CLM 4.8 at a horizontal grid resolution of 0.0625°, using the ERA-40 reanalysis data as boundary conditions. Finally, the simulated representatives are recombined to obtain a high-resolution rainfall climatology for present day climate. The methodology is also applied to ensemble simulations of three different scenarios of the global climate model ECHAM5/MPI-OM1 to derive projections of rainfall changes until 2100. Comparisons of downscaled seasonal and annual rainfall with observational data suggest that the statistical–dynamical approach is appropriate to capture the observed present-day precipitation climatology over the low lands and the first elevations of the Tienshan Mountains. On the other hand, a strong bias is found at higher altitudes, where precipitation is clearly underestimated by SDD. The application of SDD to the ECHAM5/MPI-OM1 ensemble reveals that precipitation changes by the end of the 21st century depend on the season. While for autumn an increase of seasonal precipitation is found for all simulations, a decrease in precipitation is obtained during winter for most parts of the Aksu catchment. The spread between different ECHAM5/MPI-OM1 ensemble members is strongest in spring, where trends of opposite sign are found. The largest changes in rainfall are simulated for the summer season, which also shows the most pronounced spatial heterogeneity. Most ECHAM5/MPI-OM1 realizations indicate a decrease of annual precipitation over large parts of the Tienshan, and an increase restricted to the southeast of the study area. These results provide a good basis for downscaling present-day and future rainfall distributions for hydrological purposes.

Dynamical and observational constraints on tropical Pacific sea surface temperatures at the last glacial maximum

Asynchronously coupled atmosphere and ocean general circulation model simulations are used to examine the consequences of changes in the west/east sea-surface temperature (SST) gradient across the equatorial Pacific at the last glacial maximum (LGM). Simulations forced by the CLIMAP SST for the LGM, where the west/east SST gradient across the Pacific is reduced compared to present, produce a reduction in the strength of the trade winds and a decrease in the west/east slope of the equatorial thermocline that is incompatible with thermocline depths newly inferred from foraminiferal assemblages. Stronger-than-present trade winds, and a more realistic simulation of the thermocline slope, are produced when eastern Pacific SSTs are 2°C cooler than western Pacific SSTs. Our study highlights the importance of spatial heterogeneity in tropical SSTs in determining key features of the glacial climate.

Evidence of a topographic signal in surface soil moisture derived from ENVISAT ASAR wide swath data

The susceptibility of a catchment to flooding is affected by its soil moisture prior to an extreme rainfall event. While soil moisture is routinely observed by satellite instruments, results from previous work on the assimilation of remotely sensed soil moisture into hydrologic models have been mixed. This may have been due in part to the low spatial resolution of the observations used. In this study, the remote sensing aspects of a project attempting to improve flow predictions from a distributed hydrologic model by assimilating soil moisture measurements are described. Advanced Synthetic Aperture Radar (ASAR) Wide Swath data were used to measure soil moisture as, unlike low resolution microwave data, they have sufficient resolution to allow soil moisture variations due to local topography to be detected, which may help to take into account the spatial heterogeneity of hydrological processes. Surface soil moisture content (SSMC) was measured over the catchments of the Severn and Avon rivers in the South West UK. To reduce the influence of vegetation, measurements were made only over homogeneous pixels of improved grassland determined from a land cover map. Radar backscatter was corrected for terrain variations and normalized to a common incidence angle. SSMC was calculated using change detection. To search for evidence of a topographic signal, the mean SSMC from improved grassland pixels on low slopes near rivers was compared to that on higher slopes. When the mean SSMC on low slopes was 30–90%, the higher slopes were slightly drier than the low slopes. The effect was reversed for lower SSMC values. It was also more pronounced during a drying event. These findings contribute to the scant information in the literature on the use of high resolution SAR soil moisture measurement to improve hydrologic models.

Changes in fire regimes since the Last Glacial Maximum: an assessment based on a global synthesis and analysis of charcoal data

Fire activity has varied globally and continuously since the last glacial maximum (LGM) in response to long-term changes in global climate and shorter-term regional changes in climate, vegetation, and human land use. We have synthesized sedimentary charcoal records of biomass burning since the LGM and present global maps showing changes in fire activity for time slices during the past 21,000 years (as differences in charcoal accumulation values compared to pre-industrial). There is strong broad-scale coherence in fire activity after the LGM, but spatial heterogeneity in the signals increases thereafter. In North America, Europe and southern South America, charcoal records indicate less-than-present fire activity during the deglacial period, from 21,000 to ∼11,000 cal yr BP. In contrast, the tropical latitudes of South America and Africa show greater-than-present fire activity from ∼19,000 to ∼17,000 cal yr BP and most sites from Indochina and Australia show greater-than-present fire activity from 16,000 to ∼13,000 cal yr BP. Many sites indicate greater-than-present or near-present activity during the Holocene with the exception of eastern North America and eastern Asia from 8,000 to ∼3,000 cal yr BP, Indonesia and Australia from 11,000 to 4,000 cal yr BP, and southern South America from 6,000 to 3,000 cal yr BP where fire activity was less than present. Regional coherence in the patterns of change in fire activity was evident throughout the post-glacial period. These complex patterns can largely be explained in terms of large-scale climate controls modulated by local changes in vegetation and fuel load

Potentially bioavailable metals in sediment from a tropical polymictic environment-Rio Grande Reservoir, Brazil

Background, aim and scope Although many recent studies have focused on sediment potential toxicity, few of them were performed in tropical shallow aquatic environments. Those places can suffer short-time variations, especially due to water column circulations generated by changes in temperature and wind. Rio Grande reservoir is such an example; aside from that, it suffers various anthropogenic impacts, despite its multiple uses. Materials and methods This work presents the first screening step for understanding sediment quality from Rio Grande reservoir by comparing metal content using three different sediment quality guidelines. We also aimed at verifying any possible spatial heterogeneity. Results and discussion We found spatial heterogeneity varying according to the specific metal. Results showed a tendency for metals to remain as insoluble as metal sulfide (potentially not bioavailable), since sulfide was in excess and sediment physical-chemical characteristics contribute to sulfide maintenance (low redox potential, neutral pH, low dissolved oxygen, and high organic matter content). On the other hand, metal concentrations were much higher than suggested by Canadian guidelines and regional background values, especially Cu, which raises the risk of metal remobilization in cases of water circulation. Further study steps include the temporal evaluation of AVS/SEM, a battery of bioassays and the characterization of organic compounds.

Fine-scale adaptation in a clonal sea anemone

Local adaptation in response to fine-scale spatial heterogeneity is well documented in terrestrial ecosystems. In contrast, in marine environments local adaptation has rarely been documented or rigorously explored. This may reflect real or anticipated effects of genetic homogenization, resulting from widespread dispersal in the sea. However, evolutionary theory predicts that for the many benthic species with complex life histories that include both sexual and asexual phases, each parental habitat patch should become dominated by the fittest and most competitive clones. In this study we used genotypic mapping to show that within headlands, clones of the sea anemone Actinia tenebrosa show restricted distributions to specific habitats despite the potential for more widespread dispersal. On these same shores we used reciprocal transplant experiments that revealed strikingly better performance of clones within their natal rather than foreign habitats as judged by survivorship, asexual fecundity, and growth. These findings highlight the importance of selection for fine-scale environmental adaptation in marine taxa and imply that the genotypic structure of populations reflects extensive periods of interclonal competition and site-specific selection.

Spatio-temporal analysis of house price convergence based on a demographical distance

Convergences of house prices have been studied for over three decades, but yet have been confirmed because of spatial heterogeneity and autocorrelations in house prices. A spatio-temporal approach was recently proposed to address the spatial and temporal issues related to house prices. However, most previous studies placed the focus on the spatial heterogeneity and autocorrelations from geographical locations, which neglected other spatial factors. In order to overcome this shortfall, this research argued a demographical distance, constructed by demographical structure and housing market scales, to investigate the house price convergences in Australian capital cities. The results confirmed the house price levels in Canberra, Brisbane and Perth converged to the house price level in Sydney.

Variability of sedimentary organic carbon in patchy seagrass landscapes

Seagrass ecosystems, considered among the most efficient carbon sinks worldwide, encompass a wide variety of spatial configurations in the coastal landscape. Here we evaluated the influence of the spatial configuration of seagrass meadows at small scales (metres) on carbon storage in seagrass sediments. We intensively sampled carbon stocks and other geochemical properties (δ(13)C, particle size, depositional fluxes) across seagrass-sand edges in a Zostera muelleri patchy seagrass landscape. Carbon stocks were significantly higher (ca. 20%) inside seagrass patches than at seagrass-sand edges and bare sediments. Deposition was similar among all positions and most of the carbon was from allochthonous sources. Patch level attributes (e.g. edge distance) represent important determinants of the spatial heterogeneity of carbon stocks within seagrass ecosystems. Our findings indicate that carbon stocks of seagrass areas have likely been overestimated by not considering the influence of meadow landscapes, and have important relevance for the design of seagrass carbon stock assessments.

Regional subsidies and industrial prospects of lagging regions

Large and sustained differences in economic performance across regions of developing countries have long provided motivation for fiscal incentives designed to encourage firm entry in lagging areas. Empirical evidence in support of these policies has, however, been weak at best. This paper undertakes a direct evaluation of the most prominent fiscal incentive policy in Brazil, the Fundos Constitucionais de Financiamento (Constitutional Funds). In doing so, we exploit valuable features of the Brazilian Ministry of Labor's RAIS data set to address two important elements of firm location decisions that have the potential to bias an assessment of the Funds: (i) firm “family structure” (in particular, proximity to headquarters for vertically integrated firms), and (ii) unobserved spatial heterogeneity (with the potential to confound the effects of the Funds). We find that the pull of firm headquarters is very strong relative to the Constitutional Funds for vertically integrated firms, but that, with non-parametric controls for time invariant spatial heterogeneity, the Funds provide statistically and economically significant incentives for firms in many of the targeted industries.