Analysis of precipitation climatology simulated by high resolution coupled global models over the South America

The results of coupled high resolution global models (CGCMs) over South America are discussed. HiGEM1.2 and HadGEM1.2 simulations, with horizontal resolution of ~90 and 135 km, respectively, are compared. Precipitation estimations from CMAP (Climate Prediction Center—Merged Analysis of Precipitation), CPC (Climate Prediction Center) and GPCP (Global Precipitation Climatology Project) are used for validation. HiGEM1.2 and HadGEM1.2 simulated seasonal mean precipitation spatial patterns similar to the CMAP. The positioning and migration of the Intertropical Convergence Zone and of the Pacific and Atlantic subtropical highs are correctly simulated by the models. In HiGEM1.2 and HadGEM1.2, the intensity and locations of the South Atlantic Convergence Zone are in agreement with the observed dataset. The simulated annual cycles are in phase with estimations of rainfall for most of the six regions considered. An important result is that HiGEM1.2 and HadGEM1.2 eliminate a common problem of coarse resolution CGCMs, which is the simulation of a semiannual cycle of precipitation due to the semiannual solar forcing. Comparatively, the use of high resolution in HiGEM1.2 reduces the dry biases in the central part of Brazil during austral winter and spring and in most part of the year over an oceanic box in eastern Uruguay.

Approximation by harmonic polynomials in star-shaped domains and exponential convergence of Trefftz hp-dGFEM

We study the approximation of harmonic functions by means of harmonic polynomials in two-dimensional, bounded, star-shaped domains. Assuming that the functions possess analytic extensions to a delta-neighbourhood of the domain, we prove exponential convergence of the approximation error with respect to the degree of the approximating harmonic polynomial. All the constants appearing in the bounds are explicit and depend only on the shape-regularity of the domain and on delta. We apply the obtained estimates to show exponential convergence with rate O(exp(−b square root N)), N being the number of degrees of freedom and b>0, of a hp-dGFEM discretisation of the Laplace equation based on piecewise harmonic polynomials. This result is an improvement over the classical rate O(exp(−b cubic root N )), and is due to the use of harmonic polynomial spaces, as opposed to complete polynomial spaces.

Noisy Monte Carlo: convergence of Markov chains with approximate transition kernels

Monte Carlo algorithms often aim to draw from a distribution π by simulating a Markov chain with transition kernel P such that π is invariant under P. However, there are many situations for which it is impractical or impossible to draw from the transition kernel P. For instance, this is the case with massive datasets, where is it prohibitively expensive to calculate the likelihood and is also the case for intractable likelihood models arising from, for example, Gibbs random fields, such as those found in spatial statistics and network analysis. A natural approach in these cases is to replace P by an approximation Pˆ. Using theory from the stability of Markov chains we explore a variety of situations where it is possible to quantify how ’close’ the chain given by the transition kernel Pˆ is to the chain given by P . We apply these results to several examples from spatial statistics and network analysis.

Easterly wave disturbances over Northeast Brazil: an observational analysis

This paper aims to identify the circulation associated with Easterly Wave Disturbances (EWDs) that propagate toward the Eastern Northeast Brazil (ENEB) and their impact on the rainfall over ENEB during 2006 and 2007 rainy seasons (April–July). The EWDs identification and trajectory are analyzed using an automatic tracking technique (TracKH). The EWDs circulation patterns and their main features were obtained using the composite technique. To evaluate the TracKH efficiency, a validation was done by comparing the EWDs number tracked against observed cases obtained from an observational analysis. The mean characteristics of EWDs are 5.5-day period, propagation speed of ~9.5 m·s−1, and a 4500 km wavelength. A synoptic analysis shows that between days −2 d and 0 d, the low level winds presented cyclonic relative vorticity and convergence anomalies both in 2006 and 2007. The EWDs signals are strongest at low levels. The EWDs propagation is associated with relative humidity and precipitation positive anomalies and OLR and omega negative anomalies. The EWDs tracks are seen over all ENEB and their lysis occurs between the ENEB and marginally inside the continent. The tracking captured 71% of EWDs in all periods, indicating that an objective analysis is a promising method for EWDs detection.

