The importance of vertical velocity variability for estimates of the indirect aerosol effects

The activation of aerosols to form cloud droplets is dependent upon vertical velocities whose local variability is not typically resolved at the GCM grid scale. Consequently, it is necessary to represent the subgrid-scale variability of vertical velocity in the calculation of cloud droplet number concentration. This study uses the UK Chemistry and Aerosols community model (UKCA) within the Hadley Centre Global Environmental Model (HadGEM3), coupled for the first time to an explicit aerosol activation parameterisation, and hence known as UKCA-Activate. We explore the range of uncertainty in estimates of the indirect aerosol effects attributable to the choice of parameterisation of the subgrid-scale variability of vertical velocity in HadGEM-UKCA. Results of simulations demonstrate that the use of a characteristic vertical velocity cannot replicate results derived with a distribution of vertical velocities, and is to be discouraged in GCMs. This study focuses on the effect of the variance (σw2) of a Gaussian pdf (probability density function) of vertical velocity. Fixed values of σw (spanning the range measured in situ by nine flight campaigns found in the literature) and a configuration in which σw depends on turbulent kinetic energy are tested. Results from the mid-range fixed σw and TKE-based configurations both compare well with observed vertical velocity distributions and cloud droplet number concentrations. The radiative flux perturbation due to the total effects of anthropogenic aerosol is estimated at −1.9 W m−2 with σw = 0.1 m s−1, −2.1 W m−2 with σw derived from TKE, −2.25 W m−2 with σw = 0.4 m s−1, and −2.3 W m−2 with σw = 0.7 m s−1. The breadth of this range is 0.4 W m−2, which is comparable to a substantial fraction of the total diversity of current aerosol forcing estimates. Reducing the uncertainty in the parameterisation of σw would therefore be an important step towards reducing the uncertainty in estimates of the indirect aerosol effects. Detailed examination of regional radiative flux perturbations reveals that aerosol microphysics can be responsible for some climate-relevant radiative effects, highlighting the importance of including microphysical aerosol processes in GCMs.

Analysis of the regional pattern of sea level change due to ocean dynamics and density changes for 1993-2099 in observations and CMIP5 AOGCMs

Predictions of twenty-first century sea level change show strong regional variation. Regional sea level change observed by satellite altimetry since 1993 is also not spatially homogenous. By comparison with historical and pre-industrial control simulations using the atmosphere–ocean general circulation models (AOGCMs) of the CMIP5 project, we conclude that the observed pattern is generally dominated by unforced (internal generated) variability, although some regions, especially in the Southern Ocean, may already show an externally forced response. Simulated unforced variability cannot explain the observed trends in the tropical Pacific, but we suggest that this is due to inadequate simulation of variability by CMIP5 AOGCMs, rather than evidence of anthropogenic change. We apply the method of pattern scaling to projections of sea level change and show that it gives accurate estimates of future local sea level change in response to anthropogenic forcing as simulated by the AOGCMs under RCP scenarios, implying that the pattern will remain stable in future decades. We note, however, that use of a single integration to evaluate the performance of the pattern-scaling method tends to exaggerate its accuracy. We find that ocean volume mean temperature is generally a better predictor than global mean surface temperature of the magnitude of sea level change, and that the pattern is very similar under the different RCPs for a given model. We determine that the forced signal will be detectable above the noise of unforced internal variability within the next decade globally and may already be detectable in the tropical Atlantic.

An Ensemble Smoother with Error Estimates

A smoother introduced earlier by van Leeuwen and Evensen is applied to a problem in which real obser vations are used in an area with strongly nonlinear dynamics. The derivation is new , but it resembles an earlier derivation by van Leeuwen and Evensen. Again a Bayesian view is taken in which the prior probability density of the model and the probability density of the obser vations are combined to for m a posterior density . The mean and the covariance of this density give the variance-minimizing model evolution and its errors. The assumption is made that the prior probability density is a Gaussian, leading to a linear update equation. Critical evaluation shows when the assumption is justified. This also sheds light on why Kalman filters, in which the same ap- proximation is made, work for nonlinear models. By reference to the derivation, the impact of model and obser vational biases on the equations is discussed, and it is shown that Bayes’ s for mulation can still be used. A practical advantage of the ensemble smoother is that no adjoint equations have to be integrated and that error estimates are easily obtained. The present application shows that for process studies a smoother will give superior results compared to a filter , not only owing to the smooth transitions at obser vation points, but also because the origin of features can be followed back in time. Also its preference over a strong-constraint method is highlighted. Further more, it is argued that the proposed smoother is more efficient than gradient descent methods or than the representer method when error estimates are taken into account

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%.

