Turbulence, Complexity, and Solar Flares

The issue of predicting solar flares is one of the most fundamental in physics, addressing issues of plasma physics, high-energy physics, and modelling of complex systems. It also poses societal consequences, with our ever-increasing need for accurate space weather forecasts. Solar flares arise naturally as a competition between an input (flux emergence and rearrangement) in the photosphere and an output (electrical current build up and resistive dissipation) in the corona. Although initially localised, this redistribution affects neighbouring regions and an avalanche occurs resulting in large scale eruptions of plasma, particles, and magnetic field. As flares are powered from the stressed field rooted in the photosphere, a study of the photospheric magnetic complexity can be used to both predict activity and understand the physics of the magnetic field. The magnetic energy spectrum and multifractal spectrum are highlighted as two possible approaches to this.Comment: 2 figure

Comparison of mixed layer heights from airborne high spectral resolution lidar, ground-based measurements, and the WRF-Chem model during CalNex and CARES

The California Research at the Nexus of Air Quality and Climate Change (CalNex) and Carbonaceous Aerosol and Radiative Effects Study (CARES) field campaigns during May and June 2010 provided a data set appropriate for studying the structure of the atmospheric boundary layer (BL). The NASA Langley Research Center (LaRC) airborne high spectral resolution lidar (HSRL) was deployed to California onboard the NASA LaRC B-200 aircraft to aid in characterizing aerosol properties during these two field campaigns. Measurements of aerosol extinction (532 nm), backscatter (532 and 1064 nm), and depolarization (532 and 1064 nm) profiles during 31 flights, many in coordination with other research aircraft and ground sites, constitute a diverse data set for use in characterizing the spatial and temporal distribution of aerosols, as well as the depth and variability of the daytime mixed layer (ML) height. The paper describes the modified Haar wavelet covariance transform method used to derive the ML heights from HSRL backscatter profiles. HSRL ML heights are validated using ML heights derived from two radiosonde profile sites during CARES. Comparisons between ML heights from HSRL and a Vaisala ceilometer operated during CalNex were used to evaluate the representativeness of a fixed measurement over a larger region. In the Los Angeles basin, comparisons of ML heights derived from HSRL measurements and ML heights derived from the ceilometer result in a very good agreement (mean bias difference of 10 m and correlation coefficient of 0.89) up to 30 km away from the ceilometer site, but are essentially uncorrelated for larger distances, indicating that the spatial variability of the ML height is significant over these distances and not necessarily well captured by limited ground stations. The HSRL ML heights are also used to evaluate the performance in simulating the temporal and spatial variability of ML heights from the Weather Research and Forecasting Chemistry (WRF-Chem) community model. When compared to aerosol ML heights from HSRL, thermodynamic ML heights from WRF-Chem were underpredicted in the CalNex and CARES regions, shown by a bias difference value of −157 m and −29 m, respectively. Better agreement over the Central Valley than in mountainous regions suggests that some variability in the ML height is not well captured at the 4 km grid resolution of the model. A small but significant number of cases have poor agreement when WRF-Chem consistently overestimates the ML height in the late afternoon. Additional comparisons with WRF-Chem aerosol mixed layer heights show no significant improvement over thermodynamic ML heights, confirming that any differences between measurement and model are not due to the methodology of ML height determination

Feature extraction using extrema sampling of discrete derivatives for spike sorting in implantable upper-limb neural prostheses

Next generation neural interfaces for upper-limb (and other) prostheses aim to develop implantable interfaces for one or more nerves, each interface having many neural signal channels that work reliably in the stump without harming the nerves. To achieve real-time multi-channel processing it is important to integrate spike sorting on-chip to overcome limitations in transmission bandwidth. This requires computationally efficient algorithms for feature extraction and clustering suitable for low-power hardware implementation. This paper describes a new feature extraction method for real-time spike sorting based on extrema analysis (namely positive peaks and negative peaks) of spike shapes and their discrete derivatives at different frequency bands. Employing simulation across different datasets, the accuracy and computational complexity of the proposed method are assessed and compared with other methods. The average classification accuracy of the proposed method in conjunction with online sorting (O-Sort) is 91.6%, outperforming all the other methods tested with the O-Sort clustering algorithm. The proposed method offers a better tradeoff between classification error and computational complexity, making it a particularly strong choice for on-chip spike sorting