Geochemical correlations between effusive and explosive silicic volcanics in the Saraykent region (Yozgat), central Anatolia,Turkey

Two main volcanic events are distinguished between Saraykent and Akçakışla in the Yozgat province of central Anatolia: (1) early Late Cretaceous–Palaeocene effusive activity, that produced a sequence of intermediate to felsic ‘basal lavas’; and (2) marginally later Palaeocene explosive activity that formed a series of covering ignimbrite flows. Due to their close temporal and spatial relation, geochemical comparisons were made between the silicic members of the lavas and ignimbrites, to identify chemical groups and their relative petrogenesis. The basal lavas range from calc‐alkaline basaltic andesites to dominant rhyolites. Based on trace element correlations three main geochemical groups were identified: the Akçakışla rhyolites (present as domes); Akçakışla rhyodacites‐dacites (lava flows); and Ozan‐Saraykent rhyolites (lava flows). Large‐ion lithophile elements have been mobile in all the groups, but mainly in the Akçakışla rhyolites. Rare earth element (REE) patterns show marked similarity between the Ozan and Saraykent basal lavas. The Akçakışla dome rhyolites are more fractionated with lower La_N/Yb_N ratios (c.10), whereas the Akçakışla basal lavas have much higher La_N/Yb_N ratios (c.30). The chemical coherence and petrographic similarities between the Saraykent and Ozan lavas suggest a single suite related via fractionation. Three geochemical groups were also established for the ignimbrites: Saraykent ignimbrite; Bağlıca ignimbrite‐Toklu‐Kızıldağ crystal tuffs; and Keklikpınar ignimbrite. The ignimbrites, like the basal lavas, display a pronounced depletion in Ba on ORG‐normalized plots. Relative to the basal lavas, chondrite‐normalized patterns for the ignimbrites are different in displaying negative Eu anomalies that indicate feldspar fractionation. The lack of geochemical overlap or coherence between any of the lava and ignimbrite groups suggests that they represent distinct eruptive events and are not related in any simple volcanic development and cogenetic sense. Two geochemical features are common to all the volcanic rock groups: (1) the presence of a Nb‐Ta anomaly, which is generally accepted as a crustal signature; and (2) the relatively low Y abundances which appear characteristic for the region as a whole. These fundamental features of the local silicic volcanism largely reflect source composition and effects. Copyright © 2002 John Wiley & Sons, Ltd.     

Solar Spectral Irradiance Variations in 240 – 1600 nm During the Recent Solar Cycles 21 – 23

Regular solar spectral irradiance (SSI) observations from space that simultaneously cover the UV, visible (vis), and the near-IR (NIR) spectral region began with SCIAMACHY aboard ENVISAT in August 2002. Up to now, these direct observations cover less than a decade. In order for these SSI measurements to be useful in assessing the role of the Sun in climate change, records covering more than an eleven-year solar cycle are required. By using our recently developed empirical SCIA proxy model, we reconstruct daily SSI values over several decades by using solar proxies scaled to short-term SCIAMACHY solar irradiance observations to describe decadal irradiance changes. These calculations are compared to existing solar data: the UV data from SUSIM/UARS, from the DeLand & Cebula satellite composite, and the SIP model (S2K+VUV2002); and UV-vis-IR data from the NRLSSI and SATIRE models, and SIM/SORCE measurements. The mean SSI of the latter models show good agreement (less than 5%) in the vis regions over three decades while larger disagreements (10 – 20%) are found in the UV and IR regions. Between minima and maxima of Solar Cycles 21, 22, and 23, the inferred SSI variability from the SCIA proxy is intermediate between SATIRE and NRLSSI in the UV. While the DeLand & Cebula composite provide the highest variability between solar minimum and maximum, the SIP/Solar2000 and NRLSSI models show minimum variability, which may be due to the use of a single proxy in the modeling of the irradiances. In the vis-IR spectral region, the SCIA proxy model reports lower values in the changes from solar maximum to minimum, which may be attributed to overestimations of the sunspot proxy used in modeling the SCIAMACHY irradiances. The fairly short timeseries of SIM/SORCE shows a steeper decreasing (increasing) trend in the UV (vis) than the other data during the descending phase of Solar Cycle 23. Though considered to be only provisional, the opposite trend seen in the visible SIM data challenges the validity of proxy-based linear extrapolation commonly used in reconstructing past irradiances.     

