Vegetation-groundwater dynamics at a former uranium mill site following invasion of a biocontrol agent: A time series analysis of Landsat normalized difference vegetation index data

Because groundwater recharge in dry regions is generally low, arid and semiarid environments have been considered well-suited for long-term isolation of hazardous materials (e.g., radioactive waste). In these dry regions, water lost (transpired) by plants and evaporated from the soil surface, collectively termed evapotranspiration (ET), is usually the primary discharge component in the water balance. Therefore, vegetation can potentially affect groundwater flow and contaminant transport at waste disposal sites. We studied vegetation health and ET dynamics at a Uranium Mill Tailings Radiation Control Act (UMTRCA) disposal site in Shiprock, New Mexico, where a floodplain alluvial aquifer was contaminated by mill effluent. Vegetation on the floodplain was predominantly deep-rooted, non-native tamarisk shrubs (Tamarix sp.). After the introduction of the tamarisk beetle (Diorhabda sp.) as a biocontrol agent, the health of the invasive tamarisk on the Shiprock floodplain declined. We used Landsat normalized difference vegetation index (NDVI) data to measure greenness and a remote sensing algorithm to estimate landscape-scale ET along the floodplain of the UMTRCA site in Shiprock prior to (2000–2009) and after (2010–2018) beetle establishment. Using groundwater level data collected from 2011 to 2014, we also assessed the role of ET in explaining seasonal variations in depth to water of the floodplain. Growing season scaled NDVI decreased 30% (p < .001), while ET decreased 26% from the pre- to post-beetle period and seasonal ET estimates were significantly correlated with groundwater levels from 2011 to 2014 (r² = .71; p = .009). Tamarisk greenness (a proxy for health) was significantly affected by Diorhabda but has partially recovered since 2012. Despite this, increased ET demand in the summer/fall period might reduce contaminant transport to the San Juan River during this period.     

The nature and origin of pigments in black opal from Lightning Ridge,New South Wales,Australia

Abstract

Black opal (opal-AG) owes its dark coloration to a fine-grained pigment commonly inferred to be mainly carbon, yet chemical compositions for black opals suggest there may be additional components. Here we search for such components in pigment concentrates prepared by dissolving black opal nodules (nobbies) from Lightning Ridge (NSW) in hydrofluoric acid, using electron microscopy (scanning electron microscopy, transmission electron microscopy), X-ray diffraction and laser-ablation ICP-MS. The results demonstrate the presence of sulfides—predominantly pyrite and chalcopyrite, with minor galena and Ti-oxide phases, as additional components of the pigment. ATR-FTIR analysis indicates the presence of C=O and C–H groups, consistent with an organic origin. Transmission electron microscopy images of pigment show variously deformed, originally spherical ～100?nm particles rich in sulfide and carbon, which are interpreted as thin coatings of pigment on now dissolved opaline silica spheres. Laser-ablation ICP-MS analysis identifies remnant silica in pigment concentrates, which may be interpreted as opaline silica surviving HF treatment protected as inclusions in sulfides. When examined within the context of petrographic observations from more than 1000 opal nodules (nobbies) at Lightning Ridge, these new results suggest that pigment carbon and sulfides in the nodules formed microbially under initially anoxic groundwater conditions, within pre-existing cavities concurrently being filled with silica sol ultimately derived from chemical weathering of feldspar-rich volcaniclastic sediment. Intensely black pigment layers observed at the floor of many nodules indicate settling of dark, high-density (sulfide–Ti-oxide-rich) pigment within cavities, with the implication that sulfate-reducing bacterial (SRB) activity commences early during the silica sol-gel ripening process. Microbial activity may persist until after the cavity has completely filled with the silica sol, as illustrated by abundant black opals with uniformly distributed pigment. Pigment formed at this stage may no longer be able to settle out within the ripening and increasingly viscous silica gel, thus forming pigmentation throughout the opal cavity. The existence of ‘amber’, pigment-poor opal with intensely black basal pigment layers is interpreted as signalling a lack of sulfate to sustain further SRB activity, or a change to more oxidising conditions, possibly related to interaction with surface waters within a downward-penetrating weathering front. A change in redox conditions would shut off activity of SRB and thus sulfide pigment production and allow development of aerobic microbial activity as described by others.     

Ein Drei-Bereiche-Modell zur Beschreibung der Adsolubilisation von aromatischen Verbindungen an tensidbelegtem Silikagel

