Modeling the distribution of Populus euphratica in the Heihe River Basin,an inland river basin in an arid region of China

Populus euphratica is a dominant tree species in riparian Tugai forests and forms a natural barrier that maintains the stability of local oases in arid inland river basins. Despite being critical information for local environmental protection and recovery, establishing the specific spatial distribution of P. euphratica has rarely been attempted via precise and reliable species distribution models in such areas. In this research, the potential geographic distribution of P. euphratica in the Heihe River Basin was simulated with MaxEnt software based on species occurrence data and 29 environmental variables. The result showed that in the Heihe River Basin, 820 km~2 of land primarily distributed along the banks of the lower reaches of the river is a suitable habitat for P. euphratica. We built other MaxEnt models based on different environmental variables and another eight models employing different mathematical algorithms based on the same 29 environmental variables to demonstrate the superiority of this method.MaxEnt based on 29 environmental variables performed the best among these models, as it precisely described the essential characteristics of the distribution of P. euphratica forest land. This study verified that MaxEnt can serve as an effective tool for species distribution in extremely arid regions with sufficient and reliable environmental variables. The results suggest that there may be a larger area of P. euphratica forest distribution in the study area and that ecological conservation and management of P.euphratica should prioritize suitable habitat. This research provides valuable insights for the conservation and management of degraded P. euphratica riparian forests.     

Variation of <Emphasis Type="Italic">Q</Emphasis> value before and after the 1999 Xiuyan,Liaoning Province, <Emphasis Type="Italic">M</Emphasis>=5.4 earthquake on the basis of analysis on attenuation dispersion of P waves

A method for determining medium quality factor is developed on the basis of analyzing the attenuation dispersion of the arrived first period P wave. In order to enhance signal to noise ratio, improve the resolution in measurement and reduce systematic error we applied the data resampling technique. The group velocity delay of P wave was derived by using an improved multi-filtering method. Based on a linear viscoelastic relaxation model we deduced the medium quality factor Q _m, and associated error with 95% confidence level. Applying the method to the seismic record of the Xiuyan M=5.4 earthquake sequences we obtained the following result: (1) High Q _m started to appear from Nov. 9, 1999. The events giving the deduced high Q _m value clustered in a region with their epicenter distances being between 32 and 46 km to the Yingkou station. This Q _m versus distance observation obviously deviates from the normal trend of Q _m linearly increasing with distance. (2) The average Q _m before the 29 Dec. 1999 M=5.4 earthquake is 460, while the average Q _m between the M=5.4 event and the 12 Jan. 2000 M=5.1 earthquake is 391, and the average Q _m after the M=5.1 event is 204.     

Critical frequencies <Emphasis Type="Italic">foF</Emphasis>2 as an indicator of trends in thermospheric dynamics

The time variations in three parameters during the last decades are considered. R(foF2) is the correlation coefficient between the nighttime and daytime values of foF2 for the same day. Stable trends are found for the minimum (R(foF2)(max)) and maximum (R(foF2)(min)) values of R(foF2) during a year. The foF2(night)/foF2(day) ratio demonstrates both, negative and positive trends, and the trend sign depends on the inclination I and declination D of the magnetic field. The correlation coefficient r(h, fo) between foF2 and the 100 hP level in the stratosphere demonstrates a decrease (in the years of maximum and minimum solar activity) from the 1980s to the 1990s. The trends in all three groups of data are considered under the assumption of long-term changes in the circulation in the upper atmosphere.     

Behavior of Parameters of Enhancements in the Nighttime Electron Concentration in the Ionospheric <Emphasis Type="Italic">F</Emphasis>2 Layer

We have analyzed the behavior of the F2 layer parameters during nighttime periods of enhanced electron concentration by the results of vertical sounding of the ionosphere carried out with five-minute periodicity in Almaty (76°55′ E, 43°15′ N) in 2001–2012. The results are obtained within the frameworks of the unified concept of different types of ionospheric plasma disturbances manifested as variations in the height and half-thickness of the layer accompanied by an increase and decrease of N _m F2 at the moments of maximum compression and expansion of the layer. A good correlation is found between height h _Am , which corresponds to the maximum increase, and layer peak height h _m F, while h _Am is always less than h _m F. The difference between h _Am and h _m F linearly increases with increasing h _m F. Whereas the difference is ~38 km for h _m F = 280 km, it is ~54 km for h _m F = 380 km. Additionally, the correlation is good between the increase in the electron concentration in the layer maximum ΔN _m and the maximum enhancement at the fixed height ΔN; the electron concentration enhancement in the layer maximum is about two to three times lower than its maximum enhancement at the fixed height.     