The Interaction Between Coronal Mass Ejections and Streamers: A Statistical View over 15 Years (1996 – 2010)

We report on the statistical analysis of the interaction between coronal mass ejections (CMEs) and streamers based on 15 years (from 1996 to 2010 inclusive) of observation of the solar corona with the LASCO-C2 coronagraph. We used synoptic maps and improved the method of analysis of past investigations by implementing an automatic detection of both CMEs and streamers. We identified five categories of interaction based on photometric and geometric variations between the pre- and post-CME streamers: “brightening”, “dimming”, “emergence”, “disappearance”, and “deviation”. A sixth category, “no change”, included all cases where none of the above variations is observed. A “global set” of 21?242 CMEs was considered as well as a subset of the 10 % brightest CMEs (denoted “top-ten”) and three typical periods of solar activity: minimum, intermediate, and maximum. We found that about half of the global population of CMEs are not associated with streamers, whereas 93 % of the 10 % brightest CMEs are associated. When there is a CME-streamer association, approximately 95 % of the streamers experience a change, either geometric or radiometric. The “no change” category therefore amounts to approximately 5 %, but this percentage varies from 1?–?2 % during minimum to 7?–?8 % during intermediate periods of activity; values of 3?–?5 % are recorded during maximum. Emergences and disappearances of streamers are not dominant processes; they constitute 16?–?17 % of the global set and 23 % (emergence) and 28 % (disappearance) of the “top-ten” set. Streamer deviations are observed for 57 % and 70 % of, respectively, the global set and “top-ten” CMEs. The cases of dimming and brightening are roughly equally present and each case constitutes approximately 30?–?35 % of either set, global or “top-ten”.     

The Solar Ultraviolet Spectrum Estimated Using the Mg <Emphasis Type="SmallCaps">ii</Emphasis> Index and Ca <Emphasis Type="SmallCaps">ii</Emphasis> K Disk Activity

As part of a program to estimate the solar spectrum back to the early twentieth century, we have generated fits to UV spectral irradiance measurements from 1 – 410 nm. The longer wavelength spectra (150 – 410 nm) were fit as a function of two solar activity proxies, the Mg ii core-to-wing ratio, or Mg ii index, and the total Ca ii K disk activity derived from ground based observations. Irradiance spectra at shorter wavelengths (1 – 150 nm) where used to generate fits to the Mg ii core-to-wing ratio alone. Two sets of spectra were used in these fitting procedures. The fits at longer wavelengths (150 to 410 nm) were derived from the high-resolution spectra taken by the Solar Ultraviolet Spectral Irradiance Monitor (SUSIM) on the Upper Atmospheric Research Satellite (UARS). Spectra measured by the Solar EUV Experiment (SEE) instrument on the Thermosphere Ionosphere Mesosphere Energetics and Dynamics (TIMED) satellite were used for the fits at wavelengths from 1 to 150 nm. To generate fits between solar irradiance and solar proxies, this study uses the above irradiance data, the NOAA composite Mg ii index, and daily Ca ii K disk activity determined from images measured by Big Bear Solar Observatory (BBSO). In addition to the fitting coefficients between irradiance and solar proxies, other results from this study include an estimated relationship between the fraction of the disk with enhanced Ca ii K activity and the Mg ii index, an upper bound of the average solar UV spectral irradiance during periods where the solar disk contains only regions of the quiet Sun, as was believed to be present during the Maunder Minimum, as well as results indicating that slightly more than 60% of the total solar irradiance (TSI) variability occurs between 150 and 400 nm.     

Vesta,Vestoids, and the howardite,eucrite, diogenite group: Relationships and the origin of spectral differences

