Species Richness-Environment Relationships of European Arthropods at Two Spatial Grains: Habitats and Countries

We study how species richness of arthropods relates to theories concerning net primary productivity, ambient energy, water-energy dynamics and spatial environmental heterogeneity. We use two datasets of arthropod richness with similar spatial extents (Scandinavia to Mediterranean), but contrasting spatial grain (local habitat and country). Samples of ground-dwelling spiders, beetles, bugs and ants were collected from 32 paired habitats at 16 locations across Europe. Species richness of these taxonomic groups was also determined for 25 European countries based on the Fauna Europaea database. We tested effects of net primary productivity (NPP), annual mean temperature (T), annual rainfall (R) and potential evapotranspiration of the coldest month (PETmin) on species richness and turnover. Spatial environmental heterogeneity within countries was considered by including the ranges of NPP, T, R and PETmin. At the local habitat grain, relationships between species richness and environmental variables differed strongly between taxa and trophic groups. However, species turnover across locations was strongly correlated with differences in T. At the country grain, species richness was significantly correlated with environmental variables from all four theories. In particular, species richness within countries increased strongly with spatial heterogeneity in T. The importance of spatial heterogeneity in T for both species turnover across locations and for species richness within countries suggests that the temperature niche is an important determinant of arthropod diversity. We suggest that, unless climatic heterogeneity is constant across sampling units, coarse-grained studies should always account for environmental heterogeneity as a predictor of arthropod species richness, just as studies with variable area of sampling units routinely consider area.

Limitations of Remotely-sensed Aerosol as a Spatial Proxy for Fine Particulate Matter

Recent research highlights the promise of remotely-sensed aerosol optical depth (AOD) as a proxy for ground-level PM2.5. Particular interest lies in the information on spatial heterogeneity potentially provided by AOD, with important application to estimating and monitoring pollution exposure for public health purposes. Given the temporal and spatio-temporal correlations reported between AOD and PM2.5 , it is tempting to interpret the spatial patterns in AOD as reflecting patterns in PM2.5 . Here we find only limited spatial associations of AOD from three satellite retrievals with PM2.5 over the eastern U.S. at the daily and yearly levels in 2004. We then use statistical modeling to show that the patterns in monthly average AOD poorly reflect patterns in PM2.5 because of systematic, spatially-correlated error in AOD as a proxy for PM2.5 . Furthermore, when we include AOD as a predictor of monthly PM2.5 in a statistical prediction model, AOD provides little additional information to improve predictions of PM2.5 when included in a model that already accounts for land use, emission sources, meteorology and regional variability. These results suggest caution in using spatial variation in AOD to stand in for spatial variation in ground-level PM2.5 in epidemiological analyses and indicate that when PM2.5 monitoring is available, careful statistical modeling outperforms the use of AOD.

INVESTIGATING THE IMPACTS OF STREAMBED HABITAT HETEROGENEITY ON ECOSYSTEM STRUCTURE AND PROCESSES USING BASIC AND APPLIED PERSPECTIVES

Stream restoration often focuses on increasing habitat heterogeneity to reverse ecosystem degradation. However, the connection between heterogeneity and ecosystem structure and processes is poorly understood. We looked to investigate this interaction from both applied and basic science perspectives. For the applied study, we examined two culvert replacements designed to mimic natural stream channels, to see if they were better at maintaining ecosystem processes within as well as upstream and downstream of culverts compared to non-replaced culverts. We measured three ecosystem processes (nutrient uptake, hydrologic characteristics, and coarse particulate organic matter retention) and found that stream simulation culvert restoration improved organic matter retention within culverts, and that there were no differences in processes measured upstream and downstream of both restoration designs. Our results suggest that measurements of ecosystem processes are more likely to show a response to restoration if they match the scale of the restoration activity. For the basic science study, we quantified the longitudinal spatial heterogeneity of physical and biofilm characteristics at microhabitat to segment scales on streams with different streambed variability. We found that all physical characteristics and biofilm characteristics were spatially independent at the macro-habitat scale and greater. Together, these studies demonstrate the importance of scale in ecological interactions and the value of incorporating considerations of scale into restoration activities.

