Effects of changes in canopy interception on stream runoff response and recovery following clear-cutting of a Japanese coniferous forest in Fukuroyamasawa Experimental Watershed in Japan

Understanding changes in evapotranspiration during forest regrowth is essential to predict changes of stream runoff and recovery after forest cutting. Canopy interception (Ic) is an important component of evapotranspiration, however Ic changes and the impact on stream runoff during regrowth after cutting remains unclear due to limited observations. The objective of this study was to examine the effects of Ic changes on long-term stream runoff in a regrowth Japanese cedar and Japanese cypress forest following clear-cutting. This study was conducted in two 1-ha paired headwater catchments at Fukuroyamasawa Experimental Watershed in Japan. The catchments were 100% covered by Japanese coniferous plantation forest, one of which was 100% clear-cut in 1999 when the forest was 70 years old. In the treated catchment, annual runoff increased by 301 mm/year (14% of precipitation) the year following clear-cutting, and remained 185 mm/year (7.9% of precipitation) higher in the young regrowth forest for 12–14 years compared to the estimated runoff assuming no clear-cutting. The Ic change was −358 mm/year (17% of precipitation) after cutting and was −168 mm/year (6.7% of precipitation) in the 12–14 years old regrowth forest compared to the observed Ic during the pre-cutting period. Stream runoff increased in all seasons, and the Ic change was the main fraction of evapotranspiration change in all seasons throughout the observation period. These results suggest that the change in Ic accounted for most of the runoff response following forest cutting and the subsequent runoff recovery in this coniferous forest.     

Water quality data for Fall Creek,New York,USA: 1972–1995

The 33 086 ha mixed land use Fall Creek watershed in upstate New York is part of the Great Lakes drainage system. Results from more than 3500 water samples are available in a data set that compiles flow data and measurements of various water quality analytes collected between 1972 and 1995 in all seasons and under all flow regimes in Fall Creek and its tributaries. Data is freely accessible at https://ecommons.cornell.edu/handle/1813/8148 and includes measurements of suspended solids, pH, alkalinity, calcium, magnesium, potassium, sodium, chloride, nitrate nitrogen (NO3-N), sulphate sulphur (SO4-S), phosphorus (P) fractions molybdate reactive P (MRP) and total dissolved P (TDP), percent P in sediment, and ammonium nitrogen (NH4-N). Methods, sub-watershed areas, and coordinates for sampling sites are also included. The work represented in this data set has made important scientific contributions to understanding of hydrological and biogeochemical processes that influence loading in mixed use watersheds and that have an impact on algal productivity in receiving water bodies. In addition, the work has been foundational for important regulatory and management decisions in the region.     

Modelling the effect of changing precipitation inputs on deep soil water utilization

Forests in the Southeastern United States are predicted to experience future changes in seasonal patterns of precipitation inputs as well as more variable precipitation events. These climate change‐induced alterations could increase drought and lower soil water availability. Drought could alter rooting patterns and increase the importance of deep roots that access subsurface water resources. To address plant response to drought in both deep rooting and soil water utilization as well as soil drainage, we utilize a throughfall reduction experiment in a loblolly pine plantation of the Southeastern United States to calibrate and validate a hydrological model. The model was accurately calibrated against field measured soil moisture data under ambient rainfall and validated using 30% throughfall reduction data. Using this model, we then tested these scenarios: (a) evenly reduced precipitation; (b) less precipitation in summer, more in winter; (c) same total amount of precipitation with less frequent but heavier storms; and (d) shallower rooting depth under the above 3 scenarios. When less precipitation was received, drainage decreased proportionally much faster than evapotranspiration implying plants will acquire water first to the detriment of drainage. When precipitation was reduced by more than 30%, plants relied on stored soil water to satisfy evapotranspiration suggesting 30% may be a threshold that if sustained over the long term would deplete plant available soil water. Under the third scenario, evapotranspiration and drainage decreased, whereas surface run‐off increased. Changes in root biomass measured before and 4 years after the throughfall reduction experiment were not detected among treatments. Model simulations, however, indicated gains in evapotranspiration with deeper roots under evenly reduced precipitation and seasonal precipitation redistribution scenarios but not when precipitation frequency was adjusted. Deep soil and deep rooting can provide an important buffer capacity when precipitation alone cannot satisfy the evapotranspirational demand of forests. How this buffering capacity will persist in the face of changing precipitation inputs, however, will depend less on seasonal redistribution than on the magnitude of reductions and changes in rainfall frequency.     

