ALBERT

Description: Predicting future thaw slump activity requires a sound understanding of the atmospheric drivers and geomorphic controls on mass wasting across a range of time scales. On sub-seasonal time scales, sparse measurements indicate that mass wasting at active slumps is often limited by the energy available for melting ground ice, but other factors such as rainfall or the formation of an insulating veneer may also be relevant. To study the sub-seasonal drivers, we derive topographic changes from single-pass radar interferometric data acquired by the TanDEM-X satellite (12 m resolution). The high vertical precision (around 30 cm), frequent observations (11 days) and large coverage (5000 km2) allow us to track volume losses as drivers such as the available energy change during summer in two study regions. We find that thaw slumps in the Tuktoyaktuk coastlands, Canada, are not energy limited in June, as they undergo limited mass wasting (height loss of around 0 cm/day) despite the ample available energy, indicating the widespread presence of an insulating snow or debris veneer. Later in summer, height losses generally increase (around 3 cm/day), but they do so in distinct ways. For many slumps, mass wasting tracks the available energy, a temporal pattern that is also observed at coastal yedoma cliffs on the Bykovsky Peninsula, Russia. However, the other two common temporal trajectories are asynchronous with the available energy, as they track strong precipitation events or show a sudden speed-up in late August, respectively. The observed temporal patterns are poorly related to slump characteristics like the slump area. The contrasting temporal behaviour of nearby thaw slumps highlights the importance of complex local and temporally varying controls on mass wasting.

Description: A suite of oxygenated volatile organic compounds (OVOCs – acetaldehyde, acetone, propanal, butanal and butanone) were measured concurrently in the surface water and atmosphere of the South China Sea and Sulu Sea in November 2011. A strong correlation was observed between all OVOC concentrations in the surface seawater along the entire cruise track, except for acetaldehyde, suggesting similar sources and sinks in the surface ocean. Additionally, several phytoplankton groups, such as haptophytes or pelagophytes, were also correlated to all OVOCs indicating that phytoplankton may be an important source for marine OVOCs in the South China and Sulu Seas. Humic and protein like fluorescent dissolved organic matter (FDOM) components seemed to be additional precursors for butanone and acetaldehyde. The atmospheric OVOC mixing ratios were relative high compared with literature values, suggesting the coastal region of North Borneo as a local hot spot for atmospheric OVOCs. The flux of atmospheric OVOCs was largely into the ocean for all 5 gases, with a few important exceptions near the coast of Borneo. The calculated amount of OVOCs entrained into the ocean seemed to be an important source of OVOCs to the surface ocean. When the fluxes were out of the ocean, marine OVOCs were found to be enough to control the local measured OVOC distribution in the atmosphere. Based on our model calculations, at least 0.4 ppb of marine derived acetone and butanone can reach the upper troposphere, where they may have an important influence on hydrogen oxide radical formation over the western Pacific Ocean.

Description: The Extrapolar SWIFT model is a fast ozone chemistry scheme for interactive calculation of the extrapolar stratospheric ozone layer in coupled general circulation models (GCMs). In contrast to the widely used prescribed ozone, the SWIFT ozone layer interacts with the model dynamics and can respond to atmospheric variability or climatological trends. The Extrapolar SWIFT model employs a repro-modelling approach, where algebraic functions are used to approximate the numerical output of a full stratospheric chemistry and transport model (ATLAS). The full model solves a coupled chemical differential equations system with 55 initial and boundary conditions (mixing ratio of various chemical species and atmospheric parameters). Hence the rate of change of ozone over 24  h is a function of 55 variables. Using covariances between these variables, we can find linear combinations in order to reduce the parameter space to the following nine basic variables: latitude, pressure altitude, temperature, local ozone column, mixing ratio of ozone and of the ozone depleting families (Cly, Bry, NOy and HOy). We will show that these 9 variables are sufficient to characterize the rate of change of ozone. An automated procedure fits a polynomial function of fourth degree to the rate of change of ozone obtained from several simulations with the ATLAS model. One polynomial function is determined per month which yields the rate of change of ozone over 24 h. A key aspect for the robustness of the Extrapolar SWIFT model is to include a wide range of stratospheric variability in the numerical output of the ATLAS model, also covering atmospheric states that will occur in a future climate (e.g. temperature and meridional circulation changes or reduction of stratospheric chlorine loading). For validation purposes, the Extrapolar SWIFT model has been integrated into the ATLAS model replacing the full stratospheric chemistry scheme. Simulations with SWIFT in ATLAS have proven that the systematic error is small and does not accumulate during the course of a simulation. In the context of a 10 year simulation, the ozone layer, simulated by SWIFT, shows a stable annual cycle, with inter-annual variations comparable to the ATLAS model. The application of Extrapolar SWIFT requires the evaluation of polynomial functions with 30–100 terms. Nowadays, computers can calculate such polynomial functions at thousands of model grid points in seconds. SWIFT provides the desired numerical efficiency and computes the ozone layer 104 times faster than the chemistry scheme in the ATLAS CTM.

Description: Ice retreat in the eastern Eurasian Arctic is a consequence of atmospheric and oceanic processes and regional feedback mechanisms acting on the ice cover, both in winter and summer. A correct representation of these processes in numerical models is important, since it will improve predictions of sea ice anomalies along the Northeast Passage and beyond. In this study, we highlight the importance of winter ice dynamics for local summer sea ice anomalies in thickness, volume and extent. By means of airborne sea ice thickness surveys made over pack ice areas in the south-eastern Laptev Sea, we show that years of offshore-directed sea ice transport have a thinning effect on the late-winter sea ice cover. To confirm the preconditioning effect of enhanced offshore advection in late winter on the summer sea ice cover, we perform a sensitivity study using a numerical model. Results verify that the preconditioning effect plays a bigger role for the regional ice extent. Furthermore, they indicate an increase in volume export from the Laptev Sea as a consequence of enhanced offshore advection, which has far-reaching consequences for the entire Arctic sea ice mass balance. Moreover we show that ice dynamics in winter not only preconditions local summer ice extent, but also accelerate fast-ice decay.

Description: Climate trends in the Antarctic region remain poorly characterized, owing to the brevity and scarcity of direct climate observations and the large magnitude of interannual to decadal-scale climate variability. Here, within the framework of the PAGES Antarctica2k working group, we build an enlarged database of ice core water stable isotope records from Antarctica, consisting of 112 records. We produce both unweighted and weighted isotopic (δ18O) composites and temperature reconstructions since 0 CE, binned at 5- and 10-year resolution, for seven climatically distinct regions covering the Antarctic continent. Following earlier work of the Antarctica2k working group, we also produce composites and reconstructions for the broader regions of East Antarctica, West Antarctica and the whole continent. We use three methods for our temperature reconstructions: (i) a temperature scaling based on the δ18O–temperature relationship output from an ECHAM5-wiso model simulation nudged to ERA-Interim atmospheric reanalyses from 1979 to 2013, and adjusted for the West Antarctic Ice Sheet region to borehole temperature data, (ii) a temperature scaling of the isotopic normalized anomalies to the variance of the regional reanalysis temperature and (iii) a composite-plus-scaling approach used in a previous continent-scale reconstruction of Antarctic temperature since 1 CE but applied to the new Antarctic ice core database. Our new reconstructions confirm a significant cooling trend from 0 to 1900 CE across all Antarctic regions where records extend back into the 1st millennium, with the exception of the Wilkes Land coast and Weddell Sea coast regions. Within this long-term cooling trend from 0 to 1900 CE, we find that the warmest period occurs between 300 and 1000 CE, and the coldest interval occurs from 1200 to 1900 CE. Since 1900 CE, significant warming trends are identified for the West Antarctic Ice Sheet, the Dronning Maud Land coast and the Antarctic Peninsula regions, and these trends are robust across the distribution of records that contribute to the unweighted isotopic composites and also significant in the weighted temperature reconstructions. Only for the Antarctic Peninsula is this most recent century-scale trend unusual in the context of natural variability over the last 2000 years. However, projected warming of the Antarctic continent during the 21st century may soon see significant and unusual warming develop across other parts of the Antarctic continent. The extended Antarctica2k ice core isotope database developed by this working group opens up many avenues for developing a deeper understanding of the response of Antarctic climate to natural and anthropogenic climate forcings. The first long-term quantification of regional climate in Antarctica presented herein is a basis for data–model comparison and assessments of past, present and future driving factors of Antarctic climate.

Description: Whereas ice cores from high-accumulation sites in coastal Antarctica clearly demonstrate annual layering, it is debated whether a seasonal signal is also preserved in ice cores from lower-accumulation sites further inland and particularly on the East Antarctic Plateau. In this study, we examine 5 m of early Holocene ice from the Dome Fuji (DF) ice core at a high temporal resolution by continuous flow analysis. The ice was continuously analysed for concentrations of dust, sodium, ammonium, liquid conductivity, and water isotopic composition. Furthermore, a dielectric profiling was performed on the solid ice. In most of the analysed ice, the multi-parameter impurity data set appears to resolve the seasonal variability although the identification of annual layers is not always unambiguous. The study thus provides information on the snow accumulation process in central East Antarctica. A layer counting based on the same principles as those previously applied to the NGRIP (North Greenland Ice core Project) and the Antarctic EPICA (European Project for Ice Coring in Antarctica) Dronning Maud Land (EDML) ice cores leads to a mean annual layer thickness for the DF ice of 3.0 ± 0.3 cm that compares well to existing estimates. The measured DF section is linked to the EDML ice core through a characteristic pattern of three significant acidity peaks that are present in both cores. The corresponding section of the EDML ice core has recently been dated by annual layer counting and the number of years identified independently in the two cores agree within error estimates. We therefore conclude that, to first order, the annual signal is preserved in this section of the DF core. This case study demonstrates the feasibility of determining annually deposited strata on the central East Antarctic Plateau. It also opens the possibility of resolving annual layers in the Eemian section of Antarctic ice cores where the accumulation is estimated to have been greater than in the Holocene. © Author(s) 2015.

Description: The widely used detailed SNOWPACK model has undergone constant development over the years. A notable recent extension is the introduction of a Richards equation (RE) solver as an alternative for the bucket-type approach for describing water transport in the snow and soil layers. In addition, continuous updates of snow settling and new snow density parameterizations have changed model behavior. This study presents a detailed evaluation of model performance against a comprehensive multiyear data set from Weissfluhjoch near Davos, Switzerland. The data set is collected by automatic meteorological and snowpack measurements and manual snow profiles. During the main winter season, snow height (RMSE: 〈 4.2 cm), snow water equivalent (SWE, RMSE: 〈 40 mm w.e.), snow temperature distributions (typical deviation with measurements: 〈 1.0 °C) and snow density (typical deviation with observations: 〈 50 kg m−3) as well as their temporal evolution are well simulated in the model and the influence of the two water transport schemes is small. The RE approach reproduces internal differences over capillary barriers but fails to predict enough grain growth since the growth routines have been calibrated using the bucket scheme in the original SNOWPACK model. However, the agreement in both density and grain size is sufficient to parameterize the hydraulic properties successfully. In the melt season, a pronounced underestimation of typically 200 mm w.e. in SWE is found. The discrepancies between the simulations and the field data are generally larger than the differences between the two water transport schemes. Nevertheless, the detailed comparison of the internal snowpack structure shows that the timing of internal temperature and water dynamics is adequately and better represented with the new RE approach when compared to the conventional bucket scheme. On the contrary, the progress of the meltwater front in the snowpack as detected by radar and the temporal evolution of the vertical distribution of melt forms in manually observed snow profiles do not support this conclusion. This discrepancy suggests that the implementation of RE partly mimics preferential flow effects.

Description: A still open question is how equilibrium warming in response to increasing radiative forcing (equilibrium climate sensitivity S) is depending on background climate. We here bring paleo-data based evidence on the state-dependency of S by using CO2 proxy data together with model-based reconstruction of land ice albedo over the last 5 million years. We find that the land-ice albedo forcing depends non-linearly on the background climate, while any non-linearity of CO2 radiative forcings depends on the CO2 data set used. Over the last 2 million years the combined S_[CO2,LI] from CO2 and land-ice albedo forcing is state-dependent and during interglacials at least twice as high as during glacials, thus CO2 doubling leads to an interglacial warming of 5 K. In the Pliocene data uncertainties prevents a well-supported calculation, but our analysis suggests that S_[CO2,LI] during a land-ice free northern hemisphere was smaller than during interglacials of the Pleistocene.

Description: Thermokarst lakes are typical features of the northern permafrost ecosystems, and play an important role in the thermal exchange between atmosphere and subsurface. The objective of this study is to describe the main thermal processes of the lakes and to quantify the heat exchange with the underlying sediments. The thermal regimes of five lakes located within the continuous permafrost zone of northern Siberia (Lena River Delta) were investigated using hourly water temperature and water level records covering a 3-year period (2009–2012), together with bathymetric survey data. The lakes included thermokarst lakes located on Holocene river terraces that may be connected to Lena River water during spring flooding, and a thermokarst lake located on deposits of the Pleistocene Ice Complex. Lakes were covered by ice up to 2 m thick that persisted for more than 7 months of the year, from October until about mid-June. Lake-bottom temperatures increased at the start of the ice-covered period due to upward-directed heat flux from the underlying thawed sediment. Prior to ice break-up, solar radiation effectively warmed the water beneath the ice cover and induced convective mixing. Ice break-up started at the beginning of June and lasted until the middle or end of June. Mixing occurred within the entire water column from the start of ice break-up and continued during the ice-free periods, as confirmed by the Wedderburn numbers, a quantitative measure of the balance between wind mixing and stratification that is important for describing the biogeochemical cycles of lakes. The lake thermal regime was modeled numerically using the FLake model. The model demonstrated good agreement with observations with regard to the mean lake temperature, with a good reproduction of the summer stratification during the ice-free period, but poor agreement during the ice-covered period. Modeled sensitivity to lake depth demonstrated that lakes in this climatic zone with mean depths 〉 5 m develop continuous stratification in summer for at least 1 month. The modeled vertical heat flux across the bottom sediment tends towards an annual mean of zero, with maximum downward fluxes of about 5 W m−2 in summer and with heat released back into the water column at a rate of less than 1 W m−2 during the ice-covered period. The lakes are shown to be efficient heat absorbers and effectively distribute the heat through mixing. Monthly bottom water temperatures during the ice-free period range up to 15 °C and are therefore higher than the associated monthly air or ground temperatures in the surrounding frozen permafrost landscape. The investigated lakes remain unfrozen at depth, with mean annual lake-bottom temperatures of between 2.7 and 4 °C.

Description: Deformation of ice in continental sized ice sheets determines the flow behavior of ice towards the sea. Basal dislocation glide is assumed to be the dominant deformation mechanism in the creep deformation of natural ice, but non-basal glide is active as well. Knowledge of what types of deformation mechanisms are active in polar ice is critical in predicting the response of ice sheets in future warmer climates and its contribution to sea level rise, because the activity of deformation mechanisms depends critically on deformation conditions (such as temperature) as well as on the material properties (such as grain size). One of the methods to study the deformation mechanisms in natural materials is Electron Backscattered Diffraction (EBSD). We obtained ca. 50 EBSD maps of five different depths from a Greenlandic ice core (NEEM). The step size varied between 8 and 25 micron depending on the size of the deformation features. The size of the maps varied from 2000 to 10000 grid point. Indexing rates were up to 95%, partially by saving and reanalyzing the EBSP patterns. With this method we can characterize subgrain boundaries and determine the lattice rotation configurations of each individual subgrain. Combining these observations with arrangement/geometry of subgrain boundaries the dislocation types can be determined, which form these boundaries. Three main types of subgrain boundaries have been recognized in Antarctic (EDML) ice core (Weikusat et al. 2010, 2011). Here, we present the first results obtained from EBSD measurements performed on the NEEM ice core samples from the last glacial period, focusing on the relevance of dislocation activity of the possible slip systems. Preliminary results show that all three subgrain types, recognized in the EDML core, occur in the NEEM samples. In addition to the classical boundaries made up of basal dislocations, subgrain boundaries made of non-basal dislocations are also common.