Alcohol compounds in <Emphasis Type="Italic">Azolla imbricata</Emphasis> and potential source implication for marine sediments

This study investigated the composition of long-chain alkyl diols, triols, sec-alcohols, hydroxyl acids, and other hydroxylated compounds in Azolla imbricata and compared the organic alcohol components of Azolla filiculoides, Azolla microphylla, and South China Sea (SCS) sediments in order to investigate the possible indication of Azolla being the biological source of diols and triols in SCS sediment. Large amounts of diols, monohydroxy acids, and sec-alcohols with internal hydroxy groups at ω20 were detected in the three types of Azolla. Among these, 1,ω20-diol and ω20-hydroxy acid exhibited strong even-odd predominance distribution, whereas ω20-sec-alcohol exhibited strong odd-even predominance distribution. In addition, small amounts of diols, triols, and dihydroxy acids with internal hydroxy groups at 9, 10 or ω9, ω10 were detected, among which the chain length of C₂₉ was predominate. Compounds having similar structures as those in Azolla reflected a similar biosynthetic pathway: ω20-hydroxy acid exhibiting even-odd predominance distribution is decarboxylated to ω20-sec-alcohol exhibiting odd-even predominance distribution and converted to 1, ω20-diol with even-odd predominance distribution by acyl reduction; ω9, ω10-hydroxy acid is converted to 1,20,21(1, ω9, ω10)-triol by acyl reduction, and then converted to 9,10-diol by hydrogenation and dehydration. The alcohol components in A. imbricata were clearly not the biological source of 1,13/1,14/1,15-C_{28, 30, 32} diols and 1,3,4-C_27-29 triols in the SCS sediment. Certain marine diatoms might be the source of 1,14-C_{28, 30} diol in inshore sediment, but the biological source of diols and triols in the SCS sediment requires further investigation.     

Carbon assimilation characteristics of plants in oasis-desert ecotone and their response to CO<Subscript>2</Subscript> enrichment

Six species of more than 20-year-old desert woody plants in the oasis-desert ecotone were selected for study. The results showed that: (1) in different growing seasons δ¹³C values of assimilating organ varied between -14%. and -16%. for Haloxylon ammodendron (HA),-14%. - -15%. for Calligonum mongolicum (CM) and -25%. - -28%. for Caragana korshinskii (CK), Nitraria sphaerocarpa (NS) and Hedysarum scoparium (HS). (2) The net photosynthetic rate (Pn) of HA and CM was significantly higher than those of the other species. With the decrease in Pn for the six species, their intercellular CO₂ concentration increased, but stomatal limitation value decreased under the intensive light. At the same time, the photochemical efficiency of PS II dropped to different degrees. (3) The CO₂ enrichment experiment demonstrated that, Pn of HA and CM increased to different extent under 450 μmol/mol, but their Pn reduced or approximated to the current condition under 650 μmol/mol. Under 450 μmol/mol the efficiency for solar energy utilization of CK and HS significantly reduced and under 650 μmol/mol their respiration rate exceeded photosynthesis rate. It can be concluded that HA and CM have some function of pathway for C₄ but the other three species have the function for C₃. The decline in their Pn is mainly caused by non-stomatal factors. HA, CM, CK and HS exhibited photoinhibition, which disappeared in a short time. This is a kind of positive readjustment to adapting to the desert environment. HA and CM can adapt to the high CO₂ environment, but CK and HS cannot. With the rise in atmospheric CO₂ concentration and climate warming, the latter two species in the oasis-desert ecotone may be gradually degraded or even disappear.     

An experimental study on variation of <Emphasis Type="Italic">b</Emphasis> value during deformation of rock samples with different structural models

A preliminary study of b value of rocks with two kinds of structural models has been made on the base of a new acoustic emission recording system. It shows that b value of the sample decreases obviously when the sample with compressive en echelon faults changes into a tensile one after interchange occurs between stress axis σ ₁ and σ ₂. A similar experiment is observed when the sample with tensile en echelon faults changes into that with a bend fault after two segments of the en echelon fault linking up. These facts indicate that the variation of b value may contain the information of the regional dominant structural model. Therefore, b-value analyses could be a new method for studying regional dominant structural models.     

