Modelling ecological systems in a changing world

Perspectives on computational modeling of biological systems and the significance of the SysMod community

Motivation

In recent years, applying computational modeling to systems biology has caused a substantial surge in both discovery and practical applications and a significant shift in our understanding of the complexity inherent in biological systems.

Results

In this perspective article, we briefly overview computational modeling in biology, highlighting recent advancements such as multi-scale modeling due to the omics revolution, single-cell technology, and integration of artificial intelligence and machine learning approaches. We also discuss the primary challenges faced: integration, standardization, model complexity, scalability, and interdisciplinary collaboration. Lastly, we highlight the contribution made by the Computational Modeling of Biological Systems (SysMod) Community of Special Interest (COSI) associated with the International Society of Computational Biology (ISCB) in driving progress within this rapidly evolving field through community engagement (via both in person and virtual meetings, social media interactions), webinars, and conferences.

Availability and implementation

Additional information about SysMod is available at https://sysmod.info.

Climatic niche evolution and niche conservatism of Nymphaea species in Africa, South America, and Australia

BACKGROUND

Interest in the evolution of climatic niches, particularly in understanding the potential adaptive responses of species under climate change, has increased both theoretically and within macroecological studies. These studies have provided valuable insights into how climatic traits of species influence their niche evolution. In this study, we aim to investigate whether niche conservatism plays a role in the species diversification of Nymphaea, a group of aquatic plants with a cosmopolitan distribution that is facing severe habitat loss. We will use climatic models and phylogenetic data for 23 species to reconstruct Nymphaea's niche evolution, measure niche overlap, and assess disparity through time while testing for evolutionary models.

RESULTS

There was a lot of overlap in niches both within and between groups, especially for species that can be found in many places. The breadth and peaks of the niche profile varied depending on the bioclimatic variables, which suggested that the species evolved differently to cope with changes in climate. The analysis also showed that evolutionary changes happened across the phylogeny, with weak to moderate signals. The morphological disparity index (MDI) values indicated that there were disparities within subclades over time but not between or among them. Niche reconstruction and evolution analysis revealed both convergent and divergent evolution among various variables. For example, N. immutabilis, N. atrans, N. violancea, and N. nouchali evolved towards intermediate temperatures for bio2 and bio3 (isothermity) while moving towards extreme temperatures for bio8 and bio9 (wettest and driest average quarterly temperatures).

CONCLUSION

Our study will improve our understanding of how changes in climatic niches are potentially driving the evolution of Nymphaea. It has significant scientific implications for the limits, assemblages, evolution, and diversification of species. This information is crucial for the ongoing efforts of conservation and management, particularly considering the inevitable effects of climate change.

A density functional theory for ecology across scales

Ecology lacks a holistic approach that can model phenomena across temporal and spatial scales, largely because of the challenges in modelling systems with a large number of interacting constituents. This hampers our understanding of complex ecosystems and the impact that human interventions (e.g., deforestation, wildlife harvesting and climate change) have on them. Here we use density functional theory, a computational method for many-body problems in physics, to develop a computational framework for ecosystem modelling. Our methods accurately fit experimental and synthetic data of interacting multi-species communities across spatial scales and can project to unseen data. As the key concept we establish and validate a cost function that encodes the trade-offs between the various ecosystem components. We show how this single general modelling framework delivers predictions on par with established, but specialised, approaches for systems from predatory microbes to territorial flies to tropical tree communities. Our density functional framework thus provides a promising avenue for advancing our understanding of ecological systems.

Strategic roadmap to assess forest vulnerability under air pollution and climate change

