ALBERT

Description: A Lagrangian advection scheme with shape matrix (LASM) for solving advection problems Geoscientific Model Development, 7, 2951-2968, 2014 Author(s): L. Dong, B. Wang, and L. Liu A new Lagrangian advection scheme with shape matrix (LASM) is proposed to take advantage of the extreme low numerical diffusion of the Lagrangian methods. The tracer is discretized into finite parcels, which move along the downstream trajectories. Different from other Lagrangian schemes, the parcel shape is simulated explicitly by a linear transformation matrix. By doing so, the aliasing error in the Lagrangian schemes is largely reduced without introducing substantial interparcel mixing in the pure advection stage, because the flow information will be respected when remapping tracer density onto the fixed model grids. An adaptive interparcel mixing algorithm is constructed to ensure the validity of the linear approximation of the parcel shape, where the mixing is only triggered when it is necessary and resembles the physical mixing. The total tracer mass on the parcels is conserved exactly. The new scheme is validated by using several test cases.

Description: Technical Note: Reducing the spin-up time of integrated surface water–groundwater models Hydrology and Earth System Sciences, 18, 5169-5179, 2014 Author(s): H. Ajami, J. P. Evans, M. F. McCabe, and S. Stisen One of the main challenges in the application of coupled or integrated hydrologic models is specifying a catchment's initial conditions in terms of soil moisture and depth-to-water table (DTWT) distributions. One approach to reducing uncertainty in model initialization is to run the model recursively using either a single year or multiple years of forcing data until the system equilibrates with respect to state and diagnostic variables. However, such "spin-up" approaches often require many years of simulations, making them computationally intensive. In this study, a new hybrid approach was developed to reduce the computational burden of the spin-up procedure by using a combination of model simulations and an empirical DTWT function. The methodology is examined across two distinct catchments located in a temperate region of Denmark and a semi-arid region of Australia. Our results illustrate that the hybrid approach reduced the spin-up period required for an integrated groundwater–surface water–land surface model (ParFlow.CLM) by up to 50%. To generalize results to different climate and catchment conditions, we outline a methodology that is applicable to other coupled or integrated modeling frameworks when initialization from an equilibrium state is required.

Description: Extending periodic eddy covariance latent heat fluxes through tree sapflow measurements to estimate long-term total evaporation in a peat swamp forest Hydrology and Earth System Sciences Discussions, 11, 13607-13661, 2014 Author(s): A. D. Clulow, C. S. Everson, M. G. Mengistu, J. S. Price, A. Nickless, and G. P. W. Jewitt A combination of measurement and modelling was used to find a pragmatic solution to estimate the annual total evaporation (ET) from the rare and indigenous Nkazana Peat Swamp Forest (PSF) on the east coast of Southern Africa to improve the water balance estimates within the area. Total evaporation was measured during three window periods (between seven and nine days each) using an eddy covariance (EC) system on a telescopic mast above the forest canopy. Sapflow of an understory and an emergent tree was measured using a low maintenance heat pulse velocity system for an entire hydrological year (October 2009 to September 2010). An empirical model was derived, describing the relationship between the observed ET of the Nkazana PSF measured during two of the window periods ( R 2 = 0.92 and 0.90) which, overlapped with sapflow measurements, thereby providing hourly estimates of predicted ET of the Nkazana PSF for a year, totalling 1125 mm (while rainfall was 650 mm). In building the empirical model, it was found that including the understory tree sapflow provided no benefit to the model performance. In addition, the observed emergent tree sapflow relationship with observed ET between the two field campaigns was consistent and could be represented by a single empirical model ( R 2 = 0.90; RMSE = 0.08 mm). During the window periods of EC measurement, no single meteorological variable was found to describe the Nkazana PSF ET satisfactorily. However, in terms of evaporation models, the hourly FAO56 Penman–Monteith equation best described the observed ET from EC during the August 2009 ( R 2 = 0.75), November 2009 ( R 2 = 0.85) and March 2010 ( R 2 = 0.76) field campaigns, compared to the Priestley–Taylor model ( R 2 = 0.54, 0.74 and 0.62 during the respective field campaigns). From the empirical model of ET and the FAO56 Penman–Monteith equation, a monthly crop factor ( K c ) was derived for the Nkazana PSF providing a method of estimating long-term swamp forest ET from meteorological data. The monthly crop factor indicated two distinct periods. From February to May, it was between 1.2 and 1.4 compared with June to January, when the crop factor was 0.8 to 1.0. The derived monthly K c values were verified as accurate (to one significant digit) using historical data measured at the same site, also using EC, from a~previous study. The measurements provided insights into the microclimate within a subtropical peat swamp forest and the contrasting sapflow of emergent and understory trees. They showed that expensive, high maintenance equipment can be used during manageable window periods in conjunction with low maintenance systems, dedicated to individual trees, to derive a model to estimate long-term ET over remote heterogeneous forests. In addition, the contrast in ET and rainfall emphasises the reliance of the Nkazana PSF on groundwater.