Description: Ice is a common mineral at the Earth’s surface. How much ice is stored in the Greenland and Antarctic ice sheets depends on its mechanical properties. Therefore properties of ice directly impact on human society through its role in controlling sea level. The bulk behaviour of large ice masses is the result of the behaviour of the ensemble of individual ice grains. This is strongly influenced by the viscoplastic anisotropy of these grains and their lattice orientation. Numerical modelling provides a better insight into the mechanics of ice from the micro to the ice sheet scale. We present numerical simulations that predict the microstructural evolution of an aggregate of pure ice grains at different strain rates. We simulate co-axial deformation and dynamic recrystallization up to large strain using a full-field approach. The crystal plasticity code (Lebensohn et al., 2009) is used to calculate the response of a polycrystalline aggregate that deforms by purely dislocation glide, applying a Fast Fourier Transform (FFT). This code is coupled with the ELLE microstructural modelling platform to include intracrystalline recovery, as well as grain boundary migration driven by the reduction of surface and strain energies. The results show a strong effect of recrystallization on the final microstructure, producing larger and more equiaxed grains, with smooth boundaries. This effect does not significantly modify the single-maximum pattern of c-axes that are distributed at a low angle to the shortening direction. However, in experiments with significant recrystallization the a-axes rotate towards the elongation axis at the same time as the c-axes rotate towards the compression axis. If slip systems on prismatic and/or pyramidal planes are active, it is thought that a-axes gradually concentrate with depth (Miyamoto, 2005). The bulk activity of the slip systems is different depending on the relative activity of deformation versus recrystallization: the non-basal slip systems are more active at high strain in experiments with dynamic recrystallization compared to those experiments with low recrystallization activity.

Description: The ice sheets in Greenland and Antarctica contain a significant amount of air within their upper approximately thousand meters and air hydrates below. While this air is still in exchange with the atmosphere in the permeable firn, the gas is entrapped at the firn-ice transition at 60 – 120 m depth. Understanding the dominant deformation mechanisms is essential to interpret paleo-atmosphere records and to allow a more realistic model of ice sheet dynamics. Recent research shows how the presence of air bubbles can significantly influence microdynamical processes such as grain growth and grain boundary migration (Azuma et al., 2012, Roessiger et al., 2014). Therefore, numerical modelling was performed focussing on the mechanical properties of ice with air inclusions and the implications of the presence of bubbles on recrystallisation. The full-field crystal plasticity code of Lebensohn (2001), using a Fast Fourier Transform (FFT), was coupled to the 2D numerical microstructural modelling platform Elle, following the approach by Griera et al. (2013), and used to simulate dynamic recrystallization of pure ice (Montagnat et al., 2013). FFT calculates the viscoplastic response of polycrystalline and polyphase materials that deform by dislocation glide, takes into account the mechanical anisotropy of ice and calculates dislocation densities using the local gradient of the strain-rate field. Incorporating a code for polyphase grain boundary migration driven by surface and internal strain energy reduction, based on the methodology of Becker et al. (2008) and Roessiger et al. (2014), now also enables us to model deformation of ice with air bubbles. The presence of bubbles leads to an increase in strain localization, which reduces the bulk strength of the bubbly ice. In the absence of dynamic recrystallisation, air bubbles quickly collapse at low strains and spherical to elliptical bubble shapes are only maintained when recrystallisation is activated. Our modelling confirms that strain-induced grain boundary migration already occurs in the uppermost levels of ice sheets (Kipfstuhl et al. 2009, Weikusat et al. 2009).

Description: Ice cores are the only climate archives incorporating paleo-atmosphere as individual gas inclusions, enabling the extraction and analysis of the contained gasses. A firm understanding of the processes involved is mandatory for a reliable interpretation of the gas records. One prominent process is the transition from air bubbles to crystalline air hydrates, which is known to influence, at least temporarily, the gas mixing ratios by diffusion and fractionation. This transition is still not understood completely and the existing theories do not explain the large diversity of observed hydrate morphologies. Raman tomographic measurements using the AWI cryo-Raman system provide 3D reconstructions of air hydrate morphologies. The results show complex growth structures that emphasize the importance of crystallography, microstructure and ice rheology for the hydrate formation process. Accurate hydrate volumes can be calculated from the 3D objects, improving the estimates of total gas contents.

Description: The Asian monsoon system is an important tipping element in Earth's climate with a large impact on human societies in the past and present. In light of the potentially severe impacts of present and future anthropogenic climate change on Asian hydrology, it is vital to understand the forcing mechanisms of past climatic regime shifts in the Asian monsoon domain. Here we use novel recurrence network analysis techniques for detecting episodes with pronounced non-linear changes in Holocene Asian monsoon dynamics recorded in speleothems from caves distributed throughout the major branches of the Asian monsoon system. A newly developed multi-proxy methodology explicitly considers dating uncertainties with the COPRA (COnstructing Proxy Records from Age models) approach and allows for detection of continental-scale regime shifts in the complexity of monsoon dynamics. Several epochs are characterised by non-linear regime shifts in Asian monsoon variability, including the periods around 8.5–7.9, 5.7–5.0, 4.1–3.7, and 3.0–2.4 ka BP. The timing of these regime shifts is consistent with known episodes of Holocene rapid climate change (RCC) and high-latitude Bond events. Additionally, we observe a previously rarely reported non-linear regime shift around 7.3 ka BP, a timing that matches the typical 1.0–1.5 ky return intervals of Bond events. A detailed review of previously suggested links between Holocene climatic changes in the Asian monsoon domain and the archaeological record indicates that, in addition to previously considered longer-term changes in mean monsoon intensity and other climatic parameters, regime shifts in monsoon complexity might have played an important role as drivers of migration, pronounced cultural changes, and the collapse of ancient human societies.

Description: Thermokarst lakes are important emitters of methane, a potent greenhouse gas. However, accurate estimation of methane flux from thermokarst lakes is difficult due to their remoteness and observational challenges associated with the heterogeneous nature of ebullition (bubbling). We used multi-temporal high-resolution (9–11 cm) aerial images of an interior Alaskan thermokarst lake, Goldstream Lake, acquired 2 and 4 days following freeze-up in 2011 and 2012, respectively, to characterize methane ebullition seeps and to estimate whole-lake ebullition. Bubbles impeded by the lake ice sheet form distinct white patches as a function of bubbling rate vs. time as ice thickens. Our aerial imagery thus captured in a single snapshot the ebullition events that occurred before the image acquisition. Image analysis showed that low-flux A- and B-type seeps are associated with low brightness patches and are statistically distinct from high-flux C-type and Hotspot seeps associated with high brightness patches. Mean whole-lake ebullition based on optical image analysis in combination with bubble-trap flux measurements was estimated to be 174 ± 28 and 216 ± 33 mL gas m−2 d−1 for the years 2011 and 2012, respectively. A large number of seeps demonstrated spatio-temporal stability over our two-year study period. A strong inverse exponential relationship (R2 ≥ 0.79) was found between percent surface area of lake ice covered with bubble patches and distance from the active thermokarst lake margin. Our study shows that optical remote sensing is a powerful tool to map ebullition seeps on lake ice, to identify their relative strength of ebullition and to assess their spatio-temporal variability.

Description: Source attribution and process analysis for atmospheric mercury in eastern China simulated by CMAQ-Hg Atmospheric Chemistry and Physics, 15, 8767-8779, 2015 Author(s): J. Zhu, T. Wang, J. Bieser, and V. Matthias The contribution from different emission sources and atmospheric processes to gaseous elemental mercury (GEM), gaseous oxidized mercury (GOM), particulate bound mercury (PBM) and mercury deposition in eastern China were quantified using the Community Multi-scale Air Quality (CMAQ-Hg) modeling system run with a nested domain. Natural sources (NAT) and six categories of anthropogenic mercury sources (ANTH) including cement production (CEM), domestic life (DOM), industrial boilers (IND), metal production (MET), coal-fired power plants (PP) and traffic (TRA) were considered for source apportionment. NAT were responsible for 36.6 % of annual averaged GEM concentration, which was regarded as the most important source for GEM in spite of obvious seasonal variation. Among ANTH, the influence of MET and PP on GEM were most evident especially in winter. ANTH dominated the variations of GOM and PBM concentrations with contributions of 86.7 and 79.1 %, respectively. Among ANTH, IND were the largest contributor for GOM (57.5 %) and PBM (34.4 %) so that most mercury deposition came from IND. The effect of mercury emitted from out of China was indicated by a 〉 30 % contribution to GEM concentration and wet deposition. The contributions from nine processes – consisting of emissions (EMIS), gas-phase chemical production/loss (CHEM), horizontal advection (HADV), vertical advection (ZADV), horizontal advection (HDIF), vertical diffusion (VDIF), dry deposition (DDEP), cloud processes (CLDS) and aerosol processes (AERO) – were calculated for process analysis with their comparison in urban and non-urban regions of the Yangtze River delta (YRD). EMIS and VDIF affected surface GEM and PBM concentrations most and tended to compensate each other all the time in both urban and non-urban areas. However, DDEP was the most important removal process for GOM with 7.3 and 2.9 ng m −3 reduced in the surface of urban and non-urban areas, respectively, in 1 day. The diurnal profile variation of processes revealed the transportation of GOM from urban area to non-urban areas and the importance of CHEM/AERO in higher altitudes which partly caused diffusion of GOM downwards to non-urban areas. Most of the anthropogenic mercury was transported and diffused away from urban areas by HADV and VDIF and increased mercury concentrations in non-urban areas by HADV. Natural emissions only influenced CHEM and AERO more significantly than anthropogenic. Local emissions in the YRD contributed 8.5 % more to GEM and ~ 30 % more to GOM and PBM in urban areas compared to non-urban areas.

Description: Seasonal variation of secondary organic aerosol tracers in Central Tibetan Plateau Atmospheric Chemistry and Physics, 15, 8781-8793, 2015 Author(s): R.-Q. Shen, X. Ding, Q.-F. He, Z.-Y. Cong, Q.-Q. Yu, and X.-M. Wang Secondary organic aerosol (SOA) affects the earth's radiation balance and global climate. High-elevation areas are sensitive to global climate change. However, at present, SOA origins and seasonal variations are understudied in remote high-elevation areas. In this study, particulate samples were collected from July 2012 to July 2013 at the remote Nam Co (NC) site, Central Tibetan Plateau and analyzed for SOA tracers from biogenic (isoprene, monoterpenes and β-caryophyllene) and anthropogenic (aromatics) precursors. Among these compounds, isoprene SOA (SOA I ) tracers represented the majority (26.6 ± 44.2 ng m −3 ), followed by monoterpene SOA (SOA M ) tracers (0.97 ± 0.57 ng m −3 ), aromatic SOA (SOA A ) tracer (2,3-dihydroxy-4-oxopentanoic acid, DHOPA, 0.25 ± 0.18 ng m −3 ) and β-caryophyllene SOA tracer (β-caryophyllenic acid, 0.09 ± 0.10 ng m −3 ). SOA I tracers exhibited high concentrations in the summer and low levels in the winter. The similar temperature dependence of SOA I tracers and isoprene emission suggested that the seasonal variation of SOA I tracers at the NC site was mainly influenced by the isoprene emission. The ratio of high-NO x to low-NO x products of SOA I (2-methylglyceric acid to 2-methyltetrols) was highest in the winter and lowest in the summer, due to the influence of temperature and relative humidity. The seasonal variation of SOA M tracers was impacted by monoterpenes emission and gas-particle partitioning. During the summer to the fall, temperature effect on partitioning was the dominant process influencing SOA M tracers' variation; while the temperature effect on emission was the dominant process influencing SOA M tracers' variation during the winter to the spring. SOA M tracer levels did not elevate with increased temperature in the summer, probably resulting from the counteraction of temperature effects on emission and partitioning. The concentrations of DHOPA were 1–2 orders of magnitude lower than those reported in the urban regions of the world. Due to the transport of air pollutants from the adjacent Bangladesh and northeastern India, DHOPA presented relatively higher levels in the summer. In the winter when air masses mainly came from northwestern India, mass fractions of DHOPA in total tracers increased, although its concentrations declined. The SOA-tracer method was applied to estimate secondary organic carbon (SOC) from these four precursors. The annual average of SOC was 0.22 ± 0.29 μgC m −3 , with the biogenic SOC (sum of isoprene, monoterpenes and β-caryophyllene) accounting for 75 %. In the summer, isoprene was the major precursor with its SOC contributions of 81 %. In the winter when the emission of biogenic precursors largely dropped, the contributions of aromatic SOC increased. Our study implies that anthropogenic pollutants emitted in the Indian subcontinent could be transported to the TP and have an impact on SOC over the remote NC.

Description: Influence of biomass burning on CCN number and hygroscopicity during summertime in the eastern Mediterranean Atmospheric Chemistry and Physics Discussions, 15, 21539-21582, 2015 Author(s): A. Bougiatioti, S. Bezantakos, I. Stavroulas, N. Kalivitis, P. Kokkalis, G. Biskos, N. Mihalopoulos, A. Papayannis, and A. Nenes This study investigates the CCN activity and hygroscopic properties of particles influenced by biomass burning in the eastern Mediterranean. Air masses sampled were subject to a range of atmospheric processing (several hours up to 3 days). Values of the hygroscopicity parameter, κ, were derived from cloud condensation nuclei (CCN) measurements and a Hygroscopic Tandem Differential Mobility Analyzer (HTDMA). An Aerosol Chemical Speciation Monitor (ACSM) was also used to determine the chemical composition and mass concentration of non-refractory components of the submicron aerosol fraction. During fire events, the increased organic content (and lower inorganic fraction) of the aerosol decreases the hygroscopicity parameter, κ, for all particle sizes. The reason, however, for this decrease was not the same for all size modes; smaller particle sizes appeared to be richer in less hygroscopic, less CCN-active components due to coagulation processes while larger particles become less hygroscopic during the biomass burning events due to condensation of less hygroscopic gaseous compounds. In addition, smaller particles exhibited considerable chemical dispersion (where hygroscopicity varied up to 100 % for particles of same size); larger particles, however, exhibited considerably less dispersion owing to the effects of aging and retained high levels of CCN activity. These conclusions are further supported by the observed mixing state determined by the HTDMA measurements. ACSM measurements indicate that the bulk composition reflects the hygroscopicity and chemical nature of the largest particles and a large fraction of the CCN concentrations sampled. Based on Positive Matrix Factorization (PMF) analysis of the organic ACSM spectra, CCN concentrations follow a similar trend with the BBOA component, with enhancements of CCN in biomass burning plumes ranging between 65 and 150 %, for supersaturations ranging between 0.2 and 0.7 %. Using multilinear regression, we determine the hygroscopicity of the prime organic aerosol components (BBOA, OOA-BB and OOA); it is found that the total organic hygroscopicity is very close to the inferred hygroscopicity of the oxygenated organic aerosol components. Finally, the transformation of freshly-emitted biomass burning (BBOA) to more oxidized organic aerosol (OOA-BB) can result in a two-fold increase of the inferred organic hygroscopicity. Almost 10 % of the total aerosol hygroscopicity is related to the two biomass burning components (BBOA and OOA-BB), which in turn contribute almost 35 % to the fine-particle organic water of the aerosol. This is important as organic water can contribute to the atmospheric chemistry and the direct radiative forcing.

Description: Measurement-based direct radiative effect by brown carbon over Indo-Gangetic Plain Atmospheric Chemistry and Physics Discussions, 15, 21583-21606, 2015 Author(s): A. Arola, G. L. Schuster, M. R. A. Pitkänen, O. Dubovik, H. Kokkola, A. V. Lindfors, T. Mielonen, T. Raatikainen, S. Romakkaniemi, S. N. Tripathi, and H. Lihavainen The importance of light absorbing organic aerosols, often called brown carbon (BrC), has become evident in recent years. However, there are relatively few measurement-based estimates for the direct radiative effect of BrC so far. In those earlier studies, the AErosol RObotic NETwork (AERONET) measured Aerosol Absorption Optical Depth (AAOD) and Absorption Angstrom Exponent (AAE) have been exploited. However, these two pieces of information are clearly not sufficient to separate properly carbonaceous aerosols from dust, while imaginary indices of refraction would contain more and better justified information for this purpose. This is first time that the direct radiative effect (DRE) of BrC is estimated by exploiting the AERONET-retrieved imaginary indices. We estimated it for four sites in Indo-Gangetic Plain (IGP), Karachi, Lahore, Kanpur and Gandhi College. We found a distinct seasonality, which was generally similar among all the sites, but with slightly different strengths. The monthly warming effect up to 0.5 W m -2 takes place during spring season. On the other hand, BrC results in overall cooling effect in the winter season, which can reach levels close to −1W m -2 . We then estimated similarly also DRE of black carbon and total aerosol, in order to assess the relative significance of BrC radiative effect in the radiative effects of other components. Even though BrC impact seems minor in this context, we demonstrated that it is not insignificant and moreover that it is crucial to perform spectrally resolved radiative transfer calculations to obtain good estimates for DRE of BrC.