Estimation of tree transpiration and response of tree conductance to meteorological variables in desert-oasis system of Northwest China

Transpirations of three dominated tree species, namely Mongol Scotch Pine (Pinus sylvestris var. mongolica Litvin), White elm (Ulmus pumila) and Gansu Poplar (Populus gansuensis Wang et Yang) in oasis shelter forest (Linze site) and of two dominated tree species, namely Euphrates Poplar (Populus euphratica Oliv.) and Russia olive (Elaeagnus angustifolia Linn.) in lowland desert (Erjinaqi site) have been estimated using measured sapflow in summer, autumn and winter, 2002 and in spring, 2003. An ENVIS System was used for each site to measure microclimate variables, soil moisture and sapflow every half an hour, and the study time scale is one day. In the 104 days of observation during the growing season at the Linze site, the average daily sapflow of Gansu Poplar is 9.93 L· d^-1, and the average transpiration per unit leaf area is 1.99 mm · d^-1. For White elm tree, the daily average sapflow is 4.08 L· d^-1, while the daily average transpiration per unit leaf area is 0.49 mm · d^-1. The values for Mongol Scotch Pine are 3.91 L· d^-1 and 0.25 mm · d^-1, respectively. In the total 73 days of observation during the growing season at the Erjinaqi site, the daily average sapflows of Russia olive and Euphrates Poplar are 12.1 and 20.97 L· d^-1, respectively, and the average transpirations per unit leaf area are 0.22 and 0.31 mm · d^-1, respectively. In the observation period of the growing season, tree conductances of Mongol Scotch Pine, White elm, Gansu Poplar or Russia olive show an exponential relationship with the daily average air temperature or vapour pressure deficit, but the relationship is not so obvious between tree conductance and global radiation. The transpiration process of each tree species is affected by all the observed four environmental variables. The response of tree conductance to different climatic factors changes with tree species. The effect of the same factor to the same tree species is also variable in different growing stages. The sapflow of every tree species is relatively large in later spring to early summer, and low in summer, and then reaches its largest value in later September. In the mid-November, the sapflow is relatively large, especially the deciduous tree species. This may be characteristic of the tree species in Arid Regions of Northwest China.     

Sampling uncertainties and source <Emphasis Type="Italic">b</Emphasis> likelihood for the Gutenberg-Richter <Emphasis Type="Italic">b</Emphasis> value from the Aki-Utsu method

The Aki-Utsu method of Gutenberg-Richter (G-R) b value estimation is often misapplied so that estimations not using the G-R histogram are often meaningless because they are not based on adequate samples. We propose a method to estimate the likelihood Pr(b?b _m , N, M ₁, M ₂) that an observed b _m estimate, based on a sample of N magnitudes within an [M ₁????≤?ΔM/2,?M ₂?+?ΔM/2) range, where ΔM?=?0.1 is the usual rounding applied to magnitudes, is due to a “true” source b value, b, and use these likelihoods to estimate source b ranges corresponding to various confidence levels. As an example of application of the method, we estimate the b values before and after the occurrence of a 7.4-magnitude earthquake in the Mexican subduction zone, and find a difference of 0.82 between them with 100% confidence that the b values are different.     

<Emphasis Type="Italic">Q</Emphasis><Subscript>c</Subscript> and <Emphasis Type="Italic">Q</Emphasis><Subscript>S</Subscript> wave attenuation of South African earthquakes

Quality factor Q, which describes the attenuation of seismic waves with distance, was determined for South Africa using data recorded by the South African National Seismograph Network. Because of an objective paucity of seismicity in South Africa and modernisation of the seismograph network only in 2007, I carried out a coda wave decay analysis on only 13 tectonic earthquakes and 7 mine-related events for the magnitude range 3.6?≤?M _L ?≤?4.4. Up to five seismograph stations were utilised to determine Q _c for frequencies at 2, 4, 8 and 16 Hz resulting in 84 individual measurements. The constants Q ₀ and α were determined for the attenuation relation Q _c(f)?=?Q ₀ f ^α . The result was Q ₀?=?396?±?29 and α?=?0.72?±?0.04 for a lapse time of 1.9*(t _s???t ₀) (time from origin time t ₀ to the start of coda analysis window is 1.9 times the S-travel time, t _s) and a coda window length of 80 s. This lapse time and coda window length were found to fit the most individual frequencies for a signal-to-noise ratio of at least 3 and a minimum absolute correlation coefficient for the envelope of 0.5. For a positive correlation coefficient, the envelope amplitude increases with time and Q _c was not calculated. The derived Q _c was verified using the spectral ratio method on a smaller data set consisting of nine earthquakes and one mine-related event recorded by up to four seismograph stations. Since the spectral ratio method requires absolute amplitudes in its calculations, site response tests were performed to select four appropriate stations without soil amplification and/or signal distortion. The result obtained for Q _S was Q ₀?=?391?±?130 and α?=?0.60?±?0.16, which agrees well with the coda Q _c result.     