Spatial variability of recent sedimentation in Lake Ohrid (Albania/Macedonia)

Abstract. Lake Ohrid is likely of Pliocene age and thus commonly referred to as the oldest existing lake in Europe. In this study spatial variability of recent sediment composition is assessed using >50 basin wide distributed surface sediment samples. Analysis of biogeochemical bulk parameters, selected metals, pigment concentrations as well as grain size distributions revealed a significant spatial heterogeneity in surface sediment composition. It implies that sedimentation in Lake Ohrid is controlled by an interaction of multiple natural and anthropogenic factors and processes. Major factors controlling surface sediment composition are related to differences in geological catchment characteristics, anthropogenic land use, and a counterclockwise rotating surface water current. In some instances processes controlling sediment composition also seem to impact distribution patterns of biodiversity, which suggests a common interaction of processes responsible for both patterns.

Temporal and small-scale spatial variation in grassland productivity, biomass quality, and nutrient limitation

Characterization of spatial and temporal variation in grassland productivity and nutrition is crucial for a comprehensive understanding of ecosystem function. Although within-site heterogeneity in soil and plant properties has been shown to be relevant for plant community stability, spatiotemporal variability in these factors is still understudied in temperate grasslands. Our study aimed to detect if soil characteristics and plant diversity could explain observed small-scale spatial and temporal variability in grassland productivity, biomass nutrient concentrations, and nutrient limitation. Therefore, we sampled 360 plots of 20 cm × 20 cm each at six consecutive dates in an unfertilized grassland in Southern Germany. Nutrient limitation was estimated using nutrient ratios in plant biomass. Absolute values of, and spatial variability in, productivity, biomass nutrient concentrations, and nutrient limitation were strongly associated with sampling date. In April, spatial heterogeneity was high and most plots showed phosphorous deficiency, while later in the season nitrogen was the major limiting nutrient. Additionally, a small significant positive association between plant diversity and biomass phosphorus concentrations was observed, but should be tested in more detail. We discuss how low biological activity e.g., of soil microbial organisms might have influenced observed heterogeneity of plant nutrition in early spring in combination with reduced active acquisition of soil resources by plants. These early-season conditions are particularly relevant for future studies as they differ substantially from more thoroughly studied later season conditions. Our study underlines the importance of considering small spatial scales and temporal variability to better elucidate mechanisms of ecosystem functioning and plant community assembly.

Spatial and temporal precipitation patterns of the Ebro River Basin, Spain

The Ebro River Basin, with around 85 000 km2 and located in NE Spain, is characterized by the high spatial heterogeneity of its geology, topography, climatology and land use. Rainfall is one of the most important climatic variables studied owing to its non-homogenous behaviour in event and intensity, which creates drought, water runoff and soil erosion with negative environmental and social consequences. In this work we characterized the rainfall variability pattern in the Ebro River Basin using universal multifractal (UM) analysis, which estimates the concentration of the data around the precipitation average (C1, codimension average), the degree of multiscaling behaviour in time (? index) and the maximum probable singularity in the rainfall distribution ( s). A spatial and temporal analysis of the UM parameters is applied to study the possible changes. With this porpoise, 60 daily rainfall series were selected from 132 synthetic series generated by Luna and Balairón (AEMet). These daily rainfall series present a length of 60 years, from 1950 to 2009. Each one of them was subdivided (1950?1970 and 1980?2009) to analyse the difference between the two periods. The range of variation of precipitation amounts and the frequency of dry events between both periods are discussed, as well as the evolution of the UM parameters through the years.

Heterogeneity of atmosfpheric ammonia at the landscape scale and consequences for environmental impact assessment.

