Automatic artefact removal from event-related potentials via clustering

This paper outlines a method for automatic artefact removal from multichannel recordings of event-related potentials (ERPs). The proposed method is based on, firstly, separation of the ERP recordings into independent components using the method of temporal decorrelation source separation (TDSEP). Secondly, the novel lagged auto-mutual information clustering (LAMIC) algorithm is used to cluster the estimated components, together with ocular reference signals, into clusters corresponding to cerebral and non-cerebral activity. Thirdly, the components in the cluster which contains the ocular reference signals are discarded. The remaining components are then recombined to reconstruct the clean ERPs.

Observing Forbush decreases in cloud at Shetland

Meteorological measurements from Lerwick Observatory, Shetland (60°09′N, 1°08′W), are compared with short-term changes in Climax neutron counter cosmic ray measurements. For transient neutron count reductions of 10–12%, broken cloud becomes at least 10% more frequent on the neutron minimum day, above expectations from sampling. This suggests a rapid timescale (1 day) cloud response to cosmic ray changes. However, larger or smaller neutron count reductions do not coincide with cloud responses exceeding sampling effects. Larger events are too rare to provide a robust signal above the sampling noise. Smaller events are too weak to be observed above the natural variability.

The global atmospheric electrical circuit and climate

Evidence is emerging for physical links among clouds, global temperatures, the global atmospheric electrical circuit and cosmic ray ionisation. The global circuit extends throughout the atmosphere from the planetary surface to the lower layers of the ionosphere. Cosmic rays are the principal source of atmospheric ions away from the continental boundary layer: the ions formed permit a vertical conduction current to flow in the fair weather part of the global circuit. Through the (inverse) solar modulation of cosmic rays, the resulting columnar ionisation changes may allow the global circuit to convey a solar influence to meteorological phenomena of the lower atmosphere. Electrical effects on non-thunderstorm clouds have been proposed to occur via the ion-assisted formation of ultra-fine aerosol, which can grow to sizes able to act as cloud condensation nuclei, or through the increased ice nucleation capability of charged aerosols. Even small atmospheric electrical modulations on the aerosol size distribution can affect cloud properties and modify the radiative balance of the atmosphere, through changes communicated globally by the atmospheric electrical circuit. Despite a long history of work in related areas of geophysics, the direct and inverse relationships between the global circuit and global climate remain largely quantitatively unexplored. From reviewing atmospheric electrical measurements made over two centuries and possible paleoclimate proxies, global atmospheric electrical circuit variability should be expected on many timescales

Experimental determination of layer cloud edge charging from cosmic ray ionisation

The cloud-air transition zone at stratiform cloud edges is an electrically active region where droplet charging has been predicted. Cloud edge droplet charging is expected from vertical flow of cosmic ray generated atmospheric ions in the global electric circuit. Experimental confirmation of stratiform cloud edge electrification is presented here, through charge and droplet measurements made within an extensive layer of supercooled stratiform cloud, using a specially designed electrostatic sensor. Negative space charge up to 35 pC m−3 was found in a thin (<100 m) layer at the lower cloud boundary associated with the clear air-cloud conductivity gradient, agreeing closely with space charge predicted from the measured droplet concentration using ion-aerosol theory. Such charge levels carried by droplets are sufficient to influence collision processes between cloud droplets.

Cloud optical depth retrievals from the Aerosol Robotic Network (AERONET) cloud mode observations

Cloud optical depth is one of the most poorly observed climate variables. The new “cloud mode” capability in the Aerosol Robotic Network (AERONET) will inexpensively yet dramatically increase cloud optical depth observations in both number and accuracy. Cloud mode optical depth retrievals from AERONET were evaluated at the Atmospheric Radiation Measurement program’s Oklahoma site in sky conditions ranging from broken clouds to overcast. For overcast cases, the 1.5 min average AERONET cloud mode optical depths agreed to within 15% of those from a standard ground‐based flux method. For broken cloud cases, AERONET retrievals also captured rapid variations detected by the microwave radiometer. For 3 year climatology derived from all nonprecipitating clouds, AERONET monthly mean cloud optical depths are generally larger than cloud radar retrievals because of the current cloud mode observation strategy that is biased toward measurements of optically thick clouds. This study has demonstrated a new way to enhance the existing AERONET infrastructure to observe cloud optical properties on a global scale.

