Fourier Amplitude Research Articles

Empirical ground-motion models for horizontal Fourier amplitude spectra from fixed-effects and mixed-effects analyses of the NGA-West2 database

This article presents the development of ground-motion models (GMMs) of Fourier amplitude spectra for frequencies of 0.1–20 Hz and its potential extrapolation to 100 Hz at near-source distances using the effective amplitude spectrum (EAS) ordinates developed for the NGA-West2 project and the metadata and functional form from our previous NGA-West2 GMMs. We developed the GMMs using three different approaches to study the impact of including random effects in the model development: (1) fixed-effects regression (i.e. no random effects), (2) mixed-effects regression with events as a random effect, and (3) mixed-effects regression with both events and sites as random effects. Goodness-of-fit metrics show that the GMMs were improved with the addition of each random effect. We found the variance components other than the between-site standard deviation to be magnitude-dependent, which we estimated using Bayesian inference to incorporate uncertainty in the random-effects and within-group variability. As a result, our aleatory standard deviations are larger than residual-based standard deviations that ignore uncertainty in these terms. The GMMs predict a slight spectral sag at intermediate frequencies and large magnitudes that becomes more pronounced with increasing magnitude, consistent with other empirical analyses and ground-motion simulations available in the literature. We recommend use of the GMM that includes both event and site terms because of its appropriate modeling of repeatable effects. We present the other GMMs to demonstrate how each model is improved as random effects are added and to facilitate their comparison with other models that use a similar random-effects structure.

Vibration signal characteristics of debris flow initiation and motion experiments by considering clay content

Micro accelerometers are crucial tools for monitoring and providing early warning of geological disasters. This study focuses on the potential debris flow in Wangjiayuan, Fang County, Hubei Province, aiming to explore the relationship between clay content and debris flow motion characteristics. Physical model experiments were designed with varying clay content to simulate debris flow motion, and micro accelerometers were employed to monitor and record vibration acceleration changes during debris flow initiation and motion. Frequency spectra were generated using Fourier transform to analyze the motion characteristics in the frequency domain. The results showed that the predominant period of debris flow initiation ranged from 15 Hz to 29 Hz, while the period during motion in the middle section of the flow area ranged from 36 Hz to 54 Hz. Higher clay particle content was associated with a decrease in both the frequency range and mean of the predominant period, as well as a reduction in the Fourier amplitude extremum. These findings provide a foundation for monitoring and predicting debris flows, and the vibration signal characteristics of debris flow initiation and motion can be used to infer the clay content, offering insights into early warning systems for such geological hazards.

The effect of freezing temperature on the attenuation characteristics of stress waves in water-saturated granite in freezing environments

The investigation of stress wave propagation in frozen rock masses is of great significance for dynamic stability analysis of rock masses in cold regions. This paper investigates the attenuation of stress waves in water-saturated granite with different freezing temperatures. The water-saturated granite was frozen at different temperatures (−2°C, −10°C, −20°C, −30°C) and then impacted with different pendulum angles at those temperatures. The amplitude attenuation ratio, wave velocity and spectrum of stress waves in frozen granite were investigated. The effects of freezing temperature on the propagation coefficients (attenuation coefficient and wave number) and phase velocity of stress waves were investigated. Furthermore, the effects of the pendulum angle on the amplitude attenuation ratio, wave velocity, propagation coefficient and phase velocity at different freezing temperatures were discussed. The results show that the amplitude attenuation ratio of stress waves decreases as the temperature decreases. The wave velocity increases as the temperature decreases. The Fourier amplitudes of both the incident wave and the reflected wave decrease as the temperature decreases. The attenuation coefficient and wave number decrease sharply and then decrease slightly as the temperature decreases. The phase velocity increases as the temperature decreases. The amplitude attenuation ratio, wave velocity, phase velocity and propagation coefficients approximately keep constant under different pendulum angles. The present investigation can be applied to predict the amplitude attenuation and energy dissipation of stress waves within rock masses in freezing environments.

A regional site response analysis of the Petobo area after the 2018.09.28 Palu–Donggala Indonesia earthquake using an equivalent-linear model

