Methane fluxes (FCH4) vary significantly across wetland ecosystems due to complex mechanisms, challenging accurate estimations. The interactions among environmental drivers, while crucial in regulating FCH4, have not been well understood. Here, the interactive effects of six environmental drivers on FCH4 were first analyzed using 396,322 half-hourly measurements from 22 sites across various wetland types and climate zones. Results reveal that soil temperature, latent heat turbulent flux, and ecosystem respiration primarily exerted direct effects on FCH4, while air temperature and gross primary productivity mainly exerted indirect effects by interacting with other drivers. Significant spatial variability in FCH4 regulatory mechanisms was highlighted, with different drivers demonstrated varying direct, indirect, and total effects among sites. This spatial variability was then linked to site-specific annual-average air temperature (17.7%) and water table (9.0%) conditions, allowing the categorization of CH4 sources into four groups with identified critical drivers. An improved estimation approach using a random forest model with three critical drivers was consequently proposed, offering accurate FCH4 predictions with fewer input requirements. By explicitly accounting for environmental interactions and interpreting spatial variability, this study enhances our understanding of the mechanisms regulating CH4 emissions, contributing to more efficient modeling and estimation of wetland FCH4.

Over the past century, water temperatures in many streams across the Pacific Northwest (PNW) have steadily risen, shrinking endangered salmonid habitats. The warming of PNW stream reaches can be further accelerated by wildfires burning forest stands that provide shade to streams. However, previous research on the effect of wildfires on stream water temperatures has focused on individual streams or burn events, limiting our understanding of the diversity in post-fire thermal responses across PNW streams. To bridge this knowledge gap, we assessed the impact of wildfires on daily summer water temperatures across 31 PNW stream sites, where 10-100% of their riparian area burned. To ensure robustness of our results, we employed multiple approaches to characterize and quantify fire effects on post-fire stream water temperature changes. Averaged across the 31 burned sites, wildfires corresponded to a 0.3 - 1°C increase in daily summer water temperatures over the subsequent three years. Nonetheless, post-fire summer thermal responses displayed extensive heterogeneity across burned sites where the likelihood and rate of a post-fire summer water temperature warming was higher for stream sites with greater proportion of their riparian area burned under high severity. Also, watershed features such as basin area, post-fire weather, bedrock permeability, pre-fire riparian forest cover, and winter snowpack depth were identified as strong predictors of the post-fire summer water temperature responses across burned sites. Our study offers a multi-site perspective on the effect of wildfires on summer stream temperatures in the PNW, providing insights that can inform freshwater management efforts beyond individual streams and basins.

Humans have profoundly altered phosphorus (P) cycling across scales. Agriculturally driven changes (e.g., excessive P-fertilization and manure addition), in particular, have resulted in pronounced P accumulations in soils, often known as \"soil legacy P.\" These legacy P reserves serve as persistent and long-term nonpoint sources, inducing downstream eutrophication and ecosystem services degradation. While there is considerable scientific and policy interest in legacy P, its fine-scale spatial heterogeneity, underlying drivers, and scales of variance remain unclear. Here we present an extensive field sampling (150-m interval grid) and analysis of 1438 surface soils (0-15 cm) in 2020 for two typical subtropical grassland types managed for livestock production: Intensively managed (IM) and Semi-natural (SN) pastures. We ask the following questions: (1) What is the spatial variability, and are there hotspots of soil legacy P? (2) Does soil legacy P vary primarily within pastures, among pastures, or between pasture types? (3) How does soil legacy P relate to pasture management intensity, soil and geographic characteristics? and (4) What is the relationship between soil legacy P and aboveground plant tissue P concentration? Our results showed that three measurements of soil legacy P (total P, Mehlich-1, and Mehlich-3 extractable P representing labile P pools) varied substantially across the landscape. Spatial autoregressive models revealed that soil organic matter, pH, available Fe and Al, elevation, and pasture management intensity were crucial predictors for spatial patterns of soil P, although models were more reliable for predicting total P (68.9%) than labile P. Our analysis further demonstrated that total variance in soil legacy P was greater in IM than SN pastures, and intensified pasture management rescaled spatial patterns of soil legacy P. In particular, after controlling for sample size, soil P was extremely variable at small scales, with variance diminished as spatial scale increased. Our results suggest that broad pasture- or farm-level best management practices may be limited and less efficient, especially for more IM pastures. Rather, management to curtail soil legacy P and mitigate P loading and losses should be implemented at fine scales designed to target spatially distinct P hotspots across the landscape.

