Spatiotemporal variations of actual evapotranspiration over the Lake Selin Co and surrounding small lakes (Tibetan Plateau) during 2003–2012

Actual evapotranspiration(ET_a) over the Tibetan Plateau(TP) is an important component of the water cycle,and greatly influences the water budgets of the TP lake basins.Quantitative estimation of ET_a within lake basins is fundamental to physically understanding ET_a changes,and thus will improve the understanding of the hydro logical processes and energy balance throughout the lake basins.In this study,the spatiotemporal dynamic changes of ET_α within the Lake Selin Co(the TP's largest lake) and its surrounding small lakes and land area during 2003-2012 are examined at the basin scale.This was carried out using the well-established Water and Energy Budget-based Distributed Hydrological Model(WEB-DHM) for the land area,the Penman method for the water area when unfrozen,and a simple sublimation estimation approach for the water area when frozen.The relationships between ET_a changes and controlling factors are also discussed.Results indicate that the simulated land ET_a from the WEB-DHM reasonably agrees with the estimated ET_a values from the nonlinear complementary relationship model using appropriately calibrated parameter values at a point scale.Land ET_a displayed a non-significant increase of 7.03 mm year~(-1),and largely depends on precipitation.For the water area,the combined effects of reduced wind speed and net radiation offset the effect of rising temperature and vapor pressure deficit,and contributed to a non-significant decrease in evaporation of 4.17 mm year~(-1).Sensitivity analysis shows that vapor pressure deficit and wind speed are the most sensitive variables to the changes of evaporation from the water area.     

Altitudinal effect of soil <Emphasis Type="Italic">n</Emphasis>-alkane <Emphasis Type="Italic">δ</Emphasis>D values on the eastern Tibetan Plateau and their increasing isotopic fractionation with altitude

Stable isotope paleoaltimetry has provided unprecedented insights into the topographic histories of many of the world’s highest mountain ranges. However, on the Tibetan Plateau (TP), stable isotopes from paleosols generally yield much higher paleoaltitudes than those based on fossils. It is therefore essential when attempting to interpret accurately this region’s paleoaltitudes that the empirical calibrations of local stable isotopes and the relations between them are established. Additionally, it is vital that careful estimations be made when estimate how different isotopes sourced from different areas may have been influenced by different controls. We present here 29 hydrogen isotopic values for leaf wax-derived n-alkanes (i.e., _δD_wax values, and abundance-weighted average δD values of C₂₉ and C₃₁) in surface soils, as well as the δD values of soil water (δD_sw) samples (totaling 22) from Mount Longmen (LM), on the eastern TP (altitude ~0.8–4.0 km above sea level (asl), a region climatically affected by the East Asian Monsoon (EAM). We compared our results with published data from Mount Gongga (GG). In addition, 47 river water samples, 55 spring water samples, and the daily and monthly summer precipitation records (from May to October, 2015) from two precipitation observation stations were collected along the GG transect for δD analysis. LM soil _δD_wax values showed regional differences and responded strongly to altitude, varying from?160‰ to?219‰, with an altitudinal lapse rate (ALR) of?18‰ km^?1 (R ²=0.83; p<0.0001; n=29). These _δD_wax values appeared more enriched than those from the GG transect by ~40‰. We found that both the climate and moisture sources led to the differences observed in soil _δD_wax values between the LM and GG transects. We found that, as a general rule, ε _wax/rw, ε _wax/p and ε _wax/sw values (i.e., the isotopic fractionation of _δD_wax corresponding to δD_rw, δD_p and δD_sw) increased with increasing altitude along both the LM and GG transects (up to 34‰and 50‰, respectively). Basing its research on a comparative study of _δD_wax, δD_p, δD_rw(δD_springw) and δD_sw, this paper discusses the effects of moisture recycling, glacier-fed meltwater, relative humidity (RH), evapotranspiration (ET), vegetation cover, latitude, topography and/or other factors on ε _wax/p values. Clearly, if ε _wax-p values at higher altitudes are calculated using smaller ε _wax-p values from lower altitudes, the calculated paleowaterδD_p values are going to be more depleted than the actual δD values, and any paleoaltitude would therefore be overestimated.     

Evapotranspiration from citrus trees growing in sandy soilunder drip irrigation with saline water

An experiment on evapotranspiration from citrus trees under irrigation with saline waterwas carried out for 4 months. Two lysimeters planted with a citrus tree in the green house wereused. One lysimeter was irrigated with saline water (NaCl and CaCl2 of 2000 mg/L equivalence,EC = 3.8 dS/m, SAR = 5.9) and the other was irrigated with freshwater using drip irrigation. Theapplied irrigation water was 1.2 times that of the evapotranspiration on the previous day.Evapotranspiration was calculated as the change in lysimeter weight recorded every 30 minutes.The lysimeters were filled with soil with 95.8% sand. The results of the experiment were as follows.(i) The evapotranspiration from citrus tree was reduced after irrigation with saline water. Theevapotranspiration returns to normal after leaching. However it takes months to exhaust the saltfrom the tree. ( ii ) To estimate the impact of irrigation with saline water on the evapotranspirationfrom citrus trees, the reduction coefficient due to salt stress (Ks) was used in this experiment.Evapotranspiration under irrigation with saline water (ETs) can be calculated from evapotranspira-tion under irrigation with freshwater (ET) by the equation ETs = Ks× ET. Ks can be expressed as afunction of ECsw. (iii) The critical soil-water electrical conductivity (ECsw) is 9.5 dS/m, beyondwhich adverse effects on evapotranspiration begin to appear. If ECsw can be controlled at below9.5 dS/m, saline water can be safely used for irrigation.     

Comparison between gene expression programming and traditional models for estimating evapotranspiration under hyper arid Conditions

Gene Expression Programming (GEP) was used to develop new mathematical equations for estimating daily reference evapotranspiration (ET _ref) for the Kingdom of Saudi Arabia. The daily climatic variables were collected by 13 meteorological stations from 1980 to 2010. The GEP models were trained on 65% of the climatic data and tested using the remaining 35%. The generalised Penman-Monteith model was used as a reference target for evapotranspiration (ET) values, with h _c varies from 5 to 105 cm with increment of a centimetre. Eight GEP models have been compared with four locally calibrated traditional models (Hargreaves-Samani, Irmak, Jensen-Haise and Kimberly-Penman). The results showed that the statistical performance criteria values such as determination coefficients (R ²) ranged from as low as 64.4% for GEP-MOD1, where the only parameters included (maximum, minimum, and mean temperature and crop height), to as high as 95.5% for GEP-MOD8 with which all climatic parameters included (maximum, minimum and mean temperature; maximum, minimum and mean humidity; solar radiation; wind speed; and crop height). Moreover, an interesting founded result is that the solar radiation has almost no effect on ET _ref under the hyper arid conditions. In contrast, the wind speed and plant height have a great positive impact in increasing the accuracy of calculating ET _ref. Furthermore, eight GEP models have obtained better results than the locally calibrated traditional ET _ref equations.     

On Creating Global Gridded Terrestrial Water Budget Estimates from Satellite Remote Sensing

The increasing availability and reliability of satellite remote sensing products [e.g., precipitation (P), evapotranspiration (ET), and the total water storage change (TWSC)] make it feasible to estimate the global terrestrial water budget at fine spatial resolution. In this study, we start from a reference water budget dataset that combines all available data sources, including satellite remote sensing, land surface model (LSM) and reanalysis, and investigate the roles of different non-satellite remote sensing products in closing the terrestrial water budget through a sensitivity analysis by removing/replacing one or more categories of products during the budget estimation. We also study the differences made by various satellite products for the same budget variable. We find that the gradual removal of non-satellite data sources will generally worsen the closure errors in the budget estimates, and remote sensing retrievals of P, ET, and TWSC together with runoff (R) from LSM give the worst closure errors. The gauge-corrected satellite precipitation helps to improve the budget closure (4.2–9 % non-closure errors of annual mean precipitation) against using the non-gauge-corrected precipitation (7.6–10.4 % non-closure errors). At last, a data assimilation technique, the constrained Kalman filter, is applied to enforce the water balance, and it is found that the satellite remote sensing products, though with worst closure, yield comparable budget estimates in the constrained system to the reference data. Overall, this study provides a first comparison between the water budget closure using the satellite remote sensing products and a full combination of remote sensing, LSM, and reanalysis products on a quasi-global basis. This study showcases the capability and potential of the satellite remote sensing in closing the terrestrial water budget at fine spatial resolution if properly constrained.     

Crustal Structure Along Sunda-Banda Arc Transition Zone from Teleseismic Receiver Functions

We analyzed receiver function of teleseismic events recorded at twelve Indonesian-GEOFON (IA-GE) broadband stations using nonlinear Neighbourhood Algorithm (NA) inversion and H-k stacking methods to estimate crustal thickness, V _p /V _s ratios and S-wave velocity structure along Sunda-Banda arc transition zone. We observed crustal thickness of 34–37 km in Timor Island, which is consistent with the previous works. The thick crust (> 30 km) is also found beneath Sumba and Flores Islands, which might be related to the arc-continent collision causing the thickened crust. In Timor and Sumba Islands, we observed high V _p /V _s ratio (> 1.84) with low velocity zone that might be associated with the presence of mafic and ultramafic materials and fluid filled fracture zone. The high V _p /V _s ratio observed at Sumbawa and Flores volcanic Islands might be an indication of partial melt related to the upwelling of hot asthenosphere material through the subducted slab.     

Photosynthetic characteristics of dominant tree species and canopy in the broadleaved Korean pine forest of Changbai Mountains

Based on the light-photosynthesis response measurement at leaf level, combined with over-and under-canopy eddy covariance measurements, research on photosynthetic characteristics of single trees and forest canopy was conducted. The relationship between light intensity and photosynthetic rates for leaves and canopy can be well fitted by a non-rectangular hyperbola model. Mongolian oak presented a high light compensation point, L _cp (28 μmol·m^?2·s^?1), a light saturation point L _sp (>1800 μmol·m^?2·s^?1), and a maximal net photosynthetic rate P _max (9.96 μmol·m^?2·s^?1), which suggest that it is a typical heliophilous plant. Mono maple presented the highest apparent quantum efficiency α (0.066) but the lowest, L _cp (16 μmol·m^?2·s^?1), L _sp (≈800 μmol·m^?2·s^?1), and P _max (4.51 μmol·m^?2·s^?1), which suggest that it is heliophilous plant. Korean pine showed the lowest α value but a higher P _max, which suggest that it is a semi-heliophilous plant. At the canopy level, the values of both α and P _max approached the upper limit of reported values in temperate forests, while L _cp was within the lower limit. Canopy photosynthetic characteristics were well consistent with those of leaves. Both showed a high ability to photosynthesize. However, environmental stresses, especially high vapor pressure deficits, could significantly reduce the photosynthetic ability of leaves and canopy.     

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.     

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.     

Shear wave velocity-based evaluation and design of stone column improved ground for liquefaction mitigation

The evaluation and design of stone column improvement ground for liquefaction mitigation is a challenging issue for the state of practice. In this paper, a shear wave velocity-based approach is proposed based on the well-defined correlations of liquefaction resistance (CRR)-shear wave velocity (V _s)-void ratio (e) of sandy soils, and the values of parameters in this approach are recommended for preliminary design purpose when site specific values are not available. The detailed procedures of pre- and post-improvement liquefaction evaluations and stone column design are given. According to this approach, the required level of ground improvement will be met once the target V _s of soil is raised high enough (i.e., no less than the critical velocity) to resist the given earthquake loading according to the CRR-V _s relationship, and then this requirement is transferred to the control of target void ratio (i.e., the critical e) according to the V _s-e relationship. As this approach relies on the densification of the surrounding soil instead of the whole improved ground and is conservative by nature, specific considerations of the densification mechanism and effect are given, and the effects of drainage and reinforcement of stone columns are also discussed. A case study of a thermal power plant in Indonesia is introduced, where the effectiveness of stone column improved ground was evaluated by the proposed V _s-based method and compared with the SPT-based evaluation. This improved ground performed well and experienced no liquefaction during subsequent strong earthquakes.     

On the impact of wind on the development of wave field during storm Britta

The observation of extreme waves at FINO 1 during storm Britta on the 1st November 2006 has initiated a series of research studies regarding the mechanisms behind. The roles of stability and the presence of the open cell structures have been previously investigated but not conclusive. To improve our understanding of these processes, which are essential for a good forecast of similarly important events offshore, this study revisits the development of storm Britta using an atmospheric and wave coupled modeling system, wind and wave measurements from ten stations across the North Sea, cloud images and Synthetic Aperture Radar (SAR) data. It is found here that a standard state-of-the-art model is capable of capturing the important characteristics of a major storm like Britta, including the storm path, storm peak wind speed, the open cells, and peak significant wave height (H _s ) for open sea. It was also demonstrated that the impact of the open cells has negligible contribution to the development of extreme H _s observed at FINO 1. At the same time, stability alone is not sufficient in explaining the development of extreme H _s . The controlling conditions for the development of Britta extreme H _s observed at FINO 1 are the persistent strong winds and a long and undisturbed fetch over a long period.     

Estimation of Seismic Site Coefficient and Seismic Microzonation of Imphal City,India, Using the Probabilistic Approach

Seismic site coefficients (F _s ) for Imphal city have been estimated based on 700 synthetically generated earthquake time histories through stochastic finite fault method, considering various combinations of magnitudes and fault distances that may affect Imphal city. Seismic hazard curves and Uniform Hazard Response Spectra (UHRS) are presented for Imphal city. F _s have been estimated based on site response analyses through SHAKE-91 for a period range of engineering interest (PGA to 3.0 s), for 5% damping. F _s were multiplied by UHRS values to obtain surface level spectral acceleration with 2 and 10% probability of exceedance in 50 year (～2500 and ～500 year) return period. Comparison between predicted mean surface level response spectra and IS-1893 code shows that spectral acceleration value is higher for longer periods (i.e., >1.0 s), for ～500 year return period, and lower for periods shorter than 0.2 s for ～2500 year return period.     

Homologous temperature of olivine: Implications for creep of the upper mantle and fabric transitions in olivine

The homologues temperature of a crystalline material is defined as T/T_m, where T is temperature and T_m is the melting(solidus) temperature in Kelvin. It has been widely used to compare the creep strength of crystalline materials. The melting temperature of olivine system,(Mg,Fe)_2SiO_4, decreases with increasing iron content and water content, and increases with confining pressure. At high pressure, phase transition will lead to a sharp change in the melting curve of olivine. After calibrating previous melting experiments on fayalite(Fe_2SiO_4), the triple point of fayalite-Fe_2SiO_4 spinel-liquid is determined to be at 6.4 GPa and 1793 K. Using the generalized means, the solidus and liquidus of dry olivine are described as a function of iron content and pressure up to 6.4 GPa. The change of T/T_m of olivine with depth allows us to compare the strength of the upper mantle with different thermal states and olivine composition. The transition from semi-brittle to ductile deformation in the upper mantle occurs at a depth where T/T_m of olivine equals 0.5. The lithospheric mantle beneath cratons shows much smaller T/T_m of olivine than orogens and extensional basins until the lithosphere-asthenosphere boundary where T/T_m 0.66, suggesting a stronger lithosphere beneath cratons. In addition, T/T_m is used to analyze deformation experiments on olivine. The results indicate that the effect of water on fabric transitions in olivine is closely related with pressure. The hydrogen-weakening effect and its relationship with T/T_m of olivine need further investigation. Below 6.4 GPa(200 km), T/T_m of olivine controls the transition of dislocation glide from [100] slip to [001] slip. Under the strain rate of 10~(-12)–10~(-15) s~(-1) and low stress in the upper mantle, the [100](010) slip system(A-type fabric) becomes dominant when T/T_m 0.55–0.60. When T/T_m 0.55–0.60, [001] slip is easier and low T/T_m favors the operation of [001](100) slip system(C-type fabric). This is consistent with the widely observed A-type olivine fabric in naturally deformed peridotites, and the C-type olivine fabric in peridotites that experienced deep subduction in ultrahigh-pressure metamorphic terranes. However, the B-type fabric will develop under high stress and relatively low T/T_m. Therefore, the homologues temperature of olivine established a bridge to extrapolate deformation experiments to rheology of the upper mantle. Seismic anisotropy of the upper mantle beneath cratons should be simulated using a four-layer model with the relic A-type fabric in the upper lithospheric mantle, the B-type fabric in the middle layer, the newly formed A- or B-type fabric near the lithosphere-asthenosphere boundary, and the asthenosphere dominated by diffusion creep below the Lehmann discontinuity. Knowledge about transition mechanisms of olivine fabrics is critical for tracing the water distribution and mantle flow from seismic anisotropy.     

A way of estimating the characteristic slip displacement

During the ruptures of an earthquake,the strain energy.△E,.will be transferred into,at least,three parts,i.e..the seismic radiation energy(E_s),fracture energy(E_g),and frictional energy(E_f),that is,△E = E_s + E_g + E_f.Friction,which is represented by a velocity- and state-dependent friction law by some researchers,controls the three parts.One of the main parameters of the law is the characteristic slip displacement.D_c.It is significant and necessary to evaluate the reliable value of D_c from observed and inverted seismic data.Since D_c controls the radiation efficiency.η_R = E_s/(E_s+ E_g),the value of η_r is a good constraint of estimating D_c.Integrating observed data and inverted results of source parameters from recorded seismograms.the values of E_s and E_g of an earthquake can be measured,thus leading to the value of η_R.The constraint used to estimate the reliable value of D_c will be described in this work.An example of estimates of D_c.based on the observed and inverted values of source parameters of the September 20,1999 M_S 7.6 Chi-Chi(Ji-Ji).Taiwan region,earthquake will be presented.     

G: Fracture energy,friction and dissipation in earthquakes

Recent estimates of fracture energy G ^′ in earthquakes show a power-law dependence with slip u which can be summarized as G ^′ ∝ u ^a where a is a positive real slightly larger than one. For cracks with sliding friction, fracture energy can be equated to G _f : the post-failure integral of the dynamic weakening curve. If the dominant dissipative process in earthquakes is friction, G ^′ and G _f should be comparable and show a similar scaling with slip. We test this hypothesis by analyzing experiments performed on various cohesive and non-cohesive rock types, under wet and dry conditions, with imposed deformation typical of seismic slip (normal stress of tens of MPa, target slip velocity > 1 m/s and fast accelerations ≈ 6.5 m/s²). The resulting fracture energy G _f is similar to the seismological estimates, with G _f and G ^′ being comparable over most of the slip range. However, G _f appears to saturate after several meters of slip, while in most of the reported earthquake sequences, G ^′ appears to increase further and surpasses G _f at large magnitudes. We analyze several possible causes of such discrepancy, in particular, additional off-fault damage in large natural earthquakes.     

<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.     

Seismotectonic studies of the pleistoseist area of the <Emphasis Type="Italic">M</Emphasis> <Subscript> <Emphasis Type="Italic">s</Emphasis> </Subscript> = 6.9 Ilin-Tass earthquake in northeast Yakutia

The complex seismotectonic studies of the pleistoseist area of the Ilin-Tas earthquake (M_s = 6.9), one of the strongest seismic events ever recorded by the regional seismic network in northeastern Russia, are carried out. The structural tectonic position, morphotectonic features of present-day topography, active faults, and types of Cenozoic deformations of the epicentral zone are analyzed. The data of the instrumental observations are summarized, and the manifestations of the strong seismic events in the Yana–Indigirka segment of the Cherskii seismotectonic zone are considered. The explanation is suggested for the dynamical tectonic setting responsible for the Andrei-Tas seismic maximum. This setting is created by the influence of the Kolyma–Omolon indenter, which intrudes into the Cherskii seismotectonic zone from the region of the North American lithospheric plate and forms the main seismogenic structures of the Yana–Indigirka segment in the frontal zone (the Ilin-Tas anticlinorium). The highest seismic potential is noted in the Andrei- Tas block—the focus of the main tectonic impacts from the Kolyma–Omolon superterrane. The general trend of this block coincides with the orientation of the major axis of isoseismal ellipses (azimuth 50°–85°), which were determined from the observations of macroseismic effects on the ground after the Uyandina (M_s = 5.6), Andrei-Tas (M_s = 6.1), and Ilin-Tas (M_s = 6.9) earthquakes.     

Experimental Seismic Event-screening Criteria at the Prototype International Data Center

—?Experimental seismic event-screening capabilities are described, based on the difference of body-and surface-wave magnitudes (denoted as M _s :m _b ) and event depth. These capabilities have been implemented and tested at the prototype International Data Center (PIDC), based on recommendations by the IDC Technical Experts on Event Screening in June 1998. Screening scores are presented that indicate numerically the degree to which an event meets, or does not meet, the M _s :m _b and depth screening criteria. Seismic events are also categorized as onshore, offshore, or mixed, based on their 90% location error ellipses and an onshore/offshore grid with five-minute resolution, although this analysis is not used at this time to screen out events.¶Results are presented of applications to almost 42,000 events with m _b ?≥?3.5 in the PIDC Standard Event Bulletin (SEB) and to 121 underground nuclear explosions (UNE's) at the U.S. Nevada Test Site (NTS), the Semipalatinsk and Novaya Zemlya test sites in the Former Soviet Union, the Lop Nor test site in China, and the Indian, Pakistan, and French Polynesian test sites. The screening criteria appear to be quite conservative. None of the known UNE's are screened out, while about 41 percent of the presumed earthquakes in the SEB with m _b ?≥?3.5 are screened out. UNE's at the Lop Nor, Indian, and Pakistan test sites on 8 June 1996, 11 May 1998, and 28 May 1998, respectively, have among the lowest M _s :m _b scores of all events in the SEB.¶To assess the validity of the depth screening results, comparisons are presented of SEB depth solutions to those in other bulletins that are presumed to be reliable and independent. Using over 1600 events, the comparisons indicate that the SEB depth confidence intervals are consistent with or shallower than over 99.8 percent of the corresponding depth estimates in the other bulletins. Concluding remarks are provided regarding the performance of the experimental event-screening criteria, and plans for future improvements, based on recent recommendations by the IDC Technical Experts on Event Screening in May 1999.     

Experimental study on the shear behavior of the interface between cushion materials and the concrete raft

Cushion is a layer of granular materials between the raft and the ground. The shear behavior of the interface between the cushion and the raft may influence the seismic performance of the superstructure. In order to quantify such influences, horizontal shear tests on the interfaces between different cushion materials and concrete raft under monotonic and cyclic loading were carried out. The vertical pressure P_v, material type and cushion thickness h_c were taken as variables. Conclusions include: 1) under monotonic loading, P_v is the most significant factor; the shear resistance P_(hmax) increases as P_v increases, but the normalized factor of resistance μ_n has an opposite tendency; 2) for the materials used in this study, μ_n varies from 0.40 to 0.70, the interface friction angle δ_s varies from 20° to 35°, while u_(max) varies from 3 mm to 15 mm; 3) under cyclic loading, the interface behavior can be abstracted as a "three-segment" back-bone curve, the main parameters include μ_n, the displacement u_1 and stiffness K_1 of the elastic stage, the displacement u_2 and stiffness K_2 of the plastic stage; 4) by observation and statistical analysis, the significance of different factors, together with values of K_1, K_2 and μ_n have been obtained.     

Effects of temporally correlated infiltration on water flow in an unsaturated–saturated system

Effects of temporally correlated infiltration on water flow in an unsaturated–saturated system were investigated. Both white noise and exponentially correlated infiltration processes were considered. The moment equations of the pressure head (ψ) were solved numerically to obtain the variance and autocorrelation functions of ψ at 14 observation points. Monte Carlo simulations were conducted to verify the numerical results and to estimate the power spectrum of ψ (S _ψψ ). It was found that as the water flows through the system, the variance of the ψ (\( \sigma_{\psi }^{2} \)) were damped by the system: the deeper in the system, the smaller the \( \sigma_{\psi }^{2} \), and the larger the correlation timescale of the infiltration process (λ _I ), the larger the \( \sigma_{\psi }^{2} \). The unsaturated–saturated system gradually filters out the short-term fluctuations of ψ and the damping effect is most significant in the upper part of the system. The fluctuations of ψ is non-stationary at early time and becomes stationary as time progresses: the larger the value of λ _I , the longer the non-stationary period. The correlation timescale of the ψ (λ _ψ ) increases with depth and approaches a constant value at depth: the larger the value of λ _I , the larger the value of λ _ψ . The results of the estimated S _ψψ is consistent with those of the variance and autocorrelation function.