Using Models of the Deep Structure of the Earth’s Crust and Upper Mantle Based on GOCE Gravity Data in Metallogenic Analysis

Izvestiya, Atmospheric and Oceanic Physics - Gravimetric data obtained by the GOCE spacecraft contributes to the development of global models of the deep structure of the Earth’s crust and...     

Natural environmental changes and human impact reflected in sediments of a high alpine lake in Switzerland

From the high alpine Sägistalsee (1935 m a.s.l.), 13.50 m of continuously laminated sediments comprising the last 9050 years, were analyzed. Even though Sägistalsee is a high elevation site, human-induced environmental changes start as early as 4300 cal. BP and leave a clearly detectable signal in the mineralogy of the sediments, which is much stronger than the signal from natural environmental changes that occurred before this time. Variations in the physical and mineralogical sediment properties of this clastic sequence reflect erosional changes in the catchment, where almost pure limestone contrasts with carbonaceous, quartz-bearing marl, and shist. The calcite/quartz (Cc/Qz) ratio was found to be most indicative of these changes, which occurred around AD 1850 and at 650, 2000, 3700, and 6400 cal. BP. The first four are interpreted as erosion events, which are related to human-induced changes in the vegetation cover and land use. We associate them to the recent development of tourism and grazing, the medieval intensification of pasturing, Roman forest clearance, and Bronze Age forest clearance, respectively. The Cc/Qz-ratio increases significantly within less than 100 years during these events, reflecting the erosion of unweathered or poorly weathered soils. The time intervals in between are characterized by a gradually decreasing Cc/Qz-ratio and reflect the stabilization or formation of new soils. Only the change at 6400 cal. BP, which represents the initial gradual stabilization of the catchment, is related to the immigration of Picea abies.     

Interrelationship between geochemical elements of sediment and groundwater at Samrak Park Delta of Nakdong River Basin in Korea: multivariate statistical analyses and artificial neural network approaches

This research was conducted at Samrak Park Delta of Nakdong River Basin in Busan Metropolitan City, Korea. The main objective of this study was to evaluate the interrelationship of geochemical elements in sediments and groundwater through multivariate statistical analyses and a multilayer perceptron artificial neural network model. The mean concentrations of chemical elements were Si (46%), Fe (16.9%), Al (15.7%), K (7.5%) and Ca (4.5%) in sediments, and Na (8650 mg/L), Mg (999 mg/L), Ca (432 mg/L), K (293 mg/L) and Cl (17,640 mg/L) in groundwater, respectively. The principal component analysis produced 3 kinds of factors with the variances of 63.37, 27.02 and 9.62%, respectively. It is suggested that the chemical components of sediments and groundwater were mainly originated from source rocks and seawater intrusion, with the minor impacts of irrigation and industry. Cluster analysis also showed that chemical elements were mainly controlled by the natural geogenic sources and seawater intrusion. Multilayer perceptron of artificial neural network (ANN) presented the good interrelationship between sediment and groundwater. The determination coefficients (R ²) between ANN predicted values and observed values in groundwater showed the high values of 0.61–0.97 except Mg, Mn and Sr. It is revealed that the chemical components of sediment and groundwater were derived from local geological origin and from the minor impact of anthropogenic sources. Multivariate analyses and ANN contributed to the identification of the mutual relationship between the geochemical elements of sediment and those of groundwater.     

Karst groundwater protection in the Kupa River catchment area and sustainable development

One of the most significant water resources in the Republic of Croatia is the catchment area of the Kupa River, located in the region bordering the Republic of Slovenia. About 88% of the total amount of water in this catchment originates in Croatia and just 12% from Slovenia; therefore, the largest part of the catchment area (about 1000 km²) is on the Croatian side of the border. It is a typical karst area of the Dinarides with aquifers characterized by a relatively rapid water exchange, high groundwater flow velocities and aquifers open to human impact from the surface. Consequently, the aquifers are highly vulnerable and at risk. Due to the availability of large quantities of high-quality spring water (about 6 m³/s), the entire area has a strategic importance within the context of any future development strategy pertaining to the western part of Croatia. The catchment area on the Croatian side was investigated using a wide range of research methods that included a classical hydrogeological approach, the detailed hydrologic calculation of water balance to the hydrogeochemical analyses and modelling. The objective was to determine protection zones and protection measures for the whole area. The difficulties are increased due to the fact that the karst catchment area is crossed by major traffic corridors, oil pipelines and a railway and that many settlements and a highly developed wood industry are present. The combination of protecting water resources with adequate prevention measures and necessary remedial activities that should satisfy the very strict requirements necessary for the protection of the karst aquifers while still allowing for present and future human activities is difficult – but not impossible – to achieve. One good example is the present highway with a closed dewatering system and waste water treatment before the water passes into the karst underground system.     

Analysis of snowpack dynamics during the spring melt season for a sub‐alpine site using point measurements and numerical modeling

