New buoy observation system for tsunami and crustal deformation

We have developed a new system for real-time observation of tsunamis and crustal deformation using a seafloor pressure sensor, an array of seafloor transponders and a Precise Point Positioning (PPP ) system on a buoy. The seafloor pressure sensor and the PPP system detect tsunamis, and the pressure sensor and the transponder array measure crustal deformation. The system is designed to be capable of detecting tsunami and vertical crustal deformation of ±8 m with a resolution of less than 5 mm. A noteworthy innovation in our system is its resistance to disturbance by strong ocean currents. Seismogenic zones near Japan lie in areas of strong currents like the Kuroshio, which reaches speeds of approximately 5.5 kt (2.8 m/s) around the Nankai Trough. Our techniques include slack mooring and new acoustic transmission methods using double pulses for sending tsunami data. The slack ratio can be specified for the environment of the deployment location. We can adjust slack ratios, rope lengths, anchor weights and buoy sizes to control the ability of the buoy system to maintain freeboard. The measured pressure data is converted to time difference of a double pulse and this simple method is effective to save battery to transmit data. The time difference of the double pulse has error due to move of the buoy and fluctuation of the seawater environment. We set a wire-end station 1,000 m beneath the buoy to minimize the error. The crustal deformation data is measured by acoustic ranging between the buoy and six transponders on the seafloor. All pressure and crustal deformation data are sent to land station in real-time using iridium communication.     

Near-field tsunami amplification factors in the Kii Peninsula,Japan for Dense Oceanfloor Network for Earthquakes and Tsunamis (DONET)

We investigated the correlation between coastal and offshore tsunami heights by using data from the Dense Oceanfloor Network for Earthquakes and Tsunamis (DONET) observational array of ocean-bottom pressure gauges in the Nankai trough off the Kii Peninsula, Japan. For near-field earthquakes, hydrostatic pressure changes may not accurately indicate sea surface fluctuations, because ocean-bottom pressure gauges are simultaneously displaced by crustal deformation due to faulting. To avoid this problem, we focused on the average waveform of the absolute value of the hydrostatic pressure changes recorded at all the DONET stations during a tsunami. We conducted a Monte Carlo tsunami simulation that revealed a clear relationship between the average waveforms of DONET and tsunami heights at the coast. This result indicates the possibility of accurate real-time prediction of tsunamis by use of arrays of ocean-bottom pressure gauges.     

Detectability analysis of interplate fault slips in the Nankai subduction thrust using seafloor observation instruments

To help the decision making regarding where to locate new observation instruments on the seafloor, we examined the detectability of interplate earthquakes and slow slips in the Nankai subduction thrust in Japan using seafloor observation instruments. Here, the detectability is defined as the smallest magnitude of the interplate fault slip detected by the assumed observation points based on crustal deformation simulation. In the detectability analyses, we considered the effect of sensor drifts that are particularly associated with seafloor observations. In addition, we introduced high-resolution three-dimensional (3D) finite element modeling of crustal deformation to consider the effect of the topography and 3D heterogeneous crustal structure around the Nankai Trough. The results of the detectability analyses show that introducing new seafloor stations for tilt observation in the Nankai region should increase the detectability of small- or medium-sized interplate earthquakes and slow slips significantly. Based on the obtained results, we also discuss the advantage of both the existing and the new observation instruments in detecting interplate fault slips.

     

Improvement of earthquake locations with the Marine Cable Hosted Observatory (MACHO) offshore NE Taiwan

