Precipitation diurnal cycle and summer climatology assessment over South America: An evaluation of Regional Climate Model version 3 simulations

Regional Climate Model version 3 (RegCM3) simulations of 17 summers (1988-2004) over part of South America south of 5 degrees S were evaluated to identify model systematic errors. Model results were compared to different rainfall data sets (Climate Research Unit (CRU), Climate Prediction Center (CPC), Global Precipitation Climatology Project (GPCP), and National Centers for Environmental Prediction (NCEP) reanalysis), including the five summers mean (1998-2002) precipitation diurnal cycle observed by the Tropical Rainfall Measuring Mission (TRMM)-Precipitation Radar (PR). In spite of regional differences, the RegCM3 simulates the main observed aspects of summer climatology associated with the precipitation (northwest-southeast band of South Atlantic Convergence Zone (SACZ)) and air temperature (warmer air in the central part of the continent and colder in eastern Brazil and the Andes Mountains). At a regional scale, the main RegCM3 failures are the underestimation of the precipitation in the northern branch of the SACZ and some unrealistic intense precipitation around the Andes Mountains. However, the RegCM3 seasonal precipitation is closer to the fine-scale analyses (CPC, CRU, and TRMM-PR) than is the NCEP reanalysis, which presents an incorrect north-south orientation of SACZ and an overestimation of its intensity. The precipitation diurnal cycle observed by TRMM-PR shows pronounced contrasts between Tropics and Extratropics and land and ocean, where most of these features are simulated by RegCM3. The major similarities between the simulation and observation, especially the diurnal cycle phase, are found over the continental tropical and subtropical SACZ regions, which present afternoon maximum (1500-1800 UTC) and morning minimum (0900-1200 UTC). More specifically, over the core of SACZ, the phase and amplitude of the simulated precipitation diurnal cycle are very close to the TRMM-PR observations. Although there are amplitude differences, the RegCM3 simulates the observed nighttime rainfall in the eastern Andes Mountains, over the Atlantic Ocean, and also over northern Argentina. The main simulation deficiencies are found in the Atlantic Ocean and near the Andes Mountains. Over the Atlantic Ocean the convective scheme is not triggered; thus the rainfall arises from the grid-scale scheme and therefore differs from the TRMM-PR. Near the Andes, intense (nighttime and daytime) simulated precipitation could be a response of an incorrect circulation and topographic uplift. Finally, it is important to note that unlike most reported bias of global models, RegCM3 does not trigger the moist convection just after sunrise over the southern part of the Amazon.

A Lagrangian identification of major sources of moisture over Central Brazil and La Plata Basin

[1] This work examines the main sources of moisture over Central Brazil and La Plata Basin during the year through a new Lagrangian diagnosis method which identifies the humidity contributions to the moisture budget over a region. This methodology computes budgets of evaporation minus precipitation by calculating changes in the specific humidity along back-trajectories for the previous 10 d. The origin of all air masses residing over each region was tracked during a period of 5 years (2000-2004). These regions were selected because they coincide with two centers of action of a known dipole precipitation variability mode observed in different temporal scales (from intra seasonal up to inter decadal timescales) and are related to the climatic variability of the South American Monsoon System. The results suggested the importance of the tropical south Atlantic as a moisture source for Central Brazil, and of recycling for La Plata basin. It seems that the Tropical South Atlantic plays an important role as a moisture source for Central Brazil and La Plata basin along the year, particularly during the austral summer. The north Atlantic is also an additional source for both regions during the austral summer.

Placement of Kuhlmanniodendron Fiaschi & Groppo in Lindackerieae (Achariaceae, Malpighiales) confirmed by analyses of rbcL sequences, with notes on pollen morphology and wood anatomy

The phylogenetic placement of Kuhlmanniodendron Fiaschi & Groppo (Achariaceae) within Malpighiales was investigated with rbcL sequence data. This genus was recently created to accommodate Carpotroche apterocarpa Kuhlm., a poorly known species from the rainforests of Espirito Santo, Brazil. One rbcL sequence was obtained from Kuhlmanniodendron and analyzed with 73 additional sequences from Malpighiales, and 8 from two closer orders, Oxalidales and Celastrales, all of which were available at Genbank. Phylogenetic analyses were carried out with maximum parsimony and Bayesian inference; bootstrap analyses were used in maximum parsimony to evaluate branch support. The results confirmed the placement of Kuhlmanniodendron together with Camptostylus, Lindackeria, Xylotheca, and Caloncoba in a strongly supported clade (posterior probability = 0.99) that corresponds with the tribe Lindackerieae of Achariaceae (Malpighiales). Kuhlmanniodendron also does not appear to be closely related to Oncoba (Salicaceae), an African genus with similar floral and fruit morphology that has been traditionally placed among cyanogenic Flacourtiaceae (now Achariaceae). A picrosodic paper test was performed in herbarium dry leaves, and the presence of cyanogenic glycosides, a class of compounds usually found in Achariaceae, was detected. Pollen morphology and wood anatomy of Kuhlmanniodendron were also investigated, but both pollen (3-colporate and microreticulate) and wood, with solitary to multiple vessels, scalariform perforation plates and other features, do not seem to be useful to distinguish this genus from other members of the Achariaceae and are rather common among the eudicotyledons as a whole. However, perforated ray cells with scalariform plates, an uncommon wood character, present in Kuhlmanniodendron are similar to those found in Kiggelaria africana (Pangieae, Achariaceae), but the occurrence of such cells is not mapped among the angiosperms, and it is not clear how homoplastic this character could be.