PMMoV has been widely used to normalize the concentration of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA, influenza, and respiratory syncytial virus (RSV) to account for variations in the fecal content of wastewater. PMMoV is also used as an internal RNA recovery control for wastewater-based epidemiology (WBE) tests. While potentially useful for the interpretation of WBE data, previous studies have suggested that PMMoV concentration can be affected by various physico-chemical characteristics of wastewater. There is also the possibility that laboratory methods, particularly the variability in centrifugation steps to remove supernatant from pellets can cause PMMoV variability. The goal of this study is to improve our understanding of the main drivers of PMMoV variability by assessing the relationship between PMMoV concentration, the physico-chemical characteristics of wastewater, and the methodological approach for concentrating wastewater samples. We analyzed 24-hour composite wastewater samples collected from the influent stream of three wastewater treatment plants (WWTPs) located in the City of Toronto, Ontario, Canada. Samples were collected 3 to 5 times per week starting from the beginning of March 2021 to mid-July 2023. The influent flow rate was used to partition the data into wet and dry weather conditions. Physico-chemical characteristics (e.g., total suspended solids (TSS), biological oxygen demand (BOD), alkalinity, electrical conductivity (EC), and ammonia (NH3)) of the raw wastewater were measured, and PMMoV was quantified. Spatial and temporal variability of PMMoV was observed throughout the study period. PMMoV concentration was significantly higher during dry weather conditions. Multiple linear regression analysis demonstrates that the number and type of physico-chemical parameters that drive PMMoV variability are site-specific, but overall BOD and alkalinity were the most important predictors. Differences in PMMoV concentration for a single WWTP between two different laboratory methods, along with a weak correlation between pellet mass and TSS using one method may indicate that differences in sample concentration and subjective subsampling bias could alter viral recovery and introduce variability to the data.

Human-driven multiple pressures impact freshwater ecosystems worldwide, reducing biodiversity, and impacting ecosystem functioning and services provided to human societies. Multi-metric indices (MMIs) are suitable tools for tracking the effects of anthropogenic pressures on freshwater ecosystems because they incorporate various biological metrics responding to multiple pressures at different levels of biological organization. However, the performance and applicability of MMIs depend on their metrics\' selection and their calibration against natural environmental gradients. In this study, we aimed to unravel i) how incorporating functional trait-based metrics affects the performance of MMIs, ii) how disentangling the natural environmental gradients from anthropogenic pressures effects affects the performance of MMIs, and iii) how the performance of MMIs developed using a metric performance-driven approach compares with MMIs developed using an index performance-driven approach. We carried out a field survey measuring abiotic and biotic variables at 53 sites in the Karun River basin (Iran) in 2018. For functional trait-based metrics, we used 15 macroinvertebrate traits and calculated community-weighted mean trait values and functional diversity indices. We used random forest modeling to account for the effect of natural environmental gradients on each metric. Based on our results, incorporating functional traits increased the MMI performance significantly and facilitated ecological interpretation of MMIs. Both taxonomic and functional components of macroinvertebrate assemblages co-varied strongly with natural environmental gradients, and accounting for these covariations improved the performance of MMIs. Finally, we found that index performance-driven MMIs performed better in terms of precision, bias, sensitivity, and responsiveness than metric performance-driven MMIs.