A new <Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript> estimation parameter for use in earthquake early warning systems

We propose a method that employs the squared displacement integral (ID2) to estimate earthquake magnitudes in real time for use in earthquake early warning (EEW) systems. Moreover, using τ _c and P _d for comparison, we establish formulas for estimating the moment magnitudes of these three parameters based on the selected aftershocks (4.0 ≤ M _s  ≤ 6.5) of the 2008 Wenchuan earthquake. In this comparison, the proposed ID2 method displays the highest accuracy. Furthermore, we investigate the applicability of the initial parameters to large earthquakes by estimating the magnitude of the Wenchuan M _s 8.0 mainshock using a 3-s time window. Although these three parameters all display problems with saturation, the proposed ID2 parameter is relatively accurate. The evolutionary estimation of ID2 as a function of the time window shows that the estimation equation established with ID2 ^Ref determined from the first 8-s of P wave data can be directly applicable to predicate the magnitudes of 8.0. Therefore, the proposed ID2 parameter provides a robust estimator of earthquake moment magnitudes and can be used for EEW purposes.     

Precursor times of abnormal <Emphasis Type="Italic">b</Emphasis>-values prior to mainshocks

Seismic observations exhibit the presence of abnormal b-values prior to numerous earthquakes. The time interval from the appearance of abnormal b-values to the occurrence of mainshock is called the precursor time. There are two kinds of precursor times in use: the first one denoted by T is the time interval from the moment when the b-value starts to increase from the normal one to the abnormal one to the occurrence time of the forthcoming mainshock, and the second one denoted by T _p is the time interval from the moment when the abnormal b-value reaches the peak one to the occurrence time of the forthcoming mainshock. Let T* be the waiting time from the moment when the abnormal b-value returned to the normal one to the occurrence time of the forthcoming mainshock. The precursor time, T (usually in days), has been found to be related to the magnitude, M, of the mainshock expected in a linear form as log(T)?=?q?+?rM where q and r are the coefficient and slope, respectively. In this study, the values of T, T _p , and T* of 45 earthquakes with 3?≤?M?≤?9 occurred in various tectonic regions are compiled from or measured from the temporal variations in b-values given in numerous source materials. The relationships of T and T _p , respectively, versus M are inferred from compiled data. The difference between the values of T and T _p decreases with increasing M. In addition, the plots of T*/T versus M, T* versus T, and T* versus T-T* will be made and related equations between two quantities will be inferred from given data.     

Resources and environmental systems management under synchronic interval uncertainties

Resources and environmental systems management (RESM) is challenged by the synchronic effects of interval uncertainties in the related practices. The synchronic interval uncertainties are misrepresented as random variables, fuzzy sets, or interval numbers in conventional RESM programming techniques including stochastic programming. This may lead to ineffectiveness of resources allocation, high costs of recourse measures, increased risks of unreasonable decisions, and decreased optimality of system profits. To fill the gap of few corresponding studies, a synchronic interval linear programming (SILP) method is proposed in this study. The proposition of interval sets and interval functions and coupling them with linear programming models lead to development of an SILP model for RESM. This enables incorporation of interval uncertainties in resource constraints and synchronic interval uncertainties in the programming objective into the optimization process. An analysis of the distribution-independent geometric properties of the feasible regions of SILP models results in proposition of constraint violation likelihoods. The tradeoff between system optimality and constraint violation is analyzed. The overall optimality of SILP systems under synchronic intervalness is quantified through proposition of integrally optimal solutions. Integration of these efforts leads to a violation-constrained interval integral method for optimization of RESM systems under synchronic interval uncertainties. Comparisons with selected existing methods reveal the effectiveness of SILP at eliminating negativity of synchronic intervalness, enabling risk management of and achieving overall optimality of RESM systems, and enhancing the reliability of optimization techniques for RESM problems. The exploited framework for analyzing synchronic interval uncertainties in RESM systems is helpful for addressing synchronisms of other uncertainties such as randomness or fuzziness and avoiding the resultant decision mistakes and disasters due to neglecting them.     

