Noise reduction and detection of weak, coherent signals through phase-weighted stacks

We present a new tool for efficient incoherent noise reduction for array data employing complex trace analysis. An amplitude-unbiased coherency measure is designed based on the instantaneous phase, which is used to weight the samples of an ordinary, linear stack. The result is called the phase-weighted stack (PWS) and is cleaned from incoherent noise. PWS thus permits detection of weak but coherent arrivals. The method presented can easily be extended to phase-weighted cross-correlations or be applied in the τ p domain. We illustrate and discuss the advantages and disadvantages of PWS in comparison with other coherency measures and present examples. We further show that our non-linear stacking technique enables us to detect a weak lower-mantle P -to- S conversion from a depth of approximately 840 km on array data. Hints of an 840 km discontinuity have been reported; however, such a discontinuity is not yet established due to the lack of further evidence.     

High-resolution structures of the Landers fault zone inferred from aftershock waveform data

High-frequency body waves recorded by a temporary seismic array across the surface rupture trace of the 1992 Landers, California, earthquake were used to determine fault-zone structures down to the seismogenic depth. We first developed a technique to use generalized ray theory to compute synthetic seismograms for arbitrarily oriented tabular low-velocity fault-zone models. We then generated synthetic waveform record sections of a linear array across a vertical fault zone. They show that both arrival times and waveforms of P and S waves vary systematically across the fault due to transmissions and reflections from boundaries of the low-velocity fault zone. The waveform characteristics and arrival-time patterns in the record sections allow us to locate the boundaries of the fault zone and to determine its P - and S -wave velocities independently as well as its depth extent. Therefore, the trade-off between the fault-zone width and velocities can be avoided. Applying the method to the Landers waveform data reveals a low-velocity zone with a width of 270–360 m and a 35–60 per cent reduction in P and S velocities relative to the host rock. The analysis suggests that the low-velocity zone extends to a depth of ∼7 km. The western boundary of the low-velocity zone coincides with the observed main surface rupture trace.     

Determination and analysis of long-wavelength transition zone structure using SS precursors

Global mapping of 410 and 660 km discontinuity topography and transition zone thickness has proven to be a powerful tool for constraining mantle chemistry, dynamics and mineralogy. Numerous seismic and mineral physics studies suggest that the 410 km discontinuity results from the phase change of olivine to wadsleyite and the 660 km discontinuity results from the phase change of ringwoodite to perovskite and magnesiowustite. Underside reflections of the 410 and 660 km discontinuities arrive as precursors to SS . With the recent development of a semi-automated method of determining SS arrivals, we have more than tripled the Flanagan and Shearer (1998a) data set of handpicked SS waveforms. We are able to increase resolution by stacking waveforms in 5° rather than 10° radius bins as well as increasing data coverage significantly in the southern hemisphere. The resulting SS-S410S and SS-S660S times are heavily influenced by upper-mantle velocity structure. We perform a joint inversion for discontinuity topography and velocity heterogeneity as well as performing a simple velocity correction to the precursor differential times and find little difference between the two methods. The 660 km discontinuity topography and transition zone thickness are correlated with velocities in the transition zone whereas the 410 km discontinuity topography is not. In addition, the 410 km discontinuity topography is not correlated with the 660 km discontinuity topography, rather anticorrelated, as expected due to the opposite signs of the Clapeyron slopes of their respective phase changes. These results suggest that, whereas the topography of 660 km discontinuity could be dominated by thermal effects, the topography of the 410 km discontinuity is likely dominated by compositional effects. In addition, unlike previous studies which find less topography on the 410 km discontinuity than on the 660 km discontinuity, our 410 and 660 km topography have similar amplitudes.     

Traveltime measurements from noise correlation: stability and detection of instrumental time-shifts