Conventional monitoring and mapping approaches are laborious, expensive, and time-consuming because they need a large number of data and consequently extensive sampling and experimental operations. Therefore, due to the growing concern about the potential of contamination of soils and agricultural products with heavy metals (HMs), a field experiment was conducted on 77 farm lands in an area of 2300 ha in the southeast of Shiraz (Iran) to investigate the source of metal contamination in the soils and vegetables and to model spatial distribution of HMs (iron, Fe; manganese, Mn; copper, Cu; zinc, Zn; cadmium, Cd; nickel, Ni, and lead, Pb) over the region using geographic information system (GIS) and geostatistical (Ordinary Kriging, OK) approaches and compare the results with deterministic approaches (Inverse Distance Weighting, IDW with different weighting power). Furthermore, some ecological and health risks indices including Pollution index (PI), Nemerow integrated pollution index (NIPI), pollution load index (PLI), degree of contamination (Cdeg), modified contamination degree (mCd), PIaverage and PIvector for soil quality, multi-element contamination (MEC), the probability of toxicity (MERMQ), the potential ecological index (RI), total hazard index (THI) and total carcinogenic risk index (TCR) based on ingestion, inhalation, and dermal exposure pathways for adults and children respectively for analyzing the noncarcinogenic and carcinogenic risks were calculated. Experimental semivariogram of the mentioned HMs were calculated and theoretical models (i.e., exponential, spherical, Gaussian, and linear models) were fitted in order to model their spatial structures and to investigate the most representative models. Moreover, principal component analysis (PCA) and cluster analysis (CA) were used to identify sources of HMs in the soils. Results showed that IDW method was more efficient than the OK approach to estimate the properties and HMs contents in the soils and plants. The estimated daily intake of metals (DIM) values of Pb and Ni exceeded their safe limits. In addition, Cd was the main element responsible for ecological risk. The PIave and PIvector indices showed that soil quality in the study area is not suitable. According to mCd values, the soils classified as ultra-high contaminated for Cu and Cd, extremely high for Zn and Pb, very high, high, and very low degree of contamination for Ni, Mn, and Fe, respectively. 36, 60, and 4 % of the sampling sites had high, medium, and low risk levels with 49, 21, and 9 % probability of toxicity, respectively. The maximum health risk index (HRI) value of 20.42 with extremely high risk for children was obtained for Ni and the HI for adults and children were 0.22 and 1.55, respectively. The THI values of Pb and Cd were the highest compared to the other HMs studied, revealing a possible non-cancer risk in children associated with exposure to these metals. The routes of exposure with the greatest influence on the THI and TCR indices were in the order of ingestion > inhalation > dermal. Therefore, ingestion, as the main route of exposure, is the route of greatest contribution to health risks. PCA analysis revealed that Fe, Mn, Cu, and Ni may originate from natural sources, while Fe was appeared to be controlled by fertilizer, and Cu primarily coming from pesticide, while Cd and Pb were mainly associated with the anthropogenic contamination, atmospheric depositions, and terrific in the urban soils. While, Zn mainly originated from fertilization. Findings are vital for developing remediation approaches for controlling the contaminants distribution as well as for monitoring and mapping the quality and health of soil resources.

The contents of eight carbonaceous subfractions were determined by simultaneously collecting PM2.5 samples from four sites in different functional areas of Tianjin in 2021. The results showed that the organic carbon (OC) concentration was 3.7 μg·m-3 to 4.4 μg·m-3, and the elemental carbon (EC) concentration was 1.6 μg·m-3 to 1.7 μg·m-3, with the highest OC concentration in the central urban area. There was no significant difference in EC concentration. The concentration of PM2.5 showed the distribution characteristics of the surrounding city>central city>peripheral area. The OC/EC minimum ratio method was used to estimate the concentrations of secondary organic carbon (SOC) in PM2.5, and the results showed that the secondary pollution was more prominent in the surrounding city, with SOC accounting for 48.8%. The correlation between carbon subcomponents in each functional area showed the characteristics of the peripheral area>central area>surrounding area, all showing the strongest correlation between EC1 and OC2 and EC1 and OC4. By including the carbon component concentration into the positive definite matrix factorization (PMF) model for source apportionment, the results showed that road dust sources(9.7%-23.5%), coal-combustion sources (10.2%-13.3%), diesel vehicle exhaust (12.6%-20.2%)and gasoline vehicle exhaust (18.9%-38.8%)were the main sources of carbon components in PM2.5 in Tianjin. The pollution sources of carbon components were different in different functional areas, with the central city and peripheral areas mainly affected by gasoline vehicle exhaust; the surrounding city was more prominently affected by the secondary pollution and diesel vehicle exhaust.

