Development of an Information and Database System for the IERS: Status and outlook

The focus of this paper is to present the concepts for a new dynamic, i.e. database-driven information system for the International Earth Rotation and Reference Systems Service (IERS). Its goal is to coordinate the data and information flow between the participating institutions. Therefore, all relevant data and products of the IERS should be archived to guarantee their long-term availability. At the same time the consistency and timeliness of the contents of the information system can be guaranteed by managing the information in a database. Additionally, the metadata of all products (including the IERS publications) should be modelled in a database to allow the users to search for specific data or topics with respect to space, time and content. In order to be able to link related data the heterogeneous formats of the products have to be transformed into a common format. The eXtensible Markup Language (XML) is going to be used, to perform this ambitious task. The usage of XML not only allows to link related data but also to easily exchange data and to present them in different ways like Web pages, pdf, etc. The system will be completed by several tools to manage and coordinate the tasks of the Central Bureau (CB) and by an information system with respect to IERS-related topics. The standardised database tools will also allow an easy exchange of information with and links to other databases within the IGGOS project to realize a powerful instrument to serve the Earth observing system.     

Analysis of change of red tide species in Yodo River estuary by the numerical ecosystem model

Occurrence number of red tides in Osaka Bay in Japan is more than 20 cases every year. Diatom red tide was dominant in Osaka Bay, but the non-diatom red tide was dominant in early 1990s. Therefore, the material cycling in Yodo River estuary in Osaka Bay during August from 1991 to 2000 was analyzed by using the numerical ecosystem model and field observation data to clarify the reasons of change in red tide species. Year-to-year variation in calculated concentration ratio of diatom to non-diatom corresponds to the variation in observed ratio of red tide days of diatom to non-diatom. Limiting nutrient of primary production is phosphate over the period. Diatom dominated from 1991 to 1993, but it was difficult for non-diatom to grow due to the limitation by physical condition. Non-diatom was able to grow because of good physical and nutrient conditions from 1994 to 1996. And diatom dominated again under the good physical condition, and phosphorus supply was not enough for non-diatom to grow from 1998 to 2000. Phosphate concentration in the lower layer of Yodo River estuary was important to the variation in red tide species in the upper layer of Yodo River estuary.     

Discussion on the energy partition of volcanic eruption

Energy of volcanic eruption is divided into various ways, such as lava flows, ash production, explosion, ground vibration, air waves and so on. In order to elucidate the cruptive mechanism, energy partition of volcanic eruption in unit time is considered to be a key to solve the problem. In this short paper, partition of eruption energy was calculated for typical three basaltic volcanoes. Ratio of tremor energy to thermal energy is almost the same for both Mihara and Nyragongo in spite of their different eruptive mechanism. Tremor energy was found to be small for Kilauea Volcano compared with the above two volcanoes. This was attributed to the different source mechanism of tremor. Partition of eruption energy was calculated for Also Volcano which showed phreatic explosion in 1958. Kinetic energy is almost ten times of thermal energy. Ratio of energy converted to elastic waves to thermal energy is found to be equivalent to the ratio of energy of elastic waves to the chemical energy in case of explosion of chemical explosives.     

一种适用于大地震临震预测的地声监测系统

介绍了震前地声的历史记载以及震前地声的产生机制和特征．还有如何利用震前地声观测结果对地震发震的三要素进行预测。并针对震前地声监测所存在的问题、困难和发展趋势,提出了一种适用于大地震临震预测的地声监测系统,该系统可详尽、完整地记录震前地声的强度、频度、频谱,并可通过在地震易发区域的大规模密集布点,组成这一区域的震前地声监测网络,以监测大地震孕育过程,以及对大地震的三要素进行预测预报。     

Use of a shared tuned mass damper (STMD) to reduce vibration and pounding in adjacent structures

The dynamic response of tall civil structures due to earthquakes is very important to civil engineers. Structures exposed to earthquakes experience vibrations that are detrimental to their structural components. Structural pounding is an additional problem that occurs when buildings experience earthquake excitation. This phenomena occurs when adjacent structures collide from their out‐of‐phase vibrations. Many energy dissipation devices are presently being used to reduce the system response. Tuned mass dampers (TMD) are commonly used to improve the response of structures. The stiffness and damping properties of the TMD are designed to be a function of the natural frequency of the building to which it is connected. This research involves attaching adjacent structures with a shared tuned mass damper (STMD) to reduce both the structures vibration and probability of pounding. Because the STMD is connected to both buildings, the problem of tuning the STMD stiffness and damping parameters becomes an issue. A design procedure utilizing a performance function is used to obtain the STMD parameters to result in the best overall system response. Copyright © 2001 John Wiley & Sons, Ltd.     