We test the feasibility of using Green's functions extracted from records of ambient seismic noise to monitor temporal changes in the Earth crust properties by repeated measurements at regional distances. We use about 11 yr of continuous recordings to extract surface waves between three pairs of stations in California. The correlations are computed in a moving 1-month window and we analyse the temporal evolution of measured interstation traveltimes. The comparison of the arrival times in the positive and negative correlation time of Rayleigh and Love waves allows us to separate time-shifts associated with any form of physical change in the medium, those resulting from clock drift or other instrumental errors, and those due to change in the localization of the noise sources. This separation is based on the principle of time symmetry. When possible, we perform our analysis in two different period bands: 5–10 and 10–20 s. The results indicate that significant instrumental time errors (0.5 s) are present in the data. These time-shifts can be measured and tested by closure relation and finally corrected independently of any velocity model. The traveltime series show a periodic oscillation that we interpret as the signature of the seasonal variation of the region of origin of the seismic noise. Between 1999 and 2005, the final arrival time fluctuations have a variance of the order of 0.01 s. This allows us to measure interstation traveltimes with errors smaller than 0.3 per cent of the interstation traveltime and smaller than 1 per cent of the used wave period. This level of accuracy was not sufficient to detect clear physical variation of crustal velocity during the considered 11 yr between the three stations in California. Such changes may be more easily detectable when considering pairs of stations more closely located to each other and in the vicinity of tectonically active faults or volcanoes.     

Teleseismic Sn: a guided wave in the mantle

The short-period seismic phase Sn has been interpreted by Stephens & Isacks as a lid wave' in which the seismic energy is constrained to the uppermost few tens of kilometres of the mantle. We have extended their normal-mode interpretation for structures both with and without low-velocity zones (LVZ) in the upper mantle. We have used spherical, anelastic models of the Earth. For a model of an oceanic mantle with a LVZ, we agree that Sn is a lid wave for sources above 200–250 km, if only the onset of Sn is considered. The later portions of the Sn wave train sample the structure as deeply as the 420-km discontinuity. For deeper foci, the pseudo-lid wave does not appear to be generated; even the onset of Sn samples the deeper mantle structure. For a model of a continental mantle without a LVZ, in general, sources at all depths above the 420-km discontinuity appear to generate teleseismic Sn which samples the entire mantle as deeply as the discontinuity and which travels with a velocity significantly greater than the lid velocity. Thus the velocity of Sn may be an important diagnostic to determine whether or not a LVZ exists in the upper mantle.     

类型变更的相邻气象观测站的日气温资料整合

利用日气温资料,分析了因气象观测站类型变更而致的气候资料连续性和均一性问题及其影响,提出了订正相邻观测站日气温并将因观测站类型变更而致的相邻站不连续日气温资料整合的方法。主要结论有:(1)气象观测站类型变更而致的气候数据不连续既影响资料的均一性,也影响气候变化研究的结果。在考虑或不考虑测站类型变更(仅1980年以后)时,对我国过去60年1月气温趋势变化的估计结果差别达6.0%。(2)虽然海拔差异明显影响相邻站气温,但城乡差别等下垫面及其周边环境差异因素的作用极为显著;可导致两个地理位置相近的测站最大月气温差别超过0.50℃。(3)利用相邻测站的月气温差异进行各月的日气温订正可以消除海拔、台站下垫面及其周边环境差异对观测资料的均一性影响;使订正后的序列能更好地反映出气候的年际变化特征;从而可为我国正在开展的气候区划新方案、气候变化对区划影响及冷暖期环境格局变化等研究工作提供更均一的气候观测基础资料。     

Assessment of ski condition reliability in the Spanish and Andorran Pyrenees for the second half of the 20th century

