Review: Simulation-optimization models for the management and monitoring of coastal aquifers

The literature on the application of simulation-optimization approaches for management and monitoring of coastal aquifers is reviewed. Both sharp- and dispersive-interface modeling approaches have been applied in conjunction with optimization algorithms in the past to develop management solutions for saltwater intrusion. Simulation-optimization models based on sharp-interface approximation are often based on the Ghyben-Herzberg relationship and provide an efficient framework for preliminary designs of saltwater-intrusion management schemes. Models based on dispersive-interface numerical models have wider applicability but are challenged by the computational burden involved when applied in the simulation-optimization framework. The use of surrogate models to substitute the physically based model during optimization has been found to be successful in many cases. Scalability is still a challenge for the surrogate modeling approach as the computational advantage accrued is traded-off with the training time required for the surrogate models as the problem size increases. Few studies have attempted to solve stochastic coastal-aquifer management problems considering model prediction uncertainty. Approaches that have been reported in the wider groundwater management literature need to be extended and adapted to address the challenges posed by the stochastic coastal-aquifer management problem. Similarly, while abundant literature is available on simulation-optimization methods for the optimal design of groundwater monitoring networks, applications targeting coastal aquifer systems are rare. Methods to optimize compliance monitoring strategies for coastal aquifers need to be developed considering the importance of monitoring feedback information in improving the management strategies.     

Clusters of Galaxies and the Cosmic Web with Square Kilometre Array

The intra-cluster and inter-galactic media that pervade the large scale structure of the Universe are known to be magnetized at sub-micro Gauss to micro Gauss levels and to contain cosmic rays. The acceleration of cosmic rays and their evolution along with that of magnetic fields in these media is still not well understood. Diffuse radio sources of synchrotron origin associated with the Intra-Cluster Medium (ICM) such as radio halos, relics and mini-halos are direct probes of the underlying mechanisms of cosmic ray acceleration. Observations with radio telescopes such as the Giant Metrewave Radio Telescope, the Very Large Array and the Westerbork Synthesis Radio Telescope have led to the discoveries of about 80 such sources and allowed detailed studies in the frequency range 0.15–1.4 GHz of a few. These studies have revealed scaling relations between the thermal and non-thermal properties of clusters and favour the role of shocks in the formation of radio relics and of turbulent re-acceleration in the formation of radio halos and mini-halos. The radio halos are known to occur in merging clusters and mini-halos are detected in about half of the cool-core clusters. Due to the limitations of current radio telescopes, low mass galaxy clusters and galaxy groups remain unexplored as they are expected to contain much weaker radio sources. Distinguishing between the primary and the secondary models of cosmic ray acceleration mechanisms requires spectral measurements over a wide range of radio frequencies and with high sensitivity. Simulations have also predicted weak diffuse radio sources associated with filaments connecting galaxy clusters. The Square Kilometre Array (SKA) is a next generation radio telescope that will operate in the frequency range of 0.05–20 GHz with unprecedented sensitivities and resolutions. The expected detection limits of SKA will reveal a few hundred to thousand new radio halos, relics and mini-halos providing the first large and comprehensive samples for their study. The wide frequency coverage along with sensitivity to extended structures will be able to constrain the cosmic ray acceleration mechanisms. The higher frequency (>5 GHz) observations will be able to use the Sunyaev–Zel’dovich effect to probe the ICM pressure in addition to tracers such as lobes of head–tail radio sources. The SKA also opens prospects to detect the ‘off-state’ or the lowest level of radio emission from the ICM predicted by the hadronic models and the turbulent re-acceleration models.     

Classification of rocks radionuclide data using machine learning techniques

The aim of this study is to assess the performance of linear discriminate analysis, support vector machines (SVMs) with linear and radial basis, classification and regression trees and random forest (RF) in the classification of radionuclide data obtained from three different types of rocks. Radionuclide data were obtained for metamorphic, sedimentary and igneous rocks using gamma spectroscopic method. A P-type high-purity germanium detector was used for the radiometric study. For analysis purpose, we have determined activity concentrations of ²³²Th, ²²⁶Ra and ⁴⁰K radionuclides, published elsewhere (Rafique et al. in Russ Geol Geophys 55:1073–1082, 2014), in different rock samples and built the classification model after pre-processing the data using three times tenfold cross-validation. Using this model, we have classified the new samples into known categories of sedimentary, igneous and metamorphic. The statistics depicts that RF and SVM with radial kernel outperform as compared to other classification methods in terms of error rate, area under the curve and with respect to other performance measures.     

