Using Sequential Self-Calibration Method to Identify Conductivity Distribution: Conditioning on Tracer Test data

An iterative inverse method, the sequential self-calibration method, is developed for mapping spatial distribution of a hydraulic conductivity field by conditioning on nonreactive tracer breakthrough curves. A streamline-based, semi-analytical simulator is adopted to simulate solute transport in a heterogeneous aquifer. The simulation is used as the forward modeling step. In this study, the hydraulic conductivity is assumed to be a deterministic or random variable. Within the framework of the streamline-based simulator, the efficient semi-analytical method is used to calculate sensitivity coefficients of the solute concentration with respect to the hydraulic conductivity variation. The calculated sensitivities account for spatial correlations between the solute concentration and parameters. The performance of the inverse method is assessed by two synthetic tracer tests conducted in an aquifer with a distinct spatial pattern of heterogeneity. The study results indicate that the developed iterative inverse method is able to identify and reproduce the large-scale heterogeneity pattern of the aquifer given appropriate observation wells in these synthetic cases.     

运用STOMP数值模拟器进行多空介质中尿素水解以及方解石沉淀数值模拟

遍布于美国能源部所辖野外场地的地下放射性核素和痕量金属污染对美国能源部的长期监管职能提出了最严峻的挑战。目前较有前景的一个原地固化污染物方法是利用尿素酶催化尿素水解,进而推动方解石沉淀及与之共同沉淀。运用STOMP数值模拟器实施一系列反应性溶质数值模拟来研究系统中水流、运移、反应与反应导致的介质性质变化之间的紧密的非线性耦合,并在研究中系统性增加反应系统的复杂度和物理-化学的非均质性以考察模拟效果。通过调整反应速率表达式中尿素酶浓度和沉淀反应常数,对尿素水解和方解石沉淀速率进行了敏感度分析。结果显示:尿素水解速率是矿物沉淀总量的决定性因素,沉淀物的空间分布既取决于水解速率也取决于沉淀速率;在一维水流柱6个空隙体积的时间段模拟中,观察到了5%的最大空隙度减少;当二维模拟中考虑一个低渗透性包裹体时,被改变的流场重新对生成的沉淀组成物进行了分布,形成了扭曲的反应锋面。模拟结果也显示沉淀在低渗透性区域边界上形成,这说明位于这个区域里面的污染物可能被囊裹而从流线中被隔离起来。     

Polarized Raman spectra of beryl and bazzite

The polarized Raman spectra of four different beryl crystals were studied at room temperature in the range from 30 to 4000 cm^-1. The spectra show significant differences between the samples studied, and corrections are proposed for the reference Raman spectra of beryl previously reported by Adams and Gardner (1974). Type II water is observed in two crystals; the corresponding symmetric Raman stretching band at 3595 cm^-1 is extremely strong for an impurity (about 20% of the strongest beryl lattice mode). Another, sharper, band of similar intensity at 3605 cm^-1 could possibly originate from a hydroxyl stretching mode. Additional weaker bands are observed around 1600 cm^-1 and 3600–3750 cm^-1. The first polarized Raman spectra of bazzite are presented and discussed.     

Observations of the convective plume of a lake under cold-air advective conditions

Moderating effects of Lake Apopka, Florida on downwind surface temperatures were evaluated under cold-air advective conditions. Point temperature measurements north and south of the lake and data obtained from a thermal scanner flown at 1.6 km indicate that surface temperatures directly downwind may be higher than surrounding surface temperatures by as much as 5 °C under conditions of moderate winds (~4 m s^–1). No substantial temperature effects were observed with surface wind speed less than 1 m s^–1. Fluxes of sensible and latent heat from Lake Apopka were calculated from measurements of lake temperature, net radiation, relative humidity and air temperature above the lake. Bulk transfer coefficients and the Bowen ratio were calculated and found to be in agreement with reported data for non-advective conditions.IFAS Journal Series No. 1006.     

Use of high spectral resolution remote sensing to determine leaf palatability of eucalypt trees for folivorous marsupials

In this study we assess the feasibility of remotely measuring canopy biochemistry, and thus the potential for conducting large-scale mapping of habitat quality. A number of studies have found nutrient composition of eucalypt foliage to be a major determinant of the distribution of folivorous marsupials. More recently it has been demonstrated that a specific group of secondary plant chemicals, the diformylphloroglucinols (DFPs), are the most important feeding deterrents, and are thus vital determinants of habitat quality. We report on the use of laboratory spectroscopy to attempt to identify one such DFP, sideroxylonal-A, in the foliage of Eucalyptus melliodora, one of the few eucalypt species browsed by folivorous marsupials. Reflectance spectra were obtained for freeze-dried, ground leaves using near infrared spectroscopy (NIRS) and for both oven-dried and fresh whole leaves using a laboratory-based (FieldSpec) spectroradiometer. Modified partial least squares (MPLS) regression was used to develop calibration equations for sideroxylonal-A concentration based on the reflectance spectra transformed as both the first and second difference of absorbance (Log 1/R). The predictive ability of the calibration equations was assessed using the standard error of calibration statistic (SECV). Coefficients of determination (r²) were highest for the ground leaf spectra (0.98), followed by the fresh leaf and dry leaf spectra (0.94 and 0.87, respectively). When applied to independent validation sub-sets, sideroxylonal-A was most accurately predicted from the ground leaf spectra (r² = 0.94), followed by the dry leaf and fresh leaf spectra (0.72 and 0.53, respectively). Two spectral regions, centred on 674 nm and 1394 nm, were found to be highly correlated with sideroxylonal-A concentration for each of the three spectral data sets studied. Results from this study suggest that calibration equations derived from modified partial least squares regression may be used to predict sideroxylonal-A concentration, and hence leaf palatability, of Eucalyptus melliodora trees, thereby indicating that the remote estimation of habitat quality of eucalypt forests for marsupial folivores is feasible.     