Description: Contrail life cycle and properties from 1 year of MSG/SEVIRI rapid-scan images Atmospheric Chemistry and Physics, 15, 8739-8749, 2015 Author(s): M. Vázquez-Navarro, H. Mannstein, and S. Kox The automatic contrail tracking algorithm (ACTA) – developed to automatically follow contrails as they age, drift and spread – enables the study of a large number of contrails and the evolution of contrail properties with time. In this paper we present a year's worth of tracked contrails, from August 2008 to July 2009 in order to derive statistically significant mean values. The tracking is performed using the 5 min rapid-scan mode of the Spinning Enhanced Visible and Infrared Imager (SEVIRI) on board the Meteosat Second Generation (MSG) satellites. The detection is based on the high spatial resolution of the images provided by the Moderate Resolution Imaging Spectroradiometer on board the Terra satellite (Terra/MODIS), where a contrail detection algorithm (CDA) is applied. The results show the satellite-derived average lifetimes of contrails and contrail-cirrus along with the probability density function (PDF) of other geometric characteristics such as mean coverage, distribution and width. In combination with specifically developed algorithms (RRUMS; Rapid Retrieval of Upwelling irradiance from MSG/SEVIRI and COCS (Cirrus Optical properties derived from CALIOP and SEVIRI), explained below) it is possible to derive the radiative forcing (RF), energy forcing (EF), optical thickness (τ) and altitude of the tracked contrails. Mean values here retrieved are duration, 1 h; length, 130 km; width, 8 km; altitude, 11.7 km; optical thickness, 0.34. Radiative forcing and energy forcing are shown for land/water backgrounds in day/night situations.

Description: Accounting for the effects of sastrugi in the CERES clear-sky Antarctic shortwave angular distribution models Atmospheric Measurement Techniques, 8, 3163-3175, 2015 Author(s): J. Corbett and W. Su The Cloud and the Earth's Radiant Energy System (CERES) instruments on NASA's Terra, Aqua and Soumi NPP satellites are used to provide a long-term measurement of Earth's energy budget. To accomplish this, the radiances measured by the instruments must be inverted to fluxes by the use of a scene-type-dependent angular distribution model (ADM). For permanent snow scenes over Antarctica, shortwave (SW) ADMs are created by compositing radiance measurements over the full viewing zenith and azimuth range. However, the presence of small-scale wind blown roughness features called sastrugi cause the BRDF (bidirectional reflectance distribution function) of the snow to vary significantly based upon the solar azimuth angle and location. This can result in monthly regional biases between −12 and 7.5 Wm −2 in the inverted TOA (top-of-atmosphere) SW flux. The bias is assessed by comparing the CERES shortwave fluxes derived from nadir observations with those from all viewing zenith angles, as the sastrugi affect fluxes inverted from the oblique viewing angles more than for the nadir viewing angles. In this paper we further describe the clear-sky Antarctic ADMs from Su et al. (2015). These ADMs account for the sastrugi effect by using measurements from the Multi-Angle Imaging Spectro-Radiometer (MISR) instrument to derive statistical relationships between radiance from different viewing angles. We show here that these ADMs reduce the bias and artifacts in the CERES SW flux caused by sastrugi, both locally and Antarctic-wide. The regional monthly biases from sastrugi are reduced to between −5 and 7 Wm −2 , and the monthly-mean biases over Antarctica are reduced by up to 0.64 Wm −2 , a decrease of 74 %. These improved ADMs are used as part of the Edition 4 CERES SSF (Single Scanner Footprint) data.

Description: Space-borne observation of methane from atmospheric infrared sounder version 6: validation and implications for data analysis Atmospheric Measurement Techniques Discussions, 8, 8563-8597, 2015 Author(s): X. Xiong, F. Weng, Q. Liu, and E. Olsen Atmospheric Methane (CH 4 ) is generated as a standard product in recent version of the hyperspectral Atmospheric Infrared Sounder (AIRS-V6) aboard NASA's Aqua satellite at the NASA Goddard Earth Sciences Data and Information Services Center (NASA/GES/DISC). Significant improvements in AIRS-V6 was expected but without a thorough validation. This paper first introduced the improvements of CH 4 retrieval in AIRS-V6 and some characterizations, then presented the results of validation using ~ 1000 aircraft profiles from several campaigns spread over a couple of years and in different regions. It was found the mean biases of AIRS CH 4 at layers 343–441 and 441–575 hPa are −0.76 and −0.05 % and the RMS errors are 1.56 and 1.16 %, respectively. Further analysis demonstrates that the errors in the spring and in the high northern latitudes are larger than in other seasons or regions. The error is correlated with Degree of Freedoms (DOFs), particularly in the tropics or in the summer, and cloud amount, suggesting that the "observed" spatiotemporal variation of CH 4 by AIRS is imbedded with some artificial impact from the retrieval sensitivity in addition to its variation in reality, so the variation of information content in the retrievals needs to be taken into account in data analysis of the retrieval products. Some additional filtering (i.e. rejection of profiles with obvious oscillation as well as those deviating greatly from the norm) for quality control is recommended for the users to better utilize AIRS-V6 CH 4 , and their implementation in the future versions of the AIRS retrieval algorithm is under consideration.

Description: Impact of aerosols on the OMI tropospheric NO 2 retrievals over industrialized regions: how accurate is the aerosol correction of cloud-free scenes via a simple cloud model? Atmospheric Measurement Techniques Discussions, 8, 8385-8437, 2015 Author(s): J. Chimot, T. Vlemmix, J. P. Veefkind, J. F. de Haan, and P. F. Levelt The Ozone Monitoring Instrument (OMI) instrument has provided daily global measurements of tropospheric NO 2 for more than a decade. Numerous studies have drawn attention to the complexities related to measurements of tropospheric NO 2 in the presence of aerosols. Fine particles affect the OMI spectral measurements and the length of the average light path followed by the photons. However, they are not explicitly taken into account in the current OMI tropospheric NO 2 retrieval chain. Instead, the operational OMI O 2 -O 2 cloud retrieval algorithm is applied both to cloudy scenes and to cloud free scenes with aerosols present. This paper describes in detail the complex interplay between the spectral effects of aerosols, the OMI O 2 -O 2 cloud retrieval algorithm and the impact on the accuracy of the tropospheric NO 2 retrievals through the computed Air Mass Factor (AMF) over cloud-free scenes. Collocated OMI NO 2 and MODIS Aqua aerosol products are analysed over East China, in industrialized area. In addition, aerosol effects on the tropospheric NO 2 AMF and the retrieval of OMI cloud parameters are simulated. Both the observation-based and the simulation-based approach demonstrate that the retrieved cloud fraction linearly increases with increasing Aerosol Optical Thickness (AOT), but the magnitude of this increase depends on the aerosol properties and surface albedo. This increase is induced by the additional scattering effects of aerosols which enhance the scene brightness. The decreasing effective cloud pressure with increasing AOT represents primarily the absorbing effects of aerosols. The study cases show that the actual aerosol correction based on the implemented OMI cloud model results in biases between −20 and −40 % for the DOMINO tropospheric NO 2 product in cases of high aerosol pollution (AOT ≥ 0.6) and elevated particles. On the contrary, when aerosols are relatively close to the surface or mixed with NO 2 , aerosol correction based on the cloud model results in overestimation of the DOMINO tropospheric NO 2 product, between 10 and 20 %. These numbers are in line with comparison studies between ground-based and OMI tropospheric NO 2 measurements under conditions with high aerosol pollution and elevated particles. This highlights the need to implement an improved aerosol correction in the computation of tropospheric NO 2 AMFs.

Description: Sensitivity of thermal infrared sounders to the chemical and micro-physical properties of UTLS secondary sulphate aerosols Atmospheric Measurement Techniques Discussions, 8, 8439-8481, 2015 Author(s): P. Sellitto and B. Legras Monitoring upper tropospheric-lower stratospheric (UTLS) secondary sulphate aerosols and their chemical and micro-physical properties from satellite nadir observations is crucial to better understand their formation and evolution processes and then to estimate their impact to the UTLS chemistry, and on regional and global radiative balance. Here we present a study aimed at the evaluation of the sensitivity of thermal infrared (TIR) satellite nadir observations to the chemical composition and the size distribution of idealized UTLS sulphate aerosol layers. The extinction properties of sulphuric acid/water droplets, for different sulphuric acid mixing ratios and temperatures, are systematically analysed. The extinction coefficients are derived by means of a Mie code, using refractive indexes taken from the GEISA (Gestion et Étude des Informations Spectroscopiques Atmosphériques: Management and Study of Spectroscopic Information) spectroscopic database and log-normal size distributions with different effective radii and number concentrations. IASI (Infrared Atmospheric Sounding Interferometer) pseudo-observations are generated using forward radiative transfer calculations performed with the 4A (Automatized Atmospheric Absorption Atlas) radiative transfer model, to estimate the impact of the extinction of idealized aerosol layers, at typical UTLS conditions, on the brightness temperature spectra observed by this satellite instrument. We found a marked and typical spectral signature of these aerosol layers between 700 and 1200 cm −1 , due to the absorption bands of the sulphate and bi-sulphate ions and the undissociated sulphuric acid, with the main absorption peaks at 1170 and 905 cm −1 . The dependence of the aerosol spectral signature to the sulphuric acid mixing ratio, and effective number concentration and radius, as well as the role of interferring parameters like the ozone, sulphur dioxide, carbon dioxide and ash absorption, and temperature and water vapour profile uncertainties, are analyzed and critically discussed.

Description: Global cloud top height retrieval using SCIAMACHY limb spectra: model studies and first results Atmospheric Measurement Techniques Discussions, 8, 8295-8352, 2015 Author(s): K.-U. Eichmann, L. Lelli, C. von Savigny, H. Sembhi, and J. P. Burrows Cloud top heights (CTH) were retrieved for the period 1 January 2003 to 7 April 2012 using height-resolved limb spectra measured with the Scanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIAMACHY) on board ENVISAT (ENVIronmental SATellite). In this study, we tested the sensitivity of the colour index method used in the retrieval code SCODA (SCIAMACHY Cloud Detection Algorithm) and the accuracy of the retrieved CTHs in comparison to other methods. Sensitivity studies using the radiative transfer model SCIATRAN showed that the method is capable of generally detecting cloud tops down to about 5 km and very thin cirrus clouds even up to the tropopause. Volcanic particles can also be detected that occasionally reach the lower stratosphere. Low clouds at 2–3 km can only be retrieved under very clean atmospheric conditions, as light scattering of aerosols interferes with the cloud retrieval. Upper tropospheric ice clouds are detectable for cloud optical depths down to about τ N = 0.005, which is in the subvisual range. The detection sensitivity decreases towards the surface. An optical thickness of roughly 0.1 was the lower detection limit for water cloud top heights at 5 km. This value is much lower than thresholds reported for the passive cloud detection in nadir viewing direction. Comparisons with SCIAMACHY nadir cloud top heights, calculated with the Semi-Analytical CloUd Retrieval Algorithm (SACURA), showed a good agreement in the global cloud field distribution. But only opaque clouds (τ N 〉 5) are detectable with the nadir passive retrieval technique in the UV-visible and infrared wavelength range. So due to the frequent occurrence of thin and sub-visual cirrus clouds in the tropics, large cloud top height deviations were detected between both viewing geometries. Also the land/sea contrast seen in nadir retrievals was not detected in limb mode. Co-located cloud top height measurements of the limb viewing Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) on ENVISAT for the period from January 2008 to March 2012 were compared, showing good agreement to within 1 km, which is smaller than the vertical field of view of both instruments. Lower stratospheric aerosols from volcanic eruptions occasionally interfered with the cloud retrieval and inhibited detection of tropospheric clouds. Examples of the impact of these events are shown for the volcanoes Kasatochi in August 2008, Sarychev Peak in June 2009, and Nabro in June 2010. Long-lasting aerosol layers were detected after these events in the Northern Hemisphere down to the tropics. Particle top heights up to about 22 km were retrieved in 2009, when the enhanced lower stratospheric aerosol layer persisted for about 7 months. Up to about 82 % of the Northern hemispheric lower stratosphere between 30° and 70° was covered by scattering particles in August 2009 and nearly half in October 2008.

Description: Operation of the Airmodus A11 nano Condensation Nucleus Counter at various inlet pressures, various operation temperatures and design of a new inlet system Atmospheric Measurement Techniques Discussions, 8, 8483-8508, 2015 Author(s): J. Kangasluoma, A. Franchin, J. Duplissy, L. Ahonen, F. Korhonen, M. Attoui, J. Mikkilä, K. Lehtipalo, J. Vanhanen, M. Kulmala, and T. Petäjä Measuring sub-3 nm particles outside of controlled laboratory conditions is a challenging task, as many of the instruments are operated at their limits and are subjected to changing ambient conditions. In this study, we advance the current understanding on the operation of Airmodus A11 nano Condensation Nucleus Counter (nCNC), which consists of a A10 Particle Size Magnifier (PSM) and A20 condensation particle counter (CPC). We explore the effect of the inlet line pressure on the measured particle concentration. We identify two different regions inside the instrument where supersaturation of working fluid can take place. We show the possibility of varying the cut-off of the instrument from 1 to 6 nm, a wider size range than the one usually covered by the PSM. We also present a new inlet system, which allows automated measurements of the background, minimizes the diffusion losses in the sampling line and is equipped with an electrostatic filter to remove ions. Finally, our view of the guidelines for optimal use of the Airmodus nCNC are provided.

Description: Looking beyond stratification: a model-based analysis of the biological drivers of oxygen depletion in the North Sea Biogeosciences Discussions, 12, 12543-12610, 2015 Author(s): F. Große, N. Greenwood, M. Kreus, H. J. Lenhart, D. Machoczek, J. Pätsch, L. A. Salt, and H. Thomas The problem of low oxygen conditions, often referred to as hypoxia, occurs regularly in the North Sea, a temperate European shelf sea. Stratification represents a major process regulating the seasonal dynamics of bottom oxygen. However, lowest oxygen conditions in the North Sea do not occur in the regions of strongest stratification. This suggests that stratification is an important prerequisite for hypoxia, but that the complex interaction between hydrodynamics and the biological processes drives its development. In this study we use the ecosystem model HAMSOM-ECOHAM5 to provide a general characteristic of the different North Sea oxygen regimes, and to quantify the impact of the different physical and biological factors driving the oxygen dynamics below the thermocline and in the bottom layer. We show that the North Sea can be subdivided into three different regimes in terms of oxygen dynamics: (1) a highly productive, non-stratified coastal regime, (2) a productive, seasonally stratified regime with a small sub-thermocline volume, and (3) a productive, seasonally stratified regime with a large sub-thermocline volume, with regime 2 being highly susceptible to hypoxic conditions. Our analysis of the different processes driving the oxygen development reveals that inter-annual variations in the oxygen conditions are caused by variations in primary production, while spatial differences can be attributed to differences in stratification and water depth. In addition, we show that benthic bacteria represent the main oxygen consumers in the bottom layer, consistently accounting for more than 50 % of the overall consumption. By providing these valuable insights, we show that ecosystem models can be a useful tool for the interpretation of observations and the estimation of the impact of anthropogenic drivers on the North Sea oxygen conditions.

Description: Impact of the oceanic geothermal heat flux on a glacial ocean state Climate of the Past Discussions, 11, 3597-3624, 2015 Author(s): M. Ballarotta, F. Roquet, S. Falahat, Q. Zhang, and G. Madec The oceanic geothermal heating (OGH) has a significant impact on the present-day ocean state, but its role during glacial periods, when the ocean circulation and stratification were different from those of today, remains poorly known. In the present study, we analyzed the response of the glacial ocean to OGH, by comparing ocean simulations of the Last Glacial Maximum (LGM, ∼ 21 ka ago) including or not geothermal heating. We found that applying the OGH warmed the Antarctic Bottom Waters (AABW) by ∼ 0.4 °C and increased the abyssal circulation by 15 to 30 % north of 30° S in the deep Pacific and Atlantic basins. The geothermally heated deep waters were then advected toward the Southern Ocean where they upwelled to the surface due to the Ekman transport. The extra heat transport towards Antarctica acted to reduce the amount of sea ice contributing to the freshening of the whole AABW overturning cell. The global amount of salt being conserved, this bottom freshening induced a salinification of the North Atlantic and North Pacific surface and intermediate waters, contributing to the deepening of the North Atlantic Deep Water. This indirect mechanism is responsible for the largest observed warming, found in the North Atlantic deep western boundary current between 2000 and 3000 m (up to 2 °C). The characteristic time scale of the ocean response to the OGH corresponds to an advective time scale (associated with the overturning of the AABW cell) rather than a diffusive time scale. The OGH might facilitate the transition from a glacial to an inter-glacial state but its effect on the deep stratification seems insufficient to drive alone an abrupt climate change.

