Concentration–discharge relationships reflect chemostatic characteristics of US catchments

Concentration–discharge relationships have been widely used as clues to the hydrochemical processes that control runoff chemistry. Here we examine concentration–discharge relationships for solutes produced primarily by mineral weathering in 59 geochemically diverse US catchments. We show that these catchments exhibit nearly chemostatic behaviour; their stream concentrations of weathering products such as Ca, Mg, Na, and Si typically vary by factors of only 3 to 20 while discharge varies by several orders of magnitude. Similar patterns are observed at the inter‐annual time scale. This behaviour implies that solute concentrations in stream water are not determined by simple dilution of a fixed solute flux by a variable flux of water, and that rates of solute production and/or mobilization must be nearly proportional to water fluxes, both on storm and inter‐annual timescales. We compared these catchments' concentration–discharge relationships to the predictions of several simple hydrological and geochemical models. Most of these models can be forced to approximately fit the observed concentration–discharge relationships, but often only by assuming unrealistic or internally inconsistent parameter values. We propose a new model that also fits the data and may be more robust. We suggest possible tests of the new model for future studies. The relative stability of concentration under widely varying discharge may help make aquatic environments habitable. It also implies that fluxes of weathering solutes in streams, and thus fluxes of alkalinity to the oceans, are determined primarily by water fluxes. Thus, hydrology may be a major driver of the ocean‐alkalinity feedback regulating climate change. Copyright © 2009 John Wiley & Sons, Ltd.     

Biogeochemical impact of a model western iron source in the Pacific Equatorial Undercurrent

Trace element distributions in the source waters of the Pacific Equatorial Undercurrent (EUC) show the existence of elevated total acid-soluble iron concentrations. This region has been suggested to contribute enough bioavailable iron to regulate interannual and interglacial variability in biological productivity downstream in the high-nitrate low-chlorophyll upwelling zone of the eastern equatorial Pacific. We investigated the advection and first-order biogeochemical impact of an imposed, data-based iron maximum in the western Pacific EUC using an ecosystem model forced by a global dynamical model. We imposed two source profiles of iron constrained by total acid-soluble iron measurements. Though the data for total acid-soluble iron included both dissolved and acid-soluble particulate iron species, we treated all of the total acid-soluble iron as if it was dissolved and bioavailable. A deeper (270 m) source was centered in the density horizon of the observed iron maximum and a shallower (180 m) source was located in the core of our model's EUC, where a dissolved iron maximum has been frequently postulated. These source runs were compared with a control run that contained no specific source of iron associated with the EUC. In the source runs elevated iron concentrations were simulated in the EUC across its entire zonal span, evident as a subsurface plume of dissolved iron slightly below the core of the EUC. In the control run there was no iron maximum associated with the EUC. Upwelling of iron-replete water in the central and eastern equatorial Pacific increased integrated primary productivity in the Wyrtki box (180°W:90°W, 5°S:5°N, 0:200 m) by 41% and 66% for the deeper and shallower iron perturbation, respectively. The source runs increased the realism of the zonal extent of HNLC conditions and the meridional distributions of biological productivity, relative to the control run. However, in the source simulations surface chlorophyll concentrations were too high by a factor of two and maximum surface nitrate concentrations were too low, relative to climatologies. The relative abundance of diatoms roughly doubled upon the input of additional iron, exceeding field observations. Though biogeochemical data are limited and we did not adjust parameters to optimize the model fits to observations, these results suggest that acid-soluble particulate iron supplied to the EUC in the western equatorial Pacific is unlikely to be entirely bioavailable.     

Simulation of the Magnetothermal Instability

In many magnetized, dilute astrophysical plasmas, thermal conduction occurs almost exclusively parallel to magnetic field lines. In this case, the usual stability criterion for convective stability, the Schwarzschild criterion, which depends on entropy gradients, is modified. In the magnetized long mean free path regime, instability occurs for small wavenumbers when (∂ P/∂z) (∂ ln T/∂ z) > 0, which we refer to as the Balbus criterion. We refer to the convective-type instability that results as the magnetothermal instability (MTI). We use the equations of MHD with anisotropic electron heat conduction to numerically simulate the linear growth and nonlinear saturation of the MTI in plane-parallel atmospheres that are unstable according to the Balbus criterion. The linear growth rates measured from the simulations are in excellent agreement with the weak field dispersion relation. The addition of isotropic conduction, e.g. radiation, or strong magnetic fields can damp the growth of the MTI and affect the nonlinear regime. The instability saturates when the atmosphere becomes isothermal as the source of free energy is exhausted. By maintaining a fixed temperature difference between the top and bottom boundaries of the simulation domain, sustained convective turbulence can be driven. MTI-stable layers introduced by isotropic conduction are used to prevent the formation of unresolved, thermal boundary layers. We find that the largest component of the time-averaged heat flux is due to advective motions as opposed to the actual thermal conduction itself. Finally, we explore the implications of this instability for a variety of astrophysical systems, such as neutron stars, the hot intracluster medium of galaxy clusters, and the structure of radiatively inefficient accretion flows. J. M. Stone: Program in Applied and Computational Mathematics, Princeton University, Princeton, NJ 08544     

