The study area spans 34,484 square kilometers across the cities of Linfen and Yuncheng, a region long dominated by coal mining, chemical production, and heavy industry. During the winter heating season, cold, dry Siberian air masses settle over the basin, wind speeds drop, and the surrounding mountains trap pollutants close to their sources. Coal-fired heating and open straw burning add to the burden, making winter the season with the most severe and complex pollution episodes. It is precisely this combination of topography, meteorology, and intensive human activity that made southern Shanxi an ideal natural laboratory for asking whether the arrangement of urban landscapes can measurably shape air quality and carbon output.

To answer that question, the researchers assembled an unusually detailed dataset. Carbon emissions came from the Open Data Inventory for Anthropogenic Carbon Dioxide, a 1-kilometer gridded product built by disaggregating nighttime lights, power plant emissions, and sectoral energy use. Particulate matter concentrations for PM2.5 and PM10 were drawn from 66 ground monitoring stations run by Shanxi’s environmental authorities between November 2022 and March 2023, then interpolated into continuous pollution surfaces using Bayesian kriging—a method that quantifies its own uncertainty and copes well with unevenly spaced stations. Cross-validation showed the interpolation performed reliably, with coefficients of determination of 0.714 for PM2.5 and 0.705 for PM10.

The landscape itself was characterized at remarkable resolution. Sentinel-2A imagery at 10 meters, processed with atmospheric correction and object-based classification, divided the land into green land, wetland, cropland, water body, and impervious surface, achieving an overall classification accuracy of 83.2 percent. For the vertical dimension of the city, the team turned to Gaofen-2 imagery at 0.8 meters, exploiting the geometric relationship between building shadows and true building heights—validated with field measurements—to calculate metrics such as floor area ratio, building coverage ratio, average height, and a newly constructed building rugosity index. Green space was described with eight patch-level indices capturing coverage, patch density, shape, contiguity, and connectivity, supplemented by satellite-derived net primary productivity as a measure of vegetation’s carbon-sequestration capacity.

The resulting maps revealed a striking and somewhat counterintuitive geography. Carbon emissions clustered intensely in the industrial heartlands—Yaodu District, Yanhu District, and Hejin City—where maximum emission intensity reached 168,486 tonnes per square kilometer. Particulate matter, by contrast, followed a southwest-to-northeast gradient, highest around Yanhu District and lowest in the greener, better-ventilated northeastern counties such as Gu, Anze, and Fushan. Crucially, the two pollution types overlapped locally but diverged overall: Yanhu District was a hotspot for both, yet Yaodu District, a core carbon-emission zone, did not show correspondingly peak particulate levels. County-level Pearson analysis confirmed only weak positive correlations between emissions and PM2.5 (r = 0.422) or PM10 (r = 0.362), a finding the authors attribute in part to Yaodu’s service-heavy economy, higher penetration of natural gas and electricity, centralized heating, and strict emission controls in the urban core.

When the team correlated seventeen landscape indices against emissions and pollution, a clear pattern emerged. Architectural intensity was consistently bad news: floor area ratio showed a strong positive correlation with carbon emissions (r = 0.716 in winter), and the building rugosity index correlated strongly with both PM2.5 (r = 0.845) and PM10 (r = 0.759). The mechanism is twofold. Denser, taller construction directly increases energy demand and emissions, while rougher urban canopies modify airflow, restrict ventilation pathways, create recirculation zones, and prolong the time that particles linger in the air people breathe. Green landscapes told the opposite story. The percentage of vegetated land correlated negatively with PM2.5 at a remarkable r = −0.930, and patch cohesion, largest patch index, and net primary productivity all tracked lower pollution and lower emissions. Vegetation intercepts and deposits particles, stabilizes carbon uptake, and forms ecological corridors that aid dispersion. Notably, carbon emissions proved more sensitive to vegetation productivity, while particulate concentrations responded more strongly to the spatial morphology of green patches and the structural form of buildings.

Because averages can hide local variation, the researchers applied geographically weighted regression, which fits a separate set of coefficients at every location rather than one global average. The models explained the spatial variation well—adjusted R-squared values of 0.769 for carbon emissions, 0.680 for PM2.5, and 0.789 for PM10—and the local coefficients revealed pronounced spatial heterogeneity, with strong effects concentrated in the northern, southwestern, and northeastern parts of the region and weaker effects in the center and southeast. Spatial autocorrelation analysis using Moran’s I with 999 Monte Carlo simulations confirmed significant clustering for both variables, with particulate pollution showing stronger spatial clustering than carbon emissions, consistent with the regional transport of aerosols across county boundaries.

One of the study’s most elegant analyses focused on Yuncheng Salt Lake, a unique inland saline ecosystem ringed by thirty concentric 1-kilometer buffer zones extending 30 kilometers from the lake’s edge. Carbon emissions and particulate concentrations both declined steadily outward from the lake, tracking the fading intensity of building development. In the near zone, where floor area ratios and building coverage remained high, the link between construction and pollution was strongest; in the far zone, it essentially vanished, apart from a localized emissions peak at 29 kilometers traced to the industrial built-up area of Yongji City rather than a data artifact. Meanwhile, the negative regulatory effect of green space grew stronger with distance from the lake: fragmented green patches near the shore offered limited protection, while the large, connected, highly productive green landscapes of the outer suburbs delivered the strongest pollution-blocking and carbon-sinking benefits.

The practical implications are direct. Rather than treating smog and carbon as a single problem, the authors argue for differentiated landscape planning tuned to each county’s dominant burden. Areas prioritizing particulate control, such as Yonghe and Xi counties, should build protective forest belts and ventilation corridors along mountains and river valleys while strictly regulating building forms and floor area ratios. Areas prioritizing carbon reduction, such as Wenxi and Daning counties, should expand ecological green space and vegetation productivity while curbing construction intensity. Regions needing both, including Ruicheng, Yongji, Linyi, and Hejin, should weave core carbon-sink green spaces, dust-retention buffers, and ventilation corridors into a single ecological network alongside compact, low-carbon building layouts. Around the salt lake itself, the team proposes a gradient framework: strict lakeside development restrictions and greenway connections in the near zone, large ecological patches and industrial-residential separation in the middle zone, and a reinforced mountain-and-farmland ecological barrier in the far zone. The authors caution that their analysis covers only one winter heating season, includes only building and green space landscape types, and cannot fully disentangle meteorological and socioeconomic drivers—but the core message stands. In China’s coal country, and in resource-based regions across the Yellow River Basin, the cheapest combined attack on smog and carbon may already be growing in the ground: well-connected, productive green landscapes, deployed where the data say they will work hardest.

Subject of Research: Driving mechanisms of atmospheric particulate matter and carbon emissions under the influence of landscape patterns in southern Shanxi, China

Article Title: Analysis of the driving mechanisms of atmospheric particulate matter and carbon emissions under the influence of landscape patterns in southern Shanxi, China

Article References: Ren, W., Zhao, J., & Hu, F. (2026). Analysis of the driving mechanisms of atmospheric particulate matter and carbon emissions under the influence of landscape patterns in southern Shanxi, China. Frontiers in Environmental Science, 14, Article 1944157. https://doi.org/10.3389/fenvs.2026.1944157

Image Credits: AI Generated

DOI: 10.3389/fenvs.2026.1944157

Keywords: air pollution, carbon emissions, PM2.5, landscape patterns, urban green space, geographically weighted regression, remote sensing, Shanxi, Yellow River Basin, Yuncheng Salt Lake, winter heating season, spatial heterogeneity

The research team, led by scientists from the Ningxia Academy of Agricultural and Forestry Sciences and collaborating institutions in China and Egypt, collected 487 topsoil samples from a 4,670-square-kilometer irrigation district spanning the cities of Yinchuan, Shizuishan, Wuzhong, and Zhongwei. Using a stratified random design that balanced farmland sites with adjacent forest and grassland plots, they measured concentrations of seven potentially hazardous elements: arsenic, cadmium, chromium, copper, nickel, lead, and zinc. Every sample was digested under strict EPA protocols and analyzed with graphite furnace atomic absorption spectrophotometry and inductively coupled plasma atomic emission spectrometry, with certified reference materials confirming recoveries between 96 and 105 percent.