Cytochrome-b variation in Apis mellifera samples and its association with COI-COII patterns

Five Mbo I (Mbo-A, Mbo-M, Mbo-C(1), Mbo-C(2) and Mbo-C(3)) and Hinf I (Hinf-1 to Hinf-5) patterns were observed in Apis mellifera samples after restriction of a 485 bp fragment of the mitochondrial cytochrome-b (cyt-b) gene. Associating the cyt-b Restriction fragment length polymorphism (RFLP) pattern of each sample to its respective previously established COI-COII (Dra I sites) pattern, five restriction patterns (Mbo-C(1), Mbo-C(2), Mbo-C(3), Hinf-1 and Hinf-4) were observed in samples of maternal origin associated to the evolutionary branch C. No deletions or insertions were observed and the nucleotide substitution rate was estimated at 5.4%. Higher nucleotide diversity was observed among the branch C-haplotypes when compared with A and M lineages. Further studies are needed to confirm if the cyt-b + COI-COII haplotypes help to assign certain phylogeographic patterns to the branch C and to clarify phylogenetic relationships among A. mellifera subspecies.

SGLT1 protein expression in plasma membrane of acinar cells correlates with the sympathetic outflow to salivary glands in diabetic and hypertensive rats

Salivary gland dysfunction is a feature in diabetes and hypertension. We hypothesized that sodium-glucose cotransporter 1 (SGLT1) participates in salivary dysfunctions through a sympathetic- and protein kinase A (PKA)-mediated pathway. In Wistar-Kyoto (WKY), diabetic WKY (WKY-D), spontaneously hypertensive (SHR), and diabetic SHR (SHR-D) rats, PKA/SGLT1 proteins were analyzed in parotid and submandibular glands, and the sympathetic nerve activity (SNA) to the glands was monitored. Basal SNA was threefold higher in SHR (P < 0.001 vs. WKY), and diabetes decreased this activity (similar to 50%, P < 0.05) in both WKY and SHR. The catalytic subunit of PKA and the plasma membrane SGLT1 content in acinar cells were regulated in parallel to the SNA. Electrical stimulation of the sympathetic branch to salivary glands increased (similar to 30%, P < 0.05) PKA and SGLT1 expression. Immunohistochemical analysis confirmed the observed regulations of SGLT1, revealing its location in basolateral membrane of acinar cells. Taken together, our results show highly coordinated regulation of sympathetic activity upon PKA activity and plasma membrane SGLT1 content in salivary glands. Furthermore, the present findings show that diabetic- and/or hypertensive-induced changes in the sympathetic activity correlate with changes in SGLT1 expression in basolateral membrane of acinar cells, which can participate in the salivary glands dysfunctions reported by patients with these pathologies.

A non-autonomous strongly damped wave equation: Existence and continuity of the pullback attractor

In this paper we consider the strongly damped wave equation with time-dependent terms u(tt) - Delta u - gamma(t)Delta u(t) + beta(epsilon)(t)u(t) = f(u), in a bounded domain Omega subset of R(n), under some restrictions on beta(epsilon)(t), gamma(t) and growth restrictions on the nonlinear term f. The function beta(epsilon)(t) depends on a parameter epsilon, beta(epsilon)(t) -> 0. We will prove, under suitable assumptions, local and global well-posedness (using the uniform sectorial operators theory), the existence and regularity of pullback attractors {A(epsilon)(t) : t is an element of R}, uniform bounds for these pullback attractors, characterization of these pullback attractors and their upper and lower semicontinuity at epsilon = 0. (C) 2010 Elsevier Ltd. All rights reserved.

Complexity and anisotropy in host morphology make populations less susceptible to epidemic outbreaks

One of the challenges in epidemiology is to account for the complex morphological structure of hosts such as plant roots, crop fields, farms, cells, animal habitats and social networks, when the transmission of infection occurs between contiguous hosts. Morphological complexity brings an inherent heterogeneity in populations and affects the dynamics of pathogen spread in such systems. We have analysed the influence of realistically complex host morphology on the threshold for invasion and epidemic outbreak in an SIR (susceptible-infected-recovered) epidemiological model. We show that disorder expressed in the host morphology and anisotropy reduces the probability of epidemic outbreak and thus makes the system more resistant to epidemic outbreaks. We obtain general analytical estimates for minimally safe bounds for an invasion threshold and then illustrate their validity by considering an example of host data for branching hosts (salamander retinal ganglion cells). Several spatial arrangements of hosts with different degrees of heterogeneity have been considered in order to separately analyse the role of shape complexity and anisotropy in the host population. The estimates for invasion threshold are linked to morphological characteristics of the hosts that can be used for determining the threshold for invasion in practical applications.