Although it is an integral part of global change, most of the research addressing the effects of climate change on forests have overlooked the role of environmental pollution. Similarly, most studies investigating the effects of air pollutants on forests have generally neglected the impacts of climate change. We review the current knowledge on combined air pollution and climate change effects on global forest ecosystems and identify several key research priorities as a roadmap for the future. Specifically, we recommend (1) the establishment of much denser array of monitoring sites, particularly in the South Hemisphere; (2) further integration of ground and satellite monitoring; (3) generation of flux-based standards and critical levels taking into account the sensitivity of dominant forest tree species; (4) long-term monitoring of N, S, P cycles and base cations deposition together at global scale; (5) intensification of experimental studies, addressing the combined effects of different abiotic factors on forests by assuring a better representation of taxonomic and functional diversity across the ~73,000 tree species on Earth; (6) more experimental focus on phenomics and genomics; (7) improved knowledge on key processes regulating the dynamics of radionuclides in forest systems; and (8) development of models integrating air pollution and climate change data from long-term monitoring programs.

Ozone modelling and mapping for risk assessment: An overview of different approaches for human and ecosystems health

Tropospheric ozone (O₃) is one of the most concernedair pollutants dueto its widespread impacts on land vegetated ecosystems and human health. Ozone is also the third greenhouse gas for radiative forcing. Consequently, it should be carefully and continuously monitored to estimate its potential adverse impacts especially inthose regions where concentrations are high. Continuous large-scale O₃ concentrations measurement is crucial but may be unfeasible because of economic and practical limitations; therefore, quantifying the real impact of O₃over large areas is currently an open challenge. Thus, one of the final objectives of O₃ modelling is to reproduce maps of continuous concentrations (both spatially and temporally) and risk assessment for human and ecosystem health. We here reviewedthe most relevant approaches used for O₃ modelling and mapping starting from the simplest geo-statistical approaches andincreasing in complexity up to simulations embedded into the global/regional circulation models and pro and cons of each mode are highlighted. The analysis showed that a simpler approach (mostly statistical models) is suitable for mappingO₃concentrationsat the local scale, where enough O₃concentration data are available. The associated error in mapping can be reduced by using more complex methodologies, based on co-variables. The models available at the regional or global level are used depending on the needed resolution and the domain where they are applied to. Increasing the resolution corresponds to an increase in the prediction but only up to a certain limit. However, with any approach, the ensemble models should be preferred.

Bridging Theories for Ecosystem Stability Through Structural Sensitivity Analysis of Ecological Models in Equilibrium

Ecologists are challenged by the need to bridge and synthesize different approaches and theories to obtain a coherent understanding of ecosystems in a changing world. Both food web theory and regime shift theory shine light on mechanisms that confer stability to ecosystems, but from different angles. Empirical food web models are developed to analyze how equilibria in real multi-trophic ecosystems are shaped by species interactions, and often include linear functional response terms for simple estimation of interaction strengths from observations. Models of regime shifts focus on qualitative changes of equilibrium points in a slowly changing environment, and typically include non-linear functional response terms. Currently, it is unclear how the stability of an empirical food web model, expressed as the rate of system recovery after a small perturbation, relates to the vulnerability of the ecosystem to collapse. Here, we conduct structural sensitivity analyses of classical consumer-resource models in equilibrium along an environmental gradient. Specifically, we change non-proportional interaction terms into proportional ones, while maintaining the equilibrium biomass densities and material flux rates, to analyze how alternative model formulations shape the stability properties of the equilibria. The results reveal no consistent relationship between the stability of the original models and the proportionalized versions, even though they describe the same biomass values and material flows. We use these findings to critically discuss whether stability analysis of observed equilibria by empirical food web models can provide insight into regime shift dynamics, and highlight the challenge of bridging alternative modelling approaches in ecology and beyond.

Invasive potential of tropical fruit flies in temperate regions under climate change

Tropical fruit flies are considered among the most economically important invasive species detected in temperate areas of the United States and the European Union. Detections often trigger quarantine and eradication programs that are conducted without a holistic understanding of the threat posed. Weather-driven physiologically-based demographic models are used to estimate the geographic range, relative abundance, and threat posed by four tropical tephritid fruit flies (Mediterranean fruit fly, melon fly, oriental fruit fly, and Mexican fruit fly) in North and Central America, and the European-Mediterranean region under extant and climate change weather (RCP8.5 and A1B scenarios). Most temperate areas under tropical fruit fly propagule pressure have not been suitable for establishment, but suitability is predicted to increase in some areas with climate change. To meet this ongoing challenge, investments are needed to collect sound biological data to develop mechanistic models to predict the geographic range and relative abundance of these and other invasive species, and to put eradication policies on a scientific basis.

