Sensitivity analysis of a cyanobacterial growth and movement model under two different flow regimes

Bloom-forming and toxin-producing cyanobacteria remain a persistent nuisance across the world. Modelling cyanobacterial behaviour in freshwaters is an important tool for understanding their population dynamics and predicting the location and timing of the bloom events in lakes, reservoirs and rivers. A new deterministic–mathematical model was developed, which simulates the growth and movement of cyanobacterial blooms in river systems. The model focuses on the mathematical description of the bloom formation, vertical migration and lateral transport of colonies within river environments by taking into account the major factors that affect the cyanobacterial bloom formation in rivers including light, nutrients and temperature. A parameter sensitivity analysis using a one-at-a-time approach was carried out. There were two objectives of the sensitivity analysis presented in this paper: to identify the key parameters controlling the growth and movement patterns of cyanobacteria and to provide a means for model validation. The result of the analysis suggested that maximum growth rate and day length period were the most significant parameters in determining the population growth and colony depth, respectively.

Periodic time series modeling of environmental proxy records with guaranteed positive growth rate estimation

Identifying a periodic time-series model from environmental records, without imposing the positivity of the growth rate, does not necessarily respect the time order of the data observations. Consequently, subsequent observations, sampled in the environmental archive, can be inversed on the time axis, resulting in a non-physical signal model. In this paper an optimization technique with linear constraints on the signal model parameters is proposed that prevents time inversions. The activation conditions for this constrained optimization are based upon the physical constraint of the growth rate, namely, that it cannot take values smaller than zero. The actual constraints are defined for polynomials and first-order splines as basis functions for the nonlinear contribution in the distance-time relationship. The method is compared with an existing method that eliminates the time inversions, and its noise sensitivity is tested by means of Monte Carlo simulations. Finally, the usefulness of the method is demonstrated on the measurements of the vessel density, in a mangrove tree, Rhizophora mucronata, and the measurement of Mg/Ca ratios, in a bivalve, Mytilus trossulus.

Stratospheric water vapor and climate: sensitivity to the representation in radiation codes

There has been considerable interest in the climate impact of trends in stratospheric water vapor (SWV). However, the representation of the radiative properties of water vapor under stratospheric conditions remains poorly constrained across different radiation codes. This study examines the sensitivity of a detailed line-by-line (LBL) code, a Malkmus narrow-band model and two broadband GCM radiation codes to a uniform perturbation in SWV in the longwave spectral region. The choice of sampling rate in wave number space (Δν) in the LBL code is shown to be important for calculations of the instantaneous change in heating rate (ΔQ) and the instantaneous longwave radiative forcing (ΔFtrop). ΔQ varies by up to 50% for values of Δν spanning 5 orders of magnitude, and ΔFtrop varies by up to 10%. In the three less detailed codes, ΔQ differs by up to 45% at 100 hPa and 50% at 1 hPa compared to a LBL calculation. This causes differences of up to 70% in the equilibrium fixed dynamical heating temperature change due to the SWV perturbation. The stratosphere-adjusted radiative forcing differs by up to 96% across the less detailed codes. The results highlight an important source of uncertainty in quantifying and modeling the links between SWV trends and climate.

Sensitivity analysis of an ocean carbon cycle model in the North Atlantic: an investigation of parameters affecting the air-sea CO2 flux, primary production and export of detritus

The sensitivity of the biological parameters in a nutrient-phytoplankton-zooplankton-detritus (NPZD) model in the calculation of the air-sea CO2 flux, primary production and detrital export is analysed. We explore the effect on these outputs of variation in the values of the twenty parameters that control ocean ecosystem growth in a 1-D formulation of the UK Met Office HadOCC NPZD model used in GCMs. We use and compare the results from one-at-a-time and all-at-a-time perturbations performed at three sites in the EuroSITES European Ocean Observatory Network: the Central Irminger Sea (60° N 40° W), the Porcupine Abyssal Plain (49° N 16° W) and the European Station for Time series in the Ocean Canary Islands (29° N 15° W). Reasonable changes to the values of key parameters are shown to have a large effect on the calculation of the air-sea CO2 flux, primary production, and export of biological detritus to the deep ocean. Changes in the values of key parameters have a greater effect in more productive regions than in less productive areas. The most sensitive parameters are generally found to be those controlling well-established ocean ecosystem parameterisations widely used in many NPZD-type models. The air-sea CO2 flux is most influenced by variation in the parameters that control phytoplankton growth, detrital sinking and carbonate production by phytoplankton (the rain ratio). Primary production is most sensitive to the parameters that define the shape of the photosynthesis-irradiance curve. Export production is most sensitive to the parameters that control the rate of detrital sinking and the remineralisation of detritus.

