Flow Path Length Research Articles

Identifying Subsurface Connectivity From Observations: Experimentation With Equifinality Defines Both Challenges and Pathways to Progress

ABSTRACTLinkages between landscapes and streams are increasingly described in terms of hydrological connectivity. The ability to effectively distinguish different patterns of water movement through catchments makes connectivity particularly interesting to both scientists and practical water managers. Hydrometric data (groundwater levels, soil moisture and streamflow) are often employed to infer the connection between the landscape and its drainage network. Such observational data, however, are insufficient to infer subsurface connectivity in humid settings with perennial stream flow, due to the risk of equifinality. To quantify how much subsurface flow patterns can differ and still be consistent (equifinal) with comprehensive observations of hillslope groundwater levels and stream runoff (the hydrometric data), this study used a modelling experiment based on a well‐characterised field site. Particle‐tracking simulations at different flow rates defined the water flow paths and transit times of two virtual hillslopes that differed profoundly in the vertical distribution of the saturated hydraulic conductivity. Even though the simulated weekly stream flows and groundwater levels were similar (i.e., the hillslopes were hydrometrically equifinal) particle velocities and water ages at specific locations along these hillslopes differed by orders of magnitude. Flow path lengths and catchment transit times varied up to several 100%. The hillslope‐ and stream‐based metrics used to describe connectivity also varied with stream flow rates. These results underline the need to recognise the risks for equifinality when inferring subsurface connectivity from hydrometric observations alone, even when those observations are comprehensive. The results also highlight the value of model simulations for quantifying the uncertainty in the inferred connectivity, targeting the best sampling locations/times to reduce this uncertainty with tracer data and better understanding the way connectivity influences stream chemistry.

Realization of Integrated Regional Ecological Management Based on Ecosystem Service Supply and Demand Flow Networks: An Example from a Dominant Mineral Resources Development Area

Understanding the flow processes and pattern optimization of ecosystem services (ESs) supply and demand is crucial for integrated regional ecological management. However, the understanding of the flow process of ESs at the 1 km grid scale is still limited, especially in areas dominated by mineral resource development. The landscape in these areas has undergone significant changes due to mining activities. It is urgent to construct a regional management model that integrates the flow of ecosystem services and mine restoration. This study developed a framework that links ecosystem service flows (ESFs) and ecological security patterns (ESP) based on multi-source ecological monitoring data, constructed an ES supply-demand flow network through the flow properties, and determined the sequence and optimization strategies for mine rehabilitation to achieve integrated regional management. The results show that, except for food production (FP), other services were in surplus overall, mostly in synergistic relationships, but the spatial distribution of their supply and demand was not coordinated. Surplus areas were located mainly in the eastern woodlands, and deficit areas were located in the northwestern production agglomeration centers, suggesting that areas of supply-demand imbalance can be mitigated through ecological integration. Among these, water yield (WY) had a small number of sources and sinks and is limited in area range. Habitat quality (HQ) sources and sinks had the largest area coverage and the highest number. The distribution of ESF corridors, influenced by factors such as the number of sources and sinks, flow characteristics, and spatial resistance, varied significantly. HQ exhibited a more uniform distribution range, while WY had a longer average length of flow path. Overlaying ecological and mining factors, we identified ecological strategic spots, important supply areas, beneficiary areas, and mine priority restoration areas to further optimize the overall layout and rationally allocate the intrinsic structure of the patches based on ES supply and demand.

Antecedent Hydrologic Conditions Reflected in Stream Lithium Isotope Ratios During Storms

AbstractAntecedent hydrological conditions are recorded through the evolution of dissolved lithium isotope signatures (Li) by juxtaposing two storm events in an upland watershed subject to a Mediterranean climate. Discharge and Li are negatively correlated in both events, but mean Li ratios and associated ranges of variation are distinct between them. We apply a previously developed reactive transport model (RTM) for the site to these event‐scale flow perturbations, but observed shifts in stream Li are not reproduced. To reconcile the stability of the subsurface solute weathering profile with our observations of dynamic stream Li signatures, we couple the RTM to a distribution of fluid transit times that evolve based on storm hydrographs. The approach guides appropriate flux‐weighting of fluid from the RTM over a range of flow path lengths, or equivalently fluid residence times. This flux‐weighted RTM approach accurately reproduces dynamic storm Li‐discharge patterns distinguished by the antecedent conditions of the watershed.

