Heterogeneous surface displacement pattern at the Hatchobaru geothermal field inferred from SAR interferometry time-series

We estimated surface displacements using persistent scatterer SAR interferometry (PS-InSAR) around the Hatchobaru geothermal field, Japan, from 18 ALOS/PALSAR images acquired from July 2007 to December 2010. Generally, geothermal fields, covered with natural targets such as rocky terrain and vegetation, have been one of the difficult targets for PS-InSAR analysis. However, we applied space adaptive filtering to increase the number of pixels for measuring surface displacement. The results of our analysis demonstrate ground subsidence with decaying velocity over the observation period around the geothermal field. The spatial pattern of ground subsidence includes sharp boundaries of subsidence that can be interpreted as fault traces. We demonstrated the usefulness of PS-InSAR analysis with the space adaptive filtering to estimate surface displacements with high spatial resolution and high spatial density around a geothermal field.     

Classification of forest land attributes using multi-source remotely sensed data

The aim of the study was to (1) examine the classification of forest land using airborne laser scanning (ALS) data, satellite images and sample plots of the Finnish National Forest Inventory (NFI) as training data and to (2) identify best performing metrics for classifying forest land attributes. Six different schemes of forest land classification were studied: land use/land cover (LU/LC) classification using both national classes and FAO (Food and Agricultural Organization of the United Nations) classes, main type, site type, peat land type and drainage status. Special interest was to test different ALS-based surface metrics in classification of forest land attributes. Field data consisted of 828 NFI plots collected in 2008–2012 in southern Finland and remotely sensed data was from summer 2010. Multinomial logistic regression was used as the classification method. Classification of LU/LC classes were highly accurate (kappa-values 0.90 and 0.91) but also the classification of site type, peat land type and drainage status succeeded moderately well (kappa-values 0.51, 0.69 and 0.52). ALS-based surface metrics were found to be the most important predictor variables in classification of LU/LC class, main type and drainage status. In best classification models of forest site types both spectral metrics from satellite data and point cloud metrics from ALS were used. In turn, in the classification of peat land types ALS point cloud metrics played the most important role. Results indicated that the prediction of site type and forest land category could be incorporated into stand level forest management inventory system in Finland.     

Impact of sub-pixel heterogeneity on modelled brightness temperature for an agricultural region

Knowledge of sub-pixel heterogeneity, particularly at the passive microwave scale, can improve the brightness temperature (and ultimately the soil moisture) estimation. However, the impact of surface heterogeneity (in terms of soil moisture, soil temperature and vegetation water content) on brightness temperature in an agricultural setting is relatively unknown. The Soil Moisture Active Passive Validation Experiment 2012 (SMAPVEX12) provided an opportunity to evaluate sub-pixel heterogeneity at the scale of a Soil Moisture Ocean Salinity (SMOS) or the Soil Moisture Active Passive (SMAP) radiometer footprint using field measured data. The first objective of this study was to determine if accounting for surface heterogeneity reduced the error between estimated brightness temperature (T_b) and T_b measured by SMOS. It was found that when accounting for variation in surface soil moisture, temperature and vegetation water content within the pixel footprint, the error between the modelled T_b and the measured T_b was less than if a homogeneous pixel were modelled. The correlation between the surface parameters and the error associated with not accounting for surface heterogeneity were investigated. It was found that there was low to moderate correlation between the error and the coefficient of variance associated with the measured soil moisture, soil temperature and vegetation volumetric water content during the field campaign. However, it was found that the correlations changed depending on the stage of vegetation growth and the amount of time following a precipitation event. At the start of the field campaign (following a precipitation event), there was strong correlation between the error and all three surface parameters (r ≥ 0.75). Following a precipitation event close to the middle of the field campaign (during which there was rapid growth in vegetation), there was strong correlation between the error and the variability in vegetation water content (r = 0.89), moderate correlation with soil moisture (r = 0.61) and low correlation with soil temperature (r = 0.26).     

扬子西缘乡城-丽江结合带燕山期斑岩Mo多金属矿床成矿系统

位于扬子陆块西缘与西南三江造山带结合部位的乡城-丽江地区,北起四川乡城经云南格咱南至丽江地区,在印支期斑岩铜矿带上新发现叠加了燕山期斑岩Mo多金属成矿作用,形成乡城-丽江斑岩Mo矿带。受印支期古特提斯洋盆闭合后地壳缩短与加厚的影响,燕山晚期下地壳发生拆沉作用,导致后碰撞型花岗质岩浆侵位,发育了Mo多金属成矿作用,并构成斑岩成矿系统。岩石地球化学特征表明,燕山晚期含矿斑岩具高硅(SiO_2=66.29%~79.36%)、高碱(K_2O+Na_2O=5.07%~9.24%)、富钾(K_2O/Na_2O=0.71~2.13)的特点,属于高钾钙碱性岩石系列。岩石富集轻稀土元素(LREE),具有负δEu异常;微量元素具富集大离子亲石元素K、U、Th、Rb,亏损高场强元素Nb、Ta、P、Ti的特征;大离子亲石元素的富集和Nb、Ta等高场强元素的亏损,表明形成这些岩体的岩浆主要来自地壳,具有造山带花岗岩的地球化学特征。区内代表性含矿斑岩相似的地球化学组成及分布特征表明,乡城-丽江结合带燕山晚期的成矿斑岩是由同源岩浆分异演化而来。成矿系统的分析表明,燕山晚期主要含矿斑岩成矿物质以壳源为主,但具有少量地幔物质的加入,成岩成矿物质来源具有相似或一致的源区特征。在成矿元素的组合上,由成矿斑岩体向外带表现出W、Mo→W、Mo、Cu→Pb、Zn、Ag的演化趋势。研究表明,本区燕山晚期花岗岩浆的侵入及Mo多金属成矿作用并不是独立或个别的成矿事件,而是纵跨义敦岛弧、甘孜-理塘结合带及扬子西缘的带状成矿活动,属于区内与燕山晚期岩浆侵入作用相关的统一斑岩成矿系统。     

