A semi-analytical model for coupled flow in stress-sensitive multi-scale shale reservoirs with fractal characteristics

A large number of nanopores and complex fracture structures in shale reservoirs results in multi-scale flow of oil. With the development of shale oil reservoirs, the permeability of multi-scale media undergoes changes due to stress sensitivity, which plays a crucial role in controlling pressure propagation and oil flow. This paper proposes a multi-scale coupled flow mathematical model of matrix nanopores, induced fractures, and hydraulic fractures. In this model, the micro-scale effects of shale oil flow in fractal nanopores, fractal induced fracture network, and stress sensitivity of multi-scale media are considered. We solved the model iteratively using Pedrosa transform, semi-analytic Segmented Bessel function, Laplace transform. The results of this model exhibit good agreement with the numerical solution and field production data, confirming the high accuracy of the model. As well, the influence of stress sensitivity on permeability, pressure and production is analyzed. It is shown that the permeability and production decrease significantly when induced fractures are weakly supported. Closed induced fractures can inhibit interporosity flow in the stimulated reservoir volume (SRV). It has been shown in sensitivity analysis that hydraulic fractures are beneficial to early production, and induced fractures in SRV are beneficial to middle production. The model can characterize multi-scale flow characteristics of shale oil, providing theoretical guidance for rapid productivity evaluation.

Coupled multiphysical model for investigation of influence factors in the application of microbially induced calcite precipitation

The study presents a comprehensive coupled thermo-bio-chemo-hydraulic(T-BCH)modeling framework for stabilizing soils using microbially induced calcite precipitation(MICP).The numerical model considers relevant multiphysics involved in MICP,such as bacterial ureolytic activities,biochemical reactions,multiphase and multicomponent transport,and alteration of the porosity and permeability.The model incorporates multiphysical coupling effects through well-established constitutive relations that connect parameters and variables from different physical fields.It was implemented in the open-source finite element code OpenGeoSys(OGS),and a semi-staggered solution strategy was designed to solve the couplings,allowing for flexible model settings.Therefore,the developed model can be easily adapted to simulate MICP applications in different scenarios.The numerical model was employed to analyze the effect of various factors,including temperature,injection strategies,and application scales.Besides,a TBCH modeling study was conducted on the laboratory-scale domain to analyze the effects of temperature on urease activity and precipitated calcium carbonate.To understand the scale dependency of MICP treatment,a large-scale heterogeneous domain was subjected to variable biochemical injection strategies.The simulations conducted at the field-scale guided the selection of an injection strategy to achieve the desired type and amount of precipitation.Additionally,the study emphasized the potential of numerical models as reliable tools for optimizing future developments in field-scale MICP treatment.The present study demonstrates the potential of this numerical framework for designing and optimizing the MICP applications in laboratory-,prototype-,and field-scale scenarios.

Effects of asymmetric coupling and boundary on the dynamic behaviors of a random nearest neighbor coupled system

This study investigates the dynamical behaviors of nearest neighbor asymmetric coupled systems in a confined space.First, the study derivative analytical stability and synchronization conditions for the asymmetrically coupled system in an unconfined space, which are then validated through numerical simulations. Simulation results show that asymmetric coupling has a significant impact on synchronization conditions. Moreover, it is observed that irrespective of whether the system is confined, an increase in coupling asymmetry leads to a hastened synchronization pace. Additionally, the study examines the effects of boundaries on the system's collective behaviors via numerical experiments. The presence of boundaries ensures the system's stability and synchronization, and reducing these boundaries can expedite the synchronization process and amplify its effects. Finally, the study reveals that the system's output amplitude exhibits stochastic resonance as the confined boundary size increases.

Multi-Physics Coupled Acoustic-Mechanics Analysis and Synergetic Optimization for a Twin-Fluid Atomization Nozzle