The impact of fragmentation and density regulation on forest succession in the Atlantic rain forest

The Atlantic Rain Forest, an important biodiversity hot spot, has faced severe habitat loss since the last century which has resulted in a highly fragmented landscape with a large number of small forest patches (<100 ha). For conservation planning it is essential to understand how current and future forest regeneration depends on ecological processes, fragment size and the connection to the regional seed pool. We have investigated the following questions by applying the forest growth simulation model FORMIND to the situation of the Atlantic Forest in the state of Sao Paulo, SE Brazil: (1) which set of parameters describing the local regeneration and level of density regulation can reproduce the biomass distribution and stem density of an old growth forest in a reserve? (2) Which additional processes apart from those describing the dynamics of an old growth forest, drive forest succession of small isolated fragments? (3) Which role does external seed input play during succession? Therefore, more than 300 tree species have been classified into nine plant functional types (PFTs), which are characterized by maximum potential height and shade tolerance. We differentiate between two seed dispersal modes: (i) local dispersal, i.e. all seedlings originated from fertile trees within the simulated area and (ii) external seed rain. Local seed dispersal has been parameterized following the pattern oriented approach, using biomass estimates of old growth forest. We have found that moderate density regulation is essential to achieve coexistence for a broad range of regeneration parameters. Considering the expected uncertainty and variability in the regeneration processes it is important that the forest dynamics are robust to variations in the regeneration parameters. Furthermore, edge effects such as increased mortality at the border and external seed rain have been necessary to reproduce the patterns for small isolated fragments. Overall, simulated biomass is much lower in the fragments compared to the continuous forest, whereas shade tolerant species are affected most strongly by fragmentation. Our simulations can supplement empirical studies by extrapolating local knowledge on edge effects of fragments to larger temporal and spatial scales. In particular our results show the importance of external seed rain and therefore highlight the importance of structural connectivity between regenerating fragments and mature forest stands. (C) 2009 Elsevier B.V. All rights reserved.

Nonparametric specification tests for conditional duration models

This paper deals with the testing of autoregressive conditional duration (ACD) models by gauging the distance between the parametric density and hazard rate functions implied by the duration process and their non-parametric estimates. We derive the asymptotic justification using the functional delta method for fixed and gamma kernels, and then investigate the finite-sample properties through Monte Carlo simulations. Although our tests display some size distortion, bootstrapping suffices to correct the size without compromising their excellent power. We show the practical usefulness of such testing procedures for the estimation of intraday volatility patterns.

Central limit theorem for asymmetric kernel functionals

Asymmetric kernels are quite useful for the estimation of density functions with bounded support. Gamma kernels are designed to handle density functions whose supports are bounded from one end only, whereas beta kernels are particularly convenient for the estimation of density functions with compact support. These asymmetric kernels are nonnegative and free of boundary bias. Moreover, their shape varies according to the location of the data point, thus also changing the amount of smoothing. This paper applies the central limit theorem for degenerate U-statistics to compute the limiting distribution of a class of asymmetric kernel functionals.

The aliasing effect, the Fejer Kernel and temporally aggregated long memory processes

This paper derives the spectral density function of aggregated long memory processes in light of the aliasing effect. The results are different from previous analyses in the literature and a small simulation exercise provides evidence in our favour. The main result point to that flow aggregates from long memory processes shall be less biased than stock ones, although both retain the degree of long memory. This result is illustrated with the daily US Dollar/ French Franc exchange rate series.

Estimadores do tipo núcleo para Variávei s I.I.D. com espaço de estados geral

In this work, the paper of Campos and Dorea [3] was detailed. In that article a Kernel Estimator was applied to a sequence of random variables with general state space, which were independent and identicaly distributed. In chapter 2, the estimator´s properties such as asymptotic unbiasedness, consistency in quadratic mean, strong consistency and asymptotic normality were verified. In chapter 3, using R software, numerical experiments were developed in order to give a visual idea of the estimate process