Description: Bio-optical variables in the Northern Gulf of California were analyzed using in situ and satellite data obtained during a cruise in June 2008. The study area was divided into three bio-optical regions: Upper Gulf (UG), Northern Gulf (NG), and Great Isles (GI). Each region was characterized according to phytoplankton pigment concentration, phytoplankton and nonpigmented material spectral absorption coefficients, and spectral reflectance. Observed patterns were an indication of the shift in bio-optical conditions from north to south going from turbid and eutrophic waters to mesotrophic ones. Although there was a good agreement between satellite and in situ Chla (RMSE ±33%), an overestimation of in situ Chla was observed. This was partly explained by the presence of nonalgal particles, as well as the influence of desert and continental aerosols, which is generally overcorrected in the standard processing. The UG and NG could be considered as Case  2 waters, but they did exhibit different bio-optical characteristics. This implies that both biological and optical properties should be invoked to better understand water reflectance variability in the study region and its implications for the remote sensing of Chla and biogeochemical processes.

Description: Consensus of continuous-time single-integrator multiagent systems with inherent nonlinear dynamics and measurement noises is addressed in this paper. The consensus controller is developed for directed interaction topologies. Each agent’s control input is constructed based on its own state and its neighbors’ states corrupted by noises. The control input contains a time-varying consensus gain in order to attenuate the noises. Conditions for ensuring mean square convergence under noisy measurement and asymptotic convergence in the noise-free case are derived. Finally, some simulations were carried out to show the effectiveness of our control law and how it can solve the consensus problem.

Description: Stratospheric ozone depletion from future nitrous oxide increases Atmospheric Chemistry and Physics, 14, 12967-12982, 2014 Author(s): W. Wang, W. Tian, S. Dhomse, F. Xie, J. Shu, and J. Austin We have investigated the impact of the assumed nitrous oxide (N 2 O) increases on stratospheric chemistry and dynamics using a series of idealized simulations with a coupled chemistry-climate model (CCM). In a future cooler stratosphere the net yield of NO y from N 2 O is shown to decrease in a reference run following the IPCC A1B scenario, but NO y can still be significantly increased by extra increases of N 2 O over 2001–2050. Over the last decade of simulations, 50% increases in N 2 O result in a maximal 6% reduction in ozone mixing ratios in the middle stratosphere at around 10 hPa and an average 2% decrease in the total ozone column (TCO) compared with the control run. This enhanced destruction could cause an ozone decline in the first half of this century in the middle stratosphere around 10 hPa, while global TCO still shows an increase at the same time. The results from a multiple linear regression analysis and sensitivity simulations with different forcings show that the chemical effect of N 2 O increases dominates the N 2 O-induced ozone depletion in the stratosphere, while the dynamical and radiative effects of N 2 O increases are overall insignificant. The analysis of the results reveals that the ozone depleting potential of N 2 O varies with the time period and is influenced by the environmental conditions. For example, carbon dioxide (CO 2 ) increases can strongly offset the ozone depletion effect of N 2 O.

Description: Corrigendum to "Air-sea exchange and gas-particle partitioning of polycyclic aromatic hydrocarbons in the Mediterranean" published in Atmos. Chem. Phys., 14, 8905–8915, 2014 Atmospheric Chemistry and Physics, 14, 12965-12965, 2014 Author(s): M. D. Mulder, A. Heil, P. Kukučka, J. Klánová, J. Kuta, R. Prokeš, F. Sprovieri, and G. Lammel No abstract available.

Description: Mechanistic study of secondary organic aerosol components formed from nucleophilic addition reactions of methacrylic acid epoxide Atmospheric Chemistry and Physics, 14, 12951-12964, 2014 Author(s): A. W. Birdsall, C. R. Miner, L. E. Mael, and M. J. Elrod Recently, methacrylic acid epoxide (MAE) has been proposed as a precursor to an important class of isoprene-derived compounds found in secondary organic aerosol (SOA): 2-methylglyceric acid (2-MG) and a set of oligomers, nitric acid esters, and sulfuric acid esters related to 2-MG. However, the specific chemical mechanisms by which MAE could form these compounds have not been previously studied with experimental methods. In order to determine the relevance of these processes to atmospheric aerosol, MAE and 2-MG have been synthesized and a series of bulk solution-phase experiments aimed at studying the reactivity of MAE using nuclear magnetic resonance (NMR) spectroscopy have been performed. The present results indicate that the acid-catalyzed MAE reaction is more than 600 times slower than a similar reaction of an important isoprene-derived epoxide, but is still expected to be kinetically feasible in the atmosphere on more acidic SOA. The specific mechanism by which MAE leads to oligomers was identified, and the reactions of MAE with a number of atmospherically relevant nucleophiles were also investigated. Because the nucleophilic strengths of water, sulfate, alcohols (including 2-MG), and acids (including MAE and 2-MG) in their reactions with MAE were found to be of similar magnitudes, it is expected that a diverse variety of MAE + nucleophile product species may be formed on ambient SOA. Thus, the results indicate that epoxide chain reaction oligomerization will be limited by the presence of high concentrations of non-epoxide nucleophiles (such as water); this finding is consistent with previous environmental chamber investigations of the relative humidity dependence of 2-MG-derived oligomerization processes and suggests that extensive oligomerization may not be likely on ambient SOA because of other competitive MAE reaction mechanisms.