Fine particulate matter produced during the rapid industrialization over the past decades can cause significant harm to human health.Twin-fluid atomization technology is an effective means of controlling fine particulate matter pollution.In this paper,the influences of the main parameters on the droplet size,effective atomization range and sound pressure level(SPL)of a twin-fluid nozzle(TFN)are investigated,and in order to improve the atomization performance,a multi-objective synergetic optimization algorithm is presented.A multi-physics coupled acousticmechanics model based on the discrete phase model(DPM),large eddy simulation(LES)model,and Ffowcs Williams-Hawkings(FW-H)model is established,and the numerical simulation results of the multi-physics coupled acoustic-mechanics method are verified via experimental comparison.Based on the analysis of the multi-physics coupled acoustic-mechanics numerical simulation results,the effects of the water flow on the characteristics of the atomization flow distribution were obtained.A multi-physics coupled acoustic-mechanics numerical simulation result was employed to establish an orthogonal test database,and a multi-objective synergetic optimization algorithm was adopted to optimize the key parameters of the TFN.The optimal parameters are as follows:A gas flow of 0.94 m^(3)/h,water flow of 0.0237 m^(3)/h,orifice diameter of the self-excited vibrating cavity(SVC)of 1.19 mm,SVC orifice depth of 0.53 mm,distance between SVC and the outlet of nozzle of 5.11 mm,and a nozzle outlet diameter of 3.15 mm.The droplet particle size in the atomization flow field was significantly reduced,the spray distance improved by 71.56%,and the SPL data at each corresponding measurement point decreased by an average of 38.96%.The conclusions of this study offer a references for future TFN research.

Seismic performance evaluation of hybrid coupled shear wall system with shear and flexural fuse-type steel coupling beams

Replaceable flexural and shear fuse-type coupling beams are used in hybrid coupled shear wall(HCSW)systems,enabling concrete buildings to be promptly recovered after severe earthquakes.This study aimed to analytically evaluate the seismic behavior of flexural and shear fuse beams situated in short-,medium-and high-rise RC buildings that have HCSWs.Three building groups hypothetically located in a high seismic hazard zone were studied.A series of 2D nonlinear time history analyses was accomplished in OpenSees,using the ground motion records scaled at the design basis earthquake level.It was found that the effectiveness of fuses in HCSWs depends on various factors such as size and scale of the building,allowable rotation value,inter-story drift ratio,residual drift quantity,energy dissipation value of the fuses,etc.The results show that shear fuses better meet the requirements of rotations and drifts.In contrast,flexural fuses dissipate more energy,but their sectional stiffness should increase to meet other requirements.It was concluded that adoption of proper fuses depends on the overall scale of the building and on how associated factors are considered.

Coupled DEM-FDM analyses of the effects of falling rock’s shape and impact angle on response of granular cushion and rock shed

Rock shed is an effective protection measure against rockfall.To investigate the influences of falling rock’s shape and impact angle on the impact effect of the cushioned rock shed,a modeling approach for a rock shed with a cushion layer using PFC-FLAC.The granular cushion is modeled as an aggregate of discrete non-cohesion particles,while the concrete plate and the beam are modeled as zones.The falling rock with different sphericities and impact angles is modeled as a rigid assembly.The numerical model is validated by comparing the simulation results with experimental and numerical results from previous literature.This model is applied to analyze the effects of rock shape and impact angle on the dynamic interaction effects between falling rock and cushioned rock shed,including the impact force,transmitted bottom force,penetration depth,and plate deflection.The numerical results show that the variation in the falling rock’s shape has different effects on the falling rock with different impact angles.These findings could support rock shed design by revealing the limitations of the assumptions in the past research,which may result in unsafe rock sheds for some rockfall cases.

Stress-corrosion coupled damage localization induced by secondary phases in bio-degradable Mg alloys:phase-field modeling

In this study,a phase-field scheme that rigorously obeys conservation laws and irreversible thermodynamics is developed for modeling stress-corrosion coupled damage(SCCD).The coupling constitutive relationships of the deformation,phase-field damage,mass transfer,and electrostatic field are derived from the entropy inequality.The SCCD localization induced by secondary phases in Mg is numerically simulated using the implicit iterative algorithm of the self-defined finite elements.The quantitative evaluation of the SCCD of a C-ring is in good agreement with the experimental results.To capture the damage localization,a micro-galvanic corrosion domain is defined,and the buffering effect on charge migration is explored.Three cases are investigated to reveal the effect of localization on corrosion acceleration and provide guidance for the design for resistance to SCCD at the crystal scale.

Electrochemical reduction of carbon dioxide to produce formic acid coupled with oxidative conversion of biomass

