On the Influence of Sea Surface Height Variability on Satellite Altimeter Derived Gravity

The effect of sea surface height (SSH) variability is one of the primary factors that limit the accuracy and resolution of altimeter-derived gravity values. We propose a method to estimate the influence of variation of the sea surface height on the accuracy of satellite-derived gravity by simulation technique, with a case study around Indonesian waters. Wederived an Indonesian marine gravity map using the Geosat-geodetic mission (GM). Since most of the area studied is located around coastal and shallow areas, the measurement of SSH of this area is less accurate. To obtain a distribution of SSH variability over the study area, Topex/Poseidon (T/P) data were first processed and assessed. Processing 52 cycles of the Topex/Poseidon data, the root mean square (RMS) of SSH variability for each cycle was found to vary from 1 to 179 cm. Further, for the purpose of estimating the accuracy of altimeter-derived gravity, we derived several levels of Gaussian noise, computed simulation data by adding the Gaussian noise to Geosat data, and determined simulated gravity maps. Based on the distribution of RMS values from T/P data and standard deviation (STD) differences between the simulated and the original gravity maps, we estimated the accuracy of the gravity map. Around Indonesian waters, the accuracy of the gravity map influenced by SSH variation was estimated to be within the range 0.8~93 mgal.     

Reaction-induced grain boundary cracking and anisotropic fluid flow during prograde devolatilization reactions within subduction zones

Devolatilization reactions during prograde metamorphism are a key control on the fluid distribution within subduction zones. Garnets in Mn-rich quartz schist within the Sanbagawa metamorphic belt of Japan are characterized by skeletal structures containing abundant quartz inclusions. Each quartz inclusion was angular-shaped, and showed random crystallographic orientations, suggesting that these quartz inclusions were trapped via grain boundary cracking during garnet growth. Such skeletal garnet within the quartz schist formed related to decarbonation reactions with a positive total volume change (?V _t > 0), whereas the euhedral garnet within the pelitic schists formed as a result of dehydration reaction with negative ?V _t values. Coupled hydrological–chemical–mechanical processes during metamorphic devolatilization reactions were investigated by a distinct element method (DEM) numerical simulation on a foliated rock that contained reactive minerals and non-reactive matrix minerals. Negative ?V _t reactions cause a decrease in fluid pressure and do not produce fractures within the matrix. In contrast, a fluid pressure increase by positive ?V _t reactions results in hydrofracturing of the matrix. This fracturing preferentially occurs along grain boundaries and causes episodic fluid pulses associated with the development of the fracture network. The precipitation of garnet within grain boundary fractures could explain the formation of the skeletal garnet. Our DEM model also suggests a strong influence of reaction-induced fracturing on anisotropic fluid flow, meaning that dominant fluid flow directions could easily change in response to changes in stress configuration and the magnitude of differential stress during prograde metamorphism within a subduction zone.     

Effects of Earth's curvature and radial heterogeneity in dislocation studies:Case studies of the 2008 Wenchuan earthquake and the 2004 Sumatra earthquake

Recently,effects of Earth's curvature and radial heterogeneity on coseismic deformations are often investigated based on the 2004 Sumatra earthquake.However,such effects are strongly related to earthquake types.As a low dip angle event,the 2004 Sumatra earthquake is not a good seismic case for such a topic since the effects for moderate dip angle events are much bigger.In this study,the half-space and spherical dislocation theories are used,respectively,to calculate coseismic displacements caused by the 2008 Wenchuan earthquake and the 2004 Sumatra earthquake.Effects of Earth's curvature and stratification are investigated through the discrepancies of results calculated using the two dislocation theories.Results show that the effects of Earth's curvature and stratification for the 2008 Wenchuan earthquake are much larger than those for the 2004 Sumatra earthquake.Ignoring the effects will cause errors up to 100%-200% in far field displacements for a moderate dip angle event like the 2008 Wenchuan earthquake.Such great effects are much bigger than those conclusions of previous studies.Besides,comparison with observations verifies that spherical dislocation theories yield better results than half-space ones in far fields.     