Response of Pinus taiwanensis growth to climate changes at its southern limit of Daiyun Mountain,mainland China Fujian Province

The Pinus taiwanensis trees of a Chinese endemic species form pure forests at infertile sites in humid subtropical China, which can aid soil reservation and ecological resiliency at such sites. Herein, we used dendrochronological methods to investigate the growth patterns and their relationship with climate by analyzing 158 cores of 79 P. taiwanensis trees at 4 sites in their southernmost distribution in mainland China at Daiyun Mountain in Fujian Province. Tree growths at sites with favorable conditions generally showed an age-related growth trend, which decreased from approximately 5 mm to 0.5 mm in about a century. Trees in moderately stressed environments established high growth in their juvenile periods but were highly susceptible to environmental stresses such as a sharp growth decline in the 1990s. The temperature in February and the moisture in July are the major limiting factors for most of the tree growths, except for a few extremely stressed P. taiwanensis trees. The growth of the dwarf P. taiwanensis trees, with all ring diameters of the first 20 years less than 4 mm, is mainly limited by the dry climate in May. Under the potential future warming trend, drought stress can be particularly threatening to these dwarf P. taiwanensis trees, which are likely to be the germplasm resources for this species in this region.     

Spatial inhomogeneity of variations in the <Emphasis Type="Italic">F</Emphasis><Subscript>2</Subscript> layer critical frequencies under the conditions of winter 1995

The daily samples of the hourly measurements of the foF2 critical frequency, obtained on January 5–21, 1995, at the midlatitude and high-latitude automated ionospheric stations (geographic latitude higher than 60°), are considered. The {fo} sets are transformed into the {δfoF2} sets of relative variations, for which asymmetry (A) and excess (E) are calculated. The selected stations are grouped into 20 pairs of automated ionospheric stations (AISs) located at distances of 200–10 000 km from one another. Sign estimates of the cross-correlation coefficients between the sets of 16 A and E values for different pairs of stations are used. Two types of structures of the statistical invariant spatial distribution are established: the structures with a scale of about 300 km, invariant with respect to latitude, and with a scale of about 6000 km (for only high latitudes).     

Influence of aquatic plants on the hydrogen isotope composition of sedimentary long-chain n-alkanes in the Lake Qinghai region,Qinghai-Tibet Plateau

The hydrogen isotopic composition(δD) of leaf wax long-chain n-alkanes(C27, C29, and C31) from lacustrine sediments has been widely applied to reconstruct terrestrial paleoclimatic and paleohydrological changes. However, few studies have addressed whether the aquatic-derived n-alkanes can affect the δD values of lake sedimentary long-chain n-alkanes, which are usually regarded as a recorder of the terrestrial hydrological signals. Here we systematically investigated δD values of long-chain n-alkanes from modern aquatic plants, both near-shore and off-shore surface sediments, surrounding terrestrial plant litters, as well as river water and lake water in Lake Qinghai and its satellite lakes on the northeastern Qinghai-Tibet Plateau. Our data showed that(i) δD values of long-chain n-alkanes from aquatic plants varied from-184‰ to-132‰ for n-C27, from-183‰ to-138‰ for n-C29, and from-189‰ to-130‰ for n-C31, respectively, with no significant differences among the three n-alkanes homologues;(ii) δD values of long-chain n-alkanes from aquatic plants were generally more positive than those from surrounding terrestrial plants, possibly because that they recorded the D-enrichment of lake water in this semi-arid region;(iii) δD values of long-chain n-alkanes from surface sediments showed significant differences among the three n-alkanes homologues, due to the larger aquatic input of n-C27 to the sedimentary lipid pool than that of n-C31, and(iv) n-C27 δD values of near-shore aquatic plants and near-shore sediments are more negative than those from off-shore as a result of lower δD values of near-shore lake water. Our findings indicate that in this region(i) the offset between sedimentary n-C27 and n-C31 δD values(ΔδDC27-C31) could potentially be used to evaluate if sedimentary long-chain n-alkanes are derived from a single source;(ii) while δD values of n-C27 may be influenced by lake water hydrological changes, sedimentary n-C31 is derived predominantly from terrestrial plants and thus its δD can serve as a relatively reliable indicator for terrestrial paleoclimatic and paleohydrological reconstructions.     