The results paint a picture of contamination that is anything but uniform. Cadmium, the most spatially erratic of the seven metals, showed a coefficient of variation of roughly 56 percent, while arsenic and lead varied by more than a third across the landscape. Mean cadmium concentrations reached 0.32 milligrams per kilogram — more than three times the regional background of 0.097 — and just above the national screening threshold of 0.30 milligrams per kilogram set by China’s GB 15618-2018 standard. Lead, zinc, chromium, copper, and nickel all exceeded their respective background values as well, though most remained comfortably below regulatory limits, suggesting the contamination is real but not yet at crisis levels.

To untangle what was driving this patchwork of pollution, the researchers deployed two complementary computational tools. The first, an Optimal Parameters-based Geographical Detector (OPGD), quantifies how much of the spatial variation in metal concentrations can be explained by each of eleven environmental variables — from soil pH and texture to elevation, vegetation cover, and population density. The second, Positive Matrix Factorization (PMF), is a receptor model that mathematically deconstructs the measured concentrations into distinct source signatures, assigning each metal a proportional origin.

The OPGD results delivered the study’s most striking revelation. Total nitrogen, a proxy for long-term fertilizer and manure inputs, explained 79.6 percent of the spatial variability in the metal dataset — the highest Q-statistic of any factor tested. For chromium alone, nitrogen accounted for a Q-value of 0.796, while copper and zinc followed closely at 0.612 and 0.627 respectively. Soil texture, particularly the relative proportions of sand, silt, and clay, emerged as the second most important control, reflecting how fine particles adsorb and retain metal ions. Organic matter, elevation, and population density played secondary roles, and vegetation cover contributed surprisingly little.

But the story grew more intricate when the researchers examined factor interactions. Clay content alone explained only 13 percent of arsenic’s spatial pattern, yet when combined with soil organic matter, the joint explanatory power leapt to 81 percent — a nonlinear enhancement that suggests arsenic is held in place by a coupled mechanism of organic complexation and adsorption onto fine particles. Similarly, the interaction between elevation and organic matter explained 89 percent of chromium’s variability, pointing to terrain-driven redistribution working in tandem with organic retention. These synergies reveal that no single environmental lever controls metal fate; instead, soil chemistry operates as an interlocking system.

The PMF analysis then answered the question of where the metals actually come from. Five distinct source profiles emerged, together explaining more than 90 percent of the total variance. Cadmium was overwhelmingly tied to agricultural inputs — phosphate fertilizers and pesticides — accounting for 77 percent of its factor profile. Zinc and copper clustered around industrial emissions, with loadings of 56 and 38 percent respectively, concentrated near manufacturing hubs in Shizuishan and Wuzhong. Chromium and nickel traced back to geogenic parent material, reflecting the natural lithology of the Yellow River alluvial plain. Lead, with a striking 89 percent loading on a single factor, pointed squarely at traffic residues — a legacy of historical leaded gasoline alongside modern brake and tire wear. Arsenic, at 79 percent of a mixed anthropogenic factor, suggested combined inputs from irrigation water and diffuse industrial discharge.

When the team overlaid these source maps with the OPGD-derived risk layers, a composite hazard index revealed four contamination zones across the irrigation district. Low-risk areas, dominated by forest and grassland, sat far from human influence. Moderate-risk transitional zones gave way to high-risk farmland where cadmium and lead levels were elevated, and very-high-risk clusters appeared near industrial and urban centers where fertilizer use, traffic emissions, and industrial discharge converged. Cross-correlation between the independently derived PMF factor scores and OPGD risk layers yielded coefficients between 0.62 and 0.81, providing strong statistical validation that the two models were capturing the same underlying reality.

The spatial pattern of risk was anything but random. Hotspots concentrated in the central and lower western reaches of the irrigation district, with a prominent secondary cluster in the southeast. Arsenic and lead risks were largely confined to the central-western section, while nickel — the most broadly distributed contaminant — formed three distinct clusters spanning the central-western, central-eastern, and southeastern regions. Cadmium, chromium, and zinc hotspots clustered in the upstream portion of the irrigation system, where water enters the network first and presumably carries the highest dissolved loads.

What makes this study particularly consequential is what it implies about the future of food security in irrigated agroecosystems worldwide. The Yellow River Basin irrigates roughly 15 percent of China’s farmland and supports nearly 12 percent of its population, and the researchers note that comparable systems — the Indo-Gangetic Plain, the Nile Delta, the Murray–Darling Basin — face similar pressures from fertilizer dependence and industrial emissions. The finding that total nitrogen, a manageable agronomic variable, exerts the strongest statistical control over metal distribution suggests that precision fertilization, reduced over-application, and the adoption of low-cadmium phosphate sources could meaningfully slow contamination before it reaches irreversible thresholds. The authors caution that their hazard maps reflect contamination intensity rather than bioavailability or ecological toxicity, and they recommend future work incorporating DTPA-extractable metal fractions, soil enzyme assays, and human health risk models to translate spatial patterns into actionable exposure estimates. For now, the message is clear: the nitrogen farmers apply to boost yields is also writing a chemical record in the soil, and reading that record may be the key to keeping the Yellow River’s breadbasket productive for generations to come.

Subject of Research: Spatial drivers and source apportionment of potentially hazardous elements in irrigated Yellow River Basin soils

Article Title: Spatial heterogeneity, source contributions, and hazard mapping of potentially hazardous elements in irrigated soils of the Yellow River Basin, Northwest China

Article References: Xu, H., Liu, X., Chen, S., Luo, L., Wenfan, W., Ma, L., AbdelRahman, M. A. E., & Sabra, D. M. (2026). Spatial heterogeneity, source contributions, and hazard mapping of potentially hazardous elements in irrigated soils of the Yellow River Basin, Northwest China. Frontiers in Soil Science, 6, Article 1966899. https://doi.org/10.3389/fsoil.2026.1966899

Image Credits: AI Generated

DOI: 10.3389/fsoil.2026.1966899

Keywords: heavy metals, soil contamination, Yellow River Basin, cadmium, fertilizer, PMF, OPGD, spatial heterogeneity, irrigated agriculture, source apportionment, hazard mapping, Ningxia

The research, led by Tianyu Wu and Zaida Luthey-Schulten together with colleagues including Marie-Christin Spindler, Abner T. Apsley, Emmy Earnest, Zane R. Thornburg, and Julia Mahamid, focuses on the galactose utilization network of budding yeast, Saccharomyces cerevisiae. This circuit, often called the galactose switch, has been a workhorse of quantitative biology for half a century. When glucose is absent and galactose is present, a cascade of transcriptional activation flips a set of genes on, allowing the cell to import and metabolize the alternative sugar. Because its components are so well characterized, the galactose switch is a natural test bed for asking how much spatial organization actually contributes to regulatory dynamics.

At the heart of the new work is a modeling framework that departs sharply from the well-mixed tradition. The team built spatially resolved models that integrate experimentally derived intracellular architectures, including the organization of the yeast chromosomes, the endoplasmic reticulum, and distinct spatial populations of ribosomes, the molecular machines that translate messenger RNA into protein. Electron microscopy and biochemical constraints guided the construction of these geometries, so that the simulated cell reflects what is actually observed in real yeast rather than an idealized compartment.

Computationally, the authors adopted a hybrid strategy. Gene expression is handled by a reaction–diffusion master equation, a stochastic formalism in which the cellular volume is subdivided into small lattice sites. Molecules diffuse between neighboring sites, and chemical reactions fire only when the participating species occupy the same site or adjacent sites. This enforces what modelers call locality: a transcription factor cannot activate a gene it has not physically reached, and a messenger RNA cannot be translated by a ribosome that is nowhere nearby. In contrast, the faster metabolic and transport processes are captured by ordinary differential equations, which keeps the hybrid scheme computationally tractable without sacrificing spatial realism where it counts most.

The results offer a nuanced picture of which features of cellular geography matter and which do not. When the team varied the geometry of the chromosomes within the nucleus, the output of Gal2p, the galactose permease responsible for importing galactose across the plasma membrane, was barely affected in the present model. Chromosome positioning, in other words, appears not to be a dominant determinant of how much permease the network produces, at least under the conditions examined.

The picture changes dramatically when translation is considered. Gal2p is a membrane protein, and membrane proteins in yeast are typically inserted into the endoplasmic reticulum before being trafficked to the plasma membrane. When the researchers restricted the translation of GAL2 messenger RNA to the population of ribosomes that are physically bound to the ER surface, they observed the largest decrease in Gal2p abundance across all the spatial scenarios they tested. The explanation lies in the journey itself: protein produced at the ER must be delivered through the secretory pathway to the cell surface, and this spatial separation between the site of synthesis and the site of function imposes a real cost on the amount of functional permease that arrives at the membrane.