On causes and impacts of land subsidence in Bandung Basin, Indonesia

The Bandung Basin is a large intra-montane basin surrounded by volcanic highlands, in western Java, Indonesia, inhabited by more than seven million people. The basin, an area of about 2,300 km², is a highland plateau at approximately 650–700 m above sea level and is surrounded by up to 2,400 m high Late Tertiary and Quaternary volcanic terrain. Based on the results of nine GPS surveys conducted since 2000 up to 2011, it was shown that several locations in the Bandung Basin have experienced land subsidence, with an average rate of about ?8 cm/year and can go up to about ?23 cm/year in certain locations. A hypothesis has been proposed by several studies that land subsidence observed in several locations in the Bandung Basin has been caused mainly by excessive groundwater extraction. It is found that there is a strong correlation between the rates of groundwater level lowering with the GPS-derived rates of land subsidence in several locations in Bandung Basin. The GPS results in this study detected significant subsidence in the textile industry area, where very large volumes of groundwater are usually extracted. The impact of land subsidence in Bandung can be seen in several forms, mainly in the cracking and damage of houses, buildings and infrastructure. Land subsidence also aggravates the flooding in Bandung Basin, which has brought huge economic losses and deteriorated the quality of life and environment in the affected areas.     

Three-Dimensional Subsurface Structure Model of Kansai International Airport by Integration of Borehole Data and Seismic Profiles

The three-dimensional subsurface structure model around Kansai International Airport (KIX) has been developed based on the geologically and geotechnically investigated results of a large amount of borehole data for estimation of subsidence. The model consists of the alternation of marine clay and coarse deposits. The many seismic surveys and borehole drillings were carried out around the KIX before the constructing the airport. The stratigraphy of the model was renewed by the KIX18-1, which was about 1,300 m long drilling core and was drilled near the 2nd runway of the KIX from 2006 to 2007. In this study, the subsurface geological model was revised by integration of renewed borehole data and seismic profiles to incorporate complex warping structure and heterogenic lateral variation.     

Detection of white spot syndrome virus (WSSV) ofPenaeus chinensis by in situ hybridization

White Spot Syndrome Virus (WSSV) was purified from hemolymph of infected shrimp. After nucleic acid extraction from the purified virus particles, EcoR I-digested fragments of the WSSV genome were cloned; three of these fragments were used as non-radioactive probes labeled with DIG-11-dUTP. The probes hybridized in situ, with sections located in the nuclei of all WSSV-infected tissues. The virus was detected in the gill, stomach, epidermis, and connective tissue and so on, but not detected in healthy shrimp tissues and epithelial cells of hepatopancreatic tubules of diseased shrimp. This work was supported by SRF for ROCS, SEM and the National 863 Project (Grant 819-Q-08) and Project under major State Basic Research Development Program (Grant G1999012002).     

Development of a GPGPU-parallelized hybrid finite-discrete element method for modeling rock fracture

The hybrid finite-discrete element method (FDEM) is widely used for engineering applications, which, however, is computationally expensive and needs further development, especially when rock fracture process is modeled. This study aims to further develop a sequential hybrid FDEM code formerly proposed by the authors and parallelize it using compute unified device architecture (CUDA) C/C++ on the basis of a general-purpose graphics processing unit (GPGPU) for rock engineering applications. Because the contact detection algorithm in the sequential code is not suitable for GPGPU parallelization, a different contact detection algorithm is implemented in the GPGPU-parallelized hybrid FDEM. Moreover, a number of new features are implemented in the hybrid FDEM code, including the local damping technique for efficient geostatic stress analysis, contact damping, contact friction, and the absorbing boundary. Then, a number of simulations with both quasi-static and dynamic loading conditions are conducted using the GPGPU-parallelized hybrid FDEM, and the obtained results are compared both quantitatively and qualitatively with those from either theoretical analysis or the literature to calibrate the implementations. Finally, the speed-up performance of the hybrid FDEM is discussed in terms of its performance on various GPGPU accelerators and a comparison with the sequential code, which reveals that the GPGPU-parallelized hybrid FDEM can run more than 128 times faster than the sequential code if it is run on appropriate GPGPU accelerators, such as the Quadro GP100. It is concluded that the GPGPU-parallelized hybrid FDEM developed in this study is a valuable and powerful numerical tool for rock engineering applications.     