Electrochemical reduction of CO_(2)(CO_(2)RR)has become a research hot spot in recent years in the context of carbon neutrality.HCOOH is one of the most promising products obtained by electrochemical reduction of CO_(2) due to its high energy value as estimated by market price per energy unit and wide application in chemical industry.Biomass is the most abundant renewable resource in the natural world.Coupling biomass oxidative conversion with CO_(2)RR driven by renewable electricity would well achieve carbon negativity.In this work,we comprehensively reviewed the current research progress on CO_(2)RR to produce HCOOH and coupled system for conversion of biomass and its derivatives to produce value-added products.Sn-and Bi-based electrocatalysts are discussed for CO_(2)RR with regards to the structure of the catalyst and reaction mechanisms.Electro-oxidation reactions of biomass derived sugars,alcohols,furan aldehydes and even polymeric components of lignocellulose were reviewed as alternatives to replace oxygen evolution reaction(OER)in the conventional electrolysis process.It was recommended that to further improve the efficiency of the coupled system,future work should be focused on the development of more efficient and stable catalysts,careful design of the electrolytic cells for improving the mass transfer and development of environment-friendly processes for recovering the formed formate and biomass oxidation products.

Effect of coil and chamber structure on plasma radial uniformity in radio frequency inductively coupled plasma

Enhancing plasma uniformity can be achieved by modifying coil and chamber structures in radio frequency inductively coupled plasma(ICP)to meet the demand for large-area and uniformly distributed plasma in industrial manufacturing.This study utilized a two-dimensional self-consistent fluid model to investigate how different coil configurations and chamber aspect ratios affect the radial uniformity of plasma in radio frequency ICP.The findings indicate that optimizing the radial spacing of the coil enhances plasma uniformity but with a reduction in electron density.Furthermore,optimizing the coil within the ICP reactor,using the interior point method in the Interior Point Optimizer significantly enhances plasma uniformity,elevating it from 56%to 96%within the range of the model sizes.Additionally,when the chamber aspect ratio k changes from 2.8 to 4.7,the plasma distribution changes from a center-high to a saddleshaped distribution.Moreover,the plasma uniformity becomes worse.Finally,adjusting process parameters,such as increasing source power and gas pressure,can enhance plasma uniformity.These findings contribute to optimizing the etching process by improving plasma radial uniformity.

Customized modulation on plasma uniformity by non-uniform magnetic field in capacitively coupled plasma

A two-dimensional fluid model based on COMSOL Multiphysics is developed to investigate the modulation of static magnetic field on plasma homogeneity in a capacitively coupled plasma(CCP)chamber. To generate a static magnetic field, direct current is applied to a circular coil located at the top of the chamber. By adjusting the magnetic field's configuration, which is done by altering the coil current and position, both the plasma uniformity and density can be significantly modulated. In the absence of the magnetic field, the plasma density exhibits an inhomogeneous distribution characterized by higher values at the plasma edge and lower values at the center. The introduction of a magnetic field generated by coils results in a significant increase in electron density near the coils. Furthermore, an increase in the sets of coils improves the uniformity of the plasma. By flexibly adjusting the positions of the coils and the applied current,a substantial enhancement in overall uniformity can be achieved. These findings demonstrate the feasibility of using this method for achieving uniform plasma densities in industrial applications.

Contrasts of bimodal tropical instability waves(TIWs)-induced wind stress perturbations in the Pacific Ocean among observations,ocean models,and coupled climate models

The coupling between wind stress perturbations and sea surface temperature(SST)perturbations induced by tropical instability waves(TIWs)in the Pacific Ocean has been revealed previously and proven crucial to both the atmosphere and ocean.However,an overlooked fact by previous studies is that the loosely defined“TIWs”actually consist of two modes,including the Yanai wave-based TIW on the equator(hereafter eTIW)and the Rossby wave-based TIW off the equator(hereafter vTIW).Hence,the individual feedbacks of the wind stress to the bimodal TIWs remain unexplored.In this study,individual coupling relationships are established for both eTIW and v TIW,including the relationship between the TIW-induced SST perturbations and two components of wind stress perturbations,and the relationship between the TIW-induced wind stress perturbation divergence(curl)and the downwind(crosswind)TIW-induced SST gradients.Results show that,due to different distributions of eTIW and vTIW,the coupling strength induced by the eTIW is stronger on the equator,and that by the vTIW is stronger off the equator.The results of any of eTIW and vTIW are higher than those of the loosely defined TIWs.We further investigated how well the coupling relationships remained in several widely recognized oceanic general circulation models and fully coupled climate models.However,the coupling relationships cannot be well represented in most numerical models.Finally,we confirmed that higher resolution usually corresponds to more accurate simulation.Therefore,the coupling models established in this study are complementary to previous research and can be used to refine the oceanic and coupled climate models.