Application of the SHETRAN basin‐scale,landslide sediment yield model to the Llobregat basin,Spanish Pyrenees

The first application of the SHETRAN basin‐scale, landslide erosion and sediment yield model is carried out for a major landsliding event in the upper 505 km² of the Llobregat basin, in the eastern Spanish Pyrenees, in November 1982. The model simulates the spatial distribution of shallow landslides and their sediment yield. Acknowledging uncertainty in the model parameter evaluation, the aim of the application was not to reproduce the observed occurrence of landslides as accurately as possible with one simulation, but to bracket the observed pattern with several simulations representing uncertainty in the key input conditions. Bounds on the landslide simulations were thus determined as a function of uncertainty in the vegetation root cohesion (used in the model factor of safety calculations). The resulting upper bound considerably overestimates the observed pattern (17 000 landslides compared with an observation of around 700), but it reproduces several of the principal clusters in the observed pattern. The lower bound contains around 500 landslides. The sediment yield estimates (2670–14 630 t km^?2) are comparable to measurements elsewhere in the Pyrenees for extreme events. The results demonstrate an ability to simulate the basin‐scale landslide response to a rainfall event and the resulting sediment yield. They also highlight the need for further research in setting the uncertainty bounds and in avoiding large overestimates of landslide occurrence arising in part from a current inability to model small‐scale controls for a basin of the given size. Copyright © 2006 John Wiley & Sons, Ltd.     

Abundance and degree of residency of humpback dolphins Sousa plumbea in Mossel Bay,South Africa

Indian Ocean humpback dolphins Sousa plumbea inhabit nearshore waters from South Africa to eastern India. Humpback dolphins are vulnerable to conservation threats due to their naturally small population sizes and use of nearshore habitats, where human activities are highest. We investigated the abundance and residency of this species inhabiting Mossel Bay, South Africa, using photographic mark-recapture. Data were collected during 81 surveys in Mossel Bay between 2011 and 2013. Open population modelling using the POPAN parameterisation produced a ‘super-population’ estimate of 125 individuals (95% CI: 61–260) and within-year estimates of between 33 and 86 individuals (2011: 71 [95% CI: 30–168]; 2012: 33 [15–73], 32 [15–70]; 2013: 46 [20–108]). Although less appropriate, closed capture models were also run for comparison with previous studies in the region and generated similar, but slightly smaller, population estimates within each year. We compared our catalogue with opportunistic data collected from East London, Plettenberg Bay, De Hoop and Gansbaai. The only catalogue matches attained were between Plettenberg Bay (n = 44 identified) and Mossel Bay (n = 67 identified), separated by 140?km. Population exchange was moderate, with nine individuals resighted in multiple years between these two areas. This study supports previous findings of long-range movements for this species and provides a baseline from which to assess future impacts on the population.     

Cluster analysis on the bulk elemental compositions of Antarctic stony meteorites

Remote sensing observations by recent successful missions to small bodies have revealed the difficulty in classifying the materials which cover their surfaces into a conventional classification of meteorites. Although reflectance spectroscopy is a powerful tool for this purpose, it is influenced by many factors, such as space weathering, lighting conditions, and surface physical conditions (e.g., particle size and style of mixing). Thus, complementary information, such as elemental compositions, which can be obtained by X‐ray fluorescence (XRF) and gamma‐ray spectrometers (GRS), have been considered very important. However, classifying planetary materials solely based on elemental compositions has not been investigated extensively. In this study, we perform principal component and cluster analyses on 12 major and minor elements of the bulk compositions of 500 meteorites reported in the National Institute of Polar Research (NIPR), Japan database. Our unique approach, which includes using hierarchical cluster analysis, indicates that meteorites can be classified into about 10 groups purely by their bulk elemental compositions. We suggest that Si, Fe, Mg, Ca, and Na are the optimal set of elements, as this set has been used successfully to classify meteorites of the NIPR database with more than 94% accuracy. Principal components analysis indicates that elemental compositions of meteorites form eight clusters in the three‐dimensional space of the components. The three major principal components (PC1, PC2, and PC3) are interpreted as (1) degree of differentiations of the source body (i.e., primitive versus differentiated), (2) degree of thermal effects, and (3) degree of chemical fractionation, respectively.     