This paper presents a site-specific seismic ground response evaluation through convolution–deconvolution analysis in the Balaroa–Petobo area during the 2018 Palu–Donggala Indonesia earthquake. The equivalent-linear ground response analysis for the earthquake time history recorded at Balaroa was carried out using DEEPSOIL software. The results of the analysis indicate that the EW component of the earthquake motion was amplified more severely (amax) than was the NS component, as it propagated to the Petobo surface. The amplification of the bedrock motion on the Petobo surface was more serious than that on the Balaroa surface, which appears to be due to the differences in the subsurface stratification and material properties of the two sites. The Fourier spectrum and response spectra also showed greater maximum spectral accelerations (Sa,max) and maximum Fourier amplitudes (Af) at the Petobo site than at the Balaroa site. The frequency of surface soil both the Petobo and Balaroa sites computed by using comparison between response spectra analysis and the local modes analysis VS/4*H was indicated the potential decline of surface soil stiffness at Petobo area appear to account for the structural damage and liquefaction flow slides during the 2018 incident.

Innovative segmentation technique for aerial power lines via amplitude stretching transform

Accurate segmentation of power line targets helps quickly locate faults, evaluate line conditions, and provides key image data support and analysis for the safe and stable operation of the power system.The aerial power line in segmentation due to the target is small, and the imaging reflected energy is weak, so the Unmanned Aerial Vehicle (UAV) aerial power line image is very susceptible to the interference of the environment line elements and noise, resulting in the detection of the power line target in the image of the defective, intermittent, straight line interferences and other low accuracy and real-time efficiency is not high. For this reason, this paper designs a pure amplitude stretching kernel function to form a Fourier amplitude vector field and uses this amplitude vector field to implement the stretching transformation of the amplitude field of the aerial power line image, so that the angular field after the Fourier inverse transformation can better react to the spatial domain line targets, and finally, after the Relative Total Variation (RTV) processing, the power line can be well detected. The proposed algorithm is compared with the main power line segmentation algorithms, such as Region Convolutional Neural Networks(R-CNN) and Phase Stretch Transform(PST). The average values of evaluation indicators PPA, MMPA and MMIoU of the image segmentation results of the proposed algorithm reach 0.96, 0.96 and 0.95 respectively, and the average time lag of detection is less than 0.2s, indicating that the accuracy and real-time performance of the segmentation results of the proposed algorithm are significantly better than those of the above algorithms.

Shaking Table Tests of Underground Substation Considering Soil–Structure–Equipment Interaction

ABSTRACTTo investigate the seismic dynamic response characteristics of underground substations, considering the interaction among soil, structure, and equipment, a model soil field shaking table test with a similarity ratio of 1:25 was conducted. Detailed analysis was conducted on the acceleration amplification factors, Fourier amplitude spectra, and displacements of soil, structure, and equipment during the experimental process. The experimental study revealed the following findings: It was concluded that the seismic excitation had a significant impact on the acceleration response of equipment near the edges of the structure, demonstrating an amplification effect. Both the acceleration amplification factors of the station structure and the site generally decreased with increasing input seismic intensity. As the seismic intensity increased, this decreasing trend gradually weakened. Additionally, the peak values of acceleration Fourier spectra often decrease with increasing burial depth. The examination of the postearthquake model system revealed the presence of surface cracks in the soil, which exhibited varying degrees of severity. These cracks were observed both parallel and perpendicular to the direction of excitation, with a significant concentration of soil fissures occurring directly above the structural elements.

Mapping near-surface S-wave attenuation in the onshore southeastern flank of the Red Sea Rifting system in Saudi Arabia

This study aims to quantify the site attenuation parameter, kappa (κ), in the southeastern flank of the Red Sea rifting system. To determine κ, a least-squares best-fit method was applied to three-component seismograms, allowing for the modelling of the Fourier amplitude spectra of S-waves. The dataset consisted of 53 earthquakes, ranging in local magnitudes from 3.0 to 4.8 ML. The estimates of κ for 17 sites in the southeastern flank of the Red Sea region revealed average values ranging from 0.026 ± 0.01s to 0.053 ± 0.01s. A statistical analysis of κ dependence on epicentral distance, near-surface geology, and earthquake magnitude demonstrated the dependence of κ on distance and surface geology rather than earthquake magnitudes. The spatial distribution of the parameter κ, which includes both path and site effects, shows higher values at soft rock sites in the eastern region than hard rock sites in the western area. Owing to surface geological conditions, the site-dependent parameter κo exhibited average values of 0.002 ± 0.003s and 0.024 ± 0.002s for the hard and soft rocks, respectively. The results may be used to improve seismic hazard assessments and earthquake ground motion simulations.

