Validation of the GRAPE single view aerosol retrieval for ATSR-2 and insights into the long term global AOD trend over the ocean

The Global Retrieval of ATSR Cloud Parameters and Evaluation (GRAPE) project has produced a global data-set of cloud and aerosol properties from the Along Track Scanning Radiometer-2 (ATSR-2) instrument, covering the time period 1995�2001. This paper presents the validation of aerosol optical depths (AODs) over the ocean from this product against AERONET sun-photometer measurements, as well as a comparison to the Advanced Very High Resolution Radiometer (AVHRR) optical depth product produced by the Global Aerosol Climatology Project (GACP). The GRAPE AOD over ocean is found to be in good agreement with AERONET measurements, with a Pearson's correlation coefficient of 0.79 and a best-fit slope of 1.0±0.1, but with a positive bias of 0.08±0.04. Although the GRAPE and GACP datasets show reasonable agreement, there are significant differences. These discrepancies are explored, and suggest that the downward trend in AOD reported by GACP may arise from changes in sampling due to the orbital drift of the AVHRR instruments.

Fingerprints of changes in annual and seasonal precipitation from CMIP5 models over land and ocean

By comparing annual and seasonal changes in precipitation over land and ocean since 1950 simulated by the CMIP5 (Coupled Model Intercomparison Project, phase 5) climate models in which natural and anthropogenic forcings have been included, we find that clear global-scale and regional-scale changes due to human influence are expected to have occurred over both land and ocean. These include moistening over northern high latitude land and ocean throughout all seasons and over the northern subtropical oceans during boreal winter. However we show that this signal of human influence is less distinct when considered over the relatively small area of land for which there are adequate observations to make assessments of multi-decadal scale trends. These results imply that extensive and significant changes in precipitation over the land and ocean may have already happened, even though, inadequacies in observations in some parts of the world make it difficult to identify conclusively such a human fingerprint on the global water cycle. In some regions and seasons, due to aliasing of different kinds of variability as a result of sub sampling by the sparse and changing observational coverage, observed trends appear to have been increased, underscoring the difficulties of interpreting the apparent magnitude of observed changes in precipitation.

Simultaneous lidar and EOS MLS measurements, and modeling, of a rare polar ozone filament event over Mauna Loa Observatory, Hawaii

In mid-March 2005, a rare lower stratospheric polar vortex filamentation event was observed simultaneously by the JPL lidar at Mauna Loa Observatory, Hawaii, and by the EOS MLS instrument onboard the Aura satellite. The event coincided with the beginning of the spring 2005 final warming. On 16 March, the filament was observed by lidar around 0600 UT between 415 K and 455 K, and by MLS six hours earlier. It was seen on both the lidar and MLS profiles as a layer of enhanced ozone, peaking at 1.7 ppmv in a region where the climatological values are usually around or below 1 ppmv. Ozone profiles measured by lidar and MLS were compared to profiles from the Chemical Transport Model MIMOSA-CHIM. The agreement between lidar, MLS, and the model is excellent considering the difference in the sampling techniques. MLS was also able to identify the filament at another location north of Hawaii.

Statistical uncertainty of changes in winter storms over the North Atlantic and Europe in an ensemble of transient climate simulations

Winter storms are among the most important natural hazards affecting Europe. We quantify changes in storm frequency and intensity over the North Atlantic and Europe under future climate scenarios in terms of return periods (RPs) considering uncertainties due to both sampling and methodology. RPs of North Atlantic storms' minimum central pressure (CP) and maximum vorticity (VOR) remain unchanged by 2100 for both the A1B and A2 scenarios compared to the present climate. Whereas shortened RPs for VOR of all intensities are detected for the area between British Isles/North-Sea/western Europe as early as 2040. However, the changes in storm VOR RP may be unrealistically large: a present day 50 (20) year event becomes approximately a 9 (5.5) year event in both A1B and A2 scenarios by 2100. The detected shortened RPs of storms implies a higher risk of occurrence of damaging wind events over Europe.

Vertical distribution and radiative effects of mineral dust and biomass burning aerosol over West Africa during DABEX

This paper presents measurements of the vertical distribution of aerosol extinction coefficient over West Africa during the Dust and Biomass-burning Aerosol Experiment (DABEX)/African Monsoon Multidisciplinary Analysis dry season Special Observing Period Zero (AMMA-SOP0). In situ aircraft measurements from the UK FAAM aircraft have been compared with two ground-based lidars (POLIS and ARM MPL) and an airborne lidar on an ultralight aircraft. In general, mineral dust was observed at low altitudes (up to 2 km), and a mixture of biomass burning aerosol and dust was observed at altitudes of 2–5 km. The study exposes difficulties associated with spatial and temporal variability when intercomparing aircraft and ground measurements. Averaging over many profiles provided a better means of assessing consistent errors and biases associated with in situ sampling instruments and retrievals of lidar ratios. Shortwave radiative transfer calculations and a 3-year simulation with the HadGEM2-A climate model show that the radiative effect of biomass burning aerosol was somewhat sensitive to the vertical distribution of aerosol. In particular, when the observed low-level dust layer was included in the model, the absorption of solar radiation by the biomass burning aerosols increased by 10%. We conclude that this absorption enhancement was caused by the dust reflecting solar radiation up into the biomass burning aerosol layer. This result illustrates that the radiative forcing of anthropogenic absorbing aerosol can be sensitive to the presence of natural aerosol species.