Experimental examination on the heterogeneity parameter <Emphasis Type="Italic">C</Emphasis><Subscript>v</Subscript> of earthquake precursors

Two rock samples with different structures and materials were deformed under a biaxial loading system, and multipoint strain measurements were performed for each sample. The distribution of strain anomalies during the deformation and the instability process were analyzed by using C _v value put forward by WANG Xiao-qing and CHEN Xue-zhong, et al, a parameter to describe the heterogeneous distribution of earthquake precursors, so as to examine the method of C _v value and to explore its physical meaning experimentally. The result shows that the change of C _v value is correlated to the change of deformation characteristics and is an effective parameter to describe the heterogeneity of precursor distribution. C _v value increases firstly and then decreases before the instability, and the instability occurs when C _v value decreases to the level before increasing. This indicates that C _v value may be a useful parameter for earthquake prediction.     

Frictional and viscous effects on the nucleation phase of an earthquake

We study the frictional and viscous effects on earthquake nucleation, especially for the nucleation phase, based on a one-degree-of-freedom spring-slider model with friction and viscosity. The frictional and viscous effects are specified by the characteristic displacement, U _c, and viscosity coefficient, η, respectively. Simulation results show that friction and viscosity can both lengthen the natural period of the system and viscosity increases the duration time of motion of the slider. Higher viscosity causes a smaller amplitude of lower velocity motion than lower viscosity. A change of either U _c (under large η) or η (under large U _c) from a large value (U _ch for U _c and η _h for η) to a small one (U _cl for U _c and η _l for η) in two stages during sliding can result in a clear nucleation phase prior to the P-wave. The differences δU _c = U _ch ? U _cl and δη = η _h ? η _l are two important factors in producing a nucleation phase. The difference between the nucleation phase and the P-wave increases with either δU _c or δη. Like seismic observations, the peak amplitude of P-wave, which is associated with the earthquake magnitude, is independent upon the duration time of nucleation phase. A mechanism specified with a change of either η or U _c from a larger value to a smaller one due to temporal variations in pore fluid pressure and temperature in the fault zone based on radiation efficiency is proposed to explain the simulation results and observations.     

Statistics for sample splitting for the calibration and validation of hydrological models

Hydrological models have been widely applied in flood forecasting, water resource management and other environmental sciences. Most hydrological models calibrate and validate parameters with available records. However, the first step of hydrological simulation is always to quantitatively and objectively split samples for use in calibration and validation. In this paper, we have proposed a framework to address this issue through a combination of a hierarchical scheme through trial and error method, for systematic testing of hydrological models, and hypothesis testing to check the statistical significance of goodness-of-fit indices. That is, the framework evaluates the performance of a hydrological model using sample splitting for calibration and validation, and assesses the statistical significance of the Nash–Sutcliffe efficiency index (E_f), which is commonly used to assess the performance of hydrological models. The sample splitting scheme used is judged as acceptable if the E_f values exceed the threshold of hypothesis testing. According to the requirements of the hierarchical scheme for systematic testing of hydrological models, cross calibration and validation will help to increase the reliability of the splitting scheme, and reduce the effective range of sample sizes for both calibration and validation. It is illustrated that the threshold of E_f is dependent on the significance level, evaluation criteria (both regarded as the population), distribution type, and sample size. The performance rating of E_f is largely dependent on the evaluation criteria. Three types of distributions, which are based on an approximately standard normal distribution, a Chi square distribution, and a bootstrap method, are used to investigate their effects on the thresholds, with two commonly used significance levels. The highest threshold is from the bootstrap method, the middle one is from the approximately standard normal distribution, and the lowest is from the Chi square distribution. It was found that the smaller the sample size, the higher the threshold values are. Sample splitting was improved by providing more records. In addition, outliers with a large bias between the simulation and the observation can affect the sample values of E_f, and hence the output of the sample splitting scheme. Physical hydrology processes and the purpose of the model should be carefully considered when assessing outliers. The proposed framework in this paper cannot guarantee the best splitting scheme, but the results show the necessary conditions for splitting schemes to calibrate and validate hydrological models from a statistical point of view.