Some Biases in Sampling Multilevel Piezometers for Volatile Organics

Multilevel piezometers are cost-effective monitoring devices for determining the three-dimensional distribution of solutes in ground water. Construction includes flexible tubing (plastic or Teflon®). Their sampling is subject to a number of'potential biases, particularly: (1) losses of volatile organic solutes via volatilization, (2) sorption onto the flexible tubing of the piezometers, (3) leaching of organics from this tubing, and (4) collection of unrepresentative samples due to inadequate piezometer flushing. It is shown that these biases are minimal or are easily controlled in most situations.
Another source of bias has been recognized. Organic solutes present in ground water above the screened level can penetrate the flexible plastic or Teflon tubing and contaminate the sampled water being drawn through this tubing. Laboratory tests and field results indicate this transmission causes low organic contaminant concentrations to be erroneously attributed to ground water which is free of such contaminants. The transmitted organics apparently desorb from the plastic tubing during flushing of even 40 piezometer volumes.
Recognition of this transmission problem provides for a better interpretation of existing organic contaminant distribution data. Caution is advised when considering the use of these monitoring devices in organic solute contaminant studies.     

Fraile Detection of Solar $$-Modes by Fossat et al.

The internal gravity modes of the Sun are notoriously difficult to detect, and the claimed detection of gravity modes presented by Fossat et al. (Astron. Astrophys.604, A40, 2017) is thus very exciting. Given the importance of these modes for understanding solar structure and dynamics, the results must be robust. While Fossat et al. described their method and parameter choices in detail, the sensitivity of their results to several parameters was not presented. Therefore, we test the sensitivity of the results to a selection of the parameters. The most concerning result is that the detection vanishes when we adjust the start time of the 16.5-year velocity time-series by a few hours. We conclude that this reported detection of gravity modes is extremely fragile and should be treated with utmost caution.     

Outgoing long wave radiation variability from IR satellite data prior to major earthquakes

Our analysis of the continuous outgoing long wave earth radiation (OLR) indicates anomalous variations prior to a number of medium to large earthquakes. The most recent analysis of OLR is from the M9.0 Sumatra Andaman Islands mega trust event. We compared the reference fields for December 2001 to 2004 and found OLR anomalous values, > 80 W/m², (2σ) within the epicentral area on Dec 21, 2004, 5 days before the event. We used the NOAA/IR daily (one degree) and monthly (two and half degree) gridded data to differentiate between the global and seasonal variability and the transient local anomalies. The cause of such anomalies is not fully understood; one possible explanation is the existence of thermal outgoing radiation as a result of near ground air ionization and latent heat change due to change of air humidity and temperature. This phenomenon is hypothesized to be part of a relationship between tectonic stresses, electrochemical and thermodynamic processes in the atmosphere and increasing mid IR flux, all part of a family of electromagnetic (EM) phenomena related to earthquake activity. The time scale of the observed variations is a few weeks before the onset of the seismic event. In comparison with several years of data, the observed time-series preceding the earthquake had unusually high OLR. The OLR anomaly corresponds to a large area of ground coverage and coincides with the main epicentral zone. The significance of these observations is explored using data from most recent East Asian earthquake swarm of December 2004 and three other earthquakes.     

Estimating uncertainty associated with water stages from a single SAR image