Description: Woodland Survey of Great Britain 1971–2001 Earth System Science Data, 7, 203-214, 2015 Author(s): C. M. Wood, S. M. Smart, and R. G. H. Bunce The Woodland Survey of Great Britain is a unique data set, consisting of a detailed range of ecological measurements at a national scale, covering a time span of 30 years. A set of 103 woods spread across Britain were first surveyed in 1971, which were again surveyed in 2000–2003 (for convenience referred to subsequently as the "2001 survey"). Standardised methods of describing the trees, shrubs, ground flora, soils and general habitats present were used for both sets of surveys. The sample of 1648 plots spread through 103 woodland sites located across Britain makes it probably the most extensive quantitative ecological woodland survey undertaken in Britain; it is also notable for the range of sites that have been revisited after such a long interval. The data set provides a unique opportunity to explore the effects of a range of potential drivers of woodland change that operated between 1971 and 2001. The data set is available in four discrete parts, which have been assigned the following DOIs: doi:10.5285/4d93f9ac-68e3-49cf-8a41-4d02a7ead81a (Kirby et al., 2013b), doi:10.5285/d6409d40-58fe-4fa7-b7c8-71a105b965b4 (Kirby et al., 2013d), doi:10.5285/fb1e474d-456b-42a9-9a10-a02c35af10d2 (Kirby et al., 2013c), doi:10.5285/2d023ce9-6dbe-4b4f-a0cd-34768e1455ae (Kirby et al., 2013a).

Description: Estimating greenhouse gas emissions from travel – a GIS-based study Geographica Helvetica, 70, 185-192, 2015 Author(s): S. Kuonen Conferences, meetings and congresses are an important part of today's economic and scientific world. But the environmental impact, especially from greenhouse gas emissions associated with travel, can be extensive. Anthropogenic greenhouse gas (GHG) emissions account for the warming of the atmosphere and oceans. This study draws on the need to quantify and reduce greenhouse gas emissions associated with travel activities and aims to give suggestions for organizers and participants on possible ways to reduce greenhouse gas emissions, demonstrated on the example of the European Geography Association (EGEA) Annual Congress 2013 in Wasilkow, Poland. The lack of a comprehensive methodology for the estimation of greenhouse gas emissions from travel led to an outline of a methodology that uses geographic information systems (GIS) to calculate travel distances. The calculation of travel distances in GIS is adapted from actual transportation infrastructure, derived from the open-source platform OpenStreetMap. The methodology also aims to assess the possibilities to reduce GHG emissions by choosing different means of transportation and a more central conference location. The results of the participants of the EGEA congress, who shared their travel data for this study, show that the total travel distance adds up to 238 000 km, with average travel distance of 2429 km per participant. The travel activities of the participants in the study result in total GHG emissions of 39 300 kg CO 2 -eq including both outward and return trip. On average a participant caused GHG emissions of 401 kg CO 2 -eq. In addition, the analysis of the travel data showed differences in travel behaviour depending on the distance between conference site and point of origin. The findings on travel behaviour have then been used to give an estimation of total greenhouse gas emissions from travel for all participants of the conference, which result in a total amount of 79 711 kg CO 2 -eq. The potential for reducing greenhouse gas emissions by substituting short flights with train rides and car rides with bus and train rides is limited. Only 6 % of greenhouse gas emissions could be saved by applying these measures. Further considerable savings could only be made by substituting longer flights (32.6 %) or choosing a more central conference location (26.3 %).

Description: Singular vector based targeted observations of chemical constituents: description and first application of the EURAD-IM-SVA Geoscientific Model Development Discussions, 8, 6267-6307, 2015 Author(s): N. Goris and H. Elbern Measurements of the large dimensional chemical state of the atmosphere provide only sparse snapshots of the state of the system due to their typically insufficient temporal and spatial density. In order to optimize the measurement configurations despite those limitations, the present work describes the identification of sensitive states of the chemical system as optimal target areas for adaptive observations. For this purpose, the technique of singular vector analysis (SVA), which has been proved effective for targeted observations in numerical weather predication, is implemented into the chemical transport model EURAD-IM (EURopean Air pollution and Dispersion – Inverse Model) yielding the EURAD-IM-SVA. Besides initial values, emissions are investigated as critical simulation controlling targeting variables. For both variants, singular vectors are applied to determine the optimal placement for observations and moreover to quantify which chemical compounds have to be observed with preference. Based on measurements of the airship based ZEPTER-2 campaign, the EURAD-IM-SVA has been evaluated by conducting a comprehensive set of model runs involving different initial states and simulation lengths. Since the considered cases are restricted in terms of considered chemical compounds and selected areas, they allow for a retracing of the results and a confirmation of their correctness. Our analysis shows that the optimal placement for observations of chemical species is not entirely determined by mere transport and mixing processes. Rather, a combination of initial chemical concentrations, chemical conversions, and meteorological processes determine the influence of chemical compounds and regions. We furthermore demonstrate that the optimal placement of observations of emission strengths is highly dependent on the location of emission sources and that the benefit of including emissions as target variables outperforms the value of initial value optimisation with growing simulation length. The obtained results confirm the benefit of considering both initial values and emission strengths as target variables and of applying the EURAD-IM-SVA for measurement decision guidance with respect to chemical compounds.

Description: A global scale mechanistic model of the photosynthetic capacity Geoscientific Model Development Discussions, 8, 6217-6266, 2015 Author(s): A. A. Ali, C. Xu, A. Rogers, R. A. Fisher, S. D. Wullschleger, N. G. McDowell, E. C. Massoud, J. A. Vrugt, J. D. Muss, J. B. Fisher, P. B. Reich, and C. J. Wilson Although plant photosynthetic capacity as determined by the maximum carboxylation rate (i.e., V c, max25 ) and the maximum electron transport rate (i.e., J max25 ) at a reference temperature (generally 25 °C) is known to vary substantially in space and time in response to environmental conditions, it is typically parameterized in Earth system models (ESMs) with tabulated values associated to plant functional types. In this study, we developed a mechanistic model of leaf utilization of nitrogen for assimilation (LUNA V1.0) to predict the photosynthetic capacity at the global scale under different environmental conditions, based on the optimization of nitrogen allocated among light capture, electron transport, carboxylation, and respiration. The LUNA model was able to reasonably well capture the observed patterns of photosynthetic capacity in view that it explained approximately 55 % of the variation in observed V c, max25 and 65 % of the variation in observed J max25 across the globe. Our model simulations under current and future climate conditions indicated that V c, max25 could be most affected in high-latitude regions under a warming climate and that ESMs using a fixed V c, max25 or J max25 by plant functional types were likely to substantially overestimate future global photosynthesis.

Description: Uncertainties in calculating precipitation climatology in East Asia Hydrology and Earth System Sciences Discussions, 12, 7765-7783, 2015 Author(s): J. Kim and S. K. Park This study examines the uncertainty in calculating the fundamental climatological characteristics of precipitation in the East Asia region from multiple fine-resolution gridded analysis datasets based on in-situ rain gauge observations. Five observation-based gridded precipitation datasets are used to derive the long-term means, standard deviations in lieu of interannual variability and linear trends over the 28-year period from 1980 to 2007. Both the annual and summer (June–July–August) mean precipitation is examined. The agreement amongst these precipitation datasets are examined using multiple metrics including the signal-to-noise ratio (SNR) defined as the ratio between long-term means and the corresponding standard deviations, and Taylor diagrams which allows examinations of the pattern correlation, the standard deviation, and the centered root mean square error. It is found that the five gauge-based precipitation analysis datasets agree well in the long-term mean and interannual variability in most of the East Asia region including eastern China, Manchuria, South Korea, and Japan, which are densely populated and have fairly high density observation networks. The regions of large inter-dataset variations include Tibetan Plateau, Mongolia, northern Indo-China, and North Korea. The regions of large uncertainties are typically lightly populated and are characterized by severe terrain and/or extreme high elevations. Unlike the long-term mean and interannual variability, agreements between datasets in the linear trend is weak, both for the annual and summer mean values. In most of the East Asia region, the SNR for the linear trend is below 0.5, i.e., the inter-dataset variability exceeds the multi-data ensemble mean. The uncertainty in the spatial distribution of long-term means among these datasets occurs both in the spatial pattern and variability, but the uncertainty for the interannual variability and time trend is much larger in the variability than in the pattern correlation. Thus, care must be taken in using long-term trends calculated from gridded precipitation analysis data for climate studies over the East Asia region.

Description: Qualitative soil moisture assessment in semi-arid Africa – the role of experience and training on inter-rater reliability Hydrology and Earth System Sciences, 19, 3505-3516, 2015 Author(s): M. Rinderer, H. C. Komakech, D. Müller, G. L. B. Wiesenberg, and J. Seibert Soil and water management is particularly relevant in semi-arid regions to enhance agricultural productivity. During periods of water scarcity, soil moisture differences are important indicators of the soil water deficit and are traditionally used for allocating water resources among farmers of a village community. Here we present a simple, inexpensive soil wetness classification scheme based on qualitative indicators which one can see or touch on the soil surface. It incorporates the local farmers' knowledge on the best soil moisture conditions for seeding and brick making in the semi-arid environment of the study site near Arusha, Tanzania. The scheme was tested twice in 2014 with farmers, students and experts (April: 40 persons, June: 25 persons) for inter-rater reliability, bias of individuals and functional relation between qualitative and quantitative soil moisture values. During the test in April farmers assigned the same wetness class in 46 % of all cases, while students and experts agreed on about 60 % of all cases. Students who had been trained in how to apply the method gained higher inter-rater reliability than their colleagues with only a basic introduction. When repeating the test in June, participants were given improved instructions, organized in small subgroups, which resulted in a higher inter-rater reliability among farmers. In 66 % of all classifications, farmers assigned the same wetness class and the spread of class assignments was smaller. This study demonstrates that a wetness classification scheme based on qualitative indicators is a robust tool and can be applied successfully regardless of experience in crop growing and education level when an in-depth introduction and training is provided. The use of a simple and clear layout of the assessment form is important for reliable wetness class assignments.

Description: Nonlinear effects of locally heterogeneous hydraulic conductivity fields on regional stream–aquifer exchanges Hydrology and Earth System Sciences Discussions, 12, 7727-7764, 2015 Author(s): J. Zhu, C. L. Winter, and Z. Wang Computational experiments are performed to evaluate the effects of locally heterogeneous conductivity fields on regional exchanges of water between stream and aquifer systems in the Middle Heihe River Basin (MHRB) of northwestern China. The effects are found to be nonlinear in the sense that simulated discharges from aquifers to streams are systematically lower than discharges produced by a base model parameterized with relatively coarse effective conductivity. A similar, but weaker, effect is observed for stream leakage. The study is organized around three hypotheses: (H1) small-scale spatial variations of conductivity significantly affect regional exchanges of water between streams and aquifers in river basins, (H2) aggregating small-scale heterogeneities into regional effective parameters systematically biases estimates of stream–aquifer exchanges, and (H3) the biases result from slow-paths in groundwater flow that emerge due to small-scale heterogeneities. The hypotheses are evaluated by comparing stream–aquifer fluxes produced by the base model to fluxes simulated using realizations of the MHRB characterized by local (grid-scale) heterogeneity. Levels of local heterogeneity are manipulated as control variables by adjusting coefficients of variation. All models are implemented using the MODFLOW simulation environment, and the PEST tool is used to calibrate effective conductivities defined over 16 zones within the MHRB. The effective parameters are also used as expected values to develop log-normally distributed conductivity ( K ) fields on local grid scales. Stream-aquifer exchanges are simulated with K fields at both scales and then compared. Results show that the effects of small-scale heterogeneities significantly influence exchanges with simulations based on local-scale heterogeneities always producing discharges that are less than those produced by the base model. Although aquifer heterogeneities are uncorrelated at local scales, they appear to induce coherent slow-paths in groundwater fluxes that in turn reduce aquifer–stream exchanges. Since surface water–groundwater exchanges are critical hydrologic processes in basin-scale water budgets, these results also have implications for water resources management.

Description: Near–surface air temperature and snow skin temperature comparison from CREST-SAFE station data with MODIS land surface temperature data Hydrology and Earth System Sciences Discussions, 12, 7665-7687, 2015 Author(s): C. L. Pérez Díaz, T. Lakhankar, P. Romanov, J. Muñoz, R. Khanbilvardi, and Y. Yu Land Surface Temperature (LST) is a key variable (commonly studied to understand the hydrological cycle) that helps drive the energy balance and water exchange between the Earth's surface and its atmosphere. One observable constituent of much importance in the land surface water balance model is snow. Snow cover plays a critical role in the regional to global scale hydrological cycle because rain-on-snow with warm air temperatures accelerates rapid snow-melt, which is responsible for the majority of the spring floods. Accurate information on near-surface air temperature ( T -air) and snow skin temperature ( T -skin) helps us comprehend the energy and water balances in the Earth's hydrological cycle. T -skin is critical in estimating latent and sensible heat fluxes over snow covered areas because incoming and outgoing radiation fluxes from the snow mass and the air temperature above make it different from the average snowpack temperature. This study investigates the correlation between MODerate resolution Imaging Spectroradiometer (MODIS) LST data and observed T -air and T -skin data from NOAA-CREST-Snow Analysis and Field Experiment (CREST-SAFE) for the winters of 2013 and 2014. LST satellite validation is imperative because high-latitude regions are significantly affected by climate warming and there is a need to aid existing meteorological station networks with the spatially continuous measurements provided by satellites. Results indicate that near-surface air temperature correlates better than snow skin temperature with MODIS LST data. Additional findings show that there is a negative trend demonstrating that the air minus snow skin temperature difference is inversely proportional to cloud cover. To a lesser extent, it will be examined whether the surface properties at the site are representative for the LST properties within the instrument field of view.

Description: Postdisaster search and rescue is an important application of ultrawideband (UWB) radar systems, which mainly detect trapped victims by their respiratory-motion response. The development of a respiratory-motion detection (RMD) algorithm that can eliminate nonstationary clutter and noise is a challenging task for the application. A new algorithm is proposed to deal with the task in this letter. It uses the multichannel singular spectrum analysis (MSSA) technique to reconstruct the respiratory-motion response detected by a UWB radar. During the reconstruction, the periodicity and range interrelation characteristics of the response are exploited to adaptively identify signal subspaces. The performance of the algorithm is verified both by simulated and real data. The results show its improved performance over the reference algorithms, e.g., a singular-value-decomposition-based algorithm. The adaptive-MSSA-based RMD algorithm has great promise not only in practical use but also for future research of UWB-radar-based human being remote sensing.

Description: Land–sea contamination observed in Soil Moisture and Ocean Salinity (SMOS) brightness temperature images is found to have two main contributions: the floor error inherent of image reconstruction and a multiplicative error either in the antenna temperature or in the visibility samples measured by the correlator. The origin of this last one is traced down to SMOS calibration parameters to yield a simple correction scheme, which is validated against several geophysical scenarios. Autoconsistency rules in interferometric synthesis together with redundant and complementary calibration procedures provide a robust SMOS calibration scheme.

Description: This letter proposes a novel algorithm, which is based on the generalized method of moments (GMM), for the estimation and correction of phase errors induced in synthetic aperture radar (SAR) imagery. The GMM algorithm is used to replace the original phase-estimation kernel in the basic structure of the phase-gradient-autofocus algorithm. Since this novel algorithm does not require the observed signal to be a certain distribution model, it is able to estimate arbitrary phase errors. The GMM algorithm has the ability of estimating range-dependent phase errors, which makes it an efficient estimator. As a result, higher accuracy of the estimated phase errors and a better focused image can be achieved. Excellent results have been obtained in autofocusing and imaging experiments on real SAR data.

Description: Strong clutter reflections of terrain and marine surfaces obscure the contrast between the target-of-interest and clutter (terrain and marine surface reflections) in synthetic aperture radar (SAR) images and consequently hinder the efficiency of image interpretation and analysis. To overcome this problem, this letter proposes an efficient clutter suppression method in SAR images, which is named shedding irrelevant patterns (SIP). The essence is to construct a regression function that can suppress clutter and preserve the target patterns concurrently. We assume that the clutter is irrelevant to the target-of-interest and distinguishable in patterns in terms of image-pixel distribution and intensity (spatial information). Experimental results show the efficiency of the proposed method in both clutter suppression and target pattern preservation.