Modelling transport pathways of faults with low hydraulic connectivity in mudstones with low swelling capacity

Faults in some deep mudstones have poor hydraulic connectivity owing to high normal stress on the fault planes. Designing a method for modelling solute transport pathways in such faults/fractures using available data is a critical issue vis-à-vis the safety assessment of radioactive waste disposal. In this study, faults in deep siliceous mudstones with low swelling capacity are investigated using cross-hole hydraulic and tracer tests between two boreholes. The water pressures observed during the hydraulic test are reproduced via hydraulic simulations, assuming a one-dimensional (1D) flow channel with a length 8–80 times the linear distance (4.5 m) between the two borehole test sections. This flow channel length indicates a tortuosity remarkably higher than that previously reported in discrete fracture network simulations and laboratory experiments using a single fracture (i.e. 1–5). Advection–dispersion simulations using a pipe model with a tortuosity of 8–80 reproduces a time series of tracer concentrations observed during the tracer test. The radius of the pipe simulated in the tracer test was 7–39 times the flow channel radius estimated via the hydraulic test, which is in agreement with the ratio of 5–20 previously reported based on in situ tracer and hydraulic tests of faulted/fractured rocks. The simulated longitudinal dispersivity is 1/9–1/20 of the estimated flow path lengths, in agreement with the laboratory and field tracer study values. These results indicate that transport pathways in faults with low hydraulic connectivity can be modelled using a highly tortuous 1D pipe flow path.

Physical design for microfluidic biochips considering actual volume management and channel storage

In recent years, microfluidic biochips have been widely applied in various fields of human society. The optimization design of system-architecture based on continuous-flow microfluidic biochips has been widely studied. However, most previous work was based on the traditional chip architecture with dedicated storage, which not only limits the performance of biochips but also increases their manufacturing costs. In order to improve the execution efficiency and reduce the manufacturing cost, a distributed channel-storage architecture can be used to temporarily cache intermediate fluids in idle flow channels. Under this architecture, careful consideration of the volume management of the fluid to be cached is a prerequisite for ensuring the reliability of bioassay results. However, the existing work has not considered the volume management of the fluid to be cached in detail. This may cause the volume of the fluid to not match the capacity of the storage channel, which can contaminate other fluids and lead to incorrect bioassay results or increase the manufacturing cost of biochips due to long storage channels. In this paper, we propose a physical design method for microfluidic biochips that considers the actual volume of fluid while utilizing distributed channel storage. We address this problem by taking a placement and routing co-design strategy throughout the iterative process of the simulated annealing algorithm. Experimental results under multiple benchmarks show that the proposed method can effectively reduce the completion time of bioassays, minimize the flow path length, and decrease the number of intersections.

Incorporating flowpaths as an explicit measure of river-floodplain connectivity to improve estimates of floodplain sediment deposition

Variation in floodplain topography can lead to gradual flooding and increase river-floodplain connectivity. We show that incorporating flowpaths as an explicit measure of river-floodplain connectivity can improve estimates of floodplain sediment deposition. We focus on the floodplain of the South River, downstream of Waynesboro, Virginia, where measurements of mercury accumulation have been used to estimate decadal-scale sedimentation rates. We developed a two-dimensional Hydrologic Engineering Center's River Analysis System (2D HEC-RAS) hydrodynamic model and used simulated model results with sediment deposition data to create regression models describing sedimentation across the floodplain. All of our statistical models incorporated a flowpath length from the location on the floodplain downstream to the riverbank as an explicit measure of river-floodplain connectivity that improved our estimates of floodplain sediment deposition (r2 = 0.514). We applied our best regression model to our hydrodynamic model results to create a map of floodplain sedimentation rate and discuss differences of three separate sections of floodplain. We found that floodplains with variable topography had wider, bimodal probability distribution functions (PDFs) of sedimentation rate (aggregated spatially) than floodplains without this topographic relief (with narrower log-normal PDFs). Our work highlights how floodplain topography and river-floodplain connectivity affect sedimentation rates and can help inform the development of floodplain sediment budgets.

