Spatial Congruence of Biodiversity and Ecosystem Services: From Foundational Theory to Biomedical Applications

Abstract

This article synthesizes current research on the spatial relationships between biodiversity and ecosystem services (ES), a field critical for informed conservation planning and sustainable resource management. We explore the foundational theory behind B-ES relationships, reviewing evidence that shows variable spatial congruence—from high overlap to poor correlation—ac different services and regions. The article provides a comprehensive overview of methodological approaches for spatial prioritization and hotspot identification, alongside a critical examination of common challenges such as inappropriate proxies and scale-dependency. Furthermore, we present validation studies and comparative analyses that assess the real-world effectiveness of different conservation strategies. Finally, the conclusion discusses the implications of these ecological findings for biomedical and clinical research, particularly in the context of drug discovery from natural compounds and the role of ecosystem services in supporting human health.

The Foundation: Unraveling the Spatial Links Between Biodiversity and Ecosystem Services

Biodiversity, ecosystem functioning, and ecosystem services form a critical hierarchical chain that supports environmental and human health [1]. Biodiversity represents the variety of life at genetic, species, and ecosystem levels, while ecosystem functions encompass the biological, geochemical, and physical processes that occur within ecosystems [2]. These functions, in turn, underpin final ecosystem services—the directly consumed benefits that humanity derives from nature [3] [4]. Understanding these interrelationships is fundamental for spatial congruence analysis, which examines how these layers align geographically to inform conservation and management decisions [5] [6].

The conceptual framework in Figure 1 illustrates the hierarchical organization from biodiversity to human wellbeing, guided by management actions. This framework is essential for analyzing spatial congruence between biodiversity and ecosystem services.

Figure 1. Conceptual hierarchy from biodiversity to human wellbeing.

Defining the Core Conceptual Framework

Biodiversity: The Foundation of Life's Variability

Biodiversity encompasses the full spectrum of biological variety, defined by the Convention on Biological Diversity as "the variability among living organisms from all sources including, inter alia, terrestrial, marine and other aquatic ecosystems and the ecological complexes of which they are part; this includes diversity within species, between species and of ecosystems" [4]. This multifaceted concept operates across three primary levels:

• Genetic diversity: The variation of genetic information within and between populations of species [7]
• Species diversity: The variety of species within a habitat or region, including richness, evenness, and taxonomic differences [8]
• Ecosystem diversity: The diversity of habitats, biotic communities, and ecological processes within a region [7]

Biodiversity metrics extend beyond simple species counts to include functional traits that influence ecosystem processes, and phylogenetic relationships that reflect evolutionary history [4]. Recent research has revealed a universal "core-to-transition" organization of biodiversity across biogeographical regions, with gradients in richness, endemicity, and biota overlap forming consistent spatial patterns [9].

Ecosystem Functions and Processes: The Mechanisms of Nature

Ecosystem functioning refers to "the joint effects of all processes (fluxes of energy and matter) that sustain an ecosystem" through biological activities over time and space [3]. These processes include:

• Energy flows: Primary production, respiration, and trophic transfers
• Biogeochemical cycling: Nutrient cycling, decomposition, and water purification
• Biological processes: Pollination, seed dispersal, and pest regulation
• Physical processes: Soil formation, sediment retention, and water regulation [3]

These functions are regulated by both abiotic factors (temperature, moisture, pH) and biotic drivers (species traits, interactions) operating across multiple spatial and temporal scales [1] [3]. The relationship between biodiversity and ecosystem function (B-EF) often exhibits non-additive effects, where the loss of key species can disproportionately impact functional capacity due to complementary resource use and facilitative interactions [3].

Final Ecosystem Services: Benefits to Humanity

Final ecosystem services represent the subset of ecosystem components and processes that are directly enjoyed, consumed, or used to yield human wellbeing [4]. They are distinguished from "intermediate services" (which are essentially ecosystem functions) by their direct contribution to human benefits [4]. The Millennium Ecosystem Assessment categorized final services into four primary types:

• Provisioning services: Material outputs from ecosystems (food, water, timber, medicinal resources)
• Regulating services: Benefits from regulation of ecosystem processes (climate regulation, water purification, pollination)
• Cultural services: Non-material benefits (aesthetic, spiritual, educational, recreational)
• Supporting services: underlying processes necessary for producing all other services (primary production, nutrient cycling) [3]

Table 1: Quantitative Evidence of Biodiversity's Role in Delivering Final Ecosystem Services

Final Ecosystem Service	Quantitative Measure	Biodiversity Linkage	Scale
Food Provision	>75% of global food crops rely on pollinators [2]	Pollinator diversity increases crop yield stability [10]	Global
Climate Regulation	Forests store ~2.6 billion tonnes of CO₂ annually [2]	Plant diversity enhances carbon sequestration [1]	Global
Water Purification	Wetlands provide 75% of global freshwater resources [2]	Microbial and plant diversity enhances filtration [3]	Regional
Medicinal Resources	>50% of modern medicines from natural sources [2]	Genetic diversity provides biochemical variety [7]	Global
Disease Regulation	>75% of emerging infectious diseases are zoonotic [2]	Biodiversity reduces pathogen transmission [1]	Landscape

Experimental Protocols for Spatial Congruence Analysis

Protocol 1: Assessing Biodiversity and Ecosystem Service Relationships

Objective: Quantify spatial relationships between biodiversity metrics and final ecosystem service delivery using modeling approaches [6].

Workflow:

• Biodiversity Assessment:
  • Collect species occurrence data from field surveys, museum records, or satellite observations [9] [6]
  • Model species distributions using Maximum Entropy (MaxEnt) modeling with environmental predictors [6]
  • Calculate biodiversity indices (species richness, functional diversity, phylogenetic diversity) using Zonation software or similar prioritization tools [6]

• Ecosystem Service Quantification:

  • Map final ecosystem services using the InVEST (Integrated Valuation of Ecosystem Services and Tradeoffs) model suite [6]
  • Validate model outputs with field measurements (e.g., water quality monitoring, soil sampling) [6]
  • Apply random forest machine learning algorithms to analyze nonlinear relationships between biodiversity and ecosystem services [6]

• Spatial Congruence Analysis:

  • Conduct spatial overlay analysis to identify areas of high/low congruence
  • Calculate correlation coefficients between biodiversity and ecosystem service layers
  • Use partial dependence plots (PDP) to visualize marginal effects of biodiversity on services [6]

Application Note: This protocol was successfully applied in the Wujiang River Basin, revealing that water-related services (water supply, water conservation) showed the strongest positive correlations with biodiversity, contributing 66.24% and 44.72% respectively to biodiversity patterns [6].

Protocol 2: Measuring Biodiversity-Ecosystem Function Relationships

Objective: Experimentally quantify how biodiversity components regulate specific ecosystem functions under field conditions [10] [3].

Workflow:

• Experimental Design:
  • Establish sampling plots across biodiversity gradients (natural or manipulated)
  • Measure key ecosystem functions (e.g., decomposition, productivity, nutrient cycling)
  • Quantify biodiversity across multiple dimensions (taxonomic, functional, phylogenetic)

• Field Measurements:

  • Primary Productivity: Measure aboveground biomass harvesting and root sampling
  • Decomposition: Use standardized litter bags with mass loss measurements
  • Nutrient Cycling: Analyze soil nutrient pools, microbial biomass, and process rates
  • Pollination: Conduct pollinator exclusion experiments and fruit set monitoring [10]

• Statistical Analysis:

  • Fit biodiversity-ecosystem function (BEF) relationships using generalized additive models
  • Test for non-additive effects via interaction terms in mixed-effects models
  • Identify trait-based mechanisms using structural equation modeling [3]

Application Note: This approach revealed in alpine meadows that moderate grazing significantly increased water use efficiency (measured via leaf δ13C), demonstrating how biodiversity mediates ecosystem responses to management [1].

The Scientist's Toolkit: Essential Research Reagents and Solutions

Table 2: Key Research Reagent Solutions for Biodiversity and Ecosystem Service Analysis

Tool/Category	Specific Examples	Function/Application	Key References
Biodiversity Assessment	Zonation, MaxEnt, Infomap	Spatial prioritization, species distribution modeling, bioregion delineation	[9] [6]
Ecosystem Service Modeling	InVEST model suite, ARIES	Quantifying and mapping final ecosystem services	[6]
Functional Trait Databases	TRY Plant Trait Database, EltonTraits	Providing species-level functional trait data	[4]
Genetic Analysis	DNA barcoding kits, metabarcoding primers	Assessing genetic and microbial diversity	[1]
Remote Sensing Data	Landsat, Sentinel-2, MODIS	Landscape-scale biodiversity and function monitoring	[9] [6]
Statistical Analysis	R packages (vegan, FD, randomForest)	Analyzing biodiversity-ecosystem service relationships	[6] [4]

Analytical Framework for Spatial Congruence Research

The relationship between biodiversity and final ecosystem services exhibits complex spatial patterns that require sophisticated analytical frameworks. Research in the Wujiang River Basin demonstrated that different ecosystem services show varying spatial congruency with biodiversity, with water-related services showing the strongest alignment [6]. The random forest and partial dependence plot machine learning models revealed nonlinear relationships, emphasizing that biodiversity's influence on services is context-dependent [6].

Table 3: Spatial Congruence Patterns Between Biodiversity and Ecosystem Services

Congruence Pattern	Description	Management Implications	Identified In
High Congruence	Biodiversity hotspots coincide with multiple ecosystem service hotspots	Priority areas for conservation with co-benefits	[5] [6]
Trade-off Areas	High biodiversity but low service provision, or vice versa	Require careful planning to balance objectives	[5]
Spatial Mismatch	Protected areas not optimally located for both biodiversity and services	Need for systematic conservation planning	[5]
Scale Dependency	Congruence varies across spatial scales	Multi-scale assessment necessary	[4]
Service-Specific Relationships	Varying biodiversity-service correlations across service types	Service-specific management strategies	[6]

The integration of metacommunity and meta-ecosystem perspectives is crucial for spatial congruence analysis, as connectivity between habitats influences both biodiversity patterns and ecosystem service flows [3]. Recent findings of a universal "core-to-transition" organization in biogeographical regions [9] provide a template for predicting how biodiversity and services may be distributed across landscapes.

Defining the core concepts of biodiversity, ecosystem functions, and final ecosystem services provides the foundation for spatial congruence analysis in ecological research. The experimental protocols and analytical frameworks presented here enable researchers to quantify these relationships and identify areas where conservation can maximize joint benefits. As the field advances, integrating trait-based approaches, multi-scale analyses, and mechanistic understanding of B-EF relationships will enhance our ability to predict and manage the complex interconnections between nature's diversity and human wellbeing.

A core pursuit in biogeography and conservation science is understanding the spatial relationships between biodiversity and ecosystem services. The central debate hinges on spatial congruence—the degree to which areas of high biodiversity overlap with areas that provide critical benefits to human well-being. Some studies suggest strong overlaps, advocating for efficient conservation planning, while others find weak or mismatched patterns, indicating potential trade-offs.

This document provides application notes and protocols for analyzing these spatial relationships, enabling researchers to quantify congruence and interpret its implications for environmental management and policy. The content is framed within a broader thesis on biodiversity and ecosystem service spatial congruence analysis research, offering standardized methods for researchers, scientists, and drug development professionals whose work intersects with natural products and ecological impacts.

Key Concepts and Analytical Framework

Defining Spatial Congruence and Chorotypes

In biogeography, groups of species with congruent geographical distributions, potentially caused by common historical and ecological spatial processes, constitute chorotypes [11]. The degree of spatial range congruence is a critical, explicit parameter in this analysis. The Spatial Congruence Analysis (SCAN) method identifies chorotypes by mapping direct and indirect spatial relationships among species, using a quantifiable measure of congruence [11].

When applied to the relationship between biodiversity and ecosystem services, the concept of a "chorotype" can be extended to describe areas where multiple ecosystem services and biodiversity metrics show spatially congruent patterns.

The Quantitative Measure of Congruence

Determining spatial congruence is non-trivial because ranges may differ in position, area, and shape. A generic spatial index, analogous to the Jaccard index, calculates congruence between two species (or metrics) as the product of the area of overlap weighted by the relative area of each [11]. The formula for this congruence index is presented in the experimental protocols below.

Experimental Protocols for Congruence Analysis

Protocol 1: Spatial Congruence Analysis (SCAN) for Chorotype Identification

Purpose: To identify groups of species, ecosystem services, or biodiversity metrics that share statistically significant congruent geographical distributions.

Principles: The method uses a one-layered network where vertices represent species (or metrics) and edges represent pairwise spatial congruence estimates. The network is analyzed for each reference species separately by an algorithm that searches for spatial relationships [11].

Workflow:

• Data Preparation: Obtain polygon (shapefile) range maps for all species or spatial data layers for ecosystem services and biodiversity metrics. Data should be in the same projection and coordinate system.
• Calculate Pairwise Congruence: For every pair of entities, compute the spatial congruence index: Congruence (A,B) = (Area of Overlap) / (Area of A) * (Area of Overlap) / (Area of B) This yields a value between 0 (no overlap) and 1 (perfect congruence) [11].
• Set Congruence Threshold: Define a reference value of congruence as an explicit numerical parameter (e.g., 90%). This threshold determines the minimum congruence required to establish a spatial link.
• Network Analysis: Initiate the analysis with a reference entity. Connect it to all other entities that meet or exceed the congruence threshold.
• Iterative Grouping: For each connected entity, repeat the comparison process, accruing more entities to the group. Continue until no new entities are found that meet the congruence threshold, forming a "closed list" or "partial chorotype."
• Multi-Threshold Analysis: Repeat steps 3-5 at progressively lower congruence thresholds (e.g., 85%, 80%). The set of all partial chorotypes derived from a particular reference entity, across all thresholds, constitutes the final "chorotype."

Visualization: The following workflow diagram illustrates the SCAN protocol:

Protocol 2: Asset-Based Impact Assessment for Corporate Biodiversity Footprinting

Purpose: To quantify the direct impacts of corporate physical assets (e.g., manufacturing facilities, mines) on ecosystem services and biodiversity, enabling transparent ESG reporting and within-sector comparisons.

Principles: This approach leverages high-spatial-resolution, global ecosystem service models and satellite imagery to map the footprints of physical assets and assess their impact relative to a baseline of potential natural vegetation [12]. It is analogous to Scope 1 greenhouse gas emissions accounting.

Workflow:

• Asset Mapping: Compile a global database of corporate physical assets with precise geolocation. Use high-resolution satellite imagery (e.g., Sentinel-2, Landsat) to map individual asset footprints accurately [12].
• Global Baseline Mapping: Create global maps for selected ecosystem service and biodiversity metrics under a scenario where all land is returned to its potential natural vegetation. Example metrics include:
  • Ecosystem Services: Coastal risk reduction, sediment retention, nitrogen retention, nature access.
  • Biodiversity: Species richness, habitat for Red List species, habitat for endemic species, Key Biodiversity Areas (KBAs) [12].
• Impact Quantification: For each asset, assume that development results in the complete loss of ecosystem service or biodiversity values within its footprint, compared to the natural baseline. The impact is the value of the metric at that location multiplied by the area of the footprint.
• Aggregation: Scale impacts from individual assets to the corporate level by summing the impacts of all assets owned by a company.
• Normalization: Normalize total impacts by company revenue or other operational metrics to enable fair cross-company comparisons.

Visualization: The following diagram outlines the asset-based impact assessment process:

Quantitative Data on Spatial Congruence

Variation in Corporate Impacts on Nature

Analysis of over 585,000 physical assets belonging to 2,173 global companies revealed significant variation in impacts on ecosystem services and biodiversity, both between and within sectors [12].

Table 1: Sector-Level Variation in Direct Impacts on Nature (Average per Company)

Sector	Relative Impact on Ecosystem Services	Relative Impact on Biodiversity	Key High-Impact Metrics
Utilities	High	High	Nitrogen retention, Sediment retention
Real Estate	High	Medium-High	Coastal risk reduction, Species richness
Materials	High	High	Sediment retention, Red List species habitat
Financials*	Medium-High	Medium-High	(Varies based on asset types in portfolio)
Energy	Medium	Medium	Nitrogen retention, KBAs
Consumer Staples	Low-Medium	Low-Medium	Nature access, Species richness

*Financial sector impacts are derived from their holdings of physical assets [12].

Table 2: Comparative Spatial Congruence of Ecosystem Services and Biodiversity

Metric	Geographic Concentration (Land area in 90th percentile for ≥3 metrics)	Key Global Hotspots
Ecosystem Services	0.5% (~700,000 km²)	Coastal zones, areas near streams/rivers, India, China, Europe, Midwestern US, East Africa, global urban areas [12].
Biodiversity	3.2% of total land area	Amazon, Southeast Asia, Guinean Forests of West Africa, Eastern Afromontane region [12].

The Scientist's Toolkit: Research Reagent Solutions

Table 3: Essential Materials and Tools for Spatial Congruence Research

Item	Function/Benefit
Global Species Range Maps (e.g., IUCN Red List, BirdLife International)	Provides standardized polygon data for species distributions, essential for calculating range congruences [11].
High-Resolution Satellite Imagery (e.g., Sentinel-2, Landsat)	Enables precise mapping of asset footprints and land-use change, moving beyond point locations to accurate area calculations [12].
Process-Based Ecosystem Service Models (e.g., InVEST)	Open-source, spatially explicit models that predict the provision and value of ecosystem services under different land-use scenarios [12].
R Statistical Environment with sf package	A computational platform for performing spatial analyses, including the implementation of methods like SCAN for calculating congruence indices and managing shapefiles [11].
Corporate Asset Databases (e.g., Orbis, FactSet)	Provides the foundational list and location (latitude/longitude) of physical assets owned by companies for large-scale impact assessments [12].

Application Note: Comparative Analysis of Spatial Congruence

This application note synthesizes methodologies and findings from two pivotal case studies examining the spatial relationships between biodiversity and ecosystem services. The South African national-scale assessment and the Baiyangdian Basin study in China provide complementary insights for researchers conducting spatial congruence analysis. Both studies reveal complex, scale-dependent relationships with significant implications for conservation planning and ecosystem management.

Study Contexts and Key Findings

Table 1: Comparative Overview of Case Study Characteristics

Characteristic	South Africa Case Study	Baiyangdian Basin Case Study
Geographic Scale	National (1.22 million km²)	Watershed (31,200 km²)
Primary Ecosystem Services Analyzed	Water flow, carbon storage, soil accumulation, nutrient retention	Water yield, nutrient export, carbon storage, habitat quality
Biodiversity Metrics	Species richness, biodiversity priority areas	Habitat quality (as proxy for biodiversity)
Spatial Congruence Finding	Moderate overlap; strongest in grasslands and savannas [13]	Significant trade-offs between regulating and supporting services [14]
Key Management Implication	Ecosystem services can provide additional rationale for biodiversity conservation in specific biomes [13]	Targeted spatial strategies needed to address service-specific trade-offs [14]

Protocol: Spatial Congruence Analysis Framework

Stage 1: Data Collection and Preparation

• Ecosystem Service Quantification: Utilize standardized models (e.g., InVEST) to map service provision spatially [14]
• Biodiversity Data Compilation: Incorporate species richness data, habitat quality maps, and existing conservation priority areas [13] [15]
• Spatial Resolution Standardization: Process all spatial data to common resolution and projection system (e.g., 100m grid cells) [14]

Stage 2: Hotspot Identification and Analysis

• Service-Specific Hotspots: Apply top-richest cells method (e.g., highest 10-30% of values) for individual ecosystem services [16]
• Multi-Service Hotspots: Identify areas of overlapping high provision for multiple services using spatial overlay techniques [15]
• Congruence Assessment: Calculate spatial correlation coefficients and overlap percentages between biodiversity priorities and ecosystem service hotspots [13]

Stage 3: Trade-off and Synergy Analysis

• RMSE Method Implementation: Quantify trade-off strength using Root Mean Square Error method to capture uneven rates of change between ecosystem services [14]
• Spatial Heterogeneity Mapping: Apply geographical detector and MGWR models to identify spatially varying relationships between services [14]

Experimental Workflow Visualization

Spatial Congruence Analysis Workflow

Key Research Reagent Solutions

Table 2: Essential Analytical Tools for Spatial Congruence Research

Tool/Category	Specific Examples	Research Application
Spatial Analysis Platforms	ArcGIS, QGIS, GRASS	Core spatial data processing, overlay analysis, and mapping [14]
Ecosystem Service Models	InVEST, ARIES, SOLVES	Quantification and mapping of ecosystem service provision [14]
Statistical Analysis Tools	R, Python, MATLAB	Correlation analysis, RMSE calculations, regression modeling [14]
Prioritization Software	Marxan, Zonation	Systematic conservation planning incorporating biodiversity and services [16]
Geospatial Data	Land cover maps, DEM, soil data, climate data	Input parameters for ecosystem service models [14]

Protocol: Integrated Spatial Assessment Methodology

Field Data Collection Specifications

• Land Use/Land Cover Validation: Conduct ground-truthing with GPS-enabled field surveys at minimum 50 validation points per major land cover class
• Biodiversity Sampling: Implement standardized plot-based surveys for species richness assessment, following protocols appropriate to taxa [13]
• Hydrological Measurements: For watershed studies, establish monitoring stations for water quantity and quality parameters [14]

Ecosystem Service Modeling Protocol

Water Yield Assessment

• Model: InVEST Seasonal Water Yield Model [14]
• Required Parameters: Precipitation, evapotranspiration, soil depth, plant available water content, land use/land cover, watershed boundaries [14]
• Validation Method: Compare modeled outputs with stream gauge measurements where available

Carbon Storage Assessment

• Model: InVEST Carbon Storage and Sequestration Model [14]
• Pool Estimation: Aboveground biomass, belowground biomass, soil organic matter, dead organic matter
• Data Sources: Field measurements, literature values for ecosystem-specific carbon densities

Spatial Congruence Analysis Protocol

Hotspot Identification

• Classification Method: Apply quantile-based classification (e.g., top 20% of values) for each ecosystem service and biodiversity metric [16]
• Sensitivity Analysis: Test multiple threshold levels (10%, 20%, 30%) to assess robustness of findings [16]

Overlap Analysis

• Spatial Correlation: Calculate Pearson or Spearman correlation coefficients between biodiversity and ecosystem service layers [13]
• Percent Overlap: Quantify the percentage of biodiversity priority areas that coincide with ecosystem service hotspots [13] [15]

Trade-off Analysis Visualization

Ecosystem Service Trade-off Drivers

Data Integration and Synthesis Protocol

• Multi-Scale Analysis: Conduct analyses at multiple spatial scales (local, watershed, regional) to assess scale dependencies [13] [14]
• Uncertainty Propagation: Document and quantify uncertainty sources throughout the analytical chain
• Stakeholder Validation: Present preliminary findings to local experts and stakeholders for contextual interpretation [15]

Application to Conservation Planning

Priority Area Selection Framework

The case studies demonstrate that spatial congruence between biodiversity and ecosystem services is partial rather than complete. In South Africa, approximately 30% of ecosystem service hotspots were located within terrestrial critical biodiversity areas [15], suggesting that conservation strategies focused solely on biodiversity would miss significant ecosystem service values. The Baiyangdian study revealed that trade-offs intensified in specific geographical contexts, particularly in southeast watershed areas where built-up lands are concentrated [14].

Management Recommendations

• Integrated Prioritization: Combine biodiversity and ecosystem service criteria in systematic conservation planning [16]
• Spatially Explicit Interventions: Target specific management practices to address location-specific trade-offs [14]
• Policy Integration: Mainstream congruence analysis into land use planning and natural resource management decisions [13] [15]

This protocol provides a standardized methodology for analyzing spatial relationships between biodiversity and ecosystem services, enabling comparative studies across different regions and ecosystems. The integration of these approaches strengthens conservation planning by simultaneously addressing biodiversity protection and ecosystem service maintenance.

Spatial congruence analysis represents a critical methodology for identifying priority areas where high biodiversity and essential ecosystem services co-occur. Understanding these synergies is paramount for researchers and policymakers aiming to optimize conservation strategies that deliver dual benefits for ecological integrity and human wellbeing. The complex, scale-dependent relationships between biodiversity and ecosystem services necessitate sophisticated spatial analysis techniques to disentangle patterns and drivers across landscapes [17]. This protocol provides a comprehensive framework for conducting such analyses, enabling the identification of true win-win areas that can inform everything from local conservation planning to corporate environmental impact assessments [12].

Recent research has demonstrated that while positive relationships often exist between biodiversity attributes and ecosystem services, these associations are highly complex and service-dependent [18]. The spatial patterns exhibit significant scale variations, with different factors dominating at different scales—anthropogenic influences typically controlling relationships at finer scales (e.g., 12 km) while physical environmental factors dominate at broader scales (e.g., 83 km) [17]. This protocol addresses these complexities through a multi-scale approach that captures the hierarchical nature of biodiversity-ecosystem service relationships.

Application Notes: Key Principles and Considerations

Foundational Concepts