A Three-site Model To Describe the Adsolubilization of Aromatic Compounds at Surfactant Covered Silica Surfactant modified minerals are undergoing increasing interest because of their potential use as environmental sorbents. In this work, partition coefficients of several aromatic compounds between the aqueous phase and hexadecyltrimethylammoniumbromide covered silica were determined in batch experiments and by a chromatographic method. The surfactant coverage density and the pH were varied in the experiments. It turns out that changing the pH led to different partition coefficients at the same coverage density. The dependency of the partitioning on pH is explained by the formation of different types of surfactant aggregates at the silica surface. The effect of the formation of different types of aggregates is modified by the polarity and the functional groups of the aromatic compounds. We suggest a three-site model to explain this observation. In the three-site model, three regions in the surfactant aggregates are distinguished in which the aromatic compounds can be sorbed into the surfactant aggregates.     

Ground-Motion Scaling in the Western Alps

In order to empirically obtain the scaling relationships for the high-frequency ground motion in the Western Alps (NW Italy), regressions are carried out on more than 7500 seismograms from 957 regional earthquakes. The waveforms were selected from the database of 6 three-component stations of the RSNI (Regional Seismic network of Northwestern Italy). The events, M _W ranging between 1.2 and 4.8, were recorded within a hypocentral distance of 200 km during the time period: 1996–2001. The peak ground velocities are measured in selected narrow-frequency bands, between 0.5 and 14 Hz. Results are presented in terms of a regional attenuation function for the vertical ground motion, a set of vertical excitation terms at the reference station STV2 (hard-rock), and a set of site terms (vertical and horizontal), all relative to the vertical component of station STV2.The regional propagation of the ground motion is modeled after quantifying the expected duration of the seismic motion as a function of frequency and hypocentral distance. A simple functional form is used to take into account both the geometrical and the anelastic attenuation: a multi-variable grid search yielded a quality factor Q(f) = 310f ^0.20, together with a quadri-linear geometrical spreading at low frequency. A simpler, bi-linear geometrical spreading seems to be more appropriate at higher frequencies (f > 1.0 Hz). Excitation terms are matched by using a Brune spectral model with variable, magnitude-dependent stress drop: at M _w 4.8, we used Δσ = 50 MPa. A regional distance-independent attenuation parameter is obtained (κ₀ = 0.012 s) by modelling the average spectral decay at high frequency of small earthquakes.In order to predict the absolute levels of ground shaking in the region, the excitation/attenuation model is used through the Random Vibration Theory (RVT) with a stochastic point-source model. The expected peak-ground accelerations (PGA) are compared with the ones derived by Ambraseys et al. (1996) for the Mediterranean region and by Sabetta and Pugliese (1996) for the Italian territory.     

Target-oriented imaging using extended image volumes: a low-rank factorization approach

Imaging in geological challenging environments has led to new developments, including the idea of generating reflection responses by means of interferometric redatuming at a given target datum in the subsurface, when the target datum lies beneath a complex overburden. One way to perform this redatuming is via conventional model-based wave-equation techniques. But those techniques can be computationally expensive for large-scale seismic problems since the number of wave-equation solves is equal to two times the number of sources involved during seismic data acquisition. Also conventional shot-profile techniques require lots of memory to save full subsurface extended image volumes. Therefore, we can only form subsurface image volumes in either horizontal or vertical directions. To exploit the information hidden in full subsurface extended image volumes, we now present a randomized singular value decomposition-based approach built upon the matrix probing scheme, which takes advantage of the algebraic structure of the extended imaging system. This low-rank representation enables us to overcome both the computational cost associated with the number of wave-equation solutions and memory usage due to explicit storage of full subsurface extended image volumes employed by conventional migration methods. Experimental results on complex geological models demonstrate the efficacy of the proposed methodology and allow practical reflection-based extended imaging for large-scale five-dimensional seismic data.     

A GIS based approach to back trajectory analysis for the source apportionment of aerosol constituents and its first application

A new method has been developed to combine back trajectory statistics with a detailed land cover analysis. It provides numeric proxies for the residence times of sampled air masses above certain land cover classes (marine, natural vegetation, agricultural lands, urban areas, and bare areas), as well as further meteorological parameters (mean trajectory length, solar radiation along trajectory, and local height of the boundary mixing layer). The method has been implemented into a GIS-enabled database system to allow for an efficient processing of large datasets with low computational demands. A principal component analysis was performed on a dataset including the modelled residence times, the modelled meteorological parameters, some measured meteorological parameters (wind speed and temperature), and the concentrations of 10 particle constituents (inorganic ions and organic and elemental carbon) in 5 particle size ranges for 29 winter- and summertime samples at an urban background site in Leipzig, Germany. Six principal components could be extracted which together explained about 80% of the total variance in the dataset. The factors could be attributed to the influence of meteorology to continental background pollution, secondary formation processes in polluted air masses, wood burning, aged sea-salt, local traffic, and long-range transported crustal material. The modelled residence times and the meteorological parameters were generally consistent with the existing knowledge of specific particle sources and thereby facilitated and strengthened the interpretation of the factors. Moreover, they allowed for a clear distinction between continental background pollution and secondary formation processes, which has not been possible in previous source apportionment studies. The results demonstrate that the combined usage of back trajectory, land cover, and meteorological data by the presented method yields valuable additional information on the history of sampled air masses, which can improve the quality of source apportionment of atmospheric aerosol constituents.     

