Throughfall partitioning by trees

Although we know that rainfall interception (the rain caught, stored, and evaporated from aboveground vegetative surfaces and ground litter) is affected by rain and throughfall drop size, what was unknown until now is the relative proportion of each throughfall type (free throughfall, splash throughfall, canopy drip) beneath coniferous and broadleaved trees. Based on a multinational data set of >120 million throughfall drops, we found that the type, number, and volume of throughfall drops are different between coniferous and broadleaved tree species, leaf states, and timing within rain events. Compared with leafed broadleaved trees, conifers had a lower percentage of canopy drip (51% vs. 69% with respect to total throughfall volume) and slightly smaller diameter splash throughfall and canopy drip. Canopy drip from leafless broadleaved trees consisted of fewer and smaller diameter drops (D_{50_DR}, 50th cumulative drop volume percentile for canopy drip, of 2.24 mm) than leafed broadleaved trees (D_{50_DR} of 4.32 mm). Canopy drip was much larger in diameter under woody drip points (D_{50_DR} of 5.92 mm) than leafed broadleaved trees. Based on throughfall volume, the percentage of canopy drip was significantly different between conifers, leafed broadleaved trees, leafless broadleaved trees, and woody surface drip points (p ranged from <0.001 to 0.005). These findings are partly attributable to differences in canopy structure and plant surface characteristics between plant functional types and canopy state (leaf, leafless), among other factors. Hence, our results demonstrating the importance of drop‐size‐dependent partitioning between coniferous and broadleaved tree species could be useful to those requiring more detailed information on throughfall fluxes to the forest floor.     

Quantifying threats to groundwater resources in the Republic of Maldives Part II: Recovery from tsunami marine overwash events

Marine overwash events are among the most serious short‐term threats to groundwater supply of small coral islands. During such events, seawater can inundate small islands partially or completely, causing salinization of the aquifer. A comprehensive knowledge of freshwater lens recovery is essential for water planners on these islands. In this study, a numerical modelling approach is used to quantify recovery of the freshwater lens on 4 islands of the Maldives after a tsunami‐induced overwash event similar to that experienced from the Indian Ocean earthquake in December 2004. The islands vary in size (0.2 to 10.1 km²) and span the climatic regions of the Maldives. A tested 3‐dimensional SEAWAT groundwater model for each island is used to simulate the recovery process. Recharge rates from historical rainfall data and from global climate models are imposed on each island during the post‐overwash recovery period. The effect of groundwater pumping on lens recovery also is examined. Results show abrupt decrease in fresh groundwater volumes for each island, followed by recovery that is significantly influenced by island size and recharge patterns. Overall, salinization is more widespread on small islands (<1 km²), but recovery is more rapid than for large islands. Between 50% and 90% of lens recovery occurs after 2 years for small islands (<1 km²) whereas only 35% and 55% for large islands. Imposing pumping rates required to sustain the local population lengthened the recovery time between 5% and 15%, with smaller islands having the higher percentage. However, the governing factor on recovery time is the spatial extent of land surface inundation by the overwash event, with wave height and duration of the event having a negligible impact. A strong relationship exists between required recovery time and island surface area, thereby providing a method to determine recovery time for other atoll islands not investigated in this study with similar geologic structure. Our results can be used to aid in managing water resources during the post‐overwash period.     

The Wasatch Environmental Observatory: A mountain to urban research network in the semi-arid western US

