The Lagrangian particle dispersion model FLEXPART was originally designed for calculating long-range and mesoscale dispersion of air pollutants from point sources, such that occurring after an ...accident in a nuclear power plant. In the meantime, FLEXPART has evolved into a comprehensive tool for atmospheric transport modeling and analysis at different scales. A need for further multiscale modeling and analysis has encouraged new developments in FLEXPART. In this paper, we present a FLEXPART version that works with the Weather Research and Forecasting (WRF) mesoscale meteorological model. We explain how to run this new model and present special options and features that differ from those of the preceding versions. For instance, a novel turbulence scheme for the convective boundary layer has been included that considers both the skewness of turbulence in the vertical velocity as well as the vertical gradient in the air density. To our knowledge, FLEXPART is the first model for which such a scheme has been developed. On a more technical level, FLEXPART-WRF now offers effective parallelization, and details on computational performance are presented here. FLEXPART-WRF output can either be in binary or Network Common Data Form (NetCDF) format, both of which have efficient data compression. In addition, test case data and the source code are provided to the reader as a Supplement. This material and future developments will be accessible at http://www.flexpart.eu.
The CARES campaign was conducted during June, 2010 in the vicinity of Sacramento, California to study aerosol formation and aging in a region where anthropogenic and biogenic emissions regularly mix. ...Here, we describe measurements from an Aerodyne High Resolution Aerosol Mass Spectrometer (AMS), an Ionicon Proton Transfer Reaction Mass Spectrometer (PTR-MS), and trace gas detectors (CO, NO, NO sub(x)) deployed on the G-1 research aircraft to investigate ambient gas- and particle-phase chemical composition. AMS measurements showed that the particle phase is dominated by organic aerosol (OA) (85% on average) with smaller concentrations of sulfate (5%), nitrate (6%) and ammonium (3%) observed. PTR-MS data showed that isoprene dominated the biogenic volatile organic compound concentrations (BVOCs), with monoterpene concentrations generally below the detection limit. Using two different metrics, median OA concentrations and the slope of plots of OA vs. CO concentrations (i.e., Delta OA/ Delta CO), we contrast organic aerosol evolution on flight days with different prevailing meteorological conditions to elucidate the role of anthropogenic and biogenic emissions on OA formation. Airmasses influenced predominantly by biogenic emissions had median OA concentrations of 2.2 mu g m super(-3) and near zero Delta OA/ Delta CO. Those influenced predominantly by anthropogenic emissions had median OA concentrations of 4.7 mu g m super(-3) and Delta OA/ Delta CO ratios of 35-44 mu g m super(-3) ppmv. But, when biogenic and anthropogenic emissions mixed, OA levels were enhanced, with median OA concentrations of 11.4 mu g m super(-3) and Delta OA/ Delta CO ratios of 77-157 mu g m super(-3) ppmv. Taken together, our observations show that production of OA was enhanced when anthropogenic emissions from Sacramento mixed with isoprene-rich air from the foothills. After considering several anthropogenic/biogenic interaction mechanisms, we conclude that NO sub(x) concentrations play a strong role in enhancing SOA formation from isoprene, though the chemical mechanism for the enhancement remains unclear. If these observations are found to be robust in other seasons and in areas outside of Sacramento, regional and global aerosol modules will need to incorporate more complex representations of NO sub(x)-dependent SOA mechanisms and yields into their algorithms. Ultimately, accurately predicting OA mass concentrations and their effect on radiation balance will require a mechanistically-based treatment of the interactions of biogenic and anthropogenic emissions.
We report on several features in the energy spectrum from an ultralow-noise germanium detector operated deep underground. By implementing a new technique able to reject surface events, a number of ...cosmogenic peaks can be observed for the first time. We discuss an irreducible excess of bulklike events below 3 keV in ionization energy. These could be caused by unknown backgrounds, but also dark matter interactions consistent with DAMA/LIBRA. It is not yet possible to determine their origin. Improved constraints are placed on a cosmological origin for the DAMA/LIBRA effect.
Fifteen months of cumulative CoGeNT data are examined for indications of an annual modulation, a predicted signature of weakly interacting massive particle (WIMP) interactions. Presently available ...data support the presence of a modulated component of unknown origin, with parameters prima facie compatible with a galactic halo composed of light-mass WIMPs. Unoptimized estimators yield a statistical significance for a modulation of ∼2.8σ, limited by the short exposure.
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.
An Aerodyne high resolution time-of-flight aerosol mass spectrometer (HR-ToF-AMS) was deployed during the Carbonaceous Aerosols and Radiative Effects Study (CARES) that took place in northern ...California in June 2010. We present results obtained at Cool (denoted as the T1 site of the project) in the foothills of the Sierra Nevada Mountains, where intense biogenic emissions are periodically mixed with urban outflow transported by daytime southwesterly winds from the Sacramento metropolitan area. During this study, the average mass loading of submicrometer particles (PM1) was 3.0 μg m−3, dominated by organics (80%) and sulfate (9.9%). The organic aerosol (OA) had a nominal formula of C1H1.38N0.004OM0.44, thus an average organic mass-to-carbon (OM/OC) ratio of 1.70. Two distinct oxygenated OA factors were identified via Positive matrix factorization (PMF) of the high-resolution mass spectra of organics. The more oxidized MO-OOA (O/C = 0.54) was interpreted as a surrogate for secondary OA (SOA) influenced by biogenic emissions whereas the less oxidized LO-OOA (O/C = 0.42) was found to represent SOA formed in photochemically processed urban emissions. LO-OOA correlated strongly with ozone and MO-OOA correlated well with two 1st generation isoprene oxidation products (methacrolein and methyl vinyl ketone), indicating that both SOAs were relatively fresh. A hydrocarbon like OA (HOA) factor was also identified, representing primary emissions mainly due to local traffic. On average, SOA (= MO-OOA + LO-OOA) accounted for 91% of the total OA mass and 72% of the PM1 mass observed at Cool. Twenty three periods of urban plumes from T0 (Sacramento) to T1 (Cool) were identified using the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem). The average PM1 mass loading was considerably higher in urban plumes than in air masses dominated by biogenic SOA. The change in OA mass relative to CO (ΔOA/ΔCO) varied in the range of 5-196 μg m−3 ppm−1, reflecting large variability in SOA production. The highest ΔOA/ΔCO was reached when air masses were dominated by anthropogenic emissions in the presence of a high concentration of biogenic volatile organic compounds (BVOCs). This ratio, which was 97 μg m−3 ppm−1 on average, was much higher than when urban plumes arrived in a low BVOC environment (~36 μg m−3 ppm−1) or during other periods dominated by biogenic SOA (35 μg m−3 ppm−1). These results demonstrate that SOA formation is enhanced when anthropogenic emissions interact with biogenic precursors.
A claim for evidence of dark matter interactions in the DAMA experiment has been recently reinforced. We employ a new type of germanium detector to conclusively rule out a standard isothermal ...galactic halo of weakly interacting massive particles as the explanation for the annual modulation effect leading to the claim. Bounds are similarly imposed on a suggestion that dark pseudoscalars might lead to the effect. We describe the sensitivity to light dark matter particles achievable with our device, in particular, to next-to-minimal supersymmetric model candidates.
The performance of the Weather Research and Forecasting regional model with chemistry (WRF-Chem) in simulating the spatial and temporal variations in aerosol mass, composition, and size over ...California is quantified using the extensive meteorological, trace gas, and aerosol measurements collected during the California Nexus of Air Quality and Climate Experiment (CalNex) and the Carbonaceous Aerosol and Radiative Effects Study (CARES) conducted during May and June of 2010. The overall objective of the field campaigns was to obtain data needed to better understand processes that affect both climate and air quality, including emission assessments, transport and chemical aging of aerosols, aerosol radiative effects. Simulations were performed that examined the sensitivity of aerosol concentrations to anthropogenic emissions and to long-range transport of aerosols into the domain obtained from a global model. The configuration of WRF-Chem used in this study is shown to reproduce the overall synoptic conditions, thermally driven circulations, and boundary layer structure observed in region that controls the transport and mixing of trace gases and aerosols. Reducing the default emissions inventory by 50% led to an overall improvement in many simulated trace gases and black carbon aerosol at most sites and along most aircraft flight paths; however, simulated organic aerosol was closer to observed when there were no adjustments to the primary organic aerosol emissions. We found that sulfate was better simulated over northern California whereas nitrate was better simulated over southern California. While the overall spatial and temporal variability of aerosols and their precursors were simulated reasonably well, we show cases where the local transport of some aerosol plumes were either too slow or too fast, which adversely affects the statistics quantifying the differences between observed and simulated quantities. Comparisons with lidar and in situ measurements indicate that long-range transport of aerosols from the global model was likely too high in the free troposphere even though their concentrations were relatively low. This bias led to an over-prediction in aerosol optical depth by as much as a factor of 2 that offset the under-predictions of boundary-layer extinction resulting primarily from local emissions. Lowering the boundary conditions of aerosol concentrations by 50% greatly reduced the bias in simulated aerosol optical depth for all regions of California. This study shows that quantifying regional-scale variations in aerosol radiative forcing and determining the relative role of emissions from local and distant sources is challenging during `clean' conditions and that a wide array of measurements are needed to ensure model predictions are correct for the right reasons. In this regard, the combined CalNex and CARES data sets are an ideal test bed that can be used to evaluate aerosol models in great detail and develop improved treatments for aerosol processes.
A comparison between observed aerosol optical properties from the MILAGRO field campaign, which took place in the Mexico City Metropolitan Area (MCMA) during March 2006, and values simulated by the ...Weather Research and Forecasting (WRF-Chem) model, reveals large differences. To help identify the source of the discrepancies, data from the MILAGRO campaign are used to evaluate the "aerosol chemical to aerosol optical properties" module implemented in the full chemistry version of the WRF-Chem model. The evaluation uses measurements of aerosol size distributions and chemical properties obtained at the MILAGRO T1 site. These observations are fed to the module, which makes predictions of various aerosol optical properties, including the scattering coefficient, Bscat; the absorption coefficient, Babs; and the single-scattering albedo, ϖ0; all as a function of time. Values simulated by the module are compared with independent measurements obtained from a photoacoustic spectrometer (PAS) at a wavelength of 870 nm. Because of line losses and other factors, only "fine mode" aerosols with aerodynamic diameters less than 2.5 μm are considered here. Over a 10-day period, the simulations of hour-by-hour variations of Bscat are not satisfactory, but simulations of Babs and ϖ0 are considerably better. When averaged over the 10-day period, the computed and observed optical properties agree within the uncertainty limits of the measurements and simulations. Specifically, the observed and calculated values are, respectively: (1) Bscat, 34.1±5.1 Mm−1 versus 30.4±3.4 Mm−1; (2) Babs, 9.7±1.0 Mm−1 versus 11.7±1.2 Mm−1; and (3) ϖ0, 0.78±0.05 and 0.74±0.03. The discrepancies in values of ϖ0 simulated by the full WRF-Chem model thus cannot be attributed to the "aerosol chemistry to optics" module. The discrepancy is more likely due, in part, to poor characterization of emissions near the T1 site, particularly black carbon emissions.