大气垂直折射与时序垂直角值的谱分析

本文介绍了对时序垂直角观测值进行谱分析的方法和所得到的结论。使以往对于大气折射影响及其变化规律的定性认识得以定量化。     

Predicting ground-motion variations at the Turkey-flat test site,California

Results of a single group participating in an international experiment are analyzed. The experiment served to verify computational predictions of the ground-motion variations due to near-surface geological effects at a site established for that purpose by the California Department of Conservation. Based on an acceleration record at a rock location, and geotechnical model of medium, records at the other locations of a nearby sedimentary deposit were predicted. A 2-D finite-difference sensitivity analysis suggested that the lateral wave-propagation effects are negligibly small, and locally 1-D computations are sufficient for the present site. Those computations are compared with observations not available to the authors during the blind prediction. Peak accelerations, peak velocities and RMS accelerations were predicted with errors less than 159%, 114% and 62%, respectively. Maxima of the response spectra were fitted within a factor of 2. The predicted and observed Husid's plots (i.e., the normalized cumulative plots of the acceleration squared) have the correlation coefficients 0.98. The detected misfits do not show any simple relation to the instrument location, component, frequency, or time.     

In-Situ Observation of Fluid Mud in the North Passage of Yangtze Estuary, China

Observations of fluid mud were made in the lower North Passage of the Yangtze Estuary in February 2000, on 10 -11 August 2000, on 30 - 31 August 2000 (after two strong typhoons), on 21 - 24 August 2000 (neap tide) and on 3 -6 September 2000 (mean tide) respectively. In situ data show that the fluid mud in this area consists of fine cohesive sediment (median size 7.23 μm). The formation and movement of fluid mud varied during the neap-spring and flood-ebb tidal cycle. Observations suggest that fluid mud phenomena in this area may be categorised in a three-fold manner as slack water, storm and saltwedge features. The thickness of the fluid mud layer of slack water during the neap tide ranged from 0.2 to 0.96 m, whereas during the mean tide, the thickness ranged from 0.17 to 0.73 m, and the thickness of the fluid mud layer was larger during slack water than at the flood peak. Shoals cover an area of 800 km^2 with a water depth smaller than 5 m. Erosion of these extensive intertidal mudflats due to storm action provides an abundant sediment source. This is particularly significant in this estuary when the tidal level is lower than 5 m. The lower North Passage is a typical zone of saltwater wedging, so the saltwedge fluid mud has the most extensive spatial range in the estuary.     

利用GPS组合观测值探测与修复周跳算法实现

本文在目前常用的周跳探测和修复方法的基础上,根据相位减伪距组合,Melbourne-Wübbena组合,电离层残差组合三种GPS组合观测值优缺点互补,将这三种组合观测值有机结合产生可靠的探测与修复周跳算法并用程序实现,通过实例分析证明此算法能够探测和修复任意大小的周跳。     

Long-term seismic noise investigations on Shikotan Island: First results

In 2003–2004, long-term seismic noise observations were launched on Shikotan Island (Lesser Kuril Range) based on the “Shikotan” dormant regional seismic station. The geological and geophysical data on the registration area are reported. Information about the equipment and its technical specifications is given. The precursors to the strongest local earthquakesthat occurred in the Shikotan Island region in January 2005–March 2007 are identified.     

Estimation of covariance matrix from geochemical data with observations below detection limits

Multivariate statistical analyses have been extensively applied to geochemical measurements to analyze and aid interpretation of the data. Estimation of the covariance matrix of multivariate observations is the first task in multivariate analysis. However, geochemical data for the rare elements, especially Ag, Au, and platinum-group elements, usually contain observations the below detection limits. In particular, Instrumental Neutron Activation Analysis (INAA) for the rare elements produces multilevel and possibly extremely high detection limits depending on the sample weight. Traditionally, in applying multivariate analysis to such incomplete data, the observations below detection limits are first substituted, for example, each observation below the detection limit is replaced by a certain percentage of that limit, and then the standard statistical computer packages or techniques are used to obtain the analysis of the data. If a number of samples with observations below detection limits is small, or the detection limits are relatively near zero, the results may be reasonable and most geological interpretations or conclusions are probably valid. In this paper, a new method is proposed to estimate the covariance matrix from a dataset containing observations below multilevel detection limits by using the marginal maximum likelihood estimation (MMLE) method. For each pair of variables, sayY andZ whose observations containing below detection limits, the proposed method consists of three steps: (i) for each variable separately obtaining the marginal MLE for the means and the variances, , , , and forY andZ: (ii) defining new variables by and and lettingA=C+D andB=C–D, and obtaining MLE for variances, and forA andB; (iii) estimating the correlation coefficient _YZ by and the covariance _YZ by . The procedure is illustrated by using a precious metal geochemical data set from the Fox River Sill, Manitoba, Canada.     

Analysis of COMPTEL gamma-ray burst locations and spectra

In its first three years of operation, the COMPTEL instrument on theCompton Gamma-Ray Observatory has measured the locations (mean accuracy 1°) and spectra (0.75-30 MeV) of 18 gamma-ray bursts and continues to observe new events at a rate of 1/month. With good angular resolution and sensitivity at MeV energies, the growing COMPTEL burst catalog is an important new piece of evidence in the on-going GRB mystery. The COMPTEL burst locations are consistent with an isotropic distribution of sources, yet the spatial coincidence of two of the bursts indicates the possibility of repetition. The COMPTEL burst spectra are in most cases consistent with a single power law model with spectral index in the range 2–3. However, two bursts show evidence of a spectral break in the MeV range. Measurement of rapid variability at MeV energies in the stronger bursts provides evidence that either the sources are nearby (within the Galaxy) or the gamma-ray emission is relativistically beamed. We present an overview of analysis results obtained from the COMPTEL burst catalog concentrating on the search for burst repetition and the implications of highly variable MeV emission.     