Description: A concern in hyperspectral image classification is the high number of required training samples. When traditional classifiers are applied, feature reduction (FR) techniques are the most common approaches to deal with this problem. Subspace-based classifiers, which are developed based on high-dimensional space characteristics, are another way to handle the high dimension of hyperspectral images. In this letter, a novel subspace-based classification approach is proposed and compared with basic and improved subspace-based classifiers. The proposed classifier is also compared with traditional classifiers that are accompanied by an FR technique and the well-known support vector machine classifier. Experimental results prove the efficiency of the proposed method, especially when a limited number of training samples are available. Furthermore, the proposed method has a very high level of automation and simplicity, as it has no parameters to be set.

Description: To improve the spatial density of measurement points of persistent-scatterer interferometry, distributed scatterer (DS) should be considered and processed. An important procedure in DS interferometry is the phase triangulation (PT). This letter introduces two modified PT algorithms (i.e., equal-weighted PT and coherence-weighted PT) and analyzes the mathematical relations between different published PT methods (i.e., the maximum-likelihood phase estimator, least squares estimator, and eigendecomposition-based phase estimators). The analysis shows that the above five PT methods share very similar mathematical forms with different weight values in the estimation procedure.

Description: Sparse representation-based classifier and its variants have been widely adopted for hyperspectral image (HSI) classification recently. However, sparse representation is unstable so that similar features might obtain significantly different sparse codes. Despite the instability, we find that the sparse codes follow a class-dependent distribution under the structured dictionary consisting of training samples from all classes. Based on this observation, a novel discriminative feature, sparse code histogram (SCH), is developed for HSI classification. By counting the SCH of each sample from the sparse codes of its spatial neighbors, we can statistically obtain the distribution pattern of sparse codes of the class to which the sample belongs, and then treat the SCH as a new feature for classification. To reduce the possible outliers among the neighbors, a shape-adaptive neighborhood extractor is also employed to enhance the stability of the histogram feature. Experimental results demonstrate that SCH enjoys a strong discriminative power, which can achieve notably better performance than several state-of-the-art methods for HSI classification with limited training samples.

Description: Does drought alter hydrological functions in forest soils? An infiltration experiment Hydrology and Earth System Sciences Discussions, 12, 7689-7725, 2015 Author(s): K. F. Gimbel, H. Puhlmann, and M. Weiler The water cycle is expected to change in future and severely affect precipitation patterns across central Europe and in other parts of the world, leading to more frequent and severe droughts. Usually, it is assumed that system properties, like soil properties, remain stable and will not be affected by drought events. To study if this assumption is appropriate, we address the effects of drought on the infiltration behavior of forest soils using dye tracer experiments on six sites in three regions across Germany, which were forced into drought conditions. The sites cover clayey, loamy and sandy textured soils. In each region, we compared a deciduous and a coniferous forest stand to address differences between the main tree species. The results of the dye tracer experiments show clear evidence for changes in infiltration behavior at the sites. The infiltration changed at the clayey plots from regular and homogeneous flow to fast preferential flow. Similar behavior was observed at the loamy plots, where large areas in the upper layers remained dry, displaying signs of strong water repellency. This was confirmed by WDPT tests, which revealed, in all except one plot, moderate to severe water repellency. Water repellency was also accountable for the change of regular infiltration to fingered flow in the sandy soils. The results of this study suggest that the "drought-history" or generally the climatic conditions in the past of a soil are more important than the actual antecedent soil moisture status regarding hydrophobicity and infiltration behavior; and also, that drought effects on infiltration need to be considered in hydrological models to obtain realistic predictions concerning water quality and quantity in runoff and groundwater recharge.

Description: Several detection statistics have been proposed for detecting fine ground disturbances between two synthetic aperture radar (SAR) images, such as vehicle tracks. The standard method involves estimating a local correlation coefficient between images. Other methods have been proposed using various statistical hypothesis tests. One of these alternative methods is a generalized likelihood ratio test (GLRT), which compares a full-correlation image model to a no-correlation image model. In this letter, we expand the GLRT to polarimetric SAR data and derive the appropriate GLRT detection statistics. Additionally, we explore relaxing the equal variance/equal polarimetric covariance assumptions used in previous results and find improved performance on macroscopic scene changes.

Description: For downward-looking linear array 3-D synthetic aperture radar, the resolution in cross-track direction is much lower than the ones in range and azimuth. Hence, superresolution reconstruction algorithms are desired. Since the cross-track signal to be reconstructed is sparse in the object domain, compressive sensing algorithm has been used. However, the imaging processing on the 3-D scene brings large computational loads, which renders challenges in both data acquisition and processing. To cover this shortage, truncated singular value decomposition is utilized to reconstruct a reduced-redundancy spatial measurement matrix. The proposed algorithm provides advantages in terms of computational time while maintaining the quality of the scene reconstructions. Moreover, our results on uniform linear array are generally applicable to sparse nonuniform linear array. Superresolution properties and reconstruction accuracies are demonstrated using simulations under the noise and clutter scenarios.

Description: This letter proposes a signal processing method of passive bistatic radar (PBR) exploiting an uncooperative radar as an illuminator. Compared with other opportunity illuminators, the transmitting signal of a radar usually has a better ambiguity function, which leads to a higher range resolution. Two channels are needed in PBR system. The reference channel is used to estimate radar signal parameters and reconstruct directly propagated signal. The surveillance channel is used to receive scattered wave. An array antenna and a simultaneous multibeam algorithm are necessary in the surveillance channel due to the flexible beam scanning of the uncooperative radar. The procedure of the proposed method is explained in detail, which is then followed by a field experiment. Preliminary results from the field experiment show that the proposed method can be applied to target angle and bistatic range measurement successfully.

Description: In order to achieve 3-D imaging with an airborne down-looking linear-array synthetic aperture radar (LASAR), a uniform virtual antenna array may be obtained by aperture synthesis of the cross-track sparse multiple-input–multiple-output array. However, the actual 3-D imaging quality is unavoidably degraded by errors in the virtual element position. In this letter, we investigate the effects of these errors on the forms and the degrees of image quality degradation by decomposing the error-related stochastic processes via an orthogonal transform based on discrete Legendre polynomials. It should be noted that these analyses are helpful for designing a LASAR system and providing a reference for specifying the requisite precision of measurement devices and calibration methods. Finally, we briefly consider the use of calibration methods to eliminate the effects of errors.

Description: In problems where labeled data are scarce, semisupervised learning (SSL) techniques are an attractive framework that can exploit both labeled and unlabeled data. These approaches typically rely on a smoothness assumption such that examples that are similar in input space should also be similar in label space. In many domains, such as remotely sensed hyperspectral image (HSI) classification, the data violate this assumption. In response, we propose a general method by which a neighborhood graph used in SSL is learned using binary classifiers that are trained to predict whether a pair of pixels shares the same label. Working within the framework of semisupervised neural networks (SSNNs), we show that our approach improves on the performance of the SSNN on two HSI data sets.

Description: In this letter, new models for the spatial correlation of sea clutter texture and intensity are proposed as improved versions of current power law models or exponential decay model. The models for texture have three unknown parameters, and thus can be called triparametric models. The structure of the models is a weighted sum of two components, which can describe the decaying process of the correlation coefficient with spatial lags, as well as the periodic behavior due to the existence of transient coherent structures in sea clutter. Unknown parameters are optimized by the nonlinear least square fit method. Models for sea clutter intensity can be obtained through a linear transform for uncorrelated speckle based on the compound-Gaussian representation of sea clutter. The proposed models are validated and compared with current models using S- and C-band measured sea clutter data. Analysis results indicate the effectiveness of the proposed models in that they can describe the behavior of spatial correlation coefficients with higher accuracy.

Description: Terrestrial laser scanning (TLS) has become a popular tool for acquiring source data points which can be used to construct digital elevation models (DEMs) for a wide number of applications. A TLS point cloud often has a very fine spatial resolution, which can represent well the spatial variation of a terrain surface. However, the uncertainty in DEMs created from this relatively new type of source data is not well understood, which forms the focus of this letter. TLS survey data representing four terrain surfaces of different characteristics were used to explore the effects of surface complexity and typical TLS data density (in terms of data point spacing) on DEM accuracy. The spatial variation in TLS data can be decomposed into parts corresponding to the signal of spatial variation (of terrain surfaces) and noise due to measurement error. We found a linear relation between the DEM error and the typical TLS data spacings considered (30–100 mm) which arises as a function of the interpolation error, and a constant contribution from the propagated data noise. This letter quantifies these components for each of the four surfaces considered and shows that, for the interpolation method considered here, higher density sampling would not be beneficial.

Description: A novel way to estimate the live fuel moisture content (LFMC) was explored from the ratio of canopy water content (CWC) and foliage dry biomass (FDB). The CWC was estimated using the PROSAIL (PROSPECT + SAIL) radiative transfer model from the Landsat 8 product. A weak constraint 4-D variational data assimilation method was employed to assimilate the temporally estimated leaf area index into a soil-water-atmosphere-plant (SWAP) model for optimizing the model control variables. Then, the SWAP model was reinitialized with this optimum set of control variables, and better prediction of FDB was obtained. Results showed that a high accuracy level was achieved for the estimated CWC ( $R^{2}=0.91$ , $mbox{RMSE}=84.74 mbox{g/m}^2$ ) and FDB ( $R^2=0.88$ , $mbox{RMSE}=48.54 mbox{g/m} ^2$ ) when compared with in situ measured values. However, the accuracy level of estimated LFMC was poor ( $R^2=0.59$ , $mbox{RMSE} =30.85%$ ) . Further analyses find that the estimated LFMC is reliable for low LFMC but challenged for high LFMC, which indicates that the presented method still makes sense to the assessment of wildfire risk since the wildfire generally occurs when the vegetation is in low LFMC condition.

Description: In this letter, we present the use of experimental human micro-Doppler signature data gathered by a multistatic radar system to discriminate between unarmed and potentially armed personnel walking along different trajectories. Different ways of extracting suitable features from the spectrograms of the micro-Doppler signatures are discussed, particularly empirical features such as Doppler bandwidth, periodicity, and others, and features extracted from singular value decomposition (SVD) vectors. High classification accuracy of armed versus unarmed personnel (between 90% and 97% depending on the walking trajectory of the people) can be achieved with a single SVD-based feature, in comparison with using four empirical features. The impact on classification performance of different aspect angles and the benefit of combining multistatic information is also evaluated in this letter.

Description: Automatic urban area detection in remote sensing images is an important application in the field of earth observation. Most of the existing methods employ feature classifiers and thereby contain a data training process. Moreover, some methods cannot detect urban areas in complex scenes accurately. This letter proposes an automatic urban area detection method that uses multiple features that have different resolutions. First, a downsampled low-resolution image is used to segment the candidate area. After the corner points of the urban area are extracted, a weighted Gaussian voting matrix technique is employed to integrate the corner points into the candidate area. Then, the edge features and homogeneous region are extracted by using the original high-resolution image. Using these results as the input, the processes of guided filtering and contrast enhancement can finally detect accurately the urban areas. This method combines multiple features, such as corner, edge, and regional characteristics, to detect the urban areas. The experimental results show that the proposed method has better detection accuracy for urban areas than the existing algorithms.

Description: In marine sciences, time series are often nonlinear and nonstationary. Adequate and specific methods are needed to analyze such series. In this letter, an application of the empirical mode decomposition method (EMD) associated to the Hilbert spectral analysis (HSA) is presented. Furthermore, EMD-based time-dependent intrinsic correlation (TDIC) analysis is applied to consider the correlation between two nonstationary time series. Four temperature time series obtained from automatic measurements in nearshore waters of the Réunion island are considered, recorded every 10 min from July 2011 to January 2012. The application of the EMD on these series and the estimation of their power spectra using the HSA are illustrated. The authors identify low-frequency tidal waves and display the pattern of correlations at different scales and different locations. By TDIC analysis, it was concluded that the high-frequency modes have small correlation, whereas the trends are perfectly correlated.

Description: Detecting vehicles in aerial images provides important information for traffic management and urban planning. Detecting the cars in the images is challenging due to the relatively small size of the target objects and the complex background in man-made areas. It is particularly challenging if the goal is near-real-time detection, i.e., within few seconds, on large images without any additional information, e.g., road database and accurate target size. We present a method that can detect the vehicles on a 21-MPixel original frame image without accurate scale information within seconds on a laptop single threaded. In addition to the bounding box of the vehicles, we extract also orientation and type (car/truck) information. First, we apply a fast binary detector using integral channel features in a soft-cascade structure. In the next step, we apply a multiclass classifier on the output of the binary detector, which gives the orientation and type of the vehicles. We evaluate our method on a challenging data set of original aerial images over Munich and a data set captured from an unmanned aerial vehicle (UAV).

Description: We compare five slope correction methods developed by Walter et al. , Montes et al. , Schleppi et al. , España et al. , and Gonsamo et al. (referred to as WAL, MON, SCH, ESP, and GON, respectively) using artificial fisheye pictures simulated by graphics software and a lookup table (LUT) retrieval method. The LUT is built by simulating the directional gap fraction as a function of leaf area index (LAI) and average leaf inclination angle (ALIA) using the Poisson law. LAI and ALIA estimates correspond to the case of the LUT that provides the lowest root-mean-square error between the observed gap fractions after slope correction and the simulated ones. Three LAI values (1.5, 3.5, and 5.5), four ALIA values (26.8°, 45°, 57.5°, and 63.2°), and three slope angles (0°, 20°, and 50°) constituted 36 samples of random scenes. ESP is recommended because its results are accurate and independent on the leaf angle distribution (LAD), while GON only performs well for spherical LAD. The three other methods present less good performances with underestimation or overestimation of LAI and/or ALIA depending on the LAD, and the recommended order for them is MON, SCH, and WAL.

Description: In this letter, an improved phase correlation (PC) method based on 2-D plane fitting and the maximum kernel density estimator (MKDE) is proposed, which combines the idea of Stone's method and robust estimator MKDE. The proposed PC method first utilizes a vector filter to minimize the noise errors of the phase angle matrix and then unwraps the filtered phase angle matrix by the use of the minimum cost network flow unwrapping algorithm. Afterward, the unwrapped phase angle matrix is robustly fitted via MKDE, and the slope coefficients of the 2-D plane indicate the subpixel shifts between images. The experiments revealed that the improved method can effectively avoid the impact of outliers on the phase angle matrix during the plane fitting and is robust to aliasing and noise. The matching accuracy can reach 1/50th of a pixel using simulated data. The real image sequence tracking experiment was also undertaken to demonstrate the effectiveness of the proposed PC method with a registration accuracy of root-mean-square error better than 0.1 pixels.

Description: Accurately mapping forest carbon density by combining sample plots and remotely sensed images has become popular because this method provides spatially explicit estimates. However, mixed pixels often impede the improvement of the estimation. In this letter, regression modeling and spectral unmixing analysis were integrated to improve the estimation of forest carbon density for the You County of Hunan, China, using Landsat Thematic Mapper images. Linear spectral unmixing with and without a constraint (LSUWC and LSUWOC) and nonlinear spectral unmixing (NSU) were compared to derive the fractions of five endmembers, particularly forests. Stepwise regression, logistic regression, and polynomial regression (PR) with and without the forest fraction used as an independent variable and the product of the forest fraction image and the map from the best model without the forest fraction were compared. The models were developed using 56 sample plots, and their results were validated using 26 test plots. The decomposition of mixed pixels was assessed using higher spatial resolution SPOT images and a corresponding land cover map. The results showed that 1) LSUWC more accurately estimated the endmember fractions than LSUWOC and NSU, 2) PR had the greatest estimation accuracy of forest carbon, and 3) combining regression modeling and spectral unmixing increased the estimation accuracy by 31%–39%, and introducing the forest fraction into the regressions performed better than the product of forest fraction image and the results from PR without the fraction. This implied that the integrations provided great potential in reducing the impacts of mixed pixels in mapping forest carbon.

Description: This letter proposes a multiresolution technique to address the high computational cost in remote sensing image registration. The scale-invariant feature transform is applied to detect keypoints and descriptors, and then, global information combined with descriptors is utilized to establish keypoint mappings. Keypoints are first classified according to their octaves. Then, in the lowest resolution, the keypoints of the largest octave are mapped with descriptors and the global information, giving an initial affine transformation $T_0$ . In the next octave, the keypoints of the second largest octave are mapped by employing $T_0$ to narrow the space of matching keypoints. By this means, the process of establishing keypoint correspondences is conducted from one resolution (octave) to the next as the obtained transformation gets finer until we get to the highest resolution. Due to the high computational expense of computing global information, the proposed technique is important for aligning large-size remote sensing imagery. Experimental results show that the proposed method can achieve a comparable registration accuracy but with a less computational cost than directly building keypoint mappings on images of large size.