In this study temporal trends of 14 climate and snow parameters related to ski conditions were analyzed for 11 ski stations located in the central Pyrenees (Spain and Andorra). We also investigated whether there was a temporal association for the analyzed parameters, such that the occurrence in a particular year of good (or bad) climate or snow conditions as represented by one parameter was similarly reflected by the other parameters. The lack of reliable climate and snow measurements was overcome by the use of simulated climate data retrieved from a high resolution hindcast simulation available for the period 1960–2006. These data were also used as inputs for an energy and mass snow energy model to obtain snow series. The results showed trends in ski reliability parameters for the 1960–2006 period. The number of days having a snowpack deeper than 30 cm and 100 cm showed declines at low and mid altitudes. The start of the ski season appears progressively delayed for all stations, and the ski season shortened. The frequency of rainy days increased at 3 stations and decreased at 8, while the frequency of days having heavy snowfall increased at 8 stations and declined at 3. Days having potential for snowmaking declined at all stations. The number of days having a wind-chill < −20 °C also decreased markedly, as overall did the number of days having a wind speed greater than 80th percentile. The main findings from the assessment of temporal associations between climate and snow parameters were positive correlations between snow depth and windy conditions. Seasons having a higher frequency of very cold days had a lower frequency of heavy snowfall and rainy days. Thus, the adverse effects on the ski industry of lesser snow availability may have been partially negated by the occurrence of fewer days of closure because of high winds, or other adverse meteorological factors.     

P-Wave Velocities in the Upper Mantle Beneath Fennoscandia and Western Russia

A detailed and extensive record section constructed from recordings at the NORSAR array of presumed explosions in continental Russia exhibits two distinct ( T , Δ) triplications. The reliable identification of these upper mantle travel-time branches is possible because of the dense areal sampling of the NORSAR configuration. A simple upper mantle P- velocity model which can account adequately for the data involves velocity discontinuities at depths of 420 km and 690 km, and fairly uniform velocity gradients elsewhere. For this model, the first arrival branch for Δ≤ 21° extends as a second arrival to a distance of about 33°, at which distance it is terminated by the 420-km discontinuity. Rays bottoming between depths of 420 and 690 km span the distance range 16°≤Δ≤ 28°, and give first arrivals in the range 21°≤, Δ, 24°. Rays which penetrate the 690-km discontinuity give rise to secondary arrivals in the range 19°≤Δ≤ 25°, and first arrivals for distances Δ≤ 25°.     

A lower mantle S-wave triplication and the shear velocity structure of D"

Summary. A lower mantle S-wave triplication detected with short- and long-period WWSSN and CSN recordings indicates a substantial shear velocity discontinuity near 280 km above the core–mantle boundary. The triplication can be observed in rotated SH seismograms from intermediate and deep focus events throughout the distance range from 70° to 95°. Three distinct source region–receiver array combinations that have been investigated in detail demonstrate consistent travel time and relative amplitude behaviour of the triplication, with slight systematic shifts in the triplication indicating up to 40 km variations in the depth of the discontinuity. Modelling of the observations with synthetic seismograms produced with the Cagniard de Hoop and reflectivity methods constrains the shear velocity increase to be 235 ± 0.25 per cent, comparable to upper mantle discontinuities. Short-period observations indicate that the velocity increase may be a sharp first-order discontinuity, or may extend over a transition zone no more than 50 km thick. The shear velocity gradient below the discontinuity, within the D" layer, is not well-constrained by the SH data, but slightly positive or near zero velocity gradients are consistent with the long-period amplitude ratios of ScSH/SH .     

The lithosphere–asthenosphere boundary beneath the western United States

S receiver functions from 67 broad-band seismic stations in the western United States clearly reveal the existence of a mantle discontinuity with velocity reduction downward, which we interpret as the lithosphere–asthenosphere boundary (LAB). The average depth of the LAB is ∼70 km. The boundary is relatively sharp with an overall sharpness of less than 20 km. The boundary is more prominent south of the Mendocino Triple Junction, where the Farallon Plate has completely subducted. This may indicate partial melts at the base of the lithosphere caused by the upwelling of the asthenospheric flow through the slab window. A double low velocity zone is observed at base of the lithosphere beneath southern Sierra Nevada, implying a second melting zone at a depth of ∼100 km, well correlated with previous studies of lithospheric delamination in the area.     

Upper mantle heterogeneities beneath Central Europe