Abstract— Spectra of asteroid 4 Vesta and 21 small (estimated diameters less than 10 km) asteroids with Vesta‐like spectral properties (Vestoids) were measured at visible and near‐infrared wavelengths (～0.44 to ～1.65 μm). All of the measured small asteroids (except for 2579 Spartacus) have reflectance spectra consistent with surface compositions similar to eucrites and howardites and consistent with all being derived from Vesta. None of the observed asteroids have spectra similar to diogenites. We find no spectral distinction between the 15 objects tabulated as members of the Vesta dynamical family and 6 of the 7 sampled “non‐family” members that reside just outside the semi‐major axis (a), eccentricity (e), and inclination (i) region of the family. The spectral consistency and close orbital (a‐e‐i) match of these “non‐family” objects to Vesta and the Vesta family imply that the true bounds of the family extend beyond the subjective cut‐off for membership. Asteroid 2579 Spartacus has a spectrum consistent with a mixture of eucritic material and olivine. Spartacus could contain olivine‐rich material from Vesta's mantle or may be unrelated to Vesta altogether. Laboratory measurements of the spectra of eucrites show that samples having nearly identical compositions can display a wide range of spectral slopes. Finer particle sizes lead to an increase in the slope, which is usually referred to as reddening. This range of spectral variation for the best‐known meteoritic analogs to the Vestoids, regardless of whether they are actually related to each other, suggests that the extremely red spectral slopes for some Vestoids can be explained by very fine‐grained eucritic material on their surfaces.     

Stratigraphy and geochemistry of pillow basalts within the ophiolitic mélange of the Izmir–Ankara–Erzincan suture zone: implications for the geotectonic character of the northern branch of Neotethys

The Izmir-Ankara-Erzincan suture zone of Turkey is a broad zone of ophiolitic mélange containing numerous fragmented blocks ranging in age from Triassic to Cretaceous. Stratigraphic sequences for various mélange units are compared, together with the geochemistry of associated basaltic pillow lavas of Cretaceous age. A review of geochemical data for the pillow lavas demonstrate: (a) a dominant group of alkalic basalts with enriched incompatible elements, variable Zr/Y and Zr/Nb ratios, and (b) a range of tholeiitic basalts with slightly depleted to mildly enriched (normalized) rare earth patterns, (La/Yb)_N 0.4-3.0, and generally low Zr/Y ratios. The alkalic basalts can be chemically matched to Pacific Ocean seamounts, although the close association of red radiolarites and cherts suggests that many basalts represent the margins of such structures, rather than the main seamount edifice. Nd-Sr isotope data are typical for ocean island basalts and represent an admixture of a dominant EM-1 source and a depleted MORB-like source. Enhanced δ¹⁸O compositions are a consequence of submarine alteration and not crustal contamination. Tholeiitic compositions have affinities with both N- and E-type MORB, although most are probably representative of tholeiitic ocean islands. Overall the basalts are mainly representative of structures built on the ocean floor, rather than the oceanic crust itself, being scraped off the subducting crust and preserved in the mélange of the accretionary wedge.     

Comparing single-phase,non-Newtonian approaches with experimental results: Validating flume-scale mud and debris flow in HEC-RAS

Post-fire debris flows and tailing impoundment failures destroy lives and property. These geologic hazards – and other similar processes – fall on a continuum between classic Newtonian flood analyses and geotechnical stability analyses. The US Army Corps of Engineers (USACE) is developing a non-Newtonian library (DebrisLib) that includes a suite of rheological and clastic approaches to hyper-concentrated, mudflow and debris flow dynamics. The Hydrologic Engineering Center (HEC) has implemented these non-Newtonian methods into the widely used, public-domain open-channel hydraulics and morphodynamic software, HEC-RAS (river analysis system). This work presents part of the verification and validation of these non-Newtonian approaches, applying several rheological equations to published laboratory results high-concentration flume experiments. This study tested the linear Bingham model as well as the turbulent and Bagnold quadratic terms of the O'Brien equation. HEC-RAS also includes the non-linear Herschel–Bulkley (HB) approach, which quantifies shear-thickening or shear-thinning processes. The study used these non-Newtonian models in HEC-RAS to simulate 10 of Parsons et al.’s (2001) flume experiments, which measured the snout and plug velocity of fluids with high solid concentrations (C_v = 68–74%) and a broad range of material gradations (d₅₀ = 0.05–1 mm, d₁₅ = 0.006–0.1 mm). The experiments also measured and back-calculated Bingham and HB parameters of the materials, finding HB powers between 0.45 and 1.25 (i.e. fluids that are dilatant, pseudo-plastic and visco-plastic). The rheological models incorporated into DebrisLib and implemented in HEC-RAS reproduce experimental data well for most experiments. The Bingham model generated a plug velocity root-mean-square error (RMSE) of 0.21 m/s using standard flow parameters and Parsons et al.’s calibrated parameters, a substantial improvement over the unmodified shallow water flow equations (RMSE 0.77 m/s). Experiments with strong snout effects tended to generate higher residuals, especially in the snout velocity. The RMSE associated with the O'Brien equation was larger with the Parsons et al. fit parameters, but similar (0.23 m/s) with measured parameters. The turbulent parameter was the largest (often the dominant) parameter in most O'Brien simulations, with the dispersive stress only proving significant for the coarsest material. DebrisLib had to use a modified version of the dispersive term to simulate these concentrations. Both the 2D depth-averaged shallow water equation (SWE) and diffusion wave equation were used to simulate the experiments. The best results were obtained with the SWE with horizontal mixing.     