Secular tidal changes in lunar orbit and Earth rotation

Small tidal forces in the Earth–Moon system cause detectable changes in the orbit. Tidal energy dissipation causes secular rates in the lunar mean motion n, semimajor axis a, and eccentricity e. Terrestrial dissipation causes most of the tidal change in n and a, but lunar dissipation decreases eccentricity rate. Terrestrial tidal dissipation also slows the rotation of the Earth and increases obliquity. A tidal acceleration model is used for integration of the lunar orbit. Analysis of lunar laser ranging (LLR) data provides two or three terrestrial and two lunar dissipation parameters. Additional parameters come from geophysical knowledge of terrestrial tides. When those parameters are converted to secular rates for orbit elements, one obtains dn/dt = \(-25.97\pm 0.05 ''/\)cent\(^{2}\), da/dt = 38.30 ± 0.08 mm/year, and di/dt = ?0.5 ± 0.1 \(\upmu \)as/year. Solving for two terrestrial time delays and an extra de/dt from unspecified causes gives \(\sim \) \(3\times 10^{-12}\)/year for the latter; solving for three LLR tidal time delays without the extra de/dt gives a larger phase lag of the N2 tide so that total de/dt = \((1.50 \pm 0.10)\times 10^{-11}\)/year. For total dn/dt, there is \(\le \)1 % difference between geophysical models of average tidal dissipation in oceans and solid Earth and LLR results, and most of that difference comes from diurnal tides. The geophysical model predicts that tidal deceleration of Earth rotation is \(-1316 ''\)/cent\(^{2}\) or 87.5 s/cent\(^{2}\) for UT1-AT, a 2.395 ms/cent increase in the length of day, and an obliquity rate of 9 \(\upmu \)as/year. For evolution during past times of slow recession, the eccentricity rate can be negative.     

Evidence for high‐temperature fractionation of lithium isotopes during differentiation of the Moon

Lithium isotope and abundance data are reported for Apollo 15 and 17 mare basalts and the LaPaz low‐Ti mare basalt meteorites, along with lithium isotope data for carbonaceous, ordinary, and enstatite chondrites, and chondrules from the Allende CV3 meteorite. Apollo 15 low‐Ti mare basalts have lower Li contents and lower δ⁷Li (3.8 ± 1.2‰; all uncertainties are 2 standard deviations) than Apollo 17 high‐Ti mare basalts (δ⁷Li = 5.2 ± 1.2‰), with evolved LaPaz mare basalts having high Li contents, but similar low δ⁷Li (3.7 ± 0.5‰) to Apollo 15 mare basalts. In low‐Ti mare basalt 15555, the highest concentrations of Li occur in late‐stage tridymite (>20 ppm) and plagioclase (11 ± 3 ppm), with olivine (6.1 ± 3.8 ppm), pyroxene (4.2 ± 1.6 ppm), and ilmenite (0.8 ± 0.7 ppm) having lower Li concentrations. Values of δ⁷Li in low‐ and high‐Ti mare basalt sources broadly correlate negatively with ¹⁸O/¹⁶O and positively with ⁵⁶Fe/⁵⁴Fe (low‐Ti: δ⁷Li ≤4‰; δ⁵⁶Fe ≤0.04‰; δ¹⁸O ≥5.7‰; high‐Ti: δ⁷Li >6‰; δ⁵⁶Fe >0.18‰; δ¹⁸O <5.4‰). Lithium does not appear to have acted as a volatile element during planetary formation, with subequal Li contents in mare basalts compared with terrestrial, martian, or vestan basaltic rocks. Observed Li isotopic fractionations in mare basalts can potentially be explained through large‐degree, high‐temperature igneous differentiation of their source regions. Progressive magma ocean crystallization led to enrichment in Li and δ⁷Li in late‐stage liquids, probably as a consequence of preferential retention of ⁷Li and Li in the melt relative to crystallizing solids. Lithium isotopic fractionation has not been observed during extensive differentiation in terrestrial magmatic systems and may only be recognizable during extensive planetary magmatic differentiation under volatile‐poor conditions, as expected for the lunar magma ocean. Our new analyses of chondrites show that they have δ⁷Li ranging between ?2.5‰ and 4‰. The higher δ⁷Li in planetary basalts than in the compilation of chondrites (2.1 ± 1.3‰) demonstrates that differentiated planetary basalts are, on average, isotopically heavier than most chondrites.