We examined the consequences of the spatial heterogeneity of atmospheric ammonia (NH3) by measuring and modelling NH3 concentrations and deposition at 25 m grid resolution for a rural landscape containing intensive poultry farming, agricultural grassland, woodland and moorland. The emission pattern gave rise to a high spatial variability of modelled mean annual NH3 concentrations and dry deposition. Largest impacts were predicted for woodland patches located within the agricultural area, while larger moorland areas were at low risk, due to atmospheric dispersion, prevailing wind direction and low NH3 background. These high resolution spatial details are lost in national scale estimates at 1 km resolution due to less detailed emission input maps. The results demonstrate how the spatial arrangement of sources and sinks is critical to defining the NH3 risk to semi-natural ecosystems. These spatial relationships provide the foundation for local spatial planning approaches to reduce environmental impacts of atmospheric NH3.

Ross River virus disease clusters and spatial relationship with mosquito biting exposure in Redland Shire, southern Queensland, Australia

The spatial heterogeneity in the risk of Ross River virus (family Togaviridae, genus Alphavirus, RRV) disease, the most common mosquito-borne disease in Australia, was examined in Redland Shire in southern Queensland, Australia. Disease cases, complaints from residents of intense mosquito biting exposure, and human population data were mapped using a geographic information system. Surface maps of RRV disease age-sex standardized morbidity ratios and mosquito biting complaint morbidity ratios were created. To determine whether there was significant spatial variation in disease and complaint patterns, a spatial scan analysis method was used to test whether the number of cases and complaints was distributed according to underlying population at risk. Several noncontiguous areas in proximity to productive saline water habitats of Aedes vigilax (Skuse), a recognized vector of RRV, had higher than expected numbers of RRV disease cases and complaints. Disease rates in human populations in areas which had high numbers of adult Ae. vigilax in carbon dioxide- and octenol-baited light traps were up to 2.9 times those in areas that rarely had high numbers of mosquitoes. It was estimated that targeted control of adult Ae. vigilax in these high-risk areas could potentially reduce the RRV disease incidence by an average of 13.6%. Spatial correlation was found between RRV disease risk and complaints from residents of mosquito biting. Based on historical patterns of RRV transmission throughout Redland Shire and estimated future human population growth in areas with higher than average RRV disease incidence, it was estimated that RRV incidence rates will increase by 8% between 2001 and 2021. The use of arbitrary administrative areas that ranged in size from 4.6 to 318.3 km2, has the potential to mask any small scale heterogeneity in disease patterns. With the availability of georeferenced data sets and high-resolution imagery, it is becoming more feasible to undertake spatial analyses at relatively small scales.

Do predators aggregate in response to pest density in agroecosystems? Assessing within-field spatial patterns

1. The spatial heterogeneity of predator populations is an important component of ecological theories pertaining to predator-prey dynamics. Most studies within agricultural fields show spatial correlation (positive or negative) between mean predator numbers and prey abundance across a whole field over time but generally ignore the within-field spatial dimension. We used explicit spatial mapping to determine if generalist predators aggregated within a soybean field, the size of these aggregations and if predator aggregation was associated with pest aggregation, plant damage and predation rate. 2. The study was conducted at Gatton in the Lockyer Valley, 90 km west of Brisbane, Australia. Intensive sampling grids were used to investigate within-field spatial patterns. The first row of each grid was located in a lucerne field (10 m from interface) and the remaining rows were in an adjacent soybean field. At each point on the grid the abundance of foliage-dwelling and ground-dwelling pests and predators was measured, predation rates [using sentinel Helicoverpa armigera (Hubner) egg cards] and plant damage were estimated. Eight grids were sampled across two summer cropping seasons (2000/01, 2001/02). 3. Predators exhibited strong spatial patterning with regions of high and low abundance and activity within what are considered to be uniform soybean fields. Ground-dwelling and foliage-dwelling predators were often aggregated in patches approximately 40 m across. 4. Lycosidae (wolf spiders) displayed aggregation and were consistently more abundant within the lucerne, with a decreasing trap catch with distance from the lucrene/soybean interface. This trend was consistent between subsequent grids in a single field and between fields. 5. The large amount of spatial variability in within-field arthropod abundance (pests and predators) and activity (egg predation and plant damage) indicates that whole field averages were misleading. This result has serious implications for sampling of arthropod abundance and pest management decision-making based on scouting data. 6. There was a great deal of temporal change in the significant spatial patterns observed within a field at each sampling time point during a single season. Predator and pest aggregations observed in these fields were generally not stable for the entire season. 7. Predator aggregation did not correlate consistently with pest aggregation, plant damage or predation rate. Spatial patterns in predator abundance were not associated consistently with any single parameter measured. The most consistent positive association was between foliage-dwelling predators and pests (significant in four of seven grids). Inferring associations between predators and prey based on an intensive one-off sampling grid is difficult, due to the temporal variability in the abundance of each group. 8. Synthesis and applications. This study demonstrated that generalist predator populations are rarely distributed randomly and field edges and adjacent crops can have an influence on within-field predator abundance. This must be considered when estimating arthropod (pest and predator) abundance from a set of samples taken at random locations within a field.