Design of a 130 MW Axial Turbine Operating with a Supercritical Carbon Dioxide Mixture for the SCARABEUS Project

Supercritical carbon dioxide (sCO2) can be mixed with dopants such as titanium tetrachloride (TiCl4), hexafluoro-benzene (C6F6), and sulphur dioxide (SO2) to raise the critical temperature of the working fluid, allowing it to condense at ambient temperatures in dry solar field locations. The resulting transcritical power cycles have lower compression work and higher thermal efficiency. This paper presents the aerodynamic flow path design of a utility-scale axial turbine operating with an 80–20% molar mix of CO2 and SO2. The preliminary design is obtained using a mean line turbine design method based on the Aungier loss model, which considers both mechanical and rotor dynamic criteria. Furthermore, steady-state 3D computational fluid dynamic (CFD) simulations are set up using the k-ω SST turbulence model, and blade shape optimisation is carried out to improve the preliminary design while maintaining acceptable stress levels. It was found that increasing the number of stages from 4 to 14 increased the total-to-total efficiency by 6.3% due to the higher blade aspect ratio, which reduced the influence of secondary flow losses, as well as the smaller tip diameter, which minimised the tip clearance losses. The final turbine design had a total-to-total efficiency of 92.9%, as predicted by the CFD results, with a maximum stress of less than 260 MPa and a mass flow rate within 1% of the intended cycle’s mass flow rate. Optimum aerodynamic performance was achieved with a 14-stage design where the hub radius and the flow path length are 310 mm and 1800 mm, respectively. Off-design analysis showed that the turbine could operate down to 88% of the design reduced mass flow rate with a total-to-total efficiency of 80%.

Simulation of Water Flow Path Length (WFPL) and Water Film Depth (WFD) for Wide Expressway Asphalt Pavement

This paper simulates actual rainfall conditions and raindrops flowing to form a water flow path (WFP) on the pavement surface of the wide expressway. Then, the different linear combination conditions, including longitudinal slope (LS), transverse slope (superelevation, TS), gradual change rate of TS, and pavement width (PW), were simulated and analyzed. The results show that (1) the influence of each linear index on the maximum water film path length (WFPLmax) and maximum water flow depth (WFDmax) differs (according to the absolute values of Beta, LS has the greatest influence on WFPLmax, and PW has the greatest influence on the WFDmax for both straight-line and circular-curve sections); (2) when the design value of LS is between 1.1% and 4%, the WFDmax can be effectively reduced by lowering the value of LS; (3) in the case of a high design value of LS, it can be considered to increase the TS of the pavement arch from 2% to 2.5% to effectively reduce the WFPLmax, and the wider PW, the better the reducing effect; (4) while widening the expressway, adjusting the TS from 2% to 2.5% can effectively offset the increasing effect of PW on the WFDmax. This research aims to fill the research gap in the simulation of runoff characteristics of wide expressway asphalt pavements and to improve the alignment design of expressways from the drainage perspective for the improvement of driving safety.

Landscape structures regulate the contrasting response of recession along rainfall amounts

Abstract. Streamflow recession, shaped by hydrological processes, runoff dynamics, and catchment storage, is heavily influenced by landscape structure and rainstorm characteristics. However, our understanding of how recession relates to landscape structure and rainstorm characteristics remains inconsistent, with limited research examining their combined impact. This study examines this interplay in shaping recession responses upon 291 sets of recession parameters obtained through the decorrelation process. The data originate from 19 subtropical mountainous rivers and cover events with a wide spectrum of rainfall amounts. Key findings indicate that the recession coefficient (a) increases while the exponent (b) decreases with the L/G ratio (the median of ratios between flow-path length and gradient), suggesting that longer and gentler hillslopes facilitate flow accumulation and aquifer connectivity, ultimately reducing nonlinearity. Additionally, in large catchments, the exponent (b) increases with increasing rainfall due to greater landscape heterogeneity. Conversely, in small catchments, it declines with rainfall, indicating that these catchments have less landscape heterogeneity and thus reduced runoff heterogeneity. Our findings underscore the necessity for further validation of how L/G and drainage area regulate recession responses to varying rainfall levels across diverse regions.