In order to improve the locating capability for offshore earthquakes and tsunamis monitored off northeastern Taiwan, a cable-based ocean bottom seismographic observatory named “Marine Cable Hosted Observatory” (MACHO) was constructed and began operation at the end of 2011. The installed instruments of the observatory include a broadband seismometer, a strong-motion seismometer and a pressure gauge. In addition, various scientific instruments could be deployed for other purposes as well. At present, the seismic data are transmitted in real-time via a fiber cable, and integrated into the current inland seismographic network in Taiwan. The ocean bottom station has contributed to provide high quality seismic data already. According to observations from January 2012 to June 2013, there were a total of 15,168 earthquakes recorded by the system. By using the data from the ocean bottom station, the number of relocated earthquakes with an azimuth gap less than 180 degrees substantially increase about 34 %. Meanwhile, the root–mean–square of the time residual, the error in epicenter, and the error in depth of the earthquake locations decrease. Therefore, the implementation of MACHO has the advantage of extending the coverage of existing the Taiwan seismic network to the offshore, providing more accurate and real-time seismic data for offshore earthquakes monitoring. The results show that MACHO is crucial and necessary for monitoring seismic activities in northeastern Taiwan.     

Assimilation of simulated float data in Lagrangian coordinates

We implement an approach for the accurate assimilation of Lagrangian data into regional general ocean circulation models. The forward model is expressed in Lagrangian coordinates and simulated float data are incorporated into the model via four-dimensional variational data assimilation. We show that forward solutions computed in Lagrangian coordinates are reliable for time periods of up to 100 days with phase speeds of 1 m/s and deformation radius of 35 km. The position and depth of simulated floats are assimilated into the viscous, Lagrangian shallow water equations. The weights for the errors in the model and data are varied and the assimilation results react appropriately. We show the effect of different spatial and temporal samplings of float data on all Lagrangian trajectories in the computational domain. At the end of the assimilation period, results from the Lagrangian shallow water equations could be interpolated and used as initial and boundary conditions in an Eulerian general ocean circulation model.     

The impact of mean dynamic topography on a sea-level anomaly assimilation in the South China Sea based on an eddy-resolving model

The sea-level anomaly (SLA) from a satellite altimeter has a high accuracy and can be used to improve ocean state estimation by assimilation techniques. However, the lack of an accurate mean dynamic topography (MDT) is still a bothersome issue in an ocean data assimilation. The previous studies showed that the errors in MDT have significant impacts on assimilation results, especially on the time-mean components of ocean states and on the time variant parts of states via nonlinear ocean dynamics. The temporal-spatial differences of three MDTs and their impacts on the SLA analysis are focused on in the South China Sea (SCS). The theoretical analysis shows that even for linear models, the errors in MDT have impacts on the SLA analysis using a sequential data assimilation scheme. Assimilation experiments, based on EnOI scheme and HYCOM, with three MDTs from July 2003 to June 2004 also show that the SLA assimilation is very sensitive to the choice of different MDTs in the SCS with obvious differences between the experimental results and observations in the centre of the SCS and in the vicinity of the Philippine Islands. A new MDT for assimilation of SLA data in the SCS was proposed. The results from the assimilation experiment with this new MDT show a marked reduction (increase) in the RMSEs (correlation coefficient) between the experimental and observed SLA. Furthermore, the subsurface temperature field is also improved with this new MDT in the SCS.     

Altimeter Signal-to-Noise for Deep Ocean Processes in Operational Systems

The ocean signal for this study is the sea surface height due to the slowly varying (greater than 5-day) ocean processes, which are predominantly the deep ocean mesoscale. These processes are the focus of present assimilation systems for monitoring and predicting ocean circulation due to ocean fronts and eddies and the associated environmental changes that impact real time activities in areas with depths greater than about 200 m. By this definition, signal-to-noise may be estimated directly from altimeter data sets through a crossover point analysis. The RMS variability in crossover differences is due to instrument noise, errors in environmental corrections to the satellite observation, and short time period oceanic variations. The signal-to-noise ratio indicates that shallow areas are typically not well observed due to the high frequency fluctuations. Many deep ocean areas also contain significant high frequency variability such as the subpolar latitudes, which have large atmospheric pressure systems moving through, and these in turn generate large errors in the inverse barometer correction. Understanding the spatial variations of signal to noise is a necessary prerequisite for correct assimilation of the data into operational systems.     