Seismic activity prediction using computational intelligence techniques in northern Pakistan

Earthquake prediction study is carried out for the region of northern Pakistan. The prediction methodology includes interdisciplinary interaction of seismology and computational intelligence. Eight seismic parameters are computed based upon the past earthquakes. Predictive ability of these eight seismic parameters is evaluated in terms of information gain, which leads to the selection of six parameters to be used in prediction. Multiple computationally intelligent models have been developed for earthquake prediction using selected seismic parameters. These models include feed-forward neural network, recurrent neural network, random forest, multi layer perceptron, radial basis neural network, and support vector machine. The performance of every prediction model is evaluated and McNemar’s statistical test is applied to observe the statistical significance of computational methodologies. Feed-forward neural network shows statistically significant predictions along with accuracy of 75% and positive predictive value of 78% in context of northern Pakistan.     

柴达木盆地昆特依盐湖沉积特征及其盐类资源评价

昆特依盐湖位于柴达木盆地西北部,是一个相对独立的沉积盆地,可反映盆地西北部的演化历史,对研究盆地内盐湖演化有重要意义。基于研究区已有的研究成果,总结和评价了昆特依盐湖沉积与盐类矿产资源,并得到以下结论:1)昆特依湖自1 925～774 ka期间以暖湿型气候为主,其后呈干冷与暖湿气候交替出现,最早成盐期为1150 ka,最终于30.6 ka成为干盐湖。昆特依盆地早期沉积的细碎屑沉积物沉积速率(0.22 mm/a)慢于晚期沉积的石盐层(0.33 mm/a);2)昆特依盐湖液体KCl储量为1.32×10~8 t,MgCl储量为2.79×10~8 t,NaCl储量为24.15×10~8 t;3)估算出赋存于昆特依盐湖石盐矿中的锂资源量达39×10~4 t,同时估算出液体中的锂资源量约有74×10~4 t,二者总量达百万吨;4)昆特依矿田的成因属典型的内陆湖相沉积,成矿控制因素以构造作用为主,其次为气候及物质来源。关于昆特依盐湖沉积特征与盐类矿产资源特征的研究,有利于对本地区地质演化更全面的认识,对合理开发利用盐类资源有重要意义。     

Magnetopause boundary structure deduced from the high-time resolution particle experiment on the Equator-S spacecraft

An electrostatic analyser (ESA) onboard the Equator-S spacecraft operating in coordination with a potential control device (PCD) has obtained the first accurate electron energy spectrum with energies &7 eV-100 eV in the vicinity of the magnetopause. On 8 January, 1998, a solar wind pressure increase pushed the magnetopause inward, leaving the Equator-S spacecraft in the magnetosheath. On the return into the magnetosphere approximately 80 min later, the magnetopause was observed by the ESA and the solid state telescopes (the SSTs detected electrons and ions with energies &20–300 keV). The high time resolution (3 s) data from ESA and SST show the boundary region contains of multiple plasma sources that appear to evolve in space and time. We show that electrons with energies &7 eV–100 eV permeate the outer regions of the magnetosphere, from the magnetopause to &6Re. Pitch-angle distributions of &20–300 keV electrons show the electrons travel in both directions along the magnetic field with a peak at 90° indicating a trapped configuration. The IMF during this interval was dominated by Bx and By components with a small Bz.     

Origin of Boron in the Gas Hure Salt Lake of NorthwesternQaidam Basin,China: Evidence from Hydrochemistry andBoron Isotopes

The Gas Hure Salt Lake(GHSL) in the northwestern Qaidam Basin, western China, is rich in boron(B)resources, but its B-resource origin is hardly known. Hydrochemical compositions and B isotope characteristics of different waters were collected around the GHSL, including the river water, stream water, spring water, salt-lake brine,intercrystalline brine, well water, drilling brine, and solar pond brine. The hydrochemical signatures suggest that silicates,carbonates and evaporates are the main B-bearing rocks during the water dynamic. The reservoir estimation of B resources shows that the Kulamulekesay River(KLMR) and the Atekan River(ATKR) contribute annually 18.3 tons and 22.84 tons of B, respectively, with a total amount of 11.72 × 10~4 tons of B during the past 5.7 ka. In comparison with the known B reservoir(32.96 × 10~4 tons) in the GHSL, a significant amount of B in the GHSL was probably recharged from deep fluids and sediments around the GHSL. The B concentration and B-enrichment degree are shaped by the evaporation process,which are highly elevated at the carnallite and bischofite stages.     

Aggregate dust model to study the polarization properties of comet C/1996 B2 Hyakutake