Anomalous diffusion model in Fourier space describing time correlation of shear Alfvénic turbulence

Time-correlation function in shear Alfvénic turbulence is examined from the point of view of stochastic dynamics in Fourier space. The complex random oscillator model for Elsasser variables, which has been used to discuss the violation of the Taylor hypothesis, is revisited. By using Fourier phase diffusion obeying the scaled Brownian motion, the generalized random oscillator model is derived. Auto- and cross-correlations given by the resultant model indicate that the characteristics (parameters) of anomalous diffusion of Fourier phase are closely related to the residual energy, while the cross-helicity is directly given by the energy portion among Elsasser variables. The stochastic uncertainty of the Fourier amplitude is also incorporated by taking the approach of the stochastic growth theory into account, resulting in the generalized model including terms that mimic both local and non-local interactions.

Parametric sensitivity analysis of the friction coefficient of coated textured surface in the rotary vane actuator end seal

This paper utilizes the extended Fourier amplitude sensitivity test (EFAST) method to analyze the effects of varying speeds on the friction coefficient of the coated textured end seal surface in rotary vane actuators. The analysis focuses on the first-order, total, and coupled sensitivities of root mean square (RMS) roughness, texture depth, and texture area ratio. The accuracy and validity of the numerical simulations are experimentally verified. At low speeds (20 rpm), RMS roughness shows the highest first-order sensitivity, 1.247 and 2.181 times greater than texture depth and texture area ratio. Using Support Vector Regression and Particle Swarm Optimization, the optimal parameters were determined: an RMS roughness ofs 0.1012 μm, a texture area ratio of 77.99%, and a texture depth of 2.9601 μm. The model can also predict the friction coefficient for various parameter combinations at different speeds.

Analysis of Depth Spatial Variations of Accelerograms Recorded from the Center for Engineering Strong Motion Data: Numerical Simulation for Multipoint Ground Motions

Studies on depth spatial variations are rare because of the limited number of measured recordings of the subsurface along the depth direction from the ground surface. However, these depth variations can directly influence the seismic dynamic response of underground structures or the foundations of space structures through multi-support excitation. In this study, depth spatial variations, including phase and amplitude, are analyzed to describe ground motions along different depths based on surface-earthquake time histories. The observed recordings from the Center for Engineering Strong Motion Data are used to investigate the spatial variation characteristics of the phase and amplitude in the depth direction, represented by the lagged coherence and Fourier amplitude spectrum ratio, respectively. A corresponding model is established under various site conditions to simulate multipoint correlated ground motions along the depth from the surface. Finally, a numerical simulation of multipoint ground motions is conducted to verify the application of the established models. This study provides an effective tool for describing spatial variations in both phase and amplitude, which are crucial for simulating ground motions at various depths.

An identification for channel mislabel of strong motion records based on Siamese neural network

Strong motion records are first-hand data for studying the seismic response of sites or engineering structures, and their objectivity is crucial for the credibility of the results in earthquake engineering and engineering seismology. However, domestic and international earthquake data may be mislabeled between horizontal and vertical channels. This issue is typically addressed by manually comparing the similarity between the three components of strong motion records, which is inherently subjective and inefficient in identification. To achieve the intelligent recognition of massive records, this study used 14,983 sets of ground motion records with significant differences between horizontal and vertical components from the NGA-West2 database. A Siamese neural network preliminarily distinguished the similarity between the acceleration waveform and the three components of the Fourier amplitude spectrum (FAS) of ground motion records. Combined with manual identification, an efficient and accurate method for identifying vertical components in ground motion records was proposed, and applied to verify the channel directions of the strong motion records in Strong Motion Network in China. It was found that 308 sets of records from 170 stations were suspected of mislabeling vertical and horizontal components. This advancement significantly enhances the objectivity of strong motion records. This proposed method holds potential for remote maintenance of strong motion stations, verifying the channels of strong motion instruments, and mitigating the negative impact of channel confusion on research results.

Variability and uncertainty of variance components in fixed-effects and mixed-effects ground-motion models