Description: Megacity emission plume characteristics in summer and winter investigated by mobile aerosol and trace gas measurements: the Paris metropolitan area Atmospheric Chemistry and Physics, 14, 12931-12950, 2014 Author(s): S.-L. von der Weiden-Reinmüller, F. Drewnick, Q. J. Zhang, F. Freutel, M. Beekmann, and S. Borrmann For the investigation of megacity emission plume characteristics mobile aerosol and trace gas measurements were carried out in the greater Paris region in July 2009 and January–February 2010 within the EU FP7 MEGAPOLI project (Megacities: Emissions, urban, regional and Global Atmospheric POLlution and climate effects, and Integrated tools for assessment and mitigation). The deployed instruments measured physical and chemical properties of sub-micron aerosol particles, gas phase constituents of relevance for urban air pollution studies and meteorological parameters. The emission plume was identified based on fresh pollutant (e.g., particle-bound polycyclic aromatic hydrocarbons, black carbon, CO 2 and NO x ) concentration changes in combination with wind direction data. The classification into megacity influenced and background air masses allowed a characterization of the emission plume during summer and winter environmental conditions. On average, a clear increase of fresh pollutant concentrations in plume compared to background air masses was found for both seasons. For example, an average increase of 190% (+ 8.8 ng m −3 ) in summer and of 130% (+ 18.1 ng m −3 ) in winter was found for particle-bound polycyclic aromatic hydrocarbons in plume air masses. The aerosol particle size distribution in plume air masses was influenced by nucleation and growth due to coagulation and condensation in summer, while in winter only the latter process (i.e., particle growth) seemed to be initiated by urban pollution. The observed distribution of fresh pollutants in the emission plume – its cross sectional Gaussian-like profile and the exponential decrease of pollutant concentrations with increasing distance to the megacity – are in agreement with model results. Differences between model and measurements were found for plume center location, plume width and axial plume extent. In general, dilution was identified as the dominant process determining the axial variations within the Paris emission plume. For in-depth analysis of transformation processes occurring in the advected plume, simultaneous measurements at a suburban measurement site and a stationary site outside the metropolitan area using the mobile laboratory have proven to be most useful. Organic aerosol oxidation was observed in summer, while in winter transformation processes seemed to occur at a slower rate.

Description: A regional CO 2 observing system simulation experiment for the ASCENDS satellite mission Atmospheric Chemistry and Physics, 14, 12897-12914, 2014 Author(s): J. S. Wang, S. R. Kawa, J. Eluszkiewicz, D. F. Baker, M. Mountain, J. Henderson, T. Nehrkorn, and T. S. Zaccheo Top–down estimates of the spatiotemporal variations in emissions and uptake of CO 2 will benefit from the increasing measurement density brought by recent and future additions to the suite of in situ and remote CO 2 measurement platforms. In particular, the planned NASA Active Sensing of CO 2 Emissions over Nights, Days, and Seasons (ASCENDS) satellite mission will provide greater coverage in cloudy regions, at high latitudes, and at night than passive satellite systems, as well as high precision and accuracy. In a novel approach to quantifying the ability of satellite column measurements to constrain CO 2 fluxes, we use a portable library of footprints (surface influence functions) generated by the Stochastic Time-Inverted Lagrangian Transport (STILT) model in combination with the Weather Research and Forecasting (WRF) model in a regional Bayesian synthesis inversion. The regional Lagrangian particle dispersion model framework is well suited to make use of ASCENDS observations to constrain weekly fluxes in North America at a high resolution, in this case at 1° latitude × 1° longitude. We consider random measurement errors only, modeled as a function of the mission and instrument design specifications along with realistic atmospheric and surface conditions. We find that the ASCENDS observations could potentially reduce flux uncertainties substantially at biome and finer scales. At the grid scale and weekly resolution, the largest uncertainty reductions, on the order of 50%, occur where and when there is good coverage by observations with low measurement errors and the a priori uncertainties are large. Uncertainty reductions are smaller for a 1.57 μm candidate wavelength than for a 2.05 μm wavelength, and are smaller for the higher of the two measurement error levels that we consider (1.0 ppm vs. 0.5 ppm clear-sky error at Railroad Valley, Nevada). Uncertainty reductions at the annual biome scale range from ~40% to ~75% across our four instrument design cases and from ~65% to ~85% for the continent as a whole. Tests suggest that the quantitative results are moderately sensitive to assumptions regarding a priori uncertainties and boundary conditions. The a posteriori flux uncertainties we obtain, ranging from 0.01 to 0.06 Pg C yr −1 across the biomes, would meet requirements for improved understanding of long-term carbon sinks suggested by a previous study.