IAP-25-140
Mathematical modelling of Ecosystem Resilience in a Changing Climate
Ongoing climate change threatens to instigate dramatic changes in ecosystems across the globe, prompting a renewed focus on understanding ecosystem stability and resilience. Vegetation, from forests and savannas to semi-arid drylands, is significantly impacted by extreme weather events (e.g. fires, flooding, drought) and shifts in climate on longer scales (e.g. variations in annual rainfall). Moreover, large, vegetated ecosystems are crucial in carbon capture and the global carbon cycle, as well as underpinning the livelihoods of 2.3 billion people. This project will explore a class of mathematical models of vegetative ecosystems to capture these relationships and answer questions about climate-driven impacts on vegetation.
Existing models of these ecosystems can be arranged into two broad categories: (1) Simplified phenomenological models studied by theoreticians that allow deep understanding of dynamics (Meron, 2011), or (2) Detailed process-based models that are realistic, but difficult to analyse or interpret (Stewart et al, 2014), limiting their utility for assessing ecosystem resilience. Our goal is to use mathematical coarse graining techniques to derive a new class of models that bridge the gap between approaches (1) and (2). These models should faithfully reflect underlying ecological processes but remain tractable enough to allow for detailed understanding of their dynamics (Patterson et al., 2020). This tractability will also enable us to apply resilience metrics and methods from mathematical ecology to interpret model outputs in the context of wider climactic shifts. Finally, our models and outputs will be validated by and compared to remote-sensed data. Overall, this research will enable us to better understand the forces underlying vegetation patterning and resource losses, and the signals that such patterning gives us in terms of early-warning signs of ecosystem change.
This project spans mathematical ecology, geomorphological modelling, and statistical analysis of remote-sensing data. The supervisory team brings together essential expertise spanning environmental science, geography, and mathematics. This collaborative approach is vital to capturing the complexity of the ecological processes of interest and translating research findings into practical insights, including tools such as stakeholder-facing dashboards for decision-making and ecosystem monitoring.
Methodology
Depending on the student, differing amounts of time will be spent on each facet of the project, though some level of engagement across the different approaches will be necessary. They will begin by building on the existing modelling frameworks developed in Durham to create tractable models for spatial fire spread (Patterson et al., 2020) and herbivory (Stewart et al., 2014, Tiwari et al. 2025) that can be compared to remote sensing data. Once the models are validated, we will be to quantify the likelihood of state transitions in patterned semi-arid ecosystems, develop early warning signals for these transitions, and further explore how these systems respond to ongoing climate change.
There will be an emphasis throughout on building interpretable mechanistic models, where modelling assumptions, such as time and length scales where the mechanisms are likely to be valid, are clear. We will explicitly consider a spectrum of modelling approaches to balance different approaches with the need for systematic analyses and comparisons with data. We will leverage mathematical techniques from areas like theoretical ecology to understand dynamical motifs in both simple and detailed models, such as pattern formation as a resilience mechanism for vegetation. The use of landscape-flux theory (Xu et al., 2023), as well as other tools such as bifurcation analysis and continuation, will be used to understand the variety and likelihood of ecosystem transitions expected as climate-relevant parameters vary.
We will co-create stakeholder-facing dashboards allow easy communication of management and policy-relevant lessons. These tools will be built using modern Web frameworks, such as WebGL within JavaScript, that could interface with backend servers for model execution or data access. Examples of existing projects of this form include VisualPDE.com (led by Krause) and LimbNET.embl.es.
Project Timeline
Year 1
– Literature review across mechanistic, statistical vegetation modelling,and geomorphology.
– Gain proficiency in data analysis and mathematical modelling, including compartmental process-based frameworks and partial differential equation models (PDEs).
– Learn simulation and analysis of mathematical models (dynamical systems).
– Analyse simplified ecosystem models.
Year 2
– Develop fast codes for realistic process-based models.
– Connect process-based modelling with simplified models via coarse-graining.
– Prepare dashboards and interactive visualisations for communicating findings with public and public-policy audiences.
– Prepare modelling work for publication.
Year 3
– Complete development and testing of the suite of models.
– Compare simplified and mechanistic models to remote-sensed datasets (e.g. LTER vegetation databases).
– Interface with stakeholders (e.g. public policy relevant placement, attend policy-relevant conferences) and present initial dashboards/visualisations for feedback.
– Prepare modelling-empirical comparisons for publication.
Year 3.5
– Present at a major or international conference.
– Finalise remaining work for publication and integrate papers into thesis format.
– Prepare for viva.
Training
& Skills
Scientific training
The student will learn advanced techniques in dynamical systems and partial differential equations (Patterson, Krause), essential for building and analysing complex vegetation models. Training will also include ecological and geomorphological modelling (Turnbull-Lloyd, Wainwright), and data analysis covering high-performance computing and remote-sensed vegetation modelling (Sanderson). This broad training is valuable in environmental sciences and other multidisciplinary fields.
Opportunities for career development
Students in the supervisory team’s groups have co-supervised undergraduates, taught practicals, participated in outreach, and held leadership roles. The candidate will present at conferences and co-author publications. Training in interdisciplinary communication will be highlighted, benefiting future academic, government, or industry careers, alongside placement with a governmental or NGO partner from our network of contacts, which will strengthen stakeholder engagement and provide policy experience.
Research environment
The student joins a multidisciplinary team spanning environmental science and mathematical modelling, providing broad perspectives on climate-relevant science and multidisciplinary research. Collaboration with Patterson’s and Krause’s groups involves the Applied Mathematics Section (10 faculty, 7 postdocs, 5 PhD students), which hosts weekly lunches and seminars. Involvement with the Biophysical Sciences Institute offers cross-disciplinary training, outreach, and enrichment. The student will also join weekly seminars with Newcastle’s Modelling, Evidence & Public Policy group (5 postdocs, 10 PhD students) led by Sanderson.
To contextualise research, Turnbull-Lloyd and Wainwright provide guidance on dryland ecogeomorphic processes—drawing on experience in the Mediterranean, Africa, and North America—and train in extracting vegetation and fire data from remote sensing. Dr Sanderson will guide the use of remote-sensed and existing datasets for multi-scale vegetation modelling, maximising prediction utility for sustainable management.
Given the project’s policy impact potential, the student is encouraged to join outreach events (e.g., “Pint of Science”, Durham’s “Celebrate Science”). The student will present at key mathematics and ecology conferences, including the British Ecological Society meeting, the British Applied Mathematics Colloquium, and the Society for Mathematical Biology Annual meeting.
References & further reading
Butler & Sanderson (2021) A rapid method for quantifying small‐scale vegetation patch structure to complement conventional quadrat surveys. Appl. Veg. Sci., 24, e12610
Butler & Sanderson (2022). National-scale predictions of plant assemblages via community distribution models: leveraging published data to guide future surveys. J. Appl. Ecol., 59, 1559–1571
Stewart, Parsons, Wainwright, Okin, Bestelmeyer, Fredrickson & Schlesinger (2014). Modeling emergent patterns of dynamic desert ecosystems. Ecol. Monogr., 84(3), 373–410.
Meron (2011). Modeling dryland landscapes. Math. Model. Nat. Phenom., 6(1), 163–187.
Patterson, Levin, Staver & Touboul (2020). Probabilistic foundations of spatial mean-field models in ecology and applications. SIAM J. Appl. Dyn. Syst., 19(4), 2682–2719.
Tiwari, Turnbull & Wainwright (2025). Local- and global-scale hydrological and sediment connectivity over grassland and shrubland hillslopes. J. Hydrol., 655, 132896.
Xu, Patterson, Levin & Wang (2023). Non-equilibrium early-warning signals for critical transitions in ecological systems. PNAS, 120(5), e2218663120.