Detecting infiltration and impacts of introduced water using strontium isotopes

Water introduced to surface drainages, such as agricultural and roadway runoff, mine drainage, or coalbed natural gas (CBNG)-produced water, potentially can be of environmental concern. In order to mitigate potential environmental effects, it may be important to be able to trace water discharged to the surface as it infiltrates and interacts with near-surface aquifers. We have chosen to study water withdrawn during CBNG production for isotope tracing in the hyporheic zone because it poses a variety of economic, environmental, and policy issues in the Rocky Mountain states. Ground water quality must be protected as CBNG water is added to semiarid ecosystems. Strontium (Sr) isotopes are effective fingerprints of the aquifer from which water originates. In this study, CBNG water was found to have a higher (87)Sr/(86)Sr ratio than the local alluvial aquifer water. This measurable difference allows the strontium isotope ratio and concentration to be used as tracers of CBNG water following its discharge to the surface. The dissolution and mobilization of salts from soil are an important contributor to ground water quality degradation. In the Powder River basin of Wyoming, the soils are calcium carbonate-buffered systems. The chemical similarity of strontium to calcium allows it to substitute into calcium minerals and enabled us to use strontium isotopes to identify calcium salts mobilized from the soil. Strontium isotopes are an effective monitor of the source of ions and the volume and direction of introduced water flow in the hyporheic zone.     

Model and parameter uncertainty in IDF relationships under climate change

Quantifying distributional behavior of extreme events is crucial in hydrologic designs. Intensity Duration Frequency (IDF) relationships are used extensively in engineering especially in urban hydrology, to obtain return level of extreme rainfall event for a specified return period and duration. Major sources of uncertainty in the IDF relationships are due to insufficient quantity and quality of data leading to parameter uncertainty due to the distribution fitted to the data and uncertainty as a result of using multiple GCMs. It is important to study these uncertainties and propagate them to future for accurate assessment of return levels for future. The objective of this study is to quantify the uncertainties arising from parameters of the distribution fitted to data and the multiple GCM models using Bayesian approach. Posterior distribution of parameters is obtained from Bayes rule and the parameters are transformed to obtain return levels for a specified return period. Markov Chain Monte Carlo (MCMC) method using Metropolis Hastings algorithm is used to obtain the posterior distribution of parameters. Twenty six CMIP5 GCMs along with four RCP scenarios are considered for studying the effects of climate change and to obtain projected IDF relationships for the case study of Bangalore city in India. GCM uncertainty due to the use of multiple GCMs is treated using Reliability Ensemble Averaging (REA) technique along with the parameter uncertainty. Scale invariance theory is employed for obtaining short duration return levels from daily data. It is observed that the uncertainty in short duration rainfall return levels is high when compared to the longer durations. Further it is observed that parameter uncertainty is large compared to the model uncertainty.     

对黄河流域水资源保护的探索与实践

黄河水资源贫乏,水污染严重,水资源供需矛盾日益尖锐,黄河流域水资源保护面临着诸多的困难和压力.新《水法》的颁布,给黄河水资源保护带来新的机遇和挑战,依法建立法制化、社会化的流域与区域相结合的水资源保护管理体系,完善法规,建立健全执法体系,以水功能区划管理为重心,实施入河总量控制,建立和完善重大水污染事件快速反应机制,加强水质监测机制与技术创新,加强水资源保护前期和科研工作,建设水资源保护信息管理系统,采取多部门联合治污等项措施,是做好黄河流域水资源保护工作的重要前提和任务.     

Review of vibroseis data acquisition and processing for better amplitudes: adjusting the sweep and deconvolving for the time-derivative of the true groundforce

