Understanding the initiation, propagation and evolution of water injection-induced fractures is essential for geo-energy applications. Hydromechanical stimulation experiments were conducted in a deep borehole drilled into crystalline bedrock to gain insights into these processes, involving simultaneous in-situ measurements of three-dimensional fracture displacements, injection flow rates, and water pressure in 2.4-m isolated borehole sections at 500-m depth. Three distinct sections were tested in the COSC-1 borehole (Sweden): a section of intact rock, a section with a hydraulically conductive fracture and a section with nonconductive fractures. Acoustic televiewer measurements conducted before and after the experiments confirmed the generation of new fractures. Accurate positioning of measurement tools was ensured through gamma log profiling and an innovative FFEC-based method for detecting flowing fractures. The tests revealed several transitional pressure values associated with mechanical events, with intact rock requiring the highest pressure to induce fracturing, followed by the nonconductive fracture section and the initially conductive fracture section. Following fluid injection, transient pressure decays were observed that were associated with leakage from newly generated fractures, providing insights into fracture behaviour under stimulation. Vertical displacements were predominant across the different tests, with measured displacements typically ranging from 10 to 100 µm. Fracture activation modes primarily involved the normal opening of subhorizontal fractures that were parallel to the metamorphic foliation, with some irreversible slip at higher pressures. However, a more complex scenario was observed in the test interval with previously nonconductive fractures, involving competition between the opening of subhorizontal fractures and reverse shearing of a steeply dipping fracture.

Background  Antimicrobial resistance is a major global health concern, with the environment playing a key role in its emergence and spread. Understanding the relationships between environmental factors, microbial communities, and resistance mechanisms is vital for elucidating environmental resistome dynamics. In this study, we characterized the environmental resistome of the Baltic Sea and evaluated how environmental gradients and spatial variability, alongside its microbial communities and associated functional genes, influence resistome diversity and composition across geographic regions.

Results  We analyzed the metagenomes of benthic sediments from 59 monitoring stations across a 1,150 km distance of the Baltic Sea, revealing an environmental resistome comprised of predicted antimicrobial resistance genes (ARGs) associated with resistance against 26 antibiotic classes. We observed spatial variation in its resistance profile, with higher resistome diversity in the northern regions and a decline in the dead zones and the southern areas. The combined effects of salinity and temperature gradients, alongside nutrient availability, created a complex environmental landscape that shaped the diversity and distribution of the predicted ARGs. Salinity predominantly influenced microbial communities and predicted ARG composition, leading to clear distinctions between high-saline regions and those with lower to mid-level salinity. Furthermore, our analysis suggests that microbial community composition and mobile genetic elements might be crucial in shaping ARG diversity and composition.

Conclusions  We presented that salinity and temperature were identified as the primary environmental factors influencing resistome diversity and distribution across geographic regions, with nutrient availability further shaping these patterns in the Baltic Sea. Our study also highlighted the interplay between microbial communities, resistance, and associated functional genes in the benthic ecosystem, underscoring the potential role of microbial and mobile genetic element composition in ARG distribution. Understanding how environmental factors and microbial communities modulate environmental resistomes will help predict the impact of future environmental changes on resistance mechanisms in complex aquatic ecosystems.

Background  The Baltic Sea receives petroleum hydrocarbons from various point sources. The degradation of these contaminants in the environment is typically facilitated by a variety of microorganisms that possess a range of genes and metabolic functions related to the degradation of various hydrocarbon substrates. However, our understanding of natural attenuation and the microbial capacity to degrade these contaminants within the Baltic Sea ecosystem remains limited. In this study, we compiled metagenomes from the benthic and pelagic ecosystems across the Baltic Sea to identify microorganisms and characterize their genes and metabolic functions involved in the degradation of hydrocarbon compounds.

Results  Known hydrocarbon-degrading phyla, i.e., Pseudomonadota, Myxococcota A, Actinomycetota, and Desulfobacterota, were identified within the Baltic Sea metagenome-assembled genomes (MAGs). Notably, 80% of the MAGs exhibited multiple hydrocarbon degradation gene annotations (> 10 reads per kilobase million). Aerobic degradation was the predominant pathway for hydrocarbon degradation across environmental samples. Hydrocarbon degradation gene abundances varied among samples and Baltic Sea subbasins, with long-chain alkanes and dibenzothiophene compounds being the preferred substrates. Species richness and diversity of both benthic and pelagic microorganisms positively correlated with hydrocarbon degradation gene diversity, with the pelagic ecosystem exhibiting significantly higher richness and diversity compared to the benthic ecosystem. Additionally, the composition of the hydrocarbon degradation genes across the Baltic Sea subbasins was influenced by oil spill history, with areas that experienced higher spill volumes showing lower microbial diversity, suggesting potential enrichment of specific hydrocarbon degraders. Among the environmental factors assessed, depth played a significant role in shaping the composition of genes involved in hydrocarbon degradation within the Baltic Sea.