Description: Gold nanorods (GNRs) with surface plasmon resonance peak at 1063.8 nm were fabricated and experimentally exploited as the single and combined saturable absorber (SA) in a $Q$ -switched Nd:YAG laser for the first time. In the situation using GNRs as a single SA, the maximum pulse energy of 19 $muhbox{J} $ was achieved at a pulse repetition rate of 20 kHz. However, due to the small effective modulation depth of the GNRs, microsecond pulses were generated in such a GNRs-based passively $Q$ -switched laser. A novel configuration using GNRs as a combined SA in an acoustooptic $Q$ -switched Nd:YAG laser was proposed. The small nonlinear loss modulation induced by the GNRs played a distinctive role in flexibly manipulating the nanosecond pulse waveforms. Our study demonstrated the feasibility of generating high-energy pulses by singly using GNRs-based SA in all solid-state lasers. Meanwhile, a new application mode incorporating GNRs as a combined SA was studied, providing an alternative method to take advantage of the nanomaterial-based SA with a relatively small modulation depth.

Description: Presents corrections to the paper, "Estimation of forest biomass from two-level model inversion of single-pass InSAR data" (Soja, M.J., et al.,Trans. Geosci. Remote Sens., vol. 53, no. 9, pp. 5083???5099, Sep. 2015).

Description: In this letter, we present an efficient parallel implementation of composite kernels in support vector machines (SVMs) for hyperspectral image (HSI) classification. Our implementation makes effective use of commodity graphics processing units (GPUs). Specifically, we port the calculation of composite kernels to GPUs, perform intensive computations based on NVidia's compute unified device architecture, and execute the rest of the operations related with control and small data calculations in the CPU. Our experimental results, conducted using real hyperspectral data sets and NVidia GPU platforms, indicate significant improvements in terms of computational effectiveness, achieving near-real-time performance of spatial–spectral HSI classification for the first time in the literature.

Description: Compared with airborne laser scanning, terrestrial laser scanning (TLS) offers ground-based point cloud data of trees and provides greater potential to accurately estimate tree and stand parameters. However, there is a lack of effective methods to accurately identify locations of individual trees from TLS point cloud data. It is also unknown whether the estimation accuracy of the parameters, including tree height (H), diameter at breast height (DBH), and so on, using TLS can meet the requirement of forest management and planning. In this letter, a novel method to effectively process point cloud data and further determine the locations of individual trees in a stand based on the central coordinates of point cloud data on a defined grid according to the largest DBH was developed. Moreover, a point-cloud-data-based convex hull algorithm and the cylinder method were, respectively, used to estimate DBH and H of individual trees. This study was conducted in a pure Chinese fir plantation of 45 trees located in Huang-Feng-Qiao forest farm, You County of Hunan, China. The comparison of the estimated and observed values showed that the obtained tree locations had errors of less than 20 cm, and the relative root mean square errors for the estimates of both DBH and H were less than 5%. This implies that TLS is very promising for the retrieval of tree and stand parameters in forest stands. For the applications of these methods to mixed forests with a structure of multilayer canopies, further examination is needed.

Description: Commercial, industrial, and military aerospace designs are increasingly deploying MEMS micro inertial measurement unit (MIMU) for motion control, automation, and positioning applications, such as the unmanned aerial vehicle (UAV), robot, and smart phone. On the one hand, MIMU has the merit of low cost, small size, low-power consumption, and high shock resistance, but on the other hand, low-cost MIMU is affected by systematic error caused by the instability of the drift, scaling factors, and axes misalignment, which may lead to large errors in the position and attitude’s determination from time to time. That means calibration before use is an effective way to improve the practical precision of MIMU. However, many customers have no precise turntable to calibrate the MIMU before they use it. To address these problems, this paper presents an easy self-calibration method to implement calibration of the MIMU on a common table only with an inclined surface, no precise turntable is needed. The calibration method is based on the following principles. First, the module of the output vector of the orthogonal configured three-axis accelerometers is equal to unit gravity. Second, when IMU rotates to a known gesture with a stable axis, the angles can be calculated through integration. Third, when the accelerometers’ parameters are calculated, it can act as a level datum. Furthermore, the accelerometers on the inclined surface are used to determine the rotating heading datum. Finally, after a series static positions test and rotating test, the parameters can be extracted and estimated. To demonstrate the success and the convenience of the proposed method, comparison experiments with the precision turntable have been made on an ADI’s MIMU. The calibration results show that the accuracy and precision of this method is quite equivalent with the turntable-based calibration, and the scale factors error with an order of magnitude always equal or less than $10^{-5}$ . The observed static and dynamic yaw maximum angular error in a certain period is <0.8°, the pitch maximum angular error is <0.5°, and the roll maximum angular error is <0.3°.

Description: Time–frequency packing (TFP) transmission provides the highest achievable spectral efficiency with a constrained symbol alphabet and detector complexity. In this paper, the application of the TFP technique to fiber-optic systems is investigated and experimentally demonstrated. The main theoretical aspects, design guidelines, and implementation issues are discussed, focusing on those aspects which are peculiar to TFP systems. In particular, adaptive compensation of propagation impairments, matched filtering, and maximum a posteriori probability detection are obtained by a combination of a two-dimensional equalizer and four eight-state parallel Bahl–Cocke–Jelinek–Raviv (BCJR) detectors. A novel algorithm that ensures adaptive equalization, channel estimation, and a proper distribution of tasks between the equalizer and BCJR detectors is proposed. A set of irregular low-density parity-check codes with different rates is designed to operate at low error rates and approach the spectral efficiency limit achievable by TFP at different signal-to-noise ratios. An experimental demonstration of the designed system is finally provided with five dual-polarization QPSK-modulated optical carriers, densely packed in a 100-GHz bandwidth, employing a recirculating loop to test the performance of the system at different transmission distances.

Description: We investigated nonlinear optical characteristics of Tungsten disulfide (WS 2 ) films and experimentally demonstrated their high potential for application as nonlinear saturable absorbers in passively mode-locked fiber lasers. Side polished fiber (SPF) was fabricated and WS 2 film was overlaid to provide an efficient evanescent field interaction. The WS 2 film was prepared using two methods: liquid phase exfoliation to form few-layer nano-sheets, and chemical vapor deposition (CVD) to grow uniform multilayer WS 2 on a SiO 2 substrate. Two SPF saturable absorbers were prepared by either spin coating WS 2 solution or lifting off the multilayer CVD WS 2 on SPF. An all-fiber ring cavity was built and the WS 2 film overlaid on SPF was employed as a mode locker along with Er-doped fiber as a gain medium. Using the spin-coated WS 2 SPF, stable soliton-like pulses were generated with a spectral width of 5.6 nm and 467 fs pulse duration. The fiber laser cavity containing CVD WS2 SPF generated a transform-limited soliton pulse train with a spectral width of 8.23 nm and a pulse duration of 332 fs. Our study confirmed a high potential of WS 2 film as a novel 2-D nonlinear optical material for laser applications.

Description: We propose and experimentally demonstrate an on-chip all-optical differential-equation solver capable of solving second-order ordinary differential equations (ODEs) characterizing continuous-time linear time-invariant systems. The photonic device is implemented by a self-coupled microresonator on a silicon-on-insulator platform with mutual coupling between the cavity modes. Owing to the mutual mode coupling within the same resonant cavity, the resonance wavelengths induced by different cavity modes are self-aligned, thus avoiding precise wavelength alignment and unequal thermal wavelength drifts as in the case of cascaded resonators. By changing the mutual mode coupling strength, the proposed device can be used to solve second-order ODEs with tunable coefficients. System demonstration using the fabricated device is carried out for 10-Gb/s optical Gaussian and super-Gaussian input pulses. The experimental results are in good agreement with theoretical predictions of the solutions, which verify the feasibility of the fabricated device as a tunable second-order photonic ODE solver.

Description: We demonstrate a 100 Gb/s short reach system using a multicarrier transmitter based on a gain switched monolithically integrated laser. An optical comb source with 12.5-GHz free spectral range is achieved by gain-switching an integrated passive feedback laser. The 100 Gb/s wavelength division multiplexed, single sideband, direct detection, orthogonal frequency division multiplexed (WDM-SSB-DD-OFDM) system operates over 25 km standard single mode fiber exhibiting a spectral efficiency of 1.8 b/s/Hz. Receiver sensitivity of –14.2 dBm is achieved after 25 km transmission. Performance optimization with phase and amplitude precompensation is employed to improve the SSB OFDM modulation thereby reducing the interchannel interference and overcoming the power fading induced by the optical filter. We also present a theoretical analysis of the SSB-OFDM modulation.

Description: In this paper, narrow-band emission lines are generated by means of two random distributed feedback fiber laser schemes. Spectral line-widths as narrow as 3.2 pm have been measured, which significantly improves previous reported results. The laser is analyzed with the aim of obtaining a spectral line-width as narrow as possible. Additionally a variation of this setup for multi-wavelength operation is also validated. Both schemes present a simple topology that use a combination of phase-shifted fiber Bragg gratings and regular fiber Bragg gratings as filtering elements.

Description: We report on the fabrication and characterization of 19-cell hypocycloid-shape Kagome fibers with core size larger than 100 μm. These inhibited coupling fibers present low propagation loss (100 dB/km) over broad transmission range with low chromatic dispersion combined with ultra-low power overlap with silica surround, making them an efficient solution for ultra-high power laser handling, ultra-fast laser delivery, and plasma photonics applications.

Description: The iterative demodulation and decoding algorithm introduced in 2005 by Colavolpe, Barbieri, and Caire to cope with channels affected by phase noise needs pilot symbols to bootstrap. However, pilot symbols reduce the spectral efficiency of the system and, consequently, system's throughput. The aim of this paper is to show that trellis-based demodulation can be used to bootstrap the iterative process without the need of pilot symbols. Also, the complexity issue of trellis-based demodulation is addressed in this paper. The result is that the performance of iterative demodulation and decoding after the iterations is virtually unaffected by complexity reduction, provided that the reduced-complexity demodulator guarantees cycle-slip-free operation. From the numerical results presented in this paper, we show that cycle-slip-free operation can be achieved with substantial complexity reduction also for phase noise associated with linewidths of practical interest.

Description: The optical pulse evolution in a highly nonlinear normal dispersion-increasing fiber has been considered, both experimentally and theoretically. It was found that large spectral broadening in tapered waveguides could occur without temporal instabilities and impose the linear frequency modulation, i.e., chirp, required for high-quality pulse compression. The pedestal-free pulses have been demonstrated after dechirping in a standard single-mode fiber.

Description: This work proposes a novel localized surface plasmon resonance (LSPR) biochemical sensor featuring high sensitivity and a high resolution. The sensor was divided into two subcomponents according to their distinct functions; namely, single-mode fiber and metal array. Single-mode fibers located on the left and right sides of the sensors function as the input and output for optical fiber signals. A metal array comprising an arrangement of cylindrical nanometal particles served as the detection area of the sensor. To effectively reduce the memory capacity and calculation time, two innovative techniques (i.e., object meshing and boundary meshing) were integrated with the finite element method. With the area of the triangular elements used as a basis, the object boundary, small object, medium object, and large objects were meshed at a ratio of 1:8:160:1600. The improved numerical simulation methods and six design procedures were adopted to develop and analyze the proposed LSPR biochemical sensor. The results show that the novel LSPR biochemical sensor outperformed two current high-performance biochemical sensors and provided additional advantages such as short length (approximately 430 μm), high resolution (approximately –120 dB), and high sensitivity (approximately 127 604 nm/RIU).

Description: We propose an ultra-broadband super light absorber by integrating different-sized tapered hyperbolic metamaterial (HMM) waveguides, each of which has a different and wide absorption band due to broadband slow-light response, into a unit cell. We numerically demonstrate that such an absorber is superior to a single-sized HMM absorber in terms of absorption bandwidth, while maintaining a comparable absorption efficiency. A three different-sized HMM absorber presents the capability of working with an ultra-wide frequency band ranging from 1 to 30 THz, which is much larger than previously proposed absorbers working in the same spectral region. Such a design shows great promise for a broad range of applications such as thermal emitters, photovoltaics, optical-chemical energy harvesting, and stealth technology, where ultra-wideband absorption is in very high demand.

Description: This paper proposes a ring-based integrated wireless optical network architecture and an associated protocol that involves Ethernet passive optical network (EPON) and long term evolution (LTE) wireless network. The architecture along with the proposed protocol is instrumental toward the reduction of handover delay. The proposed ring-based EPON architecture facilitates the implementation of the X2 interface for LTE network by enabling the optical network units of the EPON backhaul to directly communicate with each other. The work further discusses about an open access network architecture where a single EPON can be used by multiple mobile service providers without compromising information security. Extensive simulations have been carried out to evaluate the performance of the proposed network. An analytical model has been introduced to calculate the queuing delay experienced by the X2 interface. The model has been validated with the simulations results.

Description: High spectrum efficiency and fast restoration speed are highly desired for survivable elastic optical networks (EONs). In this paper, we take the advantages of failure-independent path-protecting preconfigured cycles (FIPP p -cycles) and investigate how to realize spectrum efficient resilience design with them. We first study the problem of offline service provisioning with FIPP p -cycles. We formulate an integer linear programming model and prove that the problem is $mathcal{NP}$ hard. Then, several time-efficient heuristics are designed for FIPP p -cycle formulation and related routing, modulation format, and spectrum assignment. Extensive simulations on offline provisioning verify that the heuristics can obtain near-optimal solutions. Next, we consider online service provisioning with FIPP p -cycles in dynamic EONs. In order to overcome the decrease of protection efficiency during dynamic network operation, we propose a p -cycle reconfiguration scheme to reoptimize protection structures on-the-fly. Simulation results demonstrate that the proposed algorithms can improve spectrum efficiency and reduce bandwidth blocking probability simultaneously.

Description: Photo-thermo-refractive (PTR) glasses are the class of polyfunctional materials that combine the properties of several monofunctional materials such as photo-refractive, laser, luminescent, and plasmonic ones. Based on PTR glasses, various diffractive holographic elements and photonic devices were developed in both the volume and fiber versions. In this paper, the fabrication of optical planar waveguides on PTR glass by low-temperature ion exchange is reported for the first time. Planar waveguides were fabricated through substituting the sodium ions in glass by silver, potassium, rubidium, and cesium ones from the nitrate melts. The silver waveguides were shown to have the largest depth (27 μm) and reveal no birefringence. For the silver waveguides, an increase in the refractive index is caused by differences in the polarizabilities of cations exchanged. The maximum increment of the refractive index was observed in the cesium waveguides (0.0512). An increase in the refractive index and also appearing the birefringence in potassium, rubidium, and cesium waveguides are found to be due to the compressive mechanical stresses and their relaxation. The potentialities of the ion exchange technology for fabricating, in PTR glasses, planar gradient waveguides with low losses (0.5 dB/cm) are demonstrated, the potentialities extending the application field of PTR glasses in photonics.

Description: State determination of catalytic converters based on an ultra-wideband communication system Journal of Sensors and Sensor Systems, 4, 255-262, 2015 Author(s): I. Motroniuk, R. Stöber, and G. Fischerauer A novel microwave-based approach for monitoring the state of aftertreatment systems such as diesel particulate filters (DPFs), three-way catalytic converters (TWCs), and selective catalytic reduction (SCR) catalysts is proposed. The volume inside the metallic housing of the DPF, TWC, or SCR is considered as a wireless communication channel between two terminals of a communication system. It is shown that, depending on the transmission channel characteristics, the properties of the catalyst, such as the catalyst state, can be inferred. This is done by means of an ultra-wideband (UWB) measurement and the subsequent evaluation and processing of the waveform in the time and frequency domains.

Description: A comparison of chemical mechanisms using tagged ozone production potential (TOPP) analysis Atmospheric Chemistry and Physics, 15, 8795-8808, 2015 Author(s): J. Coates and T. M. Butler Ground-level ozone is a secondary pollutant produced photochemically from reactions of NO x with peroxy radicals produced during volatile organic compound (VOC) degradation. Chemical transport models use simplified representations of this complex gas-phase chemistry to predict O 3 levels and inform emission control strategies. Accurate representation of O 3 production chemistry is vital for effective prediction. In this study, VOC degradation chemistry in simplified mechanisms is compared to that in the near-explicit Master Chemical Mechanism (MCM) using a box model and by "tagging" all organic degradation products over multi-day runs, thus calculating the tagged ozone production potential (TOPP) for a selection of VOCs representative of urban air masses. Simplified mechanisms that aggregate VOC degradation products instead of aggregating emitted VOCs produce comparable amounts of O 3 from VOC degradation to the MCM. First-day TOPP values are similar across mechanisms for most VOCs, with larger discrepancies arising over the course of the model run. Aromatic and unsaturated aliphatic VOCs have the largest inter-mechanism differences on the first day, while alkanes show largest differences on the second day. Simplified mechanisms break VOCs down into smaller-sized degradation products on the first day faster than the MCM, impacting the total amount of O 3 produced on subsequent days due to secondary chemistry.