An introduction to agent-based models as an accessible surrogate to field-based research and teaching

There are many barriers to fieldwork including cost, time, and physical ability. Unfortunately, these barriers disproportionately affect minority communities and create a disparity in access to fieldwork in the natural sciences. Travel restrictions, concerns about our carbon footprint, and the global lockdown have extended this barrier to fieldwork across the community and led to increased anxiety about gaps in productivity, especially among graduate students and early-career researchers. In this paper, we discuss agent-based modeling as an open-source, accessible, and inclusive resource to substitute for lost fieldwork during COVID-19 and for future scenarios of travel restrictions such as climate change and economic downturn. We describe the benefits of Agent-Based models as a teaching and training resource for students across education levels. We discuss how and why educators and research scientists can implement them with examples from the literature on how agent-based models can be applied broadly across life science research. We aim to amplify awareness and adoption of this technique to broaden the diversity and size of the agent-based modeling community in ecology and evolutionary research. Finally, we discuss the challenges facing agent-based modeling and discuss how quantitative ecology can work in tandem with traditional field ecology to improve both methods.

A guide to ecosystem models and their environmental applications

Applied ecology has traditionally approached management problems through a simplified, single-species lens. Repeated failures of single-species management have led us to a new paradigm - managing at the ecosystem level. Ecosystem management involves a complex array of interacting organisms, processes and scientific disciplines. Accounting for interactions, feedback loops and dependencies between ecosystem components is therefore fundamental to understanding and managing ecosystems. We provide an overview of the main types of ecosystem models and their uses, and discuss challenges related to modelling complex ecological systems. Existing modelling approaches typically attempt to do one or more of the following: describe and disentangle ecosystem components and interactions; make predictions about future ecosystem states; and inform decision making by comparing alternative strategies and identifying important uncertainties. Modelling ecosystems is challenging, particularly when balancing the desire to represent many components of an ecosystem with the limitations of available data and the modelling objective. Explicitly considering different forms of uncertainty is therefore a primary concern. We provide some recommended strategies (such as ensemble ecosystem models and multi-model approaches) to aid the explicit consideration of uncertainty while also meeting the challenges of modelling ecosystems.

Measurement Uncertainty in Ecological and Environmental Models

In many applied cases of ecological and environmental modeling there is sizeable variation in the independent variables as a result of measurement and sampling errors. This uncertainty may lead to biased predictions. It is possible to avoid this problem by increased sampling and by modeling the errors using hierarchical modeling.

Advancing an interdisciplinary framework to study seed dispersal ecology

Although dispersal is generally viewed as a crucial determinant for the fitness of any organism, our understanding of its role in the persistence and spread of plant populations remains incomplete. Generalizing and predicting dispersal processes are challenging due to context dependence of seed dispersal, environmental heterogeneity and interdependent processes occurring over multiple spatial and temporal scales. Current population models often use simple phenomenological descriptions of dispersal processes, limiting their ability to examine the role of population persistence and spread, especially under global change. To move seed dispersal ecology forward, we need to evaluate the impact of any single seed dispersal event within the full spatial and temporal context of a plant's life history and environmental variability that ultimately influences a population's ability to persist and spread. In this perspective, we provide guidance on integrating empirical and theoretical approaches that account for the context dependency of seed dispersal to improve our ability to generalize and predict the consequences of dispersal, and its anthropogenic alteration, across systems. We synthesize suitable theoretical frameworks for this work and discuss concepts, approaches and available data from diverse subdisciplines to help operationalize concepts, highlight recent breakthroughs across research areas and discuss ongoing challenges and open questions. We address knowledge gaps in the movement ecology of seeds and the integration of dispersal and demography that could benefit from such a synthesis. With an interdisciplinary perspective, we will be able to better understand how global change will impact seed dispersal processes, and potential cascading effects on plant population persistence, spread and biodiversity.