A priori error analysis of space–time Trefftz discontinuous Galerkin methods for wave problems

We present and analyse a space–time discontinuous Galerkin method for wave propagation problems. The special feature of the scheme is that it is a Trefftz method, namely that trial and test functions are solution of the partial differential equation to be discretised in each element of the (space–time) mesh. The method considered is a modification of the discontinuous Galerkin schemes of Kretzschmar et al. (2014) and of Monk & Richter (2005). For Maxwell’s equations in one space dimension, we prove stability of the method, quasi-optimality, best approximation estimates for polynomial Trefftz spaces and (fully explicit) error bounds with high order in the meshwidth and in the polynomial degree. The analysis framework also applies to scalar wave problems and Maxwell’s equations in higher space dimensions. Some numerical experiments demonstrate the theoretical results proved and the faster convergence compared to the non-Trefftz version of the scheme.

Impact of a Fram Strait cyclone on ice edge, drift, divergence, and concentration: Possibilities and limits of an observational analysis

This study aims at the determination of a Fram Strait cyclone track and of the cyclone’s impact on ice edge, drift, divergence, and concentration. A 24 h period on 13–14 March 2002 framed by two RADARSAT images is analyzed. Data are included from autonomous ice buoys, a research vessel, Special Sensor Microwave Imager (SSM/I) and QuikSCAT satellite, and the operational European Centre for Medium-Range Weather Forecasts (ECMWF) model. During this 24 h period the cyclone moved northward along the western ice edge in the Fram Strait, crossed the northern ice edge, made a left-turn loop with 150 km diameter over the sea ice, and returned to the northern ice edge. The ECMWF analysis places the cyclone track 100 km too far west over the sea ice, a deviation which is too large for representative sea ice simulations. On the east side of the northward moving cyclone, the ice edge was pushed northward by 55 km because of strong winds. On the rear side, the ice edge advanced toward the open water but by a smaller distance because of weaker winds there. The ice drift pattern as calculated from the ice buoys and the two RADARSAT images is cyclonically curved around the center of the cyclone loop. Ice drift divergence shows a spatial pattern with divergence in the loop center and a zone of convergence around. Ice concentration changes as retrieved from SSM/I data follow the divergence pattern such that sea ice concentration increased in areas of divergence and decreased in areas of convergence.

Non-stationary frequency analysis of extreme daily rainfall in Sao Paulo, Brazil

This work is an assessment of frequency of extreme values (EVs) of daily rainfall in the city of Sao Paulo. Brazil, over the period 1933-2005, based on the peaks-over-threshold (POT) and Generalized Pareto Distribution (GPD) approach. Usually. a GPD model is fitted to a sample of POT Values Selected With a constant threshold. However. in this work we use time-dependent thresholds, composed of relatively large p quantities (for example p of 0.97) of daily rainfall amounts computed from all available data. Samples of POT values were extracted with several Values of p. Four different GPD models (GPD-1, GPD-2, GPD-3. and GDP-4) were fitted to each one of these samples by the maximum likelihood (ML) method. The shape parameter was assumed constant for the four models, but time-varying covariates were incorporated into scale parameter of GPD-2. GPD-3, and GPD-4, describing annual cycle in GPD-2. linear trend in GPD-3, and both annual cycle and linear trend in GPD-4. The GPD-1 with constant scale and shape parameters is the simplest model. For identification of the best model among the four models WC used rescaled Akaike Information Criterion (AIC) with second-order bias correction. This criterion isolates GPD-3 as the best model, i.e. the one with positive linear trend in the scale parameter. The slope of this trend is significant compared to the null hypothesis of no trend, for about 98% confidence level. The non-parametric Mann-Kendall test also showed presence of positive trend in the annual frequency of excess over high thresholds. with p-value being virtually zero. Therefore. there is strong evidence that high quantiles of daily rainfall in the city of Sao Paulo have been increasing in magnitude and frequency over time. For example. 0.99 quantiles of daily rainfall amount have increased by about 40 mm between 1933 and 2005. Copyright (C) 2008 Royal Meteorological Society