CAS-ESM2.0 Successfully Reproduces Historical Atmospheric CO_(2) in a Coupled Carbon−Climate Simulation

The atmospheric carbon dioxide(CO_(2))concentration has been increasing rapidly since the Industrial Revolution,which has led to unequivocal global warming and crucial environmental change.It is extremely important to investigate the interactions among atmospheric CO_(2),the physical climate system,and the carbon cycle of the underlying surface for a better understanding of the Earth system.Earth system models are widely used to investigate these interactions via coupled carbon-climate simulations.The Chinese Academy of Sciences Earth System Model version 2(CAS-ESM2.0)has successfully fixed a two-way coupling of atmospheric CO_(2)with the climate and carbon cycle on land and in the ocean.Using CAS-ESM2.0,we conducted a coupled carbon-climate simulation by following the CMIP6 proposal of a historical emissions-driven experiment.This paper examines the modeled CO_(2)by comparison with observed CO_(2)at the sites of Mauna Loa and Barrow,and the Greenhouse Gases Observing Satellite(GOSAT)CO_(2)product.The results showed that CAS-ESM2.0 agrees very well with observations in reproducing the increasing trend of annual CO_(2)during the period 1850-2014,and in capturing the seasonal cycle of CO_(2)at the two baseline sites,as well as over northern high latitudes.These agreements illustrate a good ability of CAS-ESM2.0 in simulating carbon-climate interactions,even though uncertainties remain in the processes involved.This paper reports an important stage of the development of CAS-ESM with the coupling of carbon and climate,which will provide significant scientific support for climate research and China’s goal of carbon neutrality.

Abnormal transition of the electron energy distribution with excitation of the second harmonic in low-pressure radio-frequency capacitively coupled plasmas

The self-excited second harmonic in radio-frequency capacitively coupled plasma was significantly enhanced by adjusting the external variable capacitor.At a lower pressure of 3 Pa,the excitation of the second harmonic caused an abnormal transition of the electron energy probability function,resulting in abrupt changes in the electron density and temperature.Such changes in the electron energy probability function as well as the electron density and temperature were not observed at the higher pressure of 16 Pa under similar harmonic changes.The phenomena are related to the influence of the second harmonic on stochastic heating,which is determined by both amplitude and the relative phase of the harmonics.The results suggest that the self-excited high-order harmonics must be considered in practical applications of lowpressure radio-frequency capacitively coupled plasmas.

Different El Niño Flavors and Associated Atmospheric Teleconnections as Simulated in a Hybrid Coupled Model

A previously developed hybrid coupled model(HCM)is composed of an intermediate tropical Pacific Ocean model and a global atmospheric general circulation model(AGCM),denoted as HCMAGCM.In this study,different El Niño flavors,namely the Eastern-Pacific(EP)and Central-Pacific(CP)types,and the associated global atmospheric teleconnections are examined in a 1000-yr control simulation of the HCMAGCM.The HCMAGCM indicates profoundly different characteristics among EP and CP El Niño events in terms of related oceanic and atmospheric variables in the tropical Pacific,including the amplitude and spatial patterns of sea surface temperature(SST),zonal wind stress,and precipitation anomalies.An SST budget analysis indicates that the thermocline feedback and zonal advective feedback dominantly contribute to the growth of EP and CP El Niño events,respectively.Corresponding to the shifts in the tropical rainfall and deep convection during EP and CP El Niño events,the model also reproduces the differences in the extratropical atmospheric responses during the boreal winter.In particular,the EP El Niño tends to be dominant in exciting a poleward wave train pattern to the Northern Hemisphere,while the CP El Niño tends to preferably produce a wave train similar to the Pacific North American(PNA)pattern.As a result,different climatic impacts exist in North American regions,with a warm-north and cold-south pattern during an EP El Niño and a warm-northeast and cold-southwest pattern during a CP El Niño,respectively.This modeling result highlights the importance of internal natural processes within the tropical Pacific as they relate to the genesis of ENSO diversity because the active ocean–atmosphere coupling is allowed only in the tropical Pacific within the framework of the HCMAGCM.

A color image encryption scheme based on a 2D coupled chaotic system and diagonal scrambling algorithm

