Explore the vast range of titles available.

The detection of causal interactions is of great importance when inferring complex ecosystem functional and structural networks for basic and applied research. Convergent cross mapping (CCM) based on nonlinear state-space reconstruction made substantial progress about network inference by measuring how well historical values of one variable can reliably estimate states of other variables. Here we investigate the ability of a developed optimal information flow (OIF) ecosystem model to infer bidirectional causality and compare that to CCM. Results from synthetic datasets generated by a simple predator-prey model, data of a real-world sardine-anchovy-temperature system and of a multispecies fish ecosystem highlight that the proposed OIF performs better than CCM to predict population and community patterns. Specifically, OIF provides a larger gradient of inferred interactions, higher point-value accuracy and smaller fluctuations of interactions and α -diversity including their characteristic time delays. We propose an optimal threshold on inferred interactions that maximize accuracy in predicting fluctuations of effective α -diversity, defined as the count of model-inferred interacting species. Overall OIF outperforms all other models in assessing predictive causality (also in terms of computational complexity) due to the explicit consideration of synchronization, divergence and diversity of events that define model sensitivity, uncertainty and complexity. Thus, OIF offers a broad ecological information by extracting predictive causal networks of complex ecosystems from time-series data in the space-time continuum. The accurate inference of species interactions at any biological scale of organization is highly valuable because it allows to predict biodiversity changes, for instance as a function of climate and other anthropogenic stressors. This has practical implications for defining optimal ecosystem management and design, such as fish stock prioritization and delineation of marine protected areas based on derived collective multispecies assembly. OIF can be applied to any complex system and used for model evaluation and design where causality should be considered as non-linear predictability of diverse events of populations or communities.

Non-pharmaceutical interventions (NPIs) including resource allocation, risk communication, social distancing and travel restriction, are mainstream actions to control the spreading of Coronavirus disease 2019 (COVID-19) worldwide. Different countries implemented their own combinations of NPIs to prevent local epidemics and healthcare system overloaded. Portfolios, as temporal sets of NPIs have various systemic impacts on preventing cases in populations. Here, we developed a probabilistic modeling framework to evaluate the effectiveness of NPI portfolios at the macroscale. We employed a deconvolution method to back-calculate incidence of infections and estimate the effective reproduction number by using the package EpiEstim. We then evaluated the effectiveness of NPIs using ratios of the reproduction numbers and considered them individually and as a portfolio systemically. Based on estimates from Japan, we estimated time delays of symptomatic-to-confirmation and infection-to-confirmation as 7.4 and 11.4 days, respectively. These were used to correct surveillance data of other countries. Considering 50 countries, risk communication and returning to normal life were the most and least effective yielding the aggregated effectiveness of 0.11 and − 0.05 that correspond to a 22.4% and 12.2% reduction and increase in case growth. The latter is quantified by the change in reproduction number before and after intervention implementation. Countries with the optimal NPI portfolio are along an empirical Pareto frontier where mean and variance of effectiveness are maximized and minimized independently of incidence levels. Results indicate that implemented interventions, regardless of NPI portfolios, had distinct incidence reductions and a clear timing effect on infection dynamics measured by sequences of reproduction numbers. Overall, the successful suppression of the epidemic cannot work without the non-linear effect of NPI portfolios whose effectiveness optimality may relate to country-specific socio-environmental factors.