• Spatial Congruence: Defined as the geographic overlap between areas of high biodiversity value and areas of high ecosystem service provision. True congruence occurs when these areas overlap significantly beyond chance expectations.
• Scale Dependence: Relationships between biodiversity and ecosystem services vary across spatial scales, necessitating multi-scale analyses to avoid misleading conclusions [17].
• Metric Selection: Biodiversity and ecosystem services must be measured using multiple complementary indicators, as single metrics often fail to capture the full complexity of these relationships [12].

Practical Implementation Contexts

The methodologies outlined in this protocol apply to multiple research and application contexts:

• Conservation Planning: Identifying priority areas that simultaneously protect biodiversity and maintain critical ecosystem services.
• Corporate Impact Assessment: Evaluating the effects of business operations on nature using standardized, comparable metrics [12].
• Policy Development: Informing spatial planning decisions and biodiversity offset strategies.
• Research Applications: Advancing understanding of the functional relationships between biodiversity components and ecosystem service flows.

Experimental Protocols

Phase I: Data Acquisition and Preparation

Protocol 1.1: Spatial Data Collection

Purpose: To acquire foundational spatial datasets for biodiversity and ecosystem services assessment.

Procedure:

• Biodiversity Data Compilation
  • Obtain species distribution data from global databases (GBIF, IUCN Red List)
  • Acquire habitat maps and ecosystem classifications
  • Collect data on protected areas and key biodiversity areas
  • Compile metrics for species richness, threatened species, endemic species, and key biodiversity areas [12]

• Ecosystem Services Quantification

  • Select appropriate ES metrics based on research objectives (see Table 1)
  • Utilize process-based models for ES mapping (e.g., InVEST, ARIES)
  • Employ biophysical indicators for ES quantification (see Table 2)
  • Validate model outputs with field measurements and literature values

• Environmental Covariate Collection

  • Compile topographic variables (elevation, slope, aspect)
  • Obtain climate data (temperature, precipitation, seasonality)
  • Collect soil properties (texture, pH, organic matter)
  • Acquire land use/land cover classifications
  • Gather socio-economic data where relevant

Quality Control: Perform cross-validation with independent datasets and sensitivity analysis on model parameters.

Protocol 1.2: Data Pre-processing

Purpose: To standardize spatial datasets for integrated analysis.

Procedure:

• Spatial Alignment
  • Establish common spatial extent for all datasets
  • Resample to consistent spatial resolution (e.g., 30m, 100m, 1km)
  • Apply consistent coordinate reference system
  • Verify spatial alignment accuracy

• Data Transformation
  • Apply appropriate normalization techniques where needed
  • Conduct multicollinearity assessment for predictor variables
  • Implement spatial autocorrelation analysis
  • Create standardized measurement units across datasets

Troubleshooting: Address spatial mismatch through uncertainty analysis and documentation of limitations.

Phase II: Spatial Analysis and Congruence Identification

Protocol 2.1: Multi-scale Spatial Analysis

Purpose: To identify spatial patterns and relationships across multiple scales.

Procedure:

• Factorial Kriging Analysis (FKA)
  • Fit a linear model of co-regionalization (LMC) to the data
  • Decompose total variation into multiple spatial components
  • Identify key spatial scales (ranges) of variability
  • Calculate spatial components for each scale identified [17]

• Spectral Analysis

  • Perform wavelet analysis to detect scale-specific patterns
  • Apply Fourier transforms for periodic pattern identification
  • Implement multi-resolution analysis for hierarchical decomposition

• Cross-scale Correlation Analysis

  • Calculate correlation coefficients between biodiversity and ES at each scale
  • Identify scale-dependent relationships and thresholds
  • Document consistency or variation in relationships across scales

Technical Notes: FKA requires substantial computational resources for large datasets; consider cloud computing for extensive study areas.

Protocol 2.2: Congruence Quantification

Purpose: To measure and map the spatial overlap between biodiversity and ecosystem services.

Procedure:

• Hotspot Identification
  • Classify areas into percentile-based categories (e.g., top 10%, 25%)
  • Identify areas of high value for both biodiversity and ecosystem services
  • Calculate spatial congruence indices
  • Generate congruence maps visualizing overlap patterns

• Statistical Significance Testing

  • Implement randomization tests to assess significance of overlap
  • Calculate null distributions through spatial permutation
  • Determine probability values for observed congruence patterns
  • Account for spatial autocorrelation in significance testing

• Uncertainty Propagation

  • Quantify uncertainty from source data through error propagation
  • Perform sensitivity analysis on classification thresholds
  • Generate confidence intervals for congruence metrics
  • Document limitations and assumptions

Validation: Use independent datasets and case studies to validate congruence patterns.

Data Presentation

Table 1: Ecosystem Service Metrics for Spatial Congruence Analysis

Ecosystem Service Category	Specific Metric	Measurement Unit	Data Sources	Spatial Resolution	Temporal Resolution
Water Quality Regulation	Nitrogen retention	kg N retained/ha/year	MODIS, Landsat	30m-1km	Annual
	Phosphorus retention	kg P retained/ha/year	Soil surveys, DEM	30m-100m	Annual
Water Flow Regulation	Sediment retention	tons sediment retained/ha/year	USLE/RUSLE models	30m	Annual
Crop Production	Yield potential	kg/ha/year	Agricultural statistics	100m-1km	Annual
Nature Access	Population within 1-hour travel	Number of people	Population data, travel time	100m	5-year updates
Coastal Protection	Coastal risk reduction	People protected from flooding	Coastal models, population	30m-100m	Static

Table 2: Biodiversity Metrics for Spatial Congruence Analysis

Biodiversity Dimension	Specific Metric	Measurement Unit	Data Sources	Spatial Resolution	Key Considerations
Taxonomic Diversity	Species richness	Number of species/unit area	GBIF, IUCN	1km-10km	Sampling bias correction
Conservation Priority	Red List species habitat	Presence/quality of habitat	IUCN Red List	100m-1km	Habitat specificity
Endemism	Endemic species richness	Number of endemic species	Regional assessments	1km-100km	Scale dependency
Key Areas	Key Biodiversity Areas	Presence/overlap with KBAs	KBA Secretariat	Polygon-based	Categorical data

Spatial Overlap Category	Ecosystem Services	Biodiversity	Interpretation	Management Implications
High Overlap Areas	0.5% of land area in 90th percentile for ≥3 services	3.2% of land area in 90th percentile for ≥3 metrics	Limited geographic congruence	Highest priority for conservation
Sectoral Impact Variation	Utilities, real estate, materials have largest impacts	Financial sectors also have large impacts	Cross-sectoral variation	Sector-specific guidelines needed
Spatial Concentration	High heterogeneity, varies over small distances	Greater concentration in tropics	Different spatial patterns	Requires high-resolution data

Visualization

Spatial Analysis Workflow

Multi-scale Analytical Framework

The Scientist's Toolkit: Research Reagent Solutions

Table 4: Essential Analytical Tools and Platforms

Tool/Platform Category	Specific Solution	Primary Function	Application Context	Data Compatibility
Spatial Analysis Software	ArcGIS 10.2+	Geographic information system	Data processing, mapping, analysis	Multi-format, raster/vector
	QGIS 3.0+	Open-source GIS alternative	Cost-effective spatial analysis	Multi-format, raster/vector
Statistical Programming	R 3.0.3+ with tidytext, topicmodels	Statistical computing, text mining	Data analysis, topic modeling [19]	CSV, spatial formats
	Python with geospatial libraries	Custom analysis, machine learning	Advanced spatial statistics	Multiple data formats
Ecosystem Service Modeling	InVEST 3.0+	Ecosystem service mapping	ES quantification, scenario analysis	Spatial data inputs
	ARIES	Artificial intelligence for ES	Rapid ES assessment	Cloud-based integration
Biodiversity Assessment	Species Population Models	Population viability analysis	Species-specific impact assessment	Field data, habitat maps
	IBAT	Integrated Biodiversity Assessment	Screening, risk assessment	Global biodiversity datasets
Satellite Imagery	Sentinel-2	Land cover classification	Change detection, habitat mapping	10m resolution
	Landsat 8/9	Historical land use analysis	Long-term trend assessment	30m resolution
Corporate Assessment	Global Asset Impact Platform	Corporate impact quantification	ESG reporting, impact assessment [12]	Company asset data

Understanding the mechanisms through biodiversity influences ecosystem functioning is a cornerstone of ecology, with profound implications for conservation, ecosystem service management, and environmental policy. The relationship extends beyond simple counts of species to encompass how the identity, traits, and spatial distributions of organisms govern ecological processes. Recent research has increasingly highlighted that spatial organization is not merely a backdrop but a fundamental driver of these biodiversity-ecosystem function (BEF) relationships. A groundbreaking 2025 study examining over 30,000 species of vertebrates, invertebrates, and plants revealed a universal core-to-transition organization within biogeographical regions [9]. This organization reflects ordered gradients in species richness, range size, endemicity, and biogeographical transitions, which in turn are shaped by environmental filters acting on regional hotspots and permeable biogeographical boundaries [9]. This spatial perspective provides a critical framework for analyzing the mechanistic underpinnings of BEF relationships, suggesting that general forces operating across the tree of life shape how biodiversity is organized and consequently how it functions [9].

Complementing this macro-scale perspective, a 2025 mechanistic partitioning of biodiversity effects demonstrated that ecosystem stability—encompassing both resistance and resilience to perturbations like drought—is governed by distinct selection and complementarity effects [20]. Biodiversity can enhance ecosystem resistance both by making species assemblages less sensitive to perturbations and through selection effects that favor less sensitive species during ecological assembly [20]. Similarly, ecosystem resilience is influenced by a complex interplay of species-level resilience traits and selection effects that benefit either highly resilient species or those with compensatory resistance-recovery trade-offs [20]. Together, these advances underscore that the mechanistic underpinnings of BEF relationships operate across spatial and organizational scales, from regional species pools to local species interactions and trait combinations.

Key Mechanistic Pathways and Experimental Evidence

Primary Mechanisms Linking Biodiversity to Ecosystem Function

Research has identified several fundamental mechanistic pathways through which biodiversity governs ecosystem functioning. The selection effect, also known as the sampling effect, occurs when diverse communities have a higher probability of containing and becoming dominated by species with particularly influential traits that drive ecosystem processes [20]. In contrast, complementarity effects emerge when species coexist and utilize resources in different ways, leading to more complete resource utilization and increased total productivity [20] [1]. These complementarity effects can result from various mechanisms including temporal niche partitioning (species being active at different times), spatial heterogeneity (species utilizing different microhabitats), and facilitative interactions where one species enhances the environment for others [1].

The interplay of these mechanisms generates ecosystem stability, which encompasses both resistance (the ability to withstand perturbation) and resilience (the speed of recovery after disturbance) [20]. The 2025 partitioning of net biodiversity effects revealed surprising complexity in these relationships—for instance, resilience can be enhanced not only by selection for highly resilient species but also through selection effects where changes in species resistance indirectly benefit less resilient species through competitive release or other community dynamics [20]. In grassland ecosystems subjected to drought, this mechanistic partitioning precisely explained how functional complementarity and differences in life-history strategies allow more diverse plant communities to both resist and swiftly recover from major disturbances [20].

Table 1: Core Mechanisms Linking Biodiversity to Ecosystem Function

Mechanism	Functional Principle	Ecosystem Outcome	Experimental Evidence
Selection Effect	Dominance by high-functioning species	Increased primary productivity; Enhanced resistance	Grassland drought experiments [20]
Complementarity Effect	Niche differentiation and facilitation	More complete resource use; Increased nutrient cycling	Multi-trophic diversity studies [1]
Resistance Selection	Assembly favors disturbance-tolerant species	Maintained function during perturbation	Drought modeling in grasslands [20]
Resilience Selection	Recovery benefits to resilient species	Faster functional recovery post-disturbance	Ecosystem stability partitioning [20]
Spatial Insurance	Cross-habitat subsidies and dispersal	Stable ecosystem services at landscape scales	Mangrove landscape studies [21]

Spatial Congruence as a Framework for BEF Analysis

The Spatial Congruence Analysis (SCAN) method provides a powerful analytical framework for detecting biogeographical patterns relevant to BEF research [22]. SCAN identifies chorotypes—groups of species with significantly congruent geographical distributions—by mapping direct and indirect spatial relationships among species using an explicitly defined congruence threshold [22]. This method quantifies range congruence between species without assumptions about historical homology, instead focusing on the role of strict or relaxed congruence limits in pattern recognition [22]. The spatial congruence index (CS) is calculated as: CSab = (Oab/Aa) × (Oab/Ab), where Oab is the area of overlap between species a and b, and Aa and Ab are their respective range areas [22].

Application of SCAN to South American bird distributions revealed that only a small portion of range overlaps is biogeographically meaningful, and chorotypes can vary from simple groups of highly congruent species to complex patterns with numerous alternative component species and spatial configurations [22]. These patterns offer insights about possible processes driving BEF relationships at different degrees of spatial congruence, with metrics such as congruence depth, species richness, and ratio between common and total areas allowing detailed comparisons between patterns across regions and taxa [22]. This approach is particularly valuable for identifying gradients of species distributions expanding from core areas, which aligns with the recently described core-to-transition organization of biodiversity [9].

Application Notes: Protocol for Spatial BEF Analysis

Protocol 1: Spatial Congruence Analysis of Biogeographical Patterns

Purpose: To identify groups of species with congruent distributions (chorotypes) that may represent functional units in ecosystem processes, using the Spatial Congruence Analysis (SCAN) method [22].

Materials and Software:

• Species distribution data (polygon shapefiles)
• R statistical environment with sf package [22]
• GIS software for visualization

Procedure:

• Data Preparation: Compile polygon range maps for all target species. Ensure consistent projection and coordinate reference systems.
• Calculate Pairwise Congruence: For each species pair, compute the Spatial Congruence Index (CS) using the formula: CSab = (Oab/Aa) × (Oab/Ab), where Oab is the overlap area, and Aa and Ab are the respective range areas [22].
• Set Congruence Threshold (CT): Begin with high congruence (e.g., 90%) and progressively lower the threshold in increments (e.g., 5% decreases).
• Identify Direct and Indirect Congruences: For each CT, identify species pairs with CS ≥ CT (direct congruence). Then identify indirectly connected species through chains of direct connections [22].
• Form Partial Chorotypes: For a reference species at a given CT, iteratively add all directly and indirectly congruent species until no new species are added, creating a "closed list" or partial chorotype [22].
• Compile Comprehensive Chorotype: Repeat for multiple CT values to capture the full set of partial chorotypes for each reference species, representing the complete chorotype across congruence levels.
• Calculate Chorotype Metrics: For each chorotype, compute:
  • Congruence Depth: Range of CT values at which the chorotype persists
  • Species Richness: Number of component species
  • Common Area Ratio: Ratio between commonly shared area and total chorotype area [22]

Interpretation: Complex chorotypes with numerous alternative species compositions and wide congruence depth suggest stronger historical or environmental constraints on community assembly, which may translate to more consistent BEF relationships. Simple chorotypes with narrow congruence depth may indicate more stochastic assemblages with potentially more variable ecosystem functioning.

Protocol 2: Quantifying Biodiversity Effects on Ecosystem Stability

Purpose: To partition net biodiversity effects on ecosystem resistance and resilience following the methodology established in recent mechanistic studies [20].

Experimental Design:

• Establish Diversity Gradient: Create experimental communities representing a gradient of species richness and functional diversity.
• Apply Perturbation: Implement a standardized disturbance (e.g., drought, heating, nutrient pulse) with appropriate control treatments.
• Monitor Response Trajectories: Measure key ecosystem processes (e.g., productivity, nutrient cycling) before, during, and after perturbation to quantify resistance and resilience.

Resistance Calculation: Resistance = 1 - (│ΔPerturbed - ΔControl│ / ΔControl), where Δ represents the change in ecosystem process rates during disturbance [20].

Resilience Calculation: Resilience = (Recovery ratepost-disturbance) / (Deviation from pre-disturbance baseline) [20].

Mechanistic Partitioning:

• Quantify Selection Effects: Compare the performance of species in mixture versus monoculture to identify cases where dominant species drive ecosystem responses.
• Identify Complementarity Effects: Calculate the extent to which species mixtures outperform monoculture expectations based on species-specific contributions.
• Partition Resistance and Resilience Pathways: Analyze how species traits and interactions contribute specifically to resistance during disturbance versus recovery afterward [20].

Statistical Analysis: Use structural equation modeling to test pathways linking biodiversity components to resistance and resilience mechanisms, and ultimately to ecosystem stability metrics.

Research Reagent Solutions: Essential Tools for BEF Research

Table 2: Key Research Materials and Analytical Tools for BEF Studies

Category	Specific Tool/Method	Research Application	Key References
Spatial Analysis	SCAN (Spatial Congruence Analysis)	Identifying chorotypes and biogeographical patterns	[22]
Bioregion Delineation	Infomap Bioregion algorithm	Detecting complex gradients in species distribution	[9] [22]
Ecosystem Stability Metrics	Resistance-Resilience Partitioning	Quantifying biodiversity effects on stability components	[20]
Field Monitoring	Harmonised biodiversity protocols (Biodiversa+) Standardized variables: species richness, abundance, biomass	Cross-scale comparable data collection	[23]
Taxonomic Reference	GBIF Backbone Taxonomy	Standardized species identification	[23]
Ecosystem Classification	IUCN Global Ecosystem Typology, EUNIS habitat classification	Consistent habitat definition	[23]

Conceptual Framework and Visual Synthesis

The relationship between biodiversity and ecosystem functioning operates through interconnected mechanisms that span spatial and temporal scales. The following diagram synthesizes the primary mechanistic pathways and their spatial context:

Figure 1: Mechanistic pathways linking biodiversity to ecosystem functioning through spatial organization.

The spatial organization of biodiversity creates a template upon which mechanistic processes operate. As illustrated in Figure 2, the core-to-transition structure of biogeographical regions establishes gradients in biodiversity aspects that influence how ecological functions are distributed across landscapes:

Figure 2: Core-to-transition organization of biodiversity aspects across biogeographical regions.

The mechanistic underpinnings of biodiversity-ecosystem functioning relationships are fundamentally spatial in nature. The recently identified core-to-transition organization of biodiversity across biogeographical regions provides a macro-scale template that constrains and guides local BEF relationships [9]. Simultaneously, mechanistic partitioning of biodiversity effects reveals that selection and complementarity operate in distinct ways to enhance ecosystem resistance and resilience [20]. The SCAN methodology offers a powerful tool for quantifying the spatial congruence patterns that underlie these relationships [22].

For researchers and conservation practitioners, these advances enable more predictive understanding of how biodiversity loss—particularly non-random extinction correlated with spatial patterns—will impact ecosystem functioning and service provision. The protocols and analytical frameworks presented here facilitate standardized assessment of BEF relationships across taxa and ecosystems, supporting evidence-based conservation prioritization and ecosystem management. Future research should focus on integrating these spatial and mechanistic perspectives to develop unified models that predict ecosystem responses to global change across scales, from local species interactions to continental biogeographical patterns.

Mapping and Prioritization: Methods for Spatial Analysis of Ecosystem Services

Spatial prioritization has emerged as a critical precursor to effective ecosystem service (ES) management, enabling decision-makers to identify locations where resource investment will yield the greatest return for both biodiversity conservation and human well-being [24]. This framework addresses the fundamental challenge of limited conservation resources by providing a systematic approach for identifying priority areas where place-based ES management should occur [24]. Unlike traditional conservation prioritization, spatial prioritization of ES must specifically account for human beneficiaries, the availability of alternatives to service provision, and the capacity to meet beneficiary demand [24]. The growing integration of ES considerations into environmental policies, particularly in rapidly developing regions [25], underscores the practical importance of robust spatial prioritization frameworks. This protocol outlines the key components and methodologies for implementing a comprehensive spatial prioritization approach that balances ecological and socio-economic objectives, with particular relevance to biodiversity and ecosystem service spatial congruence analysis research.

Theoretical Foundation: Core Components for ES Prioritization

The Ecosystem Service Cascade Framework

The ecosystem service cascade framework provides a conceptual foundation for spatial prioritization by illustrating how ecological structures and processes transform into ecosystem functions, then into ES, and finally into benefits that contribute to human well-being [26]. This framework distinguishes between individual components of the process—ecosystem properties and integrity, ES, human well-being, and exchange-value—thereby avoiding confusion regarding their relationships [26]. In spatial prioritization applications, the cascade framework helps researchers identify which components of the ES delivery process are being measured and mapped, ensuring a comprehensive approach that connects ecological supply to human benefits [26] [27].

Essential Spatial Prioritization Components

A comprehensive spatial prioritization framework for ES management must integrate four key components that distinguish it from traditional conservation prioritization [24]:

• Supply of Services: The biophysical capacity of an ecosystem to provide services
• Beneficiary Demand: The level of need or desire for ES by human populations
• Human-Derived Alternatives: The availability of technological or engineered substitutes for natural ES provision
• Site Dependency and Scale: The degree to which service delivery depends on specific locations and the spatial extent of service flow

These components interact within complex socio-ecological systems, creating trade-offs and synergies that must be explicitly considered in spatial planning [28] [24]. For instance, in Brazil's agricultural landscapes, trade-off analyses have revealed extreme contrasts between agricultural expansion and the conservation of biodiversity and ecosystem services (BES) [28].

Table 1: Key Components of Spatial Prioritization for Ecosystem Services

Component	Description	Measurement Approaches	Data Requirements
ES Supply	Biophysical capacity of ecosystems to provide services	Biophysical modeling, remote sensing, field surveys	Land cover maps, species distributions, ecosystem process data
ES Demand	Level of need or use of services by beneficiaries	Population data, stakeholder surveys, consumption patterns	Demographic data, visitor counts, household surveys
Alternative Availability	Presence of human-engineered substitutes for ES	Infrastructure inventories, market analyses, technological assessments	Utility service areas, market prices, replacement cost data
Spatial Flow	Movement of services from supply to demand areas	Connectivity modeling, distance decay functions, network analysis	Landscape connectivity, hydrological flow, transportation networks

Data Requirements and Assessment Methods

Foundational Spatial Data

Spatial prioritization requires the integration of diverse datasets representing both ecological and social systems. Essential data layers include:

• Land Use/Land Cover (LULC): High-resolution LULC data forms the foundation for assessing ecosystem service supply [29] [30]. The China 30 m resolution land use dataset exemplifies the granular data needed for accurate ES assessment [30].
• Species Distributions: Data on target species distributions, particularly threatened, endemic, or keystone species, is crucial for integrating biodiversity conservation with ES management [31]. Systematic conservation planning approaches identify 200 or more target species in comprehensive assessments [31].
• Socio-economic Data: Population density, income levels, infrastructure development, and market accessibility help quantify demand for ES and opportunity costs of conservation [32].
• Biophysical Parameters: Climate data, soil characteristics, topographic information, and hydrological patterns enable modeling of ecosystem processes underpinning service provision [31] [30].

Ecosystem Service Assessment Methods

Multiple methodological approaches exist for quantifying and mapping ecosystem services:

• Unit Value Transfer Methods: The value equivalent factor method, pioneered by Costanza et al. and refined by Xie et al., applies standardized ES values to land cover classes [29] [30]. This approach benefits from straightforward calculation and minimal data requirements but can overlook regional variation in service provision [29].
• Biophysical Modeling: Process-based models like InVEST (Integrated Valuation of Ecosystem Services and Tradeoffs) simulate the biophysical processes underlying service provision, such as water yield, sediment retention, and carbon storage [31]. These models provide more spatially explicit results but require extensive parameterization [31].
• Social Valuation Methods: The Social Values for Ecosystem Services (SolVES) model integrates survey data on perceived social values with environmental variables to map non-material ES benefits [33]. This approach effectively captures cultural services like aesthetics, recreation, and spiritual values [33].

Table 2: Ecosystem Service Assessment Methods for Spatial Prioritization

Method Category	Specific Methods	Best Applications	Limitations
Value Transfer	Equivalent factor method, value coefficients	Regional screening assessments, comparative analyses	Insensitive to local context, assumes uniform value
Biophysical Modeling	InVEST, ARIES, SWAT	Watershed management, climate regulation services	Data intensive, complex parameterization
Social Preference	SolVES, participatory mapping, willingness-to-pay	Cultural services, urban planning, recreation management	Subject to response biases, limited spatial extrapolation
Habitat Suitability	Species Distribution Models (MaxEnt), habitat connectivity	Biodiversity-focused conservation, corridor design	Species-specific, limited to known habitat associations

Integrated Methodological Workflow

Systematic Conservation Planning Approach

The following workflow outlines a comprehensive methodology for spatial prioritization of ecosystem services:

Figure 1: Spatial Prioritization Workflow for ES Management

Detailed Experimental Protocols

Protocol 1: Priority Area Identification Using Systematic Conservation Planning

Purpose: To identify spatially explicit priority areas for ecosystem service management that balance species conservation and ES provision.

Materials and Software:

• GIS software (ArcGIS, QGIS)
• Systematic conservation planning tools (C-Plan, Zonation, Marxan)
• Species distribution modeling software (MaxEnt)
• Ecosystem service modeling tools (InVEST)
• Land use/land cover datasets
• Species occurrence records
• Socio-economic data

Procedure:

• Define Conservation Targets: Identify 150-200 target species based on IUCN Red List status, national protection level, endemism, economic value, and keystone status in ecosystems [31].
• Model Species Distributions: Use MaxEnt or similar species distribution models to predict potential habitat ranges for target species based on occurrence records and environmental variables [31].
• Map Ecosystem Services: Employ InVEST models to quantify and map key ecosystem services (carbon storage, water yield, sediment retention) [31].
• Calculate Irreplaceability: Use C-Plan or similar systematic conservation planning software to calculate irreplaceability (Ir) values for each planning unit based on species characteristics and conservation targets [31].
• Normalize and Integrate Datasets: Normalize species irreplaceability values and ecosystem service values to a common scale (0-1) using min-max normalization [31].
• Identify Priority Areas: Overlay normalized datasets using equal weighting or stakeholder-defined weights to derive comprehensive priority conservation areas [31].
• Conduct Gap Analysis: Compare identified priority areas with existing protected area networks to identify conservation gaps [31].

Analysis: Priority conservation areas typically cover 7-10% of the total study area while protecting a disproportionate share of biodiversity and ecosystem service values [31].

Protocol 2: Social Values Assessment for Cultural Services

Purpose: To quantify and map social values for ecosystem services, particularly cultural services, to inform spatial planning decisions.

Materials and Software:

• Social Values for Ecosystem Services (SolVES) model
• Survey tools for data collection
• GIS software with spatial analysis capabilities
• Environmental variable datasets (elevation, slope, distance to water, land cover)

Procedure:

• Survey Design: Develop a survey instrument that allows respondents to allocate points or virtual currency among different social value types (aesthetic, biodiversity, cultural, recreational, spiritual, therapeutic) [33].
• Data Collection: Administer surveys to representative stakeholders, collecting demographic information and spatial preferences for ecosystem service values [33].
• Spatial Explicit Response: Ask respondents to identify specific locations they associate with different social values using participatory mapping approaches [33].
• Environmental Variable Preparation: Compile relevant environmental variables (elevation, slope, distance to water, land cover types) that may influence social value distributions [33].
• Model Run: Execute SolVES model to generate value maps and quantitative metrics for each social value type [33].
• Hotspot Analysis: Conduct spatial autocorrelation analysis to identify significant clusters (hotspots and coldspots) of social values [33].
• Relationship Analysis: Examine correlations between different social value types and between social values and environmental variables [33].

Analysis: Research indicates that aesthetic values typically cover the largest spatial extent among social values, while spiritual and therapeutic values exhibit more limited distributions [33]. Significant spatial clustering patterns emerge, with respondents who prioritize aesthetic values also tending to appreciate biodiversity, recreational, spiritual, and therapeutic values [33].

The Scientist's Toolkit: Research Reagent Solutions

Table 3: Essential Tools and Data Sources for Spatial Prioritization Research

Tool/Data Category	Specific Solutions	Function	Access Information
Spatial Planning Software	C-Plan, Marxan, Zonation	Identify priority areas using systematic algorithms	Open source or commercial licensing
Species Distribution Modeling	MaxEnt, BIOMOD	Predict potential species habitats from occurrence data	Open source with Java requirement
Ecosystem Service Modeling	InVEST, ARIES, SolVES	Quantify and map ecosystem service supply and value	Open source, various platforms
Social Survey Platforms	Qualtrics, Survey123, KoboToolbox	Collect and manage social preference data	Commercial and open source options
Land Use/Land Cover Data	CNLUCC, CORINE, MODIS	Foundation for ES assessment and change analysis	Various spatial and temporal resolutions
Opportunity Cost Data	EULOCC1K, national statistics	Estimate economic trade-offs of conservation actions	European layer available at 1km resolution [32]

Implementation Considerations and Analytical Framework

Addressing Spatial and Temporal Dynamics

Effective spatial prioritization must account for the dynamic nature of socio-ecological systems. Key considerations include:

• Land Use Change Analysis: Utilize land use dynamics indices and transition matrices to understand historical patterns and project future scenarios [29]. The single dynamic attitude (Kt) and comprehensive dynamic attitude (St) metrics quantify rates of land use change over specific periods [29].
• ESV Trajectory Assessment: Monitor ecosystem service value (ESV) fluctuations over time, typically using 20-year periods with 5-year intervals to capture trends [29] [30]. In Xi'an, China, ESV increased by 938.8 million yuan from 2000-2020, with high-value areas concentrated in forested regions and river systems [29].
• Fragmentation Impacts: Assess how ecosystem fragmentation affects service delivery, as perimeter-area ratios and patch proximity correlate with declines in ES provision over decades [25].

Priority Setting and Gap Analysis

Figure 2: Conservation Priority and Gap Analysis Framework

Integration with Policy and Planning

Spatial prioritization outcomes must ultimately inform decision-making processes:

• Policy Alignment: Assess how identified priorities align with existing environmental policies and protected area networks [25] [31]. In Ethiopia, research revealed limited integration of ES concepts into agricultural policies, which prioritized provisioning services over regulating and cultural services [25].
• Stakeholder Engagement: Incorporate participatory approaches throughout the prioritization process to enhance legitimacy and implementation success [26] [24].
• Ecological Compensation Mechanisms: Develop compensation priority scores (e.g., ECPS - Ecological Compensation Priority Score) based on the ratio of non-market ESV to GDP per unit area to guide fiscal transfer mechanisms [30]. In Xizang, the northwestern Qiangtang Plateau desert ecological zone exhibited the highest compensation priority, with a theoretical compensation amount of approximately 1.6 trillion CNY in 2020 [30].

This comprehensive framework for spatial prioritization in ecosystem service management provides researchers and practitioners with a structured approach to balance biodiversity conservation with human well-being across diverse landscapes. The integration of ecological data, social values, economic considerations, and policy contexts enables more effective and sustainable resource management decisions.

The identification of biodiversity hotspots is a foundational step in conservation planning, enabling the prioritization of areas with exceptional ecological value. Within the context of analysing spatial congruence between biodiversity and ecosystem services, precise hotspot identification is critical. Techniques range from straightforward metric-based rankings to complex spatial clustering algorithms, each with distinct advantages and applications. These methods allow researchers to transform raw spatial data on species distributions, ecosystem service flows, or phylogenetic endemism into actionable maps that guide resource allocation and policy decisions. The choice of technique directly influences which areas are designated for conservation action, making a thorough understanding of methodological nuances essential for robust spatial congruence analysis.

Core Hotspot Identification Techniques

Top Richest Cells (Quantile Method)

The Top Richest Cells method, also known as the quantile approach, is one of the most direct techniques for hotspot identification. It involves ranking all grid cells in a study area from high to low based on the value of a target biodiversity or ecosystem service variable (e.g., species richness, carbon storage capacity). Hotspots are then defined as the top percentile of cells, such as the richest 5%, 10%, or 20% [16]. This method is particularly useful for initial, broad-scale assessments.

• Experimental Protocol:
  • Data Preparation: Compile a georeferenced presence-absence matrix (PAM) for the taxonomic group or a spatial layer for the ecosystem service of interest at a consistent resolution (e.g., 1°×1° or 100 km x 100 km grid) [34].
  • Metric Calculation: For each cell, calculate the target metric. For biodiversity, this is often species richness (the count of species per cell). For a more nuanced view, Weighted Endemicity (WE) can be calculated, which weights species by the inverse of their range size, preventing widespread species from dominating the analysis [35].
  • Ranking and Delineation: Rank all cells based on the calculated metric. Apply a predetermined quantile threshold (e.g., the top 20% of cells) to classify hotspots [34] [35].
  • Validation: The robustness of the identified hotspots can be tested by examining their sensitivity to different quantile thresholds and comparing them to known protected areas to assess current representation [34].

Spatial Clustering Methods

Spatial clustering methods identify hotspots based on statistically significant spatial autocorrelation, where similar values are not randomly distributed but cluster together on the landscape. Unlike the quantile method, these techniques explicitly incorporate the geographic location of cells in the analysis.

• Getis-Ord Gi*: This statistic identifies spatial clusters of high values ("hot spots") and low values ("cold spots") [16] [36]. It produces a Z-score for each feature, indicating whether the observed spatial clustering of high or low values is more pronounced than one would expect by random chance.

• Anselin Local Moran's I: This Local Indicator of Spatial Association (LISA) classifies features into five categories: high-high clusters (hotspots), low-low clusters (coldspots), high-low outliers, low-high outliers, and statistically non-significant features [36].

• Experimental Protocol for Getis-Ord Gi*:

  • Input Data: Use a grid layer with a calculated biodiversity or ecosystem service value for each cell.
  • Define Conceptualization of Spatial Relationships: Choose a method to define how cells influence each other (e.g., "Inverse Distance" or "Zone of Indifference").
  • Define Distance Threshold: Set a critical distance or a fixed number of neighbours beyond which spatial influence is zero.
  • Run Analysis: Calculate the Gi* statistic for each cell in the study area. The resulting Z-score and p-value indicate the significance of clustering.
  • Interpretation: Cells with high values and a high Z-score (significant positive spatial autocorrelation) are identified as hotspots. As demonstrated in a study of the European rabbit on Lemnos Island, this can effectively reveal core areas of high density [36].

Threshold and Intensity Methods

The Threshold Method defines hotspots based on an expert-defined biophysical value, independent of the overall distribution of the data. For example, a hotspot for soil accumulation might be defined as areas with a soil depth ≥0.8 m and ≥70% litter cover [16]. This method is highly objective but requires robust, pre-existing scientific justification for the threshold.

The Intensity Method is similar to the quantile approach but focuses on the absolute provision of a service or level of diversity. It defines hotspots as "areas which provide large proportions of a particular service" [16]. For instance, a hotspot could be defined as the set of cells that collectively contribute to the top 10% of the total service provision in the region.

Richness-Based Methods for Multiple Services

When moving from single to multiple ecosystem services or biodiversity facets, the Richness Method is commonly applied. It defines hotspots as areas where multiple services or diversity metrics overlap [16]. In this approach, a hotspot for three services would be a cell that is a hotspot for Service A and Service B and Service C. This method directly supports congruence analysis by pinpointing areas of co-occurrence, which are high priorities for integrated conservation strategies.

Table 1: Comparative Summary of Core Hotspot Identification Techniques

Technique	Core Principle	Key Advantages	Key Limitations	Best Use-Case
Top Richest Cells	Ranks cells by a single metric and selects the top quantile.	Simple, intuitive, and easy to implement.	Ignores spatial configuration; sensitive to scale and tied values.	Initial, broad-scale assessments and single-metric prioritization.
Spatial Clustering	Identifies statistically significant clusters of high values.	Accounts for spatial autocorrelation; produces statistically robust results.	More computationally intensive; results sensitive to neighbourhood definition.	Identifying core areas of density and understanding spatial structure.
Threshold/Intensity	Uses absolute biophysical or provision-level thresholds.	Objective and directly linked to ecological function.	Requires robust scientific basis for thresholds; may not capture relative importance.	When clear, scientifically-defensible targets exist (e.g., minimum habitat area).
Richness (Overlap)	Identifies areas of overlap for multiple service/diversity hotspots.	Directly targets spatial congruence; maximizes multi-functionality.	Can overlook areas important for a single, critical service.	Multi-objective conservation planning and ecosystem service bundle analysis.

Advanced Considerations in Hotspot Analysis

Integrating Multiple Dimensions of Biodiversity

Effective congruence analysis requires looking beyond taxonomic species richness. Phylogenetic Diversity (PD), which captures evolutionary history, and Functional Diversity (FD), which reflects the range of ecological traits, are critical dimensions that often do not align perfectly with species richness [37]. For example, a study on simulated plant communities found that hotspots for taxonomic richness, functional richness, and phylogenetic diversity showed only moderate overlap, demonstrating that one dimension cannot serve as a proxy for another [37].

Furthermore, Traditional Ecological Knowledge (TEK) offers a vital socio-ecological dimension. Metrics such as the diversity of Indigenous plant names and the number of documented traditional uses for species can identify areas of high biocultural value. Research shows that TEK-based hotspots sometimes, but not always, overlap with those identified by Western scientific metrics, highlighting the risk of overlooking culturally significant areas if only standard metrics are used [37].

Macroevolutionary Routes to Hotspot Formation

Understanding the evolutionary processes that generate hotspots can inform their conservation. Global analyses of mammals and birds reveal distinct macroevolutionary routes:

• Tropical Realms: Hotspots (e.g., in the Neotropics and Indo-Malay) are primarily shaped by high rates of in situ speciation over the past 25 million years [35].
• Temperate Realms: Hotspots (e.g., in the Palearctic and Nearctic) accumulated diversity through high rates of immigration from surrounding regions, rather than exceptionally high speciation [35].

This fundamental difference implies that conservation strategies may need to vary: protecting the complex landscapes that drive speciation in the tropics, versus maintaining connectivity corridors that facilitate dispersal in temperate zones.

Methodological Impact on Conservation Outcomes

The choice of hotspot identification method has direct, practical consequences for conservation planning. A comparative study found that the spatial configuration of priority areas selected by four different hotspot methods (top richest cells, spatial clustering, intensity, and richness) showed very little overlap (Kappa statistics of 0.11–0.27) [16]. Furthermore, these methods differed significantly in their cost-effectiveness and compactness when compared to systematic conservation planning tools like Marxan. This underscores that the term "hotspot" is not a universal standard but a context-dependent construct, and the method should be carefully selected to match specific conservation goals.

The following workflow diagram illustrates the decision-making process for selecting and applying these techniques within a broader research framework.

Diagram 1: Hotspot Identification and Analysis Workflow

Table 2: Key Research Reagent Solutions for Hotspot Identification

Item/Category	Function & Application in Hotspot Analysis
Presence-Absence Matrix (PAM)	The foundational data structure; a binary matrix (grid cells x species) representing species distributions, from which richness and rarity metrics are derived [34].
Spatial Analysis Software	Platforms like R (with packages such as sp, sf, raster, spdep), ArcGIS (ESRI Inc.), and QGIS are essential for spatial data handling, metric calculation, and running clustering algorithms [34] [36].
Global Biodiversity Datasets	Data sources like the Global Amphibian Assessment, GBIF, and IUCN Red List spatial data provide the foundational species range maps required for analysis at broad scales [34] [35].
Phylogenetic Trees	Comprehensive, time-calibrated phylogenies for the taxonomic group of interest (e.g., mammals, birds) are necessary for calculating phylogenetic diversity metrics and diversification rates [35].
Systematic Conservation Planning Software	Tools like Marxan and Zonation provide complementary, target-based approaches for prioritization, against which hotspot method outcomes can be compared [38] [16].
Environmental Predictor Layers	Spatial data on net primary productivity (NPP), topography (e.g., Terrain Ruggedness Index), climate, and land use are used to model distributions and explain the drivers of hotspot patterns [35] [39].

Multi-Criteria Decision-Making and Scenario Analysis

Multi-criteria decision-making (MCDM) provides a structured analytical approach for evaluating complex scenarios involving multiple conflicting objectives and criteria. In the context of biodiversity and ecosystem service research, MCDM methods enable researchers and policymakers to systematically analyze spatial congruence and trade-offs across diverse ecological values. The integration of scenario analysis with MCDM creates a powerful framework for exploring potential futures under different conservation and development pathways, supporting more robust and defensible environmental decisions.

This protocol details the application of MCDM and scenario analysis specifically for spatial congruence analysis between biodiversity patterns and ecosystem service distributions. The structured approach enables researchers to identify priority areas that maximize conservation outcomes while considering complex socio-ecological trade-offs. With accelerating biodiversity loss and increasing pressure on ecosystem services, these methodologies provide critical support for spatial planning and resource allocation decisions.

Experimental Protocols

Spatial Data Collection and Preparation Protocol

Objective: To collect, process, and harmonize spatial data on biodiversity, ecosystem services, and environmental variables for subsequent MCDM analysis.

Materials and Equipment:

• Geographic Information System (GIS) software (e.g., QGIS, ArcGIS)
• Remote sensing data platforms (e.g., Copernicus, Landsat)
• Spatial data on biodiversity occurrences (GBIF, OBIS)
• Environmental variables datasets (WorldClim, SoilGrids)
• High-performance computing resources for spatial analysis

Procedure:

• Data Identification and Acquisition

  • Compile biodiversity data from global and regional databases (GBIF for terrestrial species, OBIS for marine species) [40]
  • Collect ecosystem service data using standardized models:
    • Water yield: Apply Budyko curve methodology with annual precipitation data [41]
    • Carbon sequestration: Calculate using net primary productivity (NPP) from CASA model [41]
    • Biodiversity: Assess habitat quality using InVEST habitat quality model [41]
    • Cultural services: Quantify using SolVES 3.0 model incorporating landscape metrics and survey data [41]
  • Gather environmental drivers: land use/land cover, climate, topography, soil, and socio-economic data

• Data Preprocessing and Harmonization

  • Reproject all spatial layers to a common coordinate system and resolution
  • Address missing data through appropriate imputation methods
  • Normalize data distributions where necessary
  • Rescale data to consistent spatial units (e.g., 1km² grid cells)

• Spatial Congruence Analysis

  • Calculate correlation coefficients between biodiversity and ecosystem service layers
  • Identify areas of high and low spatial congruence using hotspot analysis
  • Document all processing steps in reproducible workflows

Table 1: Essential Spatial Data Types for Biodiversity-Ecosystem Service Congruence Analysis

Data Category	Specific Variables	Source Examples	Spatial Resolution	Temporal Resolution
Biodiversity	Species occurrence records, habitat quality, species richness	GBIF, OBIS, IUCN Red List, regional monitoring networks [23] [40]	Varies (point locations to 1km grid)	Annual to decadal
Ecosystem Services	Water yield, carbon sequestration, pollination, recreation	InVEST, ARIES, CoSting Nature models [41]	30m - 1km	Annual
Environmental Drivers	Land use/land cover, climate, topography, soil characteristics	Copernicus, MODIS, WorldClim, SoilGrids	10m - 1km	Annual to decadal
Socio-economic	Population density, economic activity, infrastructure	National statistics, census data, night-time lights	100m - 1km	Annual

Ordered Weighted Averaging for Scenario Analysis

Objective: To implement Ordered Weighted Averaging (OWA) for evaluating ecosystem service trade-offs under different development-conservation scenarios.

Theoretical Background: OWA is a multi-criteria decision-making approach that enables flexible aggregation of criteria based on decision-maker risk preferences. The method orders criterion values before applying weights, allowing control over the trade-off level between criteria [41].

Procedure:

• Scenario Definition

  • Define a range of scenarios representing different conservation and development priorities (e.g., protection-focused, balanced, development-focused)
  • Establish clear weighting schemes for each scenario reflecting different policy preferences

• Criteria Standardization

  • Transform all ecosystem service maps to a common scale (0-1) using linear scaling or value functions
  • Ensure directionality is consistent (higher values always preferred)

• OWA Implementation

  • For each grid cell, rank ecosystem service values from highest to lowest
  • Apply ordered weights according to the decision attitude (risk-averse, risk-neutral, risk-seeking)
  • Calculate weighted overall score for each grid cell using the formula:
    where ωj is the jth weight and bj is the jth largest criterion value [41]

• Hotspot and Coldspot Identification

  • Apply spatial clustering analysis (e.g., Getis-Ord Gi*) to identify statistically significant hotspots and coldspots
  • Calculate area statistics for each hotspot category

• Sensitivity Analysis

  • Test robustness of results to weight variations
  • Identify scenarios where priority areas remain stable versus those where they shift significantly

MCDM Framework Implementation for Spatial Conservation Prioritization

Objective: To implement a comprehensive MCDM framework for identifying priority areas that balance biodiversity conservation and ecosystem service provision.

Materials and Equipment:

• R or Python with MCDM packages (scikit-criteria, MCDA)
• GIS software with spatial analysis capabilities
• High-performance computing for large raster operations

Procedure:

• Structuring the Decision Problem

  • Define clear objectives for spatial prioritization
  • Identify relevant criteria (biodiversity and ecosystem service indicators)
  • Establish constraints (budget, minimum protection targets)

• Weight Elicitation

  • Conduct expert surveys or stakeholder workshops to assign criterion weights
  • Consider using entropy weighting for objective weight determination [42]
  • Document rationale for weight assignments

• Alternative Evaluation

  • Generate and evaluate different spatial configurations of protected areas
  • Apply MCDM methods such as:
    • VIKOR: Balancing group utility and individual regret [42]
    • AHP: Hierarchical pairwise comparisons [43]
    • TOPSIS: Proximity to ideal solution

• Result Validation

  • Compare MCDM results with existing protected areas network
  • Conduct gap analysis to identify underrepresented ecosystems
  • Validate with field data where available

Table 2: MCDM Methods Comparison for Spatial Conservation Planning

Method	Key Features	Strengths	Limitations	Best Application Context
Ordered Weighted Averaging (OWA)	Controls trade-off level through ordered weights	Flexible attitude to risk; Explicit scenario exploration	Requires understanding of decision attitude	Scenario-based planning with clear risk preferences [41]
VIKOR	Compromise ranking minimizing individual regret	Balances majority and minority concerns; Good for conflicting criteria	Complex implementation; Requires precise performance values	Controversial decisions with stakeholder conflicts [42]
Analytic Hierarchy Process (AHP)	Pairwise comparisons in hierarchical structure	Handles qualitative judgments; Consistency checking	Subject to ranking reversal; Many pairwise comparisons	Stakeholder-driven planning with mixed criteria types [43]
TOPSIS	Distance-based approach to ideal solution	Intuitive logic; Simple computation	Sensitive to weight assignment; Normalization issues	Rapid screening of multiple spatial alternatives

Data Analysis and Interpretation

Spatial Congruence Metrics Calculation

Objective: To quantify and interpret spatial relationships between biodiversity and ecosystem services.

Procedure:

• Spatial Correlation Analysis

  • Calculate Pearson's correlation coefficients between standardized biodiversity and ecosystem service layers
  • Generate correlation matrices to visualize relationships
  • Test statistical significance of correlations using spatial autocorrelation-adjusted tests

• Spatial Overlap Quantification

  • Identify hotspots for each criterion (top 20% of values)
  • Calculate overlap coefficients:
    • Jaccard index for hotspot overlap
    • Percentage of biodiversity hotspots that coincide with ecosystem service hotspots
  • Map areas of high and low congruence

• Trade-off Analysis

  • Identify areas where high biodiversity coincides with low ecosystem services and vice versa
  • Quantify the extent and distribution of these trade-off areas
  • Analyze environmental and socio-economic correlates of trade-off patterns

Scenario Comparison and Robustness Assessment

Objective: To compare outcomes across different scenarios and identify robust priority areas.

Procedure:

• Scenario Outcomes Comparison

  • Calculate total area and spatial configuration of priority areas for each scenario
  • Compare ecosystem service provision and biodiversity representation across scenarios
  • Identify "no-regret" areas that appear as priorities across multiple scenarios

• Efficiency Analysis

  • Calculate protection efficiency metrics for different areas [41]
  • Identify areas with high ecosystem service values relative to opportunity costs
  • Generate efficiency frontiers for different conservation budgets

• Uncertainty Propagation

  • Assess sensitivity of results to data uncertainty and weight assignments
  • Use Monte Carlo simulation to propagate uncertainty through the MCDM framework
  • Identify priority areas that remain stable under uncertainty

The Scientist's Toolkit: Research Reagent Solutions

Table 3: Essential Tools for MCDM in Biodiversity and Ecosystem Service Research

Tool Category	Specific Tools/Platforms	Function	Application Notes
Spatial Data Platforms	GBIF, OBIS, Restor, IUCN Contributions for Nature [44] [40]	Biodiversity occurrence data and project tracking	Essential for baseline data; Varying spatial and taxonomic resolution
Ecosystem Service Models	InVEST, ARIES, SolVES [41]	Quantification and mapping of ecosystem services	Model selection depends on study scale and data availability
MCDM Software	R (scikit-criteria, MCDA), Python (PyDecision), Decerns	Implementing MCDM algorithms	Open-source options available; Varying learning curves
Spatial Analysis	QGIS, ArcGIS, GRASS GIS, Whitebox Tools	Spatial data processing and analysis	Critical for data preparation and result mapping
Biodiversity Monitoring Tech	eDNA kits, acoustic sensors, remote sensing platforms [45]	Novel data collection for biodiversity	Increasingly important for filling data gaps
Decision Support Systems	Marxan, Zonation, PrioritizR	Systematic conservation planning	Specialized tools for spatial prioritization

Implementation Considerations

Data Quality and Harmonization

Successful application of MCDM for spatial congruence analysis requires careful attention to data quality and harmonization. Key considerations include:

• Scale Matching: Ensure biodiversity and ecosystem service data are at compatible spatial and temporal scales
• Uncertainty Propagation: Account for uncertainty in underlying data through sensitivity analysis
• Metadata Documentation: Maintain comprehensive metadata for all datasets to support reproducibility

Stakeholder Engagement

MCDM processes benefit significantly from appropriate stakeholder engagement:

• Weight Elicitation: Engage diverse stakeholders in criterion weight assignment to reflect multiple values
• Scenario Co-development: Collaborate with stakeholders to develop meaningful and relevant scenarios
• Result Interpretation: Facilitate discussions around trade-offs and implementation constraints

Technical Implementation

• Computational Efficiency: For large spatial datasets, employ efficient computing strategies (parallel processing, cloud computing)
• Reproducibility: Implement version control and containerization for analytical workflows
• Visualization: Develop effective visualization strategies to communicate complex MCDM results to diverse audiences

The protocols outlined provide a comprehensive framework for applying MCDM and scenario analysis to spatial congruence analysis of biodiversity and ecosystem services. By following these standardized approaches, researchers can generate comparable, transparent, and defensible results to inform conservation planning and natural resource management decisions.

Accurately quantifying ecosystem services (ES) is imperative for sustainable ecosystem management and informed policy development [46]. As the field advances, sophisticated spatial models have been developed to map and value the goods and services from nature that sustain and fulfill human life [47]. This protocol provides detailed application notes for three advanced modeling approaches—InVEST, SolVES, and TEM—that enable researchers to quantify ecosystem services with high spatial explicitness. These tools are particularly valuable for assessing biodiversity and ecosystem service spatial congruence, a core challenge in ecological research [48]. Each model offers distinct capabilities: InVEST focuses on biophysical and economic quantification of multiple ES; SolVES specializes in assessing social perceptions and cultural values; while TEM provides process-based understanding of biogeochemical cycles. When integrated within a comprehensive research framework, these models facilitate a holistic understanding of ecosystem service supply, distribution, and societal valuation—critical knowledge for balancing conservation and development objectives in landscape planning and management.

Model Comparison and Selection Framework

Table 1: Comparative overview of ecosystem service quantification models

Model	Primary Developer	Core Focus	Spatial Explicitness	Key Outputs	Temporal Dynamics	Data Requirements
InVEST	Stanford Natural Capital Project	Biophysical & economic valuation of multiple ES	High - raster-based mapping	Maps of ES provision in biophysical or economic terms [47]	Multi-year snapshots; scenario projection [48]	Land cover, DEM, biodiversity, climate, soil data
SolVES	U.S. Geological Survey	Social values & cultural services	High - integrates survey data with environmental variables	Social Value Maps, Value Indexes [33]	Cross-sectional perception analysis	Survey data, participatory mapping, environmental variables
TEM	University of New Hampshire	Biogeochemical cycles & ecosystem processes	Variable - grid-based	Carbon, nitrogen, water fluxes	Continuous simulation; climate change projections	Climate, vegetation, soil, elevation data

Table 2: Ecosystem service types addressed by each model

Ecosystem Service Category	Specific Service Types	InVEST	SolVES	TEM
Provisioning	Food production, Water yield, Timber	✓	Limited	✓
Regulating	Carbon storage, Soil conservation, Water purification, Climate regulation	✓ [49] [48]	Limited	✓
Cultural	Aesthetic value, Recreation, Spiritual values	Limited	✓ [33]	Limited
Supporting	Habitat quality, Nutrient cycling, Biodiversity	✓ [48]	✓ [33]	✓

Model-Specific Application Protocols

InVEST (Integrated Valuation of Ecosystem Services and Tradeoffs) Protocol

InVEST is a suite of free, open-source software models that use production functions to quantify how changes in ecosystem structure and function affect the flows and values of ecosystem services across landscapes [47]. The models are spatially explicit, using maps as information sources and producing maps as outputs, with results provided in either biophysical or economic terms [47]. The modular design allows users to select only those ecosystem services of interest rather than modeling all services [47].

Implementation Workflow

Detailed Methodology for Key Modules

Habitat Quality Assessment

• Data Requirements: Land use/land cover (LULC) maps, threat data (location and weight), habitat sensitivity of LULC types to each threat [48]
• Processing Steps:
  • Classify LULC types and assign habitat suitability scores (0-1)
  • Identify threat sources and assign weights, maximum effective distance, and decay type
  • Calculate relative impact of threats across the landscape
  • Compute habitat quality score using the equation: [ Q{xj} = Hj \left(1 - \frac{D{xj}^z}{D{xj}^z + k^z}\right) ] where (Q{xj}) is habitat quality of pixel x in LULC j, (Hj) is habitat suitability, (D_{xj}) is total threat level, k is half-saturation constant, and z is scaling parameter
• Output Interpretation: Values range from 0 (low quality) to 1 (high quality); useful for identifying biodiversity hotspots and conservation priorities [48]

Carbon Storage Assessment

• Data Requirements: LULC maps, carbon pool estimates (aboveground, belowground, soil, dead organic matter) for each LULC type
• Processing Steps:
  • Compile or estimate four carbon pools for each LULC class
  • Sum pools to calculate total carbon storage for each grid cell: [ C{total} = C{above} + C{below} + C{soil} + C_{dead} ]
  • For scenario analysis, model changes in carbon storage based on LULC transitions
• Output Interpretation: Maps display total carbon storage (Mg/ha); valuable for climate regulation assessments and payment for ecosystem services schemes [48]

SolVES (Social Values for Ecosystem Services) Protocol

SolVES is designed to assess, quantify, and map the perceived social values of ecosystem services, with particular emphasis on cultural services [33]. The model integrates public perception surveys with environmental data to analyze the spatial distribution of social values and investigate factors influencing them [33]. SolVES operates on the principle that social values demonstrate predictable relationships with environmental features that can be modeled and mapped [33].

Implementation Workflow

Detailed Methodology for Social Value Assessment

Survey Design and Value Allocation

• Survey Implementation:
  • Administer surveys to diverse stakeholder groups (residents, visitors, managers)
  • Utilize value allocation method: provide respondents with a virtual currency (e.g., 100 points) to allocate among different social value types according to perceived importance [33]
  • Incorporate participatory mapping: respondents mark locations of valued areas on maps or GIS interfaces
• Social Value Types (based on Dalian City case study [33]):
  • Aesthetic, biodiversity, cultural, recreational, spiritual, therapeutic, educational, ecological sustainability values
  • Document respondent preferences; Dalian study showed pronounced preference for aesthetic, cultural, and biodiversity values [33]

Data Integration and Analysis

• Environmental Variable Preparation:
  • Calculate Euclidean distance to water features, trails, roads, urban areas
  • Derive topographic variables (elevation, slope) from DEM
  • Prepare land cover classifications
• Model Execution:
  • SolVES calculates Value Index for each social value type: [ VI = \frac{\text{Number of survey points valuing that type}}{\text{Total survey points}} \times \text{Average points allocated} ]
  • Generate social value maps showing spatial distribution and intensity of each value type
  • Conduct hotspot analysis using Getis-Ord Gi* statistic to identify spatial clustering [33]
• Output Interpretation: Maps reveal areas of high social value concentration; statistical analysis identifies relationships between social values and environmental features

TEM (Terrestrial Ecosystem Model) Protocol

TEM is a process-based biogeochemistry model that simulates the cycling of carbon, nitrogen, and water through terrestrial ecosystems. While less explicitly focused on ecosystem services than InVEST or SolVES, TEM provides critical understanding of fundamental ecological processes that underpin many regulating and supporting services, particularly those related to climate regulation and nutrient cycling.

Implementation Workflow

Detailed Methodology for Ecosystem Process Modeling

Data Preparation and Model Parameterization

• Input Data Requirements:
  • Climate data: temperature, precipitation, solar radiation, humidity (daily/monthly)
  • Vegetation parameters: leaf area index, carbon allocation patterns, nutrient content
  • Soil properties: texture, pH, organic matter content, bulk density
  • Land cover classification
• Model Parameterization:
  • Define ecosystem types based on vegetation and climate characteristics
  • Set biogeochemical constants for decomposition, mineralization, photosynthesis processes
  • Calibrate using field measurements of NPP, soil carbon, and flux tower data

Simulation and Output Analysis

• Core Calculations:
  • Net Primary Production (NPP) based on temperature, light, water, and nutrient availability
  • Heterotrophic respiration constrained by temperature, moisture, and substrate quality
  • Net Ecosystem Production (NEP) = NPP - Heterotrophic Respiration
  • Nitrogen mineralization and immobilization processes
  • Water balance including evapotranspiration, runoff, and soil moisture
• Output Interpretation:
  • Maps of carbon fluxes (NPP, NEP) and stocks
  • Assessments of climate change impacts on ecosystem processes
  • Projections of future ecosystem conditions under different climate scenarios

Integrated Application for Biodiversity and Ecosystem Service Congruence

Multi-Model Integration Framework

Table 3: Strategic model integration for congruence analysis

Research Objective	Primary Model	Secondary Model	Integration Approach	Congruence Metrics
Cultural service hotspots vs. biodiversity	SolVES (social values)	InVEST (habitat quality)	Spatial overlay analysis	Co-location percentage, Correlation coefficients
Biogeochemical underpinnings of ecosystem services	TEM (processes)	InVEST (services)	Output correlation analysis	Regression slopes, Variance explanation
Scenario analysis for conservation planning	InVEST (multiple services)	SolVES (social preferences)	Weighted multi-criteria evaluation	Trade-off curves, Spatial congruence indices

Case Study Application: Yunnan-Guizhou Plateau

A comprehensive study on the Yunnan-Guizhou Plateau demonstrates the integrated application of ecosystem service models [48]. Researchers employed machine learning with the PLUS model to project land use changes, then used InVEST to evaluate future ecosystem services under three scenarios (natural development, planning-oriented, and ecological priority) for 2035 [48]. This approach allowed quantification of four key services—water yield, carbon storage, habitat quality, and soil conservation—and revealed that the ecological priority scenario demonstrated the best performance across all services [48]. The study exemplifies how model integration can inform sustainable development strategies in ecologically sensitive regions.

Advanced Analytical Techniques

Spatial Congruence Analysis