ER-associated translation more broadly reduced Gal2p delivery to the plasma membrane in the model, even when the restriction was less extreme. This finding carries a conceptual punch, because it means that a purely well-mixed model would systematically overestimate the amount of transporter available to the cell. The well-mixed assumption treats synthesis and delivery as if they happened in the same chemical pot, erasing the transit steps that in a real, compartmentalized cell determine how much protein actually reaches its destination.

What makes the study more than a technical exercise is its general implication: spatial realism can qualitatively change regulatory predictions, not merely tweak them slightly. A model that ignores intracellular organization does not simply miss a small correction term; it can miss the direction of an effect entirely. For the galactose switch, the ordering of which spatial features matter is itself a prediction, ranking ER-bound translation above chromosome geometry as a determinant of permease output. Experiments that probe transporter abundance under conditions that alter ER association or trafficking could now test that hierarchy directly.

The work also speaks to a larger ambition in computational biology: the construction of whole-cell models that simulate every molecule of a living organism. Landmark efforts such as the well-known whole-cell model of a minimal bacterium demonstrated that comprehensive simulation is feasible, but they, like most gene-regulation models, leaned on well-mixed assumptions for many cellular processes. The yeast galactose study provides a template for how experimentally grounded three-dimensional structures, from cryo-electron microscopy-derived cell architecture to mapped ribosome populations, can be folded into such models, and it warns that the extra realism may not be optional if the goal is quantitative accuracy.

There is also a practical message for the growing community of researchers who model eukaryotic cells. The hybrid stochastic–deterministic scheme used here shows one way to pay the computational price of spatial resolution only where it matters, reserving fast deterministic treatment for metabolism and transport while enforcing locality for gene expression. As imaging technologies such as cryo-electron tomography, to which the Mahamid group has contributed extensively, continue to deliver high-resolution maps of organelles, genomes, and macromolecular complexes inside intact cells, the raw material for spatially resolved modeling is becoming abundant. The challenge, this study suggests, is no longer whether cellular geometry can be modeled, but which architectural features biologists can afford to ignore. For the galactose switch, the answer is now on the record: the ER-ribosome axis cannot be ignored, and the era of pretending the yeast cell is a well-stirred bag may be drawing to a close.

Subject of Research: Spatially resolved reaction–diffusion modeling of the yeast galactose network

Article Title: Spatially resolved reaction–diffusion modeling reveals effects of intracellular spatial heterogeneity on yeast galactose network dynamics

Article References: Wu, T., Spindler, M.-C., Apsley, A. T., Earnest, E., Thornburg, Z. R., Mahamid, J., & Luthey-Schulten, Z. (2026). Spatially resolved reaction–diffusion modeling reveals effects of intracellular spatial heterogeneity on yeast galactose network dynamics. PLOS Computational Biology, 22(10), e1014811. https://doi.org/10.1371/journal.pcbi.1014811

Image Credits: AI Generated

DOI: 10.1371/journal.pcbi.1014811

Keywords: yeast, galactose switch, reaction–diffusion, systems biology, gene regulation, endoplasmic reticulum, ribosomes, chromosome organization, whole-cell modeling, Saccharomyces cerevisiae, computational biology, spatial heterogeneity

The biological motivation is broad. Bacteria colonizing a surface and metastatic cancer cells spreading through an organ both face the same fundamental situation: they grow, they migrate, and they occupy patchy environments in which conditions differ sharply from one microhabitat to the next. Migration creates the structure itself, linking colonies into networks that can be critical for colonization or metastasis. When a drug is applied, it too is distributed unevenly, penetrating some sites well and others poorly. The interplay between growth, migration, and this spatial heterogeneity makes population-level responses remarkably difficult to predict. Two environments with the same spatially averaged growth rate can produce opposite fates—one collapsing under treatment, the other persisting—because connectivity and the placement of low-drug refuges matter as much as the mean.

Hu and Wood’s central contribution is a minimal growth-migration model that captures population dynamics on discrete microhabitat structures under spatial drug heterogeneity. Each microhabitat is treated as a node with its own local growth rate, determined in part by the local drug concentration, and nodes are connected by migration pathways whose strengths encode how easily cells move between sites. The model is deliberately stripped down, containing only the essential ingredients of growth and movement, which is precisely what allows the authors to extract exact, general results rather than case-by-case simulations. The question they pose is the clinically urgent one: given a particular spatial arrangement of growth rates and drug exposures, will the total population grow or decline?

The mathematical machinery that answers this question is the most technically striking part of the paper. The authors apply a kernel transformation that maps the original, possibly complicated, microhabitat structure onto an effective fully connected graph, in which every site is linked to every other. On this transformed graph they derive a new exact criterion for population response, expressed in terms of a regularized Laplacian kernel reweighted by the local growth rates. In plain terms, the criterion condenses the entire spatial problem—every node, every migration link, every local drug concentration—into a single mathematical object whose properties determine the population’s fate. If the reweighted kernel satisfies a specific condition, the population is guaranteed to decline; if it fails, growth is possible.

This criterion is not merely an abstract formula. The authors show that it connects directly to a concept from network science called forest closeness centrality, a measure of how centrally located nodes are within the connectivity structure of a graph. Through this connection they obtain analytical bounds and sufficient conditions for population growth or decline, meaning they can certify outcomes without simulating the dynamics at all. This is a significant practical advantage: instead of running many stochastic simulations of cell growth and migration under each candidate treatment plan, one can evaluate a kernel-based quantity and know immediately whether the strategy is guaranteed to push the population toward extinction.

Among the qualitative insights the framework delivers, one stands out for its counterintuitive flavor: higher structural connectivity, such as increased migration between microhabitats, generally promotes population decline under treatment. Intuition might suggest the opposite—that well-connected colonies should be more robust, since cells from drug-protected refuges could continually reseed vulnerable sites. The mathematics shows, however, that migration also drags cells out of safe havens and into drug-exposed territory, mixing the population so that no subpopulation can remain sheltered. In the effective fully connected picture, strong connectivity spreads the lethal influence of high-drug, low-growth regions throughout the network. For treatment design, this suggests a provocative principle: rather than trying to isolate and starve each lesion or colony separately, enhancing the coupling between microhabitats—so that the drug’s effect propagates—can be a route to clearance.

The framework also addresses the inverse problem: not just predicting the response to a given drug distribution, but designing the optimal spatial assignment of drug to guarantee population decline. The authors show that this optimization reduces to selecting interconnected subcores in the effective complete graph—compact, mutually reinforcing sets of microhabitats whose combined treatment overwhelms the population. This geometric characterization turns a potentially intractable continuous optimization into a combinatorial selection problem with clear structure. Moreover, the authors handle the realistic complications that clinicians and microbiologists actually face: when only some microhabitats can be targeted, or when the drug distribution is unknown, they identify strategies that nonetheless ensure decline, providing robust guarantees under partial control and uncertainty.

To demonstrate the clinical relevance of the theory, the authors apply it to a parameterised model of colorectal cancer metastasis under systemic chemotherapy. Metastatic colorectal cancer commonly seeds multiple organs, each with distinct drug penetration characteristics, and systemic chemotherapy exposes all sites to the same regimen. When the framework is applied to this multi-organ setting, it identifies the lung as a critical growth reservoir: a site where poor drug penetration allows the metastatic population to sustain itself even when the spatially averaged drug response would predict decline. This is exactly the failure mode the theory was built to expose—an average-based calculation says the population should shrink, but the spatial structure says otherwise, because one poorly penetrated sanctuary continuously replenishes the rest.

The implication of that finding is concrete and testable: lung-targeted drug delivery emerges as a route to robust clearance in this model. By improving penetration at the critical reservoir, the spatially averaged prediction and the true population-level response can be brought back into agreement, and the guarantee of decline that the kernel criterion provides can actually be met. The authors are careful to frame this as a theoretical illustration on a parameterised model rather than a treatment recommendation, but the logic is general. Any metastatic or biofilm-associated disease in which drug penetration varies across sites could be analyzed the same way, and the framework would flag which specific sites must be better targeted for the treatment to succeed.