It has become common in the development of ground-motion models (GMMs), especially using mixed-effects regression with crossed random effects, to calculate standard deviations, referred to as variance components, from sample statistics of the residuals (i.e. random effects and within-group errors) rather than using the variance components reported by a mixed-effects regression program, calculated using the same algorithms used in the mixed-effects regression, or estimated using Bayesian inference. This practice leads to underestimating the standard deviations because it does not account for the uncertainty (i.e. standard errors) associated with fitting the random effects and within-group errors during the regression analysis. In this study, we used a series of ground-motion models for Fourier amplitude spectra developed using mixed-effects regression of an Next Generation Attenuation (NGA)-West2 database to show that residual-based standard deviations can be significantly smaller than regression-based standard deviations. These differences are exacerbated for those variance components with the fewest number of observations or when the residuals are partitioned (e.g. by earthquake magnitude). Using residual-based variance components not only results in smaller standard deviations but can lead to biased inferences when attempting to compare the efficacy of a model based solely on a comparison of the values of the variance components or their related mean square errors. It can also lead to biased fixed-effects coefficients if the variance components are derived from the residuals of a mixed-effects regression rather than estimated during a regression. We show that variance components can be decomposed into various random-effect grouping factors and partitioned into subsets of predictor variables, such as magnitude, from the total residuals of a non-Bayesian mixed-effects regression program using Bayesian inference. This process does not require the development of an entire GMM using Bayesian mixed-effects regression.

Large-scale Cosmic-ray Anisotropies with 19 yr of Data from the Pierre Auger Observatory

Results are presented for the measurement of large-scale anisotropies in the arrival directions of ultra–high-energy cosmic rays detected at the Pierre Auger Observatory during 19 yr of operation, prior to AugerPrime, the upgrade of the observatory. The 3D dipole amplitude and direction are reconstructed above 4 EeV in four energy bins. Besides the established dipolar anisotropy in R.A. above 8 EeV, the Fourier amplitude of the 8–16 EeV energy bin is now also above the 5σ discovery level. No time variation of the dipole moment above 8 EeV is found, setting an upper limit to the rate of change of such variations of 0.3% yr−1 at the 95% confidence level. Additionally, the results for the angular power spectrum are shown, demonstrating no other statistically significant multipoles. The results for the equatorial dipole component down to 0.03 EeV are presented, using for the first time a data set obtained with a trigger that has been optimized for lower energies. Finally, model predictions are discussed and compared with observations, based on two source emission scenarios obtained in the combined fit of spectrum and composition above 0.6 EeV.

Two-step global sensitivity analysis of a non-local integro-differential model for Cancer-on-Chip experiments

The present work focuses on a non-local integro-differential model reproducing Cancer-on-chip experiments where tumor cells, treated with chemotherapy drugs, secrete chemical signals stimulating the immune response. The reliability of the model in reproducing the phenomenon of interest is investigated through a global sensitivity analysis, rather than a local one, to have global information about the role of parameters, and by examining potential non-linear effects in greater detail. Focusing on a region in the parameter space, the effect of 13 model parameters on the in silico outcome is investigated by considering 11 different target outputs, properly selected to monitor the spatial distribution and the dynamics of immune cells along the period of observation. In order to cope with the large number of model parameters to be investigated and the computational cost of each numerical simulation, a two-step global sensitivity analysis is performed. First, the screening Morris method is applied to rank the effect of the 13 model parameters on each target output and it emerges that all the output targets are mainly affected by the same 6 parameters. The extended Fourier Amplitude Sensitivity Test (eFAST) method is then used to quantify the role of these 6 parameters. As a result, the proposed analysis highlights the feasibility of the considered space of parameters, and indicates that the most relevant parameters are those related to the chemical field and cell-substrate adhesion. In turn, it suggests how to possibly improve the model description as well as the calibration procedure, in order to better capture the observed phenomena and, at the same time, reduce the complexity of the simulation algorithm. On one hand, the model could be simplified by neglecting cell–cell alignment effects unless clear empirical evidences of their importance emerge. On the other hand, the best way to increase the accuracy and reliability of our model predictions would be to have experimental data/information to reduce the uncertainty of the more relevant parameters.

Performance evaluation of the USGS velocity model for the San Francisco Bay Area

