Twin experiments with the equivalent weights particle filter and HadCM3

This paper investigates the use of a particle filter for data assimilation with a full scale coupled ocean–atmosphere general circulation model. Synthetic twin experiments are performed to assess the performance of the equivalent weights filter in such a high-dimensional system. Artificial 2-dimensional sea surface temperature fields are used as observational data every day. Results are presented for different values of the free parameters in the method. Measures of the performance of the filter are root mean square errors, trajectories of individual variables in the model and rank histograms. Filter degeneracy is not observed and the performance of the filter is shown to depend on the ability to keep maximum spread in the ensemble.

A Variance-Minimizing Filter for Large-Scale Applications

A truly variance-minimizing filter is introduced and its per for mance is demonstrated with the Korteweg– DeV ries (KdV) equation and with a multilayer quasigeostrophic model of the ocean area around South Africa. It is recalled that Kalman-like filters are not variance minimizing for nonlinear model dynamics and that four - dimensional variational data assimilation (4DV AR)-like methods relying on per fect model dynamics have dif- ficulty with providing error estimates. The new method does not have these drawbacks. In fact, it combines advantages from both methods in that it does provide error estimates while automatically having balanced states after analysis, without extra computations. It is based on ensemble or Monte Carlo integrations to simulate the probability density of the model evolution. When obser vations are available, the so-called importance resampling algorithm is applied. From Bayes’ s theorem it follows that each ensemble member receives a new weight dependent on its ‘ ‘distance’ ’ t o the obser vations. Because the weights are strongly var ying, a resampling of the ensemble is necessar y. This resampling is done such that members with high weights are duplicated according to their weights, while low-weight members are largely ignored. In passing, it is noted that data assimilation is not an inverse problem by nature, although it can be for mulated that way . Also, it is shown that the posterior variance can be larger than the prior if the usual Gaussian framework is set aside. However , i n the examples presented here, the entropy of the probability densities is decreasing. The application to the ocean area around South Africa, gover ned by strongly nonlinear dynamics, shows that the method is working satisfactorily . The strong and weak points of the method are discussed and possible improvements are proposed.

Analysis Scheme in the Ensemble Kalman Filter

This paper discusses an important issue related to the implementation and interpretation of the analysis scheme in the ensemble Kalman filter . I t i s shown that the obser vations must be treated as random variables at the analysis steps. That is, one should add random perturbations with the correct statistics to the obser vations and generate an ensemble of obser vations that then is used in updating the ensemble of model states. T raditionally , this has not been done in previous applications of the ensemble Kalman filter and, as will be shown, this has resulted in an updated ensemble with a variance that is too low . This simple modification of the analysis scheme results in a completely consistent approach if the covariance of the ensemble of model states is interpreted as the prediction error covariance, and there are no further requirements on the ensemble Kalman filter method, except for the use of an ensemble of sufficient size. Thus, there is a unique correspondence between the error statistics from the ensemble Kalman filter and the standard Kalman filter approach

Assimilation of Geosat Altimeter Data for the Agulhas Current Using the Ensemble Kalman Filter with a Quasigeostrophic Model

The ring-shedding process in the Agulhas Current is studied using the ensemble Kalman filter to assimilate geosat altimeter data into a two-layer quasigeostrophic ocean model. The properties of the ensemble Kalman filter are further explored with focus on the analysis scheme and the use of gridded data. The Geosat data consist of 10 fields of gridded sea-surface height anomalies separated 10 days apart that are added to a climatic mean field. This corresponds to a huge number of data values, and a data reduction scheme must be applied to increase the efficiency of the analysis procedure. Further, it is illustrated how one can resolve the rank problem occurring when a too large dataset or a small ensemble is used.

Implicit equal-weights particle filter

Filter degeneracy is the main obstacle for the implementation of particle filter in non-linear high-dimensional models. A new scheme, the implicit equal-weights particle filter (IEWPF), is introduced. In this scheme samples are drawn implicitly from proposal densities with a different covariance for each particle, such that all particle weights are equal by construction. We test and explore the properties of the new scheme using a 1,000-dimensional simple linear model, and the 1,000-dimensional non-linear Lorenz96 model, and compare the performance of the scheme to a Local Ensemble Kalman Filter. The experiments show that the new scheme can easily be implemented in high-dimensional systems and is never degenerate, with good convergence properties in both systems.

Ratios of UV, PAR and NIR components to global solar radiation measured at Botucatu site in Brazil

The relationships between the four radiant fluxes are analyzed based on a 4 year data archive of hourly and daily global ultraviolet (I(UV)), photosynthetically active-PAR (I(PAR)), near infrared (I(NIR)) and broadband global solar radiation (I(G)) collected at Botucatu, Brazil. These data are used to establish both the fractions of spectral components to global solar radiation and the proposed linear regression models. Verification results indicated that the proposed regression models predict accurately the spectral radiant fluxes at least for the Brazilian environment. Finally, results obtained in this analysis agreed well with most published results in the literature. (c) 2010 Elsevier Ltd. All rights reserved.