Empirical analysis of the Lieb-Oxford bound in ions and molecules

Universal properties of the Coulomb interaction energy apply to all many-electron systems. Bounds on the exchange-correlation energy, in particular, are important for the construction of improved density functionals. Here we investigate one such universal property-the Lieb-Oxford lower bound-for ionic and molecular systems. In recent work [J Chem Phys 127, 054106 (2007)], we observed that for atoms and electron liquids this bound may be substantially tightened. Calculations for a few ions and molecules suggested the same tendency, but were not conclusive due to the small number of systems considered. Here we extend that analysis to many different families of ions and molecules, and find that for these, too, the bound can be empirically tightened by a similar margin as for atoms and electron liquids. Tightening the Lieb-Oxford bound will have consequences for the performance of various approximate exchange-correlation functionals. (C) 2008 Wiley Periodicals Inc.

Algorithms for two-dimensional cutting stock and strip packing problems using dynamic programming and column generation

We investigate several two-dimensional guillotine cutting stock problems and their variants in which orthogonal rotations are allowed. We first present two dynamic programming based algorithms for the Rectangular Knapsack (RK) problem and its variants in which the patterns must be staged. The first algorithm solves the recurrence formula proposed by Beasley; the second algorithm - for staged patterns - also uses a recurrence formula. We show that if the items are not so small compared to the dimensions of the bin, then these algorithms require polynomial time. Using these algorithms we solved all instances of the RK problem found at the OR-LIBRARY, including one for which no optimal solution was known. We also consider the Two-dimensional Cutting Stock problem. We present a column generation based algorithm for this problem that uses the first algorithm above mentioned to generate the columns. We propose two strategies to tackle the residual instances. We also investigate a variant of this problem where the bins have different sizes. At last, we study the Two-dimensional Strip Packing problem. We also present a column generation based algorithm for this problem that uses the second algorithm above mentioned where staged patterns are imposed. In this case we solve instances for two-, three- and four-staged patterns. We report on some computational experiments with the various algorithms we propose in this paper. The results indicate that these algorithms seem to be suitable for solving real-world instances. We give a detailed description (a pseudo-code) of all the algorithms presented here, so that the reader may easily implement these algorithms. (c) 2007 Elsevier B.V. All rights reserved.

Epsilon Birnbaum-Saunders distribution family: properties and inference

In this paper we introduce a new extension for the Birnbaum-Saunder distribution based on the family of the epsilon-skew-symmetric distributions studied in Arellano-Valle et al. (J Stat Plan Inference 128(2):427-443, 2005). The extension allows generating Birnbaun-Saunders type distributions able to deal with extreme or outlying observations (Dupuis and Mills, IEEE Trans Reliab 47:88-95, 1998). Basic properties such as moments and Fisher information matrix are also studied. Results of a real data application are reported illustrating good fitting properties of the proposed model.

Confidence in heuristic solutions?

Solutions to combinatorial optimization problems frequently rely on heuristics to minimize an objective function. The optimum is sought iteratively and pre-setting the number of iterations dominates in operations research applications, which implies that the quality of the solution cannot be ascertained. Deterministic bounds offer a mean of ascertaining the quality, but such bounds are available for only a limited number of heuristics and the length of the interval may be difficult to control in an application. A small, almost dormant, branch of the literature suggests using statistical principles to derive statistical bounds for the optimum. We discuss alternative approaches to derive statistical bounds. We also assess their performance by testing them on 40 test p-median problems on facility location, taken from Beasley’s OR-library, for which the optimum is known. We consider three popular heuristics for solving such location problems; simulated annealing, vertex substitution, and Lagrangian relaxation where only the last offers deterministic bounds. Moreover, we illustrate statistical bounds in the location of 71 regional delivery points of the Swedish Post. We find statistical bounds reliable and much more efficient than deterministic bounds provided that the heuristic solutions are sampled close to the optimum. Statistical bounds are also found computationally affordable.

Quantitative Information Flow, Relations and Polymorphic Types

This paper uses Shannon's information theory to give a quantitative definition of information flow in systems that transform inputs to outputs. For deterministic systems, the definition is shown to specialise to a simpler form when the information source and the known inputs jointly determine the inputs. For this special case, the definition is related to the classical security condition of non-interference and an equivalence is established between non-interference and independence of random variables. Quantitative information flow for deterministic systems is then presented in relational form. With this presentation, it is shown how relational parametricity can be used to derive upper and lower bounds on information flows through families of functions defined in the second order lambda calculus.