DENSITY PROFILE, VELOCITY ANISOTROPY, AND LINE-OF-SIGHT EXTERNAL CONVERGENCE OF SLACS GRAVITATIONAL LENSES

Data from 58 strong-lensing events surveyed by the Sloan Lens ACS Survey are used to estimate the projected galaxy mass inside their Einstein radii by two independent methods: stellar dynamics and strong gravitational lensing. We perform a joint analysis of these two estimates inside models with up to three degrees of freedom with respect to the lens density profile, stellar velocity anisotropy, and line-of-sight (LOS) external convergence, which incorporates the effect of the large-scale structure on strong lensing. A Bayesian analysis is employed to estimate the model parameters, evaluate their significance, and compare models. We find that the data favor Jaffe`s light profile over Hernquist`s, but that any particular choice between these two does not change the qualitative conclusions with respect to the features of the system that we investigate. The density profile is compatible with an isothermal, being sightly steeper and having an uncertainty in the logarithmic slope of the order of 5% in models that take into account a prior ignorance on anisotropy and external convergence. We identify a considerable degeneracy between the density profile slope and the anisotropy parameter, which largely increases the uncertainties in the estimates of these parameters, but we encounter no evidence in favor of an anisotropic velocity distribution on average for the whole sample. An LOS external convergence following a prior probability distribution given by cosmology has a small effect on the estimation of the lens density profile, but can increase the dispersion of its value by nearly 40%.

Multivariate analysis of the energy cycle of the South American rainy season

The South American (SA) rainy season is studied in this paper through the application of a multivariate Empirical Orthogonal Function (EOF) analysis to a SA gridded precipitation analysis and to the components of Lorenz Energy Cycle (LEC) derived from the National Centers for Environmental Prediction (NCEP) reanalysis. The EOF analysis leads to the identification of patterns of the rainy season and the associated mechanisms in terms of their energetics. The first combined EOF represents the northwest-southeast dipole of the precipitation between South and Central America, the South American Monsoon System (SAMS). The second combined EOF represents a synoptic pattern associated with the SACZ (South Atlantic convergence zone) and the third EOF is in spatial quadrature to the second EOF. The phase relationship of the EOFs, as computed from the principal components (PCs), suggests a nonlinear transition from the SACZ to the fully developed SAMS mode by November and between both components describing the SACZ by September-October (the rainy season onset). According to the LEC, the first mode is dominated by the eddy generation term at its maximum, the second by both baroclinic and eddy generation terms and the third by barotropic instability previous to the connection to the second mode by September-October. The predominance of the different LEC components at each phase of the SAMS can be used as an indicator of the onset of the rainy season in terms of physical processes, while the existence of the outstanding spectral peaks in the time dependence of the EOFs at the intraseasonal time scale could be used for monitoring purposes. Copyright (C) 2009 Royal Meteorological Society

Analysis of a hyper-diverse seed dispersal network: modularity and underlying mechanisms

Mutualistic interactions involving pollination and ant-plant mutualistic networks typically feature tightly linked species grouped in modules. However, such modularity is infrequent in seed dispersal networks, presumably because research on those networks predominantly includes a single taxonomic animal group (e.g. birds). Herein, for the first time, we examine the pattern of interaction in a network that includes multiple taxonomic groups of seed dispersers, and the mechanisms underlying modularity. We found that the network was nested and modular, with five distinguishable modules. Our examination of the mechanisms underlying such modularity showed that plant and animal trait values were associated with specific modules but phylogenetic effect was limited. Thus, the pattern of interaction in this network is only partially explained by shared evolutionary history. We conclude that the observed modularity emerged by a combination of phylogenetic history and trait convergence of phylogenetically unrelated species, shaped by interactions with particular types of dispersal agents.

The Crystal Complex of Phosphofructokinase-2 of Escherichia coli with Fructose-6-phosphate KINETIC AND STRUCTURAL ANALYSIS OF THE ALLOSTERIC ATP INHIBITION