Modeling hourly and daily fractions of UV, PAR and NIR to global solar radiation under various sky conditions at Botucatu, Brazil

In this analysis, using available hourly and daily radiometric data performed at Botucatu, Brazil, several empirical models relating ultraviolet (UV), photosynthetically active (PAR) and near infrared (NIR) solar global components with solar global radiation (G) are established. These models are developed and discussed through clearness index K(T) (ratio of the global-to-extraterrestrial solar radiation). Results obtained reveal that the proposed empirical models predict hourly and daily values accurately. Finally. the overall analysis carried Out demonstrates that the sky conditions are more important in developing correlation models between the UV component and the global solar radiation. The linear regression models derived to estimate PAR and NIR components may be obtained without sky condition considerations within a maximum variation of 8%. In the case of UV, not taking into consideration the sky condition may cause a discrepancy of up to 18% for hourly values and 15% for daily values. (C) 2008 Elsevier Ltd. All rights reserved.

Lines of heavy elements and Fe II in the UV of CS 31082-001

The abundance of heavy r-elements may provide a better understanding of the r-process, and the determination of several reference r-elements should allow a better determination of a star`s age. The space UV region (lambda < 3000 angstrom) presents a large number of lines of the heavy elements, and in the case of some elements, such as Bi, Pt, Au, detectable lines are not available elsewhere. The extreme ""r-process star"" CS 31082-001 ([Fe/H] = -2.9) was observed in the space UV to determine abundances of the heaviest stable elements, using STIS on board Hubble Space Telescope.

The Diffusion Kernel Filter

A particle filter method is presented for the discrete-time filtering problem with nonlinear ItA ` stochastic ordinary differential equations (SODE) with additive noise supposed to be analytically integrable as a function of the underlying vector-Wiener process and time. The Diffusion Kernel Filter is arrived at by a parametrization of small noise-driven state fluctuations within branches of prediction and a local use of this parametrization in the Bootstrap Filter. The method applies for small noise and short prediction steps. With explicit numerical integrators, the operations count in the Diffusion Kernel Filter is shown to be smaller than in the Bootstrap Filter whenever the initial state for the prediction step has sufficiently few moments. The established parametrization is a dual-formula for the analysis of sensitivity to gaussian-initial perturbations and the analysis of sensitivity to noise-perturbations, in deterministic models, showing in particular how the stability of a deterministic dynamics is modeled by noise on short times and how the diffusion matrix of an SODE should be modeled (i.e. defined) for a gaussian-initial deterministic problem to be cast into an SODE problem. From it, a novel definition of prediction may be proposed that coincides with the deterministic path within the branch of prediction whose information entropy at the end of the prediction step is closest to the average information entropy over all branches. Tests are made with the Lorenz-63 equations, showing good results both for the filter and the definition of prediction.

Effects of N supply on the accumulation of photosynthetic pigments and photoprotectors in Gracilaria tenuistipitata (Rhodophyta) cultured under UV radiation

We have studied the effects of nitrate supply under photosynthetic active radiation (PAR) plus ultraviolet radiation (UVR) exposure on photosynthetic pigments (chlorophyll a and carotenoids), photoprotective UV screen mycosporine-like amino acids (MAAs), and photosynthetic parameters, including the maximum quantum yield (F(v)/F(m)) and electron transport rate (ETR) on the red agarophyte Gracilaria tenuistipitata. Apical tips of G. tenuistipitata were cultivated under ten different concentrations of NO(3)(-) for 7 days. It has been shown that G. tenuistipitata cultured under laboratory conditions has the ability to accumulate high amounts of MAAs following a nitrate concentration-dependent manner under PAR+UVR. Two MAAs were identified, shinorine and porphyra-334. The relative concentration of the first increased under high concentrations of nitrate, while the second one decreased. The presence of antheraxanthin is reported for the first time in this macro-algae, which also contains zeaxanthin, lutein, and beta-carotene. The accumulation of pigments, photoprotective compounds, and photosynthetic parameters of G. tenuistipitata is directly related to N availability. All variables decreased under low N supplies and reached constant maximum values with supplements higher than 0.5 mM NO(3)(-). Our results suggest a high potential to acclimation and photoprotection against stress factors (including high PAR and UVR) directly related to N availability for G. tenuistipitata.

The effects of UV radiation on photosynthesis estimated as chlorophyll fluorescence in Zygnemopsis decussata (Chlorophyta) growing in a high mountain lake (Sierra Nevada, Southern Spain)