Correcting mean-field approximations for birth-death-movement processes

On the microscale, migration, proliferation and death are crucial in the development, homeostasis and repair of an organism; on the macroscale, such effects are important in the sustainability of a population in its environment. Dependent on the relative rates of migration, proliferation and death, spatial heterogeneity may arise within an initially uniform field; this leads to the formation of spatial correlations and can have a negative impact upon population growth. Usually, such effects are neglected in modeling studies and simple phenomenological descriptions, such as the logistic model, are used to model population growth. In this work we outline some methods for analyzing exclusion processes which include agent proliferation, death and motility in two and three spatial dimensions with spatially homogeneous initial conditions. The mean-field description for these types of processes is of logistic form; we show that, under certain parameter conditions, such systems may display large deviations from the mean field, and suggest computationally tractable methods to correct the logistic-type description.

Analysis of productive performance of crop and animal production systems : an integrated analytical framework

This article presents a two-stage analytical framework that integrates ecological crop (animal) growth and economic frontier production models to analyse the productive efficiency of crop (animal) production systems. The ecological crop (animal) growth model estimates "potential" output levels given the genetic characteristics of crops (animals) and the physical conditions of locations where the crops (animals) are grown (reared). The economic frontier production model estimates "best practice" production levels, taking into account economic, institutional and social factors that cause farm and spatial heterogeneity. In the first stage, both ecological crop growth and economic frontier production models are estimated to calculate three measures of productive efficiency: (1) technical efficiency, as the ratio of actual to "best practice" output levels; (2) agronomic efficiency, as the ratio of actual to "potential" output levels; and (3) agro-economic efficiency, as the ratio of "best practice" to "potential" output levels. Also in the first stage, the economic frontier production model identifies factors that determine technical efficiency. In the second stage, agro-economic efficiency is analysed econometrically in relation to economic, institutional and social factors that cause farm and spatial heterogeneity. The proposed framework has several important advantages in comparison with existing proposals. Firstly, it allows the systematic incorporation of all physical, economic, institutional and social factors that cause farm and spatial heterogeneity in analysing the productive performance of crop and animal production systems. Secondly, the location-specific physical factors are not modelled symmetrically as other economic inputs of production. Thirdly, climate change and technological advancements in crop and animal sciences can be modelled in a "forward-looking" manner. Fourthly, knowledge in agronomy and data from experimental studies can be utilised for socio-economic policy analysis. The proposed framework can be easily applied in empirical studies due to the current availability of ecological crop (animal) growth models, farm or secondary data, and econometric software packages. The article highlights several directions of empirical studies that researchers may pursue in the future.

Multifractal characterisation and analysis of complex networks