An assessment of aerosol‐cloud interactions in marine stratus clouds based on surface remote sensing

An assessment of aerosol-cloud interactions (ACI) from ground-based remote sensing under coastal stratiform clouds is presented. The assessment utilizes a long-term, high temporal resolution data set from the Atmospheric Radiation Measurement (ARM) Program deployment at Pt. Reyes, California, United States, in 2005 to provide statistically robust measures of ACI and to characterize the variability of the measures based on variability in environmental conditions and observational approaches. The average ACIN (= dlnNd/dlna, the change in cloud drop number concentration with aerosol concentration) is 0.48, within a physically plausible range of 0–1.0. Values vary between 0.18 and 0.69 with dependence on (1) the assumption of constant cloud liquid water path (LWP), (2) the relative value of cloud LWP, (3) methods for retrieving Nd, (4) aerosol size distribution, (5) updraft velocity, and (6) the scale and resolution of observations. The sensitivity of the local, diurnally averaged radiative forcing to this variability in ACIN values, assuming an aerosol perturbation of 500 c-3 relative to a background concentration of 100 cm-3, ranges betwee-4 and -9 W -2. Further characterization of ACI and its variability is required to reduce uncertainties in global radiative forcing estimates.

Retrievals of thick cloud optical depth from the Geoscience Laser Altimeter System (GLAS) by calibration of solar background signal

Laser beams emitted from the Geoscience Laser Altimeter System (GLAS), as well as other spaceborne laser instruments, can only penetrate clouds to a limit of a few optical depths. As a result, only optical depths of thinner clouds (< about 3 for GLAS) are retrieved from the reflected lidar signal. This paper presents a comprehensive study of possible retrievals of optical depth of thick clouds using solar background light and treating GLAS as a solar radiometer. To do so one must first calibrate the reflected solar radiation received by the photon-counting detectors of the GLAS 532-nm channel, the primary channel for atmospheric products. Solar background radiation is regarded as a noise to be subtracted in the retrieval process of the lidar products. However, once calibrated, it becomes a signal that can be used in studying the properties of optically thick clouds. In this paper, three calibration methods are presented: (i) calibration with coincident airborne and GLAS observations, (ii) calibration with coincident Geostationary Opera- tional Environmental Satellite (GOES) and GLAS observations of deep convective clouds, and (iii) cali- bration from first principles using optical depth of thin water clouds over ocean retrieved by GLAS active remote sensing. Results from the three methods agree well with each other. Cloud optical depth (COD) is retrieved from the calibrated solar background signal using a one-channel retrieval. Comparison with COD retrieved from GOES during GLAS overpasses shows that the average difference between the two retriev- als is 24%. As an example, the COD values retrieved from GLAS solar background are illustrated for a marine stratocumulus cloud field that is too thick to be penetrated by the GLAS laser. Based on this study, optical depths for thick clouds will be provided as a supplementary product to the existing operational GLAS cloud products in future GLAS data releases.

Remote sensing of cloud properties using ground-based measurements of zenith radiance