Numerical Modelling Study of Subsurface Drainage of Permeable Friction Course Considering Road Geometric Designs

This study aimed to evaluate the subsurface drainage of a permeable friction course (PFC) via two-dimensional finite element analysis. To achieve the scope, PFCs with equivalent water flow paths of length values of 10, 15, 20, and 30 m and slope values of 0.5%, 2%, 4%, 6%, and 8% were modelled based on FEniCS and implemented entirely in Python programing language to extract the time for surface ponding according to a range of rainfall intensities. The results show that when the rainfall intensity and the length of equivalent water flow path of the PFC rose, the time for surface ponding decreased. For instance, with a rainfall intensity of 10 mm/h and a slope of 0.5%, when the length of equivalent water flow path increased by 20 m, the time for surface ponding dropped by 21 min. Moreover, when the slope of the equivalent water flow path and the thickness of the PFC increased, the time for surface ponding increased. For instance, with a rainfall intensity of 10 mm/h, and a PFC with an equivalent length of 10 m, when the slope increased by 16 times, the time for surface ponding increased more than two times. The current study highlights that the thickness of the PFC has the most influence on subsurface drainage. The findings of this study indicate that at high rainfall intensities, the subsurface drainage of a PFC is not sensitive to its geometric design. Further experimental investigations are needed to evaluate and validate the subsurface drainage of a PFC considering permeability, rutting, and environmental factors.

Effects of a large-scale dam structure on upstream and downstream lateral hyporheic exchange and residence time distributions – The Xinglong Water Conservancy Dam, China

Anthropogenic in-stream structures such as large dams can alter stream hydrodynamics and cause river stage fluctuations that have the potential to affect hyporheic exchange flow (HEF) patterns between rivers and groundwater. However, the development of HEF patterns and residence time distribution (RTD) in dam-regulated riverbank aquifers has mostly been studied focusing on areas downstream of large dam structures while upstream areas have often been neglected. In order to investigate the HEF and RTD patterns under the influence of a large dam in upstream and downstream areas simultaneously, groundwater level data was collected from thirteen monitoring wells distributed over three monitoring profiles both upstream and downstream of XingLong water conservancy dam (XLD), in the Han river, China. Based on field observations, a two-dimensional, horizontal numerical model was built to assess the spatiotemporal evolution of lateral hyporheic exchange flow, the spatial extent of the hyporheic zone (HZ) and groundwater RTD under the influence of the XLD. Our results demonstrate that the impact of the XLD is different in upstream and downstream aquifer regions. The storage effect induced by the XLD caused significant groundwater recharge upstream of the dam over long period of model, while significant groundwater discharge occurs downstream of the dam. As the stage of the regulated Han River varies, the HZ along the upstream river boundary significantly enhanced by the dam whereas substantial HEF downstream was observed in response to rising or high river stage. Furthermore, the operation of the XLD resulted in some substantial flow around the dam characterized by short flow path lengths and high flow velocity, which consequently resulted in the development of a significant HZ in the perimeter surrounding the XLD. Residence time distributions showed that the water in the HZ downstream of the XLD is dominated by rejuvenated (younger) groundwater, whereas the HZ upstream of the XLD comprises a mix of older and younger waters. HEF in the immediate vicinity of the dam comprises both aged and rejuvenated discharged groundwater. Our findings emphasize the need to more holistically consider river stage dynamics in river corridors where large dams are operated to better assess and understand the consequences on riverbank storage as well as flow, transport and attenuation patterns in connected floodplain aquifers and their impact on the aquatic environment.

Flow paths and wetness conditions explain spatiotemporal variation of nitrogen retention for a temperate, humid catchment