Local ice caps in Finderup Land,North Greenland,survived the Holocene Thermal Maximum

Local glaciers and ice caps (GICs) comprise only ~5.4% of the total ice volume, but account for ~14–20% of the current ice loss in Greenland. The glacial history of GICs is not well constrained, however, and little is known about how they reacted to Holocene climate changes. Specifically, in North Greenland, there is limited knowledge about past GIC fluctuations and whether they survived the Holocene Thermal Maximum (HTM, ~8 to 5 ka). In this study, we use proglacial lake records to constrain the ice‐marginal fluctuations of three local ice caps in North Greenland including Flade Isblink, the largest ice cap in Greenland. Additionally, we have radiocarbon dated reworked marine molluscs in Little Ice Age (LIA) moraines adjacent to the Flade Isblink, which reveal when the ice cap was smaller than present. We found that outlet glaciers from Flade Isblink retreated inland of their present extent from ~9.4 to 0.2 cal. ka BP. The proglacial lake records, however, demonstrate that the lakes continued to receive glacial meltwater throughout the entire Holocene. This implies that GICs in Finderup Land survived the HTM. Our results are consistent with other observations from North Greenland but differ from locations in southern Greenland where all records show that the local ice caps at low and intermediate elevations disappeared completely during the HTM. We explain the north–south gradient in glacier response as a result of sensitivity to increased temperature and precipitation. While the increased temperatures during the HTM led to a complete melting of GICs in southern Greenland, GICs remained in North Greenland probably because the melting was counterbalanced by increased precipitation due to a reduction in Arctic sea‐ice extent and/or increased poleward moisture transport.     

Simulating the Propagation of Infrasonic Waves and Estimating the Energy of the Chelyabinsk Meteoroid Explosion Observed on February 15, 2013

Results obtained from simulating the propagation of infrasonic waves from the Chelyabinsk meteoroid explosion observed on February 15, 2013, are given. The pseudodifferential parabolic equation (PDPE) method has been used for calculations. Data on infrasonic waves recorded at the IS31 station (Aktyubinsk, Kazakhstan), located 542.7 km from the likely location of the explosion, have been analyzed. Six infrasonic arrivals (isolated clearly defined pulse signals) were recorded. It is shown that the first “fast” arrival (F) corresponds to the propagation of infrasound in a surface acoustic waveguide. The rest of the arrivals (T1–T5) are thermospheric. The agreement between the results of calculations based on the PDPE method and experimental data is satisfactory. The energy E of the explosion has been estimated using two methods. One of these methods is based on the law of conservation of the acoustic pulse I, which is a product of the wave profile area S/2 of the signal under analysis and the distance to its source E _I [kt] = 1.38 × 10^–10 (I [kg/s])^1.482. The other method is based on the relation between the energy of explosion and the dominant period T of recorded signal E _T [kt] = 1.02 × (T [s]²/σ)^3/2, where σ is the dimensionless distance determining the degree of nonlinear effects during the propagation of sound along ray trajectories. According to the data, the explosion energy E _I,T ranges from 1.87 to 32 kt TNT.     

Daily-Seasonal Variation in the Number of Remotely Observed,Pulsed, ULF Electromagnetic Earthquakes and the Spatial Distribution of Their Generation Zones on the Earth’s Surface

Geomagnetism and Aeronomy - Based on data from long-term observations at two geophysical observatories, Borok and College, distantly spaced in latitude and longitude, the results of remote...     

Wetland hydroperiod classification in the western prairies using multitemporal synthetic aperture radar

