Modeling ocean primary production: sensitivity to spectral resolution of attenuation and absorption of light

Modeling the vertical penetration of photosynthetically active radiation (PAR) through the ocean, and its utilization by phytoplankton, is fundamental to simulating marine primary production. The variation of attenuation and absorption of light with wavelength suggests that photosynthesis should be modeled at high spectral resolution, but this is computationally expensive. To model primary production in global 3d models, a balance between computer time and accuracy is necessary. We investigate the effects of varying the spectral resolution of the underwater light field and the photosynthetic efficiency of phytoplankton (α∗), on primary production using a 1d coupled ecosystem ocean turbulence model. The model is applied at three sites in the Atlantic Ocean (CIS (∼60°N), PAP (∼50°N) and ESTOC (∼30°N)) to include the effect of different meteorological forcing and parameter sets. We also investigate three different methods for modeling α∗ – as a fixed constant, varying with both wavelength and chlorophyll concentration [Bricaud, A., Morel, A., Babin, M., Allali, K., Claustre, H., 1998. Variations of light absorption by suspended particles with chlorophyll a concentration in oceanic (case 1) waters. Analysis and implications for bio-optical models. J. Geophys. Res. 103, 31033–31044], and using a non-spectral parameterization [Anderson, T.R., 1993. A spectrally averaged model of light penetration and photosynthesis. Limnol. Oceanogr. 38, 1403–1419]. After selecting the appropriate ecosystem parameters for each of the three sites we vary the spectral resolution of light and α∗ from 1 to 61 wavebands and study the results in conjunction with the three different α∗ estimation methods. The results show modeled estimates of ocean primary productivity are highly sensitive to the degree of spectral resolution and α∗. For accurate simulations of primary production and chlorophyll distribution we recommend a spectral resolution of at least six wavebands if α∗ is a function of wavelength and chlorophyll, and three wavebands if α∗ is a fixed value.

Predicting spectral features in galaxy spectra from broad-band photometry

We explore the prospects of predicting emission-line features present in galaxy spectra given broad-band photometry alone. There is a general consent that colours, and spectral features, most notably the 4000 angstrom break, can predict many properties of galaxies, including star formation rates and hence they could infer some of the line properties. We argue that these techniques have great prospects in helping us understand line emission in extragalactic objects and might speed up future galaxy redshift surveys if they are to target emission-line objects only. We use two independent methods, Artificial Neural Networks (based on the ANNz code) and Locally Weighted Regression (LWR), to retrieve correlations present in the colour N-dimensional space and to predict the equivalent widths present in the corresponding spectra. We also investigate how well it is possible to separate galaxies with and without lines from broad-band photometry only. We find, unsurprisingly, that recombination lines can be well predicted by galaxy colours. However, among collisional lines some can and some cannot be predicted well from galaxy colours alone, without any further redshift information. We also use our techniques to estimate how much information contained in spectral diagnostic diagrams can be recovered from broad-band photometry alone. We find that it is possible to classify active galactic nuclei and star formation objects relatively well using colours only. We suggest that this technique could be used to considerably improve redshift surveys such as the upcoming Fibre Multi Object Spectrograph (FMOS) survey and the planned Wide Field Multi Object Spectrograph (WFMOS) survey.

Partial spectral projected gradient method with active-set strategy for linearly constrained optimization

A method for linearly constrained optimization which modifies and generalizes recent box-constraint optimization algorithms is introduced. The new algorithm is based on a relaxed form of Spectral Projected Gradient iterations. Intercalated with these projected steps, internal iterations restricted to faces of the polytope are performed, which enhance the efficiency of the algorithm. Convergence proofs are given and numerical experiments are included and commented. Software supporting this paper is available through the Tango Project web page: http://www.ime.usp.br/similar to egbirgin/tango/.

Study of Blood Porphyrin Spectral Profile for Diagnosis of Chronic Renal Failure

