Mid and late Holocene forest fires and deforestation in the subalpine belt of the Iberian range,northern Spain

The conversion of subalpine forests into grasslands for pastoral use is a well-knownphenomenon, although for most mountain areas the timing of deforestation has not been determined. The presence of charcoal fragments in soil profiles affected by shallow landsliding enabled us to date the occurrence of fires and the periods of conversion ofsubalpine forest into grasslands in the Urbión Mountains, Iberian Range, Spain. We found that the treeline in the highest parts of the northwestern massifs of the Iberian Range(the Urbión, Demanda, Neila, and Cebollera massifs) is currently between 1500 and 1600 m a.s.l., probably because of pastoral use of the subalpine belt, whereas in the past it would have reached almost the highest divides(at approximately 2100–2200 m a.s.l.). The radiocarbon dates obtained indicate that the transformation of the subalpine belt occurred during the Late Neolithic, Chalcolithic, Bronze Age, Iron Age, and Middle Ages. Forest clearing was probably moderate during fires prior to the Middle Ages, as the small size of the sheep herds and the local character of the markets only required small clearings, and therefore more limited fires. Thus, it is likely that the forest recovered burnt areas in a few decades; this suggests the management of the forest and grasslands following a slash-andburn system. During the Middle and Modern Ages deforestation and grassland expansion affected most of the subalpine belt and coincided with the increasing prevalence of transhumance, as occurred in other mountains in the Iberian Peninsula(particularly the Pyrenees). Although the occurrence of shallow landslides following deforestation between the Neolithic and the Roman Period cannot be ruled out, the most extensive shallow landsliding processes would have occurred from the Middle Ages until recent times.     

Distributed Source Coding of Hyperspectral Images Based on Three-Dimensional Wavelet

To reduce the possibility of poor efficiency and weak anti-error capability while encoding and transmitting hyperspectral images, we present a distributed source coding scheme for hyperspectral images based on three-dimensional (3D) set partitioning in hierarchical trees (SPIHT). First, the 3D wavelet transform is performed on the hyperspectral image. Thereafter, the low frequency section is regarded as the Key frame and the high frequency section as the Wyner–Ziv frame to enable independent SPIHT coding through different transmission channels. The Wyner–Ziv encoder uses Turbo channel coding to create high frequency information that reflects the details of the image with better anti-error capacity, while the low frequency information shows the main energy of the image. In this study, we used SPIHT coding to acquire a bit stream with quality scalability. Results show that the proposed scheme is more efficient during coding, while at the same time providing improved anti-error capability and quality scalability of the bit stream.     

Robotic telescopes for high energy astrophysics in Ondřejov

We report on two small aperture robotic telescopes called BART and D50 operated in Ondřejov. Both telescopes are capable of automatic observation of gamma ray burst (GRB) optical afterglows. Coordinates of GRBs are taken from alerts distributed via Internet. Telescopes observe other interesting high energy sources when there is not any alert. The smaller telescope BART has aperture D = 254 mm. The bigger telescope D50 has a primary mirror of diameter D = 500 mm. Both telescopes are controlled by free software package RTS2 and are accessible through Internet. We describe the two telescopes and related software and show some results such as our first observed optical counterpart of GRB.     

AMR Simulations of Magnetohydrodynamic Problems by the CESE Method in Curvilinear Coordinates

The objective of this paper is to present new extensions of the space – time conservation element and solution element (CESE) method for simulations of magnetohydrodynamic (MHD) problems in general curvilinear coordinates by using an adaptive mesh refinement (AMR) grid system. By transforming the governing MHD equations from the physical space (x,y,z) to the computational space (ξ,η,ζ) while retaining the form of conservation, the CESE method is established for MHD in the curvilinear coordinates. Utilizing the parallel AMR package PARAMESH, we present the first implementation of applying the AMR CESE method for MHD (AMR-CESE-MHD) in both Cartesian and curvilinear coordinates. To show the validity and capabilities of the AMR-CESE-MHD code, a suite of numerical tests in two and three dimensions including ideal MHD and resistive MHD are carried out, with two of them in both Cartesian and curvilinear coordinates. Numerical tests show that our results are highly consistent with those obtained previously by other authors, and the results under both coordinate systems confirm each other very well.     