Excess use of nitrogen (N) fertilizer has threatened the quality of drinking water and destroyed the structure and functions of aquatic ecosystems. The N retention is defined as a catchment’s capacity to withhold nitrogen from entering surface water bodies during a certain period. Previous studies have focused on N retention in response to catchment characteristics and different climates. However, little has been done to explore the spatiotemporal variation of N retention within a catchment. This study is aimed to explore the spatiotemporal variation of the N retention and the effect of the external N input heterogeneity on the retention capacity and the N export in a small temperate, humid catchment in central Germany. A spatially distributed model was built using the fully coupled surface–subsurface model HydroGeoSphere to simulate the flow and transport of nitrate. Three scenarios considering the external N input heterogeneity, which represent different land-use patterns, were created: (i) the input was uniformly distributed over the catchment, (ii) high input was associated with high retention regions (e.g. near the hillslope top), and (iii) high input was associated with low retention regions (e.g. near the stream). The results showed that N retention exhibits strong spatial heterogeneity, which is controlled by the lengths of flow paths connecting the regions and surface water body and regional wetness conditions within a catchment. N retention showed pronounced seasonality in temperate and humid climates, being higher in wetting, drying and dry periods, but lower in wet periods. The pattern is driven by distinct shifts in the dominance of deeper flow paths and fast shallow flow paths determined by the climate seasonality. Highest denitrification, lowest export fluxes and in-stream nitrate concentration were found in scenario (ii) among the scenarios. Therefore, it was suggested that agricultural activities may avoid regions near surface water bodies to improve the surface water quality. The results of this study provide insights into the nitrogen dynamics and can guide the management of surface water quality.

Attenuation of trace organic compounds along hyporheic flow paths in a lowland sandbed stream

Knowledge of flow paths in hyporheic zones of streams is essential for understanding different turnover processes of various chemicals, such as trace organic compounds (TrOCs). Due to the spatial heterogeneity of hyporheic zones, exact flow paths in situ are unknown and thus, pore water sampling along individual flow paths remains a major challenge. To this end, we developed and tested a novel method to enforce specific flow paths under field conditions. Therefore, U-shaped pipes of 40 cm total flow path length were installed in flow direction along a river cross section in the hyporheic zone of the urban lowland River Erpe (Germany). This allowed us to determine sediment properties and to collect pore water samples along the specific flow paths during 16 consecutive hours. The pore water was analyzed for electrical conductivity, redox parameters, and TrOCs. Diurnal electrical conductivity fluctuations in the surface water and their propagation to the end of the flow paths resulted in median residence times of 11.67 or 16.92 h for different specific flow paths. A clear redox zonation along the flow paths was found. Of the 18 TrOCs studied, five were strongly attenuated (i.e., relative attenuation to surface water > 75 %), three exhibited moderate attenuation (75 to 25 %), nine were only slightly or not attenuated (< 25 %), and one compound increased in concentration (relative formation > 25 %) along the entire flow paths. For TrOCs that showed strong or moderate attenuation, the attenuation occurred mainly within the first 5 to 15 cm. This work highlights (1) the importance of knowing the specific flow paths for a comprehensive understanding of TrOC turnover in the hyporheic zone and (2) the relevance of appropriate biogeochemical and hydrological conditions for the attenuation of TrOCs.

Shallow and local or deep and regional? Inferring source groundwater characteristics across mainstem riverbank discharge faces

AbstractRiverbank groundwater discharge faces are spatially extensive areas of preferential seepage that are exposed to air at low river flow. Some conceptual hydrologic models indicate discharge faces represent the spatial convergence of highly variable age and length groundwater flowpaths, while others indicate greater consistency in source groundwater characteristics. Our detailed field investigation of preferential discharge points nested across mainstem riverbank discharge faces was accomplished by: (1) leveraging new temperature‐based recursive estimation (extended Kalman Filter) modelling methodology to evaluate seasonal, diurnal, and event‐driven groundwater flux patterns, (2) developing a multi‐parameter toolkit based on readily measured attributes to classify the general source groundwater flowpath depth and flowpath length scale, and, (3) assessing whether preferential flow points across discharge faces tend to represent common or convergent groundwater sources. Five major groundwater discharge faces were mapped along the Farmington River, CT, United States using thermal infrared imagery. We then installed vertical temperature profilers directly into 39 preferential discharge points for 4.5 months to track vertical discharge flux patterns. Monthly water chemistry was also collected at the discharge points along with one spatial synoptic of stable isotopes of water and dissolved radon gas. We found pervasive evidence of shallow groundwater sources at the upstream discharge faces along a wide valley section with deep bedrock, as primarily evidenced by pronounced diurnal discharge flux patterns. Discharge flux seasonal trends and bank storage transitions during large river flow events provided further indication of shallow, local sources. In contrast, downstream discharge faces associated with near surface cross cutting bedrock exhibited deep and regional source flowpath characteristics such as more stable discharge patterns and temperatures. However, many neighbouring points across discharge faces had similar discharge flux patterns that differed in chloride and radon concentrations, indicating the additional effects of localized flowpath heterogeneity overprinting on larger scale flowpath characteristics.