Predicting future climate at high spatial and temporal resolution

Most studies on the biological effects of future climatic changes rely on seasonally aggregated, coarse-resolution data. Such data mask spatial and temporal variability in microclimate driven by terrain, wind and vegetation, and ultimately bear little resemblance to the conditions that organisms experience in the wild. Here, I present the methods for providing fine-grained, hourly and daily estimates of current and future temperature and soil moisture over decadal timescales. Observed climate data and spatially coherent probabilistic projections of daily future weather were disaggregated to hourly and used to drive empirically calibrated physical models of thermal and hydrological microclimates. Mesoclimatic effects (cold-air drainage, coastal exposure and elevation) were determined from the coarse-resolution climate surfaces using thin-plate spline models with coastal exposure and elevation as predictors. Differences between micro and mesoclimate temperatures were determined from terrain, vegetation and ground properties using energy balance equations. Soil moisture was computed in a thin upper layer and an underlying deeper layer, and the exchange of water between these layers was calculated using the van Genuchten equation. Code for processing the data and running the models is provided as a series of R packages. The methods were applied to the Lizard Peninsula, United Kingdom, to provide hourly estimates of temperature (100 m grid resolution over entire area, 1 m for a selected area) for the periods 1983-2017 and 2041-2049. Results indicated that there is a fine-resolution variability in climatic changes, driven primarily by interactions between landscape features and decadal trends in weather conditions. High-temporal resolution extremes in conditions under future climate change were predicted to be considerably less novel than the extremes estimated using seasonally aggregated variables. The study highlights the need to more accurately estimate the future climatic conditions experienced by organisms and equips biologists with the means to do so.

Predicting population responses to environmental change from individual-level mechanisms: towards a standardized mechanistic approach

Animal populations will mediate the response of global biodiversity to environmental changes. Population models are thus important tools for both understanding and predicting animal responses to uncertain future conditions. Most approaches, however, are correlative and ignore the individual-level mechanisms that give rise to population dynamics. Here, we assess several existing population modelling approaches and find limitations to both 'correlative' and 'mechanistic' models. We advocate the need for a standardized mechanistic approach for linking individual mechanisms (physiology, behaviour, and evolution) to population dynamics in spatially explicit landscapes. Such an approach is potentially more flexible and informative than current population models. Key to realizing this goal, however, is overcoming current data limitations, the development and testing of eco-evolutionary theory to represent interactions between individual mechanisms, and standardized multi-dimensional environmental change scenarios which incorporate multiple stressors. Such progress is essential in supporting environmental decisions in uncertain future conditions.

Temporal development of Drosophila embryos is highly robust across a wide temperature range

Development is a process precisely coordinated in both space and time. Spatial precision has been quantified in a number of developmental systems, and such data have contributed significantly to our understanding of, for example, morphogen gradient interpretation. However, comparatively little quantitative analysis has been performed on timing and temporal coordination during development. Here, we use Drosophila to explore the temporal robustness of embryonic development within physiologically normal temperatures. We find that development is temporally very precise across a wide range of temperatures in the three Drosophila species investigated. However, we find temperature dependence in the timing of developmental events. A simple model incorporating history dependence can explain the developmental temporal trajectories. Interestingly, history dependence is temperature-specific, with either effective negative or positive feedback at different temperatures. We also find that embryos are surprisingly robust to shifting temperatures during embryogenesis. We further identify differences between tropical and temperate species, potentially due to different mechanisms regulating temporal development that depend on the local environment. Our data show that Drosophila embryonic development is temporally robust across a wide range of temperatures. This robustness shows interesting species-specific differences that are suggestive of different sensitivity to temperature fluctuations between Drosophila species.

Bio-Based Production Systems: Why Environmental Assessment Needs to Include Supporting Systems