Description: Size distribution and optical properties of mineral dust aerosols transported in the western Mediterranean Atmospheric Chemistry and Physics Discussions, 15, 21607-21669, 2015 Author(s): C. Denjean, F. Cassola, A. Mazzino, S. Triquet, S. Chevaillier, N. Grand, T. Bourrianne, G. Momboisse, K. Sellegri, A. Schwarzenbock, E. Freney, M. Mallet, and P. Formenti This study presents in situ aircraft measurements of Saharan mineral dust transported over the western Mediterranean basin in June–July 2013 during the ChArMEx/ADRIMED (the Chemistry-Aerosol Mediterranean Experiment/Aerosol Direct Radiative Impact on the regional climate in the MEDiterranean region) airborne campaign. Dust events differing in terms of source region (Algeria, Tunisia and Morocco), time of tranport (1–5 days) and height of transport were sampled. Mineral dust were transported above the marine boundary layer, which conversely was dominated by pollution and marine aerosols. The dust vertical structure was extremely variable and characterized by either a single layer or a more complex and stratified structure with layers originating from different source regions. Mixing of mineral dust with pollution particles was observed depending on the height of transport of the dust layers. Dust layers carried higher concentration of pollution particles at intermediate altitude (1–3 km) than at elevated altitude (〉 3 km), resulting in scattering Angstrom exponent up to 2.2 within the intermediate altitude. However, the optical properties of the dust plumes remained practically unchanged with respect to values previously measured over source regions, regardless of the altitude. Moderate light absorption of the dust plumes was observed with values of aerosol single scattering albedo at 530 nm ranging from 0.90 to 1.00 ± 0.04. Concurrent calculations from the aerosol chemical composition revealed a negligible contribution of pollution particles to the absorption properties of the dust plumes that was due to a low contribution of refractory black carbon in regards to the fraction of dust and sulfate particles. This suggests that, even in the presence of moderate pollution, likely a persistent feature in the Mediterranean, the optical properties of the dust plumes could be assimilated to those of native dust in radiative transfer simulations, modeling studies and satellite retrievals over the Mediterranean. Measurements also showed that the coarse mode of mineral dust was conserved even after 5 days of transport in the Mediterranean, which contrasts with the gravitational depletion of large particles observed during the transport of dust plumes over the Atlantic. Simulations with the WRF mesoscale meteorological model highlighted a strong vertical turbulence within the dust layers that could prevent deposition of large particles during their atmospheric transport. This has important implications for the dust radiative effects due to surface dimming, atmospheric heating and cloud formation. The results presented here add to the observational dataset necessary for evaluating the role of mineral dust on the regional climate and rainfall patterns in the western Mediterranean basin.

Description: Identification of particulate organosulfates in three megacities at the middle and lower reaches of the Yangtze River Atmospheric Chemistry and Physics Discussions, 15, 21415-21448, 2015 Author(s): X. K. Wang, S. Rossignol, Y. Ma, L. Yao, M. Y. Wang, J. M. Chen, C. George, and L. Wang PM 2.5 filter samples have been collected in three megacities i.e., Wuhan (WH), Nanjing (NJ), and Shanghai (SH) at the middle and lower reaches of the Yangtze River, respectively. Analysis of those samples using an ultra-high performance liquid chromatography (UHPLC) coupled to an orbitrap mass spectrometer (MS) allowed detection of about two hundred particulate organosulfates (OSs), including dozens of nitrooxy-organosulfates, at each location. While aliphatic OSs represented more than 78 % of the detected OSs at the three locations, aromatic OSs were much less abundant. OSs with two to four isomers accounted for about 50 % of the total OSs on average in these megacities, and the percentage of OSs with six and more isomers in WH was more significant than those in SH and NJ. The average molecular weight, and the degrees of oxidation and saturation of OSs in the WH summer samples were greater than those in WH winter samples. In SH, the average molecular weight and the degree of oxidation of OSs in summer samples were greater than those in winter samples, but the degree of saturation was similar between the two seasons. In summer, the average molecular weight, and the degrees of oxidation and unsaturation of OSs were smallest in WH among the three cities. Between NJ and SH, the average molecular weight and the degree of saturation of OSs were close and the degree of oxidation of OSs in NJ was smaller. Kendrick mass defect diagrams and Van Krevelen diagrams indicated that the characteristics of identified OSs between in NJ and in SH shared better similarity. In addition, the identity and abundance of OSs in SH showed clear seasonal and diurnal variations. OSs in summer were more abundant than they were in winter due to stronger photochemical reactions in summer. The relative abundance of OSs at night was greater than that in the daytime and more nitrooxy-OSs existed at night, probably because of active NO 3 radical chemistry at night. In SH summer samples, OSs with 5 and 10 carbons (C 5 and C 10 ) were the most abundant, indicating the importance of isoprene and monoterpenes as precursors of OSs, whereas the relative abundances of OSs with 8, 9, and more than 14 carbons (C 8 , C 9 , and C 14+ ) were also high in SH winter samples, urging the need to further understand the precursors of OSs.

Description: Trace metal characterization of aerosol particles and cloud water during HCCT 2010 Atmospheric Chemistry and Physics, 15, 8751-8765, 2015 Author(s): K. W. Fomba, D. van Pinxteren, K. Müller, Y. Iinuma, T. Lee, J. L. Collett Jr., and H. Herrmann Trace metal characterization of bulk and size-resolved aerosol and cloud water samples were performed during the Hill Cap Cloud Thuringia (HCCT) campaign. Cloud water was collected at the top of Mt. Schmücke while aerosol samples were collected at two stations upwind and downwind of Mt. Schmücke. Fourteen trace metals including Ti, V, Fe, Mn, Co, Zn, Ni, Cu, As, Sr, Rb, Pb, Cr, and Se were investigated during four full cloud events (FCEs) that fulfilled the conditions of a continuous air mass flow through the three stations. Aerosol particle trace metal concentrations were found to be lower than those observed in the same region during previous field experiments but were within a similar range to those observed in other rural regions in Europe. Fe and Zn were the most abundant elements with concentration ranges of 0.2–111.6 and 1.1–32.1 ng m −3 , respectively. Fe, Mn, and Ti were mainly found in coarse mode aerosols while Zn, Pb, and As were mostly found in the fine mode. Correlation and enrichment factor analysis of trace metals revealed that trace metals such as Ti and Rb were mostly of crustal origin while trace metals such as Zn, Pb, As, Cr, Ni, V, and Cu were of anthropogenic origin. Trace metals such as Fe and Mn were of mixed origins including crustal and combustion sources. Trace metal cloud water concentration decreased from Ti, Mn, Cr, to Co with average concentrations of 9.18, 5.59, 5.54, and 0.46 μg L −1 , respectively. A non-uniform distribution of soluble Fe, Cu, and Mn was observed across the cloud drop sizes. Soluble Fe and Cu were found mainly in cloud droplets with diameters between 16 and 22 μm, while Mn was found mostly in larger drops greater than 22 μm. Fe(III) was the main form of soluble Fe especially in the small and larger drops with concentrations ranging from 2.2 to 37.1 μg L −1 . In contrast to other studies, Fe(II) was observed mainly in the evening hours, implying its presence was not directly related to photochemical processes. Aerosol–cloud interaction did not lead to a marked increase in soluble trace metal concentrations; rather it led to differences in the chemical composition of the aerosol due to preferential loss of aerosol particles through physical processes including cloud drop deposition to vegetative surfaces.

Description: Known and unknown unknowns: the application of ensemble techniques to uncertainty estimation in satellite remote sensing data Atmospheric Measurement Techniques Discussions, 8, 8509-8562, 2015 Author(s): A. C. Povey and R. G. Grainger This paper discusses a best-practice representation of uncertainty in satellite remote sensing data. An estimate of uncertainty is necessary to make appropriate use of the information conveyed by a measurement. Traditional error propagation quantifies the uncertainty in a measurement due to well-understood perturbations in a measurement and auxiliary data – known, quantified "unknowns". The underconstrained nature of most satellite remote sensing observations requires the use of various approximations and assumptions that produce non-linear systematic errors that are not readily assessed – known, unquantifiable "unknowns". Additional errors result from the inability to resolve all scales of variation in the measured quantity – unknown "unknowns". The latter two categories of error are dominant in underconstrained remote sensing retrievals and the difficulty of their quantification limits the utility of existing uncertainty estimates, degrading confidence in such data. This paper proposes the use of ensemble techniques to present multiple self-consistent realisations of a data set as a means of depicting unquantified uncertainties. These are generated using various systems (different algorithms or forward models) believed to be appropriate to the conditions observed. Benefiting from the experience of the climate modelling community, an ensemble provides a user with a more complete representation of the uncertainty as understood by the data producer and greater freedom to consider different realisations of the data.

Description: Improved stratospheric aerosol extinction profiles from SCIAMACHY: validation and sample results Atmospheric Measurement Techniques Discussions, 8, 8353-8383, 2015 Author(s): C. von Savigny, F. Ernst, A. Rozanov, R. Hommel, K.-U. Eichmann, V. Rozanov, J. P. Burrows, and L. W. Thomason Stratospheric aerosol extinction profiles have been retrieved from SCIAMACHY/Envisat measurements of limb-scattered solar radiation. The retrieval is an improved version of an algorithm presented earlier. The retrieved aerosol extinction profiles are compared to co-located aerosol profile measurements with the SAGE II solar occultation instrument at a wavelength of 525 nm. Comparisons were carried out with two versions of the SAGE II data set (version 6.2 and the new version 7.0). In a global average sense the SCIAMACHY and the SAGE II version 7.0 extinction profiles agree to within about 10 % for altitudes above 15 km. Larger relative differences (up to 40 %) are observed at specific latitudes and altitudes. We also find differences between the two SAGE II data versions of up to 40 % for specific latitudes and altitudes. Sample results on the latitudinal and temporal variability of stratospheric aerosol extinction and optical depth during the SCIAMACHY mission period are presented. The results indicate that a series of volcanic eruptions is responsible for the increase in stratospheric aerosol optical depth from 2002 to 2012. Above about 28 km altitude volcanic eruptions are found to have negligible impact in the period 2002 to 2012.

Description: The tropopause inversion layer in baroclinic life cycles experiments: the role of diabatic and mixing processes Atmospheric Chemistry and Physics Discussions, 15, 21495-21537, 2015 Author(s): D. Kunkel, P. Hoor, and V. Wirth Recent studies on the formation of a quasi-permanent layer of enhanced static stability above the thermal tropopause revealed the contributions of dynamical and radiative processes. Dry dynamics lead to the evolution of a tropopause inversion layer (TIL) which is, however, too weak compared to observations and thus diabatic contributions are required. In this study we aim to assess the importance of diabatic as well as mixing processes in the understanding of TIL formation at midlatitudes. The non-hydrostatic model COSMO is applied in an idealized mid-latitude channel configuration to simulate baroclinic life cycles. The effect of individual diabatic, i.e. related to humidity and radiation, and turbulent processes is studied first to estimate the additional contribution of these processes to dry dynamics. In a second step these processes are stepwise included in the model to increase the complexity and finally estimate the relative importance of each process. The results suggest that including turbulence leads to a weaker TIL than in a dry reference simulation. In contrast, the TIL evolves stronger when radiation is included but the temporal occurrence is still comparable to the reference. Using various cloud schemes in the model shows that latent heat release and consecutive increased vertical motions foster an earlier and stronger appearance of the TIL than in all other life cycles. Furthermore, updrafts moisten the upper troposphere and as such increase the radiative effect from water vapor. Particularly, this process becomes more relevant for maintaining the TIL during later stages of the life cycles. Increased convergence of the vertical wind induced by updrafts and by propagating and potentially dissipating inertia-gravity waves further contributes to the enhanced stability of the lower stratosphere. Furthermore, radiative feedback of ice clouds reaching up to the tropopause is identified to potentially further affect the strength of the TIL in the region of the cloud.

Description: In this letter, a novel algorithm for attitude measurement based on a 3-D electromagnetic model (3-D em-model) is proposed. The 3-D em-model is established offline based on the geometric structure of the target, and it can be used to predict the scattering features at different target attitudes. In order to measure the attitude of the air target, we design a bistatic step frequency radar system. The directions of the two radars' lines of sight (LOSs) relative to the target are acquired by matching the high-resolution range profiles (HRRPs) from the target echoes to the HRRPs generated from the 3-D em-model. Since the directions of two radars' LOSs relative to the Earth are already known, the absolute attitude of the target can be acquired. The innovative contributions of this letter are as follows: 1) A comprehensive theoretical analysis of air target attitude measurement based on its own 3-D em-model is proposed; 2) the method can be applied to different kinds of air targets such as aircraft, satellite, missile, etc.; 3) the proposed attitude measurement method does not require target motion model in advance; and 4) the proposed algorithm can be applied to any kind of step frequency waveforms. Experiments using both data predicted by a high-frequency electromagnetic code and data measured in the chamber verify the validity of the method.

Description: The Soil Moisture and Ocean Salinity (SMOS) and Soil Moisture Active Passive (SMAP) missions provide Level-1 brightness temperature (Tb) observations that are used for global soil moisture estimation. However, the nature of these Tb data differs: the SMOS Tb observations contain atmospheric and select reflected extraterrestrial (“Sky”) radiation, whereas the SMAP Tb data are corrected for these contributions, using auxiliary near-surface information. Furthermore, the SMOS Tb observations are multiangular, whereas the SMAP Tb is measured at 40° incidence angle only. This letter discusses how SMOS Tb, SMAP Tb, and radiative transfer modeling components can be aligned in order to enable a seamless exchange of SMOS and SMAP Tb data in soil moisture retrieval and assimilation systems. The aggregated contribution of the atmospheric and reflected Sky radiation is, on average, about 1 K for horizontally polarized Tb and 0.5 K for vertically polarized Tb at 40° incidence angle, but local and short-term values regularly exceed 5 K.

Description: A multi-sensor addressing method for fiber Bragg grating aided fiber loop ringdown (FBG-FLRD) sensor array is proposed. It is capable of simultaneous measurement of temperature and force. Light from a wavelength-swept fiber laser (WSFL) is modulated into pulsed light to illuminate FBG-FLRD sensor array. Based on the time sequence of pulsed lights reflected by fiber Bragg grating array, each FBG-FLRD sensor can be distinguished. The time interval between the trigger signal of the WSFL and the pulsed light reflected by FBG offers temperature variation information. By measuring the ringdown time of each FBG-FLRD, the loss of fiber ring induced by temperature and force can be obtained. And then simultaneous measurement can be implemented in real time. To prove the validity of the proposed system, a six elements array is experimentally validated with interrogation frequency of 27.5 Hz.

Description: Impact of model developments on present and future simulations of permafrost in a global land-surface model The Cryosphere, 9, 1505-1521, 2015 Author(s): S. E. Chadburn, E. J. Burke, R. L. H. Essery, J. Boike, M. Langer, M. Heikenfeld, P. M. Cox, and P. Friedlingstein There is a large amount of organic carbon stored in permafrost in the northern high latitudes, which may become vulnerable to microbial decomposition under future climate warming. In order to estimate this potential carbon–climate feedback it is necessary to correctly simulate the physical dynamics of permafrost within global Earth system models (ESMs) and to determine the rate at which it will thaw. Additional new processes within JULES, the land-surface scheme of the UK ESM (UKESM), include a representation of organic soils, moss and bedrock and a modification to the snow scheme; the sensitivity of permafrost to these new developments is investigated in this study. The impact of a higher vertical soil resolution and deeper soil column is also considered. Evaluation against a large group of sites shows the annual cycle of soil temperatures is approximately 25 % too large in the standard JULES version, but this error is corrected by the model improvements, in particular by deeper soil, organic soils, moss and the modified snow scheme. A comparison with active layer monitoring sites shows that the active layer is on average just over 1 m too deep in the standard model version, and this bias is reduced by 70 cm in the improved version. Increasing the soil vertical resolution allows the full range of active layer depths to be simulated; by contrast, with a poorly resolved soil at least 50 % of the permafrost area has a maximum thaw depth at the centre of the bottom soil layer. Thus all the model modifications are seen to improve the permafrost simulations. Historical permafrost area corresponds fairly well to observations in all simulations, covering an area between 14 and 19 million km 2 . Simulations under two future climate scenarios show a reduced sensitivity of permafrost degradation to temperature, with the near-surface permafrost loss per degree of warming reduced from 1.5 million km 2 °C −1 in the standard version of JULES to between 1.1 and 1.2 million km 2 °C −1 in the new model version. However, the near-surface permafrost area is still projected to approximately half by the end of the 21st century under the RCP8.5 scenario.