The effect of increased UV radiation on photosynthesis estimated as in vivo chlorophyll fluorescence i.e. optimal quantum yield (F(v)/F(m)) and electron transport rate (ETR) in the green filamentous alga Zygnemopsis decussata (Streptophyta, Zygnematales) growing in the high mountain lake ""La Caldera"" (Sierra Nevada, Spain) at 3050 m altitude was evaluated. Two sets of in situ experiments were conducted: (1) On July 2006, F(v)/F(m) was measured throughout the day at different depths (0.1, 0.25, 0.5 and 1 m) and in the afternoon. ETR and phenolic compounds were determined. In addition, in order to analyze the effect of UV radiation, F(v)/F(m) was determined in algae incubated for 3 days at 0.5m under three different light treatments: PAR+UVA+UVB (PAB). PAR+UVA (PA) and PAR (P). (2) On August 2007, F(v)/F(m) was determined under PAB, PA and P treatments and desiccation/rehydration conditions. F(v)/F(m) decreased in algae growing in surface waters (0.1 m) but also at 1 m depth compared to that at 0.5 in depth. The decrease of F(v)/F(m) at noon due to photoinhibition was small (less than 10%) except in algae growing at 1 m depth (44%). The maximal electron transport rate was 3.5-5 times higher in algae growing at 0.25-0.5 m respectively than that at 0.1 and 1 m depth. These results are related to the accumulation of phenolic compounds: i.e. the algae at 0.25-0.5 in presentedrespectively about a 3-5 times higher concentration of phenolic compounds than that of algae at 0.1-1 m depth. The protection mechanisms seem to be stimulated by UVB radiation, since F(v)/F(m) was higher in the presence of UVB (PAB treatment) compared to PA or P treatments. UVA exerts the main photoinhibitory effect, not Only at midday, but also in the afternoon. UVB radiation also had a protective effect in algae grown under desiccation conditions for three days. During re-hydration, the rapid increase of F(v)/F(m) (after 1 h) was higher in the UVB-grown algae than in algae grown under UVA radiation. After 5 h. F(v)/F(m) values were similar in algae submitted to desiccation/rehydration under PAB and P treatments as they were in the control (submerged algae). The combined effect of desiccation and UVA produced the greatest decrease of photosynthesis in Z. decussata. Thifs UVB, in contrast to other species, may support the recovery process. Z. decussata can acclimate to severe stress, conditions in this high mountain lake by the photoprotection mechanism induced by UVB radiation through dynamic photoinhibition and the accumulation of phenolic compounds (UV screen and antioxidant substances).

CPDs and 6-4PPs play different roles in UV-induced cell death in normal and NER-deficient human cells

Ultraviolet (UV) light generates two major DNA lesions: cyclobutane pyrimidine dimers (CPDs) and pyrimidine-(6-4)-pyrimidone photoproducts (6-4PPs), but the specific participation of these two lesions in the deleterious effects of UV is a longstanding question. In order to discriminate the precise role of unrepaired CPDs and 6-4PPs in UV-induced responses triggering cell death, human fibroblasts were transduced by recombinant adenoviruses carrying the CPD-photolyase or 6-4PP-photolyase cDNAs. Both photolyases were able to prevent UV-induced apoptosis in cells deficient for nucleotide excision repair (NER) to a similar extent, while in NER-proficient cells UV-induced apoptosis was prevented only by CPD-photolyase, with no effects observed when 6-4PPs were removed by the specific photolyase. These results strongly suggest that both CPDs and 6-4PPs contribute to UV-induced apoptosis in NER-deficient cells, while in NER-proficient cells, CPDs are the only lesions responsible for UV-killing, probably due to the rapid repair of 6-4PPs by NER. As a consequence, the difference in skin photosensitivity, including carcinogenesis, of most of the xeroderma pigmentosum patients and of normal people is probably not only a quantitative aspect, but depends on the type of DNA damage induced by sunlight and its rate of repair. (c) 2007 Elsevier B.V. All rights reserved.

Defective transcription/repair factor IN recruitment to specific UV lesions in trichothiodystrophy syndrome

Most trichothiodystrophy (TTD) patients present mutations in the xeroderma pigmentosum D (XPD) gene, coding for a subunit of the transcription/repair factor IIH (TFHH) complex involved in nucleotide excision repair (NER) and transcription. After UV irradiation, most TTD/XPD patients are more severely affected in the NER of cyclobutane pyrimidine dimers (CPD) than of 6-4-photoproducts (6-4PP). The reasons for this differential DNA repair defect are unknown. Here we report the first study of NER in response to CPDs or 6-4PPs separately analyzed in primary fibroblasts. This was done by using heterologous photorepair; recombinant adenovirus vectors carrying photolyases enzymes that repair CPD or 64PP specifically by using the energy of light were introduced in different cell lines. The data presented here reveal that some mutations affect the recruitment of TFHH specifically to CPDs, but not to 6-4PPs. This deficiency is further confirmed by the inability of TTD/XPD cells to recruit, specifically for CPDs, NER factors that arrive in a TFIIH-dependent manner later in the NER pathway. For 6-4PPs, we show that TFHH complexes carrying an NH2-terminal XPD mutated protein are also deficient in recruitment of NER proteins downstream of TFUH. Treatment with the histone deacetylase inhibitor trichostatin A allows the recovery of TFHH recruitment to CPDs in the studied TTD cells and, for COOH-terminal XPD mutations, increases the repair synthesis and survival after UV, suggesting that this defect can be partially related with accessibility of DNA damage in closed chromatin regions.