Joule heating of Io's ionosphere by unipolar induction currents

The unexpectedly large scale height of Io's ionosphere (Kliore, A., et al., 1975, Icarus24, 407–410) together with the relatively large molecular weight of the likely principal constituent, SO₂ (Pearl, J., et al., 1979, Nature280, 755–758), suggest a high ionospheric temperature. Electrical induction in Io's ionosphere due to the corotating plasma bound to the Jovian magnetosphere is one possible source for attainment of such high temperatures. Accordingly, unipolar induction models were constructed to calculate ionospheric joule heating numerically. Heating rates produced by highly simplified models lie in the range 10^?9 to 10^?8 W/m³. These heating rates are lower than those determined from uv photodissociative heating models (Kumar, S., 1980, Geophys. Res. Lett.7, 9–12) at low levels in the ionosphere but are comparable in the upper ionosphere. The low electrical heating rate throughout most of the ionosphere is due to the power limitation imposed by the Alfvén wings which complete the electrical circuit (Neubauer, F.M., 1980, J. Geophys. Res.85, 1171–1178). Contrary to the pre-Voyager calculations of Cloutier, P. A., et al. (1978, Astrophys. Space Sci.55, 93–112), our numerical results show that the J × B force density due to unipolar induction currents in the ionosphere is much less than the gravitational force density when the combined mass of the neutral species is included. The binding and coupling of the ionosphere is principally due to the relatively dense (possibly localized) neutral SO₂ atmosphere. In regions where the ions and neutrals are collisionally coupled the ionosphere will not be stripped off by the J × B forces. However at a level above that (to which the ions move by diffusion only) the charged species would be removed. Thus there appears to be no need to postulate the existence of an intrinsic Ionian magnetic field as suggested by Kivelson, M. G., et al. (79, Science 205, 491–493) and Southwood, S. J., et al. (1980, J. Geophys. Res., in press) in order to retain the observed ionosphere.     

Reach‐scale contributions of road‐surface sediment to the Honna River,Haida Gwaii,BC

Gravel road surfaces can be a major source of fine sediment to streams, yet their contribution to channel reach sediment balances remains poorly documented. To quantify the input of road surface material and to compare this input with natural sediment sources at the reach scale, suspended sediment dynamics was examined and a 16‐month sediment balance was developed for a ~35 channel‐width (approx. 425 m) reach of the Honna River, a medium‐size, road‐affected stream located in coastal British Columbia. Of the 105 ± 33 t of suspended material passing through the reach, 18 ± 6% was attributed to the road surface. The high availability of sediment on the road surface appears to limit hysteresis in road run‐off. During rainstorms that increase streamflow, road surface material composed 0.5–15% of sediment inputs during relatively dry conditions from April to the end of September and 5–70% through wetter conditions from October to the end of March, but our data do not show evidence of major sediment accumulation on the riverbed in the reach. A comparison of modelled sediment production on the road surface with observed yields from drainage channels suggests that (1) during low intensity rainfall, ditches and drainage channels may trap sediment from road run‐off, which is subsequently released during events of greater intensity, and/or (2) production models do not effectively describe processes, such as deposition or erosion of sediment in ditches, which control sediment transport and delivery. Our findings further emphasize the risk of unpaved roads in polluting river systems and highlight the continued need for careful road design and location away from sensitive aquatic environments. Copyright © 2016 John Wiley & Sons, Ltd.