The progression to end-stage renal failure is independent of the initial pathogenic mechanism. Metabolic acidosis is a common consequence of chronic renal failure that results from inadequate ammonium excretion and decreased tubular bicarbonate reabsorption. Protoporphyrin IX (PpIX) is the immediate metabolic precursor of the heme molecule. The purpose of this study was to evaluate the levels of erythrocytes protoporphyrin IX at an animal model during progressive renal disease. A total of 36 eight-week-old male Wistar rats were divided into six groups: Normal, 4 and 8 weeks after 5/6 nephrectomy (NX). Renal function was evaluated by creatinine clearance and plasma creatinine levels. The autofluorescence of erythrocytes porphyrin of healthy and NX rats was analyzed using fluorescence spectroscopy. Emission spectra were obtained by exciting the samples at 405 nm. Significant differences between normal and NX rats autofluorescence shape occurred in the 600-700 nm spectral region. A correlation was observed between emission band intensity at 635 nm and progression of renal disease.

Polarization and Spectral Shift of Benzophenone in Supercritical Water

Monte Carlo simulation and quantum mechanics calculations based on the INDO/CIS and TD-DFT methods were utilized to study the solvatochromic shift of benzophenone when changing the environment from normal water to supercritical (P = 340.2 atm and T = 673 K) condition. Solute polarization increases the dipole moment of benzophenone, compared to gas phase, by 88 and 35% in normal and supercritical conditions, giving the in-solvent dipole value of 5.8 and 4.2 D, respectively. The average number of solute-solvent hydrogen bonds was analyzed, and a large decrease of 2.3 in normal water to only 0.8 in the supercritical environment was found. By using these polarized models of benzophenone in the two different conditions of water, we performed MC simulations to generate statistically uncorrelated configurations of the solute surrounded by the solvent molecules and subsequent quantum mechanics calculations on these configurations. When changing from normal to supercritical water environment, INDO/CIS calculations explicitly considering all valence electrons of the 235 solvent water molecules resulted in a solvatochromic shift of 1425 cm(-1) for the most intense transition of benzophenone, that is, slightly underestimated in comparison with the experimentally inferred result of 1700 cm(-1). TD-B3LYP/6-311+G(2d,p) calculations on the same configurations but with benzophenone electrostatically embedded in the 320 water molecules resulted in a solvatochromic shift of 1715 cm(-1) for this transition, in very good agreement with the experimental result. When using the unpolarized model of the benzophenone, this calculated solvatochromic shift was only 640 cm(-1). Additional calculations were also made by using BHandHLYP/6-311+G(2d,p) to analyze the effect of the asymptotic decay of the exchange functional. This study indicates that, contrary to the general expectation, there is a sizable solute polarization even in the low-density regime of supercritical condition and that the inclusion of this polarization is important for a reliable description of the spectral shifts considered here.

Second-order negative-curvature methods for box-constrained and general constrained optimization

A Nonlinear Programming algorithm that converges to second-order stationary points is introduced in this paper. The main tool is a second-order negative-curvature method for box-constrained minimization of a certain class of functions that do not possess continuous second derivatives. This method is used to define an Augmented Lagrangian algorithm of PHR (Powell-Hestenes-Rockafellar) type. Convergence proofs under weak constraint qualifications are given. Numerical examples showing that the new method converges to second-order stationary points in situations in which first-order methods fail are exhibited.

Computational aspects of harmonic wavelet Galerkin methods and an application to a precipitation front propagation model

This article is dedicated to harmonic wavelet Galerkin methods for the solution of partial differential equations. Several variants of the method are proposed and analyzed, using the Burgers equation as a test model. The computational complexity can be reduced when the localization properties of the wavelets and restricted interactions between different scales are exploited. The resulting variants of the method have computational complexities ranging from O(N(3)) to O(N) (N being the space dimension) per time step. A pseudo-spectral wavelet scheme is also described and compared to the methods based on connection coefficients. The harmonic wavelet Galerkin scheme is applied to a nonlinear model for the propagation of precipitation fronts, with the front locations being exposed in the sizes of the localized wavelet coefficients. (C) 2011 Elsevier Ltd. All rights reserved.