It is the goal of remote sensing to infer information about objects or a natural process from a remote location. This invokes that uncertainty in measurement should be viewed as central to remote sensing. In this study, the uncertainty associated with water stages derived from a single SAR image for the Alzette (G.D. of Luxembourg) 2003 flood is assessed using a stepped GLUE procedure. Main uncertain input factors to the SAR processing chain for estimating water stages include geolocation accuracy, spatial filter window size, image thresholding value, DEM vertical precision and the number of river cross sections at which water stages are estimated. Initial results show that even with plausible parameter values uncertainty in water stages over the entire river reach is 2.8 m on average. Adding spatially distributed field water stages to the GLUE analysis following a one-at-a-time approach helps to considerably reduce SAR water stage uncertainty (0.6 m on average) thereby identifying appropriate value ranges for each uncertain SAR water stage processing factor. For the GLUE analysis a Nash-like efficiency criterion adapted to spatial data is proposed whereby acceptable SAR model simulations are required to outperform a simpler regression model based on the field-surveyed average river bed gradient. Weighted CDFs for all factors based on the proposed efficiency criterion allow the generation of reliable uncertainty quantile ranges and 2D maps that show the uncertainty associated with SAR-derived water stages. The stepped GLUE procedure demonstrated that not all field data collected are necessary to achieve maximum constraining. A possible efficient way to decide on relevant locations at which to sample in the field is proposed. It is also suggested that the resulting uncertainty ranges and flood extent or depth maps may be used to evaluate 1D or 2D flood inundation models in terms of water stages, depths or extents. For this, the extended GLUE approach, which copes with the presence of uncertainty in the observed data, may be adopted.     

The ilmenite/titano-magnetite assemblage: kinetics of re-equilibration

Compositions of coexisting ilmenite and titano-magnetite commonly have been used to establish the conditions of formation for igneous and metamorphic rocks (Buddington-Lindsley geothermometer/oxybarometer). Temperatures so determined for some rocks are lower than might be expected for magmatic crystallization. These low temperatures are the result of re-equilibration of the Fe-Ti oxides during slow cooling. This process is one in which titano-magnetites are oxidized and ilmenites are reduced upon cooling along a buffer curve or in the presence of a fluid of constant composition. Textural studies indicate that the reduction of the ilmenite-hematite phase may be the slower reaction consequently, the rate-controlling step of the re-equilibration.The kinetics of the reduction reactions in the two-phase region between the joins for the ilmenite-hematite and magnetite-ulvöspinel solid solution series were studied by thermogravimetric analysis at controlledf_O2. With experiments conducted near the Ni/NiO and QFM buffers, annealed ilmenite(ss) synthetic charges were reduced isothermally. The temperature range for the experiments was 900–1250°C, with total gas-flow rates of 0.35 and 1.0 cm/s. Kinetic data were collected by measuring the weight change during each experiment as a function of time. Additionally, a cooling-rate (0.5°C/hr) experiment was conducted near the QFM buffer from 1150 to 985°C starting with an ilmenite-saturated spinel phase.At 1200°C and 0.35 cm/s flow rate, the reaction mechanism is an interfacial process where mass transport of the reacting gases is the rate-controlling step. In isothermal experiments at 1100°C and lower, diffusion of reactants through a solid product layer was the controlling mechanism. The activation energy for this diffusion process is 59±6 kcal/mole. Significant enhancement of diffusion occurs along grain boundaries.Although the individual phases may re-equilibrate internally upon cooling, the Fe-Ti contents must remain unchanged for proper use of this mineral indicator-i.e. no external re-equilibration. It is imperative that the kinetics of the situation be appreciated. These preliminary kinetic experiments indicate that extreme caution must be exercised when applying this ilmenite/titano-magnetite geothermometer/oxybarometer. These phases may only truly “quench” in the parameters of formation in the case of certain hypabyssal and extrusive occurrences     

利用卫星长波长地球物理位场对东亚地区深部结构的研究