Fish ecosystems perform ecological functions that are critically important for the sustainability of marine ecosystems, such as global food security and carbon stock. During the 21st century, significant global warming caused by climate change has created pressing challenges for fish ecosystems that threaten species existence and global ecosystem health. Here, we study a coastal fish community in Maizuru Bay, Japan, and investigate the relationships between fluctuations of ST, abundance-based species interactions and salient fish biodiversity. Observations show that a local 20% increase in temperature from 2002 to 2014 underpins a long-term reduction in fish diversity (∼25%) played out by some native and invasive species (e.g. Chinese wrasse) becoming exceedingly abundant; this causes a large decay in commercially valuable species (e.g. Japanese anchovy) coupled to an increase in ecological productivity. The fish community is analyzed considering five temperature ranges to understand its atemporal seasonal sensitivity to ST changes, and long-term trends. An optimal information flow model is used to reconstruct species interaction networks that emerge as topologically different for distinct temperature ranges and species dynamics. Networks for low temperatures are more scale-free compared to ones for intermediate (15-20°C) temperatures in which the fish ecosystem experiences a first-order phase transition in interactions from locally stable to metastable and globally unstable for high temperatures states as suggested by abundance-spectrum transitions. The dynamic dominant eigenvalue of species interactions shows increasing instability for competitive species (spiking in summer due to intermediate-season critical transitions) leading to enhanced community variability and critical slowing down despite higher time-point resilience. Native competitive species whose abundance is distributed more exponentially have the highest total directed interactions and are keystone species (e.g. Wrasse and Horse mackerel ) for the most salient links with cooperative decaying species. Competitive species, with higher eco-climatic memory and synchronization, are the most affected by temperature and play an important role in maintaining fish ecosystem stability via multitrophic cascades (via cooperative-competitive species imbalance), and as bioindicators of change. More climate-fitted species follow temperature increase causing larger divergence divergence between competitive and cooperative species. Decreasing dominant eigenvalues and lower relative network optimality for warmer oceans indicate fishery more attracted toward persistent oscillatory states, yet unpredictable, with lower cooperation, diversity and fish stock despite the increase in community abundance due to non-commercial and venomous species. We emphasize how changes in species interaction organization, primarily affected by temperature fluctuations, are the backbone of biodiversity dynamics and yet for functional diversity in contrast to taxonomic richness. Abundance and richness manifest gradual shifts while interactions show sudden shift. The work provides data-driven tools for analyzing and monitoring fish ecosystems under the pressure of global warming or other stressors. Abundance and interaction patterns derived by network-based analyses proved useful to assess ecosystem susceptibility and effective change, and formulate predictive dynamic information for science-based fishery policy aimed to maintain marine ecosystems stable and sustainable.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

How can the shape of biodiversity inform us about cities’ ecoclimatic fitness and guide their development? Can we use species as the harbingers of climatic extremes? Eco-climatically sensitive species carry information about hydroclimatic change in their distribution, fitness, and preferential gradients of habitat suitability. Conversely, environmental features outside of the species’ fitness convey information on potential ecological anomalies in response to extremes to adapt or mitigate, such as through urban parks. Here, to quantify ecosystems’ fitness, we propose a novel computational model to extract multivariate functional ecological networks and their basins, which carry the distributed signature of the compounding hydroclimatic pressures on sentinel species. Specifically, we consider butterflies and their habitat suitability (HS) to infer maximum suitability gradients that are meaningful of potential species networks and flows, with the smallest hydroclimatic resistance across urban landscapes. These flows are compared to the distribution of urban parks to identify parks’ ecological attractiveness, actual and potential connectivity, and park potential to reduce hydroclimatic impacts. The ecosystem fitness index (EFI) is novelly introduced by combining HS and the divergence of the relative species abundance (RSA) from the optimal log-normal Preston plot. In Shenzhen, as a case study, eco-flow networks are found to be spatially very extended, scale-free, and clustering for low HS gradient and EFI areas, where large water bodies act as sources of ecological corridors draining into urban parks. Conversely, parks with higher HS, HS gradients, and EFIs have small-world connectivity non-overlapping with hydrological networks. Diverging patterns of abundance and richness are inferred as increasing and decreasing with HS. HS is largely determined by temperature and precipitation of the coldest quarter and seasonality, which are critical hydrologic variables. Interestingly, a U-shape pattern is found between abundance and diversity, similar to the one in natural ecosystems. Additionally, both abundance and richness are mildly associated with park area according to a power function, unrelated to longitude but linked to the degree of urbanization or park centrality, counterintuitively. The Preston plot’s richness–abundance and abundance-rank patterns were verified to reflect the stationarity or ecological meta-equilibrium with the environment, where both are a reflection of community connectivity. Ecological fitness is grounded on the ecohydrological structure and flows where maximum HS gradients are indicative of the largest eco-changes like climate-driven species flows. These flows, as distributed stress-response functions, inform about the collective eco-fitness of communities, like parks in cities. Flow-based networks can serve as blueprints for designing ecotones that regulate key ecosystem functions, such as temperature and evapotranspiration, while generating cascading ecological benefits across scales. The proposed model, novelly infers HS eco-networks and calculates the EFI, is adaptable to diverse sensitive species and environmental layers, offering a robust tool for precise ecosystem assessment and design.