An assessment of observed vertical flux divergence in long-term eddy-covariance measurements over two Midwestern forest ecosystems

Vertical divergence of CO2 fluxes is observed over two Midwestern AmeriFlux forest sites. The differences in ensemble averaged hourly CO2 fluxes measured at two heights above canopy are relatively small (0.2–0.5 μmol m−2 s−1), but they are the major contributors to differences (76–256 g C m−2 or 41.8–50.6%) in estimated annual net ecosystem exchange (NEE) in 2001. A friction velocity criterion is used in these estimates but mean flow advection is not accounted for. This study examines the effects of coordinate rotation, averaging time period, sampling frequency and co-spectral correction on CO2 fluxes measured at a single height, and on vertical flux differences measured between two heights. Both the offset in measured vertical velocity and the downflow/upflow caused by supporting tower structures in upwind directions lead to systematic over- or under-estimates of fluxes measured at a single height. An offset of 1 cm s−1 and an upflow/downflow of 1° lead to 1% and 5.6% differences in momentum fluxes and nighttime sensible heat and CO2 fluxes, respectively, but only 0.5% and 2.8% differences in daytime sensible heat and CO2 fluxes. The sign and magnitude of both offset and upflow/downflow angle vary between sonic anemometers at two measurement heights. This introduces a systematic and large bias in vertical flux differences if these effects are not corrected in the coordinate rotation. A 1 h averaging time period is shown to be appropriate for the two sites. In the daytime, the absolute magnitudes of co-spectra decrease with height in the natural frequencies of 0.02–0.1 Hz but increase in the lower frequencies (<0.01 Hz). Thus, air motions in these two frequency ranges counteract each other in determining vertical flux differences, whose magnitude and sign vary with averaging time period. At night, co-spectral densities of CO2 are more positive at the higher levels of both sites in the frequency range of 0.03–0.4 Hz and this vertical increase is also shown at most frequencies lower than 0.03 Hz. Differences in co-spectral corrections at the two heights lead to a positive shift in vertical CO2 flux differences throughout the day at both sites. At night, the vertical CO2 flux differences between two measurement heights are 20–30% and 40–60% of co-spectral corrected CO2 fluxes measured at the lower levels of the two sites, respectively. Vertical differences of CO2 flux are relatively small in the daytime. Vertical differences in estimated mean vertical advection of CO2 between the two measurement heights generally do not improve the closure of the 1D (vertical) CO2 budget in the air layer between the two measurement heights. This may imply the significance of horizontal advection. However, a reliable assessment of mean advection contributions in annual NEE estimate at these two AmeriFlux sites is currently an unsolved problem.

Two fast radiative transfer methods to improve the temporal sampling of clouds in numerical weather prediction and climate models

The high computational cost of calculating the radiative heating rates in numerical weather prediction (NWP) and climate models requires that calculations are made infrequently, leading to poor sampling of the fast-changing cloud field and a poor representation of the feedback that would occur. This paper presents two related schemes for improving the temporal sampling of the cloud field. Firstly, the ‘split time-stepping’ scheme takes advantage of the independent nature of the monochromatic calculations of the ‘correlated-k’ method to split the calculation into gaseous absorption terms that are highly dependent on changes in cloud (the optically thin terms) and those that are not (optically thick). The small number of optically thin terms can then be calculated more often to capture changes in the grey absorption and scattering associated with cloud droplets and ice crystals. Secondly, the ‘incremental time-stepping’ scheme uses a simple radiative transfer calculation using only one or two monochromatic calculations representing the optically thin part of the atmospheric spectrum. These are found to be sufficient to represent the heating rate increments caused by changes in the cloud field, which can then be added to the last full calculation of the radiation code. We test these schemes in an operational forecast model configuration and find a significant improvement is achieved, for a small computational cost, over the current scheme employed at the Met Office. The ‘incremental time-stepping’ scheme is recommended for operational use, along with a new scheme to correct the surface fluxes for the change in solar zenith angle between radiation calculations.

Differentiation in populations of the apple scab fungus Venturia inaequalis on cultivars in a mixed orchard remains over time

The ascomycete Venturia inaequalis causes annual epidemics of apple scab worldwide. Scab development is reduced in mixed cultivar orchards compared with monocultures. To use mixtures in commercial production, we need to understand how the population of scab changes in a mixed orchard and how likely a super race, with virulence factors overcoming multiple resistance factors in the mixed orchard, is to emerge and become dominant. We used short sequence repeat (SSR) markers to investigate the temporal change of scab populations in two mixed cultivar orchards in the UK to infer the likelihood of emergence of a scab super race. There were no significant differences between the populations at the two sampling times (six or seven years apart) in either of the two mixed orchards. In one of the orchards apple scab populations on different cultivars were significantly different and the differences did not diminish over time. These results suggest that it is not inevitable that a super race of V. inaequalis will become dominant during the lifetime of a commercial apple orchard.