Complex networks have been studied extensively due to their relevance to many real-world systems such as the world-wide web, the internet, biological and social systems. During the past two decades, studies of such networks in different fields have produced many significant results concerning their structures, topological properties, and dynamics. Three well-known properties of complex networks are scale-free degree distribution, small-world effect and self-similarity. The search for additional meaningful properties and the relationships among these properties is an active area of current research. This thesis investigates a newer aspect of complex networks, namely their multifractality, which is an extension of the concept of selfsimilarity. The first part of the thesis aims to confirm that the study of properties of complex networks can be expanded to a wider field including more complex weighted networks. Those real networks that have been shown to possess the self-similarity property in the existing literature are all unweighted networks. We use the proteinprotein interaction (PPI) networks as a key example to show that their weighted networks inherit the self-similarity from the original unweighted networks. Firstly, we confirm that the random sequential box-covering algorithm is an effective tool to compute the fractal dimension of complex networks. This is demonstrated on the Homo sapiens and E. coli PPI networks as well as their skeletons. Our results verify that the fractal dimension of the skeleton is smaller than that of the original network due to the shortest distance between nodes is larger in the skeleton, hence for a fixed box-size more boxes will be needed to cover the skeleton. Then we adopt the iterative scoring method to generate weighted PPI networks of five species, namely Homo sapiens, E. coli, yeast, C. elegans and Arabidopsis Thaliana. By using the random sequential box-covering algorithm, we calculate the fractal dimensions for both the original unweighted PPI networks and the generated weighted networks. The results show that self-similarity is still present in generated weighted PPI networks. This implication will be useful for our treatment of the networks in the third part of the thesis. The second part of the thesis aims to explore the multifractal behavior of different complex networks. Fractals such as the Cantor set, the Koch curve and the Sierspinski gasket are homogeneous since these fractals consist of a geometrical figure which repeats on an ever-reduced scale. Fractal analysis is a useful method for their study. However, real-world fractals are not homogeneous; there is rarely an identical motif repeated on all scales. Their singularity may vary on different subsets; implying that these objects are multifractal. Multifractal analysis is a useful way to systematically characterize the spatial heterogeneity of both theoretical and experimental fractal patterns. However, the tools for multifractal analysis of objects in Euclidean space are not suitable for complex networks. In this thesis, we propose a new box covering algorithm for multifractal analysis of complex networks. This algorithm is demonstrated in the computation of the generalized fractal dimensions of some theoretical networks, namely scale-free networks, small-world networks, random networks, and a kind of real networks, namely PPI networks of different species. Our main finding is the existence of multifractality in scale-free networks and PPI networks, while the multifractal behaviour is not confirmed for small-world networks and random networks. As another application, we generate gene interactions networks for patients and healthy people using the correlation coefficients between microarrays of different genes. Our results confirm the existence of multifractality in gene interactions networks. This multifractal analysis then provides a potentially useful tool for gene clustering and identification. The third part of the thesis aims to investigate the topological properties of networks constructed from time series. Characterizing complicated dynamics from time series is a fundamental problem of continuing interest in a wide variety of fields. Recent works indicate that complex network theory can be a powerful tool to analyse time series. Many existing methods for transforming time series into complex networks share a common feature: they define the connectivity of a complex network by the mutual proximity of different parts (e.g., individual states, state vectors, or cycles) of a single trajectory. In this thesis, we propose a new method to construct networks of time series: we define nodes by vectors of a certain length in the time series, and weight of edges between any two nodes by the Euclidean distance between the corresponding two vectors. We apply this method to build networks for fractional Brownian motions, whose long-range dependence is characterised by their Hurst exponent. We verify the validity of this method by showing that time series with stronger correlation, hence larger Hurst exponent, tend to have smaller fractal dimension, hence smoother sample paths. We then construct networks via the technique of horizontal visibility graph (HVG), which has been widely used recently. We confirm a known linear relationship between the Hurst exponent of fractional Brownian motion and the fractal dimension of the corresponding HVG network. In the first application, we apply our newly developed box-covering algorithm to calculate the generalized fractal dimensions of the HVG networks of fractional Brownian motions as well as those for binomial cascades and five bacterial genomes. The results confirm the monoscaling of fractional Brownian motion and the multifractality of the rest. As an additional application, we discuss the resilience of networks constructed from time series via two different approaches: visibility graph and horizontal visibility graph. Our finding is that the degree distribution of VG networks of fractional Brownian motions is scale-free (i.e., having a power law) meaning that one needs to destroy a large percentage of nodes before the network collapses into isolated parts; while for HVG networks of fractional Brownian motions, the degree distribution has exponential tails, implying that HVG networks would not survive the same kind of attack.

