Nitrogen availability and summer drought, but not N:P imbalance, drive carbon use efficiency of a Mediterranean tree-grass ecosystem.

All ecosystems contain both sources and sinks for atmospheric carbon (C). A change in their balance of net and gross ecosystem carbon uptake, ecosystem-scale carbon use efficiency (CUEECO), is a change in their ability to buffer climate change. However, anthropogenic nitrogen (N) deposition is increasing N availability, potentially shifting terrestrial ecosystem stoichiometry towards phosphorus (P) limitation. Depending on how gross primary production (GPP, plants alone) and ecosystem respiration (RECO, plants and heterotrophs) are limited by N, P or associated changes in other biogeochemical cycles, CUEECO may change. Seasonally, CUEECO also varies as the multiple processes that control GPP and respiration and their limitations shift in time. We worked in a Mediterranean tree-grass ecosystem (locally called 'dehesa') characterized by mild, wet winters and summer droughts. We examined CUEECO from eddy covariance fluxes over 6 years under control, +N and + NP fertilized treatments on three timescales: annual, seasonal (determined by vegetation phenological phases) and 14-day aggregations. Finer aggregation allowed consideration of responses to specific patterns in vegetation activity and meteorological conditions. We predicted that CUEECO should be increased by wetter conditions, and successively by N and NP fertilization. Milder and wetter years with proportionally longer growing seasons increased CUEECO, as did N fertilization, regardless of whether P was added. Using a generalized additive model, whole ecosystem phenological status and water deficit indicators, which both varied with treatment, were the main determinants of 14-day differences in CUEECO. The direction of water effects depended on the timescale considered and occurred alongside treatment-dependent water depletion. Overall, future regional trends of longer dry summers may push these systems towards lower CUEECO.

Macroscopic and molecular scale assessment of thermal effects on soil-water interactions of unsaturated diesel-contaminated soils.

The soil-water interactions of unsaturated diesel-contaminated soil are crucial for assessing pollution transport during thermal remediation. This paper aims to improve our understanding of this issue by measuring the matric suction of unsaturated contaminated kaolin and carrying out molecular dynamics simulations under thermal conditions. Results show that the increase in pollutant concentration could reduce the water retention capacity of diesel-contaminated kaolin due to changes in electrochemical properties and pore characteristics of samples, as well as a decrease in interfacial tension. On the other hand, pollutants formed a protective film on the kaolinite surface to act as a liquid bridge and prevent water loss at higher temperatures, as confirmed by Fourier transform infrared spectroscopy. With rising temperatures (50-60 °C), kaolin matric suction generally decreased with higher pollutant concentrations, but this trend was not very evident at lower pollution concentrations (0-10,000 mg/kg). In addition, molecular dynamics simulations were used to demonstrate the validity of these findings. The presence of pollutants might strengthen the interaction energy between kaolinite and water (for example, increasing from 276.52 kcal/mol (25 °C) and 267.95 kcal/mol (40 °C) at 8000 mg/kg to 296.54 kcal/mol (25 °C) and 292.46 kcal/mol (40 °C) at 10,000 mg/kg), thereby enhancing the water retention capacity of kaolin. In short, the study revealed that the coating of pollutants on kaolinite could act as a protective film, which binds water molecules through van der Waals and electric field forces and thereby reduces the sensitivity of water retention capacity to temperature.

Tree species explain only half of explained spatial variability in plant water sensitivity.

Spatiotemporal patterns of plant water uptake, loss, and storage exert a first-order control on photosynthesis and evapotranspiration. Many studies of plant responses to water stress have focused on differences between species because of their different stomatal closure, xylem conductance, and root traits. However, several other ecohydrological factors are also relevant, including soil hydraulics, topographically driven redistribution of water, plant adaptation to local climatic variations, and changes in vegetation density. Here, we seek to understand the relative importance of the dominant species for regional-scale variations in woody plant responses to water stress. We map plant water sensitivity (PWS) based on the response of remotely sensed live fuel moisture content to variations in hydrometeorology using an auto-regressive model. Live fuel moisture content dynamics are informative of PWS because they directly reflect vegetation water content and therefore patterns of plant water uptake and evapotranspiration. The PWS is studied using 21,455 wooded locations containing U.S. Forest Service Forest Inventory and Analysis plots across the western United States, where species cover is known and where a single species is locally dominant. Using a species-specific mean PWS value explains 23% of observed PWS variability. By contrast, a random forest driven by mean vegetation density, mean climate, soil properties, and topographic descriptors explains 43% of observed PWS variability. Thus, the dominant species explains only 53% (23% compared to 43%) of explainable variations in PWS. Mean climate and mean NDVI also exert significant influence on PWS. Our results suggest that studies of differences between species should explicitly consider the environments (climate, soil, topography) in which observations for each species are made, and whether those environments are representative of the entire species range.