Taken together, the study offers a new theoretical perspective on drug response in spatially structured populations and a practical toolkit for optimizing spatially explicit dosing strategies in heterogeneous environments. Its exact criterion, its connection to network centrality, and its guarantees under partial control elevate the work beyond simulation-based studies of spatial heterogeneity, giving researchers a way to certify outcomes analytically. As single-cell and spatial profiling technologies continue to reveal how thoroughly heterogeneous tumors and microbial communities really are, tools of this kind—ones that treat spatial structure as the central object rather than a nuisance—may become essential to the design of therapies that work not just on average, but everywhere the disease actually lives.

Subject of Research: Population-level drug response in spatially heterogeneous microhabitat structures

Article Title: Deciphering and steering population-level response under spatial drug heterogeneity on microhabitat structures

Article References: Hu, Z., & Wood, K. (2026). Deciphering and steering population-level response under spatial drug heterogeneity on microhabitat structures. PLOS Complex Systems, 3(6), e0000113. https://doi.org/10.1371/journal.pcsy.0000113

Image Credits: AI Generated

DOI: 10.1371/journal.pcsy.0000113

Keywords: spatial heterogeneity, growth-migration model, microhabitat, drug response, Laplacian kernel, network centrality, metastasis, colorectal cancer, chemotherapy, drug penetration, population dynamics, dosing strategy

The study, published in PLOS Water, set out with a deceptively simple question. If you adjust water consumption in every locality for the same city-wide drivers—climate swings, socioeconomic trends, the shared shocks of a period that included a global pandemic and severe drought stress—do the localities then behave as statistical clones of one another, differing only by a fixed offset? The additive model the authors built first assumes exactly that. It allows each locality to sit higher or lower than the city average, capturing a static spatial ranking, but it forces every locality to follow the same temporal trajectory once common covariates are accounted for. In a city as geographically and socially layered as Bogotá, where highland neighborhoods, dense commercial cores, and sprawling southern districts experience water very differently, that assumption was always going to be under strain.

To test it rigorously, the researchers turned to the machinery of modern Bayesian inference. Rather than producing a single best-fit line, Bayesian models express uncertainty as full posterior distributions over every parameter, which makes them well suited to the noisy, short, and irregularly structured data that characterize municipal water records. The authors fitted models to monthly consumption at the locality level across six years, a period long enough to span dramatic climatic and social disruption but short enough that data-hungry methods must be handled with care. The additive specification adjusted for common temporal dynamics and city-wide climatic and socioeconomic time series, then asked how much spatial structure remained in the residuals.

The answer came from a formal sensitivity analysis, and it was unambiguous. When the researchers added a structured space–time interaction—a term that lets each locality’s deviation from the common trajectory evolve over time, borrowing strength from neighboring localities through spatially correlated priors—the improvement in model fit was enormous. The Widely Applicable Information Criterion, a standard Bayesian measure of predictive fit that penalizes complexity, dropped from 2502.7 in the additive model to 604.0 in the interaction model. The Deviance Information Criterion fell in parallel, from 2500.5 to 333.9. In the vocabulary of statistical modeling, gaps of this magnitude are not marginal refinements; they signal that the simpler model is missing something fundamental about how the data are generated. Locality deviations, the authors conclude, are not adequately represented as constant throughout the study period.

What does that mean in practical terms? Imagine two neighborhoods that consumed similar volumes of water in 2019. Under the additive model, they would be expected to track each other through the pandemic lockdowns, the rationing anxieties, and the shifting rainfall patterns of the following years, differing only by a fixed amount. The interaction model says otherwise: the same shared events landed differently in different places, and those differential responses are themselves spatially organized. Neighborhoods that are geographically close tend to deviate from the city-wide pattern in similar ways, a signature consistent with localized infrastructure, microclimates, settlement histories, and social composition shaping how water demand responds to stress. The spatial structure is not a static backdrop but an active, evolving dimension of consumption.

The study also confronted a subtler statistical hazard: heavy tails. Water consumption data, like many environmental and socioeconomic series, occasionally produce extreme observations—a billing anomaly, a meter failure, a localized crisis—that a normal distribution treats as vanishingly rare. The researchers therefore compared their Gaussian models against versions using a Student-t likelihood, which has heavier tails and down-weights extreme observations. The posterior mean degrees of freedom told the story: around 3.02 in the additive model and 2.60 in the interaction model. Values this low, well below the roughly 30 of a Gaussian distribution, indicate genuinely heavy residual tails—outliers are a real feature of these data, not a modeling artifact. The Student-t likelihood was not a paranoid precaution; it was describing something true about Bogotá’s water records.

Yet here the study delivers one of its most interesting and technically consequential findings. Despite the evidence for heavy tails, the Gaussian interaction model retained substantially better information criteria than its Student-t counterpart. That apparent paradox—acknowledging fat-tailed errors yet preferring the thinner-tailed model—reflects the way information criteria balance fit against effective complexity. The Student-t likelihood, with its additional flexibility, appears to absorb structure that the interaction term would otherwise capture, blurring the space–time signal the researchers most wanted to isolate. When the authors compared the principal spatial main effects estimated under the two likelihoods, they found a correlation of 0.975. In other words, the overall spatial ordering of localities—which places consume more, which consume less—is remarkably robust to the choice of error distribution, even as the finer temporal dynamics are sensitive to it.

That robustness matters for anyone hoping to act on these results. A water utility deciding where to target conservation campaigns, where to prioritize leak detection, or how to design tariff structures needs to trust the spatial ranking of demand. The near-perfect correlation across likelihoods suggests that ranking is solid ground. But a utility trying to forecast how a specific district will respond to the next drought, or to evaluate whether an intervention changed behavior in one locality but not its neighbors, needs the space–time interaction. The additive model, however convenient, would systematically misattribute those differential dynamics, potentially crediting or blaming city-wide policies for changes that were in fact deeply local.

The methodological lesson radiates well beyond Bogotá. Cities across Latin America and the wider global South face similar conditions: administratively defined zones of uneven size and character, consumption records of moderate length, and climatic pressures—like the multi-year rainfall deficits that have strained Andean reservoirs in recent years—that do not fall uniformly across the urban landscape. The Bayesian framework used here offers a template because it handles short series gracefully, quantifies uncertainty honestly, and makes its assumptions explicit and testable. The sensitivity analysis is the crucial move. Rather than asserting a single model, the authors demonstrated that their central conclusion—spatially heterogeneous temporal dynamics—survives a deliberate stress test of its core assumptions, from the error distribution to the structure of the interaction term.

There is also a quieter, more civic implication. When consumption is modeled as a single city-wide curve, the localities that deviate most—those that conserve aggressively under stress, or those that cannot—become statistical noise. The space–time interaction model makes those deviations visible and locatable, turning them into evidence. In a period when Bogotá and many other cities are confronting the reality that water infrastructure, climate, and inequality intersect block by block, that visibility is not an academic luxury. It is the difference between managing a city as an abstraction and managing it as the mosaic of places it actually is. The maps hidden inside Bogotá’s water bills, this study shows, were always there; it took a model willing to let space and time talk to each other to reveal them.

Subject of Research: Bayesian spatio-temporal modeling of locality-level urban water consumption in Bogotá, Colombia

Article Title: Spatial structure in urban water consumption: A Bayesian space-time analysis of Bogotá (2019–2024)

Article References: Cruz-Reyes, D. L., & Ramírez Orozco, D. L. (2026). Spatial structure in urban water consumption: A Bayesian space-time analysis of Bogotá (2019–2024). PLOS Water, 5(9), e0000512. https://doi.org/10.1371/journal.pwat.0000512

Image Credits: AI Generated

DOI: 10.1371/journal.pwat.0000512

Keywords: urban water consumption, Bogotá, Bayesian modeling, spatio-temporal analysis, PLOS Water, space-time interaction, Student-t likelihood, WAIC, water management, spatial heterogeneity, climate, sensitivity analysis

The crux of the problem is spatial heterogeneity. Conventional evapotranspiration models were developed for agricultural fields and large-scale hydrological budgets, where vegetation is comparatively uniform. A bioretention basin is nothing of the sort. Within a few meters, plant height can jump from bare ground to four-meter shrubs, buildings and highway walls cast shifting shadows, and soil moisture swings dramatically between the flooded basin floor and its drier flanks. Because evapotranspiration, or ET, is governed by leaf area, solar radiation, vapor pressure deficit, and soil water availability, this patchiness propagates directly into uncertainty. Reported estimates of how much stormwater ET removes from green infrastructure range from nearly zero to almost 90 percent, a spread so wide it borders on useless for design purposes.