A Multiple Flare Scenario where the Classic Long-Duration Flare Was Not the Source of a CME

A series of flares (GOES class M, M and C) and a CME were observed in close succession on 20 January 2004 in NOAA 10540. Radio observations, which took the form of types II, III and N bursts, were associated with these events. We use the combined observations from TRACE, EIT, Hα images from Kwasan, MDI magnetograms and GOES to understand the complex development of this event. Contrary to a standard interpretation, we conclude that the first two impulsive flares are part of the CME launch process while the following long-duration event flare represents simply the recovery phase. Observations show that the flare ribbons not only separate but also shift along the magnetic inversion line so that magnetic reconnection progresses stepwise to neighboring flux tubes. We conclude that “tether cutting” reconnection in the sheared arcade progressively transforms it to a twisted flux tube, which becomes unstable, leading to a CME. We interpret the third flare, a long-duration event, as a combination of the classical two-ribbon flare with the relaxation process following forced reconnection between the expanding CME structure and neighboring magnetic fields. Electronic Supplementary Material  The online version of this article () contains supplementary material, which is available to authorized users.     

Estimating satellite phase fractional cycle biases based on Kalman filter

Phase fractional cycle biases (FCBs) originating from satellites and receivers destroy the integer nature of PPP carrier phase ambiguities. To achieve integer ambiguity resolution of PPP, FCBs of satellites are required. In former work, least squares methods are commonly adopted to isolate FCBs from a network of reference stations. However, it can be extremely time consuming concerning the large number of observations from hundreds of stations and thousands of epochs. In addition, iterations are required to deal with the one-cycle inconsistency among FCB measurements. We propose to estimate the FCB based on a Kalman filter. The large number of observations are handled epoch by epoch, which significantly reduces the dimension of the involved matrix and accelerates the computation. In addition, it is also suitable for real-time applications. As for the one-cycle inconsistency, a pre-elimination method is developed to avoid iterations and posterior adjustments. A globally distributed network consisting of about 200 IGS stations is selected to determine the GPS satellite FCBs. Observations recorded from DoY 52 to 61 in 2016 are processed to verify the proposed approach. The RMS of wide lane (WL) posterior residuals is 0.09 cycles while that of the narrow lane (NL) is about 0.05 cycles, which indicates a good internal accuracy. The estimated WL FCBs also have a good consistency with existing WL FCB products (e.g., CNES-GRG, WHU-SGG). The RMS of differences with respect to GRG and SGG products are 0.03 and 0.05 cycles. For satellite NL FCB estimates, 97.9% of the differences with respect to SGG products are within ±?0.1 cycles. The RMS of the difference is 0.05 cycles. These results prove the efficiency of the proposed approach.     

Do Spur-Throated Grasshoppers, <Emphasis Type="Italic">Melanoplus</Emphasis> spp. (Orthoptera: Acrididae), Exert Top-Down Control on Smooth Cordgrass <Emphasis Type="Italic">Spartina alterniflora</Emphasis> in Northern New England?

Recently, strong top-down (consumer) control of cordgrass (Spartina alterniflora) has been demonstrated. Here, we manipulated the densities of cordgrass consumers, acridid grasshoppers (Melanoplus bivittatus and Melanoplus femurrubrum), to examine their impact on cordgrass in the Plum Island Estuary (PIE), MA, USA. After 1 month, there was no detectable effect of grasshopper density on S. alterniflora biomass and grasshoppers at the highest densities (34 individuals per square meter) consumed only ~14% of the standing stock biomass. However, significant impacts of grasshopper density on grazing damage were seen. For example, plant damage and scarring length increased by 160% and 6,156%, respectively, at the highest grasshopper densities relative to exclusion (zero grasshoppers) densities. Plant height was significantly reduced with increasing grasshopper densities, although this may be a function of leaf tip removal instead of reduced plant growth. No other strong consumers of cordgrass (e.g., Littoraria irrorata, Prokelisia marginata) were observed in PIE and we suggest that consumer regulation of cordgrass is weak in this system.     

The activity of silica in kimberlites,revisited