Summary. An inversion of ISC travel-time data from selected earthquakes in the distance range 30°-90° to 53 stations in Central Europe has been used to model velocity down to 600 km depth. The model explains 0.1–0.2s of the residuals, as for other array studies, leaving 0.5 s unexplained as noise. The uppermost 100 km of the mantle and crust contains inhomogeneities that correlate remarkably well with the geology. This may be due to deep-seated thermal anomalies or, in some areas, to delays introduced by passage of the rays through sedimentary cover. The deeper anomalies are smaller and unrelated to those in the lithosphere, which suggests that the asthenosphere is decoupled from the rigid lithosphere. The structure at 600 km depth is again quite inhomogeneous and might be due to undulations of the 650 km discontinuity. The models show some suggestion of a high velocity slab trending from east to west beneath the Alps.     

Tectonic implications of transtensional supradetachment basin development in an extension‐parallel transfer zone: the Kocaçay Basin,western Anatolia,Turkey

The Kocaçay Basin (KÇB) is a key area in western Anatolia – a well‐known extended terrane where regional segmentation has received limited attention – for investigating strike‐slip faults kinematically linked to detachment faults. In this paper, we present results of an integrated sedimentologic, stratigraphic, and structural study of Miocene alluvial fan/fan‐delta/lacustrine deposits that accumulated in the KÇB, a NE‐trending basin with connections to the Menderes Metamorphic Core Complex (MCC). We mapped and evaluated most of the key faults in the KÇB, many for the first time, and recognised different deformation events in the study area near the E margin of the MCC. We also present field evidence for kinematic connections between low‐angle normal and strike‐slip faults which were developed in an intermittently active basement‐involved transfer zone in western Anatolia. We find that the KÇB contains a detailed record of Miocene transtensional sedimentation and volcanism that accompanied exhumation of the MCC. Structural data reveal that the basin was initially formed by transtension (D1 phase) and subsequently uplifted and deformed, probably as a result of early Pliocene wrench‐ to extension‐dominated deformation (D2 phase) overprinted by Plio‐Quaternary extensional tectonics (D3 phase). These results are consistent with progressive deformation wherein the axis of maximum extension remained in the horizontal plane but the intermediate and maximum shortening axes switched position in the vertical plane. Combining our results with published studies, we propose a new working hypothesis that the KÇB was a transtensional supradetachment basin during the Miocene. The hypothesis could provide new insights into intermittently active extension‐parallel zone of weakness in western Anatolia.These results also suggest that the termination of low‐angle normal fault systems within an extension parallel transfer zone may have resulted in a transtensional depressions which are different from classical supradetachment basins with respect to the sedimentation and deformational pattern of the basin infills.     

两高分辨率黄土剖面末次冰消期气候事件的差异探讨

高分辨率黄土剖面的获取使我们对亚轨道时间尺度上的气候事件的研究成为可能。通过对六盘山东、西部两个高分辨率黄土剖面与GRIP冰芯记录的对比研究,我们发现两剖面都捕捉到了末次冰消期内的Younger Dryas 事件与Heinrich-1事件及Blling和Allerd事件,并且它们在时间上与冰芯的结果较为一致。更为重要的是,这些亚轨道时间尺度上的气候事件在东西两剖面内存在较大差异:Heinrich-1事件在西面沙沟剖面的烙印较东面的王官剖面要深;对于Blling和Allerd事件,沙沟剖面较王官剖面的波动要多;在沙沟剖面,Younger Dryas 的结束是一种突变行为,而在王官剖面中则是一种渐变行为;此外,沙沟剖面比王官剖面表现出更大的波动幅度和更多千年尺度上的快速波动。初步推测上述差异是由高、低纬地区对它们影响的不同所造成的。     

IS THE EXTREME SNOW DEPTH RECORD FOR NORDLAND,NORWAY RELIABLE?