Substrate inhibition by ATP is a regulatory feature of the phosphofructokinases isoenzymes from Escherichia coli (Pfk-1 and Pfk-2). Under gluconeogenic conditions, the loss of this regulation in Pfk-2 causes substrate cycling of fructose-6-phosphate (fructose-6-P) and futile consumption of ATP delaying growth. In the present work, we have broached the mechanism of ATP-induced inhibition of Pfk-2 from both structural and kinetic perspectives. The crystal structure of Pfk-2 in complex with fructose-6-P is reported to a resolution of 2 angstrom. The comparison of this structure with the previously reported inhibited form of the enzyme suggests a negative interplay between fructose-6-P binding and allosteric binding of MgATP. Initial velocity experiments show a linear increase of the apparent K(0.5) for fructose-6-P and a decrease in the apparent k(cat) as a function of MgATP concentration. These effects occur simultaneously with the induction of a sigmoidal kinetic behavior (n(H) of approximately 2). Differences and resemblances in the patterns of fructose-6-P binding and the mechanism of inhibition are discussed for Pfk-1 and Pfk-2, as an example of evolutionary convergence, because these enzymes do not share a common ancestor.

STABILITY ANALYSIS WITH APPLICATIONS OF A TWO-DIMENSIONAL DYNAMICAL SYSTEM ARISING FROM A STOCHASTIC MODEL FOR AN ASSET MARKET

We analyze the stability properties of equilibrium solutions and periodicity of orbits in a two-dimensional dynamical system whose orbits mimic the evolution of the price of an asset and the excess demand for that asset. The construction of the system is grounded upon a heterogeneous interacting agent model for a single risky asset market. An advantage of this construction procedure is that the resulting dynamical system becomes a macroscopic market model which mirrors the market quantities and qualities that would typically be taken into account solely at the microscopic level of modeling. The system`s parameters correspond to: (a) the proportion of speculators in a market; (b) the traders` speculative trend; (c) the degree of heterogeneity of idiosyncratic evaluations of the market agents with respect to the asset`s fundamental value; and (d) the strength of the feedback of the population excess demand on the asset price update increment. This correspondence allows us to employ our results in order to infer plausible causes for the emergence of price and demand fluctuations in a real asset market. The employment of dynamical systems for studying evolution of stochastic models of socio-economic phenomena is quite usual in the area of heterogeneous interacting agent models. However, in the vast majority of the cases present in the literature, these dynamical systems are one-dimensional. Our work is among the few in the area that construct and study analytically a two-dimensional dynamical system and apply it for explanation of socio-economic phenomena.

A note on convergence of Peck-Shell and Green-Lin mechanisms in the Diamond-Dybvig model

We study the effects of population size in the Peck-Shell analysis of bank runs. We find that a contract featuring equal-treatment for almost all depositors of the same type approximates the optimum. Because the approximation also satisfies Green-Lin incentive constraints, when the planner discloses positions in the queue, welfare in these alternative specifications are sandwiched. Disclosure, however, it is not needed since our approximating contract is not subject to runs.

Convergence in Brazil: recent trends and long run prospects

This paper applies to the analysis of the interstate income distribution in BraziI a set of techniques that have been widely used in the current empirical literature on growth and convergence. Usual measures of dispersion in the interstate income distribution (the coefficient of variation and Theil' s index) suggest that cr-convergence was an unequivoca1 feature of the regional growth experience in BraziI, between 1970 and 1986. After 1986, the process of convergence seems, however, to have sIowed down almost to a halt. A standard growth modeI is shown to fit the regional data well and to expIain a substantial amount of the variation in growth rates, providing estimates of the speed of (conditional) J3-convergence of approximateIy 3% p.a .. Different estimates of the long run distribution implied by the recent growth trends point towards further reductions in the interstate income inequality, but also suggest that the relative per capita incomes of a significant number of states and the number of ''very poor" and "poor" states were, in 1995, already quite c10se to their steady-state values.