A novel color image encryption scheme is developed to enhance the security of encryption without increasing the complexity. Firstly, the plain color image is decomposed into three grayscale plain images, which are converted into the frequency domain coefficient matrices(FDCM) with discrete cosine transform(DCT) operation. After that, a twodimensional(2D) coupled chaotic system is developed and used to generate one group of embedded matrices and another group of encryption matrices, respectively. The embedded matrices are integrated with the FDCM to fulfill the frequency domain encryption, and then the inverse DCT processing is implemented to recover the spatial domain signal. Eventually,under the function of the encryption matrices and the proposed diagonal scrambling algorithm, the final color ciphertext is obtained. The experimental results show that the proposed method can not only ensure efficient encryption but also satisfy various sizes of image encryption. Besides, it has better performance than other similar techniques in statistical feature analysis, such as key space, key sensitivity, anti-differential attack, information entropy, noise attack, etc.

Dynamic Response Impact of Vehicle Braking on Simply Supported Beam Bridges with Corrugated Steel Webs Based on Vehicle-Bridge Coupled Vibration Analysis

A novel approach for analyzing coupled vibrations between vehicles and bridges is presented,taking into account spatiotemporal effects and mechanical phenomena resulting fromvehicle braking.Efficient modeling and solution of bridge vibrations induced by vehicle deceleration are realized using this method.The method’s validity and reliability are substantiated through numerical examples.A simply supported beam bridge with a corrugated steel web is taken as an example and the effects of parameters such as the initial vehicle speed,braking acceleration,braking location,and road surface roughness on the mid-span displacement and impact factor of the bridge are analyzed.The results show that vehicle braking significantly amplifies mid-span displacement and impact factor responses in comparison to uniform vehicular motion across the bridge.Notably,the influence of wheelto-bridge friction forces is of particular significance and cannot be overlooked.When the vehicle initiates braking near the middle of the span,both the mid-span displacement and impact factor of the bridge exhibit substantial increases,further escalating with higher braking acceleration.Under favorable road surface conditions,the midspan displacement and the impact factor during vehicle braking may exceed the design values stipulated by codes.It is important to note that road surface roughness exerts a more pronounced effect on the impact factor of the bridge in comparison to the effects of vehicle braking.

Effect of a negative DC bias on a capacitively coupled Ar plasma operated at different radiofrequency voltages and gas pressures

The effect of a negative DC bias,|V_(dc)|,on the electrical parameters and discharge mode is investigated experimentally in a radiofrequency(RF)capacitively coupled Ar plasma operated at different RF voltage amplitudes and gas pressures.The electron density is measured using a hairpin probe and the spatio-temporal distribution of the electron-impact excitation rate is determined by phase-resolved optical emission spectroscopy.The electrical parameters are obtained based on the waveforms of the electrode voltage and plasma current measured by a voltage probe and a current probe.It was found that at a low|V_(dc)|,i.e.inα-mode,the electron density and RF current decline with increasing|V_(dc)|;meanwhile,the plasma impedance becomes more capacitive due to a widened sheath.Therefore,RF power deposition is suppressed.When|V_(dc)|exceeds a certain value,the plasma changes toα–γhybrid mode(or the discharge becomes dominated by theγ-mode),manifesting a drastically growing electron density and a moderately increasing RF current.Meanwhile,the plasma impedance becomes more resistive,so RF power deposition is enhanced with|V_(dc)|.We also found that the electrical parameters show similar dependence on|V_(dc)|at different RF voltages,andα–γmode transition occurs at a lower|V_(dc)|at a higher RF voltage.By increasing the pressure,plasma impedance becomes more resistive,so RF power deposition and electron density are enhanced.In particular,theα–γmode transition tends to occur at a lower|V_(dc)|with increase in pressure.

Seismic safety assessment with non-Gaussian random processes for train-bridge coupled systems