The transition to a bio-based economy is expected to deliver substantial environmental and economic benefits. However, bio-based production systems still come with significant environmental challenges, and there is a need for assessment methods that are adapted for the specific characteristics of these systems. In this review, we investigated how the environmental aspects of bio-based production systems differ from those of non-renewable systems, what requirements these differences impose when assessing their sustainability, and to what extent mainstream assessment methods fulfil these requirements. One unique characteristic of bio-based production is the need to maintain the regenerative capacity of the system. The necessary conditions for maintaining regenerative capacity are often provided through direct or indirect interactions between the production system and surrounding “supporting” systems. Thus, in the environmental assessment, impact categories affected in both the primary production system and the supporting systems need to be included, and impact models tailored to the specific context of the study should be used. Development in this direction requires efforts to broaden the system boundaries of conventional environmental assessments, to increase the level of spatial and temporal differentiation, and to improve our understanding of how local uniqueness and temporal dynamics affect the performance of the investigated system.

Anticipating global terrestrial ecosystem state change using FLUXNET

Ecosystems can be characterized as complex systems that traverse a variety of functional and structural states in response to changing bioclimatic forcings. A central challenge of global change biology is the robust empirical description of these states and state transitions. An ecosystem's functional state can be empirically described using Process Networks (PN) that use timeseries observations to determine the strength of process-level functional couplings between ecosystem components. A globally extensive source of in-situ observations of terrestrial ecosystem dynamics is the FLUXNET eddy-covariance network that provides standardized observations of micrometeorology and carbon, water, and energy flux dynamics. We employ the LaThuile FLUXNET synthesis dataset to delineate each month's functional state for 204 sites, yielding the LaThuile PN version 1.0 database that describes the strength of an ecosystem's functional couplings from air temperature and precipitation to carbon fluxes during each site-month. Then we calculate the elasticity of these couplings to seasonal scale forcings: air temperature, precipitation, solar radiation, and phenophase. Finally, we train artificial neural networks to extrapolate these elasticities from 204 sites to the globe, yielding maps of the estimated functional elasticity of every terrestrial ecosystem's functional states to changing seasonal bioclimatic forcings. These maps provide theoretically novel resource that can be used to anticipate ecological state transitions in response to climate change and to validate process-based models of ecological change. These elasticity maps show that each ecosystem can be expected to respond uniquely to changing forcings. Tropical forests, hot deserts, savannas, and high elevations are most elastic to climate change, and elasticity of ecosystems to seasonal air temperature is on average an order of magnitude higher than elasticity to other bioclimatic forcings. We also observed a reasonable amount of moderate relationships between functional elasticity and structural state change across different ecosystems.

Simulation tools for assessment of tick suppression treatments of Rhipicephalus (Boophilus) microplus on non-lactating dairy cattle in Puerto Rico

BACKGROUND

The southern cattle fever tick (SCFT), Rhipicephalus (Boophilus) microplus, remains endemic in Puerto Rico. Systematic treatment programmes greatly reduced and even eradicated temporarily this tick from the island. However, a systemic treatment programme that includes integrated management practices for livestock against SCFT remains to be established in the island. We describe a spatially-explicit, individual-based model that simulates climate-livestock-SCFT-landscape interactions. This model was developed as an investigative tool to aid in a research project on integrated management of the SCFT that took place in Puerto Rico between 2014 and 2017. We used the model to assess the efficacy of tick suppression and probability of tick elimination when applying safer acaricides at 3-week intervals to different proportions of a herd of non-lactating dairy cattle.

RESULTS

Probabilities of eliminating host-seeking larvae from the simulated system decreased from ≈ 1 to ≈ 0 as the percentage of cattle treated decreased from 65 to 45, with elimination probabilities ≈ 1 at higher treatment percentages and ≈ 0 at lower treatment percentages. For treatment percentages between 65% and 45%, a more rapid decline in elimination probabilities was predicted by the version of the model that produced higher densities of host-seeking larvae. Number of weeks after the first acaricide application to elimination of host-seeking larvae was variable among replicate simulations within treatment percentages, with within-treatment variation increasing markedly at treatment percentages ≤ 65. Number of weeks after first application to elimination generally varied between 30 and 40 weeks for those treatment percentages with elimination probabilities ≈ 1.