To pin down what is actually happening, the team spent a growing season measuring ET directly at eleven fixed plots within the Philadelphia basin using the closed dynamic chamber method. Transparent chambers were sealed over collars set into the soil, and a laser-based gas analyzer recorded water vapor accumulation once per second, allowing the researchers to calculate fluxes from the rate of buildup. Measurements were taken five to six times per day across seven dates from June through October, deliberately spanning dry spells and post-storm conditions. After discarding about 17 percent of records for quality reasons, 302 usable flux measurements remained, each tied to a specific combination of plant stature, topographic position, and microclimate.

The spatial context came from an unusual pairing of technologies. Terrestrial LiDAR scanning produced a bare-earth digital elevation model of the basin, while drone-based photogrammetry yielded a digital surface model of everything above it; subtracting the two mapped vegetation height across the site. To capture shading from the highway wall, nearby buildings, and a billboard, the researchers used a regional aerial LiDAR dataset and the Area Solar Radiation calculator in ArcGIS to compute shade coefficients for every pixel on every day of the study. These coefficients, typically between 0.75 and 0.95, scaled the measured solar radiation to reflect how much light actually reached each patch of ground. The basin was then segmented into 967 roughly plant-sized polygons using an image-clustering algorithm, and all environmental variables were aggregated to that scale.

With these ingredients, the team built a statistical model of ET that explicitly incorporated plant height, topographic position, shade, vapor pressure deficit, solar radiation, and soil moisture, complete with interaction terms and random effects for plot and date. The model explained most of the observed variation, with fixed effects alone accounting for 65 percent of variance. Its coefficients told a coherent physiological story: ET rose strongly with solar radiation and vapor pressure deficit, was lower on the drier upper flanks of the basin, and increased with plant height, reflecting greater leaf area per unit ground. Interestingly, the effect of radiation weakened in wetter soils, likely because near-saturated conditions restrict both evaporation and plant transpiration, while the effect of moisture on ET differed between the basin floor and its slopes.

Scaling the model across the whole basin revealed daily ET ranging from zero to about 6 millimeters per day, peaking in early and midsummer and declining to roughly 1 millimeter per day by mid-autumn. The clearest spatial pattern was elevated ET beneath the tallest vegetation, particularly a cluster of shrubs hugging the highway wall. Just as telling was what the sensitivity analysis showed. Setting all plant heights to zero reduced season-wide ET by 41 percent, while assigning every segment a uniform 100-centimeter height inflated it by 60 percent; even using the basin’s mean height of 36 centimeters everywhere depressed the total by 13 percent. Shade mattered too: removing all structural shading boosted seasonal ET by 11 percent, while assuming only a quarter of incoming radiation reached the basin cut it nearly in half. Topographic position, by contrast, had little effect on the basin-wide total.

How did the conventional models fare? The team ran six widely used equations, including the two Penman-Monteith variants endorsed by the American Society of Civil Engineers and the FAO, Priestley-Taylor, Hargreaves-Samani, Matt-Shuttleworth, and Granger-Gray. Most overpredicted ET on average. The Penman-Monteith, Priestley-Taylor, and Hargreaves-Samani models consistently overestimated relative to the empirical estimates, while Matt-Shuttleworth and Granger-Gray, which were less responsive to atmospheric conditions, overpredicted at the low end of the range and underpredicted at the high end. Digging into the discrepancies, the researchers found that differences between Penman-Monteith and empirical estimates grew systematically with wind speed, explaining up to 39 percent of the bias for the ASCE variant, whereas accounting for soil-water limitation made little difference; water stress was infrequent during the study period.

The practical payoff came from testing three correction strategies. Applying shade coefficients to the solar radiation inputs improved agreement for the Penman-Monteith and Priestley-Taylor models. Multiplicative landscape coefficients, analogous to crop coefficients in agriculture and calculated here at roughly 0.5 to 1.0, also helped those same models. But the strongest results came from additive corrections: simply subtracting a regression-derived offset from Penman-Monteith and Priestley-Taylor estimates raised concordance correlation coefficients to about 0.88 or 0.89, a dramatic improvement. The caveat is that additive correction can produce negative daily ET values when uncorrected estimates are small, so it must be used with care. None of the corrections improved the already-decent Matt-Shuttleworth and Granger-Gray estimates, and the Hargreaves-Samani model remained the poorest performer regardless of adjustment, a notable finding given that some reviews have recommended it for its minimal data requirements.

The implications ripple outward from one Philadelphia basin. For designers of similar systems, the message is that Penman-Monteith and Priestley-Taylor, once corrected, can credibly represent ET in water balances used for regulatory compliance and volume-reduction credit, rather than being ignored or trusted blindly. The study also suggests a design lever: because plant height and shade dominate basin-scale ET, landscape architects can deliberately position larger plants in unshaded locations to boost atmospheric water removal, provided the coarse, fast-draining soil media can supply enough moisture between storms. Deep-rooted shrubs and trees may be particularly valuable here, drawing water from a broader depth range during dry spells.

The authors are candid about limitations. Their empirical model was calibrated on a single basin, and ET from the tallest plants, which covered only about 10 percent of the area but exerted outsized influence, required extrapolation beyond the height of the flux chambers. Microclimatic variables other than solar radiation were assumed uniform across the basin, and lacking inflow data, the team could not compute ET’s share of total stormwater removal. Still, the basin contains many features common to bioretention design everywhere, and the growing availability of high-resolution LiDAR, drone imagery, and weather data makes locally calibrated corrections increasingly feasible. As cities pour billions into green stormwater infrastructure to combat sewer overflows and climate-driven downpours, this study delivers a refreshingly actionable insight: the atmosphere’s share of the water budget is measurable, modelable, and, with a few well-chosen coefficients, finally predictable.

Subject of Research: Spatial heterogeneity in evapotranspiration estimates for urban bioretention basins

Article Title: Incorporating spatial heterogeneity into evapotranspiration estimates for bioretention basins

Article References: Caplan, J. S., Bouda, M., Salisbury, A. B., Alonzo, M., Nyquist, J. E., Toran, L., & Eisenman, S. W. (2026). Incorporating spatial heterogeneity into evapotranspiration estimates for bioretention basins. Hydrology and Earth System Sciences, 30(18), 6075-6093. https://doi.org/10.5194/hess-30-6075-2026

Image Credits: AI Generated

DOI: 10.5194/hess-30-6075-2026

Keywords: evapotranspiration, bioretention basins, green stormwater infrastructure, Penman-Monteith, spatial heterogeneity, urban hydrology, flux chambers, LiDAR, shade coefficients, landscape coefficients, Philadelphia, stormwater management

The research, led by Xinyu Song and colleagues working across institutions in China and the United States, tackles a problem that has long hampered precision oncology: most studies of drug resistance rely on samples taken from different patients at different stages, making it impossible to distinguish true temporal evolution from background variation between individuals. What makes this study remarkable is its longitudinal design. The team assembled a cohort of six patients with EGFR-mutant NSCLC who underwent serial re-biopsies at three critical time points: at baseline before any TKI treatment, after developing resistance to a first-generation EGFR-TKI, and after subsequently failing a third-generation EGFR-TKI. These consecutive samples from the same patients are extraordinarily rare, because repeated lung biopsies are invasive and many patients decline them, yet they are the only way to watch the same tumor evolve in real time within the same body.

To interrogate these precious samples, the researchers employed GeoMx Digital Spatial Profiling, a technology that bridges the gap between bulk sequencing, which averages signals across an entire tissue piece, and single-cell methods, which lose the architectural context of where each cell sits. With GeoMx, the team segmented each biopsy into tumor-enriched and stroma-enriched areas of interest, effectively separating the malignant cells themselves from the surrounding supportive tissue of fibroblasts, immune cells, and blood vessels. For each compartment, they quantified paired RNA and protein profiles, allowing them to ask not only which genes were being transcribed but also which proteins were actually being produced. This dual-layer approach matters because gene expression and protein abundance do not always correspond, and drugs typically target proteins rather than transcripts.

The first major finding was a confirmation that the spatial segmentation was biologically meaningful: tumor-enriched and stroma-enriched compartments proved consistently distinct from one another, both at the transcriptional and the protein level. This might sound like an expected result, but it validates the analytical framework and establishes that the technology reliably captures genuine biological differences rather than technical noise. More intriguing was what happened over time. The researchers found that spatial heterogeneity, the degree to which different regions of the tumor differ from one another, remained relatively stable as patients moved from baseline to first-generation TKI resistance. But between first-generation resistance and third-generation resistance, heterogeneity increased substantially, and the stroma-enriched compartments accounted for most of that expansion. In other words, the ecosystem surrounding the tumor, rather than the tumor cells alone, appears to become increasingly diverse and disorganized as the disease advances through sequential therapies.