Extensive high-speed railway(HSR)network resembled the intricate vascular system of the human body,crisscrossing mainlands.Seismic events,known for their unpredictability,pose a significant threat to both trains and bridges,given the HSR’s extended operational duration.Therefore,ensuring the running safety of train-bridge coupled(TBC)system,primarily composed of simply supported beam bridges,is paramount.Traditional methods like the Monte Carlo method fall short in analyzing this intricate system efficiently.Instead,efficient algorithm like the new point estimate method combined with moment expansion approximation(NPEM-MEA)is applied to study random responses of numerical simulation TBC systems.Validation of the NPEM-MEA’s feasibility is conducted using the Monte Carlo method.Comparative analysis confirms the accuracy and efficiency of the method,with a recommended truncation order of four to six for the NPEM-MEA.Additionally,the influences of seismic magnitude and epicentral distance are discussed based on the random dynamic responses in the TBC system.This methodology not only facilitates seismic safety assessments for TBC systems but also contributes to standard-setting for these systems under earthquake conditions.

Strongly Coupled Ag/Sn-SnO_(2)Nanosheets Toward CO_(2)Electroreduction to Pure HCOOH Solutions at Ampere‑Level Current

Electrocatalytic reduction of CO_(2) converts intermittent renewable electricity into value-added liquid products with an enticing prospect,but its practical application is hampered due to the lack of high-performance electrocatalysts.Herein,we elaborately design and develop strongly coupled nanosheets composed of Ag nanoparticles and Sn-SnO_(2) grains,designated as Ag/Sn-SnO_(2) nanosheets(NSs),which possess optimized electronic structure,high electrical conductivity,and more accessible sites.As a result,such a catalyst exhibits unprecedented catalytic performance toward CO_(2)-to-formate conversion with near-unity faradaic efficiency(≥90%),ultrahigh partial current density(2,000 mA cm^(−2)),and superior long-term stability(200 mA cm^(−2),200 h),surpassing the reported catalysts of CO_(2) electroreduction to formate.Additionally,in situ attenuated total reflection-infrared spectra combined with theoretical calculations revealed that electron-enriched Sn sites on Ag/Sn-SnO_(2)NSs not only promote the formation of*OCHO and alleviate the energy barriers of*OCHO to*HCOOH,but also impede the desorption of H*.Notably,the Ag/Sn-SnO_(2)NSs as the cathode in a membrane electrode assembly with porous solid electrolyte layer reactor can continuously produce~0.12 M pure HCOOH solution at 100 mA cm^(−2)over 200 h.This work may inspire further development of advanced electrocatalysts and innovative device systems for promoting practical application of producing liquid fuels from CO_(2).

Hydrologic Response to Future Climate Change in the Dulong-Irra-waddy River Basin Based on Coupled Model Intercomparison Project 6

Within the context of the Belt and Road Initiative(BRI)and the China-Myanmar Economic Corridor(CMEC),the Dulong-Ir-rawaddy(Ayeyarwady)River,an international river among China,India and Myanmar,plays a significant role as both a valuable hydro-power resource and an essential ecological passageway.However,the water resources and security exhibit a high degree of vulnerabil-ity to climate change impacts.This research evaluates climate impacts on the hydrology of the Dulong-Irrawaddy River Basin(DIRB)by using a physical-based hydrologic model.We crafted future climate scenarios using the three latest global climate models(GCMs)from Coupled Model Intercomparison Project 6(CMIP6)under two shared socioeconomic pathways(SSP2-4.5 and SSP5-8.5)for the near(2025-2049),mid(2050-2074),and far future(2075-2099).The regional model using MIKE SHE based on historical hydrologic processes was developed to further project future streamflow,demonstrating reliable performance in streamflow simulations with a val-idation Nash-Sutcliffe Efficiency(NSE)of 0.72.Results showed that climate change projections showed increases in the annual precip-itation and potential evapotranspiration(PET),with precipitation increasing by 11.3%and 26.1%,and PET increasing by 3.2%and 4.9%,respectively,by the end of the century under SSP2-4.5 and SSP5-8.5.These changes are projected to result in increased annual streamflow at all stations,notably at the basin’s outlet(Pyay station)compared to the baseline period(with an increase of 16.1%and 37.0%at the end of the 21st century under SSP2-4.5 and SSP5-8.5,respectively).Seasonal analysis for Pyay station forecasts an in-crease in dry-season streamflow by 31.3%-48.9%and 22.5%-76.3%under SSP2-4.5 and SSP5-8.5,respectively,and an increase in wet-season streamflow by 5.8%-12.6%and 2.8%-33.3%,respectively.Moreover,the magnitude and frequency of flood events are pre-dicted to escalate,potentially impacting hydropower production and food security significantly.This research outlines the hydrological response to future climate change during the 21st century and offers a scientific basis for the water resource management strategies by decision-makers.