Drivers of barrier island water-table fluctuations and groundwater salinization.

Barrier islands are threatened by climate change as sea-level rise and higher frequency storm surge lead to more flooding and saltwater intrusion. Vegetation plays a vital role in preventing erosion of barrier islands due to aeolian and hydrological forces. However, vegetation on barrier islands is threatened by rising water tables causing hypoxic conditions and storm-surge overwash introducing saline water to the root zone. To better protect barrier island ecosystems, it is critical to identify the relative influence of different hydrological drivers on water table elevation and salinity, and understand how this influence varies spatially and temporally. In this study, three barrier island sites were instrumented with groundwater wells monitoring water level and specific conductance. Using these data, a set of transfer function noise models were calibrated and used to determine the relative influence of hydrologic drivers including precipitation, evapotranspiration, bay and ocean water levels, and wave height on groundwater levels and specific conductance. We found that drivers of water-level change and specific conductance vary strongly among sites, depending primarily on the surface water connectivity and the geology of the island. Sites with close connection to inlets showed more salinization and responded to a larger number of drivers, while sites that were poorly connected to the ocean responded to fewer drivers.

Wheat gluten structure and (non-)covalent network formation during deep-fat frying.

The aim of this work was to investigate wheat gluten protein network structure throughout the deep-frying process and evaluate its contribution to frying-induced micro- and macrostructure development. Gluten polymerization, gluten-water interactions, and molecular mobility were assessed as a function of the deep-frying time (0 - 180 s) for gluten-water model systems of differing hydration levels (40 - 60 % moisture content). Results showed that gluten protein extractability decreased considerably upon deep frying (5 s) mainly due to glutenin polymerization by disulfide covalent cross-linking. Stronger gliadin and glutenin protein-protein interactions were attributed to the formation of covalent linkages and evaporation of water interacting with protein chains. Longer deep-frying (> 60 s) resulted in progressively lower protein extractabilities, mainly due to the loss in gliadin protein extractability, which was associated with gliadin co-polymerization with glutenin by thiol-disulfide exchange reactions. The mobility of gluten polymers was substantially reduced during deep-frying (based on the lower T2 relaxation time of the proton fraction representing the non-exchanging protons of gluten) and gluten proteins gradually transitioned from the rubbery to the glassy state (based on the increased area of said protons). The sample volume during deep-frying was strongly correlated to the reduced protein extractability (r = -0.792, p < 0.001) and T2 relaxation time of non-exchanging protons of gluten proteins (r = -0.866, p < 0.001) thus demonstrating that the extent of gluten structural expansion as a result of deep-frying is dictated both by the polymerization of proteins and the reduction in their molecular mobility.

Teasing out the Effects of Natural Stressors at Chemically Contaminated Sites.

Aquatic ecosystems are often impacted by a multitude of stressors, many of which are introduced by a combination of anthropogenic activities such as agricultural development, urbanization, damming, and industrial discharge. Determining the primary stressors responsible for ecological impairments at a site can be complex and challenging; however, it is crucial for making informed management decisions. Improper diagnosis of an impaired system can lead to misguided attempts at remediation, which can be both time consuming and costly. We focused on the development, implementation, and evaluation of methodologies that, in combination, allowed us to identify the primary stressors. These included a four-phase, weight-of-evidence (WOE) assessment including in situ Toxicity Identification and Evaluation (iTIE) testing, physicochemical and macrobenthos characterization, reciprocal sediment transplants, and laboratory and in situ toxicity testing. The contaminants of concern (COCs) at the site were elevated levels of ammonia, chloride, pH, and total dissolved solids in groundwater upwellings into a high-quality waterway. Reciprocal transplants of site sediments and nearby reference sediments and traditional benthic sampling showed impaired benthic indices and multiple stations around a contaminated industrial settling basin. Impaired stations had elevated COCs in groundwaters but exhibited a steep vertical concentration gradient, with concentrations decreasing near the sediment-surface water interface. We describe Phase 4 of the study, which focused on teasing out the role of dissolved oxygen sags in benthic macroinvertebrate responses. Extensive submerged and emergent macrophytes, algae, and cyanobacteria co-occurred at the impaired sites and increased throughout the summer. Laboratory testing suggested that ammonia and pH were possibly toxic at the sites, based on groundwater concentrations. The in situ toxicity testing, however, showed toxicity occurring even at stations with low levels of COCs concurrently with large diurnal fluxes in dissolved oxygen (DO). A final phase using a type of iTIE approach utilized limnocorrals with and without aeration and with in situ toxicity measures using Hyalella azteca. The Phase 4 assessment revealed that low DO levels were primarily responsible for impaired benthic communities, and COC upwellings were diluted at the sediment-water interface to nontoxic levels. These findings will allow for improved management decisions for more efficient and effective restoration activities. Environ Toxicol Chem 2024;43:1524-1536. © 2024 SETAC.

