Multichannel Three-Dimensional SOLA Inversion for Local Helioseismology

Inversions for local helioseismology are an important and necessary step for obtaining three-dimensional maps of various physical quantities in the solar interior. Frequently, the full inverse problems that one would like to solve prove intractable because of computational constraints. Due to the enormous seismic data sets that already exist and those forthcoming, this is a problem that needs to be addressed. To this end, we present a very efficient linear inversion algorithm for local helioseismology. It is based on a subtractive optimally localized averaging (SOLA) scheme in the Fourier domain, utilizing the horizontal-translation invariance of the sensitivity kernels. In Fourier space the problem decouples into many small problems, one for each horizontal wave vector. This multichannel SOLA method is demonstrated for an example problem in time–distance helioseismology that is small enough to be solved both in real and Fourier space. We find that both approaches are successful in solving the inverse problem. However, the multichannel SOLA algorithm is much faster and can easily be parallelized.     

Sedimentology and geochemistry of an exceptionally preserved last interglacial sapropel S5 in the Levantine Basin (Mediterranean Sea)

A combined study of lithological, geochemical and physical sediment properties is reported from a completely laminated S5 sapropel, recovered in three gravity cores (M40-4 SL67, M51-3 SL103, M51-3 SL104) from the Pliny Trench region of the eastern Mediterranean. The thickness of the studied sapropel S5 varies between 85 and 91 cm and tops most S5-sapropels in the Mediterranean. Based on optical features like color and thickness of laminae, the sapropels were subdivided into thirteen distinct lithostratigraphic zones. These zones, as well as the finer layering pattern within them, could be followed exactly among the three cores, indicating that the processes responsible for this variation acted at least on a regional scale. The sapropel sediment is characterized by exceptionally high porosity, which is strongly correlated with Si/Ca. This relationship implies that the sapropel is in essence an organic-matter rich diatomite and its exceptional thickness can be explained by preservation of diatoms forming a loosely packed sediment fabric. Compared to other S5 sapropels, the preservation of diatoms has apparently led to a twofold increase in the thickness of the sapropel layer. Relative abundances of 10 elements were determined at ultra-high resolution (0.2 mm) by XRF-scanner over the complete length of each sapropel including several cm of enclosing marl. An analysis of the chemical data indicates that the lowermost 13 cm of the sapropel is chemically more similar to the underlying marl and that the sediment chemistry shows different signals at different scales. The strongest pattern is the contrast between the sapropel and the surrounding marl, which is accentuated in elements indicative for redox conditions as well as terrigenous sediment input and productivity. Within the sapropel, a mm- to cm-scale layering is observed. The abundances of many elements are systematically linked to the pattern of these layers, indicating a common origin, related to productivity and/or terrigenous sediment and/or redox conditions. This pattern indicates a link to a regional climatic process, making the S5 sapropel horizon in M40-4 SL67, M51-3 SL103 and M51-3 SL104 a potential high-resolution archive of climatic variability during the last interglacial in the Mediterranean Sea and its adjacent landmasses.     

Impact of reservoirs and channelization on lowland river macroinvertebrates: A case study from Central Europe

In order to assess and compare the ecological impacts of channelization and shallow lowland reservoirs, macroinvertebrate communities of a lowland metapotamal river below reservoirs with epilimnial release were studied. The study was carried out in the Dyje River (Czech Republic) at five sites located from 1.5 to 22.5 km downstream of the reservoir outfall. The five sites differed in the degree of channel modification from natural muddy banks to riprap regulation. Seven samples were collected during the years 1998 and 1999 at each site using a semiquantitative method. The data were processed using multivariate analyses and methods for assessing the ecological and functional structure of communities. Altogether, 261 species of benthic macroinvertebrates were recorded including several rare and threatened taxa. Based on the results of principal component analysis (PCA), most of the variability within the species data (the first PCA axis) was explained by the degree of channel modification, from natural muddy banks with aquatic vegetation to a man-made riprap. The second axis was strongly correlated with current velocity. The sites differed in species richness, total abundances, proportion of individual functional feeding groups, pattern of the distribution of the current preference groups, and values of several biotic indexes, all of which also corresponded to the degree of channel modification. Thus, the morphological man-made modifications of the river channel were found to be the main factor affecting lowland river macroinvertebrates and their biodiversity. Our results suggest that the biggest threat to benthic macroinvertebrate diversity of lowland rivers comes from channelization. The impact of reservoirs can be completely overwhelmed by the impact of channelization, especially when muddy banks with aquatic vegetation present a substantial part of habitat diversity and significantly contribute to the total species pool.     