Development of a new method for efficiently calculating of evaporation from the phreatic aquifer in variably saturated flow modeling

In this paper, we proposed a new method that has been developed based on the surface soil moisture content (SSMC) to more efficiently calculate the groundwater evaporation in variably saturated flow modeling. In this method, the empirical formula to calculate evaporation was modified and the value of the formula varies from zero to one as a closed interval. In addition, the simulation code for calculating the groundwater evaporation based on the SSMC method was incorporated into the EOS9 module of Tough2, a variably saturated flow modeling code. Finally, two numerical tests and a case simulation were conducted to verify the feasibility and accuracy of the SSMC method. Simulation results indicate that the SSMC method is capable of appropriately simulating the characteristics of water flow in vadose zone and the amount of evaporation with the variable water table. And such results are in coincidence with the value calculated by the logistic function method, and fit well with the measured data globally rather than locally.     

Molybdenum Mass Fractions and Isotopic Compositions of International Geological Reference Materials

The double‐spike method with multi‐collector inductively coupled plasma‐mass spectrometry was used to measure the Mo mass fractions and isotopic compositions of a set of geological reference materials including the mineral molybdenite, seawater, coral, as well as igneous and sedimentary rocks. The long‐term reproducibility of the Mo isotopic measurements, based on two‐year analyses of NIST SRM 3134 reference solutions and seawater samples, was ≤ 0.07‰ (two standard deviations, 2s, n = 167) for δ^98/95Mo. Accuracy was evaluated by analyses of Atlantic seawater, which yielded a mean δ^98/95Mo of 2.03 ± 0.06‰ (2s, n = 30, relative to NIST SRM 3134 = 0‰) and mass fraction of 0.0104 ± 0.0006 μg g^?1 (2s, n = 30), which is indistinguishable from seawater samples taken world‐wide and measured in other laboratories. The comprehensive data set presented in this study serves as a reference for quality assurance and interlaboratory comparison of high‐precision Mo mass fractions and isotopic compositions.     

High Precision Sr‐Nd‐Hf‐Pb Isotopic Compositions of USGS Reference Material BCR‐2

The Lamont‐Doherty Earth Observatory radiogenic isotope group has been systematically measuring Sr‐Nd‐Pb‐Hf isotopes of USGS reference material BCR‐2 (Columbia River Basalt 2), as a chemical processing and instrumental quality control monitor for isotopic measurements. BCR‐2 is now a widely used geochemical inter‐laboratory reference material (RM), with its predecessor BCR‐1 no longer available. Recognising that precise and accurate data on RMs is important for ensuring analytical quality and for comparing data between different laboratories, we present a compilation of multiple digestions and analyses made on BCR‐2 during the first author's dissertation research. The best estimates of Sr, Nd and Hf isotope ratios and measurement reproducibilities, after filtering at the 2s level for outliers, were ⁸⁷Sr/⁸⁶Sr = 0.705000 ± 11 (2s, 16 ppm, n = 21, sixteen digestions, one outlier), ¹⁴³Nd/¹⁴⁴Nd = 0.512637 ± 13 (2s, 25 ppm, n = 27, thirteen digestions, one outlier) and ¹⁷⁶Hf/¹⁷⁷Hf = 0.282866 ± 11 (2s, 39 ppm, n = 25, thirteen digestions, no outliers). Mean Nd and Hf values were within error of those reported by Weis et al. (2006, 2007) in their studies of RMs; mean Sr values were just outside the 2s uncertainty range of both laboratories. Moreover, a survey of published Sr‐Nd‐Hf data shows that our results fall within the range of reported values, but with a smaller variability. Our Pb isotope results on acid leached BCR‐2 aliquots (n = 26, twelve digestions, two outliers) were ²⁰⁶Pb/²⁰⁴Pb = 18.8029 ± 10 (2s, 55 ppm), ²⁰⁷Pb/²⁰⁴Pb = 15.6239 ± 8 (2s, 52 ppm), ²⁰⁸Pb/²⁰⁴Pb = 38.8287 ± 25 (2s, 63 ppm). We confirm that unleached BCR‐2 powder is contaminated with Pb, and that sufficient leaching prior to digestion is required to achieve accurate values for the uncontaminated Pb isotopic compositions.     