This stromal expansion of heterogeneity carries important implications. The tumor microenvironment is not a passive backdrop; it actively shapes drug penetration, immune surveillance, and the selective pressures that drive resistance. If the stroma becomes progressively more heterogeneous after each line of therapy, it may create an increasingly complex landscape in which resistant clones can find refuge and in which immune cells become progressively less capable of recognizing and attacking malignant cells. The finding also suggests that therapeutic strategies targeting the microenvironment, such as stromal modulation or combinations with immunotherapy, might need to be timed differently depending on how many lines of TKI treatment a patient has already received.

The immune story that emerged from the longitudinal analysis was equally striking and followed what the authors describe as a biphasic trajectory. Early remodeling, occurring as tumors transitioned from baseline to first-generation TKI resistance, featured a loss of T cell activation programs, reduced neutrophil signatures, increased myeloid remodeling, and impairment of antigen presentation. Each of these changes points in the same troubling direction: toward an immune system that is progressively less able to mount an effective response against the tumor. T cell activation is the engine of anti-cancer immunity, antigen presentation is the mechanism by which tumor cells display their abnormal proteins for immune inspection, and myeloid cells, when remodeled toward immunosuppressive states, can actively suppress the immune response. The early phase of resistance thus appears to involve a coordinated dismantling of the immune recognition machinery.

The late phase of the trajectory, from first-generation to third-generation TKI resistance, was characterized instead by stromal checkpoint reprogramming, a shift in the inhibitory signaling molecules expressed within the supportive tissue surrounding the tumor. This temporal separation is significant because it suggests that resistance is not a single event but a staged process, with different biological mechanisms dominating at different points along the treatment course. A therapy designed to counteract early myeloid remodeling would need to be deployed at a different time than one aimed at late stromal checkpoint changes. The concept of stage-dependent immune remodeling, demonstrated here at spatial resolution within individual patients, provides a template for thinking about resistance as a choreographed sequence rather than a monolithic switch.

Perhaps the most tantalizing observation came from comparing tumors that later acquired the T790M mutation with those that did not. T790M is the classic gatekeeper mutation that renders first-generation EGFR-TKIs ineffective and is the very reason third-generation inhibitors such as osimertinib were developed. The researchers found that early spatial remodeling patterns showed potential differences between samples that subsequently acquired T790M and those that resisted through other mechanisms. While the authors are careful to frame this as preliminary and hypothesis-generating rather than definitive, the implication is profound: the microenvironmental changes visible early in treatment might eventually serve as predictive markers, telling clinicians which resistance pathway a tumor is likely to take before it does so. Such foresight could allow therapeutic interventions to be staged preemptively rather than reactively.

It is important to appreciate the limitations that the authors themselves acknowledge. With six patients and fifteen samples, this is an exploratory study, and its explicit purpose is to generate testable, time-aware hypotheses rather than to establish clinical guidelines. Rare consecutive re-biopsy cohorts are inherently small, and the findings will need validation in larger, independent datasets before they can influence practice. Nevertheless, the study provides something genuinely new: a preliminary spatially resolved, longitudinally matched atlas of how the tumor microenvironment changes across sequential EGFR-TKI therapies, capturing both RNA and protein layers within anatomically defined compartments. As spatial omics technologies mature and re-biopsy practices expand, atlases of this kind could become the foundation for a new generation of treatment strategies that anticipate, rather than merely react to, the evolutionary maneuvers of resistant tumors.

Subject of Research: Tumor microenvironment remodeling during sequential EGFR-TKI resistance in EGFR-mutant non-small-cell lung cancer

Article Title: Longitudinal spatial multi-omics delineates tumor microenvironment remodeling across sequential EGFR-TKIs in EGFR-mutant NSCLC

Article References: Song, X., Pan, C., Li, Z., Zhou, X., Xiong, A., Zhao, Z., Luo, H., & Zhou, C. (2026). Longitudinal spatial multi-omics delineates tumor microenvironment remodeling across sequential EGFR-TKIs in EGFR-mutant NSCLC. Journal of Translational Medicine, 24(1), Article 1212. https://doi.org/10.1186/s12967-026-08897-2

Image Credits: AI Generated

DOI: 10.1186/s12967-026-08897-2

Keywords: NSCLC, EGFR-TKIs, drug resistance, tumor microenvironment, spatial omics, digital spatial profiling, longitudinal re-biopsy, spatial heterogeneity, T790M, immune remodeling, tumor stroma, precision oncology

The research team, led by Yunyi Gu and Chengjin Chu of Sun Yat-Sen University together with colleagues from the Guangdong Academy of Forestry, Temple University, and the Guangdong Chebaling National Nature Reserve, worked in the Chebaling forest dynamic plot, a 20-hectare tract of protected mid-subtropical evergreen forest. Rather than sampling soil at a few scattered points, the researchers deployed high-resolution field sensing at 2,500 locations across the plot, capturing a detailed map of how sodium and other soil physicochemical properties vary from meter to meter. This dense spatial sampling was essential, because the team suspected that the fine-grained mosaic of soil conditions, not just their average values, might be what structures the forest above.

That suspicion proved well founded. The sodium data revealed pronounced spatial heterogeneity, with the variance in soil sodium peaking at a distance of roughly 200 meters, indicating that the forest floor is divided into broad patches of relatively sodium-rich and sodium-poor ground. When the researchers compared sodium against every other measured soil variable, sodium showed the strongest negative correlation with total tree species richness. In other words, the places in the forest where the soil carried the most sodium were the places hosting the fewest tree species, a pattern consistent with the osmotic and ionic stress that excess sodium imposes on plant cells.

The mechanism behind that stress is well understood at the physiological level. Sodium ions compete with potassium ions for uptake sites on root membranes, disrupting the enzyme function, stomatal regulation, and osmotic balance that potassium normally supports. At elevated concentrations, sodium accumulation in plant tissues causes ion toxicity, while the osmotic effect of salt in the soil solution makes it harder for roots to extract water. For most non-halophytic plants, sodium is at best a questionable nutrient and at worst a potent toxicant, which is why salinity is a major constraint on agriculture worldwide. What has been far less clear is whether these effects matter in humid, leached forest soils far from any coast, where salinity has traditionally been assumed to be a non-issue.

The subtropical forest data suggest that assumption deserves revision, but with a crucial twist: the impact of sodium depends on which fungal partners a tree relies on. The team classified every tree species in the plot by mycorrhizal type, separating arbuscular mycorrhizal (AM) species, whose roots are colonized by fungi that penetrate root cells, from ectomycorrhizal (EcM) species, whose fungal partners form a sheath around the roots, along with a residual category of other species. When the spatial correlations were computed separately for each group, a striking asymmetry emerged. Sodium was significantly and negatively correlated with richness for AM trees and for the other category, but the relationship vanished entirely for EcM trees.

This mycorrhizal split is more than a statistical curiosity. AM and EcM fungi represent two of the most ancient and widespread symbioses in terrestrial ecosystems, and they differ profoundly in how they acquire nutrients. Ectomycorrhizal fungi can enzymatically mine organic nitrogen and phosphorus directly from soil organic matter, effectively insulating their hosts from some forms of soil chemical stress, while arbuscular mycorrhizal fungi depend more heavily on the soil solution for nutrient uptake. That difference in foraging strategy could plausibly explain why EcM-associated trees appear buffered against sodium-driven filtering, whereas AM-associated species, which dominate the species pool of subtropical Chinese forests, bear the brunt of it. The result aligns with a growing body of work showing that mycorrhizal type shapes how trees respond to soil nutrient availability and can even influence soil carbon storage and global plant biogeography.

To move beyond simple correlations, the researchers turned to random forest regression, a machine-learning approach that quantifies the relative importance of many predictor variables simultaneously. The models delivered a verdict that few forest ecologists would have predicted a decade ago: soil sodium was the single most important predictor of total tree richness, of AM tree richness, and of richness in the other category, outranking soil temperature, pH, and the macronutrients that usually headline such analyses. For EcM species, by contrast, soil temperature and pH emerged as the primary correlates, reinforcing the picture of two functionally distinct tree guilds responding to different environmental filters within the same forest.