Determination of a lower bound on Earth’s climate sensitivity

Transient and equilibrium sensitivity of Earth's climate has been calculated using global temperature, forcing and heating rate data for the period 1970–2010. We have assumed increased long-wave radiative forcing in the period due to the increase of the long-lived greenhouse gases. By assuming the change in aerosol forcing in the period to be zero, we calculate what we consider to be lower bounds to these sensitivities, as the magnitude of the negative aerosol forcing is unlikely to have diminished in this period. The radiation imbalance necessary to calculate equilibrium sensitivity is estimated from the rate of ocean heat accumulation as 0.37±0.03W m^−2 (all uncertainty estimates are 1−σ). With these data, we obtain best estimates for transient climate sensitivity 0.39±0.07K (W m^−2)^−1 and equilibrium climate sensitivity 0.54±0.14K (W m^−2)^−1, equivalent to 1.5±0.3 and 2.0±0.5K (3.7W m^−2)^−1, respectively. The latter quantity is equal to the lower bound of the ‘likely’ range for this quantity given by the 2007 IPCC Assessment Report. The uncertainty attached to the lower-bound equilibrium sensitivity permits us to state, within the assumptions of this analysis, that the equilibrium sensitivity is greater than 0.31K (W m^−2)^−1, equivalent to 1.16K(3.7W m^−2)^−1, at the 95% confidence level.

Intra-crystalline protein diagenesis (IcPD) in Patella vulgata. Part II: Breakdown and temperature sensitivity

Artificial diagenesis of the intra-crystalline proteins isolated from Patella vulgata was induced by isothermal heating at 140 °C, 110 °C and 80 °C. Protein breakdown was quantified for multiple amino acids, measuring the extent of peptide bond hydrolysis, amino acid racemisation and decomposition. The patterns of diagenesis are complex; therefore the kinetic parameters of the main reactions were estimated by two different methods: 1) a well-established approach based on fitting mathematical expressions to the experimental data, e.g. first-order rate equations for hydrolysis and power-transformed first-order rate equations for racemisation; and 2) an alternative model-free approach, which was developed by estimating a “scaling” factor for the independent variable (time) which produces the best alignment of the experimental data. This method allows the calculation of the relative reaction rates for the different temperatures of isothermal heating. High-temperature data were compared with the extent of degradation detected in sub-fossil Patella specimens of known age, and we evaluated the ability of kinetic experiments to mimic diagenesis at burial temperature. The results highlighted a difference between patterns of degradation at low and high temperature and therefore we recommend caution for the extrapolation of protein breakdown rates to low burial temperatures for geochronological purposes when relying solely on kinetic data.

Feedbacks and mechanisms affecting the global sensitivity of glaciers to climate change

Mass loss by glaciers has been an important contributor to sea level rise in the past, and is projected to contribute a substantial fraction of total sea level rise during the 21st century. Here, we use a model of the world's glaciers to quantify equilibrium sensitivities of global glacier mass to climate change, and to investigate the role of changes in glacier hypsometry for long-term mass changes. We find that 21st century glacier-mass loss is largely governed by the glacier's response to 20th century climate change. This limits the influence of 21st century climate change on glacier-mass loss, and explains why there are relatively small differences in glacier-mass loss under greatly different scenarios of climate change. The projected future changes in both temperature and precipitation experienced by glaciers are amplified relative to the global average. The projected increase in precipitation partly compensates for the mass loss caused by warming, but this compensation is negligible at higher temperature anomalies since an increasing fraction of precipitation at the glacier sites is liquid. Loss of low-lying glacier area, and more importantly, eventual complete disappearance of glaciers, strongly limit the projected sea level contribution from glaciers in coming centuries. The adjustment of glacier hypsometry to changes in the forcing strongly reduces the rates of global glacier-mass loss caused by changes in global mean temperature compared to rates of mass loss when hypsometric changes are neglected. This result is a second reason for the relatively weak dependence of glacier-mass loss on future climate scenario, and helps explain why glacier-mass loss in the first half of the 20th century was of the same order of magnitude as in the second half of the 20th century, even though the rate of warming was considerably smaller.