Retrievals of a size parameter for phytoplankton and spectral light absorption by colored detrital matter from water-leaving radiances at SeaWiFS channels in a continental shelf region off Brazil

Many efforts are currently oriented toward extracting more information from ocean color than the chlorophyll a concentration. Among biological parameters potentially accessible from space, estimates of phytoplankton cell size and light absorption by colored detrital matter (CDM) would lead to an indirect assessment of major components of the organic carbon pool in the ocean, which would benefit oceanic carbon budget models. We present here 2 procedures to retrieve simultaneously from ocean color measurements in a limited number of bands, magnitudes, and spectral shapes for both light absorption by CDM and phytoplankton, along with a size parameter for phytoplankton. The performance of the 2 procedures was evaluated using different data sets that correspond to increasing uncertainties: ( 1) measured absorption coefficients of phytoplankton, particulate detritus, and colored dissolved organic matter ( CDOM) and measured chlorophyll a concentrations and ( 2) SeaWiFS upwelling radiance measurements and chlorophyll a concentrations estimated from global algorithms. In situ data were acquired during 3 cruises, differing by their relative proportions in CDM and phytoplankton, over a continental shelf off Brazil. No local information was introduced in either procedure, to make them more generally applicable. Over the study area, the absorption coefficient of CDM at 443 nm was retrieved from SeaWiFS radiances with a relative root mean square error (RMSE) of 33%, and phytoplankton light absorption coefficients in SeaWiFS bands ( from 412 to 510 nm) were retrieved with RMSEs between 28% and 33%. These results are comparable to or better than those obtained by 3 published models. In addition, a size parameter of phytoplankton and the spectral slope of CDM absorption were retrieved with RMSEs of 17% and 22%, respectively. If these methods are applied at a regional scale, the performances could be substantially improved by locally tuning some empirical relationships.

A Bayesian Analysis of Spectral ARMA Model

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Alteration of soil properties through a weathering sequence as evaluated by spectral reflectance

In the study of physical, chemical, and mineralogical data related to the weathering of soils and the quantification of their properties, remote sensing constitutes an important technique that, in addition to conventional analyses, can contribute to soil survey. The objectives of this research were to characterize and differentiate soils developed from basaltic rocks that occur in the Parana state, Brazil and to quantify soil properties based on their spectral reflectance. These observations were used to verify the relationship between the soils and reflectance with regard to weathering, organic matter (OM), and forms of Fe. From the least to the most weathered soil, we used a Typic Argiudoll (Reddish Brunizem), Rhodudalf (Terra Roxa Estruturada), and Rhodic Hapludox (Very Dark Red Latosol). The spectral reflectances between 400 and 2500 nm were obtained in the laboratory from soil samples collected at two depth increments, 0- to 20- and 40- to 60-cm, using an Infra Red Intelligent Spectroradiometer (IRIS). Correlation, regression, and discriminant estimates were used in analyzing the soil and spectral data. Results of this study indicated that soils could be separated at the soil-type level based on reflectance intensity in various absorption bands. Soil collected in the 40- to 60-cm depth appeared to have higher reflectance intensities than those from the 0- to 20-cm depth. Removal of OM from soil samples promoted higher reflectance intensity in the entire spectrum. Amorphous and crystalline Fe influenced reflectance differently. Weathering of basaltic soils was correlated with alterations in the reflectance intensities and absorption features of the spectral curves. Multivariate analysis demonstrated that this technique was efficient in the estimation of clay, silt, kaolinite, crystalline Fe, amorphous Fe, and Mg through the use of reflected energy of the soils.

Doppler ultrasonography principles and methods of evaluation of the reproductive tract in mares