In this study, we evaluated the performance of the United States Geological Survey velocity model developed for the San Francisco Bay Area (SFBA), version 21.1. The evaluation was performed through high-resolution three-dimensional physics-based ground motion simulations of seven small-magnitude earthquakes (ranging from magnitude 3.8 to 4.4) that occurred on the eastern side of the San Francisco Bay. The simulations were performed in the frequency range from 0 to 5 Hz with a minimum shear-wave velocity of 250 m/s, which allowed the capture of wave propagation effects of the near-surface soft materials that characterize local basins. Based on the direct comparison of Fourier amplitude spectra between recorded and simulated ground motions for more than 250 stations, we found that the velocity model generally performs well in the frequency range of 0.2–5 Hz. The median value of the Fourier amplitude residuals was found to be near zero for all seven earthquakes. The slight over-prediction of 0.2 log-natural units at frequencies above 3 Hz in our simulations was attributed to the potentially inaccurate representation of the source radiation pattern by a double-couple point source model, and simple representation of shallow small-scale underground structural complexity in the velocity model. Maps of spectral amplitude differences between the simulated and recorded data were used to identify areas responsible for systematic ground motion over-predictions or under-predictions. For example, while some sub-domains over soft sediments show over-prediction patterns, the block east of the Hayward fault is prone to exhibit patterns of under-prediction. These maps can be used to guide future refinements of the SFBA velocity model. Since our simulation methodology allows for the decoupling of the source and wave propagation effects, the ground motion data generated by our simulations can also be used to quantify the epistemic uncertainty due to the velocity model, in empirically based ground motion estimates for the SFBA.

Atomic-scale imaging of laser-driven electron dynamics in solids

Resolving laser-driven electron dynamics on their natural time and length scales is essential for understanding and controlling light-induced phenomena. Capabilities to reveal these dynamics are limited by challenges in interpreting wave mixing of a driving and a probe pulse, low energy resolution at ultrashort time scales and a lack of atomic-scale resolution by standard spectroscopic techniques. Here, we demonstrate how ultrafast x-ray diffraction can access fundamental information on laser-driven electronic motion in solids. We propose a method based on subcycle-resolved x-ray-optical wave mixing that allows for a straightforward reconstruction of key properties of strong-field-induced electron dynamics with atomic spatial resolution. Namely, this technique provides both phases and amplitudes of the spatial Fourier transform of optically-induced charge distributions, their temporal behavior, and the direction of the instantaneous microscopic optically-induced electron current flow. It captures the rich microscopic structures and symmetry features of laser-driven electronic charge and current density distributions.

Precise Fourier parameters of Cepheid radial velocity curves: Towards refining the Hertzsprung progression models

Context. Radial velocity (RV) curves of Classical Cepheids allow precise determination of the resonant periods, which in turn help to constrain fundamental parameters of these stars. The RV curves of Cepheids are also useful for identifying their pulsation modes and for distance determination using the parallax-of-pulsation method. Aims. The primary goal of this paper is to derive precise Fourier parameters of the RV curves for fundamental and first-overtone Galactic Cepheids. Our secondary objectives are then to analyze the progression of the Fourier parameters up to the seventh harmonic, and to propose an identification of the pulsation modes of the stars. Methods. For each star, we carefully selected RV measurements available in the literature that yield the highest precision of Fourier parameters according to the procedure that follows. We performed a Fourier decomposition of the RV curves using the unweighted least-square method and the standard deviation of the fit was used to derive the uncertainty on the Fourier parameters. We corrected for zero-point differences between datasets and RV modulations caused by binary motion. Results. With this study we have more than doubled the number of Cepheids with published RV curve Fourier parameters and with their uncertainty properly estimated. Our sample includes 178 fundamental-mode and 33 first-overtone pulsators, as well as 7 additional Cepheids whose pulsation mode is uncertain or undetermined according to our criteria. For the fundamental-mode Cepheids, the precision of the obtained low-order Fourier phases and amplitudes is about seven times and 25% better, respectively, as compared to the precision achieved in previously published Fourier parameter surveys. With highly accurate RV Fourier phases ϕ21, we are able to firmly identify V495 Cyg as a new first-overtone Cepheid and we confirm the first-overtone nature of several other stars. In particular, α UMi should be firmly classified as a first-overtone pulsator. In three objects (VY Per, AQ Pup, and QZ Nor), we find significant γ-velocity variations, which for the first two objects (and possibly for QZ Nor as well) can be attributed to the spectroscopic binarity of these stars. Finally, the analysis of the fundamental mode Fourier parameters up to seventh order reveals tight progression of Fourier phases for all pulsation periods. Conclusions. We provide new precise Fourier parameters of Cepheid RV curves determined from RV measurements available in the literature together with unpublished data. The pulsation period coverage and the precision obtained, in particular for Fourier phase ϕ21, will be useful for studying the dynamics of Cepheid pulsations with the help of hydrodynamical models. Further RV measurements from modern high-resolution spectroscopic instruments will be important to improve these results.