The 2085 km² Jordan River Basin, and its seven sub-catchments draining the Central Wasatch Range immediately east of Salt Lake City, UT, are home to an array of hydrologic, atmospheric, climatic and chemical research infrastructure that collectively forms the Wasatch Environmental Observatory (WEO). WEO is geographically nested within a wildland to urban land-use gradient and built upon a strong foundation of over a century of discharge and climate records. A 2200 m gradient in elevation results in variable precipitation, temperature and vegetation patterns. Soil and subsurface structure reflect systematic variation in geology from granitic, intrusive to mixed sedimentary clastic across headwater catchments, all draining to the alluvial or colluvial sediments of the former Lake Bonneville. Winter snowfall and spring snowmelt control annual hydroclimate, rapid population growth dominates geographic change in lower elevations and urban gas and particle emissions contribute to episodes of severe air pollution in this closed-basin. Long-term hydroclimate observations across this diverse landscape provide the foundation for an expanding network of infrastructure in both montane and urban landscapes. Current infrastructure supports both basic and applied research in atmospheric chemistry, biogeochemistry, climate, ecology, hydrology, meteorology, resource management and urban redesign that is augmented through strong partnerships with cooperating agencies. These features allow WEO to serve as a unique natural laboratory for addressing research questions facing seasonally snow-covered, semi-arid regions in a rapidly changing world and an excellent facility for providing student education and research training.     

Tracer‐based evidence of heterogeneity in subsurface flow and storage within a boreal hillslope

Runoff from boreal hillslopes is often affected by distinct soil boundaries, including the frozen boundary and the organic‐mineral boundary (OMB), where highly porous and hydraulically conductive organic material overlies fine‐grained mineral soils. Viewed from the surface, ground cover appears as a patchwork on sub‐meter scales, with thick, moss mats interspersed with lichen‐covered, silty soils with gravel inclusions. We conducted a decameter‐scale subsurface tracer test on a boreal forest hillslope in interior Alaska to quantify locations and mechanisms of transport and storage in these soils, focusing on the OMB. A sodium bromide tracer was added as a slug addition to a pit and sampled at 40 down‐gradient wells, screened primarily at the OMB and within a 7 × 12 m well field. We maintained an elevated head in the injection pit for 8.5 hr to simulate a storm. Tracer breakthrough velocities ranged from <0.12 to 0.93 m hr^?1, with the highest velocities in lichen‐covered soils. After 12 hr and cessation of the elevated head, the tracer coalesced and was only detected in thick mosses at a trough in the OMB. By 24 hr, approximately 17% of the tracer mass could be accounted for. The majority of the mass loss occurred between 4 and 12 hr, while the tracer was in contact with lichen‐covered soils, which is consistent with tracer transport into deeper flow paths via preferential flow through discrete gravelly areas. Slow breakthroughs suggest that storage and exchange also occurred in shallow soils, likely related to saturation and drainage in fine‐grained mineral soils caused by the elevated hydraulic head. These findings highlight the complex nature of storage and transmission of water and solutes from boreal hillslopes to streams and are particularly relevant given rapid changes to boreal environments related to climate change, thawing permafrost and increasing fire severity.     

Relationships between riparian evapotranspiration and groundwater depth along a semiarid irrigated river valley

Evapotranspiration (ET) from riparian vegetation can be difficult to estimate due to relatively abundant water supply, spatial vegetation heterogeneity, and interactions with anthropogenic influences such as shallower groundwater tables, increased salinity, and nonpoint source pollution induced by irrigation. In semiarid south-eastern Colorado, reliable ET estimates are scarce for the riparian corridor that borders the Arkansas River. This work investigates relationships between the riparian ecosystem along the Arkansas River and an underlying alluvial aquifer using ET estimates from remotely sensed data and modelled water table depths. Results from a calibrated, finite-difference groundwater model are used to estimate weekly water table fluctuations in the riparian ecosystem from 1999 to 2009, and estimates of ET are calculated using the Operational Simplified Surface Energy Balance (SSEBop) model with over 200 Landsat scenes covering over 30 km² of riparian ecosystem along a 70-km stretch of the river. Comparison of calculated monthly SSEBop ET to estimated alfalfa reference ET from local micrometeorological station data indicated statistically significant high linear correspondence (R² = .87). Daily calculated SSEBop ET showed statistically significant moderate linear correspondence with data from a local weighing lysimeter (R² = .59). Simulated monthly SSEBop ET values were larger in drier years compared with wetter years, and ET variability was also larger in drier years. Peak ET most commonly occurred during the month of June for all 11 years of analysis. Relationships between ET and water table depth showed that peak monthly ET was highest when groundwater depths were less than about 3 m, and ET values were significantly lower for groundwater depths greater than 3 m. Negative sample Spearman correlation highlighted riparian corridor locations where ET increased as a result of decreased groundwater depths across years with different hydroclimatic conditions. This study shows how a combination of remotely sensed riparian ET estimates and a regional groundwater model can improve our understanding of linkages between riparian consumptive use and near-river groundwater conditions influenced by irrigation return flow and different climatic drivers.     