In natural ecosystems, large-scale structure and function arise from the spatial coordination of local habitat shaped by environmental pressures. Replicating such multiscale organization in synthetic systems—where emergent form and function are governed solely by material composition—remains a fundamental challenge in design and engineering. Here, we developed an experimental, analytical, and computational framework to create self-organizing material patterns—Terraforms—using a ternary mixture of clay, water, and binder subjected to uniform mechanical pressure. These networked patterns form instantaneously at material interfaces and exhibit 3D scale-free organization, positioned between random and regular configurations, mimicking diverse ecological patterns, with planar symmetry and microtopographic fractality. Systematic variation of relative material proportions (RMPs) revealed a critical regime—50% clay, 15–20% water, and 30–35% binder—with a water-to-binder (W/B) ratio close to the Pacioli golden ratio (1.62), yielding Pareto-optimal patterns. Within this regime, network structures minimize eco-stress, maximize hygroscopic capacity, and enhance systemic stability, thereby enabling their quantification. Transitions between pattern topologies follow gradients in the W/B ratio, consistent with a Fibonacci sequence, marking metastable states between disorder and order. Identical W/B ratios can produce distinct topologies, indicating multistability under fixed composition. In Terraforms, the W/B ratio governs short- and long-range bonding, analogous to ecological interactions, while clay acts as a foundational agent, comparable to habitat-forming species. Our results show that complex, functionally optimized network structures can emerge solely from material interactions dictated by RMPs, without systemic control, biological processes, or environmental heterogeneity. These findings establish design rules for programmable, nature-inspired materials with transformative potential for eco-functional infrastructure, ecosystem restoration, and the creation of integrated natural–built environments with deliberately engineered functions — a process we term ecological terraforming.

Fish ecosystems perform ecological functions that are critically important for the sustainability of marine ecosystems, such as global food security and carbon stock. During the 21st century, significant global warming caused by climate change has created pressing challenges for fish ecosystems that threaten species existence and global ecosystem health. Here, we study a coastal fish community in Maizuru Bay, Japan, and investigate the relationships between fluctuations of ST, abundance-based species interactions and salient fish biodiversity. Observations show that a local 20% increase in temperature from 2002 to 2014 underpins a long-term reduction in fish diversity (∼25%) played out by some native and invasive species (e.g. Chinese wrasse) becoming exceedingly abundant; this causes a large decay in commercially valuable species (e.g. Japanese anchovy) coupled to an increase in ecological productivity. The fish community is analyzed considering five temperature ranges to understand its atemporal seasonal sensitivity to ST changes, and long-term trends. An optimal information flow model is used to reconstruct species interaction networks that emerge as topologically different for distinct temperature ranges and species dynamics. Networks for low temperatures are more scale-free compared to ones for intermediate (15-20°C) temperatures in which the fish ecosystem experiences a first-order phase transition in interactions from locally stable to metastable and globally unstable for high temperatures states as suggested by abundance-spectrum transitions. The dynamic dominant eigenvalue of species interactions shows increasing instability for competitive species (spiking in summer due to intermediate-season critical transitions) leading to enhanced community variability and critical slowing down despite higher time-point resilience. Native competitive species whose abundance is distributed more exponentially have the highest total directed interactions and are keystone species (e.g. Wrasse and Horse mackerel) for the most salient links with cooperative decaying species. Competitive species, with higher eco-climatic memory and synchronization, are the most affected by temperature and play an important role in maintaining fish ecosystem stability via multitrophic cascades (via cooperative-competitive species imbalance), and as bioindicators of change. More climate-fitted species follow temperature increase causing larger divergence divergence between competitive and cooperative species. Decreasing dominant eigenvalues and lower relative network optimality for warmer oceans indicate fishery more attracted toward persistent oscillatory states, yet unpredictable, with lower cooperation, diversity and fish stock despite the increase in community abundance due to non-commercial and venomous species. We emphasize how changes in species interaction organization, primarily affected by temperature fluctuations, are the backbone of biodiversity dynamics and yet for functional diversity in contrast to taxonomic richness. Abundance and richness manifest gradual shifts while interactions show sudden shift. The work provides data-driven tools for analyzing and monitoring fish ecosystems under the pressure of global warming or other stressors. Abundance and interaction patterns derived by network-based analyses proved useful to assess ecosystem susceptibility and effective change, and formulate predictive dynamic information for science-based fishery policy aimed to maintain marine ecosystems stable and sustainable.