Statistical characteristics of detectable inhomogeneities in observed meteorological time series

Statistical characteristics of detectable inhomogeneities [IHs] in more than 600 observed meteorological time series have been investigated using 16 objective homogenisation methods. Forty and 100 year long series of monthly or annual characteristics of surface air temperature, precipitation total and relative air humidity from the Czech Republic and Hungary were examined. The area of the part of the Czech observing network used here is smaller, and the density of sites is larger, than in the Hungarian network, resulting in higher spatial correlations among data in the Czech dataset relative to the Hungarian dataset. Time series with low number of gaps were supplied with interpolated data. Before homogenisation relative time series were created, using weighted averages of time series from the same geographical region as reference series. For ease of comparison, the magnitudes of the detected IHs are normalised with the standard deviation of the noise in the relative time series. Results show that observed meteorological time series usually contain large number of small IHs, and that the magnitude distribution of IHs from different data segments are surprisingly similar. Effects of different spatial coherences on the results are discussed.     

Geological interpretation of gravity profiles through the Karlovy Vary granite massif (Czech Republic)

We examined the shape of the Late Variscan Karlovy Vary granite massif located south of the Ohre/Eger graben in Northern Bohemia by reinterpretation of existing gravity data on two perpendicular profiles. The granite body of about 360 km² total outcrop size has the elongation ratio 0.35 with the major axis trending NE-SW. The SW part of the body was crossed in the nineties by the seismic profile 9HR which localized the bottom of granites in a depth of about 10 km. We used this value as a reference datum in our gravity profiles. We positioned one of our profiles along the seismic profile 9HR and the other one perpendicularly, i.e. parallel with the elongation of the outcrop surface. We interpret the shape of the main granite body in the vicinity of Karlovy Vary as a continuous desk whose floor is horizontal (or subhorizontal) and varies along its whole extension about a depth of 10 km. This thickness is approximately identical with that of the Saxothuringian nappes imaged by seismic reflection. The near surface upper contact of the granite body is mildly inclined, and outward dipping. It changes to steep sides or inward inclined contacts in deeper levels. The Lesny-Lysina (Kynžvart) massif is a separate granite body about 324 km thick, not continuously connected with the main Karlovy Vary massif. The gravity curve suggests that granites often enclose in their endocontact large blocks of country metasediments or metabasites the existence of which is partly evidenced by their outcrops outside the line of the profile. The granite body is found density-homogenous. Minor density differences between granite varieties are caused mainly by more intense hydrothermal alterations in younger suite granites. We interpret vertical conduits for the ascent of granitic magmas to be parallel to the Jáchymov-Gera and Ohře (Eger) lineaments or the Mariánské Lázně fault zone as indicated by the elongation of some outcrops. However, they are not clearly imaged from the gravity data. The effect of the depression of the Sokolov basin along the faults parallel with the Ohře (Eger) lineament is shallow and it is not indicated by any change in the floor depth of the granite body. Comparison of the seismicity distribution suggests that the hypocenters occur mostly outside of the granite bodies or near their contact with the country rock.     

Bedload hysteresis in a glacier‐fed mountain river

Sediment transport during flood events often reveals hysteretic patterns because flow discharge can peak before (counterclockwise hysteresis) or after (clockwise hysteresis) the peak of bedload. Hysteresis in sediment transport has been used in the literature to infer the degree of sediment availability. Counterclockwise and clockwise hysteresis have been in fact interpreted as limited and unlimited sediment supply conditions, respectively. Hysteresis has been mainly explored for the case of suspended sediment transport, but it was rarely reported for bedload transport in mountain streams. This work focuses on the temporal variability of bedload transport in an alpine catchment (Saldur basin, 18.6 km², Italian Alps) where bedload transport was monitored by means of an acoustic pipe sensor which detects the acoustic vibrations induced by particles hitting a 0.5m‐long steel pipe. Runoff dynamics are dominated by snowmelt in late spring/early summer, mostly by glacier melt in late summer/early autumn, and by a combination of the snow and glacier melt in mid‐summer. The results indicate that hysteretic patterns during daily discharge fluctuations are predominantly clockwise during the snowmelt period, likely due to the ready availability of unpacked sediments within the channel or through bank erosion in the lower part of the basin. On the contrary, counterclockwise hysteresis tend to be more frequent during late glacier melting period, possibly due to the time lag needed for sediment provided by the glacial and peri‐glacial area to be transported to the monitoring section. However, intense rainfall events occurring during the glacier melt period generated predominantly clockwise hysteresis, thus indicating the activation of different sediment sources. These results indicate that runoff generation processes play a crucial role on sediment supply and temporal availability in mountain streams. Copyright © 2014 John Wiley & Sons, Ltd.     