Effect of mould surface roughness adjustments to increase the flow path length in the injection moulding process

To increase the flow path length in the mould for thin-walled packaging applications or thick-walled technical components with long flow paths, a mould surface machined by electrical discharge machining is used. By varying the process parameters, the moulding compound and the roughness of the surface, the influence of the mould surface structure on the achievable flow path length is investigated. A flow meander mould is used for this purpose. A material dependent, positive influence of the surface roughness on the flow path length can be determined. However, the influence is small compared to other influences such as mould or melt temperature.

A Multimodel Framework for Quantifying Flow and Advective Transport Controlled by Earthquake-Induced Canister Failures in a Reference Case for Radioactive Waste Geological Disposal

Characterizing flow and transport for earthquake-induced shear canister failure is critical for the performance and safety assessment of radioactive waste geological disposal. The study presents a modeling framework that integrates multiple models to account for fractures produced by shear displacements, evaluate canister failures, and simulate flow and advective transport in a conceptual repository site based on a selected reference case in an offshore island in western Taiwan. The typical KBS-3 disposal concept associated with 500 realizations of the shear-induced fracture properties is employed to quantify the uncertainty of flow and advective transport in the geological disposal site. The radionuclides in canisters are assumed to migrate through the shear-induced fractures surrounding the deposition holes. The results from 500 realizations show that two types of fractures produce a high potential to destroy canisters induced by the shear displacements. The earliest canister failure time influenced by possible shear movements is 0.23 million years for the reference case. The modeling framework identifies five canisters and the associated shear-induced fractures for flow and advective transport simulations. Based on the results of the density-dependent flow fields, the particle tracking algorithm enables the calculations of flow and transport parameters, including equivalent initial flux, equivalent flow rate, path length, travel time, and flow-related transport resistance for the identified five canisters. These parameters are critical for the performance and safety assessments of buffer erosion and canister corrosion near the disposal repository and the far field of the radioactive waste disposal site.

Influencing factors and characterization methods of the anti‐clogging ability of labyrinth channel emitters

AbstractTo quantitatively evaluate the anti‐clogging performance of labyrinth channel emitters, literary data were used to study the clogging of different emitters under different conditions. The results showed that the emitter clogging comprehensive evaluation index (Ia) can be used as an inherent property of labyrinth flow channel emitters, denoted as the anti‐clogging ability characterization index, which is affected by flow channel structure parameters and dimensions. The larger Ia is, the more likely it is for the emitter to be clogged and for clogging to be more serious. The emitter anti‐clogging ability is affected by the emitter discharge (Q), flow path length (L), width (W), depth (D), cross‐sectional area (A, which equals W × D), minimum sectional size (min(D, W)), etc. The Ia value with Q, W, D, A and min(D, W) increases first decreases and then increases, and with L increases first increases and then decreasing, and a critical threshold range exists. Ia can be estimated by Q, W, D, A and min(D, W) and can be used as an effective indicator to quickly estimate emitter anti‐clogging ability and emitter configuration in drip irrigation systems. In practical applications, an emitter with x &lt; 0.47 (flow index) should be selected, and an emitter with Ia &lt; 0.45 has a better hydraulic performance and anti‐clogging ability. To obtain a better anti‐clogging ability, Q, W, D, A and min(D, W) should be 1.33–2.13 L/h, 0.71–0.90 mm, 0.55–0.71 mm, 0.33–0.68 mm2 and 0.60–0.75 mm, respectively. When L is 26.47–58.31 mm, the emitter anti‐clogging ability is relatively low.

Thermal performance of a hybrid thermal management system based on the half helical coil coupled with air jet cooling for cylindrical Lithium-ion battery