Deep structure of the Nojima Fault, southwest Japan, estimated from borehole observations of fault-zone trapped waves

To estimate the deep structure of the southern part of the Nojima Fault, southwest Japan without the influence of near-surface structures, we analyzed the Love-wave-type fault-zone trapped waves (LTWs) recorded by a borehole seismometer at 1800 m depth. We examined the polarization, dispersion, and dominant frequency of the wavetrain following the direct S-wave in each seismogram to identify the LTW. We selected eight candidates for typical LTWs from 462 records. Because the duration of the LTW increases with hypocentral distance, we infer that the low velocity fault-zone of the Nojima Fault continues towards the seismogenic depth. In addition, since the duration of the LTW increases nonlinearly with hypocentral distance, we infer that the S-wave velocity of the fault-zone increases with depth. The location of events showing the LTW indicates that the fault-zone dips to the southeast at 75° and continues to a depth of approximately 10 km. We assumed a uniform low-velocity waveguide to estimate the average structure of the fault-zone. We estimated the average width, S-wave velocity, and Q_s of the fault-zone by comparing an analytical solution of the LTW with measured data. The average width, S-wave velocity, and Q_s of the fault-zone are 150 to 290 m, 2.5 to 3.2 km/s, and 40 to 90, respectively. Hence the fault-zone structure with a larger width and smaller velocity reduction than the fault-zone model estimated by previous surface observation is more suitable to represent the average fault-zone structure of the Nojima fault. The present study also indicated that the shallow layers and/or a shallow fault-zone structure drastically changes the characteristics of the LTW recorded at the surface, and therefore cause a discrepancy in the fault-zone model between the borehole observation and surface observation.     

Observations And Modelling Of The On-Ice And Off-Ice Air Flow Over The Northern Baltic Sea

Two cases of on-ice and off-ice air flow characterizing the opposite weather situations over the ice-edge zone in the northern Baltic Sea are analysed on the basis of aircraft observations, and modelled using atwo-dimensional mesoscale model. The stable boundary layer (SBL) during theon-ice flow exhibited little thermal modification, but a low-level jet (LLJ) was generated at the 250-m high top of the SBL. In the model, the LLJ was associated with inertial oscillations in space, while the baroclinicity explained the shape of the wind profile well above the SBL. Although the observed LLJ was most pronounced over the ice, the modelling suggests that it was not generated by the ice edge but by the coastline some 400 km upwind of the ice edge, where a much more drastic change in the thermal stratification and surface roughness took place. The generation, maintenance, and strength of the LLJ were very sensitive to the parameterization of turbulent mixing in the SBL. In the case of the off-ice flow, the modification of the air mass and the development of a convective boundary layer (CBL) both over the ice and open sea were reasonably well modelled. Sensitivity runs suggested that it was essential to take into account the effects of subgrid-scale leads, a forest in the archipelago (which was crossed by the air flow), and water vapour condensationinto ice crystals. The heat flux from leads was particularly important for the heatbudget of the CBL, and the observed growth of the CBL was partly due to theeffective mixing over the rough and relatively warm forest.     

On the circulation in the East Frisian Wadden Sea: numerical modeling and data analysis

 In this paper we use a combination of numerical modeling and data analysis to gain a better understanding of the major characteristics of the circulation in the East Frisian Wadden Sea. In particular, we concentrate on the asymmetry of the tidal wave and its modulation in the coastal area, which results in a complex pattern of responses to the sea-level forcing from the North Sea. The numerical simulations are based on the 3-D primitive equation General Estuarine Transport Model (GETM) with a horizontal resolution of 200 m and terrain-following vertical coordinates. The model is forced at its open boundaries with sea-level data from an operational model for the German Bight (German Hydrographic Office). The validation data for our model simulations include time series of tidal gauge data and surface currents measured at a pile in the back-barrier basin of the Island Langeoog, as well as several ADCP transects in the Accumer Ee tidal inlet. Circulation and turbulence characteristics are investigated for typical situations driven by spring and neap tides, and the analysis is focused on dominating temporal and spatial patterns. By investigating the response of five back-barrier basins with rather different morphologies to external forcing, an attempt is made to elucidate the dominating physical balances controlling the circulation in the individual sub-basins. It is demonstrated that the friction at the seabed tends to slow down the tidal signal in the shallow water. This leads to the establishment of flood dominance in the shallow sea north of the barrier islands. South of the islands, where the water volume of the channels at low tide is smaller than the tidal prism, the asymmetry of the tidal signal is shifted towards ebb dominance, a feature which is particularly pronounced at spring tide. At the northern open boundary, the tidal wave propagating from west to east generates a sea-level difference of ∼1 m along the boundary, and thereby triggers vigorous alongshore currents. The frictional control in the model is located in the inlets, as well as along the northern boundary. The correlation between velocity and turbulent kinetic energy tends to the establishment of a net southward transport, giving theoretical support to the observed accumulation of sediments on the intertidal flats. Weak turbulence along the northern shores of the barrier islands and the small magnitude of the residual currents there promote accumulation of suspended matter in these areas, although wave action will generally counteract this effect. Received: 29 May 2002 / Accepted: 26 September 2002 Responsible Editor: Jean-Marie Beckers Acknowledgements We are indebted to S. Dick for providing the data from the operational model of BSH and to B. Flemming for the useful discussions. The topography data and Fig. 1 have been prepared in cooperation with F. Meyer. Figure 2 has been prepared by G. Brink-Spalink. We also thank for the comments from an anonymous reviewer which helped to improve our paper.