Snowpack dynamics through October 2014–June 2017 were described for a forested, sub‐alpine field site in southeastern Wyoming. Point measurements of wetness and density were combined with numerical modeling and continuous time series of snow depth, snow temperature, and snowpack outflow to identify 5 major classes of distinct snowpack conditions. Class (i) is characterized by no snowpack outflow and variable average snowpack temperature and density. Class (ii) is characterized by short durations of liquid water in the upper snowpack, snowpack outflow values of 0.0008–0.005 cm hr^?1, an increase in snowpack temperature, and average snow density between 0.25–0.35 g cm^?3. Class (iii) is characterized by a partially saturated wetness profile, snowpack outflow values of 0.005–0.25 cm hr^?1, snowpack temperature near 0 °C, and average snow density between 0.25–0.40 g cm^?3. Class (iv) is characterized by strong diurnal snowpack outflow pattern with values as high as 0.75 cm hr^?1, stable snowpack temperature near 0 °C, and stable average snow density between 0.35–0.45 g cm^?3. Class (v) occurs intermittently between Classes (ii)–(iv) and displays low snowpack outflow values between 0.0008–0.04 cm hr^?1, a slight decrease in temperature relative to the preceding class, and similar densities to the preceding class. Numerical modeling of snowpack properties with SNOWPACK using both the Storage Threshold scheme and Richards' equation was used to quantify the effect of snowpack capillarity on predictions of snowpack outflow and other snowpack properties. Results indicate that both simulations are able to predict snow depth, snow temperature, and snow density reasonably well with little difference between the 2 water transport schemes. Richards' equation more accurately simulates the timing of snowpack outflow over the Storage Threshold scheme, especially early in the melt season and at diurnal timescales.     

THE WAVELENGTH COMPONENT VECTORGRAM: A TOOL FOR RESOLVING TWO-COMPONENT FLUORESCENT MIXTURES

We consider the problem of resolving the total luminescence spectra of two-component mixtures.Non-negativity of the constituent spectra is commonly used as a criterion to accomplish this resolution,but its application seldom leads to unique answers. Sometimes, the condition must be relaxed because nosolution satisfying non-negativity exists. In these cases judgments must be made concerning which partsof the spectra will be allowed to go negative and by how much. The wavelength component vectorgramn (WCV) is presented here as a tool for assisting theexperimenter in making these judgments. Even when non-negative resolutions exist, it may happen thatresolutions with slight negative elements will more closely match the spectra of the pure constituents.The WCV can also be used as a tool for detecting these cases. The methodology is illustrated with realdata from a two-component mixture.     

POROSITY AND PORE STRUCTURE FROM ACOUSTIC WELL LOGGING DATA1

Wyllie's time-average equation and subsequent refinements have been used for over 20 years to estimate the porosity of reservoir rocks from compressional (P)-wave velocity (or its reciprocal, transit time) recorded on a sonic log. This model, while simple, needs to be more convincingly explained in theory and improved in practice, particularly by making use of shear (S)-wave velocity. One of the most important, although often ignored, factors affecting elastic velocities in a rock is pore structure, which is also a controlling factor for transport properties of a rock. Now that S-wave information can be obtained from the sonic log, it may be used with P-waves to provide a better understanding of pore structure. A new acoustic velocities-to-porosity transform based on an elastic velocity model developed by Kuster and Toksöz is proposed. Employing an approximation to an equivalent pore aspect ratio spectrum, pore structure for reservoir rocks is taken into account, in addition to total pore volume. Equidimensional pores are approximated by spheres and rounded spheroids, while grain boundary pores and flat pores are approximated by low aspect ratio cracks. An equivalent pore aspect ratio spectrum is characterized by a power function which is determined by compressional-and shear-wave velocities, as well as by matrix and inclusion properties. As a result of this more sophisticated elastic model of porous rocks and a stricter theory of elastic wave propagation, the new method leads to a more satisfactory interpretation and fuller use of seismic and sonic log data. Calculations using the new transform on data for sedimentary rocks, obtained from published literature and laboratory measurements, are presented and compared at atmospheric pressure with those estimated from the time-average equation. Results demonstrate that, to compensate for additional complexity, the new method provides more detailed information on pore volume and pore structure of reservoir rocks. Examples are presented using a realistic self-consistent averaging scheme to consider interactions between pores, and the possibility of extending the method to complex lithologies and shaly rocks is discussed.     

The effects of topographically-controlled thermal tides in the martian upper atmosphere as seen by the MGS accelerometer

Mars Global Surveyor accelerometer observations of the martian upper atmosphere revealed large variations in density with longitude during northern hemisphere spring at altitudes of 130-160 km, all latitudes, and mid-afternoon local solar times (LSTs). This zonal structure is due to tides from the surface. The zonal structure is stable on timescales of weeks, decays with increasing altitude above 130 km, and is dominated by wave-3 (average amplitude 22% of mean density) and wave-2 (18%) harmonics. The phases of these harmonics are constant with both altitude and latitude, though their amplitudes change significantly with latitude. Near the South Pole, the phase of the wave-2 harmonic changes by 90° with a change of half a martian solar day while the wave-3 phase stays constant, suggesting diurnal and semidiurnal behaviour, respectively. We use a simple application of classical tidal theory to identify the dominant tidal modes and obtain results consistent with those of General Circulation Models. Our method is less rigorous, but simpler, than the General Circulation Models and hence complements them. Topography has a strong influence on the zonal structure.