The Battery Thermal Management System (BTMS) is crucial for the efficient and safe operation of Lithium-ion batteries (LIBs). It is challenging for liquid cooling to apply to areas of battery modules with complex surface forms, such as the positive and negative cell terminals, cell holders, etc. Herein, the thermal management scheme integrated with the half helical coil and air jet impingement is proposed for cylindrical LIBs. Heat is transferred from the LIBs' sides by helical flow and dissipated from the LIBs' ends by air jet impingement. The three-dimensional computational fluid dynamics model of battery modules with cell holders is developed. The thermal performance of the LIBs at various BTMS strategies, the number of turns of the helical coil, airflow inlet velocity, inlet temperature, and contact resistance are investigated. When the impinging air jets are arranged on both ends of the cylindrical form, the maximum values of temperature difference (ΔT) of the LIBs fall from 5.7°C to 4.0°C compared with that of liquid-based BTMS. The arrangement of the half helical duct divided into two separate sections reduces the coolant flow path length, which decreases the maximum temperature (Tmax) and ΔT to 28.9°C and 3.6°C. The proposed BTMS exhibits an enhanced heat dissipation performance and provides a design guideline for developing the hybrid BTMS.

Assessment of streamwater age using water stable isotopes in a headwater catchment of the central Tibetan Plateau

Global warming has significantly impacted the hydrological processes in alpine cryosphere region. Water age is an essential descriptor of the hydrological function within a catchment. However, the mechanism of streamwater age variability remains unclear due to limited observational data and high altitudes of alpine catchment. In this study, long-term stable isotopic data on streamwater in a catchment in the central Tibetan Plateau (TP) were collected to assess the water age using the sine-wave approach and gamma distribution. Results showed that the mean streamwater age was 77 days, and that 30 % of streamwater was less than 41 days old on average. The streamwater age in this study was relatively younger than that in low-elevation natural catchments, indicating that the rapid drainage process occurs within the glacier and permafrost catchment. The fraction of young water (Fyw) of the streamwater decreased from 39 % at an upstream site to 28 % at the outlet, revealing the impact of permafrost (low Fyw: 25 %) on streamwater age. These variabilities were related to glacier and permafrost coverage, specifically in catchments with higher glacier coverage that are prone to have a lower water age. Temporally, the streamwater age was significantly influenced by precipitation, relative humidity, and glacier change and, to a lesser extent, permafrost change. Mechanically, glacier and permafrost changes influenced the water age by increasing the vertical flowpath length. This study provides new insights into the change in hydrological processes in alpine headwater catchments under global warming.

Insights into the streamwater age in the headwater catchments covered by glaciers and permafrost, Central Tibetan Plateau

Improving our understanding of streamwater age knowledge is critical for revealing the complex hydrological processes in alpine cryosphere catchments. However, few studies on water age have been conducted in alpine cryosphere catchments due to the complicated and inclement environment. In this study, the Buqu catchment, a typical alpine catchment covered by glaciers and permafrost on the central Tibetan Plateau (TP), was selected as the study area. Using the sine-wave approach and a gamma model based on the seasonal cycle of stable isotopes in water, the young water fraction (Fyw) and mean transit time (MTT) of the Buqu catchment outlet and 23 sub-catchments was estimated to comprehensively reveal the potential driving mechanism of water age variability. The streamwater MTT for the entire catchment was 107 days, and 15.1 % of the streamwater was younger than 41 days on average. The estimated water age showed significant spatial heterogeneity with shorter water ages in high-elevation and glacier catchments and longer water ages in low-elevation and non-glacier catchments. Precipitation was the primary driver for spatial variations in water age, while the thickness of the permafrost active layer may function as an intermediate hub to drive water age variability. Mechanically, the thickness of the permafrost active layer controls the water ages by modifying the flow direction and length of water flow path. Spatially, this control mechanism is indirectly driven by the elevation gradient. The TDS concentration in streamwater is significantly related to water age, thus revealing a close link between water quality and hydrology. Our findings suggest that cryosphere retreats likely alter water age, thereby slowing water circulation rates and affecting water quality security under global warming. This study provides insights into the evolution of water ages, thereby deepening our understanding of the hydrological processes and guiding the protection of water resources in alpine headwater catchments.