The ongoing global change and the increased interest in macroecological processes call for the analysis of spatially extensive data on species communities to understand and forecast distributional changes of biodiversity. Recently developed joint species distribution models can deal with numerous species efficiently, while explicitly accounting for spatial structure in the data. However, their applicability is generally limited to relatively small spatial data sets because of their severe computational scaling as the number of spatial locations increases. In this work, we propose a practical alleviation of this scalability constraint for joint species modeling by exploiting two spatial‐statistics techniques that facilitate the analysis of large spatial data sets: Gaussian predictive process and nearest‐neighbor Gaussian process. We devised an efficient Gibbs posterior sampling algorithm for Bayesian model fitting that allows us to analyze community data sets consisting of hundreds of species sampled from up to hundreds of thousands of spatial units. The performance of these methods is demonstrated using an extensive plant data set of 30,955 spatial units as a case study. We provide an implementation of the presented methods as an extension to the hierarchical modeling of species communities framework.

Understanding the role of interspecific interactions in shaping ecological communities is one of the central goals in community ecology. In fungal communities, measuring interspecific interactions directly is challenging because these communities are composed of large numbers of species, many of which are unculturable. An indirect way of assessing the role of interspecific interactions in determining community structure is to identify the species co‐occurrences that are not constrained by environmental conditions. In this study, we investigated co‐occurrences among root‐associated fungi, asking whether fungi co‐occur more or less strongly than expected based on the environmental conditions and the host plant species examined. We generated molecular data on root‐associated fungi of five plant species evenly sampled along an elevational gradient at a high arctic site. We analysed the data using a joint species distribution modelling approach that allowed us to identify those co‐occurrences that could be explained by the environmental conditions and the host plant species, as well as those co‐occurrences that remained unexplained and thus more probably reflect interactive associations. Our results indicate that not only negative but also positive interactions play an important role in shaping microbial communities in arctic plant roots. In particular, we found that mycorrhizal fungi are especially prone to positively co‐occur with other fungal species. Our results bring new understanding to the structure of arctic interaction networks by suggesting that interactions among root‐associated fungi are predominantly positive.

Joint species distribution modelling (JSDM) is a widely used statistical method that analyzes combined patterns of all species in a community, linking empirical data to ecological theory and enhancing community-wide prediction tasks. However, fitting JSDMs to large datasets is often computationally demanding and time-consuming. Recent studies have introduced new statistical and machine learning techniques to provide more scalable fitting algorithms, but extending these to complex JSDM structures that account for spatial dependencies or multi-level sampling designs remains challenging. In this study, we aim to enhance JSDM scalability by leveraging high-performance computing (HPC) resources for an existing fitting method. Our work focuses on the Hmsc R-package, a widely used JSDM framework that supports the integration of various dataset types into a single comprehensive model. We developed a GPU-compatible implementation of its model-fitting algorithm using Python and the TensorFlow library. Despite these changes, our enhanced framework retains the original user interface of the Hmsc R-package. We evaluated the performance of the proposed implementation across various model configurations and dataset sizes. Our results show a significant increase in model fitting speed for most models compared to the baseline Hmsc R-package. For the largest datasets, we achieved speed-ups of over 1000 times, demonstrating the substantial potential of GPU porting for previously CPU-bound JSDM software. This advancement opens promising opportunities for better utilizing the rapidly accumulating new biodiversity data resources for inference and prediction.

How community‐level specialization differs among groups of organisms, and changes along environmental gradients, is fundamental to understanding the mechanisms influencing ecological communities. In this paper, we investigate the specialization of root‐associated fungi for plant species, asking whether the level of specialization varies with elevation. For this, we applied DNA barcoding based on the ITS region to root samples of five plant species equivalently sampled along an elevational gradient at a high arctic site. To assess whether the level of specialization changed with elevation and whether the observed patterns varied between mycorrhizal and endophytic fungi, we applied a joint species distribution modeling approach. Our results show that host plant specialization is not environmentally constrained in arctic root‐associated fungal communities, since there was no evidence for changing specialization with elevation, even if the composition of root‐associated fungal communities changed substantially. However, the level of specialization for particular plant species differed among fungal groups, root‐associated endophytic fungal communities being highly specialized on particular host species, and mycorrhizal fungi showing almost no signs of specialization. Our results suggest that plant identity affects associated mycorrhizal and endophytic fungi differently, highlighting the need of considering both endophytic and mycorrhizal fungi when studying specialization in root‐associated fungal communities.

Remote monitoring of natural afforestation processes on abandoned agricultural lands is crucial for assessments and predictions of forest cover dynamics, biodiversity, ecosystem functions and services. In this work, we built on the general approach of combining satellite and field data for forest mapping and developed a simple and robust method for afforestation dynamics assessment. This method is based on Landsat imagery and index-based thresholding and specifically targets suitability for limited field data. We demonstrated method’s details and performance by conducting a case study for two bordering districts of Rudnya (Smolensk region, Russia) and Liozno (Vitebsk region, Belarus). This study area was selected because of the striking differences in the development of the agrarian sectors of these countries during the post-Soviet period (1991-present day). We used Landsat data to generate a consistent time series of five-year cloud-free multispectral composite images for the 1985–2020 period via the Google Earth Engine. Three spectral indices, each specifically designed for either forest, water or bare soil identification, were used for forest cover and arable land mapping. Threshold values for indices classification were both determined and verified based on field data and additional samples obtained by visual interpretation of very high-resolution satellite imagery. The developed approach was applied over the full Landsat time series to quantify 35-year afforestation dynamics over the study area. About 32% of initial arable lands and grasslands in the Russian district were afforested by the end of considered period, while the agricultural lands in Belarus’ district decreased only by around 5%. Obtained results are in the good agreement with the previous studies dedicated to the agricultural lands abandonment in the Eastern Europe region. The proposed method could be further developed into a general universally applicable technique for forest cover mapping in different growing conditions at local and regional spatial levels.