Light pollution of freshwater ecosystems: principles, ecological impacts and remedies.

Light pollution caused by artificial light at night (ALAN) is increasingly recognized as a major driver of global environmental change. Since emissions are rapidly growing in an urbanizing world and half of the human population lives close to a freshwater shoreline, rivers and lakes are ever more exposed to light pollution worldwide. However, although light conditions are critical to aquatic species, and freshwaters are biodiversity hotspots and vital to human well-being, only a small fraction of studies conducted on ALAN focus on these ecosystems. The effects of light pollution on freshwaters are broad and concern all levels of biodiversity. Experiments have demonstrated diverse behavioural and physiological responses of species, even at low light levels. Prominent examples are skyglow effects on diel vertical migration of zooplankton and the suppression of melatonin production in fish. However, responses vary widely among taxa, suggesting consequences for species distribution patterns, potential to create novel communities across ecosystem boundaries, and cascading effects on ecosystem functioning. Understanding, predicting and alleviating the ecological impacts of light pollution on freshwaters requires a solid consideration of the physical properties of light propagating in water and a multitude of biological responses. This knowledge is urgently needed to develop innovative lighting concepts, mitigation strategies and specifically targeted measures. This article is part of the theme issue 'Light pollution in complex ecological systems'.

Non-contact robotic manipulation of floating objects: exploiting emergent limit cycles.

The study of non-contact manipulation in water, and the ability to robotically control floating objects has gained recent attention due to wide-ranging potential applications, including the analysis of plastic pollution in the oceans and the optimization of procedures in food processing plants. However, modeling floating object movements can be complex, as their trajectories are influenced by various factors such as the object's shape, size, mass, and the magnitude, frequency, and patterns of water waves. This study proposes an experimental investigation into the emergence ofrobotically controlled limit cycles in the movement of floating objects within a closed environment. The objects' movements are driven by robot fins, and the experiment plan set up involves the use of up to four fins and variable motor parameters. By combining energy quantification of the system with an open-loop pattern generation, it is possible to demonstrate all main water-object interactions within the enclosed environment. A study using dynamic time warping around floating patterns gives insights on possible further studies.

Characterization of subsurface pathways contributing to freshwater salinization of urban streams using electrical and electromagnetic imaging techniques.

Salinization of inland fresh surface waters in temperate climates is a growing concern due to increasing salt inputs from sources including chloride (Cl)-containing road salt de-icers, industrial waste, and landfill leachate. Groundwater pathways play an important role in the year-round delivery of Cl to streams, but quantifying this pathway, including spatiotemporal variability and amount of Cl mass stored in the subsurface, is challenging. The objective of this study was to demonstrate, evaluate, and compare the potential applications of the geoelectrical techniques - electromagnetics (EM) and direct current (DC) resistivity - for mapping salt contamination in shallow urban groundwater and characterizing the groundwater pathways delivering Cl to urban streams. EM and DC surveys were conducted (3D mapping and 2D time-lapse) across a 20 m salt-impacted stream section and surrounding riparian zone that is located near an arterial road and parking lot. Groundwater samples and soil cores were also collected to validate the geoelectrical results. Both the EM and DC surveys detected high salt concentrations in the shallow subsurface (up to 3 m depth) near the road, parking lot, and stream; however, DC more accurately represented groundwater Cl concentrations. DC results were used to calculate the total Cl mass in the subsurface, with the spatial mass distribution used to infer the temporal variability in the subsurface salt plume. Finally, time-lapse DC showed that the highest groundwater salt concentrations existed near the stream between June and October - this is expected to contribute to the elevated salt concentrations in the stream during summer months. This study has shown that EM and DC can be useful for identifying groundwater salt concentration, storage, and transport in a non-intrusive and efficient manner, making them valuable field tools for characterizing and quantifying groundwater salt pathways to urban streams.