Wetlands represent one of the world's most biodiverse and threatened ecosystem types and were diminished globally by about two‐thirds in the 20th century. There is continuing decline in wetland quantity and function due to infilling and other human activities. In addition, with climate change, warmer temperatures and changes in precipitation and evapotranspiration are reducing wetland surface and groundwater supplies, further altering wetland hydrology and vegetation. There is a need to automate inventory and monitoring of wetlands, and as a study system, we investigated the Shepard Slough wetlands complex, which includes numerous wetlands in urban, suburban, and agricultural zones in the prairie pothole region of southern Alberta, Canada. Here, wetlands are generally confined to depressions in the undulating terrain, challenging wetlands inventory and monitoring. This study applied threshold and frequency analysis routines for high‐resolution, single‐polarization (HH) RADARSAT‐2, synthetic aperture radar mapping. This enabled a growing season surface water extent hyroperiod‐based wetland classification, which can support water and wetland resource monitoring. This 3‐year study demonstrated synthetic aperture radar‐derived multitemporal open‐water masks provided an effective index of wetland permanence class, with overall accuracies of 89% to 95% compared with optical validation data, and RMSE between 0.2 and 0.7 m between model and field validation data. This allowed for characterizing the distribution and dynamics of 4 marsh wetlands hydroperiod classes, temporary, seasonal, semipermanent, and permanent, and mapping of the sequential vegetation bands that included emergent, obligate wetland, facultative wetland, and upland plant communities. Hydroperiod variation and surface water extent were found to be influenced by short‐term rainfall events in both wet and dry years. Seasonal hydroperiods in wetlands were particularly variable if there was a decrease in the temporary or semipermanent hydroperiod classes. In years with extreme rain events, the temporary wetlands especially increased relative to longer lasting wetlands (84% in 2015 with significant rainfall events, compared with 42% otherwise).     

Analytical study of the swing-by maneuver in an elliptical system

The objective of the present paper is to derive a set of analytical equations that describe a swing-by maneuver realized in a system of primaries that are in elliptical orbits. The goal is to calculate the variations of energy, velocity and angular momentum as a function of the usual basic parameters that describe the swing-by maneuver, as done before for the case of circular orbits. In elliptical orbits the velocity of the secondary body is no longer constant, as in the circular case, but it varies with the position of the secondary body in its orbit. As a consequence, the variations of energy, velocity and angular momentum become functions of the magnitude and the angle between the velocity vector of the secondary body and the line connecting the primaries. The “patched-conics” approach is used to obtain these equations. The configurations that result in maximum gains and losses of energy for the spacecraft are shown next, and a comparison between the results obtained using the analytical equations and numerical simulations are made to validate the method developed here.     

Determination of Supermassive Black Hole Spins Based on the Standard Shakura–Sunyaev Accretion Disk Model and Polarimetric Observations

Based on spectropolarimetry for 47 type 1 active galactic nuclei observed with the 6-m BTA telescope, we have estimated the spins of the supermassive black holes at the centers of these galaxies. We have determined the spins based on the standard Shakura–Sunyaev accretion disk model. About 70% of the investigated active galactic nuclei are shown to have Kerr supermassive black holes with a dimensionless spin greater than 0.9.     

Critical Latitude in Tidal Dynamics Using the Kara Sea as an Example

It is well known that, within the linear nonviscous equations of tidal dynamics, the amplitudes of oscillations of the barotropic and baroclinic tidal velocity components unlimitedly increase when approaching the critical latitude. It is also known that the linear equations of tidal dynamics with a constant and specified vertical eddy viscosity indicate the occurrence of significant tidal velocity shears in the near-bottom layer, which are responsible for increasing the baroclinic tidal energy dissipation, the turbulent kinetic energy, and the thickness of the bottom boundary layer. The first circumstance—the growth of the amplitudes of oscillations of the barotropic and baroclinic tidal velocity components—is due to the elimination in the original equations of small terms, which are small everywhere except for the critical latitude zone. The second circumstance—the occurrence of significant tidal velocity shears—is due to the fact that internal tidal waves, which induce the dissipation of the baroclinic tidal energy and the diapycnal diffusion, are either not taken into account or described inadequately. It is suggested that diapycnal diffusion can lead to the degeneration (complete or partial) of tidal velocity shears, with all the ensuing consequences. The aforesaid is confirmed by simulation results obtained using the QUODDY-4 high-resolution three-dimensional finite-element hydrostatic model along the 66.25° E section, which passes in the Kara Sea across the critical latitude.