Since 2013, the Centers for Disease Control and Prevention (CDC) has hosted an annual influenza season forecasting challenge. The 2015–2016 challenge consisted of weekly probabilistic forecasts of multiple targets, including fourteen models submitted by eleven teams. Forecast skill was evaluated using a modified logarithmic score. We averaged submitted forecasts into a mean ensemble model and compared them against predictions based on historical trends. Forecast skill was highest for seasonal peak intensity and short-term forecasts, while forecast skill for timing of season onset and peak week was generally low. Higher forecast skill was associated with team participation in previous influenza forecasting challenges and utilization of ensemble forecasting techniques. The mean ensemble consistently performed well and outperformed historical trend predictions. CDC and contributing teams will continue to advance influenza forecasting and work to improve the accuracy and reliability of forecasts to facilitate increased incorporation into public health response efforts.

The human microbiome is an extremely complex ecosystem considering the number of bacterial species, their interactions, and its variability over space and time. Here, we untangle the complexity of the human microbiome for the Irritable Bowel Syndrome (IBS) that is the most prevalent functional gastrointestinal disorder in human populations. Based on a novel information theoretic network inference model, we detected potential species interaction networks that are functionally and structurally different for healthy and unhealthy individuals. Healthy networks are characterized by a neutral symmetrical pattern of species interactions and scale-free topology versus random unhealthy networks. We detected an inverse scaling relationship between species total outgoing information flow, meaningful of node interactivity, and relative species abundance (RSA). The top ten interacting species are also the least relatively abundant for the healthy microbiome and the most detrimental. These findings support the idea about the diminishing role of network hubs and how these should be defined considering the total outgoing information flow rather than the node degree. Macroecologically, the healthy microbiome is characterized by the highest Pareto total species diversity growth rate, the lowest species turnover, and the smallest variability of RSA for all species. This result challenges current views that posit a universal association between healthy states and the highest absolute species diversity in ecosystems. Additionally, we show how the transitory microbiome is unstable and microbiome criticality is not necessarily at the phase transition between healthy and unhealthy states. We stress the importance of considering portfolios of interacting pairs versus single node dynamics when characterizing the microbiome and of ranking these pairs in terms of their interactions (i.e., species collective behavior) that shape transition from healthy to unhealthy states. The macroecological characterization of the microbiome is useful for public health and disease diagnosis and etiognosis, while species-specific analyses can detect beneficial species leading to personalized design of pre- and probiotic treatments and microbiome engineering.

The microbiome emits informative signals of biological organization and environmental pressure that aid ecosystem monitoring and prediction. Are the many signals reducible to a habitat-specific portfolio that characterizes ecosystem health? Does an optimally structured microbiome imply a resilient microbiome? To answer these questions, we applied our novel Eco-Evo Mandala to bacterioplankton data from four habitats within the Great Barrier Reef, to explore how patterns in community structure, function and genetics signal habitat-specific organization and departures from theoretical optimality. The Mandala revealed communities departing from optimality in habitat-specific ways, mostly along structural and functional traits related to bacterioplankton abundance and interaction distributions (reflected by ϵ and λ as power law and exponential distribution parameters), which are not linearly associated with each other. River and reef communities were similar in their relatively low abundance and interaction disorganization (low ϵ and λ) due to their protective structured habitats. On the contrary, lagoon and estuarine inshore reefs appeared the most disorganized due to the ocean temperature and biogeochemical stress. Phylogenetic distances (D) were minimally informative in characterizing bacterioplankton organization. However, dominant populations, such as Proteobacteria, Bacteroidetes, and Cyanobacteria, were largely responsible for community patterns, being generalists with a large functional gene repertoire (high D) that increases resilience. The relative balance of these populations was found to be habitat-specific and likely related to systemic environmental stress. The position on the Mandala along the three fundamental traits, as well as fluctuations in this ecological state, conveys information about the microbiome’s health (and likely ecosystem health considering bacteria-based multitrophic dependencies) as divergence from the expected relative optimality. The Eco-Evo Mandala emphasizes how habitat and the microbiome’s interaction network topology are first- and second-order factors for ecosystem health evaluation over taxonomic species richness. Unhealthy microbiome communities and unbalanced microbes are identified not by macroecological indicators but by mapping their impact on the collective proportion and distribution of interactions, which regulates the microbiome’s ecosystem function.