Controls on fluvial carbon efflux from eroding peatland catchments

Global peatlands store an unparalleled proportion of total global organic carbon but it is vulnerable to erosion into fluvial systems. Fluvial networks are being recognized as areas of carbon transformation, with eroded particulate organic carbon processed to dissolved organic carbon and CO₂. Existing studies indicate biodegradation and photodegradation as key processes controlling the transformation of organic carbon in fluvial systems, with initial concentrations of dissolved organic carbon (DOC) identified as a control on the rate of carbon mineralization. This study manipulates temperature and incident light intensity to investigate carbon mineralization rates in laboratory simulations of peatland sediment transport into fluvial systems. By directly measuring gaseous CO₂ emissions from sampled stream water, the relationship of temperature and light intensity with carbon efflux is identified. In simulations where sediment (as particulate organic matter, POM) is absent, temperature is consistently the dominant factor influencing carbon efflux rates. This influence is independent of the initial DOC concentration of the water sample. In simulations where POM was added, representing a peatland river receiving eroded terrestrial sediment, initial DOC concentration predicts 79% of the variation in total gaseous carbon efflux whereas temperature and light intensity predict 12% and 3%, respectively. When sampled stream water's mineralization rates in the presence of added POM are analysed independently, removing DOC as a model variable, the dominant variable affecting CO₂ efflux is opposite for each sample. This study presents novel data suggesting peatland erosion introduces further complexity to dynamic stream systems where rates of carbon transformation processes and the influence of specific environmental variables are interdependent. Anthropogenic climate change is identified as a leading risk factor perpetuating peatland erosion; therefore, understanding the fate of terrestrial sediment in rivers and further quantifying the benefits of protecting peatland soils will be of increasing importance to carbon budgeting and ecosystem function studies.     

Assessing Data Completeness of VGI through an Automated Matching Procedure for Linear Data

New, free and fast growing spatial data sources have appeared online, based on Volunteered Geographic Information (VGI). OpenStreetMap (OSM) is one of the most representative projects of this trend. Its increasing popularity and density makes the study of its data quality an imperative. A common approach is to compare OSM with a reference dataset. In such cases, data matching is necessary for the comparison to be meaningful, and is usually performed manually at the data preparation stage. This article proposes an automated feature‐based matching method specifically designed for VGI, based on a multi‐stage approach that combines geometric and attribute constraints. It is applied to the OSM dataset using the official data from Ordnance Survey as the reference dataset. The results are then used to evaluate data completeness of OSM in several case studies in the UK.     

Smoothing/filtering LiDAR digital surface models. Experiments with loess regression and discrete wavelets

This paper reports on the smoothing/filtering analysis of a digital surface model (DSM) derived from LiDAR altimetry for part of the River Coquet, Northumberland, UK using loess regression and the 2D discrete wavelet transform (DWT) implemented in the S-PLUS and R statistical packages. The chosen method of analysis employs a simple method to generate ‘noise’ which is then added to a smooth sample of LiDAR data; loess regression and wavelet methods are then used to smooth/filter this data and compare with the original ‘smooth’ sample in terms of RMSE. Various combinations of functions and parameters were chosen for both methods. Although wavelet analysis was effective in filtering the noise from the data, loess regression employing a quadratic parametric function produced the lowest RMSE and was the most effective.     

Understanding the Geometry of Connected Fracture Flow with Multiperiod Oscillatory Hydraulic Tests

An understanding of the spatial and hydraulic properties of fast preferential flow pathways in the subsurface is necessary in applications ranging from contaminant fate and transport modeling to design of energy extraction systems. One method for the characterization of fracture properties over interwellbore scales is Multiperiod Oscillatory Hydraulic (MOH) testing, in which the aquifer response to oscillatory pressure stimulations is observed. MOH tests were conducted on isolated intervals of wells in siliciclastic and carbonate aquifers in southern Wisconsin. The goal was to characterize the spatial properties of discrete fractures over interwellbore scales. MOH tests were conducted on two discrete fractured intervals intersecting two boreholes at one field site, and a nest of three piezometers at another field site. Fracture diffusivity estimates were obtained using analytical solutions that relate diffusivity to observed phase lag and amplitude decay. In addition, MOH tests were used to investigate the spatial extent of flow using different conceptual models of fracture geometry. Results indicated that fracture geometry at both field sites can be approximated by permeable two‐dimensional fracture planes, oriented near‐horizontally at one site, and near‐vertically at the other. The technique used on MOH field data to characterize fracture geometry shows promise in revealing fracture network characteristics important to groundwater flow and transport.     