Miocene silicic volcanism in southwestern Idaho: geochronology,geochemistry, and evolution of the central Snake River Plain

New ⁴⁰Ar-³⁹Ar geochronology, bulk rock geochemical data, and physical characteristics for representative stratigraphic sections of rhyolite ignimbrites and lavas from the west-central Snake River Plain (SRP) are combined to develop a coherent stratigraphic framework for Miocene silicic magmatism in this part of the Yellowstone ‘hotspot track’. The magmatic record differs from that in areas to the west and east with regard to its unusually large extrusive volume, broad lateral scale, and extended duration. We infer that the magmatic systems developed in response to large-scale and repeated injections of basaltic magma into the crust, resulting in significant reconstitution of large volumes of the crust, wide distribution of crustal melt zones, and complex feeder systems for individual eruptive events. Some eruptive episodes or ‘events’ appear to be contemporaneous with major normal faulting, and perhaps catastrophic crustal foundering, that may have triggered concurrent evacuations of separate silicic magma reservoirs. This behavior and cumulative time-composition relations are difficult to relate to simple caldera-style single-source feeder systems and imply complex temporal-spatial development of the silicic magma systems. Inferred volumes and timing of mafic magma inputs, as the driving energy source, require a significant component of lithospheric extension on NNW-trending Basin and Range style faults (i.e., roughly parallel to the SW–NE orientation of the eastern SRP). This is needed to accommodate basaltic inputs at crustal levels, and is likely to play a role in generation of those magmas. Anomalously high magma production in the SRP compared to that in adjacent areas (e.g., northern Basin and Range Province) may require additional sub-lithospheric processes. Electronic supplementary material The online version of this article (doi: ) contains supplementary material, which is available to authorized users. This paper constitutes part of a special issue dedicated to Bill Bonnichsen on the petrogenesis and volcanology of anorogenic rhyolites.     

Experimental and modeling study on Cr(VI) transfer from soil into surface runoff

With the development of modern agriculture, large amount of fertilizer and pesticide outflow from farming land causes great waste and serious pollution to surface water and groundwater, and threatens ecological environment and even human life. In this paper, laboratory experiments are conducted to simulate adsorbed Cr(VI) transfer from soil into runoff. A two-layer in-mixing analytical model is applied to analyze laboratory experimental results. A data assimilation (DA) method via the ensemble Kalman filter (EnKF) is used to update parameters and improve modeling results. In comparison with the experimental data, DA updated modeling results are much better than those without the updating. To make predictions better, the inflation method with a constant inflation factor via DA method was used to compensate the fast decrease of ensemble spread partially related to filter inbreeding. Based on the used rainfall and relevant physical principles, the updated value of the incomplete mixing coefficient γ is about 14.0 times of the value of the incomplete mixing coefficient α in experiment 1 and about 7.4 times in experiment 2, while the difference between the flow rate of runoff and infiltration is not so large even after reaching stable infiltration condition. The results indicate the loss of Cr(VI) in soil solute is mainly due to infiltration, rather than surface runoff. With the increase of mixing layer depth, soil adsorption capacity will increase and the loss of soil solute will decrease. The study results provide information for reducing and even preventing the agricultural nonpoint source pollution.     

Potential of legacy 2D seismic data for deep targeting and structural imaging at the Neves–Corvo massive sulphide-bearing deposit,Portugal

Seismic methods are becoming an established choice for deep mineral exploration after being extensively tested and employed for the past two decades. To investigate whether the early European mineral-exploration datasets had potential for seismic imaging that was overlooked, we recovered a low-fold legacy seismic dataset from the Neves–Corvo mine site in the Iberian Pyrite Belt in southern Portugal. This dataset comprises six 4–6 km long profiles acquired in 1996 for deep targeting. Using today's industry-scale processing algorithms, the world-class, ca. 150 Mt, Lombador massive sulphide and other smaller deposits were better imaged. Additionally, we also reveal a number of shallow but steeply dipping reflections that were absent in the original processing results. This study highlights that legacy seismic data are valuable and should be revisited regularly to take advantage of new processing algorithms and the experiences gained from processing such data in hard-rock environments elsewhere. Remembering that an initial processing job in hard rock should always aim to first obtain an overall image of the subsurface and make reflections visible, and then subsequent goals of the workflow could be set to, for example understanding relative amplitude ratios. The imaging of the known mineralization implies that this survey could likely have been among one of the pioneer studies in the world that demonstrated the capability of directly imaging massive sulphide deposits using the seismic method.