Aiming at the problem of on-line damage diagnosis in structural health monitoring (SHM), an algorithm of feature extraction and damage alarming based on auto-regressive moving-average (ARMA) time series analysis i...Aiming at the problem of on-line damage diagnosis in structural health monitoring (SHM), an algorithm of feature extraction and damage alarming based on auto-regressive moving-average (ARMA) time series analysis is presented. The monitoring data were first modeled as ARMA models, while a principalcomponent matrix derived from the AR coefficients of these models was utilized to establish the Mahalanobisdistance criterion functions. Then, a new damage-sensitive feature index DDSF is proposed. A hypothesis test involving the t-test method is further applied to obtain a decision of damage alarming as the mean value of DDSF had significantly changed after damage. The numerical results of a three-span-girder model shows that the defined index is sensitive to subtle structural damage, and the proposed algorithm can be applied to the on-line damage alarming in SHM.展开更多
The first internal transcribed spacer(ITS1) of nuclear ribosomal DNA of three wild rice species and two subspecies of cultivated rice, which are distributed in China, was amplified using PCR technique and sequenced wi...The first internal transcribed spacer(ITS1) of nuclear ribosomal DNA of three wild rice species and two subspecies of cultivated rice, which are distributed in China, was amplified using PCR technique and sequenced with automated fluorescent sequencing. The sequences of ITS1 ranged from 193 bp to 218 bp in size and G/C content varied from 69.3% to 72.7%. In pairwise comparisons among the five taxa, sequence site divergence ranged from 1.5% to 10.6%. Phylogenetic analysis of ITS1 sequences using Wagner parsimony generated a single well resolved tree, which revealed that Oryza rufipogon was much more closely related to cultivated rice species than to the other two wild species. Oryza granulata was less closely related to either cultivated rice species or the other two wild species, and might be a unique and isolated taxon in the genus Oryza. The phylogenetic relationships of the three wild rice species and two cultivated rice subspecies inferred from ITS1 sequences is highly concordant with those based on the molecular evidence from isozyme, chloroplast DNA (cpDNA), mitochondrial DNA(mtDNA) and nuclear DNA (nDNA) of the genus Oryza .展开更多
Based on the changing law of municipal water demand,a trigonometric function model for short-term water demand forecast is established using the time-series analysis approach.The method for forecasting water demand du...Based on the changing law of municipal water demand,a trigonometric function model for short-term water demand forecast is established using the time-series analysis approach.The method for forecasting water demand during holidays and under unexpected events is also presented.Meanwhile,a computer software is developed.Through actual application,this method performs well and has high accuracy,so it can be applied to the daily operation of a water distribution system and lay a foundation for on-line optimal operation.展开更多
The research conducted prediction on changes of atmosphere pollution during July 9, 2014-July 22, 2014 with SPSS based on monitored data of O3 in 13 successive weeks from 6 sites in Baoding City and demonstrated predi...The research conducted prediction on changes of atmosphere pollution during July 9, 2014-July 22, 2014 with SPSS based on monitored data of O3 in 13 successive weeks from 6 sites in Baoding City and demonstrated prediction effect of ARIMA model is good by Ljung-Box Q-test and R2, and the model can be used for prediction on future atmosphere pollutant changes.展开更多
In the field of global changes, the relationship between plant phenology and climate, which reflects the response of terrestrial ecosystem to global climate change, has become a key subject that is highly concerned. U...In the field of global changes, the relationship between plant phenology and climate, which reflects the response of terrestrial ecosystem to global climate change, has become a key subject that is highly concerned. Using the moderate-resolution imaging spectroradiometer (MODIS)/enhanced vegetation index(EVI) collected every eight days during January- July from 2005 to 2008 and the corresponding remote sensing data as experimental materials, we constructed cloud-free images via the Harmonic analysis of time series (HANTS). The cloud-free images were then treated by dynamic threshold method for obtaining the vegetation phenology in green up period and its distribution pattern. And the distribution pattern between freezing disaster year and normal year were comparatively analyzed for revealing the effect of freezing disaster on vegetation phenology in experimental plot. The result showed that the treated EVI data performed well in monitoring the effect of freezing disaster on vegetation phenology, accurately reflecting the regions suffered from freezing disaster. This result suggests that processing of remote sensing data using HANTS method could well monitor the ecological characteristics of vegetation.展开更多
In this paper, we propose a cross reference method for nonlinear time series analyzing in semi blind case, that is, the dynamic equations modeling the time series are known but the corresponding parameters are not. ...In this paper, we propose a cross reference method for nonlinear time series analyzing in semi blind case, that is, the dynamic equations modeling the time series are known but the corresponding parameters are not. The tasks of noise reduction and parameter estimation which were fulfilled separately before are combined iteratively. With the positive interaction between the two processing modules, the method is somewhat superior. Some prior work can be viewed as special cases of this general framework. The simulations for noise reduction and parameter estimation of contaminated chaotic time series show improved performance of our method compared with previous work.展开更多
This paper presents a novel approach to identify and correct the gross errors in the microelectromechanical system (MEMS) gyroscope used in ground vehicles by means of time series analysis. According to the characte...This paper presents a novel approach to identify and correct the gross errors in the microelectromechanical system (MEMS) gyroscope used in ground vehicles by means of time series analysis. According to the characteristics of autocorrelation function (ACF) and partial autocorrelation function (PACF), an autoregressive integrated moving average (ARIMA) model is roughly constructed. The rough model is optimized by combining with Akaike's information criterion (A/C), and the parameters are estimated based on the least squares algorithm. After validation testing, the model is utilized to forecast the next output on the basis of the previous measurement. When the difference between the measurement and its prediction exceeds the defined threshold, the measurement is identified as a gross error and remedied by its prediction. A case study on the yaw rate is performed to illustrate the developed algorithm. Experimental results demonstrate that the proposed approach can effectively distinguish gross errors and make some reasonable remedies.展开更多
To protect trains against strong cross-wind along Qinghai-Tibet railway, a strong wind speed monitoring and warning system was developed. And to obtain high-precision wind speed short-term forecasting values for the s...To protect trains against strong cross-wind along Qinghai-Tibet railway, a strong wind speed monitoring and warning system was developed. And to obtain high-precision wind speed short-term forecasting values for the system to make more accurate scheduling decision, two optimization algorithms were proposed. Using them to make calculative examples for actual wind speed time series from the 18th meteorological station, the results show that: the optimization algorithm based on wavelet analysis method and improved time series analysis method can attain high-precision multi-step forecasting values, the mean relative errors of one-step, three-step, five-step and ten-step forecasting are only 0.30%, 0.75%, 1.15% and 1.65%, respectively. The optimization algorithm based on wavelet analysis method and Kalman time series analysis method can obtain high-precision one-step forecasting values, the mean relative error of one-step forecasting is reduced by 61.67% to 0.115%. The two optimization algorithms both maintain the modeling simple character, and can attain prediction explicit equations after modeling calculation.展开更多
The important effects of snow cover to ground thermal decades. In the most of previous research, the effects were usually regime has received much attention of scholars during the past few evaluated through the numeri...The important effects of snow cover to ground thermal decades. In the most of previous research, the effects were usually regime has received much attention of scholars during the past few evaluated through the numerical models and many important results are found. However, less examples and insufficient data based on field measurements are available to show natural cases. In the present work, a typical case study in Mohe and Beijicun meteorological stations, which both are located in the most northern tip of China, is given to show the effects of snow cover on the ground thermal regime. The spatial (the ground profile) and time series analysis in the extremely snowy winter of 2012-2013 in Heilongjiang Province are also performed by contrast with those in the winter of 2011-2012 based on the measured data collected by 63 meteorological stations, Our results illustrate the positive (warmer) effect of snow cover on the ground temperature (GT) on the daily basis, the highest difference between GT and daily mean air temperature (DGAT) is as high as 32.35℃. Moreover, by the lag time analysis method it is found that the response time of GT from 0 cm to 20 cm ground depth to the alternate change of snow depth has 10 days lag, while at 40 cm depth the response of DGAT is not significant. This result is different from the previous research by modeling, in which the resnonse denth of ground to the alteration of snow depth is far more than 40 cm.展开更多
Oil–water two-phase flow patterns in a horizontal pipe are analyzed with a 16-electrode electrical resistance tomography(ERT) system. The measurement data of the ERT are treated as a multivariate time-series, thus th...Oil–water two-phase flow patterns in a horizontal pipe are analyzed with a 16-electrode electrical resistance tomography(ERT) system. The measurement data of the ERT are treated as a multivariate time-series, thus the information extracted from each electrode represents the local phase distribution and fraction change at that location. The multivariate maximum Lyapunov exponent(MMLE) is extracted from the 16-dimension time-series to demonstrate the change of flow pattern versus the superficial velocity ratio of oil to water. The correlation dimension of the multivariate time-series is further introduced to jointly characterize and finally separate the flow patterns with MMLE. The change of flow patterns with superficial oil velocity at different water superficial velocities is studied with MMLE and correlation dimension, respectively, and the flow pattern transition can also be characterized with these two features. The proposed MMLE and correlation dimension map could effectively separate the flow patterns, thus is an effective tool for flow pattern identification and transition analysis.展开更多
Time-series analysis is important to a wide range of disciplines transcending both the physical and social sciences for proactive policy decisions. Statistical models have sound theoretical basis and have been success...Time-series analysis is important to a wide range of disciplines transcending both the physical and social sciences for proactive policy decisions. Statistical models have sound theoretical basis and have been successfully used in a number of problem domains in time series forecasting. Due to power and flexibility, Box-Jenkins ARIMA model has gained enormous popularity in many areas and research practice for the last three decades. More recently, the neural networks have been shown to be a promising alternative tool for modeling and forecasting owing to their ability to capture the nonlinearity in the data. However, despite the popularity and the superiority of ARIMA and ANN models, the empirical forecasting performance has been rather mixed so that no single method is best in every situation. In this study, a hybrid ARIMA and neural networks model to time series forecasting is proposed. The basic idea behind the model combination is to use each model’s unique features to capture different patterns in the data. With three real data sets, empirical results evidently show that the hybrid model outperforms ARIMA and ANN model noticeably in terms of forecasting accuracy used in isolation.展开更多
Electromagnetic emission(EME) is a kind of physical phenomenon accompanying the process of deformation and fracture of loaded coal and rock and it is of importance in quantitatively analyzing its characteristics.This ...Electromagnetic emission(EME) is a kind of physical phenomenon accompanying the process of deformation and fracture of loaded coal and rock and it is of importance in quantitatively analyzing its characteristics.This will reveal the process of deformation and fracture of coal and predicting dynamic disasters in coal mines.In this study,the G-P(Grassberger and Procaccia) algorithm,calculation steps of the(if only 1 dimension) correlation dimension of time series and the identification standards of chaotic signals are introduced.Furthermore,the correlation dimensions of EME and the acoustic emission(AE) signals of time series during deformation and fracture of coal bodies are calculated and analyzed.The results show that the time series of pulses number of EME and the time series of AE count rate are chaotic and that the saturation embedding dimensions of a K3 coal sample are,respectively,5 and 6.The results can be used to provide basic parameters for predicting of EME and AE time series.展开更多
Runoff change and trend of the Naoli River Basin were studied through the time series analysis using the data from the hydrological and meteorological stations. Time series of hydrological data were from 1957 to 2009 ...Runoff change and trend of the Naoli River Basin were studied through the time series analysis using the data from the hydrological and meteorological stations. Time series of hydrological data were from 1957 to 2009 for Bao′an station, from 1955 to 2009 for Baoqing station, from 1956 to 2009 for Caizuizi station and from 1978 to 2009 for Hongqiling station. The influences of climate change and human activities on runoff change were investigated, and the causes of hydrological regime change were revealed. The seasonal runoff distribution of the Naoli River was extremely uneven, and the annual change was great. Overall, the annual runoff showed a significant decreasing trend. The annual runoff of Bao′an, Baoqing, and Caizuizi stations in 2009 decreased by 64.1%, 76.3%, and 84.3%, respectively, compared with their beginning data recorded. The wet and dry years of the Naoli River have changed in the study period. The frequency of wet year occurrence decreased and lasted longer, whereas that of dry year occurrence increased. The frequency of dry year occurrence increased from 25.0%-27.8% to 83.9%-87.5%. The years before the 1970s were mostly wet, whereas those after the 1970s were mostly dry. Precipitation reduction and land use changes contributed to the decrease in annual runoff. Rising temperature and water project construction have also contributed important effects on the runoff change of the Naoli River.展开更多
文摘对带相关观测噪声和未知噪声统计的多传感器系统,用相关方法得到噪声统计在线估值器.在按分量标量加权线性最小方差最优信息融合准则下,用现代时间序列分析方法,基于滑动平均(moving average)新息模型的辨识,提出了自校正解耦融合Wiener预报器.用动态误差系统分析(dynamic error system analysis)方法证明了自校正融合Wiener预报器收敛于最优融合Wiener预报器,因而它具有渐近最优性.它的精度比每个局部自校正Wiener预报器精度都高.它的算法简单,便于实时应用.一个目标跟踪系统的仿真例子说明了其有效性.
基金The National High Technology Research and Devel-opment Program of China (863Program) (No2006AA04Z416)the National Natural Science Foundation of China (No50538020)
文摘Aiming at the problem of on-line damage diagnosis in structural health monitoring (SHM), an algorithm of feature extraction and damage alarming based on auto-regressive moving-average (ARMA) time series analysis is presented. The monitoring data were first modeled as ARMA models, while a principalcomponent matrix derived from the AR coefficients of these models was utilized to establish the Mahalanobisdistance criterion functions. Then, a new damage-sensitive feature index DDSF is proposed. A hypothesis test involving the t-test method is further applied to obtain a decision of damage alarming as the mean value of DDSF had significantly changed after damage. The numerical results of a three-span-girder model shows that the defined index is sensitive to subtle structural damage, and the proposed algorithm can be applied to the on-line damage alarming in SHM.
基金The project supported by the Grant from Presidentof Chinese Academy of Sciences
文摘The first internal transcribed spacer(ITS1) of nuclear ribosomal DNA of three wild rice species and two subspecies of cultivated rice, which are distributed in China, was amplified using PCR technique and sequenced with automated fluorescent sequencing. The sequences of ITS1 ranged from 193 bp to 218 bp in size and G/C content varied from 69.3% to 72.7%. In pairwise comparisons among the five taxa, sequence site divergence ranged from 1.5% to 10.6%. Phylogenetic analysis of ITS1 sequences using Wagner parsimony generated a single well resolved tree, which revealed that Oryza rufipogon was much more closely related to cultivated rice species than to the other two wild species. Oryza granulata was less closely related to either cultivated rice species or the other two wild species, and might be a unique and isolated taxon in the genus Oryza. The phylogenetic relationships of the three wild rice species and two cultivated rice subspecies inferred from ITS1 sequences is highly concordant with those based on the molecular evidence from isozyme, chloroplast DNA (cpDNA), mitochondrial DNA(mtDNA) and nuclear DNA (nDNA) of the genus Oryza .
基金Natural Science Foundation of China!(No.598780 30 )
文摘Based on the changing law of municipal water demand,a trigonometric function model for short-term water demand forecast is established using the time-series analysis approach.The method for forecasting water demand during holidays and under unexpected events is also presented.Meanwhile,a computer software is developed.Through actual application,this method performs well and has high accuracy,so it can be applied to the daily operation of a water distribution system and lay a foundation for on-line optimal operation.
基金Supported by Student Research Fund of Agricultural University of Hebei(cxzr2014023)Technology Fund of Agricultural University of Hebei(ZD201406)~~
文摘The research conducted prediction on changes of atmosphere pollution during July 9, 2014-July 22, 2014 with SPSS based on monitored data of O3 in 13 successive weeks from 6 sites in Baoding City and demonstrated prediction effect of ARIMA model is good by Ljung-Box Q-test and R2, and the model can be used for prediction on future atmosphere pollutant changes.
文摘In the field of global changes, the relationship between plant phenology and climate, which reflects the response of terrestrial ecosystem to global climate change, has become a key subject that is highly concerned. Using the moderate-resolution imaging spectroradiometer (MODIS)/enhanced vegetation index(EVI) collected every eight days during January- July from 2005 to 2008 and the corresponding remote sensing data as experimental materials, we constructed cloud-free images via the Harmonic analysis of time series (HANTS). The cloud-free images were then treated by dynamic threshold method for obtaining the vegetation phenology in green up period and its distribution pattern. And the distribution pattern between freezing disaster year and normal year were comparatively analyzed for revealing the effect of freezing disaster on vegetation phenology in experimental plot. The result showed that the treated EVI data performed well in monitoring the effect of freezing disaster on vegetation phenology, accurately reflecting the regions suffered from freezing disaster. This result suggests that processing of remote sensing data using HANTS method could well monitor the ecological characteristics of vegetation.
文摘In this paper, we propose a cross reference method for nonlinear time series analyzing in semi blind case, that is, the dynamic equations modeling the time series are known but the corresponding parameters are not. The tasks of noise reduction and parameter estimation which were fulfilled separately before are combined iteratively. With the positive interaction between the two processing modules, the method is somewhat superior. Some prior work can be viewed as special cases of this general framework. The simulations for noise reduction and parameter estimation of contaminated chaotic time series show improved performance of our method compared with previous work.
基金The National Natural Science Foundation of China(No.61273236)the Natural Science Foundation of Jiangsu Province(No.BK2010239)the Ph.D.Programs Foundation of Ministry of Education of China(No.200802861061)
文摘This paper presents a novel approach to identify and correct the gross errors in the microelectromechanical system (MEMS) gyroscope used in ground vehicles by means of time series analysis. According to the characteristics of autocorrelation function (ACF) and partial autocorrelation function (PACF), an autoregressive integrated moving average (ARIMA) model is roughly constructed. The rough model is optimized by combining with Akaike's information criterion (A/C), and the parameters are estimated based on the least squares algorithm. After validation testing, the model is utilized to forecast the next output on the basis of the previous measurement. When the difference between the measurement and its prediction exceeds the defined threshold, the measurement is identified as a gross error and remedied by its prediction. A case study on the yaw rate is performed to illustrate the developed algorithm. Experimental results demonstrate that the proposed approach can effectively distinguish gross errors and make some reasonable remedies.
基金Project(2006BAC07B03) supported by the National Key Technology R & D Program of ChinaProject(2006G040-A) supported by the Foundation of the Science and Technology Section of Ministry of RailwayProject(2008yb044) supported by the Foundation of Excellent Doctoral Dissertation of Central South University
文摘To protect trains against strong cross-wind along Qinghai-Tibet railway, a strong wind speed monitoring and warning system was developed. And to obtain high-precision wind speed short-term forecasting values for the system to make more accurate scheduling decision, two optimization algorithms were proposed. Using them to make calculative examples for actual wind speed time series from the 18th meteorological station, the results show that: the optimization algorithm based on wavelet analysis method and improved time series analysis method can attain high-precision multi-step forecasting values, the mean relative errors of one-step, three-step, five-step and ten-step forecasting are only 0.30%, 0.75%, 1.15% and 1.65%, respectively. The optimization algorithm based on wavelet analysis method and Kalman time series analysis method can obtain high-precision one-step forecasting values, the mean relative error of one-step forecasting is reduced by 61.67% to 0.115%. The two optimization algorithms both maintain the modeling simple character, and can attain prediction explicit equations after modeling calculation.
基金Under the auspices of National Natural Science Foundation of China(No.41471289,41301368)Natural Science Foundation of Jilin Province(No.20140101158JC)Foundation of State Key Laboratory of Remote Sensing Science(No.OFSLRSS201517)
文摘The important effects of snow cover to ground thermal decades. In the most of previous research, the effects were usually regime has received much attention of scholars during the past few evaluated through the numerical models and many important results are found. However, less examples and insufficient data based on field measurements are available to show natural cases. In the present work, a typical case study in Mohe and Beijicun meteorological stations, which both are located in the most northern tip of China, is given to show the effects of snow cover on the ground thermal regime. The spatial (the ground profile) and time series analysis in the extremely snowy winter of 2012-2013 in Heilongjiang Province are also performed by contrast with those in the winter of 2011-2012 based on the measured data collected by 63 meteorological stations, Our results illustrate the positive (warmer) effect of snow cover on the ground temperature (GT) on the daily basis, the highest difference between GT and daily mean air temperature (DGAT) is as high as 32.35℃. Moreover, by the lag time analysis method it is found that the response time of GT from 0 cm to 20 cm ground depth to the alternate change of snow depth has 10 days lag, while at 40 cm depth the response of DGAT is not significant. This result is different from the previous research by modeling, in which the resnonse denth of ground to the alteration of snow depth is far more than 40 cm.
基金Projects(61227006,61473206) supported by the National Natural Science Foundation of ChinaProject(13TXSYJC40200) supported by Science and Technology Innovation of Tianjin,China
文摘Oil–water two-phase flow patterns in a horizontal pipe are analyzed with a 16-electrode electrical resistance tomography(ERT) system. The measurement data of the ERT are treated as a multivariate time-series, thus the information extracted from each electrode represents the local phase distribution and fraction change at that location. The multivariate maximum Lyapunov exponent(MMLE) is extracted from the 16-dimension time-series to demonstrate the change of flow pattern versus the superficial velocity ratio of oil to water. The correlation dimension of the multivariate time-series is further introduced to jointly characterize and finally separate the flow patterns with MMLE. The change of flow patterns with superficial oil velocity at different water superficial velocities is studied with MMLE and correlation dimension, respectively, and the flow pattern transition can also be characterized with these two features. The proposed MMLE and correlation dimension map could effectively separate the flow patterns, thus is an effective tool for flow pattern identification and transition analysis.
文摘Time-series analysis is important to a wide range of disciplines transcending both the physical and social sciences for proactive policy decisions. Statistical models have sound theoretical basis and have been successfully used in a number of problem domains in time series forecasting. Due to power and flexibility, Box-Jenkins ARIMA model has gained enormous popularity in many areas and research practice for the last three decades. More recently, the neural networks have been shown to be a promising alternative tool for modeling and forecasting owing to their ability to capture the nonlinearity in the data. However, despite the popularity and the superiority of ARIMA and ANN models, the empirical forecasting performance has been rather mixed so that no single method is best in every situation. In this study, a hybrid ARIMA and neural networks model to time series forecasting is proposed. The basic idea behind the model combination is to use each model’s unique features to capture different patterns in the data. With three real data sets, empirical results evidently show that the hybrid model outperforms ARIMA and ANN model noticeably in terms of forecasting accuracy used in isolation.
基金Projects 50427401 supported by the National Natural Science Foundation of China2006BAK03B06 by the National Eleventh Five-Year Key Science & Technology Project of China+2 种基金the New Century Excellent Talent Program from the Ministry of Education (No.NCET-07-0799)the Fok Ying-Tong Education Foundation for Young Teachers in Higher Education Institutions of China (No.111053)the Beijing Science and Technology New Star Plan (No.2006A081)
文摘Electromagnetic emission(EME) is a kind of physical phenomenon accompanying the process of deformation and fracture of loaded coal and rock and it is of importance in quantitatively analyzing its characteristics.This will reveal the process of deformation and fracture of coal and predicting dynamic disasters in coal mines.In this study,the G-P(Grassberger and Procaccia) algorithm,calculation steps of the(if only 1 dimension) correlation dimension of time series and the identification standards of chaotic signals are introduced.Furthermore,the correlation dimensions of EME and the acoustic emission(AE) signals of time series during deformation and fracture of coal bodies are calculated and analyzed.The results show that the time series of pulses number of EME and the time series of AE count rate are chaotic and that the saturation embedding dimensions of a K3 coal sample are,respectively,5 and 6.The results can be used to provide basic parameters for predicting of EME and AE time series.
基金Under the auspices of National Natural Science Foundation of China (No. 40830535, 41001110, 41101092, 41171092)National Basic Research Program of China (No. 2010CB951304)the CAS/SAFEA (Chinese Academy of Sciences/State Administration of Foreign Experts Affairs) International Partnership Program for Creative Research Teams, Eleventh Five-Year' Key Technological Projects of Heilongjiang Province Farm Bureau (No. HNK10A-10-01, HNK10A-10-03)
文摘Runoff change and trend of the Naoli River Basin were studied through the time series analysis using the data from the hydrological and meteorological stations. Time series of hydrological data were from 1957 to 2009 for Bao′an station, from 1955 to 2009 for Baoqing station, from 1956 to 2009 for Caizuizi station and from 1978 to 2009 for Hongqiling station. The influences of climate change and human activities on runoff change were investigated, and the causes of hydrological regime change were revealed. The seasonal runoff distribution of the Naoli River was extremely uneven, and the annual change was great. Overall, the annual runoff showed a significant decreasing trend. The annual runoff of Bao′an, Baoqing, and Caizuizi stations in 2009 decreased by 64.1%, 76.3%, and 84.3%, respectively, compared with their beginning data recorded. The wet and dry years of the Naoli River have changed in the study period. The frequency of wet year occurrence decreased and lasted longer, whereas that of dry year occurrence increased. The frequency of dry year occurrence increased from 25.0%-27.8% to 83.9%-87.5%. The years before the 1970s were mostly wet, whereas those after the 1970s were mostly dry. Precipitation reduction and land use changes contributed to the decrease in annual runoff. Rising temperature and water project construction have also contributed important effects on the runoff change of the Naoli River.