We have conducted the first extensive field test of two new methods to retrieve optical properties for overhead clouds that range from patchy to overcast. The methods use measurements of zenith radiance at 673 and 870 nm wavelengths and require the presence of green vegetation in the surrounding area. The test was conducted at the Atmospheric Radiation Measurement Program Oklahoma site during September–November 2004. These methods work because at 673 nm (red) and 870 nm (near infrared (NIR)), clouds have nearly identical optical properties, while vegetated surfaces reflect quite differently. The first method, dubbed REDvsNIR, retrieves not only cloud optical depth τ but also radiative cloud fraction. Because of the 1-s time resolution of our radiance measurements, we are able for the first time to capture changes in cloud optical properties at the natural timescale of cloud evolution. We compared values of τ retrieved by REDvsNIR to those retrieved from downward shortwave fluxes and from microwave brightness temperatures. The flux method generally underestimates τ relative to the REDvsNIR method. Even for overcast but inhomogeneous clouds, differences between REDvsNIR and the flux method can be as large as 50%. In addition, REDvsNIR agreed to better than 15% with the microwave method for both overcast and broken clouds. The second method, dubbed COUPLED, retrieves τ by combining zenith radiances with fluxes. While extra information from fluxes was expected to improve retrievals, this is not always the case. In general, however, the COUPLED and REDvsNIR methods retrieve τ to within 15% of each other.

Interaction of nitrogen dioxide (NO2) with a monolayer of oleic acid at the air–water interface – A simple proxy for atmospheric aerosol

The reactions between atmospheric oxidants and organic amphiphiles at the air water interface of an aerosol droplet may affect the size and critical supersaturation required for cloud droplet formation. We demonstrate that no reaction occurs between gaseous nitrogen dioxide (1000 ppm in air) and a monolayer of an insoluble amphiphile, oleic acid (cis-9-octadecenoic acid), at the air water interface which removes material from the air water interface. We present evidence that the NO2 isomerises the cis-9-octadecenoic (oleic) acid to trans-9-octadecenoic (elaidic) acid. The study presented here is important for future and previous studies of (1) the reaction between the nitrate radical, NO3, and thin organic films as NO2 is usually present in high concentrations in these experimental systems and (2) the effect of NO2 air pollution on the unsaturated fatty acids and lipids found at the air liquid surface of human lung lining fluid.

Exploring NOy chemistry in levitated aqueous aerosol droplets

Chemistry of reactive nitrogen oxides, NOy, is crucial for our understanding of composition and properties of the Earth’s atmosphere. The proof-of-principle experiments demonstrated that we are able to study the atmospheric fate of nitrogen oxides that has significant impact on global climate and hydrological cycle, thus affecting the likelihood of local floods and acid rain.

Night-time tropospheric chemistry of the unsaturated alcohols (Z)-pent-2-en-1-ol and pent-1-en-3-ol: Kinetic studies of reactions of NO3 and N2O5 with stress-induced plant emissions

The night-time tropospheric chemistry of two stress-induced volatile organic compounds (VOCs), (Z)-pent-2-en-1-ol and pent-1-en-3-ol, has been studied at room temperature. Rate coefficients for reactions of the nitrate radical (NO3) with these pentenols were measured using the discharge-flow technique. Because of the relatively low volatility of these compounds, we employed off-axis continuous-wave cavity-enhanced absorption spectroscopy for detection of NO3 in order to be able to work in pseudo first-order conditions with the pentenols in large excess over NO3. The rate coefficients were determined to be (1.53 +/- 0.23) x 10(-13) and (1.39 +/- 0.19) x 10(-14) cm(3) molecule(-1) s(-1) for reactions of NO3 with (Z)-pent-2-en-1-ol and pent-1-en-3-ol. An attempt to study the kinetics of these reactions with a relative-rate technique, using N2O5 as source of NO3 resulted in significantly higher apparent rate coefficients. Performing relative-rate experiments in known excesses of NO2 allowed us to determine the rate coefficients for the N2O5 reactions to be (5.0 +/- 2.8) x 10(-19) cm(3) molecule(-1) s(-1) for (Z)-pent-2-en-1-ol, and (9.1 +/- 5.8) x 10(-19) cm(3) molecule(-1) s(-1) for pent-1-en-3-ol. We show that these relatively slow reactions can indeed interfere with rate determinations in conventional relative-rate experiments.

Influence of dynamics on the changes in tropical cloud radiative forcing during the 1998 El Niño