The goal of vibroseis data acquisition and processing is to produce seismic reflection data with a known spatially-invariant wavelet, preferably zero phase, such that any variations in the data can be attributed to variations in geology. In current practice the vibrator control system is designed to make the estimated groundforce equal to the sweep and the resulting particle velocity data are cross-correlated with the sweep. Since the downgoing far-field particle velocity signal is proportional to the time-derivative of the groundforce, it makes more sense to cross-correlate with the time-derivative of the sweep. It also follows that the ideal amplitude spectrum of the groundforce should be inversely proportional to frequency. Because of non-linearities in the vibrator, bending of the baseplate and variable coupling of the baseplate to the ground, the true groundforce is not equal to the pre-determined sweep and varies not only from vibrator point to vibrator point but also from sweep to sweep at each vibrator point. To achieve the goal of a spatially-invariant wavelet, these variations should be removed by signature deconvolution, converting the wavelet to a much shorter zero-phase wavelet but with the same bandwidth and signal-to-noise ratio as the original data. This can be done only if the true groundforce is known. The principle may be applied to an array of vibrators by employing pulse coding techniques and separating responses to individual vibrators in the frequency domain. Various approaches to improve the estimate of the true groundforce have been proposed or are under development; current methods are at best approximate.     

Remediation of Heterogeneous Aquifers Subject to Uncertainty

Optimal cost pump-and-treat ground water remediation designs for containment of a contaminated aquifer are often developed using deterministic ground water models to predict ground water flow. Uncertainty in hydraulic conductivity fields used in these models results in remediation designs that are unreliable. The degree to which uncertainty contributes to the reliability of remediation designs as measured by the characterization of the uncertainty is shown to differ depending upon the geologic environments of the models. This conclusion is drawn from the optimal design costs for multiple deterministic models generated to represent the uncertainty of four distinct models with different geologic environments. A multi scenario approach that includes uncertainty into the remediation design called the deterministic method for optimization subject to uncertainty (DMOU) is applied to these distinct models. It is found that the DMOU is a method for determining a remediation design subject to uncertainty that requires minimal postprocessing efforts. Preprocessing, however, is required for the application of the DMOU to unique problems. In the ground water remediation design problems, the orientation of geologic facies with respect to the orientation of flow patterns, pumping well locations, and constraint locations are shown to affect the preprocessing, the solutions to the DMOU problems, and the computational efficiency of the DMOU approach. The results of the DMOU are compared to the results of a statistical analysis of the effects of the uncertainty on remediation designs. This comparison validates the efficacy of the DMOU and illustrates the computational advantages to using the DMOU over statistical measures.     

Building over faults: A procedure for evaluating risk

This paper examines the risk from fault offset to a building constructed over the surface trace of a fault. An expression relating the magnitude of surface offset to its return period is derived using classical seismic hazard analysis. The risk from fault offset is compared with that due to ground shaking by comparing offsets and MM intensities with equal likelihoods of occurrence. The comparison of risks is useful since it enables risk from fault movement to be related to levels of community acceptance as expressed in building codes for design against ground shaking effects. Relative risk is independent of fault activity rate, but does depend on regional constants, principally those describing the attenuation of MM intensities. For a site over an isolated fault in the Wellington, New Zealand region the risk of damage from fault rupture appears to be about equal to that from ground shaking for buildings with poor earthquake resistance (e.g. unreinforced masonry structures). But for structures designed to have some resistance to ground shaking, the risk due to fault rupture is considerably greater than that due to ground shaking.     

A note on the deflection of a baroclinic current by a continental shelf

The deflection, at a step-shelf fronted coast, of a constant potential vorticity current in a reduced-gravity, inviscid model ocean is studied theoretically. The step shelf, with a depth smaller than the reservoir depth, forces the uplifting of the approaching current and causes water column foreshortening, leading to relative vorticity generation that enhances current deflection to the right (facing the coast). As a consequence, in comparison to the case of a vertical wall coast, the proportion of the transport to the right is increased. For normal incidence for a shelf-depth/reservoir-depth ratio of 0.3 and shelf width to deformation radius ratio of 1.5, more than 90% of the approaching current transport goes to the right and less than 10% to the left. In addition, the (barotropic) dynamic pressure at the coast is low to the right and high to the left (with the highest pressure at the stagnation point). In the vertical wall case, the wall pressures on the flank are equal. For oblique incidence from the left, the deflection to the left is drastically reduced. In fact, there is practically no steady-state flow diverted to the left (less than 2%) when the approach angle is greater than 60° to the left of normal. In the vertical wall case, the same angle would have to be 90° for the flow to the left to vanish, namely only when the approach current is parallel to the coast to the right.     

A Basic Water Quality Model for the River Neckar: Part 2 – model‐based analysis of the oxygen budget and scenarios