Sensitivity of OMI SO2 measurements to variable eruptive behaviour at Soufrière Hills Volcano, Montserrat

Since 2004, the satellite-borne Ozone Mapping Instrument (OMI) has observed sulphur dioxide (SO2) plumes during both quiescence and effusive eruptive activity at Soufrière Hills Volcano, Montserrat. On average, OMI detected a SO2 plume 4-6 times more frequently during effusive periods than during quiescence in the 2008-2010 period. The increased ability of OMI to detect SO2 during eruptive periods is mainly due to an increase in plume altitude rather than a higher SO2 emission rate. Three styles of eruptive activity cause thermal lofting of gases (Vulcanian explosions; pyroclastic flows; a hot lava dome) and the resultant plume altitudes are estimated from observations and models. Most lofting plumes from Soufrière Hills are derived from hot domes and pyroclastic flows. Although Vulcanian explosions produced the largest plumes, some produced only negligible SO2 signals detected by OMI. OMI is most valuable for monitoring purposes at this volcano during periods of lava dome growth and during explosive activity.

Effect of temperature and strain rate on the compressive behaviour of supramolecular polyurethane

Supramolecular polyurethanes (SPUs) possess thermoresponsive and thermoreversible properties, and those characteristics are highly desirable in both bulk commodity and value-added applications such as adhesives, shape-memory materials, healable coatings and lightweight, impact-resistant structures (e.g. protection for mobile electronics). A better understanding of the mechanical properties, especially the rate and temperature sensitivity, of these materials are required to assess their suitability for different applications. In this paper, a newly developed SPU with tuneable thermal properties was studied, and the response of this SPU to compressive loading over strain rates from 10−3 to 104 s−1 was presented. Furthermore, the effect of temperature on the mechanical response was also demonstrated. The sample was tested using an Instron mechanical testing machine for quasi-static loading, a home-made hydraulic system for moderate rates and a traditional split Hopkinson pressure bars (SHPBs) for high strain rates. Results showed that the compression stress-strain behaviour was affected significantly by the thermoresponsive nature of SPU, but that, as expected for polymeric materials, the general trends of the temperature and the rate dependence mirror each other. However, this behaviour is more complicated than observed for many other polymeric materials, as a result of the richer range of transitions that influence the behaviour over the range of temperatures and strain rates tested.

Intra-strain variations of baroreflex sensitivity in young Wistar-Kyoto rats

Purpose: To compare baroreflex sensitivity among conscious rats of the same strain. Methods: Male WKY rats (eight weeks old) were studied. Cannulas were inserted into the abdominal aortic artery through the right femoral artery to measure mean arterial pressure (MAP) and heart rate (HR). Baroreflex gain was calculated as the ratio between variation of HR in function of the MAP variation (Delta HR/Delta MAP) tested with a depressor dose of sodium nitroprusside (SNP, 50 mu g/kg, iv) and with a pressor dose of phenylephrine (PE, 8 mu g/kg, iv). We divided the rats into four groups: 1) Low bradycardic baroreflex (LB), BG between -1 and -2 bpm/mmHg tested with PE; 2) High bradycardic baroreflex (HB), BG < -2 bpm/mmHg tested with PE; 3) Low tachycardic baroreflex (LT), BG between -1 and -2 bpm/mmHg tested with SNP and; 4) High tachycardic baroreflex (HT), BG < -2 bpm/mmHg tested with SNP. Significant differences were considered for p<0.05. Results: Approximately 82% of the rats presented reduced bradycardic reflex while 22 showed attenuated tachycardic reflex. No alterations were noted regarding basal MAP and HR, tachycardic and bradycardic peak and HR range. Conclusions: There was alteration in baroreflex sensitivity among rats of the same strain. Care should be taken when interpreting studies employing WKY as a control for the SHR.