Changes of a mutualistic network over time: reanalysis over a 10-year period

We analyzed the structure of a multispecific network or interacting ants and plants bearing extrafloral nectaries recorded in 1990 and again in 2000 in La Mancha, Veracruz, Mexico. We assessed the replicability of the number of interactions found among species and also whether there had been changes in the network structure associated with appearance of new ant and plant species during. that 10-year period. Our results show that the nested topology of the network was similar between sampling dates, group dissimilarity increased, mean number of interactions for ant species increased, the frequency distribution of standardized degrees reached higher values for plant species, more ant species and fewer plant species constituted the core of the more recent network, and the presence of new ant and plant species increased while their contribution to nestedness remained the same. Generalist species (i.e., those with the most links or interactions) appeared to maintain the stability of the network because the new species incorporated into the communities were linked to this core of generalists. Camponotus planatus was the most extreme generalist ant species (the one with the most links) in both networks, followed by four other ant species; but other species changed either their position along the continuum of generalists relative to specialists or their presence or absence within the network. Even though new species moved into the area during the decade between the surveys, the overall network structure remained unmodified.

Sampling, WLS, and Mixed Models

Mixed models may be defined with or without reference to sampling, and can be used to predict realized random effects, as when estimating the latent values of study subjects measured with response error. When the model is specified without reference to sampling, a simple mixed model includes two random variables, one stemming from an exchangeable distribution of latent values of study subjects and the other, from the study subjects` response error distributions. Positive probabilities are assigned to both potentially realizable responses and artificial responses that are not potentially realizable, resulting in artificial latent values. In contrast, finite population mixed models represent the two-stage process of sampling subjects and measuring their responses, where positive probabilities are only assigned to potentially realizable responses. A comparison of the estimators over the same potentially realizable responses indicates that the optimal linear mixed model estimator (the usual best linear unbiased predictor, BLUP) is often (but not always) more accurate than the comparable finite population mixed model estimator (the FPMM BLUP). We examine a simple example and provide the basis for a broader discussion of the role of conditioning, sampling, and model assumptions in developing inference.

Feasibility of using in situ fusion for the determination of Co, Cr and Mn in Portland cement by direct solid sampling graphite furnace atomic absorption spectrometry

In situ fusion on the boat-type graphite platform has been used as a sample pretreatment for the direct determination of Co, Cr and Mn in Portland cement by solid sampling graphite furnace atomic absorption spectrometry (SS-GF AAS). The 3-field Zeeman technique was adopted for background correction to decrease the sensitivity during measurements. This strategy allowed working with up to 200 mu g of sample. The in situ fusion was accomplished using 10 mu L of a flux mixture 4.0% m/v Na(2)CO(3) + 4.0% m/v ZnO + 0.1% m/v Triton (R) X-100 added over the cement sample and heated at 800 degrees C for 20 s. The resulting mould was completely dissolved with 10 mu L of 0.1% m/v HNO(3). Limits of detection were 0.11 mu g g(-1) for Co, 1.1 mu g g(-1) for Cr and 1.9 mu g g(-1) for Mn. The accuracy of the proposed method has been evaluated by the analysis of certified reference materials. The values found presented no statistically significant differences compared to the certified values (Student`s t-test, p<0.05). In general, the relative standard deviation was lower than 12% (n = 5). (C) 2009 Elsevier B.V. All rights reserved.

Occurrence and behavior of butyltins in intertidal and shallow subtidal surface sediments of an estuarine beach under different sampling conditions

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

Economic design of (X)over-bar and R charts under Weibull shock models

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

AATS for the (X)over-bar chart with variable parameters

A Fortran computer program is given for the computation of the adjusted average time to signal, or AATS, for adaptive (X) over bar charts with one, two, or all three design parameters variable: the sample size, n, the sampling interval, h, and the factor k used in determining the width of the action limits. The program calculates the threshold limit to switch the adaptive design parameters and also provides the in-control average time to signal, or ATS.

Constrained optimization model for the design of an adaptive (X)over-bar chart

An economic-statistical model is developed for variable parameters (VP) (X) over bar charts in which all design parameters vary adaptively, that is, each of the design parameters (sample size, sampling interval and control-limit width) vary as a function of the most recent process information. The cost function due to controlling the process quality through a VP (X) over bar chart is derived. During the optimization of the cost function, constraints are imposed on the expected times to signal when the process is in and out of control. In this way, required statistical properties can be assured. Through a numerical example, the proposed economic-statistical design approach for VP (X) over bar charts is compared to the economic design for VP (X) over bar charts and to the economic-statistical and economic designs for fixed parameters (FP) (X) over bar charts in terms of the operating cost and the expected times to signal. From this example, it is possible to assess the benefits provided by the proposed model. Varying some input parameters, their effect on the optimal cost and on the optimal values of the design parameters was analysed.