Altimeter Signal-to-Noise for Deep Ocean Processes in Operational Systems

The ocean signal for this study is the sea surface height due to the slowly varying (greater than 5-day) ocean processes, which are predominantly the deep ocean mesoscale. These processes are the focus of present assimilation systems for monitoring and predicting ocean circulation due to ocean fronts and eddies and the associated environmental changes that impact real time activities in areas with depths greater than about 200 m. By this definition, signal-to-noise may be estimated directly from altimeter data sets through a crossover point analysis. The RMS variability in crossover differences is due to instrument noise, errors in environmental corrections to the satellite observation, and short time period oceanic variations. The signal-to-noise ratio indicates that shallow areas are typically not well observed due to the high frequency fluctuations. Many deep ocean areas also contain significant high frequency variability such as the subpolar latitudes, which have large atmospheric pressure systems moving through, and these in turn generate large errors in the inverse barometer correction. Understanding the spatial variations of signal to noise is a necessary prerequisite for correct assimilation of the data into operational systems.     

underwThe Hawaii-2 Observatory seismic system

The Hawaii-2 Observatory seismic system is currently transmitting high-quality seismic data from the ocean floor in the central NE Pacific Ocean through Hawaii to the IRIS Data Management Center. The system includes broad-band seismic, geophone, acoustic, and ocean current sensors. The seismic sensors are buried about 0.4 m below the ocean floor to improve coupling to the ocean bottom and to reduce noise levels. The system can be remotely calibrated, leveled and locked, and gains can be changed on command from shore. Data are temporarily stored in the seismic package for retransmission as needed to correct for transmission problems and to prevent loss of data. Data generated are valuable for studies of the Earth's structure and the dynamics of earthquakes     

The estimation of ocean bottom seismographs' coupling characteristics by means of microshock recordings

A transient technique was used for estimating the bottom-instrument response function in experiments with bottom seismographs (OBS) in deep ocean. The sharp mechanical impacts on bottom seismograph usually induced by bottom displacements under the instrument (microshocks) are suggested as rough analogues of the impulses for the bottom-instrument system transient calibration. It has been found that microshocks has usually sufficiently small duration to be used for coupling characteristics estimates. Test measurements have shown that in many cases this method makes it possible to distinguish spectral features characterising earthquakes and seismic noise wavetrains from those which are caused by coupling resonances of the OBS-sediment inter-face.     

Ocean Tides Observed from A GPS Receiver on Floating Sea Ice Near Chinese Zhongshan Station,Antarctica

Due to limit of coverage in TOPEX/Poseidon (T/P) satellite and sparseness of in-situ tide gauges around Antarctica, the accuracy of global ocean tide models in Antarctic seas is relatively poorer than in low- and mid-latitude regions. To better understand ocean tides in Prydz Bay, east Antarctica, a GPS receiver was deployed on floating sea ice to measure tide-induced ice motion in multiple campaigns. Four online Precise Point Positioning (PPP) services are used to process the GPS data in the kinematic PPP mode, and UTide software is used to separate the major tidal constituents. Comparison between results from different processing methods (relative processing solutions from Track, kinematic PPP solutions from online services) and with bottom pressure gauge (BPG) shows that, high-accuracy tidal information can be obtained from GPS observations on floating sea ice, the root-sum-square (RSS) for the eight major constituents (O1, K1, P1, Q1, M2, S2, N2, K2) is below 4 cm. We have also studied the impacts of data span and filter edge effects at daily boundaries on the accuracy of tide estimates, and found that to obtain reliable tide estimates and neglect the filter edge effects, continuous observation longer than 30 days is necessary. Our study suggests that GPS provides an independent method to estimate tides in Prydz Bay, and can be an alternative to tidal gauges, which are costly and hard to maintain in Antarctica.     