The variability of baroreflex sensitivity in juvenile, spontaneously hypertensive rats

In this study the baroreflex sensitivity of conscious, juvenile, spontaneously hypertensive rats (SHRs) was compared. The study population consisted of 19 eight-week-old male SHRs. The baroreflex sensitivity was quantified as the derivative of the variation in heart rate (HR) and the variation of mean arterial pressure (baroreflex sensitivity = Delta HR/Delta MAP). MAP was manipulated with sodium nitroprusside (SNP) and phenylephrine (PHE), administered via an inserted cannula in the right femoral vein. The SHRs were divided into four groups: (1) low bradycardic baroreflex (LB) where the baroreflex gain (BG) was between 0 and 1 bpm/mmHg with PHE; (2) high bradycardic baroreflex (HB), where the BG was < -1 bpm/mmHg with PHE; (3) low tachycardic baroreflex (LT) where the BC was between 0 and 3 bpm/mmHg with SNP; (4) high tachycardic baroreflex (HT) where the BG was > 3 bpm/mmHg with SNP. We noted that 36.8% of the rats presented with an increased bradycardic reflex, while 27.8% demonstrated an attenuated tachycardic reflex. No significant alterations were noted regarding the basal MAP and HR. There were significant differences in the baroreflex sensitivity between SHRs in the same laboratory. One should be careful when interpreting studies employing the SHR as a research model.

Baroreflex sensitivity and oxidative stress in the LDL receptor knockout mice

This study alms at observing the effect of low-density lipoprotein (LDL) receptor deficiency in cholesterol blood levels, baroreflex sensitivity (BRS), nitric oxide (NO) bioavailability, and oxidative stress. The lack of LDL receptors in mice significantly increased the cholesterol blood levels (179 +/- 35 vs. 109 +/- 13 mg/dL) in the knockout (KO) mice compared to control. There was no difference in basal mean arterial pressure and heart rate between the groups. However, in KO mice the BRS was significantly attenuated and the antioxidant enzyme activities, measured in erythrocytes and heart, were significantly decreased. On the other hand, the oxidative damage measured by chemiluminescence and carbonyls was increased, while total plasma nitrate levels were lower in KO mice, indicating a decrease in NO availability. In conclusion, these results indicate that the lack of LDL receptor increased cholesterol blood levels, induced oxidative stress and decreased BRS. (C) 2008 Elsevier GmbH. All rights reserved.

Inflammation, oxidative stress, glomerular filtration rate, and albuminuria in elderly men : a cross-sectional study

BACKGROUND: The role of inflammation and oxidative stress in mild renal impairment in the elderly is not well studied. Accordingly, we aimed at investigating the associations between estimated glomerular filtration rate (eGFR), albumin/creatinine ratio (ACR), and markers of different inflammatory pathways and oxidative stress in a community based cohort of elderly men. FINDINGS: Cystatin C-based GFR, ACR, and biomarkers of cytokine-mediated inflammation (interleukin-6, high-sensitivity C-reactive protein[CRP], serum amyloid A[SAA]), cyclooxygenase-mediated inflammation (urinary prostaglandin F2alpha [PGF2alpha]), and oxidative stress (urinary F2 isoprostanes) were assessed in the Uppsala Longitudinal Study of Adult Men(n = 647, mean age 77 years). RESULTS: In linear regression models adjusting for age, BMI, smoking, blood pressure, LDL-cholesterol, HDL-cholesterol, triglycerides, and treatment with statins, ACE-inhibitors, ASA, and anti-inflammatory agents, eGFR was inversely associated with CRP, interleukin-6, and SAA (beta-coefficient -0.13 to -0.19, p < 0.001 for all), and positively associated with urinary F2-isoprostanes (beta-coefficient 0.09, p = 0.02). In line with this, ACR was positively associated with CRP, interleukin-6, and SAA (beta- coefficient 0.09-0.12, p < 0.02 for all), and negatively associated with urinary F2-isoprostanes (beta-coefficient -0.12, p = 0.002). The associations were similar but with lower regression coefficients in a sub-sample with normal eGFR (>60 ml/min/1.73 m2, n = 514), with the exception that F2-isoprostane and SAA were no longer associated with eGFR. CONCLUSION: Our data indicate that cytokine-mediated inflammation is involved in the early stages of impaired kidney function in the elderly, but that cyclooxygenase-mediated inflammation does not play a role at this stage. The unexpected association between higher eGFR/lower albuminuria and increased F2-isoprostanes in urine merits further studies.