Conclusions  Using metagenomics, we profiled the native microorganisms associated with hydrocarbon degradation in the Baltic Sea. This knowledge will aid in understanding the natural capacities of microbial communities, potentially linked to the natural attenuation of hydrocarbon pollutants in the area. Insights into microbial degradation potential can enhance predictions of petroleum pollutant persistence and accumulation, support mitigation strategies for marine pollution, and reveal the ecological resilience of native microbial communities in marine ecosystems.

Stricter national regulations on fertilisation and improved use of manure in agriculture, as well as set-aside schemes, do have measurable effects on the riverine loads of nitrogen to the Baltic Sea, although it can take up to a decade to reach full depletion of legacies in agricultural soils. However, reaching the environmental goals in a cost-effective way, requires more targeted measures to areas where the risk for leaching is high.

Intensive agriculture has been linked to increased nitrogen loads and adverse effects on downstream aquatic ecosystems. Sustained large net nitrogen surpluses have been shown in several contexts to form legacies in soil or waters, which delay the effects of reduction measures. In this study, detailed land use and agricultural statistics were used to reconstruct the annual nitrogen surpluses in three agriculture-dominated watersheds of Denmark (600-2700 km2) with well-drained loamy soils. These surpluses and long-term hydrological records were used as inputs to the process model ELEMeNT to quantify the nitrogen stores and fluxes for 1920-2020. A multi-objective calibration using timeseries of river nitrate loads, as well as other non-conventional data sources, allowed to explore the potential of these different data to constrain the nitrogen cycling model. We found the flux-weighted nitrate concentrations in the root zone percolate below croplands, a dataset not commonly used in calibrating watershed models, to be critical in reducing parameter uncertainty. Groundwater nitrate legacies built up in all three studied watersheds during 1950-1990 corresponding to & SIM;2% of the surplus (or & SIM;1 kg N ha yr-1) before they went down at a similar rate during 1990-2015. Over the same periods active soil nitrogen legacies first accumulated by approximately 10% of the surplus (& SIM;5 kg N ha yr-1), before undergoing a commensurate reduction. Both legacies appear to have been the drivers of hysteresis in the diffuse load at the catchments' outlet and hindrances to reaching water quality goals. Results indicate that the low cropland surpluses enforced during 2008-2015 had a larger impact on the diffuse river loads than the European Union's untargeted grass set-aside policy of 1993-2008. Collectively, the measures of 1990-2015 are estimated to have reset the diffuse load regimes of the watersheds back to the situation prevailing in the 1960s.

The generalized radial flow (GRF) model in well-test analysis employs noninteger flow dimensions to represent the variation in flow area with respect to radial distance from a borehole. However, the flow dimension is influenced not only by changes in flow area, but also by permeability variations in the flow medium. In this report, the flow dimension from the combined effect of flow dimensionality and permeability/conductance variation is interpreted and referred to as apparent flow dimension (AFD). AFD is determined using the second derivative of the drawdown-time plot from pressure transient testing, which may have varied noninteger values with time. A systematic set of investigations is presented, starting from idealized channel networks in one, two and three dimensions (1D, 2D and 3D, respectively), and proceeding to a case study with a complex fracture network based on actual field data. Interestingly, a general relation between the AFD upsurge/dip and the conductance contrast between adjacent flow channels is established. The relation is derived from calculations for 1D networks but is shown to be useful even for data interpretation for more complex 2D and 3D cases. In an application to fracture network data at a real site, the presence of flow channel clusters is identified using the AFD plot. Overall, the AFD analysis is shown to be a useful tool in detecting the conductance/dimensionality changes in the flow system, and may serve as one of the different data types that can be jointly analysed for characterizing a heterogeneous flow system.