A review on the hydration properties of dietary fibers derived from food waste and their interactions with other ingredients: opportunities and challenges for their application in the food industry.

Dietary fiber (DF) significantly affects the quality attributes of food matrices. Depending on its chemical composition, molecular structure, and degree of hydration, the behavior of DF may differ. Numerous reports confirm that incorporating DF derived from food waste into food products has significant effects on textural, sensory, rheological, and antimicrobial properties. Additionally, the characteristics of DF, modification techniques (chemical, enzymatic, mechanical, thermal), and processing conditions (temperature, pH, ionic strength), as well as the presence of other components, can profoundly affect the functionalities of DF. This review aims to describe the interactions between DF and water, focusing on the effects of free water, freezing-bound water, and unfreezing-bound water on the hydration capacity of both soluble and insoluble DF. The review also explores how the structural, functional, and environmental properties of DF contribute to its hydration capacity. It becomes evident that the interactions between DF and water, and their effects on the rheological properties of food matrices, are complex and multifaceted subjects, offering both opportunities and challenges for further exploration. Utilizing DF extracted from food waste exhibits promise as a sustainable and viable strategy for the food industry to create nutritious and high-value-added products, while concurrently reducing reliance on primary virgin resources.

Tuning the water interactions of cellulose nanofibril hydrogels using willow bark extract.

Cellulose nanofibrils (CNFs) are increasingly used as precursors for foams, films and composites, where water interactions are of great importance. In this study, we used willow bark extract (WBE), an underrated natural source of bioactive phenolic compounds, as a plant-based modifier for CNF hydrogels, without compromising their mechanical properties. We found that the introduction of WBE into both native, mechanically fibrillated CNFs and TEMPO-oxidized CNFs increased considerably the storage modulus of the hydrogels and reduced their swelling ratio in water up to 5-7 times. A detailed chemical analysis revealed that WBE is composed of several phenolic compounds in addition to potassium salts. Whereas the salt ions reduced the repulsion between fibrils and created denser CNF networks, the phenolic compounds - which adsorbed readily on the cellulose surfaces - played an important role in assisting the flowability of the hydrogels at high shear strains by reducing the flocculation tendency, often observed in pure and salt-containing CNFs, and contributed to the structural integrity of the CNF network in aqueous environment. Surprisingly, the willow bark extract exhibited hemolysis activity, which highlights the importance of more thorough investigations of biocompatibility of natural materials. WBE shows great potential for managing the water interactions of CNF-based products.

Value and limitations of machine learning in high-frequency nutrient data for gap-filling, forecasting, and transport process interpretation.

High-frequency monitoring of water quality in catchments brings along the challenge of post-processing large amounts of data. Moreover, monitoring stations are often remote and technical issues resulting in data gaps are common. Machine learning algorithms can be applied to fill these gaps, and to a certain extent, for predictions and interpretation. The objectives of this study were (1) to evaluate six different machine learning models for gap-filling in a high-frequency nitrate and total phosphorus concentration time series, (2) to showcase the potential added value (and limitations) of machine learning to interpret underlying processes, and (3) to study the limits of machine learning algorithms for predictions outside the training period. We used a 4-year high-frequency dataset from a ditch draining one intensive dairy farm in the east of The Netherlands. Continuous time series of precipitation, evapotranspiration, groundwater levels, discharge, turbidity, and nitrate or total phosphorus were used as predictors for total phosphorus and nitrate concentrations respectively. Our results showed that the random forest algorithm had the best performance to fill in data-gaps, with R2 higher than 0.92 and short computation times. The feature importance helped understanding the changes in transport processes linked to water conservation measures and rain variability. Applying the machine learning model outside the training period resulted in a low performance, largely due to system changes (manure surplus and water conservation) which were not included as predictors. This study offers a valuable and novel example of how to use and interpret machine learning models for post-processing high-frequency water quality data.

Improvement of in situ hydraulic conductivity tests for riverbeds under the influence of infiltration intake: A case study of the Olawa River, SW Poland.