Assessing groundwater availability of the Maldives under future climate conditions

Groundwater resources of the Republic of the Maldives are threatened by a variety of factors including variable future rainfall patterns, continued population growth and associated pumping demands, rising sea level, and contamination from the land surface. This study assesses changes in groundwater availability due to variable rainfall patterns and sea level rise (SLR) in the coming decades, a key component of water resources management for the country. Using a suite of two‐dimensional density‐dependent groundwater flow models, time‐dependent thickness of the freshwater lens is simulated for a range of island sizes (200 to 1,100 m) during the time period of 2011 to 2050, with recharge to the freshwater lens calculated using rainfall patterns provided by general circulation models for the three distinct geographic regions of the Maldives. The effect of SLR on the freshwater lens is quantified using estimates of shoreline recession and associated decreases in island width. If rainfall is solely considered, groundwater availability is projected to increase, as lens thickness during the 2031–2050 time periods is slightly greater (1–5%) than during the 2011–2030 time period. However, including the impact of SLR indicates an overall decrease in lens thickness, with drastic decreases (60% to 100%) projected for small islands (200 m) and moderate decreases (12% to 14%) expected for 400 m islands, which accommodate one third of the national population. Similar methodologies can be used for other atoll island nations, such as the Republic of Marshall Islands, Federated States of Micronesia, and the Republic of Kiribati. For the Maldives, results from this study can be used in conjunction with population growth estimates to determine the feasibility of including groundwater in water resources planning and management for the country.     

Galaxy formation in the reionization epoch as hinted by Wide Field Camera 3 observations of the Hubble Ultra Deep Field

We present a large sample of candidate galaxies at z ≈ 7-10,selected in the Hubble Ultra Deep Field using the new observations of the Wide Field Camera 3that was recently installed on the Hubble Space Telescope.Our sample is composed of 20 z850-dropouts(four new discoveries),15 Y105-dropouts(nine new discoveries)and 20 J125-dropouts(all new discoveries).The surface densities of the z850-dropouts are close to what was predicted by earlier studies,however,those of the Y105-and J125-dropouts are quite unexpected.While no Y105-or J125-dropouts have been found at AB < 28.0 mag,their surface densities seem to increase sharply at fainter levels.While some of these candidates seem to be close to foreground galaxies and thus could possibly be gravitationally lensed,the overall surface densities after excluding such cases are still much higher than what would be expected if the luminosity function does not evolve from z～7 to 10.Motivated by such steep increases,we tentatively propose a set of Schechter function parameters to describe the luminosity functions at z ≈ 8 and 10.As compared to their counterpart at z ≈ 7,here L*decreases by a factor demanded by the existing observations,they are allowed and seem to agree with the data better than other alternatives.If these luminosity functions are still valid beyond our current detection limit,this would imply a sudden emergence of a large number of low-luminosity galaxies when looking back in time to z ≈ 10,which,while seemingly exotic,would naturally fit in the picture of the cosmic hydrogen reionization.These early galaxies could easily account for the ionizing photon budget required by the reionization,and they would imply that the global star formation rate density might start from a very high value at z ≈ 10,rapidly reach the minimum at z ≈ 7,and start to rise again towards z ≈ 6.In this scenario,the majority of the stellar mass that the universe assembled through the reionization epoch seems still undetected by current observations at z ≈ 6.     

Star clusters around recoiled black holes in the Milky Way halo