Measurement of tsunamis in the open ocean

Abstract

The spectrum of long waves, which are a background to tsunamis, is analyzed on the basis of records of near‐bottom pressure sensors obtained in the Northwest Pacific during the first and second USA‐USSR expeditions on the investigation of tsunamis in the open ocean (1975 and 1978). Instrumental trends, tidal oscillations, and quasistationary longwave noise were contained in the records. Special filters were used to pick out small waves generated by the seismicity of the ocean bottom. A decrease of noise level from 10² cm (including tides) to 10°¹ cm could be reached. The level of long‐wave noise is found to depend on the instrument's location. Minimal disturbances in the records were observed at stations situated on the edge of the continental slope. The influence of cyclones passing over the instrument's site is deduced. It shows an increase in noise level on Nyquist frequency (0.5 min°¹), which possibly is caused by the action of swell generated by the cyclone. Seismicity of the region under investigation for the second expedition (August‐October 1978) is described, and the recurrence of tsunamis is estimated, including microtsunamis. Taking into account this analysis, records were filtered and sections corresponding to probable arrivals of tsunamis from most strong earthquakes were selected. The anomalous disturbance of ocean level with a height of about 0.5 cm was found. Presumably, it was generated by an earthquake with magnitude M = 6.     

基于背景噪声互相关技术的海底地震仪的时钟校正方法在西南印度洋脊的应用

Seismic monitoring using ocean bottom seismometers(OBS) is an efficient method for investigating earthquakes in mid-ocean ridge far away from land. Clock synchronization among the OBSs is difficult without direct communication because electromagnetic signals cannot propagate efficiently in water. Time correction can be estimated through global positioning system(GPS) synchronization if clock drift is linear before and after the deployment. However, some OBSs in the experiments at the southwest Indian ridge(SWIR) on the Chinese DY125-34 cruise had not been re-synchronized from GPS after recovery. So we attempted to estimate clock drift between each station pairs using time symmetry analysis(TSA) based on ambient noise cross-correlation. We tested the feasibility of the TSA method by analyzing daily noise cross-correlation functions(NCFs) that extract from the data of another OBS experiment on the Chinese DY125-40 cruise with known clock drift and the same deployment site. The results suggest that the NCFs' travel time of surface wave between any two stations are symmetrical and have an opposite growing direction with the date. The influence of different band-pass filters,different components and different normalized methods was discussed. The TSA method appeared to be optimal for the hydrophone data within the period band of 2–5 s in dozens of km-scale interstation distances. A significant clock drift of ~2 s was estimated between OBSs sets through linear regression during a 108-d deployment on the Chinese cruise DY125-34. Time correction of the OBS by the ambient noise cross-correlation was demonstrated as a practical approach with the appropriate parameters in case of no GPS re-synchronization.     

集合最优平滑同化方法的研究与应用

海洋状态场的历史变化过程对其分布状态有重要影响。在观测资料稀疏的情况下,合理利用历史观测资料能够为海洋数据同化提供大量有效信息。然而在目前的顺序资料同化过程中,往往只同化当前时刻的观测数据,没有考虑到历史观测资料对当前状态的约束。四维变分虽然可以体现变量在时间维度的演变过程,但引入伴随方程会增加计算代价。本文基于集合最优平滑同化算法（Ensemble Optimal Smoothing,EnOS）探讨了一种在数据同化中加入历史观测资料的简易可行方案,其能够根据历史观测数据估计当前状态,并进行单点同化实验和区域同化实验来验证该方案的有效性。实验结果表明,将历史观测资料引入到同化过程中可以把控时间演变趋势,减小分析数据与真实值之间的偏差,更有效地消除数值模式误差,提高同化质量。     

Comparison of Sea Surface Heights Derived from Satellite Altimetry and from Ocean Bottom Pressure Gauges: The SW Pacific MOTEVAS Project