Tree size added a second layer of contingency. When the community was split by diameter at breast height, with small trees defined as stems under 10 centimeters and large trees at or above that threshold, sodium again dominated as the predictor for the small size class, while available phosphorus became the strongest correlate for large trees. This size-dependent pattern makes ecological sense. Young saplings and small stems have shallow, limited root systems and thin tissues, leaving them acutely vulnerable to osmotic stress and ion toxicity in sodium-enriched patches. Large, established trees, with deep roots and substantial nutrient reserves, are better equipped to tolerate sodium but remain constrained by phosphorus availability, a resource that is notoriously scarce and chemically locked away in the highly weathered soils of subtropical China. The study thus suggests that sodium acts as an early-life-stage filter, thinning the pool of recruits before trees ever reach the canopy.

The broader implications reach well beyond one forest plot. Ecologists have recently begun to argue that sodium deserves recognition as a ‘seventh macronutrient,’ one whose shortfall or excess ramifies through food webs, from ants and termites to decomposers and megaherbivores. At the same time, the anthropogenic salt cycle is accelerating: road de-icing, irrigation, mining, and land-use change are mobilizing salts into soils and waterways at rates that rival natural fluxes. If even humid subtropical forests, far from any obvious salinity source, show sodium-driven filtering of tree diversity, then rising salt loads elsewhere could quietly erode biodiversity in ecosystems currently considered safe. The Chebaling findings also suggest that conservation planning and restoration in species-rich forests may need to account for soil sodium maps, not just carbon and nutrient targets.

For the authors, the message is that the nutrient-centered framework governing plant-soil research needs expansion. Soil sodium, they conclude, is a critical yet underappreciated predictor of tree species richness, with effects strongly differentiated by mycorrhizal type and tree size. Integrating sodium into plant-soil studies, alongside the nitrogen and phosphorus paradigms of the past century, could sharpen predictions of how forest diversity will respond to environmental change. A forest’s hidden salt map, it turns out, may be quietly writing the rules of who lives where, one seedling at a time.

Subject of Research: Effects of soil sodium heterogeneity on tree species richness in a subtropical forest, mediated by mycorrhizal type and tree size class.

Article Title: Soil sodium differentially affects tree species richness depending on mycorrhizal type and size class in a subtropical forest

Article References: Gu, Y., Zhong, Y., Behm, J., Tan, H., Li, B., Liu, S., & Chu, C. (2026). Soil sodium differentially affects tree species richness depending on mycorrhizal type and size class in a subtropical forest. Plant and Soil. https://doi.org/10.1007/s11104-026-09074-6

Image Credits: AI Generated

DOI: 10.1007/s11104-026-09074-6

Keywords: soil sodium, tree species richness, subtropical forest, mycorrhizal type, arbuscular mycorrhiza, ectomycorrhiza, spatial heterogeneity, soil salinity, nutrient limitation, forest ecology, plant-soil interactions, random forest

Water regulation is one of the least visible yet most consequential services that ecosystems provide. Forested slopes intercept rainfall, slow surface runoff, and encourage infiltration into soils and aquifers, which in turn sustains baseflow in rivers during dry seasons and buffers downstream communities against floods and landslides. When forests are cleared for agriculture, plantations, or settlements, the landscape loses much of this sponge-like behavior. Rainfall that once soaked gradually into the ground instead races over compacted surfaces, carrying topsoil with it, swelling rivers quickly, and leaving less water stored for the months when it is needed most. Sulawesi, with its rugged central highlands, four peninsulas, and dense network of small watersheds, is particularly sensitive to these shifts because so many communities depend on short, steep catchments for irrigation and drinking water.

To capture this capacity spatially, the study employed a weighted spatial index built on a Simple Additive Weighting approach, a multi-criteria decision technique that combines several input layers into a single score. In this case, the layers included natural landscape characteristics, natural vegetation, and land cover data for each of the three reference years. Each land cover type was assigned a weight reflecting its contribution to water regulation, so that intact natural forest scored highest, while agricultural land and built-up areas scored progressively lower. The weighted layers were then summed to produce a continuous index of potential water regulation capacity, which was classified into categories ranging from very high to very low. This approach allowed the author to compare like with like across a 32-year span and to pinpoint where on the island the losses were concentrated.

The results reveal a clear geographic pattern of decline. Central Sulawesi and South Sulawesi contributed most strongly to the reduction in high-capacity areas, consistent with the intensive forest conversion and agricultural expansion that have characterized those provinces in recent decades. South Sulawesi, home to the densely populated southwestern lowlands and the Makassar urban corridor, has long experienced pressure from rice cultivation, cocoa and other perennial crops, and urban growth. Central Sulawesi has seen extensive land clearing linked to agriculture and, in some areas, plantation development. By contrast, Southeast Sulawesi and North Sulawesi showed partial recovery of high-capacity classes in the later period, suggesting that regrowth of secondary vegetation, protection efforts, or shifts in land management may have locally restored some regulatory function even as the island-wide trend remained negative.

Underlying these capacity changes is a fundamental shift in what covers the land. The study found that vegetated land cover types that strongly support water regulation declined from 75 percent of the island in 1990 to 54 percent in 2022. That is a loss of roughly a fifth of the island’s vegetated surface in a single generation, largely driven by conversion to agricultural and built-up areas. The magnitude of this transition matters because the relationship between vegetation and water regulation is not linear at the margins. Small losses of forest in headwater zones can disproportionately affect downstream flow regimes, and once soils are degraded, the hydrological functions they supported can take decades to recover even if vegetation returns.

The spatial heterogeneity that the study documents is itself scientifically important. Sulawesi is not losing its water regulation capacity uniformly; it is losing it in patches, corridors, and hotspots that correspond to provincial boundaries, watershed networks, and the geography of agricultural frontiers. This heterogeneity is precisely what makes the mapping approach valuable for policy. Rather than treating the island as a single unit with a single trend, the weighted index produces provincial- and watershed-level hotspots where the decline is most acute, and therefore where interventions such as forest protection, restoration, and ecosystem service-based zoning could be prioritized with the greatest expected benefit. In practical terms, a watershed whose index has slipped from high to moderate is a candidate for targeted restoration of riparian buffers and headwater forest, while a province showing widespread contraction of high-capacity classes may warrant stronger land use controls.

The study’s methodology also illustrates a broader trend in ecosystem service science: the use of relatively simple, transparent, and repeatable indices rather than data-hungry process models. Complex hydrological simulations can capture the physics of rainfall-runoff partitioning in detail, but they demand extensive calibration data that are often unavailable in tropical regions. A Simple Additive Weighting index, by contrast, requires only consistent land cover and landscape datasets, which are increasingly available from national mapping agencies and satellite programs. The trade-off is that the index measures potential capacity rather than actual water fluxes; it indicates where the landscape is structurally capable of supporting rainfall retention, baseflow, and watershed resilience, not how many cubic meters of water are being stored in any given year. For regional planning purposes, however, this potential-based framing is often exactly what is needed, because it identifies where the hydrological infrastructure of the landscape is being dismantled before the consequences appear in stream gauges.

The findings sit within a wider body of evidence on deforestation in Indonesia and its consequences. Sulawesi has been identified in previous research as a priority island for conservation because of its exceptional endemism and the rapid loss of its forests, with documented impacts on primate habitat and unique ecosystems. Studies of deforestation patterns in Sulawesi between 1990 and 2018 have shown persistent pressure on forest cover, and work on mangrove decline in South Sulawesi has highlighted how coastal ecosystem services have been degraded alongside terrestrial ones. The new analysis extends this picture into the hydrological domain, showing that the same land conversion processes that threaten biodiversity are also dismantling the island’s capacity to regulate water, with direct implications for flood risk, dry-season water supply, and agricultural stability for millions of people.

What makes the study particularly timely is the policy context. Indonesia has experimented with forest moratoria, peatland protection, and payment-for-ecosystem-service schemes, and researchers have debated how cost-effective these instruments are in reducing emissions and protecting services. Mapping exercises like this one provide the spatial evidence base that such policies require. If restoration funds are limited, they should flow to the watersheds where the loss of regulation capacity is steepest and where downstream exposure is greatest. The study’s outputs, which identify provincial- and watershed-level hotspots, are designed to support exactly that kind of prioritization, and the author notes that they could inform ecosystem service-based zoning as part of provincial spatial planning.