Satellite measurements of the radiation budget and data from the U.S. National Centers for Environmental Prediction–National Center for Atmospheric Research reanalysis are used to investigate the links between anomalous cloud radiative forcing over the tropical west Pacific warm pool and the tropical dynamics and sea surface temperature (SST) distribution during 1998. The ratio, N, of the shortwave cloud forcing (SWCF) to longwave cloud forcing (LWCF) (N = −SWCF/LWCF) is used to infer information on cloud altitude. A higher than average N during 1998 appears to be related to two separate phenomena. First, dynamic regime-dependent changes explain high values of N (associated with low cloud altitude) for small magnitudes of SWCF and LWCF (low cloud fraction), which reflect the unusual occurrence of mean subsiding motion over the tropical west Pacific during 1998, associated with the anomalous SST distribution. Second, Tropics-wide long-term changes in the spatial-mean cloud forcing, independent of dynamic regime, explain the higher values of N during both 1998 and in 1994/95. The changes in dynamic regime and their anomalous structure in 1998 are well simulated by version HadAM3 of the Hadley Centre climate model, forced by the observed SSTs. However, the LWCF and SWCF are poorly simulated, as are the interannual changes in N. It is argued that improved representation of LWCF and SWCF and their dependence on dynamical forcing are required before the cloud feedbacks simulated by climate models can be trusted.

Comments on “An ensemble cumulus convection parameterization with explicit cloud treatment”

Wagner and Graf (2010) derive a population evolution equation for an ensemble of convective plumes, an analogue with the Lotka–Volterra equation, from the energy equations for convective plumes provided by Arakawa and Schubert (1974). Although their proposal is interesting, as the present note shows, there are some problems with their derivation.

Cloud-resolving simulations of deep convection over a heated mountain

Cloud-resolving numerical simulations of airflow over a diurnally heated mountain ridge are conducted to explore the mechanisms and sensitivities of convective initiation under high pressure conditions. The simulations are based on a well-observed convection event from the Convective and Orographically Induced Precipitation Study (COPS) during summer 2007, where an isolated afternoon thunderstorm developed over the Black Forest mountains of central Europe, but they are idealized to facilitate understanding and reduce computational expense. In the conditionally unstable but strongly inhibited flow under consideration, sharp horizontal convergence over the mountain acts to locally weaken the inhibition and moisten the dry midtroposphere through shallow cumulus detrainment. The onset of deep convection occurs not through the deep ascent of a single updraft but rather through a rapid succession of thermals that are vented through the mountain convergence zone into the deepening cloud mass. Emerging thermals rise through the saturated wakes of their predecessors, which diminishes the suppressive effects of entrainment and allows for rapid glaciation above the freezing level as supercooled cloud drops rime onto preexisting ice particles. These effects strongly enhance the midlevel cloud buoyancy and enable rapid ascent to the tropopause. The existence and vigor of the convection is highly sensitive to small changes in background wind speed U0, which controls the strength of the mountain convergence and the ability of midlevel moisture to accumulate above the mountain. Whereas vigorous deep convection develops for U0 = 0 m s−1, deep convection is completely eliminated for U0 = 3 m s−1. Although deep convection is able to develop under intermediate winds (U0 = 1.5 m s−1), its formation is highly sensitive to small-amplitude perturbations in the initial flow.

Fair weather atmospheric electricity

Abstract. Not long after Franklin’s iconic studies, an atmospheric electric field was discovered in “fair weather” regions, well away from thunderstorms. The origin of the fair weather field was sought by Lord Kelvin, through development of electrostatic instrumentation and early data logging techniques, but was ultimately explained through the global circuit model of C.T.R. Wilson. In Wilson’s model, charge exchanged by disturbed weather electrifies the ionosphere, and returns via a small vertical current density in fair weather regions. New insights into the relevance of fair weather atmospheric electricity to terrestrial and planetary atmospheres are now emerging. For example, there is a possible role of the global circuit current density in atmospheric processes, such as cloud formation. Beyond natural atmospheric processes, a novel practical application is the use of early atmospheric electrostatic investigations to provide quantitative information on past urban air pollution.