A snow depth of 370 cm at Dunderlandsdalen in winter 1919–20 is the largest recorded at stations in Nordland, Norway. During the period 1895–1924, the average maximum annual value there was 150.4 cm; at other stations it ranged from 38.5 to 190.1 cm. The ratio of maxima at other stations to that at Dunderlandsdalen was particularly low in 1919–20. In Nordland generally, that winter's recorded precipitation was slightly above the 30‐year average, but at Dunderlandsdalen it exceeded the average by 34%. At all stations except Dunderlandsdalen, 26 Jan.–1 Feb. was a dry period; at Dunderlandsdalen, 51.7 mm was recorded. Only one day without snowfall was recorded at Dunderlandsdalen between early January and early February, but elsewhere there were few days with snowfall. The difference in snowfall frequency and snow depth at Dunderlandsdalen in 1919–20 from values recorded elsewhere in Nordland contrasts with the relationships in other winters between 1895 and 1924. No observations were made at Dunderlandsdalen in winter 1917–18. Two of the householders there died in 1916. A change of personnel making the observations may have been responsible for the data gap and for the anomalous 1919 data. Changes made to buildings or the recording site in 1917 or 1918 may have resulted in increased snow depths as a result of drifting. Maintaining a record of climatic extremes and their environmental consequences is important. Data must be accurate. In view of this, it would be sensible to regard the validity of the 370 cm Dunderlandsdalen maximum as doubtful.     

ScSp observed on North Island, New Zealand: implications for subducting plate structure

Seismic phase conversions provide important constraints on the layered nature of subduction zone structures. Recordings from digital stations in North Island, New Zealand, have been examined for converted ScS ‐to‐ p ( ScSp ) arrivals from deep (>150 km) Tonga–Kermadec earthquakes to image layering in the underlying Hikurangi subduction zone. Consistent P ‐wave energy prior to ScS has been identified from stations in eastern and southern North Island, where the subducted plate interface is at a depth of between 15 and 30 km. Two ScS precursors are observed. Ray tracing indicates that the initial precursor ( ScSp ₁) corresponds to conversion from the base of an 11–14 km thick subducting Pacific crust. The second precursor is interpreted as a conversion from the top of the subducting plate. The amplitude ratio, ScSp ₁: ScS , increases from 0.10 to 0.19 from northern to southern North Island. This is within the range expected from a simple first‐order velocity discontinuity at an oceanic Moho. A 1–2 km thick layer of low‐velocity sediment at the top of the subducting plate is required to explain the remaining ScSp waveform. Our results imply that the abnormally thick Hikurangi–Chatham Plateau has been subducting beneath New Zealand for at least 2.9 Myr, thus explaining the high uplift rates observed across eastern North Island.     

距今五千年左右环境恶化事件对我国新石器文化的影响及其原因的初步探讨

距今五千年左右全球普遍存在一次变化幅度较大的突发性环境恶化事件 ,此事件导致各地生态环境发生了较大的变化 ,与此同时的我国新石器文化表现为衰退和在某些地点的中断。气候环境变化是该时段影响人类文化发展变化的主要因素。环境变化对史前文化影响的可能原因是 :气候变化导致生态环境变化 ,生态环境变化导致直接影响人类生存的生物营养源的变化 ,进而导致史前文化的变化。     

CO2 seasonal variation and global change: Test global warming from another point of view

CO₂ and temperature records at Mauna Loa, Hawaii, and other observation stations show that the correlation between CO₂ and temperature is not significant. These stations are located away from big cities, and in various latitudes and hemispheres. But the correlation is significant in global mean data. Over the last five decades, CO₂ has grown at an accelerating rate with no corresponding rise in temperature in the stations. This discrepancy indicates that CO₂ probably is not the driving force of temperature change globally but only locally (mainly in big cities). We suggest that the Earth''s atmospheric concentration of CO₂ is too low to drive global temperature change. Our empirical perception of the global warming record is due to the urban heat island effect:temperature rises in areas with rising population density and rising industrial activity. This effect mainly occurs in the areas with high population and intense human activities, and is not representative of global warming. Regions far from cities, such as the Mauna Loa highland, show no evident warming trend. The global monthly mean temperature calculated by record data, widely used by academic researchers, shows R²=0.765, a high degree of correlation with CO₂. However, the R² shows much less significance (mean R²=0.024) if calculated by each record for 188 selected stations over the world. This test suggests that the inflated high correlation between CO₂ and temperature (mean R²=0.765-0.024=0.741) used in reports from the Intergovernmental Panel on Climate Change (IPCC) was very likely produced during data correction and processing. This untrue global monthly mean temperature has created a picture:human emission drives global warming.     