A bottom pressure gauge (BPG) was installed in proximity (3.7 km at closest approach) of Jason-1 and formerly TOPEX/Poseidon (T/P) ground track No. 238 at the Wusi site, located ∼ 10 km offshore off the west coast of Santo Island, Vanuatu, Southwest (SW) Pacific. Sea level variations are inferred from the bottom pressure, seawater temperature, and salinity, corrected for the measured surface atmospheric pressure. The expansion of the water column (steric increase in sea surface height, SSH) due to temperature and salinity changes is approximated by the equation of state. We compare time series of SSH derived from T/P Side B altimeter Geophysical Data Records (GDR) and Jason-1 Interim Geophysical Data Records (IGDR), with the gauge-inferred sea level variations. Since altimeter SSH is a geocentric measurement, whereas the gauge-inferred observation is a relative sea level measurement, SSH comparison is conducted with the means of both series removed in this study. In addition, high-rate (1-Hz) bottom pressure implied wave heights (H_1/3) are compared with the significant wave height (SWH) measured by Jason-1. Noticeable discrepancy is found in this comparison for high waves, however the differences do not contribute significantly to the difference in sea level variations observed between the altimeter and the pressure gauge. In situ atmospheric pressure measurements are also used to verify the inverse barometer (IB) and the dry troposphere corrections (DTC) used in the Jason IGDR. We observe a bias between the IGDR corrections and those derived from the local sensors. Standard deviations of the sea level differences between T/P and BPG is 52 mm and is 48 mm between Jason and BPG, indicating that both altimeters have similar performance at the Wusi site and that it is feasible to conduct long-term monitoring of altimetry at such a site.     

The Ocean Loading Effect in the GPS Analysis: A Case Study in the Antarctic Peninsula Region

The accuracy of GPS data analysis for the vertical component has reached a level where the vertical crustal deformation due to the ocean tidal loading should be taken into consideration. Ocean loading affects in particular the results of the GPS analysis for those observation sets covering less than 24 hours. Especially in these cases, a correction for this phenomenon should be performed during the data analysis. On the other hand, it is possible to estimate the ocean loading effect from GPS data sets. In this way it is possible to validate models for the loading deformation derived from global ocean tide models.     

Reproducibility of spatial and temporal distribution of aseismic slips in Hyuga-nada of southwest Japan

The Hyuga-nada region of southwest Japan, which is located off the east coast of Kyushu Island, may have the potential to generate great interplate earthquakes along the Nankai trough in the future. In this area, thrust earthquakes of M = 6.7–7.2 have occurred with recurrence intervals of approximately 30 years. In association with these earthquakes, possible local heterogeneities of plate coupling may be expected within 100 km from the coast in the Hyuga-nada region. We investigate numerical experiments to determine the spatial and temporal resolution of slip on the plate interface beneath the Hyuga-nada offshore region. For this purpose, we calculated synthetic displacement data from the result of numerical simulation conducted for the afterslip following an M_w 6.8 earthquake, for existing global positioning system stations on land and planned ocean floor seismic network stations. The spatial and temporal distribution of fault slip is then estimated using a Kalman filter-based inversion. The slip distribution estimated by using ocean floor stations demonstrates that the heterogeneity of plate coupling is resolved approximately within 50 km from the coastal area. This heterogeneity corresponds to the coseismic area of an M_w 6.8 earthquake with a radius of 10 km. Our study quantitatively evaluates the spatial resolution of aseismic slip in the Hyuga-nada region. Analysis based on continuous ocean floor data is useful for resolving the spatial variations of heterogeneities in plate couplings.     

An ocean data assimilation system and reanalysis of the World Ocean hydrophysical fields