Analysis of productive performance of crop production systems: an integrated analytical framework

Recent literature has argued that environmental efficiency (EE), which is built on the materials balance (MB) principle, is more suitable than other EE measures in situations where the law of mass conversation regulates production processes. In addition, the MB-based EE method is particularly useful in analysing possible trade-offs between cost and environmental performance. Identifying determinants of MB-based EE can provide useful information to decision makers but there are very few empirical investigations into this issue. This article proposes the use of data envelopment analysis and stochastic frontier analysis techniques to analyse variation in MB-based EE. Specifically, the article develops a stochastic nutrient frontier and nutrient inefficiency model to analyse determinants of MB-based EE. The empirical study applies both techniques to investigate MB-based EE of 96 rice farms in South Korea. The size of land, fertiliser consumption intensity, cost allocative efficiency, and the share of owned land out of total land are found to be correlated with MB-based EE. The results confirm the presence of a trade-off between MB-based EE and cost allocative efficiency and this finding, favouring policy interventions to help farms simultaneously achieve cost efficiency and MP-based EE.

A systematic review of mathematical models of mosquito-borne pathogen transmission: 1970-2010

Mathematical models of mosquito-borne pathogen transmission originated in the early twentieth century to provide insights into how to most effectively combat malaria. The foundations of the Ross–Macdonald theory were established by 1970. Since then, there has been a growing interest in reducing the public health burden of mosquito-borne pathogens and an expanding use of models to guide their control. To assess how theory has changed to confront evolving public health challenges, we compiled a bibliography of 325 publications from 1970 through 2010 that included at least one mathematical model of mosquito-borne pathogen transmission and then used a 79-part questionnaire to classify each of 388 associated models according to its biological assumptions. As a composite measure to interpret the multidimensional results of our survey, we assigned a numerical value to each model that measured its similarity to 15 core assumptions of the Ross–Macdonald model. Although the analysis illustrated a growing acknowledgement of geographical, ecological and epidemiological complexities in modelling transmission, most models during the past 40 years closely resemble the Ross–Macdonald model. Modern theory would benefit from an expansion around the concepts of heterogeneous mosquito biting, poorly mixed mosquito-host encounters, spatial heterogeneity and temporal variation in the transmission process.

A finite volume scheme with preconditioned Lanczos method for two-dimensional space-fractional reaction–diffusion equations

Fractional differential equations have been increasingly used as a powerful tool to model the non-locality and spatial heterogeneity inherent in many real-world problems. However, a constant challenge faced by researchers in this area is the high computational expense of obtaining numerical solutions of these fractional models, owing to the non-local nature of fractional derivatives. In this paper, we introduce a finite volume scheme with preconditioned Lanczos method as an attractive and high-efficiency approach for solving two-dimensional space-fractional reaction–diffusion equations. The computational heart of this approach is the efficient computation of a matrix-function-vector product f(A)bf(A)b, where A A is the matrix representation of the Laplacian obtained from the finite volume method and is non-symmetric. A key aspect of our proposed approach is that the popular Lanczos method for symmetric matrices is applied to this non-symmetric problem, after a suitable transformation. Furthermore, the convergence of the Lanczos method is greatly improved by incorporating a preconditioner. Our approach is show-cased by solving the fractional Fisher equation including a validation of the solution and an analysis of the behaviour of the model.