Performance assessment of nuclear waste disposal in deep crystalline bedrock demands a thorough understanding of the related flow and transport processes. Uncertainties may arise both from the selection of the conceptual model as well as the estimation of the related model parameters. Discrete fracture network (DFN) models are widely used for such modeling while channel network models (CNM) provide an alternative representation, the latter focusing on the fact that flow and transport in deep fractured media often are dominated by a small number of long preferential flow paths. This study applies the principle of channel networks, implemented in the Pychan3d simulator, to analyze the hydraulic and tracer transport behavior in a 450-m-deep fractured granite system at the aspo Hard Rock Laboratory in Sweden, where extensive site characterization data, including hydraulic and tracer test data are available. Semi-automated calibration of channel conductances to field characterization data (flow rates, drawdowns, and tracer recoveries) is performed using PEST algorithm. It was observed that an optimal CNM connectivity map for channel conductance calibration can only be developed by jointly fitting flow rates, drawdowns and tracer mass recovery values. Results from data-calibrated CNM when compared to a corresponding calibrated DFN model shows that the CNM calibrates and adapts better than a DFN model with uniform fracture surfaces. This comparative study shows the differences and uncertainties between two models as well as examines the implications of using them for long term model predictions.

Water bodies provide essential ecosystem services linked to morphometric features that might differ between natural lakes and reservoirs. We use the HydroLAKES global dataset to quantitatively compare large (area > 1 km²) reservoirs and natural lakes in terms of scaling exponents between morphometric measures (volume, area, shore length). These exponents are further compared to those expected from geometrical assumptions and constraints. Lakes cover a larger range of volumes for the same range of surface areas than reservoirs, and have a larger volume-area scaling exponent. The volume-area scaling exponent for reservoirs (but not natural lakes) and the area-shore length exponent for all water bodies follow the predictions for self-affine surfaces. Land cover and terrain influence the scaling relations more for lakes than for reservoirs. These morphometric differences may be used to model the impact of reservoirs and lakes on hydrological processes and associated ecosystem services at regional to global scales.

To characterize the coupled hydromechanical behavior of rock fractures, the step-rate injection method for fracture in-situ properties (SIMFIP) was conducted with a specialized downhole probe developed by Guglielmi et al. (2014, https://doi.org/10.1007/s00603-013-0517-1). In June 2019, a field campaign was carried out near Åre, Sweden, where the SIMFIP probe was applied in the Collisional Orogeny in the Scandinavian Caledonides-1 scientific borehole to understand the dynamics of injection-induced fracture initiation, fracture opening, and shearing due to water injection-withdrawal in a borehole interval isolated by two packers. Three intervals were investigated at ∼500 m depth: (a) an unfractured section (intact rock), (b) a section with non-conductive fractures, and (c) a section with hydraulically conductive fractures. Pressure, injection flow rate, and borehole wall displacement were simultaneously measured during the tests. In the present study, the geometry of the induced fracture and deformation of existing fractures at different time stages of the tests are determined based on a hydrologic model by using the measured pressure and flow data during each time stage of the experiment. A numerical model for the fluid flow within the fracture and the packed-off borehole interval is implemented within COMSOL Multiphysics. By matching model simulations with observed data for all three sections, estimates of the induced and propagated fractures' radius and aperture at successive time stages have been obtained in each case. We could also determine the non-linear relationship between fracture aperture and pressure for values above fracture opening pressures. The model results provide insights for the understanding of pressure-induced fracture initiation and propagation in crystalline rock.

Quick-clay landslides are common geohazards in Nordic countries and Canada. The presence of potential quick clays is confirmed using geotechnical investigations, but near-surface geophysical methods, such as seismic and resistivity surveys, can also help identify coarse-grained materials associated with the development of quick clays. We present the results of reflection seismic investigations on land and in part of the Gota River in Sweden, along which many quick-clay landslide scars exist. This is the first time that such a large-scale reflection seismic investigation has been carried out to study the subsurface structures associated with quick-clay landslides. The results also show a reasonable correlation with radio magnetotelluric and travel-time tomography models of the subsurface. Other ground geophysical data, such as high magnetic values, suggest a positive correlation with an increased thickness of the coarse-grained layer and shallower depths to the top of the bedrock and the top of the coarse-grained layer. The morphology of the river bottom and riverbanks, e.g. subaquatic landslide deposits, is shown by side-scan sonar and bathymetric data. Undulating bedrock, covered by subhorizontal sedimentary glacial and postglacial deposits, is clearly revealed. An extensive coarse-grained layer (P-wave velocity mostly between 1500 and 2500 m s(-1) and resistivity from approximately 80 to 100 Omega m) exists within the sediments and is interpreted and modelled in a regional context. Several fracture zones are identified within the bedrock. Hydrological modelling of the coarse-grained layer confirms its potential for transporting fresh water infiltrated in fractures and nearby outcrops located in the central part of the study area. The modelled groundwater flow in this layer promotes the leaching of marine salts from the overlying clays by seasonal inflow-outflow cycles and/or diffusion, which contributes to the formation of potential quick clays.