A new version of the ocean data assimilation system (ODAS) developed at the Hydrometcentre of Russia is presented. The assimilation is performed following the sequential scheme analysis–forecast–analysis. The main components of the ODAS are procedures for operational observation data processing, a variational analysis scheme, and an ocean general circulation model used to estimate the first guess fields involved in the analysis. In situ observations of temperature and salinity in the upper 1400-m ocean layer obtained from various observational platforms are used as input data. In the new ODAS version, the horizontal resolution of the assimilating model and of the output products is increased, the previous 2D-Var analysis scheme is replaced by a more general 3D-Var scheme, and a more flexible incremental analysis updating procedure is introduced to correct the model calculations. A reanalysis of the main World Ocean hydrophysical fields over the 2005–2015 period has been performed using the updated ODAS. The reanalysis results are compared with data from independent sources.     

Large aperture seismic imaging at a convergent margin: Techniques and results from the Costa Rica seismogenic zone

In March and April 1995 a cooperative German, Costa Rican, and United States research team recorded onshore-offshore seismic data sets along the Pacific margin of Costa Rica using the R/V Ewing. Off the Nicoya Peninsula we used a linear array of ocean bottom seismometers and hydrophones (OBS/H) with onshore seismometers extending across much of the isthmus. In the central area we deployed an OBS/H areal array consisting of 30 instruments over a 9 km by 35-km area and had land stations on the Nicoya Peninsula adjacent to this marine array and also extending northeast on the main Costa Rican landmass. Our goal in these experiments was to determine the crustal velocity structure along different portions of this convergent margin and to use the dense instrument deployments to create migrated reflection images of the plate boundary zone and the subducting Cocos Plate. Our specific goal in the central area was to determine whether a subducted seamount is present at the location of the 1990, M 7 earthquake off the Nicoya Peninsula and can thus be linked to its nucleation. Subsequently we have processed the data to improve reflection signals, used the data to calculate crustal velocity models, and developed several wide-aperture migration techniques, based on a Kirchhoff algorithm, to produce reflection images. Along the northern transect we used the ocean bottom data to construct a detailed crustal velocity model, but reflections from the plate boundary and top and bottom of the subducting Cocos plate are difficult to identify and have so far produced poor images. In contrast, the land stations along this same transect recorded clear reflections from the top of the subducting plate or plate boundary, within the seismogenic zone, and we have constructed a clear image from this reflector beneath the Nicoya shelf. Data from the 3-D seismic experiment suffer from high-amplitude, coherent noise (arrivals other than reflections), and we have tried many techniques to enhance the signal to noise ratio of reflected arrivals. Due to the noise, an apparent lack of strong reflections from the plate boundary zone, and probable structural complexity, the resulting 3-D images only poorly resolve the top of the subducting Cocos Plate. The images are not able to provide compelling evidence of whether there is a subducting seamount at the 1990 earthquake hypocenter. Our results do show that OBS surveys are capable of creating images of the plate boundary zone and the subducting plate well into the seismogenic zone if coherent reflections are recorded at 1.8 km instrument spacing (2-D) and 5 km inline by 1 km crossline spacing for 3-D acquisition. However, due to typical high amplitude coherent noise, imaging results may be poorer than expected, especially in unfavorable geologic settings such as our 3-D survey area. More effective noise reduction in acquisition, possibly with the use of vertical hydrophone arrays, and in processing, with advanced multiple removal and possibly depth filtering, is required to achieve the desired detailed images of the seismogenic plate boundary zone.     

Mitobs: A seismometer system for ocean-bottom earthquake studies

The MIT ocean-bottom seismometer is a free-fall, pop-up instrument capable of recording three components of seismic data on the sea floor for periods of at least one month. Data are recorded in digital format on a specially designed magnetic tape recorder. An event recording scheme and semiconductor memories assure both efficient data storage and preservation of first motion information. Sensors and recording electronics are housed in a cylindrical pressure vessel, which sits vertically atop an expendable base plate on the ocean bottom. Attached to the pressure case are three glass spheres for buoyancy. After a pre-set time interval, a motor-driven mechanical latch release frees the instrument to float to the ocean surface for recovery.