In this paper the Basic Water Quality Model (BWQM) for the central part of River Neckar is used to analyse the oxygen budget and to assess the potentials of various measures to prevent or mitigate critical dissolved oxygen (DO) declines. It is shown that the oxygen budget is mainly governed by phytoplankton dynamics. The excessive growth of algae and the sudden break down of the resulting algal blooms may cause episodic DO depressions. Therefore, to stabilise the oxygen budget in a sustainable way, eutrophication has to be controlled within the central part of River Neckar and the upstream regions. The only feasible way to reach this goal appears to be a further drastic reduction of phosphorus emissions. In addition, it is indispensable to hold the very high standards of biochemical oxygen demand and ammonium retention at the wastewater treatment plants. A worse performance of the treatment plants would dramatically aggravate critical DO declines which may be caused by algae dynamics. As long as the oxygen budget is not completely stabilised, weir and turbine aeration can be used to mitigate DO depressions. It could be shown that the potentials of these measures suffice to keep DO at a tolerable level. However, due to the long travel times in River Neckar, it is important to start aeration up to several days before the DO minimum is reached.     

Elastic models of seismic failure in the Valley of Mexico

Restricted parts of the Valley of Mexico have been plagued by swarms of small earthquakes during recent years. The cause of these events is puzzling, but they may relate to massive artificial changes in the hydrology of the lake-bed which underlies most of Mexico City. If the valley responds elastically, plane strain finite element calculations can be used to relate geologic and seismic structure to observed seismicity and overall stability of the valley. These calculations show that unstable zones and areas of local seismicity coincide if the triggering mechanism of the seismicity is related to variations of pore pressure and water load.In order to make such a connection it is necessary to ascribe the generating mechanism of the seismicity to the density and elastic modulus anomalies associated with the valley structure, and to suggest that the phenomenon that triggers this seismicity is related to changes in water content of the near surface formations.     

Use of GRACE determined secular gravity rates for glacial isostatic adjustment studies in North-America

Monthly geopotential spherical harmonic coefficients from the GRACE satellite mission are used to determine their usefulness and limitations for studying glacial isostatic adjustment (GIA) in North-America. Secular gravity rates are estimated by unweighted least-squares estimation using release 4 coefficients from August 2002 to August 2007 provided by the Center for Space Research (CSR), University of Texas. Smoothing is required to suppress short wavelength noise, in addition to filtering to diminish geographically correlated errors, as shown in previous studies. Optimal cut-off degrees and orders are determined for the destriping filter to maximize the signal to noise ratio. The halfwidth of the Gaussian filter is shown to significantly affect the sensitivity of the GRACE data (with respect to upper mantle viscosity and ice loading history). Therefore, the halfwidth should be selected based on the desired sensitivity.It is shown that increase in water storage in an area south west of Hudson Bay, from the summer of 2003 to the summer of 2006, contributes up to half of the maximum estimated gravity rate. Hydrology models differ in the predictions of the secular change in water storage, therefore even 4-year trend estimates are influenced by the uncertainty in water storage changes. Land ice melting in Greenland and Alaska has a non-negligible contribution, up to one-fourth of the maximum gravity rate.The estimated secular gravity rate shows two distinct peaks that can possibly be due to two domes in the former Pleistocene ice cover: west and south east of Hudson Bay. With a limited number of models, a better fit is obtained with models that use the ICE-3G model compared to the ICE-5G model. However, the uncertainty in interannual variations in hydrology models is too large to constrain the ice loading history with the current data span. For future work in which GRACE will be used to constrain ice loading history and the Earth's radial viscosity profile, it is important to include realistic uncertainty estimates for hydrology models and land ice melting in addition to the effects of lateral heterogeneity.     

我国科学数据共享现状的调查

本文是科学数据共享政策法规体系研究系列文章的收笔之作.撰写本文的目的,就是为说明科学数据共享政策法规体系的研究是建立在大量的基础调研之上.调查研究是谋事之基、成事之道.从实践中寻求解决现实问题的方法,才能制定出切实可行的政策、法规,调查研究是制定政策、法规的客观依据.为使制订政策、 规章制度适应社会经济发展和变化的新情况,就深入开展调查研究,认真听取各方面的意见,深入研究相应的管理规律和运行机制.本文详述了调查过程、调查方法,以佐证政策法规研究不能脱离实际.     

Analysis of thick, non-planar boundaries using the discontinuity analyser