Background: Doppler ultrasonography is a non-invasive real time pulse-wave technique recently used for the transrectal study of the reproductive system hemodynamics in large animals. This technic is based in the Doppler Effect Principle that proposes the change in frequency of a wave for an observer (red blood cells) moving relative to the source of the respective wave (ultrasonic transducer). This method had showed to be effective and useful for the evaluation of the in vivo equine reproductive tract increasing the diagnostic, monitoring, and predictive capabilities of theriogenology in mares. However, an accurate and truthful ultrasonic exam requires the previous knowledge of the Doppler ultrasonography principles. Review: In recent years, the capabilities of ultrasound flow imaging have increased enormously. The current Doppler ultrasound machines offer three methods of evaluation that may be used simultaneously (triplex mode). In B-mode ultrasound, a linear array of transducers simultaneously scans a plane through the tissue that can be viewed as a two-dimensional gray-scale image on screen. This mode is primarily used to identify anatomically a structure for its posterior evaluation using colored ultrasound modes (Color or Spectral modes). Colored ultrasound images of flow, whether Color or Spectral modes, are essentially obtained from measurements of moving red cells. In Color mode, velocity information is presented as a color coded overlay on top of a B-mode image, while Pulsed Wave Doppler provides a measure of the changing velocity throughout the cardiac cycle and the distribution of velocities in the sample volume represented by a spectral graphic. Color images conception varies according to the Doppler Frequency that is the difference between the frequency of received echoes by moving blood red cells and wave frequency transmitted by the transducer. To produce an adequate spectral graphic it is important determine the position and size of the simple gate. Furthermore, blood flow velocity measurement is influence by the intersection angle between ultrasonic pulses and the direction of moving blood-red cells (Doppler angle). Objectively colored ultrasound exam may be done on large arteries of the reproductive tract, as uterine and ovary arteries, or directly on the target tissue (follicle, for example). Mesovarium and mesometrium attachment arteries also can be used for spectral evaluation of the equine reproductive system. Subjectively analysis of the ovarian and uterine vascular perfusion must be done directly on the corpus luteum, follicular wall and uterus (endometrium and myometrium associated), respectively. Power-flow imaging has greater sensitivity to weak blood flow and independent of the Doppler angle, improving the evaluation of vessels with small diameters and slow blood flow. Conclusion: Doppler ultrasonography principles, methods of evaluation and reproductive system anatomy have been described. This knowledge is essential for the competent equipment acquisition and precise collection and analysis of colored ultrasound images. Otherwise, the reporting of inconsistent and not reproducible findings may result in the discredit of Doppler technology ahead of the scientific veterinary community.

Effects of catecholamines on volemic replacement with saline solution and the impact on heart rate variability in rabbits subjected to hemorrhage. A study by spectral analysis

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Analysis of NDVI time series using cross-correlation and forecasting methods for monitoring sugarcane fields in Brazil

Brazil is the largest sugarcane producer in the world and has a privileged position to attend to national and international market places. To maintain the high production of sugarcane, it is fundamental to improve the forecasting models of crop seasons through the use of alternative technologies, such as remote sensing. Thus, the main purpose of this article is to assess the results of two different statistical forecasting methods applied to an agroclimatic index (the water requirement satisfaction index; WRSI) and the sugarcane spectral response (normalized difference vegetation index; NDVI) registered on National Oceanic and Atmospheric Administration Advanced Very High Resolution Radiometer (NOAA-AVHRR) satellite images. We also evaluated the cross-correlation between these two indexes. According to the results obtained, there are meaningful correlations between NDVI and WRSI with time lags. Additionally, the adjusted model for NDVI presented more accurate results than the forecasting models for WRSI. Finally, the analyses indicate that NDVI is more predictable due to its seasonality and the WRSI values are more variable making it difficult to forecast.

Nerve fiber layer in glaucomatous hemifield loss: a case-control study with time- and spectral-domain optical coherence tomography