Method for prediction of landslide movements based on random forests

Prediction of landslide movements with practical application for landslide risk mitigation is a challenge for scientists. This study presents a methodology for prediction of landslide movements using random forests, a machine learning algorithm based on regression trees. The prediction method was established based on a time series consisting of 2 years of data on landslide movement, groundwater level, and precipitation gathered from the Kostanjek landslide monitoring system and nearby meteorological stations in Zagreb (Croatia). Because of complex relations between precipitations and groundwater levels, the process of landslide movement prediction is divided into two separate models: (1) model for prediction of groundwater levels from precipitation data and (2) model for prediction of landslide movements from groundwater level data. In a groundwater level prediction model, 75 parameters were used as predictors, calculated from precipitation and evapotranspiration data. In the landslide movement prediction model, 10 parameters calculated from groundwater level data were used as predictors. Model validation was performed through the prediction of groundwater levels and prediction of landslide movements for the periods from 10 to 90 days. The validation results show the capability of the model to predict the evolution of daily displacements, from predicted variations of groundwater levels, for the period up to 30 days. Practical contributions of the developed method include the possibility of automated predictions, updated and improved on a daily basis, which would be an important source of information for decisions related to crisis management in the case of risky landslide movements.     

Assessment of potential sediment contamination using screening methods (XRF,TGA/MS) taking into account principles of green chemistry,Eastern Slovakia

Mining and milling of metal ores coupled with industries have bequeathed many countries the legacy of wide distribution of metal contaminants in sediments. The aim of this study was to assess potential sediment contamination via useful screening methods (XRF, CHNS, TGA/MS). The sediments were collected from the water reservoir Krompachy Eastern Slovakia in April 2015. Within the frame of evaluation it was found that the concentrations of the study elements (Cu, Zn, As, Pb, Cr, Ni, Cd) exceeded some of the MPC, TV and IV values. Sample c was the most polluted by metals, which evident according to it’s the highest CHNS proportion as well as the highest clay and silt proportion. In the samples studied the best correlation was confirmed between weight losses in the temperature range (400–620 °C) and the following metal concentrations: Cu (r = 0.89), Zn (r = 0.88), As (r = 0.93), Hg (r = 0.83), Pb (r = 0.87). The greatest proportions of m/z 44, m/z 18 were detected at temperatures (400–620 °C) associated with decomposition of minerals such as siderite, barite, and exothermic loss of more refractory aromatic C took also place.     

Co‐evolution of riparian vegetation and channel dynamics in an aggrading braided river system,Mount Pinatubo,Philippines

Increased bank stability by riparian vegetation can have profound impacts on channel morphology and dynamics in low‐energy systems, but the effects are less clear in high‐energy environments. Here we investigate the role of vegetation in active, aggrading braided systems at Mount Pinatubo, Philippines, and compare results with numerical modeling results. Gradual reductions in post‐eruption sediment loads have reduced bed reworking rates, allowing vegetation to finally persist year‐round on the Pasig‐Potrero and Sacobia Rivers. From 2009–2011 we collected data detailing vegetation extent, type, density, and root strength. Incorporating these data into the RipRoot model and BSTEM (Bank Stability and Toe Erosion Model) shows cohesion due to roots increases from zero in unvegetated conditions to > 10·2 kPa in densely‐growing grasses. Field‐based parameters were incorporated into a cellular model comparing vegetation strength and sediment mobility effects on braided channel dynamics. The model shows both low sediment mobility and high vegetation strength lead to less active systems, reflecting trends observed in the field. The competing influence of vegetation strength versus channel dynamics is a concept encapsulated in a dimensionless ratio between timescales for vegetation growth and channel reworking known as T*. An estimated T* between 1·5 and 2·3 for the Pasig‐Potrero River suggests channels are still very mobile and likely to remain braided until aggradation rates decline further. Vegetation does have an important effect on channel dynamics, however, by focusing flow and thus aggradation into the unvegetated fraction of braidplain, leading to an aggradational imbalance and transition to a more avulsive state. The future trajectory of channel–vegetation interactions as sedimentation rates decline is complicated by strong seasonal variability in precipitation and sediment loads, driving incision and armoring in the dry season. By 2011, incision during the dry season was substantial enough to lower the water‐table, weaken existing vegetation, and allow for vegetation removal in future avulsions. Copyright © 2015 John Wiley & Sons, Ltd.