The activity of silica in a silicate liquid in equilibrium with olivine and orthopyroxene decreases with increasing pressure. In contrast, the activity of silica in an unbuffered silicate liquid changes little with pressure. Although the implications of these pressure dependencies have been considered by previous authors in terms of inferring pressures of origin of magmas, less consideration has been given to the implications of these dependencies on the evolution of the magma en route to the surface, or to the mantle through which the magma passes. In this paper, a combination of Schreinemakers’ analysis in isothermal section and calculated reactions in space is used to (a) rationalize the absence of orthopyroxene xenocrysts in kimberlites and the relative abundance of olivine “megacrysts” therein, (b) propose another reason for the paucity of xenocrystic mantle-derived carbonates in kimberlites, (c) explain why clinopyroxene is much less reactive in the kimberlite melt than is orthopyroxene, and (d) explore the implications of the relative stabilities of olivine, orthopyroxene, and clinopyroxene in kimberlitic magma for the mantle through which the magma transits.

---

Laboratory investigation and simulation of breakthrough curves in karst conduits with pools

A series of laboratory experiments are performed under various hydrological conditions to analyze the effect of pools in pipes on breakthrough curves (BTCs). The BTCs are generated after instantaneous injections of NaCl tracer solution. In order to test the feasibility of reproducing the BTCs and obtain transport parameters, three modeling approaches have been applied: the equilibrium model, the linear graphical method and the two-region nonequilibrium model. The investigation results show that pools induce tailing of the BTCs, and the shapes of BTCs depend on pool geometries and hydrological conditions. The simulations reveal that the two-region nonequilibrium model yields the best fits to experimental BTCs because the model can describe the transient storage in pools by the partition coefficient and the mass transfer coefficient. The model parameters indicate that pools produce high dispersion. The increased tailing occurs mainly because the partition coefficient decreases, as the number of pools increases. When comparing the tracer BTCs obtained using the two types of pools with the same size, the more appreciable BTC tails that occur for symmetrical pools likely result mainly from the less intense exchange between the water in the pools and the water in the pipe, because the partition coefficients for the two types of pools are virtually identical. Dispersivity values decrease as flow rates increase; however, the trend in dispersion is not clear. The reduced tailing is attributed to a decrease in immobile water with increasing flow rate. It provides evidence for hydrodynamically controlled tailing effects.     

Stream Restoration Performance and Its Contribution to the Chesapeake Bay TMDL: Challenges Posed by Climate Change in Urban Areas

In large part driven by total maximum daily load (TMDL) mandates, the restoration community in the Chesapeake Bay region has been implementing novel best management practices (BMPs) and stream restoration designs in urban areas, such as regenerative stream/stormwater conveyance (RSC) structures and stream-wetland complexes (SWCs). However, the nutrient and sediment reduction efficiencies of these novel designs are virtually unknown, and the possibility of increasing riverine flow in the Chesapeake Bay watershed associated with climate change this century necessitates an evaluation of their performance to develop and utilize those that optimize reductions in nutrient and sediment fluxes. We compare pre- and post-construction loads (total nitrogen (TN), total phosphorus (TP), and total suspended solids (TSS)) from RSCs (i.e., upland BMPs) and a SWC (i.e., stream restoration) constructed at the outflow of a highly developed watershed in the Coastal Plain physiographic province of Maryland. The largest of the two RSCs performed best in relation to expected nutrient and sediment reductions because of superior water retention capability. By the length of river reach restored, the SWC attained from 79 to 88% of its N reduction TMDL goal, but only 19 to 23 and 2.7 to 3.1% for TP and TSS, respectively; by watershed area, % attainments of TMDL goals were much lower. Results indicate that SWCs have the potential to curtail N loading from developed catchments, but additional water quality benefits may be limited. Climate change projections indicate that there will be an increased frequency of larger-volume storms that will result in an increase in stormflow runoff from urban areas, and increased pollutant loads will likely curtail potential gains made by efforts to achieve TMDL goals. Given the large-scale implementation of BMPs currently underway to accommodate the Chesapeake Bay TMDL, the restoration community needs to adopt a concerted strategy of building climate resilience into many types of urban BMPs to help attain and maintain loads at TMDL levels in anticipation of a progressively wetter climate throughout this century.