Coordination of water molecules with Na+ cations in a beryl channel as determined by polarized IR spectroscopy

Subtle variations of frequencies in the infrared (IR) absorption spectra of beryl have been predicted based on the coordination between extra-framework cations and water molecules in two orientations (referred to as type I and type II) trapped within the channel. In this study, the polarized IR spectra of hydrated synthetic beryl and natural beryl were measured to clarify the relationships between the frequencies of the absorption bands and the coordination states of type II water. Na⁺ was assumed to be the predominant cation coordinated to type II water in our samples, as determined by chemical analyses. These measurements revealed a clear quantitative linear relationship in absorbance between bands at 3,602 and 1,619 and at 3,589 and 1,631 cm⁻¹. On the basis of experimental and theoretical studies, we assigned these pairs of bands to the ν₁ and ν₂ modes of doubly coordinated type II H₂O and to singly coordinated type II H₂O, respectively. These assignments were supported by IR measurements of annealed natural beryl. We also conducted dehydration studies of natural beryl, in which two observed dehydration peaks, at 600 and 750°C, suggested the dehydration of type I and type II water, respectively.     

Improving the extreme rainfall forecast of Typhoon Morakot (2009) by assimilating radar data from Taiwan Island and mainland China

This study examined the impact of an improved initial field through assimilating ground-based radar data from mainland China and Taiwan Island to simulate the long-lasting and extreme rainfall caused by Morakot (2009). The vortex location and the subsequent track analyzed through the radial velocity data assimilation (VDA) are generally consistent with the best track. The initial humidity within the radar detecting region and Morakot’s northward translation speed can be significantly improved by the radar reflectivity data assimilation (ZDA). As a result, the heavy rainfall on both sides of Taiwan Strait can be reproduced with the joint application of VDA and ZDA. Based on sensitivity experiments, it was found that, without ZDA, the simulated storm underwent an unrealistic inward contraction after 12-h integration, due to underestimation of humidity in the global reanalysis, leading to underestimation of rainfall amount and coverage. Without the vortex relocation via VDA, the moister (drier) initial field with (without) ZDA will produce a more southward (northward) track, so that the rainfall location on both sides of Taiwan Strait will be affected. It was further found that the improvement in the humidity field of Morakot is mainly due to assimilation of high-value reflectivity (strong convection) observed by the radars in Taiwan Island, especially at Kenting station. By analysis of parcel trajectories and calculation of water vapor flux divergence, it was also found that the improved typhoon circulation through assimilating radar data can draw more water vapor from the environment during the subsequent simulation, eventually contributing to the extreme rainfall on both sides of Taiwan Strait.     

大气校正对SPOT卫星遥测水质的影响

藉由卫星遥测进行河川水质监测，目前尚没有较明确可行之方法，如何利用较为简单且适当的SPOT卫星遥测大气校正方法，正确辨识水体水质，是本研究的主要目的。利用SPOT卫星作两阶段非监督式及监督式自动分类，确认卫星影像中水质测站对应之水体样本，并将所有样本依季节分群，俾让卫星监测水体水质样本较为均质。模拟方式采用多变量回归、类神经网络及判别分析3种模式，并比较4种不同之大气校正程序。结果发现。以水质及其指标整体预测来看，类神经网络预测结果较优于多变量回归及判别分析的结果，大气校正方法以直接采用灰度值并消除最暗像元灰度值之校正方法，即可达到不错之预测结果。综合而言，以SPOT或分辨率更高之卫星光谱遥测水质是简单可行，但仍需更多数据以验证其精确度。