An evaluation of ocean color model estimates of marine primary productivity in coastal and pelagic regions across the globe

Nearly half of the earth\u27s photosynthetically fixed carbon derives from the oceans. To determine global and region specific rates, we rely on models that estimate marine net primary productivity (NPP) thus it is essential that these models are evaluated to determine their accuracy. Here we assessed the skill of 21 ocean color models by comparing their estimates of depth-integrated NPP to 1156 in situ C-14 measurements encompassing ten marine regions including the Sargasso Sea, pelagic North Atlantic, coastal Northeast Atlantic, Black Sea, Mediterranean Sea, Arabian Sea, subtropical North Pacific, Ross Sea, West Antarctic Peninsula, and the Antarctic Polar Frontal Zone. Average model skill, as determined by root-mean square difference calculations, was lowest in the Black and Mediterranean Seas, highest in the pelagic North Atlantic and the Antarctic Polar Frontal Zone, and intermediate in the other six regions. The maximum fraction of model skill that may be attributable to uncertainties in both the input variables and in situ NPP measurements was nearly 72%. On average, the simplest depth/wavelength integrated models performed no worse than the more complex depth/wavelength resolved models. Ocean color models were not highly challenged in extreme conditions of surface chlorophyll-a and sea surface temperature, nor in high-nitrate low-chlorophyll waters. Water column depth was the primary influence on ocean color model performance such that average skill was significantly higher at depths greater than 250 m, suggesting that ocean color models are more challenged in Case-2 waters (coastal) than in Case-1 (pelagic) waters. Given that in situ chlorophyll-a data was used as input data, algorithm improvement is required to eliminate the poor performance of ocean color NPP models in Case-2 waters that are close to coastlines. Finally, ocean color chlorophyll-a algorithms are challenged by optically complex Case-2 waters, thus using satellite-derived chlorophyll-a to estimate NPP in coastal areas would likely further reduce the skill of ocean color models

Polynomials for Crystal Frameworks and the Rigid Unit Mode Spectrum

To each discrete translationally periodic bar-joint framework \C in \bR^d we associate a matrix-valued function \Phi_\C(z) defined on the d-torus. The rigid unit mode spectrum \Omega(\C) of \C is defined in terms of the multi-phases of phase-periodic infinitesimal flexes and is shown to correspond to the singular points of the function z \to \rank \Phi_\C(z) and also to the set of wave vectors of harmonic excitations which have vanishing energy in the long wavelength limit. To a crystal framework in Maxwell counting equilibrium, which corresponds to \Phi_\C(z) being square, the determinant of \Phi_\C(z) gives rise to a unique multi-variable polynomial p_\C(z_1,\dots,z_d). For ideal zeolites the algebraic variety of zeros of p_\C(z) on the d-torus coincides with the RUM spectrum. The matrix function is related to other aspects of idealised framework rigidity and flexibility and in particular leads to an explicit formula for the number of supercell-periodic floppy modes. In the case of certain zeolite frameworks in dimensions 2 and 3 direct proofs are given to show the maximal floppy mode property (order $N$). In particular this is the case for the cubic symmetry sodalite framework and some other idealised zeolites.Comment: Final version with new examples and figures, and with clearer streamlined proof

On the Theory of Superfluidity in Two Dimensions

The superfluid phase transition of the general vortex gas, in which the circulations may be any non-zero integer, is studied. When the net circulation of the system is not zero the absence of a superfluid phase is shown. When the net circulation of the vortices vanishes, the presence of off-diagonal long range order is demonstrated and the existence of an order parameter is proposed. The transition temperature for the general vortex gas is shown to be the Kosterlitz---Thouless temperature. An upper bound for the average vortex number density is established for the general vortex gas and an exact expression is derived for the Kosterlitz---Thouless ensemble.Comment: 22 pages, one figure, written in plain TeX, published in J. Phys. A24 (1991) 502

Three-periodic nets and tilings: natural tilings for nets

Rules for determining a unique natural tiling that carries a given three-periodic net as its 1-skeleton are presented and justified. A computer implementation of the rules and their application to tilings for zeolite nets and for the nets of the RCSR database are described

Streaming Algorithm for Euler Characteristic Curves of Multidimensional Images