The advent of missions comprised of phased arrays of spacecraft, with separation distances ranging down to at least mesoscales, provides the scientific community with an opportunity to accurately analyse the spatial and temporal dependencies of structures in space plasmas. Exploitation of the multi-point data sets, giving vastly more information than in previous missions, thereby allows unique study of their small-scale physics. It remains an outstanding problem, however, to understand in what way comparative information across spacecraft is best built into any analysis of the combined data. Different investigations appear to demand different methods of data co-ordination. Of the various multi-spacecraft data analysis techniques developed to affect this exploitation, the discontinuity analyser has been designed to investigate the macroscopic properties (topology and motion) of boundaries, revealed by multi-spacecraft magnetometer data, where the possibility of at least mesoscale structure is considered. It has been found that the analysis of planar structures is more straightforward than the analysis of non-planar boundaries, where the effects of topology and motion become interwoven in the data, and we argue here that it becomes necessary to customise the analysis for non-planar events to the type of structure at hand. One issue central to the discontinuity analyser, for instance, is the calculation of normal vectors to the structure. In the case of planar and ‘thin’ non-planar structures, the method of normal determination is well-defined, although subject to uncertainties arising from unwanted signatures. In the case of ‘thick’, non-planar structures, however, the method of determination becomes particularly sensitive to the type of physical sampling that is present. It is the purpose of this article to firstly review the discontinuity analyser technique and secondly, to discuss the analysis of the normals to thick non-planar structures detected in magnetometer data.     

Integrated assessment for establishing an oil environmental vulnerability map: Case study for the Santos Basin region,Brazil

The growth of maritime transport and oil exploitation activities may increase the risk of oil spills. Thus, plans and actions to prevent or mitigate impacts are needed to minimize the effects caused by oil. However, tools used worldwide to support contingency plans have not been integrated, thus leading to failure in establishing priority areas. This investigation aimed to develop indices of environmental vulnerability to oil (IEVO), by combining information about environmental sensibility to oil and results of numerical modeling of spilled oil. To achieve that, a case study concerning to oil spills scenarios in a subtropical coastal area was designed, and IEVOs were calculated and presented in maps, in order to make the information about the areas’ vulnerability more easily visualized. For summer, the extension of coastline potentially affected by oil was approximately 150 km, and most of the coastline presented medium to high vulnerability. For winter, 230 km coastline would be affected, from which 75% were classified as medium to high vulnerability. Thus, IEVO maps allowed a rapid and clearer interpretation of the vulnerability of the mapped region, facilitating the planning process and the actions in response to an oil spill.     

弹性波阻抗在时移地震中的应用分析

常规声波波阻抗由于缺乏横波的信息,使得它对于流体的变化很不敏感,而弹性波阻抗反演,除了与纵波速度,密度有关之外,还与横波速度、入射角度有关,由于考虑了AVO(Amplitude Variation with Offset)效应,使得对流体和岩性的预测能力增强.本文将弹性波阻抗反演应用于时移地震中,采用AVA(Amplitude Variation with Angle)约束稀疏脉冲反演来进行弹性波阻抗反演,得到纵横波阻抗等弹性参数,通过交会图分析,将弹性参数转换为油藏参数,如预测含水饱和度的变化,通过本方法可增强对储层流体的预测,提高储层管理能力.     

Detection of point scatterers using diffraction imaging and deep learning

Diffracted waves carry high-resolution information that can help interpreting fine structural details at a scale smaller than the seismic wavelength. However, the diffraction energy tends to be weak compared to the reflected energy and is also sensitive to inaccuracies in the migration velocity, making the identification of its signal challenging. In this work, we present an innovative workflow to automatically detect scattering points in the migration dip angle domain using deep learning. By taking advantage of the different kinematic properties of reflected and diffracted waves, we separate the two types of signals by migrating the seismic amplitudes to dip angle gathers using prestack depth imaging in the local angle domain. Convolutional neural networks are a class of deep learning algorithms able to learn to extract spatial information about the data in order to identify its characteristics. They have now become the method of choice to solve supervised pattern recognition problems. In this work, we use wave equation modelling to create a large and diversified dataset of synthetic examples to train a network into identifying the probable position of scattering objects in the subsurface. After giving an intuitive introduction to diffraction imaging and deep learning and discussing some of the pitfalls of the methods, we evaluate the trained network on field data and demonstrate the validity and good generalization performance of our algorithm. We successfully identify with a high-accuracy and high-resolution diffraction points, including those which have a low signal to noise and reflection ratio. We also show how our method allows us to quickly scan through high dimensional data consisting of several versions of a dataset migrated with a range of velocities to overcome the strong effect of incorrect migration velocity on the diffraction signal.