• Hotspot Identification: Apply Getis-Ord Gi* statistic to identify statistically significant hotspots and coldspots for each ecosystem service [33]
• Overlay Analysis: Spatially intersect biodiversity priority areas with ecosystem service hotspots using GIS overlay operations
• Correlation Analysis: Calculate spatial correlation coefficients (Pearson's r) between biodiversity indices and ecosystem service values

Trade-off and Synergy Analysis

• Relationship Quantification: Use correlation analysis (Spearman's rank) to identify trade-offs (negative correlations) and synergies (positive correlations) between ecosystem services [48]
• Trade-off Curve Development: Create production possibility frontiers to visualize optimal combinations of multiple services
• Scenario Testing: Evaluate how different land-use configurations affect the relationships between biodiversity and ecosystem services

Research Reagent Solutions

Table 4: Essential research tools and data sources for ecosystem service quantification

Category	Specific Tool/Data	Specifications	Application Context	Source Examples
Spatial Data	Land Cover Maps	Minimum 30m resolution (Landsat); preferably 10m (Sentinel-2)	All models	CORINE, MODIS, ESA CCI
	Digital Elevation Model	SRTM (30m) or ALOS (12.5m)	Terrain analysis; hydrological modeling	USGS EarthExplorer
	Climate Data	Daily temperature, precipitation, solar radiation	TEM; water yield models	WorldClim, CRU
Biodiversity Data	Species Distribution Data	Point occurrences or range maps	Habitat quality assessment	GBIF, IUCN Red List
	Habitat Condition Metrics	Field-measured vegetation structure	Model validation	Field surveys
Social Data	Survey Instruments	Structured questionnaires with mapping components	SolVES applications	Custom design
	Participatory Mapping Tools	Web-based or physical maps for respondent input	Social value assessment	Maptionnaire, Google Earth
Modeling Software	InVEST	Version 3.10.0+; requires Python 3.8+	Ecosystem service assessment	Natural Capital Project
	SolVES	Version 4.0+; QGIS plugin compatible	Social value mapping	USGS
	R/Python Libraries	raster, sf, ggplot2 (R); numpy, gdal (Python)	Data processing and analysis	CRAN, PyPI

Quality Assurance and Validation Protocols

Model Validation Techniques

InVEST Validation

• Field Verification: Conduct ground-truthing of habitat quality predictions using standardized biodiversity surveys
• Comparison with Independent Data: Validate carbon storage outputs against forest inventory data; compare water yield predictions with stream gauge measurements [46]
• Sensitivity Analysis: Perform one-at-a-time sensitivity tests to identify most influential parameters

SolVES Validation

• Survey Reliability Assessment: Calculate Cronbach's alpha for survey instruments to ensure internal consistency
• Cross-validation: Use hold-out samples or k-fold cross-validation to assess prediction accuracy [33]
• Stakeholder Feedback: Present results to local experts and residents for verification and contextual interpretation

TEM Validation

• Flux Tower Comparison: Compare simulated carbon and water fluxes with eddy covariance measurements
• Soil Carbon Validation: Validate predicted soil carbon stocks against field measurements
• Remote Sensing Comparison: Compare NPP estimates with MODIS NPP products

Addressing Model Limitations

Research indicates potential mismatches between model-based assessments and stakeholder perceptions, with one study finding stakeholder ES estimates were 32.8% higher on average than model outputs [46]. This highlights the importance of:

• Transparent Communication of model assumptions and uncertainties
• Participatory Model Refinement incorporating local knowledge
• Multi-method Approaches that combine quantitative modeling with qualitative insights

The integrated application of InVEST, SolVES, and TEM provides a powerful toolkit for quantifying ecosystem services and analyzing their spatial congruence with biodiversity. By following the standardized protocols outlined in this document, researchers can generate comparable, replicable assessments across different landscapes and scales. The case studies demonstrate that this approach yields actionable insights for conservation planning, highlighting the importance of considering both biophysical and social dimensions of ecosystems. Future methodological developments should focus on enhancing temporal dynamics, improving model integration, and strengthening validation protocols through emerging technologies like environmental DNA and remote sensing.

Systematic conservation planning (SCP) provides a quantitative framework for allocating conservation resources and actions across space and time. The integration of ecosystem services (ES) into this framework represents a critical advancement, moving beyond a singular focus on biodiversity to incorporate the myriad benefits that nature provides to people. This integration allows planners to identify areas where conservation actions can simultaneously achieve biodiversity objectives and sustain the flow of vital ES, thereby broadening the justification for conservation and engaging a wider range of stakeholders [50]. Spatial conservation prioritization (SCP) software, notably Marxan, has emerged as the primary tool for operationalizing this integration, enabling decision-makers to balance multiple, often competing, objectives in complex landscapes and seascapes [51].

The central challenge lies in moving from theoretical recognition of ES importance to practical methodologies that robustly account for the unique spatial characteristics of ES in conservation plans. Unlike many biodiversity features, ES often involve complex connectivity requirements, flows between provision and demand areas, and considerations of equitable access [51]. This document provides detailed application notes and protocols for researchers and practitioners aiming to successfully integrate ES into SCP, with a specific focus on the Marxan software suite, within the broader context of biodiversity and ES spatial congruence analysis.

Theoretical Foundation: A Connectivity Typology for Ecosystem Services

A sophisticated integration of ES into SCP requires an understanding of how different ES interact with spatial patterns. A key conceptual framework classifies ES based on their specific connectivity requirements, which fundamentally influence their ideal spatial priority patterns [51].

Table 1: Connectivity Typology for Ecosystem Services in Spatial Conservation Prioritization

Connectivity Type	Description	Exemplary Ecosystem Services	Implication for Spatial Prioritization
Low Connectivity Requirements	Services that can be provided locally in small areas and/or can be transported over long distances.	Carbon sequestration, Climate regulation [51]	Spatial pattern can be governed by requirements of other features; can be input as simple supply maps.
Provision Connectivity: Aggregation & Local Area	Services that require a minimum contiguous area for provision due to ecological processes (e.g., edge effects, minimum population sizes) or administrative/logistical needs.	Recreation, Groundwater recharge, Biodiversity conservation [51]	Priority areas should be spatially aggregated; planners must define minimum patch size or use aggregation penalties.
Provision Connectivity: Regional Networks	Services that depend on large-scale spatial dynamical processes and connected networks for maintenance.	Pollination, Meta-population persistence [51]	Requires identification and prioritization of ecologically functional networks, not just individual sites.
ES Flow: Proximity of Supply & Demand	Services where the benefit is realized only when there is connectivity (e.g., proximity, hydrological flow) between areas of service provision and areas of human demand.	Pollination, Water purification, Flood regulation [51]	Prioritization must account for both the supply of the service and the location of beneficiaries.
Dispersed Supply & Equitable Access	Services that require widespread availability and accessibility across multiple administrative or social regions to ensure equitable benefits.	Recreation, Access to wild foods [51]	Priority pattern should be more dispersed to ensure access for multiple communities or jurisdictions.

This typology provides the necessary theoretical groundwork for selecting appropriate technical methods within spatial prioritization software like Marxan. The following sections translate this theory into actionable protocols.

Methodological Protocols for Spatial Data Integration

Data Acquisition and Curation

The first step involves assembling a robust spatial dataset. Key resources include:

• UN Biodiversity Lab (UNBL): Provides access to over 400 global data layers on biodiversity, climate change, and sustainable development. It is an essential source for pre-processed, authoritative global datasets which can be used for regional to national-scale analyses [52].
• ELSA Integrated Spatial Planning Tool: A tool within UNBL that combines spatial data with systematic planning tools to help identify priority areas that support goals for biodiversity, climate, and human well-being, aligned with the Kunming-Montreal Global Biodiversity Framework [52].
• NASA-Supported Projects & Resource Watch: These platforms provide extensive Earth observation data, which can be critical for modeling ES like carbon storage, water yield, and erosion control [52].

For data-deficient species or ES, techniques like "borrowing strength" can be applied, where predictive models for poorly-sampled species borrow information from closely-related or functionally-similar well-sampled species [53].

Conceptual Workflow for Integrating Ecosystem Services

The integration of ES into a Marxan-based conservation plan follows a logical sequence from goal-setting to implementation. The diagram below visualizes this workflow, highlighting the critical decision points for ES integration.

The Researcher's Toolkit: Marxan Software Suite & Reagents

Marxan is not a single tool but a family of software designed to address different types of planning challenges. Selecting the appropriate variant and supporting tools is crucial for a successful analysis.

Table 2: The Marxan Software Suite for Conservation Planning

Software / Tool	Primary Function	Application Context	Key Consideration for ES Integration
Marxan	Designs cost-efficient networks that meet biodiversity conservation targets [54] [55].	Standard protected area design; foundational SCP.	ES can be included as conservation features or used to define cost surfaces.
Marxan with Zones	Extends Marxan to assign planning units to multiple zones (e.g., strict protection, sustainable use) to meet broader objectives [55].	Complex zoning plans; balancing multiple land-uses.	Ideal for managing trade-offs between ES provision and biodiversity conservation.
Marxan with Connectivity	Incorporates ecological connectivity (e.g., larval dispersal, animal migration) into the prioritization [55].	Designing ecological networks; freshwater/marine systems.	Directly addresses "Provision Connectivity: Regional Networks" for supporting and regulating services.
Marxan Web / MaPP	Cloud-based platform for collaborative planning with streamlined data handling and visualization [55] [56].	Stakeholder-driven processes; multi-user collaboration.	Facilitates engagement with ES beneficiaries when mapping demand and supply.
CLUZ & QMarxan	QGIS plug-ins for on-screen spatial planning and simplified interface for Marxan data preparation and analysis [55].	Desktop analysis within an open-source GIS environment.	Simplifies the process of creating and managing ES data layers for Marxan.
PrioritizR	An R package for systematic conservation prioritization using integer linear programming for optimal solutions [55].	Advanced, highly customizable statistical analysis.	Offers maximum flexibility for implementing novel ES connectivity metrics and models.

Essential Research Reagents and Data Solutions

Beyond software, a successful project requires a suite of data "reagents" and parameters.

Table 3: Essential Research Reagents for Marxan Analysis

Reagent / Parameter	Function	Protocol for Ecosystem Service Integration
Conservation Features	The biodiversity elements (species, habitats) and ES to be conserved.	Represent ES as separate feature layers (e.g., carbon stock map, recreation potential map, water yield model) [51].
Feature Targets	The quantitative amount of each feature to be secured in the plan (e.g., 30% of area).	Set targets for ES based on demand, policy goals, or stakeholder input (e.g., protect 80% of current carbon stocks, secure water provision for X people) [50].
Planning Unit Cost	The cost of conserving a planning unit, which Marxan minimizes.	Use opportunity costs of forgone agriculture or timber, or direct management costs. ES benefits can be incorporated as negative costs (i.e., benefits) [51].
Boundary Length Modifier (BLM)	A parameter that controls the compactness and aggregation of the solution.	Increase BLM to create more aggregated reserves, which is critical for ES with "Provision Connectivity: Aggregation & Local Area" requirements [51].
Connectivity Layers	Spatial data representing ecological or ES flows (e.g., larval dispersal, water flow).	Used in "Marxan with Connectivity" to prioritize interconnected sites. Essential for ES with "ES Flow" and "Regional Network" connectivity types [55].

Technical Protocol: Implementing Connectivity for ES in Marxan

This protocol provides a step-by-step guide for addressing the ES connectivity typology using Marxan's technical capabilities. The workflow translates a conceptual model of an ES into specific spatial data and software configurations.

Protocol Steps:

• Service Classification: Classify the target ES according to the typology in Table 1. For example, classify recreational ES as requiring "Provision Connectivity: Aggregation & Local Area" [51].
• Data Layer Preparation:
  • For Low Connectivity ES (e.g., Carbon): Create a simple raster layer where each planning unit's value represents the service supply (e.g., megatons of carbon stored) [51].
  • For Provision Connectivity ES (e.g., Recreation): In addition to the supply layer, this classification requires a method to induce spatial aggregation. This is primarily achieved by setting a high Boundary Length Modifier (BLM) parameter in Marxan, which penalizes fragmented solutions and favors compact reserve designs [51].
  • For ES Flow (e.g., Pollination): This requires a more complex data layer. Use "Marxan with Connectivity." The connectivity input should be a layer that represents the flow between service-providing areas (e.g., natural habitat supporting pollinator populations) and beneficiary areas (e.g., agricultural fields). This could be a model of pollinator dispersal distance between habitats and crops [51] [55].
• Target Setting: Establish quantitative conservation targets for each ES feature layer. For instance, "secure 90% of the current pollination service flow to high-value farmland" [50].
• Analysis and Calibration: Run Marxan multiple times (e.g., 100-1000 repetitions) while varying key parameters like the BLM to generate a range of near-optimal solutions. Use the Marxan calibration tool (available in interfaces like CLUZ or Zonae Cogito) to systematically explore this parameter space.
• Solution Evaluation: Analyze the output, particularly the summed solution raster which indicates how frequently each planning unit was selected across all runs. Assess whether the solutions adequately capture the spatial configuration required for the ES (e.g., are recreational areas sufficiently large and contiguous?).

Application Case: Coastal Zone Planning

A study on systematic conservation planning for people and nature demonstrates a practical application. The research integrated biodiversity, coastal protection services, and recreational services, with a focus on equitable benefit sharing [50].

Experimental Protocol Excerpt:

• Objective: Identify a network of priority areas that conserve critical biodiversity and maintain the provision of key ES while ensuring benefits are equitably distributed.
• Tools Used: Marxan software suite [50].
• Data Layers:
  • Biodiversity Features: Distributions of threatened species and habitats.
  • Ecosystem Services: Maps of "coastal protection" (e.g., mangrove and coral reef extent that mitigates storm surge) and "recreation" (e.g., accessibility of natural coastal areas for tourism) [50].
  • Costs: A composite cost layer reflecting the opportunity cost of conservation.
  • Equity Metrics: Data on human population distribution and socio-economic status to evaluate the distribution of benefits from conserved ES.
• Marxan Configuration:
  • Coastal protection was treated with "Provision Connectivity" (high BLM) as continuous stretches of natural infrastructure are more effective.
  • Recreation was treated with "Dispersed Supply" logic, potentially by setting targets within multiple sub-regions to ensure equitable access.
• Output Analysis: The analysis produced maps highlighting spatial priorities and trade-offs. The results likely demonstrated areas of high spatial congruence between biodiversity and ES (win-wins), as well as areas where objectives conflicted, requiring explicit negotiation and trade-off analysis facilitated by the Marxan output [50].

This case underscores the move from theory to practice, showing how the connectivity typology and Marxan protocols can be deployed in a real-world planning scenario to deliver more holistic and socially-just conservation outcomes.

Overcoming Analytical Challenges in B-ES Spatial Analysis

Within biodiversity and ecosystem service (BES) spatial congruence analysis, a fundamental challenge is the selection of appropriate measures. The use of oversimplified or misapplied proxies and indices can lead to misleading conclusions, poor conservation decisions, and an incomplete understanding of the relationships between biodiversity and the services ecosystems provide. This pitfall arises from a frequent disconnect between the complex, multifaceted nature of ecological systems and the simplistic metrics often used to represent them [4]. This application note details the sources and implications of this pitfall and provides protocols for selecting robust, mechanistically justified proxies in spatial BES research.

The Problem: Common Pitfalls in Proxy and Index Selection

Research into biodiversity-ecosystem service (B-ES) relationships has historically relied on top-down, correlative approaches across large spatial scales, often ignoring the underlying ecological mechanisms [4]. This has led to the widespread use of inappropriate indices, which can be categorized into two main areas.

• Table 1: Common Pitfalls in Biodiversity and Ecosystem Service Proxies

Category	Common Pitfall	Representative Example	Core Issue
Ecosystem Service (ES) Proxies	Confusing ecosystem functions (EF) for final ES [4]	Using net primary productivity (an EF) as a direct proxy for the final ES of "climate regulation" [4]	The complete pathway from function to service is not established or measured.
	Using inconsistent or non-standardized ES indices [4]	Employing different metrics for "water purification" across studies, preventing synthesis [4]	Limits generalizability and meta-analysis of B-ES findings.
Biodiversity Proxies	Over-reliance on simplistic taxonomic indices [4]	Using only species richness, which ignores species identity, abundance, and functional roles [4]	Fails to capture the mechanisms by which biodiversity influences ecosystem functions and services.
	Using policy-relevant instead of ecologically-relevant measures [4]	Using threatened species richness to infer mechanisms driving carbon storage [4]	The species driving a service may not be the same as those of conservation concern.
	Focusing on taxonomically incongruent groups [4]	Linking tree diversity to game production, or mammal diversity to carbon storage [4]	The studied taxonomic group may not be the one that mechanistically drives the target ES.

A critical, cross-cutting issue is the scale mismatch. B-ES research is often conducted at large, policy-relevant scales, while the ecological functions that underpin services operate at much finer scales [4]. Using land-cover data as a coarse proxy for ES, for instance, provides a poor fit to actual ES data and obscures within-ecosystem relationships [4]. Furthermore, different taxonomic groups can respond to different environmental gradients, meaning a richness pattern for one group (e.g., birds) may not be a valid proxy for the richness patterns of understudied groups like soil microorganisms [57].

Experimental Evidence: A Case Study in Divergent Patterns

A metabarcoding study in Amazonia highlights the risks of assuming congruence across taxa [57]. The study compared Operational Taxonomic Unit (OTU) richness from soil, litter, and insect samples with the species richness of trees and birds across a west-to-east transect and different habitat types.

• Table 2: Comparison of Richness Patterns Across Taxa in Amazonia [57]

Taxonomic Group	Data Type	West-to-East Gradient	Pattern Across Habitat Types (Terra-firme > Várzea > Igapó > Campinas)
Trees	Species Richness	Declining	Yes
Birds	Species Richness	Declining	Yes
Micro-organisms & Insects	DNA-based OTU Richness	Declining	No (Pattern was decoupled; campinas showed high OTU richness)
Conclusion	Large-scale geographic gradients may be consistent, but fine-scale, habitat-level patterns are often taxon-specific and cannot be extrapolated.

This evidence demonstrates that while large-scale diversity gradients might be consistent, the patterns at the habitat scale are highly idiosyncratic. Using trees or birds as a universal proxy for all biodiversity, particularly for understudied micro-organisms, would lead to incorrect conclusions about the distribution of diversity in these specific habitats [57].

Recommended Protocols and Workflows

To overcome the pitfall of inappropriate proxies, researchers should adopt a more mechanistic, bottom-up approach. The following workflow provides a structured method for selecting robust proxies in BES spatial congruence studies.

Protocol 1: Mechanistic Proxy Selection for BES Spatial Analysis

Title: Bottom-up Workflow for Selecting Robust Biodiversity and Ecosystem Service Proxies.

Principle: Move from correlative to mechanistic understanding by linking biodiversity to final ecosystem services via the underlying ecosystem functions they drive [4].

Procedure:

• Define the Final Ecosystem Service: Clearly articulate the specific, final service being assessed (e.g., "water purification for drinking," not "water quality") [4]. This is the endpoint of the workflow.
• Identify Underlying Ecosystem Functions: List the key ecological processes that underpin the final service. For "water purification," this could include nutrient cycling, sediment retention, and pathogen removal [4].
• Identify Key Biotic Drivers: For each ecosystem function, identify the specific taxonomic or functional groups that are known to be key mechanistic drivers. For example, the diversity and composition of soil microbial communities and the root structures of riparian plants are key biotic drivers of nutrient filtration [4].
• Select Mechanistically-Linked Biodiversity Metrics: Choose biodiversity indices that accurately reflect the identified biotic drivers. Move beyond simple species richness to include:
  • Functional Diversity Indices: Functional richness (FRic), functional divergence (FDiv), community-weighted mean (CWM) traits [4].
  • Structural Diversity Metrics: For example, foliage height diversity derived from LiDAR data, which has been shown to positively influence bird diversity at specific spatial scales [58].
• Spatial Scale Alignment: Ensure the spatial scale of the biodiversity data matches the scale at which the biotic drivers influence the ecosystem function. A scale mismatch is a primary cause of proxy failure [4] [58]. Re-assess metrics if scales are not aligned.

Protocol 2: Empirical Validation of Proxy Relationships

Title: Field Validation of Proxy-ES Relationships.

Objective: To empirically test the assumed link between a selected biodiversity proxy and the target ecosystem service across the study area.

Procedure:

• Stratified Sampling: Design a field sampling strategy that captures the expected gradient in the biodiversity proxy and the associated environmental conditions.
• Co-measurement: At each sampling site, simultaneously measure:
  • The proposed biodiversity proxy (e.g., functional trait diversity, soil microbial OTU richness).
  • The target final ecosystem service (or a well-validated, direct metric thereof).
  • Relevant abiotic covariates (e.g., soil pH, topography, land-use history).
• Statistical Modeling: Use models (e.g., structural equation modeling) to test the strength and consistency of the relationship between the proxy and the service, while accounting for covariate effects. This validates whether the proxy is a reliable predictor.

• Table 3: Key Research Reagent Solutions for Robust BES Analysis

Item Name	Function / Application	Key Considerations
Functional Trait Databases (e.g., TRY, BIEN)	Provide species-level trait data to calculate functional diversity indices (FRic, FDiv, CWM) [4].	Ensure trait data is relevant to the ecosystem functions of interest (e.g., root depth for erosion control).
High-Resolution Remote Sensing Data (e.g., LiDAR, Hyperspectral)	To map structural diversity (e.g., foliage height diversity), habitat complexity, and land cover at fine spatial resolutions [58].	Resolution must be fine enough to capture the vertical and horizontal structure relevant to the biotic drivers [58].
DNA Metabarcoding Kits	To characterize the diversity of understudied but functionally critical groups like soil bacteria, fungi, and microfauna [57].	Choice of molecular marker (16S, 18S, COI) determines the taxonomic groups targeted [57]. Be aware of database gaps.
Geospatial Analysis Software (e.g., R, QGIS, ArcGIS)	To manage, analyze, and model spatial congruence between biodiversity and ES layers, while accounting for scale dependencies.	Must be capable of handling multi-scale analysis and spatial statistics.
ColorBrewer / Viz Palette	To ensure data visualizations are clear, accessible, and accurately represent quantitative patterns without distortion [59] [60].	Use color palettes appropriate to data type (qualitative, sequential, diverging) and test for colorblind safety [59].

A critical challenge in biodiversity and ecosystem service (BES) spatial congruence analysis is the scale- and context-dependent nature of the relationships we aim to measure. Scale-dependence refers to how observed ecological patterns and relationships change across spatial, temporal, and organizational scales of measurement [61]. Context-dependence acknowledges that these relationships are modulated by local environmental conditions, historical factors, and biotic interactions. For researchers investigating spatial congruence between biodiversity and ecosystem services, failing to account for these factors can lead to erroneous conclusions and ineffective conservation policies.

Empirical studies demonstrate that the direction and strength of BES relationships can shift dramatically across scales. An analysis of multi-assemblage comparisons found that rarefaction curves frequently crossed, indicating reversals in species richness rankings across spatial scales [61]. Nearly 10% of studies examining biodiversity change in response to ecological drivers showed changes in direction across scales, where a positive relationship at one scale became negative at another [61]. Furthermore, research on resource-biodiversity relationships has documented transitions from unimodal patterns at local scales (<20 km) to monotonically increasing patterns at regional scales (>4,000 km) [62]. These findings underscore the critical importance of scale-explicit approaches in spatial congruence analysis.

Theoretical Foundation

Conceptual Framework for Scale-Dependence

The theoretical basis for scale-dependence rests on several interconnected ecological principles. The multi-species contact process model illustrates how stochastic spatial mechanisms involving local dispersal, competition, and resource heterogeneity can generate scale-dependent patterns even without niche differentiation [62]. This neutral model demonstrates that unimodal resource-richness relationships at community scales emerge from trade-offs between resource availability and patch connectivity.

The transition in relationship patterns across scales occurs because regional species pools integrate across heterogeneous local conditions. At local scales, competitive exclusion may limit diversity at high resource levels, creating hump-shaped patterns. At regional scales, these local constraints are averaged out, revealing the underlying energy-richness relationship [62]. This explains why productivity-diversity relationships documented in pond ecosystems show unimodal patterns at individual pond levels but linear increases at watershed scales [62].

Implications for Congruence Analysis

The scale-dependence of individual BES relationships necessarily affects their spatial congruence. Congruence patterns measured at one scale may not hold at others due to:

• Differential scaling of various ecosystem services
• Threshold effects in service provision across scales
• Spatial mismatch in service supply and demand areas
• Cross-scale interactions in ecological processes

These factors necessitate explicit scale considerations when designing BES congruence research and interpreting results for conservation decision-making.

Quantitative Evidence

Table 1: Documented Scale-Dependent Relationships in Biodiversity and Ecosystem Service Studies

Relationship Type	Local/Community Scale Pattern	Regional/Landscape Scale Pattern	Key Driver of Transition
Productivity-Diversity [62]	Unimodal (<20 km)	Linearly increasing (>4,000 km)	Resource heterogeneity and species pool integration
Resource-Species Richness [62]	Unimodal	Monotonically increasing	Spatial aggregation of community sites
Biodiversity Change [61]	Direction varies	Direction may reverse (10% of studies)	Sampling and spatial effects

Table 2: Accuracy Improvement from Multi-Model Ensembles Across Scales

Ecosystem Service	Number of Models	Ensemble Accuracy Improvement	Validation Data Scale
Water Supply	8	14%	Weir-defined watersheds
Recreation	5	6%	National scale
Aboveground Carbon Storage	14	6%	Plot scale
Fuelwood Production	9	3%	National scale
Forage Production	12	3%	National scale

Methodological Protocols

Multi-Scale Assessment Protocol

Objective: To quantify biodiversity-ecosystem service relationships across multiple spatial scales in a consistent, comparable framework.

Materials:

• Geospatial data layers for biodiversity indicators
• Ecosystem service model outputs or measurement data
• GIS software with zonal statistics capability
• Statistical software (R, Python, or equivalent)

Procedure:

• Scale Definition: Define a nested hierarchy of spatial scales relevant to your research questions and management applications. Include at minimum:

  • Local/plot scale (0.1-1 km)
  • Community/landscape scale (1-20 km)
  • Regional scale (20-4,000 km)

• Data Harmonization: Resample all BES data to consistent resolution and extent using standardized methods:

  • For raster data, use conservative resampling techniques
  • For vector data, apply appropriate aggregation functions
  • Document all transformations and potential uncertainties

• Multi-Scale Analysis: At each predefined scale:

  • Calculate biodiversity metrics (species richness, Shannon diversity, functional diversity)
  • Quantify ecosystem service indicators (carbon storage, water yield, recreation potential)
  • Compute correlation coefficients between BES metrics
  • Perform multivariate analyses to identify covariance structures

• Cross-Scale Comparison:

  • Test for significant differences in relationship strength and direction across scales
  • Identify scale thresholds where relationship patterns shift
  • Map spatial concordance and discordance across scales

• Uncertainty Quantification:

  • Calculate confidence intervals for all scale-specific relationships
  • Propagate measurement and model errors through analyses
  • Report spatial autocorrelation effects at each scale

Ensemble Modeling Protocol for Uncertainty Reduction

Objective: To address the "certainty gap" in BES assessments by creating model ensembles that improve accuracy and provide uncertainty estimates [63].

Materials:

• Multiple ES models (e.g., ARIES, InVEST, Co$ting Nature)
• Validation datasets (empirical measurements, national statistics)
• Computational resources for model runs
• Ensemble calculation scripts (available at github.com/GlobalEnsembles)

Procedure:

• Model Selection: Identify multiple independent models for each ecosystem service of interest. Prioritize models with:

  • Different theoretical foundations
  • Alternative input data requirements
  • Published validation records

• Model Implementation: Run all selected models using consistent input data and spatial resolution across the study area.

• Ensemble Creation: Calculate multiple ensemble types:

  • Unweighted median: Median value across all models per grid cell
  • Unweighted mean: Mean value across all models
  • Weighted ensembles: Models weighted by past performance or theoretical justification

• Accuracy Validation: Compare ensemble predictions against independent validation data using:

  • Deviance metrics (primary)
  • Spearman's ρ (secondary)
  • Spatial error analysis

• Uncertainty Mapping: Calculate and map the standard error of the mean across models for each location as a proxy for local accuracy.

• Application: Use ensemble outputs rather than individual model results for BES congruence analysis, prioritizing weighted ensembles when possible.

Multi-Scale BES Assessment Workflow

Research Reagent Solutions

Table 3: Essential Methodological Tools for Scale-Explicit BES Research

Research Tool	Function	Application Example
Rarefaction Methods [61]	Standardizes diversity comparisons across unequal samples	Comparing species richness across differently-sampled areas
Hill Numbers [61]	Integrates species richness and evenness	Quantifying diversity in a scalable, mathematically coherent framework
Spatial Zonation Algorithms	Identifies priority areas across scales	Multi-scale conservation planning for BES congruence
Model Ensemble Algorithms [63]	Improves prediction accuracy	Reducing uncertainty in ES mapping across scales
Multi-Scale Validation Datasets	Assesses accuracy across scales	Testing scale-transference of BES models

Analytical Framework

Scale-Dependence in Biodiversity Relationships

Scale-Explicit Statistical Approaches

Implement statistical methods that explicitly incorporate scale:

• Multi-level models with random effects for spatial hierarchies
• Cross-scale interaction terms in regression models
• Multivariate ordination with scale as a constraint
• Spatially-weighted regression with varying kernel sizes

Scale Transition Theory Application

Apply scale transition theory to predict how local relationships aggregate to regional patterns:

• Quantify nonlinear averaging effects
• Account for spatial covariance between processes and environments
• Model how response shapes change with scale
• Predict congruence pattern shifts across administrative boundaries

Accounting for scale and context-dependence is not merely a methodological consideration but a fundamental requirement for robust BES congruence analysis. The protocols presented here provide a structured approach to navigate scale-related challenges, while the ensemble modeling strategy directly addresses critical certainty and capacity gaps in current research practice [63]. By adopting these scale-explicit frameworks, researchers can produce more accurate, applicable insights for conservation prioritization and environmental policy across administrative boundaries and ecological jurisdictions.

Ecosystem service (ES) hotspots are geographic areas that provide disproportionately high levels of ecosystem services relative to the surrounding landscape. Within the broader context of biodiversity and ecosystem service spatial congruence analysis, identifying these areas is critical for prioritizing conservation efforts and understanding the natural capital that underpins human well-being. This protocol provides a standardized methodology for defining, identifying, and comparing ES hotspots, enabling researchers to produce consistent, repeatable, and comparable results across different geographical regions and scales.

Defining Ecosystem Service Hotspots: Concepts and Criteria

The conceptual definition of an ES hotspot is an area of high service provision. Operationally, this is typically translated into a quantitative threshold. A common approach is to classify areas falling above a specific percentile of the frequency distribution of ES values as hotspots. Table 1 summarizes standard quantitative criteria for hotspot identification.

Table 1: Common Quantitative Thresholds for Defining Ecosystem Service Hotspots

Threshold Criterion	Description	Key Application/Consideration
Top X% of Values	Areas with ES supply in the highest X percentile (e.g., top 10%, 20%) are classified as hotspots [64].	A widely used and intuitive method; allows for relative comparison within a study area. The choice of percentile (e.g., 10% vs. 20%) influences the extent of the identified hotspot.
Absolute Threshold	Areas exceeding a predefined, absolute biophysical or economic value [65].	Useful for setting management goals based on specific, measurable targets (e.g., soil erosion below 5 t/ha/year). Requires established baseline data.
Statistical Outlier	Areas identified as statistical outliers (e.g., using Z-scores) from the mean ES supply.	Effective for identifying areas of exceptionally high service provision that are distinct from the norm.

The following diagram illustrates the logical workflow and decision points involved in the process of defining an ecosystem service hotspot.

Methodological Protocol for Hotspot Identification and Comparison

This section provides a detailed, step-by-step protocol for mapping ecosystem services and identifying hotspots, incorporating insights from recent large-scale assessments.

Phase I: Spatial Assessment of Ecosystem Services

Objective: To generate high-resolution, spatially explicit maps of ecosystem service supply.

• Step 1: Define the Spatial Extent and Resolution. Determine the geographical boundaries of your study area (e.g., regional, national) and the spatial resolution of the analysis. High-resolution datasets (e.g., 30m) are critical for identifying site-specific differences and avoiding the obscuring of small-scale hotspots [65].
• Step 2: Select Target Ecosystem Services. Choose which services to map based on the research objectives and relevance to the study area. Common categories include:
  • Provisioning Services: Water yield, food production, timber [66] [65].
  • Regulating Services: Net primary productivity (NPP), soil conservation, sandstorm prevention, carbon sequestration, nitrogen processing [67] [65].
• Step 3: Data Collection and Model Parameterization. Gather input data required by the chosen ecological process models. This typically involves:
  • Land Cover/Land Use Data: High-resolution maps.
  • Climate Data: Precipitation, temperature, solar radiation.
  • Topographic Data: Digital Elevation Models (DEMs), slope.
  • Soil Data: Soil type, texture, organic matter content.
  • Ground Truthing: Use in situ observations and monitoring data to calibrate and validate model parameters [65].
• Step 4: Model Execution and Mapping. Run the calibrated models (e.g., InVEST, ARIES) to simulate and map the biophysical supply of the selected ecosystem services across the study area for the specified time period(s) [65].

Phase II: Identifying and Comparing Hotspots

Objective: To apply quantitative thresholds to ES maps to identify hotspots and analyze their spatial relationships.

• Step 5: Apply Hotspot Threshold. For each ecosystem service map, apply the chosen quantitative threshold from Table 1 (e.g., the top 20% of values) to create a binary map (hotspot vs. non-hotspot).
• Step 6: Conduct Spatial Overlay Analysis. To identify areas of multiple ES provision (ES bundle hotspots) or congruence with biodiversity, use Geographic Information System (GIS) overlay techniques:
  • ES Bundle Hotspots: Overlay individual ES hotspot maps. The intersecting areas are classified as multi-ES hotspots.
  • Biodiversity-ES Congruence: Overlay the ES hotspot map(s) with a map of biodiversity hotspots. Biodiversity hotspots are defined as regions with at least 1,500 endemic vascular plant species that have lost ≥70% of their primary vegetation [68] [69] [70].
• Step 7: Change Detection Analysis (Temporal Comparison). To analyze changes in hotspot location and extent over time, repeat Steps 1-6 for different time points (e.g., 2000, 2010, 2020). Compare the hotspot maps to quantify spatial shifts, gains, and losses [65]. This can reveal how environmental changes or management interventions affect service provision.

The following workflow diagram summarizes this two-phase protocol, highlighting the key procedures and data outputs at each stage.

Experimental Protocols from Key Studies

Protocol: Nesting Experiments in a Mapped Landscape to Assess BEF Relationships

This protocol, adapted from a study on marine sediments, demonstrates how to empirically validate the functional importance of identified hotspots by measuring Biodiversity-Ecosystem Function (BEF) relationships across a naturally varying landscape [67].

• Objective: To test how BEF relationships for nitrogen cycling change across a landscape and with increasing nutrient loading.
• Methodology:
  • Landscape-Scale Survey: Conduct an intensive grid-based survey (e.g., 400 points across 300,000 m²) of the macrobenthic community composition and sediment properties [67].
  • Identification of Experimental Locations: From the survey data, select ~28 experimental locations that represent a gradient of community composition, functional diversity, and the abundance of key ecosystem engineers [67].
  • Manipulative Experiment: At each location, establish experimental plots. Manipulate an environmental stressor relevant to the system, such as sediment nitrogen load, using multiple treatment levels (e.g., control, 150 g N m⁻², 600 g N m⁻²) applied via slow-release fertilizer injections [67].
  • Functional Response Measurement: After a suitable incubation period (e.g., 7 weeks), sample the plots to measure ecosystem functions. Key indicators for nitrogen cycling included:
    • Pore water ammonium concentration (surface and deep sediments).
    • Ammonium flux across the sediment-water interface.
    • Benthic oxygen consumption.
    • Standing stock of microphytobenthos (measured as chlorophyll a) [67].
  • Spatial Extrapolation: Use the BEF relationships derived from the experiment to predict functional performance (e.g., location of hot and cold spots) across the entire mapped landscape under different environmental scenarios [67].

Protocol for High-Resolution ES Dataset Creation

A recent study creating a high-resolution dataset for China offers a protocol for large-scale, spatially explicit ES assessment [65].

• Objective: To produce a high spatial resolution (30m) dataset of multiple ecosystem services over a 20-year period.
• Methodology:
  • Model Selection and Parameterization: Use established ecological process models to simulate net primary productivity (NPP), soil conservation, sandstorm prevention, and water yield. Model parameters were calibrated based on synthesized literature, ground monitoring data, and reconstructed remote sensing data [65].
  • Data Inputs: Utilize time-series data on land cover, climate (precipitation, temperature), topography (DEM, slope), and soil properties as primary inputs for the models [65].
  • Validation: Validate the simulated ES maps against in situ observations and compare them with existing datasets to ensure accuracy [65].
  • Trend Analysis: Calculate the change in ES supply over the study period to identify regions of significant improvement or degradation [65].

The Scientist's Toolkit: Research Reagents and Essential Materials

Table 2 details key reagents, datasets, and software solutions essential for conducting ES hotspot analyses.

Table 2: Essential Research Tools for ES Hotspot Analysis

Item/Reagent	Function/Description	Application Example
GIS Software (e.g., ArcGIS, QGIS)	The primary platform for spatial data management, analysis, cartography, and overlay operations.	Used in all phases, from processing input data to executing spatial overlays and mapping final hotspots [66] [65].
Ecological Process Models (e.g., InVEST, ARIES)	Software suites containing pre-defined models to simulate and map the biophysical supply of ecosystem services.	Used in Phase I to generate maps of water yield, soil retention, and other services [66].
Slow-Release Fertilizer (e.g., Urea, Nutricote)	Used in manipulative field experiments to simulate environmental stressors like eutrophication in a controlled manner.	Applied in the BEF experiment to create gradients of sediment nitrogen loading [67].
High-Resolution Land Cover Map	A fundamental dataset representing the physical material on the Earth's surface, a key input for most ES models.	Critical for accurately simulating ecosystem processes at local scales (e.g., 30m resolution) [65].
Global Biodiversity Hotspots GIS Data	A vector dataset depicting the 36 globally recognized biodiversity hotspots.	Used as a baseline layer for analyzing spatial congruence between biodiversity and ES hotspots [68] [69] [70].
Pore Water Sampler (e.g., syringe corer)	A device for collecting interstitial water from sediments for chemical analysis.	Used to measure pore water ammonium concentrations in the BEF experiment [67].

Optimizing for Cost-Effectiveness and Compactness in Conservation Planning

Application Notes

Conservation planning aims to efficiently allocate limited resources to achieve the greatest possible ecological benefits. Two central, and often competing, objectives in this process are cost-effectiveness—maximizing conservation outcomes per unit of cost—and compactness—designating protected areas with minimal fragmentation to enhance ecological viability and reduce management expenses. The following notes detail the application of advanced analytical frameworks to balance these objectives, particularly within the context of biodiversity and ecosystem service spatial congruence.

The Strategic Value of Spatial Congruence Analysis

Spatial Congruence Analysis (SCAN) is a chorological method that identifies groups of species (chorotypes) sharing similar geographical distributions by quantifying the spatial overlap of their range polygons [22]. This method moves beyond grid-based analyses to minimize scale distortion and uses an explicit, quantifiable Spatial Congruence Index (CS). The index is calculated as: CS,ab = (Oab/Aa) × (Oab/Ab) where Oab is the area of overlap between species a and b, and Aa and Ab are their respective range areas [22]. This approach allows planners to identify core areas of high biodiversity congruence, which are prime candidates for compact and cost-effective protection.

Quantifying Conservation Gaps and Target Setting

Global analyses reveal significant gaps in current conservation efforts. For wetlands, a critical ecosystem, only 44.03% of identified global Wetland Conservation Priorities (WCPs) are within existing protected areas, leaving 55.97% unprotected [71]. Research on the Qinghai-Tibet Plateau further shows that while existing Nature Reserves (NRs) offer some protection, they fail to safeguard all critical zones for biodiversity and ecosystem services, highlighting the need for systematic gap analysis [72].

To address these gaps, science-based target-setting scenarios are essential. For global wetland conservation, three distinct scenarios have been proposed, which can be downscaled to national levels [71]:

• Conservative Target: Increases WCP coverage by 9.40%.
• Moderate Target: Increases WCP coverage by 42.40%.
• Ambitious Target: Increases WCP coverage by 55.97%.

Integrated Frameworks for Multiple Conservation Goals

A singular focus on biodiversity can overlook other critical ecological and socio-economic values. An integrative "representativeness-vulnerability" framework, which combines biodiversity, ecosystem services (ES), and their susceptibility to human activity, facilitates the identification of priority areas that deliver multiple benefits [72]. A study in the Mira River watershed, Ecuador, demonstrated a positive spatial relationship between biodiversity and soil accumulation service in 98% of subwatersheds, with biodiversity explaining up to 92% of the variance in soil accumulation supply. A 52.5% spatial overlap was found, indicating that 15% of the subwatersheds are optimal for simultaneous management of both goals, thereby optimizing management investment [73].

Empirical evidence from German grasslands confirms that cost-effectiveness optimization can achieve high levels of multiple conservation goals—including species richness, trophic interactions, and ecosystem resilience—with limited budgets. This approach can secure up to 80% of the maximum achievable ecological benefits using only 30% of the budget required for an unconstrained ecological optimization [74]. While trade-offs are inevitable, they are minimized when planning jointly for multiple goals.

Table 1: Key Quantitative Findings from Recent Conservation Planning Research

Study Focus / Region	Key Metric	Finding	Source
Global Wetland Conservation	Current Protection of Priorities	44.03% of WCPs are protected	[71]
Global Wetland Conservation	Ambitious Target Coverage	55.97% increase in WCP coverage	[71]
Grassland Conservation, Germany	Cost-Effectiveness	Achieves 80% of ecological benefits with 30% of budget	[74]
Mira Watershed, Ecuador	Biodiversity & Soil Service Overlap	52.5% spatial congruence	[73]
Mira Watershed, Ecuador	Subwatersheds for Co-Management	15% are optimal for joint management	[73]

Experimental Protocols

This section provides detailed methodologies for the key computational and analytical procedures cited in the application notes.

Protocol A: Spatial Congruence Analysis (SCAN) for Chorotype Identification

Purpose: To identify groups of species (chorotypes) with significantly congruent geographic distributions without a priori assumptions, establishing a foundation for compact reserve design [22].

Workflow Overview:

• Data Preparation: Compile species distribution data as polygon shapefiles.
• Congruence Calculation: Compute pairwise Spatial Congruence Index (CS) for all species.
• Network Construction: Build a network where nodes represent species and edges represent CS values above a threshold.
• Chorotype Delineation: For a reference species, iteratively identify directly and indirectly connected species at decreasing congruence thresholds to form "partial chorotypes."

Materials & Software:

• Input Data: Polygon shapefiles of species ranges.
• Software: R environment with sf package for geospatial operations [22].
• Output: A list of chorotypes at various congruence levels, with associated metrics (e.g., depth, richness).

Step-by-Step Procedure:

• Data Input and Validation: Load species range polygons into R. Ensure all layers are in the same projected coordinate system to ensure accurate area calculations.
• Calculate Pairwise Congruence: For every unique pair of species (a, b), calculate the CS index using the formula provided in Section 1.1. Store results in a symmetric matrix.
• Set Congruence Threshold (CT): Begin with a high CT (e.g., 95%).
• Identify Chorotypes for a Reference Species: a. Select a reference species. b. At the current CT, find all species directly congruent with the reference species (CS ≥ CT). c. From this set, find all species directly congruent with each new member, iterating until no new species are added, forming a "closed list" or partial chorotype. d. Record the partial chorotype and its member species. e. Decrease the CT by a fixed increment (e.g., 1%) and repeat steps b-d until the group fails to close or becomes overly inclusive.
• Repeat: Execute Step 4 for every species in the dataset.
• Post-Processing and Analysis: Analyze the collected chorotypes for nestedness and gradients. Calculate descriptive metrics for each chorotype, such as:
  • Congruence: The CT value for the chorotype.
  • Depth: The range of CT values over which the chorotype persists.
  • Richness: The number of species in the chorotype.

Protocol B: Integrated Conservation Prioritization using Marxan with Human Interference Factors

Purpose: To identify cost-effective and compact priority conservation areas that meet biodiversity targets while minimizing conflict with human activities [75].

Workflow Overview:

• Habitat Suitability Modeling: Use MaxEnt to model species habitat suitability.
• Cost Layer Creation: Develop a unified cost layer representing economic and social constraints.
• Conservation Prioritization: Use Marxan to generate efficient reserve networks that represent species in a compact manner while minimizing total cost.

Materials & Software:

• Species Data: Validated species occurrence points.
• Environmental Data: Bioclimatic, topographic, and soil variables (e.g., from WorldClim).
• Cost Data: Layers for land use, nighttime lights, population density, road density, etc.
• Software: MaxEnt for species distribution modeling; Marxan for spatial prioritization; ArcGIS or R for data preprocessing and post-processing.

Step-by-Step Procedure: I. Model Species Habitat Suitability with MaxEnt: 1. Data Collection: Gather species occurrence points and screen for accuracy, removing duplicates and cultivated records [75]. 2. Variable Selection: Obtain and process ecological variables (e.g., 19 bioclimatic variables, elevation, soil data). Perform multicollinearity analysis to select a final, non-correlated set for modeling. 3. Model Calibration and Run: Use the R package ENMeval or kuenm to tune MaxEnt parameters (feature classes, regularization multiplier) and avoid overfitting. Run the optimized MaxEnt model to generate habitat suitability maps for each species.

II. Develop a Unified Human Interference Factor (HIF) Cost Layer: 1. Data Acquisition: Collect spatial data representing human pressure: * Land use/cover (e.g., from Globeland30) * Nighttime light intensity * Population density * GDP density * Road density * Human activity intensity index 2. Data Standardization: Resample all layers to a consistent spatial resolution and extent. Normalize each layer to a common scale (e.g., 0-1). 3. Weighted Integration: Use the Entropy Weight Method (EWM) to assign objective weights to each factor based on its information content. Combine the weighted layers into a single HIF raster using a weighted sum, where higher values indicate higher conservation cost (greater human interference) [75].

III. Prioritize Conservation Areas with Marxan: 1. Define Planning Units: Subdivide the study area into planning units (e.g., hexagons or watershed sub-basins). 2. Set Conservation Targets: Define representation targets for each species (e.g., protect 20-30% of its suitable habitat, defined from the MaxEnt output). 3. Assign Input Values to Planning Units: * Conservation Feature Amount: For each species, calculate the area of suitable habitat in each planning unit. * Cost: Assign the mean HIF value from the unified cost layer to each planning unit. 4. Configure and Run Marxan: Set Marxan parameters, including the Boundary Length Modifier (BLM), which is crucial for controlling the compactness of the selected priority areas. A higher BLM value favors more compact solutions. Run Marxan multiple times to generate a range of near-optimal solutions. 5. Analyze Results: The primary output is the "selection frequency" of each planning unit—how often it was chosen across multiple runs. Areas with high selection frequency (e.g., 80-100) are considered irreplaceable for achieving conservation targets cost-effectively and compactly [75].

Visualization of Workflows

The following diagrams illustrate the logical relationships and sequential steps in the key protocols described above.

Figure 1. SCAN Chorotype Identification Workflow. This diagram outlines the iterative process of identifying species groups with congruent distributions at varying thresholds of spatial overlap.

Figure 2. Integrated Conservation Prioritization Workflow. This diagram shows the parallel processes of modeling ecological data and socio-economic costs, which are combined in a systematic prioritization analysis.

The Scientist's Toolkit: Research Reagent Solutions

Table 2: Essential Materials and Data Sources for Conservation Planning Analyses

Item Name / Category	Specific Source or Tool	Function in Analysis	Key Consideration
Species Data Portals	Global Biodiversity Information Facility (GBIF)	Provides primary species occurrence point data for modeling.	Data requires rigorous cleaning for geographic inaccuracies and biases.
	National Plant Specimen Resource Center (CVH)	National/regional repository for validated specimen data.	Often contains data not available on international platforms.
Environmental Data	WorldClim Database	Source of current, historical, and future bioclimatic variables.	Choice of variable set and resolution (e.g., 1 km) is critical to model performance.
	National Tibetan Plateau / Other Data Centers	Source for specialized data (e.g., soil properties, vegetation indices).
Human Impact Data	Globeland30 Land Cover	High-resolution (30m) global land use/cover data for cost modeling.	Key for quantifying habitat conversion and human footprint.
	Nighttime Light Data (e.g., VIIRS)	Proxy for human economic activity and urbanization intensity.	Effective for creating a continuous surface of development pressure.
	GRIP Global Road Database	Provides global road density data, a key driver of habitat fragmentation.
Software for Modeling	MaxEnt	Uses maximum-entropy principle to model species habitat suitability.	Requires parameter tuning (ENMeval) to avoid overfitting.
	Marxan	Spatial conservation prioritization software for optimal reserve design.	Boundary Length Modifier (BLM) is key parameter for compactness.
	R Software Environment	Platform for data cleaning, analysis (e.g., SCAN method), and visualization.	Essential for reproducible and customizable analytical workflows.
Spatial Analysis Tools	ArcGIS / QGIS	Geographic Information System for data management, preprocessing, and map production.	Industry standard for spatial data handling and cartography.

The concept of "bundling services" in ecological contexts refers to the spatial aggregation of multiple ecosystem services that repeatedly appear together across a landscape. Understanding these bundles is critical for spatial conservation planning because it reveals where managing for one service may inherently support or undermine other services and biodiversity conservation goals. Research demonstrates that biodiversity has a positive spatial relationship with key ecosystem services like soil accumulation, with biodiversity explaining up to 92% of the variance in soil accumulation service supply in some tropical watersheds [73]. However, significant trade-offs emerge when conservation priorities are unbalanced; studies show that putting too much weight on ecosystem services can have detrimental effects on biodiversity conservation [5].

Ecosystem service bundles emerge from complex social-ecological interactions that operate across multiple spatial scales. The identification and management of these bundles requires understanding both ecological/supply-side drivers and social-ecological/demand-side drivers that influence how biodiversity-ecosystem service (BES) relationships manifest across different contexts [76]. Properly defining the "service shed" – the appropriate spatial and temporal context for quantifying a service – is essential for accurate assessment of these relationships and for informing effective conservation policies [77].

Analytical Framework for Identifying Service Bundles

Conceptual Foundations

The relationship between biodiversity and ecosystem services operates through a cascade where species composition drives ecosystem functions, which in turn generate services that benefit human populations. This relationship is spatially explicit and scale-dependent, meaning that observed BES relationships can strengthen, weaken, or change directionality depending on the spatial scale of analysis [76]. The Forest Human Nexus (FHN) exemplifies this complexity, integrating forest area per capita, forest accessibility, and population density to provide a geospatial assessment of human-forest relationships [78].

Table 1: Key Concepts in Bundle Service Analysis

Concept	Definition	Spatial Consideration
Ecosystem Service Bundle	Clusters of ecosystem services that repeatedly appear together across space [79]	Represents spatial covariance of multiple services
Service Shed	Appropriate spatial and temporal context for quantifying a service flow [77]	Boundaries depend on networks connecting ecosystem supply to human beneficiaries
Spatial Congruence	Degree of spatial overlap between biodiversity and ecosystem services [73]	Enables identification of win-win management areas
Trade-off	Situation where increase in one service leads to decrease in another [5]	Often manifests when management prioritizes single services
Synergy	Situation where multiple services reinforce each other [79]	Creates opportunities for efficient multi-service management

Spatial Analysis Workflow

The following diagram illustrates the comprehensive workflow for identifying and analyzing ecosystem service bundles:

Figure 1: Ecosystem Service Bundle Analysis Workflow. This diagram outlines the sequential process for identifying and analyzing ecosystem service bundles, from data collection to management application.

Experimental Protocols for Bundle Analysis

Multi-Dimensional Data Collection Protocol

Comprehensive bundle analysis requires integrating data across multiple dimensions. The LESS (Landscape-Ecology-Society-Space) framework provides a structured approach for capturing these dimensions [79]:

• Landscape Metrics: Calculate landscape pattern indices using Fragstats 4.0 software, including total area, landscape density, landscape shape index, aggregation index, Shannon diversity index, and Shannon evenness index derived from high-resolution (0.5m) remote sensing data with land cover classification into forest, grass, water, bare land, and impervious surface [79].

• Ecological Parameters: Quantify vegetation cover using NDVI, surface temperature using LST, humidity using WET index, and surface imperviousness using NDBSI from Landsat 8 imagery. Ensure overall accuracy of ≥85% through confusion matrix validation with 50 randomly selected verification points per category [79].

• Social Dimensions: Assess population density around target areas (e.g., parks) within a 30-minute walking distance. Map distribution of public transportation stations, shopping services, restaurant services, and entertainment services to quantify social accessibility [79].

• Spatial Characteristics: Measure border frontage rate, green barrier rate, blocking wall rate, boundary waterfront rate, and boundary openness rate to capture spatial configuration effects [79].

Bundle Identification Using Self-Organizing Maps (SOM)

The Self-Organizing Map algorithm is particularly effective for identifying ecosystem service bundles from high-dimensional datasets:

• Data Normalization: Standardize all variables (landscape, ecological, social, spatial) to a common scale (0-1) to prevent dominance by variables with larger numerical ranges.

• SOM Training: Initialize a 2D neural network with hexagonal nodes. The number of nodes should be determined by the formula 5√n, where n is the number of observation units. Train the network using the batch algorithm for 100-500 iterations.

• Bundle Identification: Cluster the trained SOM nodes using k-means clustering. Determine the optimal number of clusters (bundles) through the elbow method based on within-cluster sum of squares.

• Bundle Characterization: Calculate mean values of all variables within each identified bundle to interpret the dominant service combinations and their spatial manifestations [79].

Trade-off and Synergy Analysis Protocol

Quantifying relationships between services and biodiversity requires multiple analytical approaches:

• Geographically Weighted Regression (GWR): Implement local regression modeling to analyze spatial heterogeneity in BES relationships. Use adaptive bandwidth selection to account for varying data density across the study area. The GWR model takes the form: yᵢ = β₀(uᵢ,vᵢ) + Σβₖ(uᵢ,vᵢ)xᵢₖ + εᵢ where (uᵢ,vᵢ) denotes the coordinates of each location i [79] [73].

• Spatial Congruence Analysis: Calculate overlap coefficients between biodiversity hotspots and ecosystem service provision areas using the formula: C = (A∩B)/(A∪B) × 100% where A represents areas of high biodiversity value and B represents areas of high ecosystem service supply [73].

• Temporal Trend Analysis: For longitudinal assessments, compute Sen's slope estimator to quantify monotonic trends in metrics like Forest Area per Person (FAP): Q = median((xⱼ - xᵢ)/(j - i)) for all i < j where x represents the metric value at different time points [78].

Visualization of Bundle Relationships

Conceptual Relationships in Service Bundles

The following diagram illustrates the key relationships and trade-offs within ecosystem service bundles:

Figure 2: Ecosystem Service Bundle Relationships and Trade-offs. This diagram illustrates the cascading relationships from biodiversity to human well-being and key intervention points for managing trade-offs.

Application to Conservation Planning

Priority Area Identification

Spatial congruence analysis enables identification of areas where biodiversity and ecosystem service management can be simultaneously optimized. Research in the Rio Mira watershed of Ecuador demonstrated that 52.5% overlap occurred between biodiversity and soil accumulation service, with 15% of subwatersheds identified as high priority for simultaneous management of both objectives [73]. This approach allows conservation planners to:

• Target limited resources to areas providing the greatest joint benefits
• Avoid areas where intense trade-offs would undermine conservation goals
• Design multi-objective conservation networks that address both biodiversity and human wellbeing

Protected area network assessments often reveal spatial misalignment with integrated priorities. Analyses show that protected areas are frequently not optimally located to ensure prioritization of both ecosystem services and biodiversity, necessitating systematic spatial planning to identify complementary areas [5].

Management Strategies for Different Bundle Types

Table 2: Management Approaches for Different Bundle Types

Bundle Type	Characteristics	Management Strategy	Potential Trade-offs
Landscape-Ecological	High landscape pattern indices, strong ecological effects [79]	Conservation-focused management; maintain habitat connectivity	Potential conflict with social development needs
Social-Spatial	High social accessibility, strong spatial effects [79]	Multi-use management; focus on recreational services	Potential ecological degradation from high use
Composite-Driven	Balanced landscape, ecological, social, spatial effects [79]	Integrated adaptive management; maintain balance across objectives	Requires more complex, resource-intensive governance
Biodiversity-Led	High species richness supporting ecosystem services [73]	Biodiversity-focused protection; habitat restoration	May limit intensive service utilization
Service-Rich	Multiple high-yielding ecosystem services [5]	Sustainable harvest management; optimize service flows	Potential biodiversity compromise if over-exploited

The Scientist's Toolkit: Essential Research Reagents

Table 3: Essential Research Tools for Ecosystem Service Bundle Analysis

Tool/Category	Specific Examples	Function/Analytical Purpose
GIS & Remote Sensing	Fragstats 4.0, Landsat 8, HILDA+ dataset, GHSL	Landscape pattern analysis, land cover classification, population distribution mapping [79] [78]
Statistical Analysis	Geographically Weighted Regression (GWR), Principal Component Analysis (PCA), Sen's slope estimator	Analyze spatial heterogeneity, reduce data dimensionality, quantify temporal trends [79] [78]
Bundle Identification	Self-Organizing Maps (SOM), k-means clustering	Identify ecosystem service bundles from high-dimensional data [79]
Spatial Metrics	Forest Area per Person (FAP), Forest Proximate People (FPP), Forest Human Nexus (FHN)	Quantify human-environment spatial relationships [78]
Field Validation	Object-oriented segmentation, confusion matrix validation, 50-point verification	Ground-truth remote sensing classifications and ensure ≥85% accuracy [79]
Ecosystem Service Models	ARIES framework, SEEA EA compatible models, service shed analysis	Quantify ecosystem service flows and dependencies [77]

Validation and Real-World Applications: From Corporate Impacts to Self-Recovery

Spatial congruence analysis of biodiversity and ecosystem services (BES) is a cornerstone of modern conservation planning, land-use management, and environmental policy development. The identification of areas where high-value BES overlap, known as hotspots, and areas of low provision, termed coldspots, enables policymakers to prioritize limited resources for maximum ecological and societal benefit [80] [28]. However, the spatial configuration and ultimate outcomes of these identified areas are highly sensitive to the methodological choices made during analysis. Different techniques can delineate vastly different geographies of priority, influencing conservation strategies and resource allocation [80]. This application note provides a detailed comparative analysis of prevalent hotspot identification methods, their experimental protocols, and their impact on spatial outcomes, framed within the context of a broader thesis on BES spatial congruence.

Hotspot and coldspot analyses transform complex, continuous BES data into discrete, mappable units critical for decision-making. A hotspot is typically defined as an area providing a high level of a specific service or where multiple services overlap, whereas a coldspot is an area characterized by low ecosystem service supply [80]. The choice of delineation method significantly influences the amount, spatial extent, and clustering of these areas, with implications for managing multiple distributed spots versus fewer, larger consolidated areas [80].

Table 1: Comparative Analysis of Primary Hotspot Identification Methods

Method Category	Core Principle	Key Strengths	Key Limitations	Representative Application/Context
Spatial Clustering (e.g., Getis-Ord Gi*)	Identifies statistically significant spatial clusters of high values (hot spots) and low values (cold spots) [80].	Objectively distinguishes significant clusters from random spatial patterns; provides confidence levels (e.g., p-value, z-score) [33].	Sensitive to scale and neighborhood definition; may miss important but non-clustered high-value areas.	Mapping social value clusters in Dalian, China [33]; identifying HES critical zones in the Aglar watershed, India [80].
Intensity & Summation Approach	Sums normalized values of multiple ES to create a composite map; hotspots are the highest percentiles (e.g., top 20%) [80].	Facilitates integration of multiple, disparate ES into a single, multi-service priority map.	Can lead to double-counting and mask trade-offs between individual services [80].	Used in studies analyzing synergies and trade-offs among multiple ecosystem services [80].
Expert-Based Matrix	Assigns ES capacity scores to land use/cover classes based on quantitative data and expert judgment [80].	Highly applicable in data-scarce regions; leverages existing knowledge efficiently.	Subjective and static; may not capture fine-scale, within-class variability or recent environmental changes.	Proposed as an alternative for rapid assessments where detailed spatial modeling is not feasible [80].
Quantile & Natural Jenks	Classifies ES maps into a set number of classes based on natural data breaks or quantiles [80].	Simple and intuitive; effectively highlights the highest-value areas for a single service.	Classifications are data-specific and not directly comparable across different studies or regions.	Used exclusively to extract areas with the highest flow of a single ecosystem service [80].

The selection of a method is not merely a technical decision but a consequential one that shapes the spatial narrative of conservation. For instance, a study in the Aglar watershed of the Indian Himalayan Region demonstrated that using two different weightage-based approaches on the same set of Hydrological Ecosystem Service (HES) descriptors resulted in significantly different spatial configurations of priority areas. One method classified 3.18% of the watershed as hot spots, while another identified 2.36% as hot spots, with only a fraction of the area (6.86 km²) consistently identified by both [80]. This underscores the critical impact of methodological choice on the final conservation map.

Detailed Experimental Protocols

To ensure reproducibility and rigorous spatial analysis, the following protocols outline key methodologies for hotspot identification.

Protocol 1: Social Value Hotspot Analysis using the SolVES Model

This protocol is designed to map non-material ecosystem service values, such as aesthetic or cultural value, based on survey data [33].

1. Survey Design and Data Collection:

• Objective: Collect spatially explicit perceptual data on ecosystem service values.
• Procedure: Administer a survey to a representative sample of respondents (e.g., residents, visitors). Utilize a polygon-based mapping exercise or point allocation within a defined study area. Accompany this with a questionnaire to gather demographic and preference data (e.g., value preferences for aesthetics, biodiversity, recreation).
• Data Output: A geodatabase containing respondent-derived point locations with attributed social value types and intensities.

2. Spatial Modeling with SolVES:

• Objective: Generate predictive maps of social value intensity.
• Procedure:
  • Input Preparation: Compile a suite of spatially explicit environmental variables (e.g., distance to water, elevation, slope, land cover) that are hypothesized to influence social values.
  • Model Run: Execute the SolVES model to analyze the relationship between the survey value points and the environmental variables. The model uses a maximum entropy approach to predict the relative intensity of each social value across the landscape [33].
  • Output: A series of raster maps, one for each social value type, with pixel values representing a Value Index.

3. Hotspot Delineation via Spatial Clustering:

• Objective: Identify statistically significant social value hotspots and coldspots from the SolVES output raster.
• Procedure:
  • Data Preparation: Convert the Value Index raster into a vector grid or point layer.
  • Statistical Analysis: Perform Getis-Ord Gi* spatial autocorrelation analysis. This calculates a z-score and p-value for each feature, indicating whether high or low values are clustered spatially.
  • Classification: Define hotspots as features with a high z-score and a statistically significant small p-value (e.g., p < 0.05). Coldspots are defined by a low negative z-score and a significant p-value [33].
  • Output: A polygon map of statistically significant social value hotspots and coldspots, ready for overlay and congruence analysis with biodiversity data.

Protocol 2: Hydrological Ecosystem Service (HES) Hotspot Analysis

This protocol uses physical environmental data and hydrological modeling to identify hotspots of key water-related services [80].

1. HES Quantification via Hydrological Modeling:

• Objective: Calculate the spatial explicit biophysical supply of multiple HES.
• Procedure:
  • Model Selection: Employ a calibrated hydrological model such as the Soil and Water Assessment Tool (SWAT). Required input data includes a Digital Elevation Model (DEM), soil maps, land use/land cover (LULC) data, and long-term climate data (precipitation, temperature) [80].
  • Descriptor Calculation: Run the model to simulate and map key HES descriptors such as Water Yield (WYLD), Sediment Yield (SYLD), and Base Flow (LATQ) [80].

2. Multi-Method Hotspot Identification:

• Objective: Delineate HES hotspots using different methods to assess uncertainty.
• Procedure:
  • Method 1 - Individual Descriptor Analysis: For each HES descriptor (e.g., WYLD, SYLD), normalize the values and apply a quantile classification (e.g., top 20th percentile) to define hotspots. Combine the individual hotspot maps using a weightage-based overlay (e.g., simple summation) [80].
  • Method 2 - Grouped Service Analysis: First, group individual HES descriptors into broader ecosystem service categories (e.g., group WYLD and LATQ under "Water Provision"). Calculate a composite index for each category, then apply a quantile classification to identify hotspots for the grouped service [80].
  • Uncertainty Assessment: Perform a pixel-by-pixel comparison of the hotspot maps generated by the different methods. Quantify the area of agreement and disagreement to highlight regions where hotspot identification is method-dependent [80].

3. Synergy and Trade-off Analysis:

• Objective: Understand the interactions between different HES.
• Procedure: Conduct a Pearson correlation or spatial correlation analysis between the raster layers of the different HES descriptors. A positive correlation indicates a synergistic relationship (e.g., high water yield and high base flow), while a negative correlation indicates a trade-off (e.g., high surface runoff and low water retention) [80].

The Researcher's Toolkit: Essential Reagents & Materials

Successful spatial congruence analysis relies on a suite of specialized software, data, and analytical tools.

Table 2: Key Research Reagent Solutions for Spatial Congruence Analysis

Research Reagent	Category	Primary Function	Application Context
SolVES (Social Values for Ecosystem Services) Model	Software / Model	A spatially explicit tool that integrates social survey data with environmental variables to map perceived ecosystem service values [33].	Mapping cultural, aesthetic, and recreational service hotspots in urban and natural landscapes [33].
Soil & Water Assessment Tool (SWAT)	Hydrological Model	A public-domain model that simulates the quality and quantity of surface and ground water and predicts the environmental impact of land use and management practices [80].	Quantifying provisioning and regulating HES descriptors like water yield, sediment yield, and nutrient loading at watershed scales [80].
InVEST (Integrated Valuation of Ecosystem Services and Trade-offs)	Software Suite	A family of models developed by the Natural Capital Project to map and value ecosystem services, often with lower data demands than SWAT [80].	Spatially explicit assessment of a wide range of ES, including carbon storage, habitat quality, and water purification [80].
Getis-Ord Gi* Statistic	Spatial Analysis Tool	A spatial autocorrelation statistic used to identify statistically significant clusters of high values (hot spots) and low values (cold spots) in a raster or feature dataset [33] [80].	The core analytical step for moving from a continuous value surface to a map of statistically significant hotspots/coldspots in most protocols.
GIS Software (e.g., ArcGIS, QGIS)	Platform	The primary spatial data management, visualization, and analysis environment for handling all geospatial data layers and performing overlay analyses.	Essential for all stages, from data preparation and model integration to final mapping and congruence assessment.

The spatial configuration of biodiversity and ecosystem service hotspots is not an absolute truth revealed by a single method, but a nuanced picture that emerges from the careful selection and application of analytical techniques. As demonstrated, the use of multiple methods, such as spatial clustering and intensity analysis, can reveal significant uncertainties in hotspot location and extent [80]. This comparative analysis underscores the necessity of a method-transparent, question-driven approach. Researchers must align their choice of hotspot identification protocol with the specific objectives of their study, whether it is prioritizing areas for multi-service conservation using a summation approach, understanding the drivers of perceived value with SolVES, or managing watershed services through rigorous hydrological modeling. Acknowledging and quantifying the uncertainty inherent in these methods is critical for producing robust, defensible, and effective spatial plans that can genuinely advance the integration of climate and biodiversity action on the ground [44].

The management of post-mining and quarrying landscapes is a critical component of sustainable development, responding to global policies such as the 2030 Agenda for Sustainable Development and the European Green Deal [81]. These landscapes, often initially perceived as degraded and of poor ecological value, can undergo natural recovery and provide significant ecosystem services (ES). A core challenge in this field is the validation of ecological models used to predict biodiversity and ES recovery. This application note synthesizes evidence and protocols for validating such models with field data, with a specific focus on post-mining pond ecosystems. Framed within the context of spatial congruence analysis—a method for detecting groups of species with shared geographical distributions [22]—this document provides researchers with actionable methodologies for testing and refining predictive models of ecosystem recovery.

Evidence from Post-Mining Pond Recovery Studies

Key Findings on Natural Recovery and Model Validation

Recent research on post-mining and quarry ponds in Sardinia provides a compelling case for natural self-recovery and offers a dataset for model validation. A study of 34 quarry and 14 mining ponds, abandoned between the 1960s and 1990s with no active restoration, recorded the presence of 524 animals, vascular plants, and habitats [81]. The study developed a Bioindex and an Ecosystem Services Index (ESI) to summarize findings. Two key results provide a critical test for predictive models:

• Temporal Trend: Both the Bioindex and ESI increased significantly with time since abandonment, confirming that these environments can recover naturally over time [81].
• Divergence between Biodiversity and Services: The Bioindex and ESI were poorly correlated, suggesting that an increase in species richness does not automatically translate to enhanced ecosystem service provision. This finding indicates that models focusing solely on biodiversity may require additional parameters to accurately predict ES recovery [81].

These findings underscore the necessity of validating models against multi-faceted field data that capture both biological and functional aspects of ecosystem recovery.

Quantitative Validation of Ecosystem Service Models

The validation of ES models at large scales, though rare, is feasible and has been successfully demonstrated. A continental-scale study in sub-Saharan Africa validated multiple ES models against 1,675 independent data points across 36 countries [82]. The results, summarized in the table below, provide a benchmark for model performance and highlight the importance of model complexity and the role of human demand.

Table 1: Performance of Ecosystem Service Model Ensembles in Sub-Saharan Africa [82]

Ecosystem Service	Type of Service	Number of Validation Data Points	Key Validation Finding	Impact of Model Complexity
Water Supply	Potential	736	Models showed reasonable predictive power	More complex models sometimes provided more accurate estimates
Stored Carbon	Potential	214	Models showed reasonable predictive power	More complex models sometimes provided more accurate estimates
Firewood Use	Realized	285	Human population density was as good or better predictor than ES supply models in 85% of cases	Simpler models based on demand were often sufficient
Charcoal Use	Realized	59	Human population density was as good or better predictor than ES supply models in 85% of cases	Simpler models based on demand were often sufficient
Grazing Resources	Realized	401	Human population density was as good or better predictor than ES supply models in 85% of cases	Simpler models based on demand were often sufficient

Furthermore, global ensembles of ES models for five services (including water supply, fuelwood, forage, carbon storage, and recreation) have been shown to be 2 to 14% more accurate than individual models [63]. This ensemble approach effectively addresses "certainty gaps" and "capacity gaps," especially in data-poor regions.

Experimental Protocols for Model Validation

This section outlines a detailed, transferable protocol for validating spatial models of biodiversity and ecosystem services in recovered mining ponds and similar habitats.

Field Sampling Design and Data Collection

Objective: To collect comprehensive field data on biodiversity, ecosystem services, and environmental variables for the purpose of validating and calibrating spatial ecological models.

Site Selection:

• Select a stratified random sample of post-mining and quarry ponds, ensuring representation across a gradient of key variables [81]:
  • Time since abandonment (e.g., decades: 1960s, 1980s, 1990s, 2000s)
  • Site type (mining vs. quarry pond)
  • Water permanence (temporary vs. permanent)
  • Proximity to urban and natural areas

Field Data Collection: The following data should be collected at each site:

• Biodiversity (Bioindex):
  • Animals: Conduct standardized surveys for key taxa (e.g., amphibians, insects, birds) using transects, point counts, or trapping appropriate to the target species [81].
  • Plants: Perform vegetation surveys in defined plots to record species identity and percent cover of all vascular plants [81].
  • Habitats: Map and classify distinct habitat types within the pond and its immediate riparian zone.
• Ecosystem Services (Ecosystem Service Index - ESI):
  • Quantify both services and disservices. Metrics can include [81]:
    • Habitat Provision: Based on habitat diversity and quality.
    • Water Quality Improvement: Simple water chemistry tests (e.g., pH, turbidity, heavy metals).
    • Recreation: Evidence of human use (e.g., paths, fishing gear).
  • Record disservices, such as the presence of hazardous structures or vectors for disease.
• Environmental Covariates:
  • Record site characteristics: surface area, water depth, substrate type.
  • Measure water chemistry parameters.
  • Georeference all sampling locations with high-precision GPS.

Spatial Congruence Analysis (SCAN) for Chorotype Validation

Objective: To validate model-predicted species distributions by identifying and mapping groups of species with significantly congruent geographical ranges (chorotypes) against field-observed distributions.

Methodology: The SCAN method provides a robust framework for this analysis by using direct and indirect spatial relationships among species without the scale bias associated with grids [22].

Workflow:

• Data Preparation: Convert field-surveyed species occurrence data into polygon range maps (shapefiles) for each species.
• Calculate Spatial Congruence: For each pair of species, compute the Spatial Congruence Index (CS) [22]: CS,ab = (Oab / Aa) * (Oab / Ab) where Oab is the area of overlap between species a and b, and Aa and Ab are the total areas of their respective ranges. The index ranges from 0 (no congruence) to 1 (perfect congruence).
• Set Congruence Threshold (CT): Define a quantitative congruence threshold (e.g., CS ≥ 0.75) as a parameter for identifying significant range overlap.
• Build Chorotype Network: Construct a network where species are vertices, and edges represent direct spatial congruences that meet the CT.
• Identify Indirect Congruences: Run an algorithm to find species linked through a chain of direct connections (e.g., ABC), forming "partial chorotypes" [22].
• Model Validation: Compare the chorotypes identified from your field-derived range maps against the chorotypes predicted by your spatial model. Discrepancies indicate areas where the model requires refinement.

The following diagram illustrates this analytical workflow:

Ecosystem Service Model Validation Protocol

Objective: To validate maps of predicted ecosystem services against independently collected field data.

Methodology:

• Run ES Models: Generate ES supply and use maps for your study area using one or more modeling platforms (e.g., InVEST, Co$ting Nature, WaterWorld) [82].
• Collect Independent Validation Data: This is a critical and often overlooked step [83]. Data must be independent, not used for model calibration. Examples include:
  • Water use: Municipal water withdrawal statistics or household survey data.
  • Carbon storage: Field measurements of above-ground biomass from plots.
  • Recreation: Geotagged social media photos or visitor counts.
• Perform Spatial Validation: Extract the modeled ES values at the locations of your validation data points and compute accuracy metrics (e.g., R², Root Mean Square Error, Spearman's ρ) [63] [82].
• Test Ensemble Performance: Where multiple models are available, create a simple median ensemble (taking the median value of all models for each pixel) and validate its performance. Research shows this can improve accuracy by 2-14% over individual models [63].
• Validate Realized Services: For ES that are used by people (realized ES), test if a simple model weighting potential ES by human population density outperforms a complex biophysical model [82].

The Scientist's Toolkit: Research Reagent Solutions

The following table details key datasets, models, and analytical tools essential for conducting robust model validation in this field.

Table 2: Essential Resources for Biodiversity and Ecosystem Service Model Validation

Tool / Resource	Type	Primary Function in Validation	Access / Reference
Spatial Congruence Analysis (SCAN)	Analytical Method	Identifies chorotypes (groups of species with congruent ranges) from field data for direct comparison with model outputs.	[22]
Global ES Model Ensembles	Data & Model	Provides pre-processed, median-ensemble maps for 5 key ES; serves as a benchmark for local models and fills data gaps.	[63]
InVEST Model Suite	Software Suite	Models multiple ES (e.g., carbon, water) for scenario analysis; model outputs can be validated with field data.	[82]
Co$ting Nature	Software Platform	Assesses ES supply, demand, and conservation priorities at a global scale; provides another model for ensemble validation.	[82]
WaterWorld	Software Model	A policy support tool for evaluating water-related ES; model outputs require validation with local discharge or water quality data.	[82]
Independent Validation Datasets	Data	Critical for testing model accuracy. Includes national statistics, field biomass plots, and household survey data.	[83] [82]

The validation of ecological models with field data is not merely an academic exercise but a fundamental step for credible conservation and restoration science. Evidence from post-mining ponds and large-scale ES studies confirms that with rigorous protocols—such as structured field sampling, Spatial Congruence Analysis, and ensemble modeling—researchers can effectively close the "certainty gap." This process reveals critical insights, such as the divergent recovery paths of biodiversity and ecosystem services and the dominant role of human demand in shaping realized service flows. By adopting these application notes and protocols, researchers can generate more reliable, validated evidence to support the sustainable management of recovering ecosystems worldwide.

In the contemporary landscape of corporate sustainability, impact assessments have evolved from a supplementary exercise to a core strategic function, essential for validating environmental sustainability (ES) metrics [84]. This transition is occurring within a broader business environment where, according to PwC’s 2024 Global Investor Survey, 76% of investors state that a company's ESG performance significantly influences their investment decisions [85]. The convergence of rigorous impact assessment with emerging research on biodiversity spatial organization offers a transformative opportunity to enhance the scientific validity and practical utility of corporate ecosystem service metrics. This protocol details the methodologies for integrating principles from biogeographical congruence analysis into corporate impact assessments, creating a new frontier for validating these critical metrics against established ecological patterns.

Key Concepts and Definitions

Corporate Impact Assessment: A systematic process for measuring, managing, and communicating the social and environmental effects of an organization's activities, extending beyond traditional financial metrics to demonstrate accountability and create strategic advantage [84].

ES Metrics: Quantifiable performance indicators that measure a company's environmental practices, including but not limited to greenhouse gas emissions, energy consumption, water usage, and waste management [86].

Biogeographical Core-Transition Organization: A recently identified universal spatial structure within biogeographical regions, forming ordered layers from core areas (characterized by high species richness and endemicity) to transition zones (characterized by high biota overlap and species occupancy) [9]. This structure transcends taxa and habitats, providing a template for validating corporate impact on ecosystems.

Double Materiality: An assessment approach, embedded in standards like the European Sustainability Reporting Standards (ESRS), that evaluates both the impact of a company on the environment and society, and the financial implications of sustainability issues on the company itself [86].

Quantitative ESG Metrics Framework

The foundation of robust corporate impact assessment lies in the precise definition and consistent tracking of key performance indicators. The following table synthesizes the core environmental sustainability metrics that companies should prioritize, aligned with major reporting frameworks.

Table 1: Key Environmental Sustainability Metrics for Corporate Impact Assessment

Metric Category	Specific Metrics	Reporting Frameworks	Data Requirements
Greenhouse Gas Emissions	Scope 1, 2, and 3 emissions (absolute and intensity-based) [86]	GHG Protocol, TCFD, ISSB S2 [85]	Fuel consumption, purchased energy data, supply chain engagement data
Energy Performance	Total energy consumption (MWh), Renewable energy percentage, Energy intensity [86]	GRI, SASB, TCFD [85]	Utility bills, power purchase agreements, generation records
Resource Management	Water withdrawal/consumption, Waste generation by type, Waste diversion rate [86]	GRI, ESRS [86]	Water invoices, waste hauling records, material flow tracking
Ecosystem Impact	Land use change, Impact on protected areas, Biodiversity action plans [85]	ESRS, GRI, CSRD-mandated disclosures [85]	Geospatial data, environmental impact assessments, supplier audits

Integrating Biogeographical Research into Impact Assessment

Recent ecological research reveals a universal core-to-transition organization of biodiversity across biogeographical regions [9]. This pattern, observed across terrestrial/marine vertebrates, invertebrates, and plants, demonstrates that regions comprise ordered layers from a core (high richness/endemicity) to transition zones (high biota overlap/species occupancy). This spatial congruence analysis provides a powerful, nature-based benchmark for validating corporate ES metrics.

Experimental Protocol: Spatial Congruence Analysis for ES Metric Validation

Objective: To validate corporate land and ecosystem impact metrics by testing their alignment with the predicted core-to-transition biogeographical organization.

Methodology:

• Geospatial Data Collection:

  • Obtain high-resolution species distribution data for the operational region (e.g., from GBIF or national biodiversity databases).
  • Acquire corporate geospatial data: facility boundaries, supply chain sourcing locations, and land management units.

• Biogeographical Sector Mapping:

  • Apply a k-means clustering algorithm to the region's biodiversity data based on four aspects: species richness, biota overlap, species occupancy, and endemicity [9].
  • This will identify the seven universal biogeographical sectors, mapping the core-to-transition gradient.

• Impact Metric Correlation:

  • Overlay corporate land-use and emission data onto the biogeographical sector map.
  • Statistically correlate the intensity of corporate impact (e.g., water usage, emission levels, land modification) with the position on the core-to-transition gradient.
  • Validation Criterion: Operations in core sectors should demonstrate significantly lower environmental impact intensities and stronger positive biodiversity outcomes than those in transition zones, aligning with the inherent ecological sensitivity of core areas.

Workflow Visualization:

Detailed Corporate Impact Assessment Protocol

This section provides a step-by-step experimental protocol for conducting a comprehensive corporate impact assessment that integrates biogeographical principles.

Phase 1: Materiality Assessment and Goal Setting

• Objective: Identify environmentally material issues and set SMART targets.
• Procedure:
  • Stakeholder Engagement: Conduct surveys and interviews with investors, employees, customers, and community representatives to understand their priorities [86].
  • Double Materiality Assessment: Evaluate both the impact of the company on the environment and the financial impact of environmental issues on the company, as required by frameworks like ESRS [86].
  • Target Setting: Establish Specific, Measurable, Achievable, Relevant, and Time-bound (SMART) targets. For example: "Reduce water consumption per unit of production by 15% by 2027" or "Achieve a 20% reduction in Scope 1 and 2 emissions intensity by 2026 from a 2023 baseline" [86].

Phase 2: Data Collection and Management

• Objective: Establish robust, automated systems for reliable ES data collection.
• Procedure:
  • System Implementation: Centralize data from multiple departments (operations, supply chain, HR) into a dedicated ESG software platform [85]. As noted by Gartner (2025), 80% of large firms now use such platforms to consolidate non-financial data [85].
  • Scope 3 Emissions Tracking: Engage with suppliers to collect data on indirect value chain emissions, which can account for over 70% of a company's carbon footprint [85].
  • Biodiversity Data Integration: Incorporate geospatial data on operational footprints and supply chain locations relative to sensitive ecological areas and the core-to-transition zones identified in Protocol 4.1.

Table 2: Research Reagent Solutions for Impact Assessment

Tool / Solution	Function	Application Example
ESG Data Platform	Automated consolidation and validation of non-financial data from disparate sources [85].	Tracking GHG emissions, energy use, and water consumption across global operations in real-time.
Geographic Information System (GIS)	Spatial analysis and mapping of corporate assets and supply chains against ecological data layers [9].	Overlaying facility locations with biogeographical sector maps to assess regional ecological impact.
IoT Sensors	Continuous, real-time monitoring of environmental conditions and resource flows [84].	Measuring water effluent quality or energy consumption at a production site.
AI-Powered Analytics	Detecting anomalies in ESG data, generating predictive insights, and automating GHG forecasting [85].	Predicting future emissions trends based on production forecasts and identifying data outliers for audit.

Phase 3: Impact Validation and Analysis

• Objective: Validate ES metrics against ecological benchmarks and analyze performance.
• Procedure:
  • Core-Transition Benchmarking: Execute the Spatial Congruence Analysis described in Protocol 4.1.
  • Performance Benchmarking: Compare ES metrics (e.g., carbon intensity, water efficiency) against industry peers using platforms like Bloomberg Terminal [85].
  • Outcome Measurement: Shift from measuring outputs (e.g., "trees planted") to outcomes (e.g., "improvement in ecosystem integrity measured by species richness") [84].

Phase 4: Assurance and Reporting

• Objective: Ensure data credibility and communicate performance transparently.
• Procedure:
  • Third-Party Assurance: Engage auditors (e.g., Big Four firms) for limited or reasonable assurance on ESG disclosures, as mandated by regulations like the EU's CSRD [85].
  • Integrated Reporting: Publish combined financial and ESG performance in a unified format, as demonstrated by companies like Nestlé and Apple [85].
  • Transparent Disclosure: Report balanced outcomes, including failures and unexpected outcomes, to refine strategy and build trust [87].

Integrated Assessment Workflow:

Discussion and Strategic Implications

Adopting this rigorous, science-based approach to corporate impact assessment delivers significant strategic advantages. Companies with strong ESG profiles supported by comprehensive impact reporting experience a lower cost of capital and greater operational resilience during market disruptions [84]. Furthermore, integrating biogeographical congruence analysis moves corporate sustainability from abstract metric tracking to tangible ecosystem stewardship, ensuring that environmental performance indicators are grounded in the fundamental patterns that govern life on Earth. This alignment not only mitigates regulatory and reputational risks but also positions companies as leaders in the purpose-driven economy of tomorrow [84].

Understanding the spatial relationships between biodiversity and ecosystem services (ES) is a central goal in sustainability science. This application note addresses a critical dimension of this research: quantifying how corporate activities from different industrial sectors impact ES and biodiversity, and the significant variation in these impacts both between and within sectors. Robust spatial congruence analysis reveals that economic sectors exert highly heterogeneous pressures on nature, necessitating sector-specific protocols for accurate impact assessment. These protocols enable researchers, policymakers, and corporate sustainability professionals to move beyond generic environmental assessments to spatially explicit, asset-level evaluations that account for the unique ways different industries interact with ecological systems. The frameworks outlined here support the broader thesis that effective environmental management depends on recognizing the distinct spatial fingerprints of industrial sectors on ES and biodiversity.

Quantitative Sectoral Impact Analysis

Comprehensive analysis of over 580,000 physical assets belonging to 2,173 global public companies reveals distinct sectoral patterns in ES and biodiversity impacts. The table below summarizes the average direct impacts (analogous to Scope 1 greenhouse gas emissions) for selected sectors, highlighting the variation between industrial categories [12].

Table 1: Sector-Level Impacts on Ecosystem Services and Biodiversity

Sector	Coastal Risk Reduction	Sediment Retention	Nitrogen Retention	Species Richness	Red List Species
Utilities	0.41	0.38	0.35	0.42	0.39
Real Estate	0.39	0.36	0.33	0.40	0.37
Materials	0.37	0.34	0.31	0.38	0.35
Financials	0.35	0.32	0.29	0.36	0.33
Consumer Discretionary	0.25	0.22	0.19	0.26	0.23

Note: Impact values represent indices normalized for comparison within metrics. Sectors are ranked by average total impact across metrics. Data sourced from [12].

Critically, substantial impact variation exists within sectors, often exceeding variation between sectors. For example, within the materials sector, individual company impacts can vary by an order of magnitude depending on the location and size of their physical assets [12]. This within-sector heterogeneity underscores the limitation of sector-average approaches and confirms the necessity of high-spatial-resolution, asset-level assessments for accurate impact quantification.

Experimental Protocols for Impact Assessment

Protocol 1: Corporate Asset-Level Impact Assessment

This protocol provides a standardized methodology for quantifying the direct ES and biodiversity impacts of corporate physical assets using open-source models and global datasets [12].

Research Reagent Solutions

Table 2: Essential Materials for Corporate Impact Assessment

Item	Function	Example Sources
Global Ecosystem Service Models	Quantifies baseline service provision (e.g., water purification, coastal protection)	InVEST, WaterWorld, Co$ting Nature
Biodiversity Indicator Datasets	Provides species richness, habitat for threatened/endemic species data	IUCN Red List, Key Biodiversity Areas, Map of Life
Corporate Asset Databases	Geospatial data on company physical assets (locations, types)	Truncated, Orbis, Bloomberg
Remote Sensing Imagery	Delineates precise asset footprints; monitors changes over time	Sentinel-2, Landsat, high-resolution commercial satellites
Potential Natural Vegetation Map	Establishes ecological baseline for impact counterfactual	MOSART, Simulated Potential Vegetation

Methodology

• Asset Footprint Delineation: For each corporate physical asset, establish the precise geographical footprint. For point-location data, apply size-adjusted buffering based on median asset size per activity type (e.g., mines vs. retail stores) [12].
• Baseline Establishment: Generate global maps of ES and biodiversity metrics under a potential natural vegetation scenario for four ES (coastal risk reduction, sediment retention, nitrogen retention, nature access) and four biodiversity metrics (species richness, Red List species habitat, endemic species habitat, Key Biodiversity Areas) [12].
• Impact Quantification: Assume complete ES/biodiversity value loss within asset footprints compared to the natural baseline. Calculate impact indices for each metric within each footprint [12].
• Aggregation and Normalization: Aggregate asset-level impacts to company and sector levels. Normalize impacts by revenue for cross-company comparability [12].
• Validation: Ground-truth impact assessments using high-resolution satellite imagery and, where accessible, field verification data [12].

Figure 1: Workflow for Corporate Asset-Level Impact Assessment

Protocol 2: Post-Industrial Site Assessment

This protocol assesses the recovery of ES and biodiversity in abandoned mining and quarrying sites, providing methodology for evaluating natural recovery versus active restoration outcomes [81].

Methodology

• Site Selection and Characterization: Select abandoned mining and quarry ponds across temporal gradients (e.g., 1960s, 1980s, 1990s). Record site type (quarry/mining), water presence (temporary/permanent), surface area, time since abandonment, and distance from urban/natural areas [81].
• Biodiversity Inventory: Conduct comprehensive species surveys across taxonomic groups. Record presence of animals, vascular plants, and habitat types using standardized plot or transect methods [81].
• Ecosystem Service Assessment: Document both services and disservices across 18 ES categories. Assess provisioning (e.g., water provision), regulating (e.g., water purification), and cultural services (e.g., recreation) [81].
• Index Development: Calculate two complementary indices:
  • Bioindex: Summarizes animal, plant, and habitat diversity.
  • Ecosystem Services Index (ESI): Quantifies net ES provision, accounting for disservices [81].
• Temporal Trend Analysis: Analyze how Bioindex and ESI values change with time since abandonment, comparing quarry versus mining ponds [81].

Essential Research Tools and Frameworks

Spatial Analysis Toolkit

Table 3: Research Reagent Solutions for Spatial Congruence Analysis

Tool/Model	Primary Function	Application Context
InVEST Suite	Spatially explicit ES modeling (carbon, water, habitat)	Regional ES assessment; trade-off analysis [88]
SolVES Model	Mapping social values of ES from survey data	Cultural service assessment; urban planning [33]
Geographically Weighted Regression (GWR)	Analyzing spatial non-stationarity in ES relationships	Identifying local trade-offs/synergies [88]
Self-Organizing Maps (SOM)	Identifying ES bundles (co-occurring services)	Regional zoning; spatial planning [88]
Random Forest Model	Quantifying driver contributions to ES patterns	Identifying key social-environmental factors [88]

Figure 2: Post-Industrial Site Assessment Workflow

Key Findings and Applications

Spatial congruence analysis reveals that biodiversity impacts are an insufficient proxy for ES impacts, as they exhibit distinct spatial patterns [12]. High-impact areas for multiple ES cover only 0.5% of land area, while high-impact biodiversity areas cover 3.2%, demonstrating significant spatial divergence [12].

The application of these protocols to the Beijing-Tianjin-Hebei urban agglomeration demonstrated distinct ES interactions, with four ES pairs showing synergies and six exhibiting trade-offs [88]. Factor analysis in agricultural watersheds has further identified three distinct ES groups: "forest and water synergies," "pasture and water synergies," and "crop and water quality tradeoffs" [89]. These relationships form the foundation for sector-specific management strategies that can enhance synergies while mitigating trade-offs.

For the mining sector, detailed lithium mine assessments demonstrate how high-resolution satellite imagery can track impact evolution over time, revealing substantial variation between individual mines and highlighting opportunities for impact reduction through improved siting and management [12]. In post-mining landscapes, research shows that both Bioindex and ESI increase significantly with time since abandonment, confirming natural recovery potential, though quarry ponds generally achieve higher ESI than mining ponds, suggesting a greater need for active restoration in the latter [81].

These protocols provide the scientific foundation for targeted corporate sustainability strategies, nature-related financial disclosure, and sector-specific environmental policy development, enabling stakeholders to accurately quantify and manage their unique impacts on Earth's life support systems.

Biodiversity as a Proxy for ES? Assessing Correlation and Complementarity

The relationship between biodiversity (BD) and ecosystem services (ES) represents a critical frontier in ecological research, with profound implications for conservation policy and sustainable development. This relationship is inherently spatial, yet the degree to which biodiversity patterns can serve as reliable proxies for ecosystem service provision remains a complex and context-dependent question. Understanding the correlation and complementarity between BD and ES requires frameworks that can simultaneously map both dimensions and analyze their spatial congruence.

The core-to-transition organization of biodiversity recently identified across biogeographical regions provides a fundamental spatial template for this analysis [9]. This organization reveals that biodiversity facets—including species richness, range size, endemicity, and biogeographical transitions—form ordered layers from core regions to transition zones, creating predictable spatial gradients that may correspond with ES provision. Simultaneously, research from karst forest ecosystems demonstrates that trade-offs and synergies between ES such as water yield, carbon storage, soil conservation, and biodiversity follow distinct spatial patterns influenced by environmental and anthropogenic drivers [90].

This Application Note provides detailed protocols for assessing BD-ES spatial relationships using emerging analytical frameworks, with particular emphasis on the Spatial Congruence Analysis (SCAN) method [22] and its application across different ecological contexts.

Theoretical Foundation: The Core-to-Transition Gradient

Recent research analyzing over 30,000 species across terrestrial and marine vertebrates, invertebrates, and plants has revealed a universal spatial organization of biodiversity within biogeographical regions [9]. This "core-to-transition" structure comprises seven distinct biogeographical sectors that form ordered layers reflecting consistent combinations of four fundamental biodiversity aspects:

• Species Richness: Number of species present
• Range Size: Geographical distribution of species
• Endemicity: Proportion of species unique to the region
• Biota Overlap: Degree of species mixing between adjacent regions

These sectors form gradients from regional hotspots (high richness, high endemicity, small range sizes) to transition zones (high biota overlap, lower endemicity, larger range sizes), with this organization driven by complementary environmental filters acting on species from regional hotspots and permeable biogeographical boundaries [9]. This fundamental structure provides a predictive framework for investigating how ES correlate with BD across different spatial contexts.

Table 1: Fundamental Biodiversity Aspects in Core-to-Transition Organization

Biodiversity Aspect	Core Regions	Transition Zones	Ecological Interpretation
Species Richness	High	Variable to Low	Favorable conditions for diversification and persistence
Species Range Size	Small	Large	Environmental filtering permits only tolerant species
Endemicity	High	Low	Historical isolation and unique evolutionary pressures
Biota Overlap	Low	High	Permeable boundaries allow species mixing

Quantitative Analysis of BD-ES Relationships

Documenting Trade-offs and Synergies

Research in the South China Karst forests provides compelling quantitative evidence of the complex relationships between multiple ecosystem services and biodiversity. This study employed integrated modeling approaches (InVEST and RUSLE) across a 20-year period (2000-2020) to track four key services: water yield (WY), carbon storage (CS), soil conservation (SC), and biodiversity (Bio) [90].

Table 2: Ecosystem Service Trade-offs and Synergies in Karst Forests (2000-2020)

Ecosystem Service	Net Change (%)	Relationship with BD	Primary Driver Influence
Water Yield (WY)	+13.44	Trade-off	Positive: PrecipitationNegative: Population density
Soil Conservation (SC)	+4.94	Synergy (weak)	Positive: PrecipitationNegative: Population density
Carbon Storage (CS)	-0.03	Synergy	Positive: Temperature, Vegetation CoverNegative: Population density
Biodiversity (Bio)	-0.61	-	Positive: Temperature, Vegetation CoverNegative: Population density

The analysis revealed a predominant trade-off relationship between provisioning services (water yield) and regulating/supporting services (carbon storage, biodiversity), with biodiversity showing the most consistent decline. Spatial analysis further demonstrated that these relationships exhibited significant heterogeneity across different karst landforms, with losses concentrated in karst gorges, fault basins, and middle-high mountains [90].

Relational Values and Policy Mismatch

Expert surveys from marine systems (Wadden Sea and Algoa Bay) reveal a crucial dimension of BD-ES relationships: biodiversity changes disproportionately impact non-material NCPs (Nature's Contributions to People) compared to material and regulating categories [91]. Specifically, experts identified significant impacts on:

• Experiences (healing, recreation, aesthetic enjoyment)
• Learning and inspiration
• Cultural and spiritual identities

This finding highlights a critical policy mismatch, as conservation strategies often prioritize material and regulating services while undervaluing these relational dimensions. When using biodiversity as a proxy for ES, this suggests assessment frameworks must incorporate non-material contributions that are often overlooked in traditional analyses [91].

Experimental Protocols: Spatial Congruence Analysis (SCAN)

The Spatial Congruence Analysis (SCAN) method provides a robust framework for quantifying spatial congruence among species distributions and can be adapted to analyze BD-ES relationships [22]. SCAN operates on vector-based distribution data (polygons), avoiding scale-dependent distortions associated with grid-based approaches, and explicitly uses congruence thresholds as a central parameter.

Key Innovations of SCAN:

• Uses direct range map comparisons rather than grid-cells
• Identifies both direct and indirect spatial relationships
• Explores patterns across a continuum of congruence values
• Reveals complex chorotypes from simple to highly complex

Step-by-Step Protocol

Step 1: Data Preparation and Congruence Index Calculation

• Format species range data or ES provision maps as polygon shapefiles
• Calculate pairwise Spatial Congruence Index (CS) values for all combinations using:

CS,ab = (Oab/Aa) × (Oab/Ab)

Where: Oab = area of overlap between a and b; Aa and Ab = areas of ranges a and b respectively [22]

Step 2: Network Construction and Threshold Setting

• Construct a one-layer network connecting entities (species/ES) using pairwise CS values as edges
• Set initial Congruence Threshold (CT) to maximum stringency (e.g., 100%)
• Systematically decrease CT in increments (e.g., 1%) for subsequent analyses

Step 3: Chorotype Identification

• For each reference entity at a given CT, identify all directly congruent partners (CS ≥ CT)
• Expand group to include indirectly congruent entities connected through chains of direct relationships
• Iterate until no new entities are added ("closed list"), forming a "partial chorotype"
• Repeat across all CT values to build comprehensive chorotype profiles [22]

Step 4: BD-ES Congruence Assessment

• Apply SCAN separately to biodiversity data and ES provision maps
• Compare resulting chorotypes to identify areas of spatial congruence and divergence
• Quantify congruence strength using CS values at multiple thresholds
• Analyze environmental correlates of congruent vs. non-congruent areas

Application to BD-ES Research

When applying SCAN to BD-ES relationships, researchers can:

• Identify Spatial Correlations: Detect areas where biodiversity hotspots coincide with ES hotspots
• Quantify Complementarity: Measure the degree to which BD and ES provide complementary versus redundant spatial information
• Detect Scale Dependencies: Analyze how BD-ES relationships change across spatial scales and congruence thresholds
• Inform Conservation Planning: Identify areas where biodiversity conservation simultaneously maximizes ES provision

The method's capacity to handle both direct and indirect relationships makes it particularly valuable for capturing the complex networks of interactions that characterize BD-ES relationships in heterogeneous landscapes [22].

The Scientist's Toolkit: Research Reagent Solutions

Table 3: Essential Analytical Tools for BD-ES Spatial Congruence Research

Research Tool	Function	Application Context
SCAN Algorithm [22]	Identifies chorotypes through direct/indirect spatial relationships	Quantifying spatial congruence between BD and ES distributions
Infomap Bioregion [9]	Detects biogeographical regions using network analysis	Establishing core-to-transition biodiversity framework
InVEST Model [90]	Quantifies multiple ecosystem services	Modeling water yield, carbon storage, habitat quality
RUSLE Model [90]	Estimates soil conservation service	Calculating soil erosion prevention capacity
Spectral Variation Analysis [92]	Remote sensing of plant diversity	Estimating α- and β-diversity from hyperspectral imagery
Geodetector [90]	Identifies driving factors	Analyzing environmental drivers of BD-ES relationships

Integrated Analysis Framework

The integration of these approaches enables a comprehensive assessment of BD-ES relationships across spatial scales. The core-to-transition organization provides the foundational biodiversity gradient [9], while SCAN enables precise quantification of spatial congruence [22], and integrated modeling (InVEST/RUSLE) quantifies ES provision and trade-offs [90]. This multi-framework approach addresses both correlation (statistical relationships) and complementarity (unique information content) dimensions of the BD-ES proxy question.

The question "Biodiversity as a Proxy for ES?" necessitates a nuanced, context-dependent answer. Evidence suggests that:

• Spatial Congruence is Partial: Biodiversity patterns show significant but incomplete spatial correlation with ES provision, varying by service type and biogeographical context [9] [90]
• Core-Transition Gradients Matter: BD-ES relationships differ systematically between core biodiversity regions and transition zones [9]
• Trade-offs are Pervasive: Optimizing for single ES often creates trade-offs with biodiversity and other services [90]
• Non-Material Values are Undervalued: Biodiversity serves as a particularly important proxy for non-material NCPs that are frequently overlooked in policy [91]

Future research should prioritize: (1) multi-taxon comparisons of BD-ES relationships across biogeographical sectors; (2) dynamic analyses of how BD-ES correlations shift under environmental change; and (3) development of integrated indicators that capture both BD and ES dimensions for conservation planning and policy implementation.

Conclusion

The analysis of spatial congruence between biodiversity and ecosystem services reveals a complex and context-dependent picture. There is no universal, high correlation, necessitating a service-specific and location-aware approach to conservation planning. Methodological advancements now allow for sophisticated spatial prioritization, yet significant challenges remain, including the need for appropriate proxies, a deeper mechanistic understanding of B-ES relationships, and the integration of cost and feasibility. For biomedical researchers and drug development professionals, these findings underscore a critical implication: the conservation of biodiversity-rich areas is not merely an ecological goal but a direct investment in the natural chemical libraries from which future medicines are derived. The degradation of ecosystems equates to the irreversible loss of genetic and biochemical diversity, potentially extinguishing compounds before discovery. Future research must focus on quantifying these direct links, identifying 'win-win' areas that are crucial for both ecological integrity and biomedical prospecting, and developing policies that incentivize the conservation of nature's pharmacy through robust payments for ecosystem services schemes.

References

• 1. Editorial: Biodiversity, ecosystem functions and services [https://www.frontiersin.org/journals/ecology-and-evolution/articles/10.3389/fevo.2022.1086408/full]
• 2. Biodiversity [https://www.who.int/news-room/fact-sheets/detail/biodiversity]
• 3. Linking Biodiversity, Ecosystem Functioning and Services ... [https://www.sciencedirect.com/science/article/pii/S0065250415000288]
• 4. The quest for a mechanistic understanding of biodiversity ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC4633867/]
• 5. Integrating ecosystem services within spatial biodiversity ... [https://www.sciencedirect.com/science/article/pii/S2212041620301285]
• 6. Spatial patterns and interrelationships between biodiversity ... [https://www.biodiversity-science.net/EN/10.17520/biods.2023061]
• 7. EnviroAtlas Benefit Category: Biodiversity Conservation [https://www.epa.gov/enviroatlas/enviroatlas-benefit-category-biodiversity-conservation]
• 8. Biodiversity and ecosystem services: a multilayered ... [https://www.sciencedirect.com/science/article/pii/S0169534711002424]
• 9. A general rule on the organization of biodiversity in Earth's ... [https://www.nature.com/articles/s41559-025-02724-5]
• 10. Do biodiversity-ecosystem functioning experiments inform ... [https://www.sciencedirect.com/science/article/pii/S0006320720301865]
• 11. Spatial Congruence Analysis (SCAN): A method for detecting ... [https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0245818]
• 12. An open-source approach for measuring corporate impacts ... [https://www.nature.com/articles/s43247-024-01797-7]
• 13. Spatial congruence between biodiversity and ecosystem ... [https://www.sciencedirect.com/science/article/abs/pii/S0006320708004436]
• 14. A Case Study in Baiyangdian Basin, China [https://www.mdpi.com/2073-445X/13/10/1619]
• 15. Identifying ecosystem service hotspots for environmental ... [http://www.scielo.org.za/scielo.php?script=sci_arttext&pid=S0006-82412016000200011]
• 16. Spatial prioritisation for conserving ecosystem services [https://pmc.ncbi.nlm.nih.gov/articles/PMC4722056/]
• 17. Spatial multi-scale relationships of ecosystem services [https://www.nature.com/articles/s41598-017-09863-1]
• 18. Linkages between biodiversity attributes and ecosystem ... [https://www.sciencedirect.com/science/article/pii/S2212041614000576]
• 19. Trends and gaps in biodiversity and ecosystem services ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC9666618/]
• 20. Partitioning Net Biodiversity Effects on Ecosystem ... [https://pubmed.ncbi.nlm.nih.gov/41206877/]
• 21. Mangrove biodiversity and ecosystem services [https://www.nature.com/articles/s44358-025-00103-3]
• 22. Spatial Congruence Analysis (SCAN): A method for detecting ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC8136640/]
• 23. Specifications for cross-scale inclusion of harmonised ... [https://www.biodiversa.eu/2025/10/10/specifications-for-cross-scale-inclusion-of-harmonised-protocols/]
• 24. Identifying spatial priorities for protecting ecosystem services [https://pmc.ncbi.nlm.nih.gov/articles/PMC3814920/]
• 25. Spatial analysis of multiple ecosystem service flows from ... [https://research.utwente.nl/en/publications/spatial-analysis-of-multiple-ecosystem-service-flows-from-natural/]
• 26. Analysis of integrated studies on ecosystem service ... [https://link.springer.com/article/10.1007/s11252-025-01831-2]
• 27. Analysis of integrated studies on ecosystem service ... [https://ui.adsabs.harvard.edu/abs/2025UrbEc..28..214U/abstract]
• 28. Global insights on spatial and temporal assessments of ... [https://link.springer.com/article/10.1007/s10668-025-06797-x]
• 29. Study on the spatial and temporal evolution of ecosystem ... [https://www.nature.com/articles/s41598-024-83257-y]
• 30. Assessing ecosystem service value dynamics and ... [https://www.frontiersin.org/journals/environmental-science/articles/10.3389/fenvs.2025.1648758/full]
• 31. A comprehensive approach for the spatial optimization of ... [https://link.springer.com/article/10.1007/s10980-024-02026-4]
• 32. Opportunity cost estimates for spatial conservation ... [https://www.nature.com/articles/s41597-025-05837-5]
• 33. Social valuation of urban ecosystem services using the ... [https://www.nature.com/articles/s41599-025-05946-x]
• 34. Is Rich and Rare the Common Share? Describing Biodiversity ... [https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0056073]
• 35. Multiple macroevolutionary routes to becoming a ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC6365113/]
• 36. Spatial Distribution Patterns, Environmental Drivers, and ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC11939462/]
• 37. Simulated biodiversity hotspots from traditional ecological ... [https://www.nature.com/articles/s43247-025-02309-x]
• 38. Systematic planning shows more than half of the most ... [https://www.sciencedirect.com/science/article/pii/S2351989424002403]
• 39. Spatiotemporal dynamics and driving factors of ecosystem ... [https://www.nature.com/articles/s41598-024-76838-4]
• 40. All the OBIS-related sessions at Living Data 2025 at a glance! [https://obis.org/2025/04/24/living-data-call-for-abstracts-obis-sessions/]
• 41. Multi-Criteria Decision-Making Scenario Insights into ... [https://www.mdpi.com/2073-445X/13/11/1964]
• 42. A Hybrid Multi-Criteria Decision-Making Framework for the ... [https://www.mdpi.com/2078-2489/16/6/454]
• 43. Designing an innovative Multi-Criteria Decision Making ... [https://www.nature.com/articles/s41598-025-14283-7]
• 44. action agendas need spatial data to integrate climate and ... [https://www.nature.com/articles/s44168-025-00307-5]
• 45. Highlights from the Biodiversity Monitoring Science Fair [https://www.biodiversa.eu/2025/06/19/event-connecting-people-technology-and-data-for-biodiversity-monitoring/]
• 46. a comparative assessment of ecosystem services models ... [https://www.nature.com/articles/s41598-024-76600-w]
• 47. InVEST - Stanford Natural Capital Project [https://naturalcapitalproject.stanford.edu/software/invest]
• 48. Assessment and multi-scenario prediction of ecosystem ... [https://www.frontiersin.org/journals/ecology-and-evolution/articles/10.3389/fevo.2025.1539547/full]
• 49. Progress in regulating ecosystem services research [https://www.nature.com/articles/s40494-025-01674-2]
• 50. Systematic conservation planning for people and nature [https://www.sciencedirect.com/science/article/pii/S2212041624000445]
• 51. Ecosystem services and connectivity in spatial ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC7154789/]
• 52. UN Biodiversity Lab – Providing decision makers with the best ... [https://unbiodiversitylab.org/en/]
• 53. Spatial Biodiversity Modeling [https://bgc.yale.edu/spatial-biodiversity-modeling]
• 54. Marxan: Home [https://marxansolutions.org/]
• 55. Software [https://marxansolutions.org/software/]
• 56. MARXAN Platform Changes the Nature of Conservation ... [https://www.nature.org/en-us/what-we-do/our-insights/perspectives/marxan-microsoft-conservation-planning-platform/]
• 57. The pitfalls of biodiversity proxies: Differences in richness ... [https://www.nature.com/articles/s41598-019-55490-3]
• 58. Understanding the effects of spatial scaling on ... [https://www.sciencedirect.com/science/article/pii/S1569843225000883]
• 59. How to choose colors for data visualizations [https://www.atlassian.com/data/charts/how-to-choose-colors-data-visualization]
• 60. Data Visualization: Design Considerations [https://guides.lib.berkeley.edu/data-visualization/design]
• 61. Embracing scale-dependence to achieve a deeper ... [https://pubmed.ncbi.nlm.nih.gov/30182500/]
• 62. Scale-dependence of resource-biodiversity relationships [https://www.sciencedirect.com/science/article/abs/pii/S002251930300239X]
• 63. Model ensembles of ecosystem services fill global certainty ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC10081842/]
• 64. Hotspots of biodiversity and ecosystem services [https://oneecosystem.pensoft.net/article/24719/]
• 65. A high spatial resolution dataset of ecosystem services ... [https://essd.copernicus.org/preprints/essd-2025-107/]
• 66. Spatial assessment of ecosystem services (phase 1) [https://www.peer.eu/publications/spatial-assessment-of-ecosystem-services-phase-1]
• 67. Changes in the location of biodiversity–ecosystem function ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC5394661/]
• 68. Biodiversity Hotspots Defined | CEPF [https://www.cepf.net/our-work/biodiversity-hotspots/hotspots-defined]
• 69. What Are Biodiversity Hotspots? [https://www.conservation.org/learning/biodiversity-hotspots]
• 70. What are biodiversity hotspots? [https://www.ifaw.org/journal/what-are-biodiversity-hotspots]
• 71. Global conservation priorities for wetlands and setting post ... [https://www.nature.com/articles/s43247-023-01195-5]
• 72. Integrative analysis of biodiversity, ecosystem services ... [https://www.sciencedirect.com/science/article/abs/pii/S0048969723017151]
• 73. Identifying Priority Areas for Biodiversity and Soil... [https://reference-global.com/article/10.2478/jlecol-2025-0030]
• 74. Cost-effective conservation of multiple goals in grasslands [https://www.sciencedirect.com/science/article/pii/S0006320725004380]
• 75. Species richness prediction and priority conservation ... [https://www.nature.com/articles/s41598-025-11025-7]
• 76. Key questions for understanding drivers of biodiversity ... [https://link.springer.com/article/10.1007/s10980-024-01842-y]
• 77. An overview of methods to enhance ecosystem services ... [https://aries.integratedmodelling.org/aries-hub/an-overview-of-methods-to-enhance-ecosystem-services-assessments-within-the-global-biodiversity-framework/]
• 78. A 45-year global analysis of the spatial human forest nexus [https://www.nature.com/articles/s43247-025-02514-8]
• 79. A Framework for Assessing Trade-Offs and Synergies in ... [https://www.mdpi.com/1999-4907/14/8/1614]
• 80. Identifying Conservation Priority Areas of Hydrological ... [https://www.mdpi.com/2072-4292/16/18/3409]
• 81. Changes in biodiversity and ecosystem services over time ... [https://link.springer.com/article/10.1007/s10750-025-05892-4]
• 82. A Continental-Scale Validation of Ecosystem Service Models [https://link.springer.com/article/10.1007/s10021-019-00380-y]
• 83. Ecosystem services mapping and modelling. Where is the ... [https://www.sciencedirect.com/science/article/pii/S2666683925000252]
• 84. The Strategic Value of Impact Reporting: Why It Matters ... [https://www.givingforce.com/blog/the-strategic-value-of-impact-reporting-why-it-matters-more-than-ever-in-2025]
• 85. How Do Companies Measure Their ESG Performance in ... [https://senecaesg.com/insights/how-do-companies-measure-their-esg-performance-in-2025/]
• 86. 7-Key-ESG-Metrics-Every-Company-Should-Track-2025 [https://www.fiegenbaum.solutions/en/blog/7-key-esg-metrics-2025]
• 87. The impact investment sector's open secret is spotty data [https://www.weforum.org/stories/2025/10/measuring-impact-investing/]
• 88. Urban Spatial Management and Planning Based on the ... [https://www.mdpi.com/2072-4292/17/7/1258]
• 89. Spatial interactions among ecosystem services in an ... [https://pmc.ncbi.nlm.nih.gov/articles/PMC3718150/]
• 90. Trade-offs and synergies of ecosystem services and their ... [https://www.nature.com/articles/s40494-025-01847-z]
• 91. Marine biodiversity change impacts relational values [https://www.nature.com/articles/s44183-025-00148-z]
• 92. Focus on Spatial Analysis of Biodiversity Patterns [https://iopscience.iop.org/collections/multi-240214-01]