Reducing bias and quantifying uncertainty in watershed flux estimates: the R package loadflex

Many ecological insights into the function of rivers and watersheds emerge from quantifying the flux of solutes or suspended materials in rivers. Numerous methods for flux estimation have been described, and each has its strengths and weaknesses. Currently, the largest practical challenges in flux estimation are to select among these methods and to implement or apply whichever method is chosen. To ease this process of method selection and application, we have written an R software package called loadflex that implements several of the most popular methods for flux estimation, including regressions, interpolations, and the special case of interpolation known as the period-weighted approach. Our package also implements a lesser-known and empirically promising approach called the “composite method,” to which we have added an algorithm for estimating prediction uncertainty. Here we describe the structure and key features of loadflex, with a special emphasis on the rationale and details of our composite method implementation. We then demonstrate the use of loadflex by fitting four different models to nitrate data from the Lamprey River in southeastern New Hampshire, where two large floods in 2006–2007 are hypothesized to have driven a long-term shift in nitrate concentrations and fluxes from the watershed. The models each give believable estimates, and yet they yield different answers for whether and how the floods altered nitrate loads. In general, the best modeling approach for each new dataset will depend on the specific site and solute of interest, and researchers need to make an informed choice among the many possible models. Our package addresses this need by making it simple to apply and compare multiple load estimation models, ultimately allowing researchers to estimate riverine concentrations and fluxes with greater ease and accuracy

Complexity factor of spherically anisotropic polytropes from gravitational decoupling

In this work we will analyse the complexity factor, proposed by L. Herrera, of spherically symmetric static matter distributions satisfying a polytropic equation through the gravitational decoupling method. Specifically, we will use the 2-steps GD, which is a particular case of the Extended Geometric Deformation (EGD), to obtain analytic polytropic solutions of Einstein's equations. In order to give an example, we construct a model satisfying a polytropic equation of state using Tolman IV as seed solution

Hydrostatic Equilibrium of a Perfect Fluid Sphere with Exterior Higher-Dimensional Schwarzschild Spacetime

We discuss the question of how the number of dimensions of space and time can influence the equilibrium configurations of stars. We find that dimensionality does increase the effect of mass but not the contribution of the pressure, which is the same in any dimension. In the presence of a (positive) cosmological constant the condition of hydrostatic equilibrium imposes a lower limit on mass and matter density. We show how this limit depends on the number of dimensions and suggest that $\Lambda > 0$ is more effective in 4D than in higher dimensions. We obtain a general limit for the degree of compactification (gravitational potential on the boundary) of perfect fluid stars in $D$-dimensions. We argue that the effects of gravity are stronger in 4D than in any other number of dimensions. The generality of the results is also discussed

Definitive Evidence for Order-by-Quantum-Disorder in Er2Ti2O7

Here we establish the systematic existence of a U(1) degeneracy of all symmetry-allowed Hamiltonians quadratic in the spins on the pyrochlore lattice, at the mean-field level. By extracting the Hamiltonian of Er2Ti2O7 from inelastic neutron scattering measurements, we then show that the U(1)-degenerate states of Er2Ti2O7 are its classical ground states, and unambiguously show that quantum fluctuations break the degeneracy in a way which is confirmed by experiment. This is the first definitive observation of order-by-disorder in any material. We provide further verifiable consequences of this phenomenon, and several additional comparisons between theory and experiment.Comment: 4.5 pages, 3 figures, 7.5 pages of Supplemental Material, 8 supplemental figure

Mass and Charge in Brane-World and Non-Compact Kaluza-Klein Theories in 5 Dim

In classical Kaluza-Klein theory, with compactified extra dimensions and without scalar field, the rest mass as well as the electric charge of test particles are constants of motion. We show that in the case of a large extra dimension this is no longer so. We propose the Hamilton-Jacobi formalism, instead of the geodesic equation, for the study of test particles moving in a five-dimensional background metric. This formalism has a number of advantages: (i) it provides a clear and invariant definition of rest mass, without the ambiguities associated with the choice of the parameters used along the motion in 5D and 4D, (ii) the electromagnetic field can be easily incorporated in the discussion, and (iii) we avoid the difficulties associated with the "splitting" of the geodesic equation. For particles moving in a general 5D metric, we show how the effective rest mass, as measured by an observer in 4D, varies as a consequence of the large extra dimension. Also, the fifth component of the momentum changes along the motion. This component can be identified with the electric charge of test particles. With this interpretation, both the rest mass and the charge vary along the trajectory. The constant of motion is now a combination of these quantities. We study the cosmological variations of charge and rest mass in a five-dimensional bulk metric which is used to embed the standard k = 0 FRW universes. The time variations in the fine structure "constant" and the Thomson cross section are also discussed.Comment: V2: References added, discussion extended. V3 is identical to V2, references updated. To appear in General Relativity and Gravitatio

Accelerated expansion from braneworld models with variable vacuum energy

In braneworld models a variable vacuum energy may appear if the size of the extra dimension changes during the evolution of the universe. In this scenario the acceleration of the universe is related not only to the variation of the cosmological term, but also to the time evolution of $G$ and, possibly, to the variation of other fundamental "constants" as well. This is because the expansion rate of the extra dimension appears in different contexts, notably in expressions concerning the variation of rest mass and electric charge. We concentrate our attention on spatially-flat, homogeneous and isotropic, brane-universes where the matter density decreases as an inverse power of the scale factor, similar (but at different rate) to the power law in FRW-universes of general relativity. We show that these braneworld cosmologies are consistent with the observed accelerating universe and other observational requirements. In particular, $G$ becomes constant and $\Lambda_{(4)} \approx const \times H^2$ asymptotically in time. Another important feature is that the models contain no "adjustable" parameters. All the quantities, even the five-dimensional ones, can be evaluated by means of measurements in 4D. We provide precise constrains on the cosmological parameters and demonstrate that the "effective" equation of state of the universe can, in principle, be determined by measurements of the deceleration parameter alone. We give an explicit expression relating the density parameters $\Omega_{\rho}$, $\Omega_{\Lambda}$ and the deceleration parameter $q$. These results constitute concrete predictions that may help in observations for an experimental/observational test of the model.Comment: References added, typos correcte