Using high‐resolution isotope data and alternative calibration strategies for a tracer‐aided runoff model in a nested catchment

Testing hydrological models over different spatio‐temporal scales is important for both evaluating diagnostics and aiding process understanding. High‐frequency (6‐hr) stable isotope sampling of rainfall and runoff was undertaken during 3‐week periods in summer and winter within 12 months of daily sampling in a 3.2‐km² catchment in the Scottish Highlands. This was used to calibrate and test a tracer‐aided model to assess the (a) information content of high‐resolution data, (b) effect of different calibration strategies on simulations and inferred processes, and (c) model transferability to <1‐km² subcatchment. The 6‐hourly data were successfully incorporated without loss of model performance, improving the temporal resolution of the modelling, and making it more relevant to the time dynamics of the isotope and hydrometric response. However, this added little new information due to old‐water dominance and riparian mixing in this peatland catchment. Time variant results, from differential split sample testing, highlighted the importance of calibrating to a wide range of hydrological conditions. This also provided insights into the nonstationarity of catchment mixing processes, in relation to storage and water ages, which varied markedly depending on the calibration period. Application to the nested subcatchment produced equivalent parameterization and performance, highlighting similarity in dominant processes. The study highlighted the utility of high‐resolution data in combination with tracer‐aided models, applied at multiple spatial scales, as learning tools to enhance process understanding and evaluation of model behaviour across nonstationary conditions. This helps reveal more fully the catchment response in terms of the different mechanistic controls on both wave celerites and particle velocities.     

Low carbon dioxide partial pressure in a productive subtropical lake

The purposes of this study were to assess if Lake Apopka (FL, USA) was autotrophic or heterotrophic based on the partial pressure of dissolved carbon dioxide (pCO₂) in the surface water and to evaluate factors that influence the long-term changes in pCO₂. Monthly average pH, alkalinity and other limnological variables collected between 1987 and 2006 were used to estimate dissolved inorganic carbon (DIC), pCO₂ and CO₂ flux between surface water and atmosphere. Results indicated that average pCO₂ in the surface water was 196 μatm, well below the atmospheric pCO₂. Direct measurements of DIC concentration on three sampling dates in 2009 also supported pCO₂ undersaturation in Lake Apopka. Supersaturation in CO₂ occurred in this lake in only 13% of the samples from the 20-year record. The surface-water pCO₂ was inversely related to Chl a concentrations. Average annual CO₂ flux was 28.2 g C m⁻² year⁻¹ from the atmosphere to the lake water and correlated significantly with Chl a concentration, indicating that biological carbon sequestration led to the low dissolved CO₂ concentration. Low pCO₂ and high invasion rates of atmospheric CO₂ in Lake Apopka indicated persistent autotrophy. High rates of nutrient loading and primary production, a high buffering capacity, a lack of allochthonous loading of organic matter, and the dominance of a planktivorous–benthivorous fish food web have supported long-term net autotrophy in this shallow subtropical eutrophic lake. Our results also showed that lake restoration by the means of nutrient reduction resulted in significantly lower total phosphorus (TP) and Chl a concentrations, and higher pCO₂.     

Assessing urban rainfall-runoff response to stormwater management extent