CONCLUSIONS

Explicit simulation of the spatial and temporal dynamics of off-host (host-seeking) larvae in response to control methods should be an essential element of research that involves the evaluation of integrated SCFT management programmes. This approach could provide the basis to evaluate novel control technologies and to develop protocols for their cost-effective use with other treatment methods.

Ringed seal demography in a changing climate

Climate change is affecting species' distributions and abundances worldwide. Baseline population estimates, against which future observations may be compared, are necessary if we are to detect ecological change. Arctic sea ice ecosystems are changing rapidly and we lack baseline population estimates for many ice-associated species. Provided we can detect them, changes in Arctic marine ecosystems may be signaled by changes in indicator species such as ringed seals (Pusa hispida). Ringed seal monitoring has provided estimates of survival and fertility rates, but these have not been used for population-level inference. Using matrix population models, we synthesized existing demographic parameters to obtain estimates of historical ringed seal population growth and structure in Amundsen Gulf and Prince Albert Sound, Canada. We then formalized existing hypotheses about the effects of emerging environmental stressors (i.e., earlier spring ice breakup and reduced snow depth) on ringed seal pup survival. Coupling the demographic model to ice and snow forecasts available from the Coupled Model Intercomparison Project resulted in projections of ringed seal population size and structure up to the year 2100. These projections showed median declines in population size ranging from 50% to 99%. Corresponding to these projected declines were substantial changes in population structure, with increasing proportions of ringed seal pups and adults and declining proportions of juveniles. We explored if currently collected, harvest-based data could be used to detect the projected changes in population stage structure. Our model suggests that at a present sample size of 100 seals per year, the projected changes in stage structure would only be reliably detected by mid-century, even for the most extreme climate models. This modeling process revealed inconsistencies in existing estimates of ringed seal demographic rates. Mathematical population models such as these can contribute both to understanding past population trends as well as predicting future ones, both of which are necessary if we are to detect and interpret future observations.

Harnessing Deep Learning in Ecology: An Example Predicting Bark Beetle Outbreaks

Addressing current global challenges such as biodiversity loss, global change, and increasing demands for ecosystem services requires improved ecological prediction. Recent increases in data availability, process understanding, and computing power are fostering quantitative approaches in ecology. However, flexible methodological frameworks are needed to utilize these developments towards improved ecological prediction. Deep learning is a rapidly evolving branch of machine learning, yet has received only little attention in ecology to date. It refers to the training of deep neural networks (DNNs), i.e. artificial neural networks consisting of many layers and a large number of neurons. We here provide a reproducible example (including code and data) of designing, training, and applying DNNs for ecological prediction. Using bark beetle outbreaks in conifer-dominated forests as an example, we show that DNNs are well able to predict both short-term infestation risk at the local scale and long-term outbreak dynamics at the landscape level. We furthermore highlight that DNNs have better overall performance than more conventional approaches to predicting bark beetle outbreak dynamics. We conclude that DNNs have high potential to form the backbone of a comprehensive disturbance forecasting system. More broadly, we argue for an increased utilization of the predictive power of DNNs for a wide range of ecological problems.

Ecosystem heterogeneity and diversity mitigate Amazon forest resilience to frequent extreme droughts

The impact of increases in drought frequency on the Amazon forest's composition, structure and functioning remain uncertain. We used a process- and individual-based ecosystem model (ED2) to quantify the forest's vulnerability to increased drought recurrence. We generated meteorologically realistic, drier-than-observed rainfall scenarios for two Amazon forest sites, Paracou (wetter) and Tapajós (drier), to evaluate the impacts of more frequent droughts on forest biomass, structure and composition. The wet site was insensitive to the tested scenarios, whereas at the dry site biomass declined when average rainfall reduction exceeded 15%, due to high mortality of large-sized evergreen trees. Biomass losses persisted when year-long drought recurrence was shorter than 2-7 yr, depending upon soil texture and leaf phenology. From the site-level scenario results, we developed regionally applicable metrics to quantify the Amazon forest's climatological proximity to rainfall regimes likely to cause biomass loss > 20% in 50 yr according to ED2 predictions. Nearly 25% (1.8 million km² ) of the Amazon forests could experience frequent droughts and biomass loss if mean annual rainfall or interannual variability changed by 2σ. At least 10% of the high-emission climate projections (CMIP5/RCP8.5 models) predict critically dry regimes over 25% of the Amazon forest area by 2100.