An Investigation of Digestion Methods for Trace Elements in Bauxite and Their Determination in Ten Bauxite Reference Materials Using Inductively Coupled Plasma‐Mass Spectrometry

Trace elements from samples of bauxite deposits can provide useful information relevant to the exploration of the ore‐forming process. Sample digestion is a fundamental and critical stage in the process of geochemical analysis, which enables the acquisition of accurate trace element data by ICP‐MS. However, the conventional bomb digestion method with HF/HNO₃ results in a significant loss of rare earth elements (REEs) due to the formation of insoluble AlF₃ precipitates during the digestion of bauxite samples. In this study, the digestion capability of the following methods was investigated: (a) ‘Mg‐addition’ bomb digestion, (b) NH₄HF₂ open vessel digestion and (c) NH₄F open vessel digestion. ‘Mg‐addition’ bomb digestion can effectively suppress the formation of AlF₃ and simultaneously ensure the complete decomposition of resistant minerals in bauxite samples. The addition of MgO to the bauxite samples resulted in (Mg + Ca)/Al ratios ≥ 1. However, adding a large amount of MgO leads to significant blank contamination for some transition elements (V, Cr, Ni and Zn). The NH₄HF₂ or NH₄F open vessel digestion methods can also completely digest resistant minerals in bauxite samples in a short period of time (5 hr). Unlike conventional bomb digestion with HF/HNO₃, the white precipitates and the semi‐transparent gels present in the NH₄HF₂ and NH₄F digestion methods could be efficiently dissolved by evaporation with HClO₄. Based on these three optimised digestion methods, thirty‐seven trace elements including REEs in ten bauxite reference materials (RMs) were determined by ICP‐MS. The data obtained showed excellent inter‐method reproducibility (agreement within 5% for REEs). The relative standard deviation (% RSD) for most elements was < 6%. The concentrations of trace elements in the ten bauxite RMs showed agreement with the limited certified (Li, V, Cr, Cu, Zn, Ga, Sr, Zr and Pb) and information values (Co, Ba, Ce and Hf) available. New trace element data for the ten RMs are provided, some of which for the first time.     

Production and Certification of Synthetic Selenium Isotopic Reference Materials

Due to intensive research into selenium isotopes in recent years, the increasing requirement for reliable and comparable measurement results has created a strong demand for selenium isotopic certified reference materials (iCRM) that were previously not available. To address this, eleven selenium iCRMs were developed, including ten synthetic iCRMs (GBW 04447–GBW 04456) and one natural iCRM (GBW 04457). The synthetic iCRMs were prepared with ⁷⁶Se, ⁷⁸Se, ⁸⁰Se and ⁸²Se solutions and a natural selenium solution; the natural iCRM was prepared with highly pure selenium material. The property values of isotope ratios in these iCRMs were certified by calibrated mass spectrometry with a collision cell multi‐collector ICP‐MS. The mass discrimination effect of the instrument was corrected with corresponding ⁷⁸Se/⁷⁶Se isotope mixtures and ⁸²Se/⁷⁶Se isotope mixtures, which were gravimetrically prepared with purified, isotopically enriched selenium materials. Homogeneity and stability tests were performed, and no significant influences were found. The uncertainty of the property values of the iCRMs was evaluated according to the Guide to the Expression of Uncertainty in Measurement (GUM) of ISO/BIPM and ISO Guide 35. The δ^82/76Se value of GBW 04457 relative to NIST SRM 3149 was also calculated. These iCRMs are intended for use in calibration of instruments and evaluation of methods for the determination of selenium isotope ratios.     

Twenty-five years of spectral UV-B measurements over Canada,Europe and Japan: Trends and effects from changes in ozone,aerosols, clouds,and surface reflectivity

Spectral UV records of solar irradiance at stations over Europe, Canada, and Japan were used to study long-term trends at 307.5 nm for a 25-year period, from 1992 to 2016. Ground-based measurements of total ozone, as well as satellite measurements of the Aerosol Index, the Total Cloud Cover and the surface reflectivity were also used in order to attribute the estimated changes of the UV to the corresponding changes of these factors. The present study shows that over the Northern Hemisphere, the long-term changes in UV-B radiation reaching the Earth's surface vary significantly over different locations, and that the main drivers of these variations are changes in aerosols and total ozone. At high latitudes, part of the observed changes may also be attributed to changes in the surface reflectivity. Over Japan, the UV-B irradiance at 307.5 nm has increased significantly by ∼3%/decade during the past 25 years, possibly due to the corresponding significant decrease of its absorption by aerosols. It was found that the greatest part of this increase took place before the mid-2000s. The only European station, over which UV radiation increases significantly, is that of Thessaloniki, Greece. Analysis of the clear-sky irradiance for the particular station shows increasing irradiance at 307.5 nm by ∼3.5%/decade during the entire period of study, with an increasing rate of change during the last decade, possibly again due to the decreasing absorption by aerosols.