There are also lessons here for how tropical islands elsewhere are monitored. Sulawesi’s story, a 15-percentage-point decline in high-capacity land over three decades, driven overwhelmingly by the replacement of vegetation with agriculture and settlements, is a pattern that repeats across Southeast Asia, Madagascar, the Caribbean, and many other island systems where steep terrain concentrates both water services and land use pressure. The demonstration that a weighted spatial index built from freely available land cover data can track these changes consistently over more than three decades offers a template for other data-limited regions. As climate change intensifies both droughts and extreme rainfall in the tropics, the landscapes that quietly regulate water are becoming more valuable, not less. Sulawesi’s shrinking high-capacity core is a quantitative reminder that this natural infrastructure is finite, spatially specific, and, once lost, slow and expensive to rebuild.

Subject of Research: Land cover change and water regulation ecosystem services on Sulawesi Island, Indonesia

Article Title: Spatial heterogeneity of landscapes and land cover in Sulawesi Island in relation to water regulation ecosystem services

Article References: Nursaputra, M. (2026). Spatial heterogeneity of landscapes and land cover in Sulawesi Island in relation to water regulation ecosystem services. Environmental Monitoring and Assessment, 198(10), Article 1086. https://doi.org/10.1007/s10661-026-15892-z

Image Credits: AI Generated

DOI: 10.1007/s10661-026-15892-z

Keywords: ecosystem services, water regulation, land cover change, spatial heterogeneity, Sulawesi Island, deforestation, Simple Additive Weighting, watershed resilience, forest restoration, Indonesia, landscape ecology, remote sensing

The research team, led by Lantian Zhang of the Institute of Geographic Sciences and Natural Resources Research with collaborators including Mark W. Rosenberg of Queen’s University, assembled an unusually ambitious dataset covering 1,376 counties across China, of which 379 are classified as Kashin-Beck disease endemic. Rather than treating the disease itself as the focus, the researchers used these endemic counties as a natural laboratory for studying healthy ageing under conditions of ecological fragility and long-term socioeconomic disadvantage. Their outcome measure was old-age healthy life expectancy at age 60, a metric endorsed by the World Health Organization that combines mortality data with the prevalence of good health to estimate how many years a 60-year-old can expect to live in reasonable health.

Calculating this measure for more than a thousand counties across fifteen years was itself a technical feat. The team drew on census microdata from 2005, 2010, and 2020 provided through the Peking University Research Data Center, applying the Sullivan method to construct life tables from age 60 onward. Self-rated health, in which individuals assess their own overall condition, served as the health indicator, with people classified as healthy or basically healthy counted among the healthy states. Because no nationally representative mortality survey was conducted in 2005, the researchers estimated that year’s mortality using the 2000 national survey, refined through small-area estimation models that incorporated county-level covariates such as education rates, ethnic composition, income, and population density, then scaled the results to match provincial statistics.

On the explanatory side, the researchers built a multidimensional sustainability indicator system aligned with the United Nations Sustainable Development Goals, spanning the classic three pillars of sustainable development. The ecological dimension included forest coverage, vegetation indices measured by satellite, soil conservation, air pollutants such as fine particulate matter and ozone, and carbon emissions. The economic dimension captured income, employment, industrial structure, and financial resources. The social dimension encompassed safe drinking water, sanitation, housing, road density, education, and critically, the density of health facilities, health technicians, and hospital beds. Each indicator was standardized, and dimension-specific indices were constructed using a hybrid weighting scheme that combined expert-driven analytic hierarchy process weights with data-driven entropy weights in a 70:30 ratio, a balance the team tested for robustness against alternative weightings.

The analytical engine of the study was a hierarchical Bayesian framework implemented through Integrated Nested Laplace Approximation, a computational approach that has revolutionized spatial statistics by making complex models tractable. The researchers started with a baseline model decomposing healthy life expectancy into structured and unstructured spatial effects, a common temporal trend, and region-specific deviations from that trend. They then layered on increasingly sophisticated specifications, culminating in a spatiotemporal varying coefficient model in which the regression coefficients themselves shift across both space and time, and a further extension allowing county-by-time interactions. Model selection relied on the deviance information criterion, the Watanabe-Akaike information criterion, and cross-validation, while a spatiotemporal variance partitioning index quantified how much of the variation in healthy life expectancy each sustainability dimension explained.

The headline finding is stark: social sustainability matters most, and it matters most precisely where societies are weakest. Across all counties, a one-unit increase in the standardized social sustainability index was associated with an 11.66-year difference in healthy life expectancy at 60, with a 95 percent credible interval of 10.85 to 12.48 years. Ecological and economic sustainability showed positive but far smaller associations of roughly 2.3 and 2.2 years respectively. Yet the averages conceal dramatic heterogeneity. In Kashin-Beck endemic counties, the social sustainability coefficient rose to 13.56 years, while economic sustainability showed no detectable association at all. In non-endemic counties, the pattern reversed: economic sustainability was positively linked to healthy ageing, and social effects were somewhat weaker. On the Qinghai-Tibet Plateau, endemic counties showed the strongest ecological associations anywhere in the study, at 6.15 years, suggesting that in harsh high-altitude environments the natural surroundings themselves become a meaningful determinant of healthy old age.

Perhaps the most striking result is temporal. When the researchers decomposed the variation in healthy life expectancy, temporal dynamics of the sustainability indicators explained more than 90 percent of the variance in both endemic and non-endemic areas, dwarfing spatial effects, with model residuals accounting for less than 5 percent. This dominance of change over time reflects the sweeping nationwide improvements in living conditions and healthcare that China experienced between 2005 and 2020, including the 2009 New Healthcare Reform that expanded insurance coverage and strengthened primary care. In the Kashin-Beck affected areas specifically, decades of government intervention, including selenium supplementation, population relocation, and ecological compensation, appear to have paid dividends as the cohorts affected by earlier epidemics reached old age. The study also documented a narrowing of regional health inequalities, with non-endemic counties averaging 15.96 healthy years at age 60 compared with 15.15 years in endemic counties, a significant but not overwhelming gap.

The geography of these associations shifted in revealing ways. In 2005, the strongest links between social sustainability and healthy ageing were concentrated on the Qinghai-Tibet Plateau, where geographic inaccessibility and weak infrastructure made clinics, sanitation, and schools the binding constraints on wellbeing. By 2020, those strongest associations had migrated toward more economically developed regions, a redistribution the authors interpret through the lens of evolving policy: as basic healthcare shortages were alleviated in poorer areas, the marginal importance of social conditions migrated to wealthier regions grappling with chronic disease burdens that demand long-term, integrated health and social care. Meanwhile, the influence of ecological sustainability grew steadily over the fifteen-year window, hinting that as other determinants improve, environmental quality becomes an increasingly relevant pillar of healthy ageing.

For policymakers, the message is that one size does not fit all. In regions with fragile socioeconomic foundations, investments in social infrastructure, healthcare accessibility, sanitation, and education offer the greatest returns for healthy ageing, while sustained economic development and environmental governance matter over the longer horizon. The authors are candid about limitations: self-rated health captures a subjective dimension of wellbeing, the Sullivan method relies on cross-sectional data and cannot disentangle birth cohort effects, and as an ecological analysis the study describes regional patterns rather than individual-level mechanisms. Still, by demonstrating that the sustainability-health relationship is neither fixed in place nor uniform across the landscape, the study delivers a methodological lesson as important as its substantive one. Static, aspatial analyses of healthy ageing, it suggests, are not merely incomplete; they can obscure the very dynamics that determine whether the world’s fastest-ageing regions, many of them its most environmentally fragile, can convert rising life expectancy into genuinely healthy longevity.

Subject of Research: Spatiotemporal associations between ecological, economic, and social sustainability and old-age healthy life expectancy in ecologically fragile regions of China

Article Title: Spatiotemporal effects of sustainability dimensions on old-age healthy life expectancy in ecologically fragile regions

Article References: Zhang, L., Gu, L., Yang, L., Li, H., Wang, L., & Rosenberg, M. W. (2026). Spatiotemporal effects of sustainability dimensions on old-age healthy life expectancy in ecologically fragile regions. Environmental and Sustainability Indicators, 32, Article 101546. https://doi.org/10.1016/j.indic.2026.101546

Image Credits: AI Generated

DOI: 10.1016/j.indic.2026.101546

Keywords: healthy life expectancy, population ageing, sustainability indicators, Kashin-Beck disease, Bayesian spatiotemporal modeling, Qinghai-Tibet Plateau, Loess Plateau, social sustainability, Sullivan method, spatial heterogeneity, Sustainable Development Goals, China