The observed linear polarization data of comet Hyakutake are studied at wavelengths λ=0.365 μm,λ=0.485 μm and 0.684 μm through simulations using Ballistic Particle-Cluster Aggregate and Ballistic Cluster-Cluster Aggregate aggregates of 128 spherical monomers.We first found that the size parameter of the monomer,x～1.56-1.70,turned out to be the most suitable which provides the best fits to the observed dust scattering properties at three wavelengths:λ=0.365 μm,0.485 μm and 0.684 μm.Thus,the effective radius of the aggregate(r)lies in the range 0.45 μm ≤ r≤ 0.49 μm at λ=0.365 μm;0.60μm < r < 0.66 μm at λ=0.485 μmand 0.88 μm < r < 0.94 μm at λ=0.684 μm.Now using superposition T-MATRIX code and the power-law size distribution,n(r)～r-3,the best-fitting values of complex refractive indices are calculated for the observed polarization data at the above three wavelengths.The best-fitting complex refractive indices(n,k)are found to be (1.745,0.095)at λ=0.365 μm,(1.743,0.100)at λ=0.485 μm and(1.695,0.100)at λ=0.684 μm.The refractive indices derived from the present analysis correspond to a mixture of both silicates and organics,which are in good agreement with the in situ measurement of comets by different spacecraft.     

Ammonia emission from subtropical crop land area in India

Ammonia (NH₃) emission from wheat (November to April) and rice (July to October) crops was measured using the chemiluminescence method at a subtropical agricultural area of India during 2009?C2010. Samples were collected from the canopy height during different growth stages of wheat crop to study the variations of NH₃ emission during different growth stages of the crop. Background atmospheric concentration of NH₃ was measured at 5 m height at the study site. Background NH₃ concentration was subtracted from the NH₃ concentration at crop canopy height to estimate the emission of NH₃ from crop canopy. The NH₃ emission from the wheat crop were recorded as 33.3 to 57.0; 15.3 to 29.2; 10.3 to 28.0; 8.7 to 23.9 and 13.9 to 28.9 ??g m^?2 d^?1 during sowing, crown root initiation (CRI), panicle initiation, grain filling and maturity stages of the crop respectively. The NH₃ emission followed a diurnal pattern with significant correlation with ambient temperature at different crop growth stages. Cumulative seasonal NH₃ emission to the atmosphere was accounted for the loss of ??10% of applied N-fertilizer during the wheat crop growing period. Immediate increase in NH₃ emission was recorded from rice crop, grown under temperature gradient tunnel (TGT). However, the NH₃ emission inside the TGT decreases within 3?C4 h after the N-fertilizer application. Continuous estimation of NH₃ concentration at the crop canopy inside the TGT, suggests that the NH₃ emission to the atmosphere reaches its peak within ??20 h of N-fertilizer application and continues up to 5 d following a diurnal pattern.     

Effect of soil aging on arsenic fractionation and bioaccessibility in inorganic arsenical pesticide contaminated soils

Ingestion of As – contaminated soil by children is a growing concern in former agricultural lands converted to residential or recreational land use areas. The mobility and bioavailability of As is controlled by its reactions with soil particles. The degree and strength of As retention by soil constituents may vary greatly with time. The present authors hypothesize that aging results in reduced mobility of As thereby decreasing As release and its bioavailability. The present study is aimed at evaluating the effect of aging on soil As fractionation and bioaccessibility in a temperature and humidity-controlled greenhouse setting. The design allowed the evaluation of dynamic interactions between soils, pesticides, water, and plants. Therefore, 4 soil types (Immokalee, Millhopper, Pahokee Muck, and Orelia) were selected based on their potential differences in As reactivity. The soils were amended with the pesticide Na arsenate at two rates. Rice was used as the test crop. Soil samples collected after different time periods (0, 6 months, 1 a and 3 a) were extracted for soil-As forms via a sequential extraction technique. Bioaccessible As was extracted via an in vitro gastrointestinal method. At time 0, most of the extractable As in soil was in the soluble form, resulting in high bioaccessibility. As expected, soluble and exchangeable fractions decreased with time for up to 6 months, but remained constant thereafter. After 3 a of soil–pesticide equilibration, As bioaccessibility was still high in all the soils except for the Pahokee Muck. No significant difference in As bioaccessibility was observed between the soils. Arsenic was present predominantly as As(V) with 5–10% of the total dissolved As being present as As(III). Data obtained suggest that although aging had an impact on the geochemical forms, gastric pH was the sole important factor effecting As bioaccessibility.