Comparison of Five Eulerian Air Pollution Forecasting Systems for the Summer of 1999 Using the German Ozone Monitoring Data

Eulerian state-of-the-art air pollution forecasting systems on the European scale are operated routinely by several countries in Europe. DWD and FUB, both Germany, NERI, Denmark, NILU, Norway, and SMHI, Sweden, operate some of these systems. To apply such modeling systems, e.g. for regulatory purposes according to new EU directives, an evaluation and comparison of the model systems is fundamental in order to assess their reliability. One step in this direction is presented in this study: The model forecasts from all five systems have been compared to measurements of ground level ozone in Germany. The outstanding point in this investigation is the availability of a huge amount of data – from forecasts by the different model systems and from observations. This allows for a thorough interpretation of the findings and assures the significance of the observed features. Data from more than 300 measurement stations for a 5-month period (May–September 1999) of the German monitoring networks have been used in this comparison. Different spatial and temporal statistical parameters were applied in the evaluation. Generally, it was found that the most comprehensive models gave the best results. However, the less comprehensive and computational cheaper models also produced good results. The extensive comparison made it possible to point out weak points in the different models and to describe the individual model behavior for a full summer period in a climatological sense. The comparison also gave valuable information for an assessment of individual measurement stations and complete monitoring networks in terms of the representativeness of the observation data.     

Mechanism development and modelling of tropospheric multiphase halogen chemistry: The CAPRAM Halogen Module 2.0 (HM2)

A new detailed multiphase halogen mechanism, the CAPRAM Halogen Module 2.0 (HM2), has been developed and coupled to the multiphase chemistry mechanism RACM-MIM2ext/CAPRAM 3.0n. The overall mechanism comprises 1,705 reactions including 595 reactions of the HM2. Halogen chemistry box model studies have been, for the first time, performed with a non-permanent cloud scenario for pristine open ocean regions in mid-latitudes. Moreover, detailed time-resolved reaction flux analysis has been used to investigate the multiphase halogen reaction cycles in more detail. Clouds significantly change the multiphase halogen chemical system and new reaction cycles are proposed for in-cloud conditions. While most gas phase concentrations are decreased for chlorine and iodine species, they are increased for bromine. Flux analyses determined the relative contributions of the methylene dihalides CH₂IX (X = Cl, Br, I) as the main I atom source with a contribution of about 80 % to the total iodocarbon sources. Furthermore, HOI was confirmed to be important for chlorine activation. It is shown that 25 % of the ozone loss can be attributed to halogens. VOC oxidation by halogens is important as halogens account for about 20 % of the methane oxidation and up to 80 % of the oxidation of other VOCs. In other cases, enhanced VOC and VOC oxidation product concentration levels were found. For example, 15 % of the methyl peroxyl radicals are formed after the reaction of chlorine atoms with methane or methyl hydroperoxide. In the aqueous phase, changes in the oxidation of organics do only occur for highly oxidised organics without a C-H bond. For example, over 80 % of oxalic acid are oxidised by electron transfer with Cl₂ ^? in deliquescent particles during non-cloud periods.     

Scenarios for Modeling Multiphase Tropospheric Chemistry

Besides observational data model calculations are a very importanttool for improving our understanding of multiphase chemistryin the troposphere. Before a chemical model can be used for that purposeit is necessary to show that the model does what it is intendedto do. A protocol has been developed thatcan be used as a basis for the verification of the numericsand the correct implementation of thechemical balance equations.The protocol defines meteorological parameters and initial conditionsfor a zerodimensional (box) model. Several scenarios cover the pollutedas well as the remote marine and continental boundary layer and also thefree troposphere. Calculations by different groupswith different modelsand numerical solvers demonstrate that the protocol is clear and complete.The excellent agreement between the results of all groups are a major step of verification of the participating models.The scenarios may also serve as well documented base cases forsensitivity studies.