Spatial and temporal variability of daily precipitation in Haihe River basin, 1958–2007

The seasonal variability and spatial distribution of precipitation are the main cause of flood and drought events. The study of spatial distribution and temporal trend of precipitation in river basins has been paid more and more attention. However, in China, the precipitation data are measured by weather stations (WS) of China Meteorological Administration and hydrological rain gauges (RG) of national and local hydrology bureau. The WS data usually have long record with fewer stations, while the RG data usually have short record with more stations. The consistency and correlation of these two data sets have not been well understood. In this paper, the precipitation data from 30 weather stations for 1958–2007 and 248 rain gauges for 1995–2004 in the Haihe River basin are examined and compared using linear regression, 5-year moving average, Mann-Kendall trend analysis, Kolmogorov-Smirnov test, Z test and F test methods. The results show that the annual precipitation from both WS and RG records are normally distributed with minor difference in the mean value and variance. It is statistically feasible to extend the precipitation of RG by WS data sets. Using the extended precipitation data, the detailed spatial distribution of the annual and seasonal precipitation amounts as well as their temporal trends are calculated and mapped. The various distribution maps produced in the study show that for the whole basin the precipitation of 1958–2007 has been decreasing except for spring season. The decline trend is significant in summer, and this trend is stronger after the 1980s. The annual and seasonal precipitation amounts and changing trends are different in different regions and seasons. The precipitation is decreasing from south to north, from coastal zone to inland area.     

Late Miocene – Recent exhumation of the central Himalaya and recycling in the foreland basin assessed by apatite fission-track thermochronology of Siwalik sediments, Nepal

Thermochronological analysis of detrital sediments derived from the erosion of mountain belts and contained in the sedimentary basins surrounding them allows reconstructing the long-term exhumation history of the sediment source areas. The effective closure temperature of the thermochronological system analysed determines the spatial and temporal resolution of the analysis through the duration of the lag time between closure of the system during exhumation and its deposition in the sedimentary basin. Here, we report apatite fission-track (AFT) data from 31 detrital samples collected from Miocene to Pliocene stratigraphic sections of the Siwalik Group in western and central Nepal, as well as three samples from modern river sediments from the same area, that complement detrital zircon fission-track (ZFT) and U–Pb data from the same samples presented in a companion paper. Samples from the upper part of the stratigraphic sections are unreset and retain a signal of source-area exhumation; they show spatial variations in source-area exhumation rates that are not picked up by the higher-temperature systems. More deeply buried samples have been partially reset within the Siwalik basin and provide constraints on the thermal and kinematic history of the fold-and-thrust belt itself. The results suggest that peak source-area exhumation rates have been constant at ∼1.8 km Myr⁻¹ over the last ∼7 Ma in central Nepal, whereas they ranged between 1 and ∼1.5 km Myr⁻¹ in western Nepal over the same time interval; these spatial variations may be explained by either a tectonic or climatic control on exhumation rates, or possibly a combination of the two. Increasing lag times within the uppermost part of the sections suggest an increasing component of apatites that have been recycled within the Siwalik belt and are corroborated by AFT ages of modern river sediment downstream as well as the record of the distal Bengal Fan. The most deeply buried and most strongly annealed samples record onset of exhumation of the frontal Siwaliks along the Himalayan frontal thrust at ∼2 Ma and continuous shortening at rates comparable with the present-day shortening rates from at least 0.3 Ma onward.     

Central‐Auckland rainfall, 1853–2020: Sites histories and implications for developing a long‐term rainfall record

Previous work related to central‐Auckland site Albert Park showed that sparse site histories for inner‐city weather stations can be usefully supplemented by mining searchable online historical archives (Papers Past, DigitalNZ), at least as far back as 1909. Assessment of potential is extended here back to 1853, within the context of developing a long‐term rainfall record for Auckland. Although late 19th century sources are thinner and of poorer quality, much useful information could be gleaned about the locations and exposures of early sites. Hypothesised implications for the rainfall record suggest that all early site data need to be adjusted upwards.