Insights on CDI parametric controls and dependencies using gloabal sensitivity analysis

Capacitive deionization (CDI) emerges as a promising desalination technology due to its high energy efficiency and potential for energy recovery. Its performance improvement and application extension highly require understanding how variations in multiple input parameters, as well as their interactions, lead to variations in output, which is not straightforward in the complex system of CDI. In this work we use the extended Fourier amplitude sensitivity test (eFAST) method to conduct a global sensitivity analysis (GSA) in order to examine the effects of several parameters, including material properties, operating conditions, and ionic properties, on the energy consumption and salt removal efficiency of CDIs. The research reveals that high capacitance and efficient mass transfer are crucial for CDI advancement, with a trade-off between desalination rate and energy efficiency. The interplay between input parameters is revealed by factorial design and GSA, which show that no single parameter alone determines salt adsorption capacity but that feedwater salt concentration and current density are important when considering both individual and coupled effects. The study suggests that the specific application and feedwater characteristics should be taken into consideration when deciding between the CC and CV modes of desalination. While CC mode excels in water production and energy efficiency for high salinity feedwater pretreatment, CV mode is better at treating low salinity waterr. The comprehensive parametric analysis provides valuable insights into optimizing CDI design and operating conditions to balance desalination capacity and energy efficiency.

Isotropic High‐Frequency Radiation in Near‐Fault Seismic Data

AbstractWe compare Fourier Amplitude Spectra of Fault Normal (FN) and Fault Parallel (FP) seismograms at near‐fault sites for seven strike‐slip earthquakes with moment magnitudes Mw ≥ 6. For all events we find large FN/FP ratios at low frequencies consistent with near‐fault S‐wave radiation patterns for strike‐slip earthquakes. However, the difference diminishes with increasing frequency and FN/FP is about 1 above a transition frequency. The results may reflect small tensile/isotropic components in the earthquake rupture zones that homogenize the high‐frequency radiation in different directions at near‐fault sites. The FN/FP ratios at low frequencies and transition frequencies above which FN ∼ FP vary among the analyzed earthquakes and have no clear correlation with the magnitudes. The lack of correlation may signify a characteristic scale (e.g., process zone size, duration of source time function) controlling the isotropic radiation, and/or wave propagation and other effects that mask the source effects.

Regional water cycle health evaluation and obstacle factor analysis: A case study of the Yellow River basin in Henan Province, China

Ecological and socioeconomic systems are interconnected through water cycles. Due to intensive anthropogenic activities and global climate change, the water cycle has undergone profound changes over the last few decades. This study aimed to evaluate the health status of the regional water cycle, explore the etiology of water cycle disorders, and provide suggestions for improving the health of the water cycle. Based on the Driving force-Pressure-State-Impact-Response (DPSIR) model, an index system was constructed to evaluate regional water cycle health. The water cycle was evaluated using a random forest (RF) model, which can achieve desirable fitting and generalization properties. The Yellow River Basin in Henan Province (YRBHN), China, was used as an example to demonstrate the application of this system. An Extended Fourier Amplitude Sensitivity Test (EFAST)-Cloud model was used to validate the evaluation results of the RF model. The main influencing factors in the water cycle process were then identified using the obstacle degree model, and ensuing policy recommendations were presented. The results demonstrated that the health state of the water cycle in the YRBHN shifted from unhealthy to sub-healthy during 2011–2020, with the response layer showing the most evident positive change. The health states of Zhengzhou and Luoyang were superior to those of the other cities, with average health index numbers of 3.1006 and 3.0078, respectively. Although the evaluation index numbers of the EFAST-Cloud and RF models were slightly different, their evaluation results were similar, indicating that the RF model was applicable for water cycle health evaluation. The proportion of ecological water use and the response (R) layer provided the greatest obstacle degree in the index and criterion layers, respectively. Based on the current state of the water cycle in the YRBHN, some policy implications were proposed. This study proposed a new method for regional water cycle health evaluation and decision making that provides theoretical support for regional water resource management and is of considerable practical significance.