Malaria ranks high among prevalent and ravaging infectious diseases in sub-Saharan Africa (SSA). The negative impacts, disease burden, and risk are higher among children and pregnant women as part of the most vulnerable groups to malaria in Nigeria. However, the burden of malaria is not even in space and time. This study explores the spatial variability of malaria prevalence among children under five years (U5) in medium-sized rapidly growing city of Akure, Nigeria using model-based geostatistical modeling (MBG) technique to predict U5 malaria burden at a 100 × 100 m grid, while the parameter estimation was done using Monte Carlo maximum likelihood method. The non-spatial logistic regression model shows that U5 malaria prevalence is significantly influenced by the usage of insecticide-treated nets-ITNs, window protection, and water source. Furthermore, the MBG model shows predicted U5 malaria prevalence in Akure is greater than 35% at certain locations while we were able to ascertain places with U5 prevalence > 10% (i.e. hotspots) using exceedance probability modelling which is a vital tool for policy development. The map provides place-based evidence on the spatial variation of U5 malaria in Akure, and direction on where intensified interventions are crucial for the reduction of U5 malaria burden and improvement of urban health in Akure, Nigeria.

This study addressed the scarcity of NH3 measurements in urban Europe and the diverse monitoring protocols, hindering direct data comparison. Sixty-nine datasets from Finland, France, Italy, Spain, and the UK across various site types, including industrial (IND, 8), traffic (TR, 12), urban (UB, 22), suburban (SUB, 12), and regional background (RB, 15), are analyzed to this study. Among these, 26 sites provided 5, or more, years of data for time series analysis. Despite varied protocols, necessitating future harmonization, the average NH3 concentration across sites reached 8.0 ± 8.9 μg/m3. Excluding farming/agricultural hotspots (FAHs), IND and TR sites had the highest concentrations (4.7 ± 3.2 and 4.5 ± 1.0 μg/m3), followed by UB, SUB, and RB sites (3.3 ± 1.5, 2.7 ± 1.3, and 1.0 ± 0.3 μg/m3, respectively) indicating that industrial, traffic, and other urban sources were primary contributors to NH3 outside FAH regions. When referring exclusively to the FAHs, concentrations ranged from 10.0 ± 2.3 to 15.6 ± 17.2 μg/m3, with the highest concentrations being reached in RB sites close to the farming and agricultural sources, and that, on average for FAHs there is a decreasing NH3 concentration gradient towards the city. Time trends showed that over half of the sites (18/26) observed statistically significant trends. Approximately 50 % of UB and TR sites showed a decreasing trend, while 30 % an increasing one. Meta-analysis revealed a small insignificant decreasing trend for non-FAH RB sites. In FAHs, there was a significant upward trend at a rate of 3.51[0.45,6.57]%/yr. Seasonal patterns of NH3 concentrations varied, with urban areas experiencing fluctuations influenced by surrounding emissions, particularly in FAHs. Diel variation showed differing patterns at urban monitoring sites, all with higher daytime concentrations, but with variations in peak times depending on major emission sources and meteorological patterns. These results offer valuable insights into the spatio-temporal patterns of gas-phase NH3 concentrations in urban Europe, contributing to future efforts in benchmarking NH3 pollution control in urban areas.

A major motivation for multiple atmospheric probe measurements at Uranus is the understanding of dynamic processes that create and maintain spatial variation in thermal structure, composition, and horizontal winds. But origin questions-regarding the planet\'s formation and evolution, and conditions in the protoplanetary disk-are also major science drivers for multiprobe exploration. Spatial variation in thermal structure reveals how the atmosphere transports heat from the interior, and measuring compositional variability in the atmosphere is key to ultimately gaining an understanding of the bulk abundances of several heavy elements. We review the current knowledge of spatial variability in Uranus\' atmosphere, and we outline how multiple probe exploration would advance our understanding of this variability. The other giant planets are discussed, both to connect multiprobe exploration of those atmospheres to open questions at Uranus, and to demonstrate how multiprobe exploration of Uranus itself is motivated by lessons learned about the spatial variation at Jupiter, Saturn, and Neptune. We outline the measurements of highest value from miniature secondary probes (which would complement more detailed investigation by a larger flagship probe), and present the path toward overcoming current challenges and uncertainties in areas including mission design, cost, trajectory, instrument maturity, power, and timeline.