Description: Numerical simulations of the Cordilleran ice sheet through the last glacial cycle The Cryosphere Discussions, 9, 4147-4203, 2015 Author(s): J. Seguinot, I. Rogozhina, A. P. Stroeven, M. Margold, and J. Kleman Despite more than a century of geological observations, the Cordilleran ice sheet of North America remains poorly understood in terms of its former extent, volume and dynamics. Although geomorphological evidence is abundant, its complexity is such that whole ice-sheet reconstructions of advance and retreat patterns are lacking. Here we use a numerical ice sheet model calibrated against field-based evidence to attempt a quantitative reconstruction of the Cordilleran ice sheet history through the last glacial cycle. A series of simulations is driven by time-dependent temperature offsets from six proxy records located around the globe. Although this approach reveals large variations in model response to evolving climate forcing, all simulations produce two major glaciations during marine oxygen isotope stages 4 (61.9–56.5 ka) and 2 (23.2–16.8 ka). The timing of glaciation is better reproduced using temperature reconstructions from Greenland and Antarctic ice cores than from regional oceanic sediment cores. During most of the last glacial cycle, the modelled ice cover is discontinuous and restricted to high mountain areas. However, widespread precipitation over the Skeena Mountains favours the persistence of a central ice dome throughout the glacial cycle. It acts as a nucleation centre before the Last Glacial Maximum and hosts the last remains of Cordilleran ice until the middle Holocene (6.6–6.2 ka).

Description: Stimulated Brillouin scattering during electron gyro-harmonic heating at EISCAT Annales Geophysicae, 33, 983-990, 2015 Author(s): H. Y. Fu, W. A. Scales, P. A. Bernhardt, S. J. Briczinski, M. J. Kosch, A. Senior, M. T. Rietveld, T. K. Yeoman, and J. M. Ruohoniemi Observations of secondary radiation, stimulated electromagnetic emission (SEE), produced during ionospheric modification experiments using ground-based, high-power, high-frequency (HF) radio waves are considered. The High Frequency Active Auroral Research Program (HAARP) facility is capable of generating narrowband SEE in the form of stimulated Brillouin scatter (SBS) and stimulated ion Bernstein scatter (SIBS) in the SEE spectrum. Such narrowband SEE spectral lines have not been reported using the European Incoherent Scatter (EISCAT) heater facility before. This work reports the first EISCAT results of narrowband SEE spectra and compares them to SEE previously observed at HAARP during electron gyro-harmonic heating. An analysis of experimental SEE data shows observations of emission lines within 100 Hz of the pump frequency, interpreted as SBS, during the 2012 July EISCAT campaign. Experimental results indicate that SBS strengthens as the pump frequency approaches the third electron gyro-harmonic. Also, for different heater antenna beam angles, the CUTLASS radar backscatter induced by HF radio pumping is suppressed near electron gyro-harmonics, whereas electron temperature enhancement weakens as measured by EISCAT/UHF radar. The main features of these new narrowband EISCAT observations are generally consistent with previous SBS measurements at HAARP.

Description: Optimizing the impact of temperature on bio-hydrogen production from food waste and its derivatives under no pH control using statistical modelling Biogeosciences Discussions, 12, 12823-12850, 2015 Author(s): A. Sattar, C. Arslan, C. Ji, S. Sattar, K. Yousaf, and S. Hashim The effect of temperature on bio-hydrogen production by co-digestion of sewerage sludge with food waste and its two derivatives, i.e. noodle waste and rice waste, was investigated by statistical modelling. Experimental results showed that increasing temperature from mesophilic (37 °C) to thermophilic (55 °C) was an effective mean for increasing bio-hydrogen production from food waste and noodle waste, but it caused a negative impact on bio-hydrogen production from rice waste. The maximum cumulative bio-hydrogen production of 650 mL was obtained from noodle waste under mesophilic temperature condition. Most of the production was observed during 48 h of incubation that continued till 72 h of incubation, and a decline in pH during this interval was 4.3 and 4.4 from a starting value of 7 under mesophilic and thermophilic conditions, respectively. Most of glucose consumption was also observed during 72 h of incubation and the maximum consumption was observed during the first 24 h, which was the same duration where the maximum pH drop occurred. The maximum hydrogen yields of 82.47 mL VS −1 , 131.38 mL COD −1 , and 44.90 mL glucose −1 were obtained from mesophilic food waste, thermophilic noodle waste and mesophilic rice waste respectively. The production of volatile fatty acids increased with an increase in time and temperature from food waste and noodle waste reactors whereas it decreased with temperature in rice waste reactors. The statistical modelling returned good results with high values of coefficient of determination ( R 2 ) for each waste type when it was opted for the study of cumulative hydrogen production, glucose consumption and volatile fatty acid production. The 3-D response surface plots developed by the statistical models helped a lot in developing better understanding of the impact of temperature and incubation time.

Description: Characterization and source apportionment of organic aerosol using offline aerosol mass spectrometry Atmospheric Measurement Techniques Discussions, 8, 8599-8644, 2015 Author(s): K. R. Daellenbach, C. Bozzetti, A. Křepelová, F. Canonaco, R. Wolf, P. Zotter, P. Fermo, M. Crippa, J. G. Slowik, Y. Sosedova, Y. Zhang, R.-J. Huang, L. Poulain, S. Szidat, U. Baltensperger, A. S. H. Prévôt, and I. El Haddad Field deployments of the Aerodyne Aerosol Mass Spectrometer (AMS) have significantly advanced real-time measurements and source apportionment of non-refractory particulate matter. However, the cost and complex maintenance requirements of the AMS make impractical its deployment at sufficient sites to determine regional characteristics. Furthermore, the negligible transmission efficiency of the AMS inlet for supermicron particles significantly limits the characterization of their chemical nature and contributing sources. In this study, we utilize the AMS to characterize the water-soluble organic fingerprint of ambient particles collected onto conventional quartz filters, which are routinely sampled at many air quality sites. The method was applied to 256 particulate matter (PM) filter samples (PM 1 , PM 2.5 , PM 10 ) collected at 16 urban and rural sites during summer and winter. We show that the results obtained by the present technique compare well with those from co-located online measurements, e.g. AMS or Aerosol Chemical Speciation Monitor (ACSM). The bulk recoveries of organic aerosol (60–91 %) achieved using this technique, together with low detection limits (0.8 μg of organic aerosol on the analyzed filter fraction) allow its application to environmental samples. We will discuss the recovery variability of individual hydrocarbon, oxygen containing and other ions. The performance of such data in source apportionment is assessed in comparison to ACSM data. Recoveries of organic components related to different sources as traffic, wood burning and secondary organic aerosol are presented. This technique, while subjected to the limitations inherent to filter-based measurements (e.g. filter artifacts and limited time resolution) may be used to enhance the AMS capabilities in measuring size-fractionated, spatially-resolved long-term datasets.

Description: Statistical detection and modeling of the over-dispersion of winter storm occurrence Natural Hazards and Earth System Science, 15, 1757-1761, 2015 Author(s): M. Raschke In this communication, I improve the detection and modeling of the over-dispersion of winter storm occurrence. For this purpose, the generalized Poisson distribution and the Bayesian information criterion are introduced; the latter is used for statistical model selection. Moreover, I replace the frequently used dispersion statistics by an over-dispersion parameter which does not depend on the considered return period of storm events. These models and methods are applied in order to properly detect the over-dispersion in winter storm data for Germany, carrying out a joint estimation of the distribution models for different samples.

Description: Impacts of beaver dams on hydrologic and temperature regimes in a mountain stream Hydrology and Earth System Sciences, 19, 3541-3556, 2015 Author(s): M. Majerova, B. T. Neilson, N. M. Schmadel, J. M. Wheaton, and C. J. Snow Beaver dams affect hydrologic processes, channel complexity, and stream temperature in part by inundating riparian areas, influencing groundwater–surface water interactions, and changing fluvial processes within stream systems. We explored the impacts of beaver dams on hydrologic and temperature regimes at different spatial and temporal scales within a mountain stream in northern Utah over a 3-year period spanning pre- and post-beaver colonization. Using continuous stream discharge, stream temperature, synoptic tracer experiments, and groundwater elevation measurements, we documented pre-beaver conditions in the first year of the study. In the second year, we captured the initial effects of three beaver dams, while the third year included the effects of ten dams. After beaver colonization, reach-scale (~ 750 m in length) discharge observations showed a shift from slightly losing to gaining. However, at the smaller sub-reach scale (ranging from 56 to 185 m in length), the discharge gains and losses increased in variability due to more complex flow pathways with beaver dams forcing overland flow, increasing surface and subsurface storage, and increasing groundwater elevations. At the reach scale, temperatures were found to increase by 0.38 °C (3.8 %), which in part is explained by a 230 % increase in mean reach residence time. At the smallest, beaver dam scale (including upstream ponded area, beaver dam structure, and immediate downstream section), there were notable increases in the thermal heterogeneity where warmer and cooler niches were created. Through the quantification of hydrologic and thermal changes at different spatial and temporal scales, we document increased variability during post-beaver colonization and highlight the need to understand the impacts of beaver dams on stream ecosystems and their potential role in stream restoration.

Description: Does the Budyko curve reflect a maximum power state of hydrological systems? A backward analysis Hydrology and Earth System Sciences Discussions, 12, 7821-7842, 2015 Author(s): M. Westhoff, E. Zehe, P. Archambeau, and B. Dewals Almost all catchments plot within a small envelope around the Budyko curve. This apparent behaviour suggests that organizing principles may play a role in the evolution of catchments. In this paper we applied the thermodynamic principle of maximum power as the organizing principle. In a top-down approach we derived mathematical formulations of the relation between relative wetness and gradients driving runoff and evaporation for a simple one-box model. We did this in such a way that when the conductances are optimized with the maximum power principle, the steady state behaviour of the model leads exactly to a point on the Budyko curve. Subsequently we derived gradients that, under constant forcing, resulted in a Budyko curve following the asymptotes closely. With these gradients we explored the sensitivity of dry spells and dynamics in actual evaporation. Despite the simplicity of the model, catchment observations compare reasonably well with the Budyko curves derived with dynamics in rainfall and evaporation. This indicates that the maximum power principle may be used (i) to derive the Budyko curve and (ii) to move away from the empiricism in free parameters present in many Budyko functions. Future work should focus on better representing the boundary conditions of real catchments and eventually adding more complexity to the model.

Description: Impacts of land use change and climate variations on annual inflow into Miyun Reservoir, Beijing, China Hydrology and Earth System Sciences Discussions, 12, 7785-7819, 2015 Author(s): J. K. Zheng, G. Sun, W. H. Li, X. X. Yu, C. Zhang, Y. B. Gong, and L. H. Tu Miyun reservoir, the only surface water source for Beijing city, has experienced water supply decline in recent decades. Previous studies suggest that both land use change and climate contributes the changes of water supply in this critical watershed. However, the specific causes of the decline in Miyun reservoir are debatable in a non-stationary climate in the past four decades. The central objective of this study was to quantify the separate and collective contributions of land use change and climate variability to the decreasing inflow into Miyun reservoir during 1961–2008. Different from previous studies, this work objectively identified breakpoints by analyzing the long-term historical hydrometeorology and land cover records. To effectively study the different impacts of the climate variation and land cover change during different sub-periods, annual water balance model (AWB), climate elasticity model (CEM), and rainfall–runoff model (RRM) were employed to conduct attribution analysis synthetically. We found a significant decrease in annual streamflow ( p 〈 0.01), a significant positive trend in annual potential evapotranspiration ( p 〈 0.01), and an insignificant negative trend in annual precipitation ( p 〉 0.1) during 1961–2008. Combined with historical records, we identified two breakpoints as in 1983 and 1999 for the period 1961–2008 by the sequential Mann–Kendall Test and Double Mass Curve. Climate variability alone did not explain the decrease in inflow to Miyun reservoir. Reduction of water yield was closely related to increase in evapotranspiration rates due to the expansion of forestlands and reduction in cropland and grassland, and was likely exacerbated by increased water consumption for domestic and industrial uses in the basin. Our study found that the contribution to the observed streamflow decline from land use change fell from 64–92 % during 1984–1999 to 36–58 % during 2000–2008, whereas the contribution from climate variation climbed from 8–36 % during the 1984–1999 to 42–64 % during 2000–2008. Model uncertainty analysis further demonstrated that climate warming played a dominant role in streamflow reduction in the 2000s. We conclude that future climate change and variability will further challenge the goal of water supply of Miyun reservoir to meet water demand. A comprehensive watershed management strategy needs to consider the climate variations besides vegetation management.

Description: Strange VLF bursts in northern Scandinavia: case study of the afternoon "mushroom-like" hiss on 8 December 2013 Annales Geophysicae, 33, 991-995, 2015 Author(s): J. Manninen, N. G. Kleimenova, A. Kozlovsky, I. A. Kornilov, L. I. Gromova, Y. V. Fedorenko, and T. Turunen We investigate a non-typical very low frequency (VLF) 1–4 kHz hiss representing a sequence of separated noise bursts with a strange "mushroom-like" shape in the frequency–time domain, each one lasting several minutes. These strange afternoon VLF emissions were recorded at Kannuslehto (KAN, ϕ = 67.74° N, λ = 26.27° E; L ∼ 5.5) in northern Finland during the late recovery phase of the small magnetic storm on 8 December 2013. The left-hand (LH) polarized 2–3 kHz "mushroom caps" were clearly separated from the right-hand (RH) polarized "mushroom stems" at the frequency of about 1.8–1.9 kHz, which could match the lower ionosphere waveguide cutoff (the first transverse resonance of the Earth–ionosphere cavity). We hypothesize that this VLF burst sequence could be a result of the modulation of the VLF hiss electron–cyclotron instability from the strong Pc5 geomagnetic pulsations observed simultaneously at ground-based stations as well as in the inner magnetosphere by the Time History of Events and Macroscale Interactions during Substorms mission probe (THEMIS-E; ThE). This assumption is confirmed by a similar modulation of the intensity of the energetic (1–10 keV) electrons simultaneously observed by the same ThE spacecraft. In addition, the data of the European Incoherent Scatter Scientific Association (EISCAT) radar at Tromsø show a similar quasi-periodicity in the ratio of the Hall-to-Pedersen conductance, which may be used as a proxy for the energetic particle precipitation enhancement. Our findings suggest that this strange mushroom-like shape of the considered VLF hiss could be a combined mutual effect of the magnetospheric ULF–VLF (ultra low frequency–very low frequency) wave interaction and the ionosphere waveguide propagation.

Description: Radiative and climate impacts of a large volcanic eruption during stratospheric sulfur geoengineering Atmospheric Chemistry and Physics Discussions, 15, 21837-21881, 2015 Author(s): A. Laakso, H. Kokkola, A.-I. Partanen, U. Niemeier, C. Timmreck, K. E. J. Lehtinen, H. Hakkarainen, and H. Korhonen Both explosive volcanic eruptions, which emit sulfur dioxide into the stratosphere, and stratospheric geoengineering via sulfur injections can potentially cool the climate by increasing the amount of scattering particles in the atmosphere. Here we employ a global aerosol-climate model and an earth system model to study the radiative and climate impacts of an erupting volcano during solar radiation management (SRM). According to our simulations, the radiative impacts of an eruption and SRM are not additive: in the simulated case of concurrent eruption and SRM, the peak increase in global forcing is about 40 % lower compared to a corresponding eruption into a clean background atmosphere. In addition, the recovery of the stratospheric sulfate burden and forcing was significantly faster in the concurrent case since the sulfate particles grew larger and thus sedimented faster from the stratosphere. In our simulation where we assumed that SRM would be stopped immediately after a volcano eruption, stopping SRM decreased the overall stratospheric aerosol load. For the same reasons, a volcanic eruption during SRM lead to only about 1/3 of the peak global ensemble-mean cooling compared to an eruption under unperturbed atmospheric conditions. Furthermore, the global cooling signal was seen only for 12 months after the eruption in the former scenario compared to over 40 months in the latter. In terms of the global precipitation rate, we obtain a 36 % smaller decrease in the first year after the eruption and again a clearly faster recovery in the concurrent eruption and SRM scenario. We also found that an explosive eruption could lead to significantly different regional climate responses depending on whether it takes place during geoengineering or into an unperturbed background atmosphere. Our results imply that observations from previous large eruptions, such as Mt Pinatubo in 1991, are not directly applicable when estimating the potential consequences of a volcanic eruption during stratospheric geoengineering.