Vegetation demographics in Earth System Models: A review of progress and priorities

Numerous current efforts seek to improve the representation of ecosystem ecology and vegetation demographic processes within Earth System Models (ESMs). These developments are widely viewed as an important step in developing greater realism in predictions of future ecosystem states and fluxes. Increased realism, however, leads to increased model complexity, with new features raising a suite of ecological questions that require empirical constraints. Here, we review the developments that permit the representation of plant demographics in ESMs, and identify issues raised by these developments that highlight important gaps in ecological understanding. These issues inevitably translate into uncertainty in model projections but also allow models to be applied to new processes and questions concerning the dynamics of real-world ecosystems. We argue that stronger and more innovative connections to data, across the range of scales considered, are required to address these gaps in understanding. The development of first-generation land surface models as a unifying framework for ecophysiological understanding stimulated much research into plant physiological traits and gas exchange. Constraining predictions at ecologically relevant spatial and temporal scales will require a similar investment of effort and intensified inter-disciplinary communication.

Making ecological models adequate

Critical evaluation of the adequacy of ecological models is urgently needed to enhance their utility in developing theory and enabling environmental managers and policymakers to make informed decisions. Poorly supported management can have detrimental, costly or irreversible impacts on the environment and society. Here, we examine common issues in ecological modelling and suggest criteria for improving modelling frameworks. An appropriate level of process description is crucial to constructing the best possible model, given the available data and understanding of ecological structures. Model details unsupported by data typically lead to over parameterisation and poor model performance. Conversely, a lack of mechanistic details may limit a model's ability to predict ecological systems' responses to management. Ecological studies that employ models should follow a set of model adequacy assessment protocols that include: asking a series of critical questions regarding state and control variable selection, the determinacy of data, and the sensitivity and validity of analyses. We also need to improve model elaboration, refinement and coarse graining procedures to better understand the relevancy and adequacy of our models and the role they play in advancing theory, improving hind and forecasting, and enabling problem solving and management.

Altered development, oxidative stress and DNA damage in Leptodactylus chaquensis (Anura: Leptodactylidae) larvae exposed to poultry litter

Poultry litter (PL), which is usually used as organic fertilizer, is a source of nutrients, metals, veterinary pharmaceuticals and bacterial pathogens, which, through runoff, may end up in the nearest aquatic ecosystems. In this study, Leptodactylus chaquensis at different development stages (eggs, larval stages 28 and 31 here referred to as stages I, II and III respectively) were exposed to PL test sediments as follows: 6.25% (T1), 12.5% (T2); 25% (T3); 50% (T4); 75% (T5); 100% PL (T6) and to dechlorinated water as control. Larval survival, development endpoints (growth rate -GR-, development rate -DR-, abnormalities), antioxidant enzyme activities (Catalase -CAT- and Glutathione-S-Transferase -GST-), and genotoxic effect (DNA damage index by the Comet assay) were analyzed at different times. In stage I, no egg eclosion was observed in treatments T3-T6, and 50% of embryo mortality was recorded after 24h of exposure to T2. In stages II and III, mortality in treatments T3-T6 reached 100% between 24 and 48h. In the three development stages evaluated, the DR and GR were higher in controls than in PL treatments (T1, T2), except for those T1-treated larvae of stage II. Larvae of stage I showed five types of morphological abnormalities, being diamond body shape and lateral displacement of the intestine the most prevalent in T1, whereas larvae of stages II and III presented lower prevalence of abnormalities. In stage I, CAT activity was similar to that of control (p>0.05), whereas it was higher in T1- and T2- treated larvae of stages II and III than controls (p<0.05). In stages I and III, GST activity was similar to that of controls (p>0.05), whereas it was inhibited in T1-treated larvae of stage II (p<0.05). T1- and T2-treated larvae of stages II and III caused higher DNA damage respect to controls (p<0.05), varying from medium to severe damage (comet types II, III and IV). These results showed that PL treatments altered development and growth and induced oxidative stress and DNA damage, resulting ecotoxic for L. chaquensis larvae.