有研究表明从中国北部经黄海、朝鲜半岛、日本海至日本本州岛存在一个大型断裂带,被称为“北纬40°断裂带”.为了验证该断裂带的存在性,本文通过GRACE和CHAMP卫星的位场数据中的大尺度区域重力和磁力异常来研究该区的大型构造特征.区域性研究范围为24°N~56°N,90°E~150°E,详细研究范围为32°N~42°N,122°E~132°E.2008年Taylor等利用CHAMP卫星2005年6月至12月收集的轨道数据,从测量数据中消除地核和外部产生磁场效应,经化磁极之后得到磁异常图,从中来追踪构造特征,并根据地质和地球物理资料建立了数学模型来进行解释.本文在此基础上,增加了GRACE卫星2003年10月份测量的重力异常场数据.实测差分位场数据利用高斯-勒让德求积法转换为垂向重力异常值,通过减去最新地球重力场模型（EGM96）的场值来去除波长超过1100 km的长波长异常,所得的重力异常数据和磁场数据进行比较.为了能从中得到研究区的共同构造特征,建立一个球面棱柱模型进行正演模拟、并进行重、磁异常场的波数相关性分析.研究发现在重、磁异常中均发现一个呈负相关的纬向构造特征,东亚地区“北纬40°断裂带”横穿了整个朝鲜半岛.     

Isolation of 12 microsatellite markers following a pyrosequencing procedure and cross-priming in two invasive cryptic species, Alexandrium catenella (group IV) and A. tamarense (group III) (Dinophyceae)

Alexandrium catenella (group IV) and Alexandrium tamarense (group III) (Dinophyceae) are two cryptic invasive phytoplankton species belonging to the A. tamarense species complex. Their worldwide spread is favored by the human activities, transportation and climate change. In order to describe their diversity in the Mediterranean Sea and understand their settlements and maintenances in this area, new microsatellite markers were developed based on Thau lagoon (France) samples of A. catenella and A. tamarense strains. In this study twelve new microsatellite markers are proposed. Five of these microsatellite markers show amplifications on A. tamarense and ten on A. catenella. Three of these 12 microsatellite markers allowed amplifications on both cryptic species. Finally, the haplotypic diversity ranged from 0.000 to 0.791 and 0.000 to 0.942 for A. catenella and A. tamarense respectively.     

Microstructural Inhomogeneity and Mechanical Anisotropy Associated with Bedding in Rothbach Sandstone

This study present the result of conventional triaxial tests conducted on samples of Rothbach sandstone cored parallel, oblique (at 45 degrees) and perpendicular to the bedding at effective pressures ranging from 5 to 250 MPa. Mechanical and microstructural data were used to determine the role of the bedding on mechanical strength and failure mode. We find that samples cored at 45 degrees to the bedding yield at intermediate level of differential stress between the ones for parallel and perpendicular samples at all effective pressures. Strain localization at high confining pressure (i.e., in the compactive domain) is observed in samples perpendicular and oblique to the bedding but not in samples cored parallel to the bedding. However, porosity reduction is comparable whether compactive shear bands, compaction bands or homogeneous cataclastic flow develop. Microstructural data suggest that (1) mechanical anisotropy is controlled by a preferred intergranular contact alignment parallel to the bedding and that (2) localization of compaction is controlled by bedding laminations and grain scale heterogeneity, which both prevent the development of well localized compaction features.     

Field Survey and Numerical Simulations: A Review of the 1998 Papua New Guinea Tsunami

— The Papua New Guinea (PNG) tsunami of 1998 is re-examined through a detailed review of the field survey as well as numerous numerical computations. The discussion of the field survey explores a number of possible misinterpretations of the recorded data. The survey data are then employed by a numerical model as a validation tool. A Boussinesq model and a nonlinear shallow water wave (NLSW) model are compared in order to quantify the effect of frequency dispersion on the landslide-generated tsunami. The numerical comparisons indicate that the NLSW model is a poor estimator of offshore wave heights. However, due to what appears to be depth-limited breaking seaward of Sissano spit, both numerical models are in agreement in the prediction of maximum water elevations at the overtopped spit. By comparing three different hot-start initial profiles of the tsunami wave, it is shown that the initial shape and orientation of the tsunami wave is secondary to the initial displaced water mass in regard to prediction of water elevations on the spit. These numerical results indicate that agreement between numerical prediction of runup values with field recorded values at PNG cannot be used to validate either a NLSW tsunami propagation model or a specific landslide tsunami hot-start initial condition. Finally, with the use of traditional tsunami codes, a new interpretation of the PNG runup measurements is presented.