We present an efficient algorithm to compute Euler characteristic curves of gray scale images of arbitrary dimension. In various applications the Euler characteristic curve is used as a descriptor of an image. Our algorithm is the first streaming algorithm for Euler characteristic curves. The usage of streaming removes the necessity to store the entire image in RAM. Experiments show that our implementation handles terabyte scale images on commodity hardware. Due to lock-free parallelism, it scales well with the number of processor cores. Our software---CHUNKYEuler---is available as open source on Bitbucket. Additionally, we put the concept of the Euler characteristic curve in the wider context of computational topology. In particular, we explain the connection with persistence diagrams

Smooth adiabatic evolutions with leaky power tails

Adiabatic evolutions with a gap condition have, under a range of circumstances, exponentially small tails that describe the leaking out of the spectral subspace. Adiabatic evolutions without a gap condition do not seem to have this feature in general. This is a known fact for eigenvalue crossing. We show that this is also the case for eigenvalues at the threshold of the continuous spectrum by considering the Friedrichs model.Comment: Final form, to appear in J. Phys. A; 11 pages, no figure

One-Dimensional Discrete Stark Hamiltonian and Resonance Scattering by Impurities

A one-dimensional discrete Stark Hamiltonian with a continuous electric field is constructed by extension theory methods. In absence of the impurities the model is proved to be exactly solvable, the spectrum is shown to be simple, continuous, filling the real axis; the eigenfunctions, the resolvent and the spectral measure are constructed explicitly. For this (unperturbed) system the resonance spectrum is shown to be empty. The model considering impurity in a single node is also constructed using the operator extension theory methods. The spectral analysis is performed and the dispersion equation for the resolvent singularities is obtained. The resonance spectrum is shown to contain infinite discrete set of resonances. One-to-one correspondence of the constructed Hamiltonian to some Lee-Friedrichs model is established.Comment: 20 pages, Latex, no figure

On Krein-like theorems for noncanonical Hamiltonian systems with continuous spectra: application to Vlasov-Poisson

The notions of spectral stability and the spectrum for the Vlasov-Poisson system linearized about homogeneous equilibria, f_0(v), are reviewed. Structural stability is reviewed and applied to perturbations of the linearized Vlasov operator through perturbations of f_0. We prove that for each f_0 there is an arbitrarily small delta f_0' in W^{1,1}(R) such that f_0+delta f_0$is unstable. When$f_0$ is perturbed by an area preserving rearrangement, f_0 will always be stable if the continuous spectrum is only of positive signature, where the signature of the continuous spectrum is defined as in previous work. If there is a signature change, then there is a rearrangement of f_0 that is unstable and arbitrarily close to f_0 with f_0' in W^{1,1}. This result is analogous to Krein's theorem for the continuous spectrum. We prove that if a discrete mode embedded in the continuous spectrum is surrounded by the opposite signature there is an infinitesimal perturbation in C^n norm that makes f_0 unstable. If f_0 is stable we prove that the signature of every discrete mode is the opposite of the continuum surrounding it.Comment: Submitted to the journal Transport Theory and Statistical Physics. 36 pages, 12 figure

Uniqueness in MHD in divergence form: right nullvectors and well-posedness

Magnetohydrodynamics in divergence form describes a hyperbolic system of covariant and constraint-free equations. It comprises a linear combination of an algebraic constraint and Faraday's equations. Here, we study the problem of well-posedness, and identify a preferred linear combination in this divergence formulation. The limit of weak magnetic fields shows the slow magnetosonic and Alfven waves to bifurcate from the contact discontinuity (entropy waves), while the fast magnetosonic wave is a regular perturbation of the hydrodynamical sound speed. These results are further reported as a starting point for characteristic based shock capturing schemes for simulations with ultra-relativistic shocks in magnetized relativistic fluids.Comment: To appear in J Math Phy

Synthesis and characterization of Nb2O5@C core-shell nanorods and Nb2O5nanorods by reacting Nb(OEt)5via RAPET (reaction under autogenic pressure at elevated temperatures) technique

The reaction of pentaethoxy niobate, Nb(OEt)5, at elevated temperature (800 °C) under autogenic pressure provides a chemical route to niobium oxide nanorods coated with amorphous carbon. This synthetic approach yielded nanocrystalline particles of Nb2O5@C. As prepared Nb2O5@C core-shell nanorods is annealed under air at 500 °C for 3 h (removing the carbon coating) results in neat Nb2O5nanorods. According to the TEM measurements, the Nb2O5crystals exhibit particle sizes between 25 nm and 100 nm, and the Nb2O5crystals display rod-like shapes without any indication of an amorphous character. The optical band gap of the Nb2O5nanorods was determined by diffuse reflectance spectroscopy (DRS) and was found to be 3.8 eV