Cash flow and discount rate risk decomposition and ICAPM for the US and brazilian stock markets

Esse trabalho é uma aplicação do modelo intertemporal de apreçamento de ativos desenvolvido por Campbell (1993) e Campbell e Vuolteenaho (2004) para as carteiras de Fama-French 2x3 brasileiras no period de janeiro de 2003 a abril de 2012 e para as carteiras de Fama-French 5x5 americanas em diferentes períodos. As varíaveis sugeridas por Campbell e Vuolteenaho (2004) para prever os excessos de retorno do mercado acionário americano no period de 1929 a 2001 mostraram-se também bons preditores de excesso de retorno para o mercado brasileiro no período recente, com exceção da inclinação da estrutura a termo das taxas de juros. Entretanto, mostramos que um aumento no small stock value spread indica maior excesso de retorno no futuro, comportamento que não é coerente com a explicação para o prêmio de valor sugerida pelo modelo intertemporal. Ainda, utilizando os resíduos do VAR preditivo para definir o risco de choques de fluxo de caixa e de choques nas taxas de desconto das carteiras de teste, verificamos que o modelo intertemporal resultante não explica adequadamente os retornos observados. Para o mercado norte-americano, concluímos que a abilidade das variáveis propostas para explicar os excessos de retorno do mercado varia no tempo. O sucesso de Campbell e Vuolteenaho (2004) em explicar o prêmio de valor para o mercado norte-americano na amostra de 1963 a 2001 é resultado da especificação do VAR na amostra completa, pois mostramos que nenhuma das varíaveis é um preditor de retorno estatisticamente significante nessa sub-amostra.

An improved and fast UHPLC-PDA methodology for determination of L-ascorbic and dehydroascorbic acids in fruits and vegetables: evaluation of degradation rate during storage

This study provides a versatile validated method to determine the total vitamin C content, as the sum of the contents of L-ascorbic acid (L-AA) and dehydroascorbic acid (DHAA), in several fruits and vegetables and its degradability with storage time. Seven horticultural crops from two different origins were analyzed using an ultrahigh-performance liquid chromatographic–photodiode array (UHPLC-PDA) system, equipped with a new trifunctional high strength silica (100% silica particle) analytical column (100 mm×2.1 mm, 1.7 μm particle size) using 0.1% (v/v) formic acid as mobile phase, in isocratic mode. This new stationary phase, specially designed for polar compounds, overcomes the problems normally encountered in HPLC and is suitable for the analysis of large batches of samples without L-AA degradation. In addition, it proves to be an excellent alternative to conventional C18 columns for the determination of L-AA in fruits and vegetables. The method was fully validated in terms of linearity, detection (LOD) and quantification (LOQ) limits, accuracy, and inter/intraday precision. Validation experiments revealed very good recovery rate of 96.6±4.4% for L-AA and 103.1±4.8 % for total vitamin C, good linearity with r2-values >0.999 within the established concentration range, excellent repeatability (0.5%), and reproducibility (1.6%) values. The LOD of the method was 22 ng/mL whereas the LOQ was 67 ng/mL. It was possible to demonstrate that L-AA and DHAA concentrations in the different horticulture products varied oppositely with time of storage not always affecting the total amount of vitamin C during shelf-life. Locally produced fruits have higher concentrations of vitamin C, compared with imported ones, but vegetables showed the opposite trend. Moreover, this UHPLC-PDA methodology proves to be an improved, simple, and fast approach for determining the total content of vitamin C in various food commodities, with high sensitivity, selectivity, and resolving power within 3 min of run analysis.