The vertical hydraulic conductivity (KV) of riverbeds plays a pivotal role in controlling water exchanges between the surface water and groundwater. The recent papers have focused on the spatial and temporal variability of KV (K measured in the "z" direction), while point measurement schemes seem to require more attention. The most popular measurement technique is the falling head method, in which a PVC pipe is placed in the riverbed. The most commonly used formulas to determine hydraulic conductivity are those based on the Darcy's law and its modifications. This article presents how the falling head method should be carried out and what errors need to be considered in conditions of the riverbank intake. The simultaneous measurements and analyses of the falling head vs. the riverbed head can produce satisfying outcomes and reduce errors. A full utility MS Excel spreadsheet based on an automated formula with a solver was developed.

Morphology and swelling of thin films of dialcohol xylan.

Polysaccharides are excellent network formers and are often processed into films from water solutions. Despite being hydrophilic polysaccharides, the typical xylans liberated from wood are sparsely soluble in water. We have previously suggested that an additional piece to the solubilization puzzle is modification of the xylan backbone via oxidative cleavage of the saccharide ring. Here, we demonstrate the influence of the degree of modification, i.e., degree of oxidation (DO) on xylan solubilization and consequent film formation and stability. Oxidized and reduced wood xylans (i.e., dialcohol xylans) with the highest DO (77 %) within the series exhibited the smallest hydrodynamic diameter (dh) of 60 nm in dimethylsulfoxide (DMSO). We transferred the modified xylans into films credit to their established solubility and then quantified the film water interactions. Dialcohol xylans with intermediate DOs (42 and 63 %) did not form continuous films. The films swelled slightly when subjected to humidity. However, the film with the highest DO demonstrated a significant moisture uptake that depended on the film mass and was not observed with the other modified grades or with unmodified xylan.

Experimental and simulated distribution and interaction of water in cellulose esters with alkyl chain substitutions.

This study investigated the effect of the average length of substituted side chains in different cellulose esters on water sorption and the water association mechanism. For this purpose, a set of esters with a similar total degree of substitution was selected: cellulose acetate, cellulose acetate propionate, and cellulose acetate butyrate. Dynamic vapor sorption was used to determine the effect of the side chain length on sorption, desorption, and the occurrence of water clustering. Since water association in the structure was of interest, molecular dynamics simulations were performed on cellulose acetate and cellulose acetate propionate. This study showed that cellulose acetate appears to be water-sensitive and experiences hysteresis upon water sorption, which was attributed to structural changes. The simulations also showed that water is screened out by the side chains and forms intermolecular hydrogen bonds, primarily to the carbonyl oxygen rather than the residual hydroxyl groups.

Combined Surface-Subsurface Stream Restoration Structures Can Optimize Hyporheic Attenuation of Stream Water Contaminants.

There is a design-to-function knowledge gap regarding how engineered stream restoration structures can maximize hyporheic contaminant attenuation. Surface and subsurface structures have each been studied in isolation as techniques to restore hyporheic exchange, but surface-subsurface structures have not been investigated or optimized in an integrated manner. Here, we used a numerical model to systematically evaluate key design variables for combined surface (i.e., weir height and length) and subsurface (i.e., upstream and downstream baffle plate spacing) structures. We also compared performance metrics that place differing emphasis on hyporheic flux versus transit times. We found that surface structures tended to create higher flux, shorter transit time flowpaths, whereas subsurface structures promoted moderate-flux, longer transit time flowpaths. Optimal combined surface-subsurface structures could increase fluxes and transit times simultaneously, thus providing conditions for contaminant attenuation that were many times more effective than surface or subsurface structures alone. All performance metrics were improved by the presence of an upstream plate and the absence of a downstream plate. Increasing the weir length tended to improve all metrics, whereas the optimal weir height varied based on metrics. These findings may improve stream restoration by better aligning specific restoration goals with appropriate performance metrics and hyporheic structure designs.

Distinguishing Direct Human-Driven Effects on the Global Terrestrial Water Cycle.

Population growth is increasing the pressure on water resource availability. For useful assessment and planning for societal water availability impacts, it is imperative to disentangle the direct influences of human activities in the landscape from external climate-driven influences on water flows and their variation and change. In this study we used the water balance model, a gridded global hydrological model, to quantify and distinguish human-driven change components, modified by interventions such as dams, reservoirs, and water withdrawals for irrigation, industry, and households, from climate-driven change components on four key water balance variables in the terrestrial hydrological system (evapotranspiration, runoff, soil moisture, storage change). We also analyzed emergent effect patterns in and across different parts of the world, facilitating exploration of spatial variability and regional patterns on multiple spatial scales, from pixel to global, including previously uninvestigated parts of the world. Our results show that human activities drive changes in all hydrological variables, with different magnitudes and directions depending on geographical location. The differences between model scenarios with and without human activities were largest in regions with the highest population densities. In such regions, which also have relatively large numbers of dams for irrigation, water largely tends to be removed from storage and go to feed increased runoff and evapotranspiration fluxes. Our analysis considers a more complete set of hydrological variables than previous studies and can guide further research and management planning for future hydrological and water availability trends, including in relatively data-poor parts of the world.