Estimador do Tipo Núcleo para Densidades Limites de Cadeias de Markov com Espaço de Estados Geral

In this work we studied the consistency for a class of kernel estimates of f f (.) in the Markov chains with general state space E C Rd case. This study is divided into two parts: In the first one f (.) is a stationary density of the chain, and in the second one f (x) v (dx) is the limit distribution of a geometrically ergodic chain

Estimadores do tipo núcleo para a densidade de transição de uma cadeia de markov com espaço de estados geral

In this work, we studied the strong consistency for a class of estimates for a transition density of a Markov chain with general state space E ⊂ Rd. The strong ergodicity of the estimates for the density transition is obtained from the strong consistency of the kernel estimates for both the marginal density p(:) of the chain and the joint density q(., .). In this work the Markov chain is supposed to be homogeneous, uniformly ergodic and possessing a stationary density p(.,.)

Spatial distribution of trachoma cases in the City of Bauru, State of São Paulo, Brazil, detected in 2006: defining key areas for improvement of health resources

Introduction: The objective of this study was to analyze the spatial behavior of the occurrence of trachoma cases detected in the City of Bauru, State of São Paulo, Brazil, in 2006 in order to use the information collected to set priority areas for optimization of health resources. Methods: the trachoma cases identified in 2006 were georeferenced. The data evaluated were: schools where the trachoma cases studied, data from the 2000 Census, census tract, type of housing, water supply conditions, distribution of income and levels of education of household heads. In the Google Earth® software and TerraView® were made descriptive spatial analysis and estimates of the Kernel. Each area was studied by interpolation of the density surfaces exposing events to facilitate to recognize the clusters. Results: of the 66 cases detected, only one (1.5%) was not a resident of the city's outskirts. A positive association was detected of trachoma cases and the percentage of heads of household with income below three minimum wages and schooling under eight years of education. Conclusions: The recognition of the spatial distribution of trachoma cases coincided with the areas of greatest social inequality in Bauru City. The micro-areas identified are those that should be prioritized in the rationalization of health resources. There is the possibility of using the trachoma cases detected as an indicator of performance of micro priority health programs.

Comparison of two methods (visual estimates and filming) for counts of horn flies (Haematobia irritans irritans) (L.) (Diptera : Muscidae)

Visual estimates are generally used for counts of horn flies, Haematobia irritans (L.) and play an important role as an instrument to quantify fly populations in scientific studies. In this study, horn fly counts were performed on 30 Nelore steers in the municipality of Aracatuba, SP Brazil, from January to December 1998. Flies were counted weekly by two methods: the estimate method whereby estimates of the number of flies on one side of the animal are obtained by visual observation, and the filming method whereby images of flies from both sides of the animal are recorded with a video camera. The tape was then played on a videotape recorder coupled to a television and the flies were counted on the screen. Both methods showed variations in horn fly population density during the period studied. However, significant differences (p < 0.05) were observed between the two methods with the filming method permitting the visualization of a larger number of flies than the estimate method. In addition, the filming method permitted safe and reliable counts hours after the images were taken, with the advantage that the tape can serve as an archive for random re-counts. (C) 2002 Elsevier B.V. B.V. All rights reserved.

Theoretical estimates of consumable food and probability of acquiring food in larvae of Chrysomya putoria (Diptera: Calliphoridae)

An indirect estimate of consumable food and probability of acquiring food in a blowfly species, Chrysomya putoria, is presented. This alternative procedure combines three distinct models to estimate consumable food in the context of the exploitative competition experienced by immature individuals in blowfly populations. The relevant parameters are derived from data for pupal weight and survival and estimates of density-independent larval mortality in twenty different larval densities. As part of this procedure, the probability of acquiring food per unit of time and the time taken to exhaust the food supply are also calculated. The procedure employed here may be valuable for estimations in insects whose immature stages develop inside the food substrate, where it is difficult to partial out confounding effects such as separation of faeces. This procedure also has the advantage of taking into account the population dynamics of immatures living under crowded conditions, which are particularly characteristic of blowflies and other insects as well.

Sampling density and proportion for the characterization of the variability of Oxisol attributes on different materials

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