Gravitational wave emission by coalescing black holes (BHs) kicks the remnant BH with a typical velocity of hundreds of  km s⁻¹  . This velocity is sufficiently large to remove the remnant BH from a low-mass galaxy but is below the escape velocity from the Milky Way (MW) galaxy. If central BHs were common in the galactic building blocks that merged to make the MW, then numerous BHs that were kicked out of low-mass galaxies should be freely floating in the MW halo today. We use a large statistical sample of possible merger tree histories for the MW to estimate the expected number of recoiled BH remnants present in the MW halo today. We find that hundreds of BHs should remain bound to the MW halo after leaving their parent low-mass galaxies. Each BH carries a compact cluster of old stars that populated the core of its original host galaxy. Using the time-dependent Fokker–Planck equation, we find that the present-day clusters are  ≲1 pc  in size, and their central bright regions should be unresolved in most existing sky surveys. These compact systems are distinguishable from globular clusters by their internal (Keplerian) velocity dispersion greater than 100 km s⁻¹ and their high mass-to-light ratio owing to the central BH. An observational discovery of this relic population of star clusters in the MW halo would constrain the formation history of the MW and the dynamics of BH mergers in the early Universe. A similar population should exist around other galaxies and may potentially be detectable in M31 and M33.     

Gravitational waves from scattering of stellar-mass black holes in galactic nuclei

Stellar-mass black holes (BHs) are expected to segregate and form a steep density cusp around supermassive black holes (SMBHs) in galactic nuclei. We follow the evolution of a multimass system of BHs and stars by numerically integrating the Fokker–Planck energy diffusion equations for a variety of BH mass distributions. We find that the BHs 'self-segregate', and that the rarest, most massive BHs dominate the scattering rate closest to the SMBH  (≲10⁻¹ pc)  . BH–BH binaries form out of gravitational wave emission during BH encounters. We find that the expected rate of BH coalescence events detectable by Advanced LIGO is  ∼1–10² yr⁻¹  , depending on the initial mass function of stars in galactic nuclei and the mass of the most massive BHs. We find that the actual merger rate is likely ∼10 times larger than this due to the intrinsic scatter of stellar densities in many different galaxies. The BH binaries that form this way in galactic nuclei have significant eccentricities as they enter the LIGO band (90 per cent with   e > 0.9  ), and are therefore distinguishable from other binaries, which circularize before becoming detectable. We also show that eccentric mergers can be detected to larger distances and greater BH masses than circular mergers, up to  ∼700 M_⊙  . Future ground-based gravitational wave observatories will be able to constrain both the mass function of BHs and stars in galactic nuclei.     

Compositional and lithological diversity among brecciated lunar meteorites of intermediate iron concentration

Abstract— We present new compositional data for 30 lunar stones representing about 19 meteorites. Most have iron concentrations intermediate to those of the numerous feldspathic lunar meteorites (3–7% FeO) and the basaltic lunar meteorites (17–23% FeO). All but one are polymict breccias. Some, as implied by their intermediate composition, are mainly mixtures of brecciated anorthosite and mare basalt, with low concentrations of incompatible elements such as Sm (1–3 μg/g). These breccias likely originate from points on the Moon where mare basalt has mixed with material of the FHT (Feldspathic Highlands Terrane). Others, however, are not anorthosite‐basalt mixtures. Three (17–75 μ/g Sm) consist mainly of nonmare mafic material from the nearside PKT (Procellarum KREEP Terrane) and a few are ternary mixtures of material from the FHT, PKT, and maria. Some contain mafic, nonmare lithologies like anorthositic norites, norites, gabbronorites, and troctolite. These breccias are largely unlike breccias of the Apollo collection in that they are poor in Sm as well as highly feldspathic anorthosite such as that common at the Apollo 16 site. Several have high Th/Sm compared to Apollo breccias. Dhofar 961, which is olivine gabbronoritic and moderately rich in Sm, has lower Eu/Sm than Apollo samples of similar Sm concentration. This difference indicates that the carrier of rare earth elements is not KREEP, as known from the Apollo missions. On the basis of our present knowledge from remote sensing, among lunar meteorites Dhofar 961 is the one most likely to have originated from South Pole‐Aitken basin on the lunar far side.