Assessment of cellulose interactions with water by ssNMR: 1H->13C transfer kinetics revisited.

To evaluate cellulose interactions with water, 1H->13C polarization transfer kinetics during Variable Contact Time CP-MAS NMR spectroscopy were studied and modelled using cellulose of different origins. The increase in the temporal resolution of the plot relating signal intensity to contact-time made it possible to compare different physical models for use in fitting the kinetic curve. These models involve combinations of variables, such as proton spin diffusions, that require a better understanding of their physicochemical and structural bases. To that end, hydrogen interactions were modulated by adding water, first by varying cellulose water content, second by exchanging hydroxyl protons with D2O, and last by varying the spinning rate. The results demonstrate that this approach makes it possible to probe interactions of polysaccharides with structural water, as well as to follow the evolution of the proton-proton interactions during hydration through spin diffusion times.

Understanding groundwater behaviors and exchange dynamics in a linked catchment-floodplain-lake system.

Groundwater and surface water are hydrologically interconnected systems that exhibit dynamic water, heat and mass exchanges. In this study, a conceptual framework was used to investigate groundwater behaviors and associated hydrological exchanges by combining field measurements, digital filtering and analytical approaches, exemplified by a linked catchment-floodplain-lake system (Poyang Lake, China). The results show that the hydrological regime for both groundwater and surface water exhibit a seasonal variability in the lake catchment. Topographically, the lake catchment can be divided into the mountainous baseflow, ungauged lateral groundwater and floodplain groundwater that contribute to the lake storage changes. Although groundwater flow is generally from the mountainous catchment to the lake floodplain areas due to topographic effects, precipitation provides an additional input for the shallow groundwater and is expected to enhance the groundwater dynamics in terms of spatially heterogeneous responses. The estimation indicates that about 40 % of the catchment river discharge may be coming from the mountainous baseflow (~290 × 108 m3/yr) and discharged into the lake through a surface flow pathway. The ungauged groundwater-lake interaction shows the annual discharge volume is up to 10 × 108 m3/yr and associated exchange fluxes tend to be stronger during spring-summer months (23-45 m3/s) than those of autumn-winter months (9-22 m3/s). Additionally, the floodplain groundwater-lake exchange (~9.5 × 108 m3/yr) indicates that groundwater generally receives the lake water during summer months (mean flux = 110 m3/s) and discharges into the lake during other months (90 m3/s) through a subsurface pathway. This study highlights the importance of groundwater's contributions to the surface river-lake system in terms of the flux variability and different transport pathways. The outcomes of this work will benefit future water resources management and applications by providing a methodology for predicting the groundwater hydrology of large lake-catchment systems.

Response of water-biochar interactions to physical and biochemical aging.

Biochar aging may affect the interactions of biochar with water and thus its performance as soil amendment; yet the specific mechanisms underlying these effects are poorly understood. By means of FTIR, N2 adsorption, Hg intrusion porosimetry, thermogravimetric analysis, 13C solid state nuclear magnetic resonance (NMR) and 1H NMR relaxometry, we investigated changes in the chemistry and structure of biochar as well as its interaction with water after biochar aging, both physical (simulated by ball-milling) and biochemical (simulated by co-composting). Three different porosities of biochar were examined: <5 nm, 1 µm and 10 µm diameter sizes. Physical aging caused the disappearance of the porosity at 10 µm. With biochemical aging, biochar underwent an enrichment of oxygenated functional groups either as a result of surface functionalisation processes or by the deposition of fresh organic matter layers on the surface and pores of biochar. 1H NMR relaxometry revealed that the proportion of water strongly interacting with biochar increased with both physical and biochemical aging. Although biochemical aging significantly altered the composition of biochar surface and modulates its interaction with water, 1H NMR relaxometry proved that physical aging had a relatively stronger influence on water mobility and dynamics in biochar, lowering both T1 and T2 relaxation times in the initial contact times of biochar and water.