Fine-scale climate change: modelling spatial variation in biologically meaningful rates of warming

The existence of fine-grain climate heterogeneity has prompted suggestions that species may be able to survive future climate change in pockets of suitable microclimate, termed 'microrefugia'. However, evidence for microrefugia is hindered by lack of understanding of how rates of warming vary across a landscape. Here, we present a model that is applied to provide fine-grained, multidecadal estimates of temperature change based on the underlying physical processes that influence microclimate. Weather station and remotely derived environmental data were used to construct physical variables that capture the effects of terrain, sea surface temperatures, altitude and surface albedo on local temperatures, which were then calibrated statistically to derive gridded estimates of temperature. We apply the model to the Lizard Peninsula, United Kingdom, to provide accurate (mean error = 1.21 °C; RMS error = 1.63 °C) hourly estimates of temperature at a resolution of 100 m for the period 1977-2014. We show that rates of warming vary across a landscape primarily due to long-term trends in weather conditions. Total warming varied from 0.87 to 1.16 °C, with the slowest rates of warming evident on north-east-facing slopes. This variation contributed to substantial spatial heterogeneity in trends in bioclimatic variables: for example, the change in the length of the frost-free season varied from +11 to -54 days and the increase in annual growing degree-days from 51 to 267 °C days. Spatial variation in warming was caused primarily by a decrease in daytime cloud cover with a resulting increase in received solar radiation, and secondarily by a decrease in the strength of westerly winds, which has amplified the effects on temperature of solar radiation on west-facing slopes. We emphasize the importance of multidecadal trends in weather conditions in determining spatial variation in rates of warming, suggesting that locations experiencing least warming may not remain consistent under future climate change.

BioVeL: a virtual laboratory for data analysis and modelling in biodiversity science and ecology

Background

Making forecasts about biodiversity and giving support to policy relies increasingly on large collections of data held electronically, and on substantial computational capability and capacity to analyse, model, simulate and predict using such data. However, the physically distributed nature of data resources and of expertise in advanced analytical tools creates many challenges for the modern scientist. Across the wider biological sciences, presenting such capabilities on the Internet (as “Web services”) and using scientific workflow systems to compose them for particular tasks is a practical way to carry out robust “in silico” science. However, use of this approach in biodiversity science and ecology has thus far been quite limited.

Results

BioVeL is a virtual laboratory for data analysis and modelling in biodiversity science and ecology, freely accessible via the Internet. BioVeL includes functions for accessing and analysing data through curated Web services; for performing complex in silico analysis through exposure of R programs, workflows, and batch processing functions; for on-line collaboration through sharing of workflows and workflow runs; for experiment documentation through reproducibility and repeatability; and for computational support via seamless connections to supporting computing infrastructures. We developed and improved more than 60 Web services with significant potential in many different kinds of data analysis and modelling tasks. We composed reusable workflows using these Web services, also incorporating R programs. Deploying these tools into an easy-to-use and accessible ‘virtual laboratory’, free via the Internet, we applied the workflows in several diverse case studies. We opened the virtual laboratory for public use and through a programme of external engagement we actively encouraged scientists and third party application and tool developers to try out the services and contribute to the activity.

Conclusions

Our work shows we can deliver an operational, scalable and flexible Internet-based virtual laboratory to meet new demands for data processing and analysis in biodiversity science and ecology. In particular, we have successfully integrated existing and popular tools and practices from different scientific disciplines to be used in biodiversity and ecological research.

Electronic supplementary material

The online version of this article (doi:10.1186/s12898-016-0103-y) contains supplementary material, which is available to authorized users.