We report an empirical analysis of the hydrologic response of three small, highly impervious urban watersheds to pulse rainfall events, to assess how traditional stormwater management (SWM) alters urban hydrographs. The watersheds vary in SWM coverage from 3% to 61% and in impervious cover from 45% to 67%. By selecting a set of storm events that involved a single rainfall pulse with >96% of total precipitation delivered in 60 min, we reduced the effect of differences between storms on hydrograph response to isolate characteristic responses attributable to watershed properties. Watershed-average radar rainfall data were used to generate local storm hyetographs for each event in each watershed, thus compensating for the extreme spatial and temporal heterogeneity of short-duration, intense rainfall events. By normalizing discharge values to the discharge peak and centring each hydrograph on the time of peak we were able to visualize the envelope of hydrographs for each group and to generate representative composite hydrographs for comparison across the three watersheds. Despite dramatic differences in the fraction of watershed area draining to SWM features across these three headwater tributaries, we did not find strong evidence that SWM causes significant attenuation of the hydrograph peak. Hydrograph response for the three watersheds is remarkably uniform despite contrasts in SWM, impervious cover and spatial patterns of land cover type. The primary difference in hydrograph response is observed on the recession limb of the hydrograph, and that change appears to be associated with higher storm-total runoff in the watersheds with more area draining to SWM. Our findings contribute more evidence to the work of previous authors suggesting that SWM is less effective at attenuating urban hydrographs than is commonly assumed. Our findings also are consistent with previous work concluding that percent impervious cover may have greater influence on runoff volume than percent SWM coverage.     

Evidence Of Dynamical Coupling Between The Residual Layer And The Developing Convective Boundary Layer

The diurnal cycle of the atmospheric boundary layer (ABL) hasbeen documented on 8 August 1998 in the framework of the Étude et Simulation de la QUalité de l'air en Ile-de-France (ESQUIF) experiment that took place in the Paris area. The ABL structure was documented by means of a ground-based lidar, surface meteorological stations and soundings. The interaction between the residual layer and the convective boundary layer is investigated using the collected data as well as mesoscale modelling. As opposed to the generally accepted concept, we find evidence of entrainment at the top ofthe residual layer. High temporal simulations of the 8 August 1998 casemade with the mesoscale atmospheric model Meso-NH also show evidenceof mixing at the top of the residual layer (RL). This mixing is believed to be related to the presence of convective (gravity) waves in the RL.     

Lightning evolution related to radar-derived microphysics in the 21 July 1998 EULINOX supercell storm

Results of a combined analysis of data from a C-band polarimetric Doppler radar and a 3D VHF interferometric lightning mapping system, as obtained during the European Lightning Nitrogen Oxides project (EULINOX) field campaign, are presented. For 21 July 1998, the lightning data from a supercell thunderstorm weakly indicate a tendency for a bi-level vertical distribution of lightning VHF emissions around the −15°C and −30°C temperature levels. Also, in some parts of the clouds, evidence is found for the presence of a lower positive charge center near the freezing level. However, where strong vertical motions prevail, VHF emissions are not organized in horizontal layers but in oblique or vertical regions. Correlation of VHF signals with radar quantities shows that in the growing storm, peak VHF activity is low and related to reflectivity factors around 30 dBZ, while after the mature stage, the peak VHF activity is about three times larger. The highest density of VHF signals is now found near reflectivity factors of 45 dBZ. A polarimetric hydrometeor classification indicates that during storm development, most lightning activity occurs where graupel and, secondarily, snow and small dry hail are present. In the decaying phase of the supercell hailstorm, however, most lightning VHF emissions stem from the region with hail and heavy rain. Furthermore, while the VHF signal frequency per cubic kilometer in the graupel and rain regions remains nearly constant throughout the supercell life cycle, the signal frequency in the hail region rises during storm decay.     

Idealized Mesoscale Model Simulations of Open Cellular Convection Over the Sea

The atmospheric conditions during an observed case of open cellular convection over the North Sea were simulated using the Weather Research and Forecasting (WRF) numerical model. Wind, temperature and water vapour mixing ratio profiles from the WRF simulation were used to initialize an idealized version of the model, which excluded the effects of topography, surface inhomogeneities and large-scale weather forcing. Cells with an average diameter of 17.4 km developed. Simulations both with and without a capping inversion were made, and the cell-scale kinetic energy budget was calculated for each case. By considering all sources of explicit diffusion in the model, the budgets were balanced. In comparison with previous work based on observational studies, the use of three-dimensional, gridded model data afforded the possibility of calculating all terms in the budgets, which showed that the important terms in the budgets were buoyancy, pressure balance and inter-scale transfer to subgrid scales. Cells were also composited to calculate the average cell-scale flow and each of the budget terms on two-dimensional cross-sections through the cells, parallel and perpendicular to the mean wind direction.