Purpose: To evaluate the retinal nerve fiber layer measurements with time-domain (TD) and spectral-domain (SD) optical coherence tomography (OCT), and to test the diagnostic ability of both technologies in glaucomatous patients with asymmetric visual hemifield loss. Methods: 36 patients with primary open-angle glaucoma with visual field loss in one hemifield (affected) and absent loss in the other (non-affected), and 36 age-matched healthy controls had the study eye imaged with Stratus-OCT (Carl Zeiss Meditec Inc., Dublin, California, USA) and 3 D OCT-1000 (Topcon, Tokyo, Japan). Peripapillary retinal nerve fiber layer measurements and normative classification were recorded. Total deviation values were averaged in each hemifield (hemifield mean deviation) for each subject. Visual field and retinal nerve fiber layer "asymmetry indexes" were calculated as the ratio between affected versus non-affected hemifields and corresponding hemiretinas. Results: Retinal nerve fiber layer measurements in non-affected hemifields (mean [SD] 87.0 [17.1] mu m and 84.3 [20.2] mu m, for TD and SD-OCT, respectively) were thinner than in controls (119.0 [12.2] mu m and 117.0 [17.7] mu m, P<0.001). The optical coherence tomography normative database classified 42% and 67% of hemiretinas corresponding to non-affected hemifields as abnormal in TD and SD-OCT, respectively (P=0.01). Retinal nerve fiber layer measurements were consistently thicker with TD compared to SD-OCT. Retinal nerve fiber layer thickness asymmetry index was similar in TD (0.76 [0.17]) and SD-OCT (0.79 [0.12]) and significantly greater than the visual field asymmetry index (0.36 [0.20], P<0.001). Conclusions: Normal hemifields of glaucoma patients had thinner retinal nerve fiber layer than healthy eyes, as measured by TD and SD-OCT. Retinal nerve fiber layer measurements were thicker with TD than SD-OCT. SD-OCT detected abnormal retinal nerve fiber layer thickness more often than TD-OCT.

Optic Disc Margin Anatomy in Patients with Glaucoma and Normal Controls with Spectral Domain Optical Coherence Tomography

Objective: To characterize optic nerve head (ONH) anatomy related to the clinical optic disc margin with spectral domain-optical coherence tomography (SD-OCT). Design: Cross-sectional study. Participants: Patients with open-angle glaucoma with focal, diffuse, and sclerotic optic disc damage, and age-matched normal controls. Methods: High-resolution radial SD-OCT B-scans centered on the ONH were analyzed at each clock hour. For each scan, the border tissue of Elschnig was classified for obliqueness (internally oblique, externally oblique, or nonoblique) and the presence of Bruch's membrane overhanging the border tissue. Optic disc stereophotographs were co-localized to SD-OCT data with customized software. The frequency with which the disc margin identified in stereophotographs coincided with (1) Bruch's membrane opening (BMO), defined as the innermost edge of Bruch's membrane; (2) Bruch's membrane/border tissue, defined as any aspect of either outside BMO or border tissue; or (3) border tissue, defined as any aspect of border tissue alone, in the B-scans was computed at each clock hour. Main Outcome Measures: The SD-OCT structures coinciding with the disc margin in stereophotographs. Results: There were 30 patients (10 with each type of disc damage) and 10 controls, with a median (range) age of 68.1 (42-86) years and 63.5 (42-77) years, respectively. Although 28 patients (93%) had 2 or more border tissue configurations, the most predominant one was internally oblique, primarily superiorly and nasally, frequently with Bruch's membrane overhang. Externally oblique border tissue was less frequent, observed mostly inferiorly and temporally. In controls, there was predominantly internally oblique configuration around the disc. Although the configurations were not statistically different between patients and controls, they were among the 3 glaucoma groups. At most locations, the SD-OCT structure most frequently identified as the disc margin was some aspect of Bruch's membrane and border tissue external to BMO. Bruch's membrane overhang was regionally present in the majority of patients with glaucoma and controls; however, in most cases it was not visible as the disc margin. Conclusions: The clinically perceived disc margin is most likely not the innermost edge of Bruch's membrane detected by SD-OCT. These findings have important implications for the automated detection of the disc margin and estimates of the neuroretinal rim. Financial Disclosure(s): Proprietary or commercial disclosure may be found after the references. Ophthalmology 2012;119:738-747 (C) 2012 by the American Academy of Ophthalmology.