ALBERT

Description: This paper describes the Global Climate Observing System (GCOS) Reference Upper-Air Network (GRUAN) approach to managing the transition from the Vaisala RS92 to the Vaisala RS41 as the operational radiosonde. The goal of GRUAN is to provide long-term high-quality reference observations of upper-air essential climate variables (ECVs) such as temperature and water vapor. With GRUAN data being used for climate monitoring, it is vital that the change of measurement system does not introduce inhomogeneities to the data record. The majority of the 27 GRUAN sites were launching the RS92 as their operational radiosonde, and following the end of production of the RS92 in the last quarter of 2017, most of these sites have now switched to the RS41. Such a large-scale change in instrumentation is unprecedented in the history of GRUAN and poses a challenge for the network. Several measurement programs have been initiated to characterize differences in biases, uncertainties, and noise between the two radiosonde types. These include laboratory characterization of measurement errors, extensive twin sounding studies with RS92 and RS41 on the same balloon, and comparison with ancillary data. This integrated approach is commensurate with the GRUAN principles of traceability and deliberate redundancy. A 2-year period of regular twin soundings is recommended, and for sites that are not able to implement this, burden-sharing is employed such that measurements at a certain site are considered representative of other sites with similar climatological characteristics. All data relevant to the RS92–RS41 transition are archived in a database that will be accessible to the scientific community for external scrutiny. Furthermore, the knowledge and experience gained regarding GRUAN's RS92–RS41 transition will be extensively documented to ensure traceability of the process. This documentation will benefit other networks in managing changes in their operational radiosonde systems. Preliminary analysis of the laboratory experiments indicates that the manufacturer's calibration of the RS41 temperature and humidity sensors is more accurate than for the RS92, with uncertainties of

Description: Trend estimates with different signs are reported in the literature for lower stratospheric water vapour considering the time period between the late 1980s and 2010. The NOAA (National Oceanic and Atmospheric Administration) frost point hygrometer (FPH) observations at Boulder (Colorado, 40.0° N, 105.2° W) indicate positive trends (about 0.1 to 0.45 ppmv decade−1). On the contrary, negative trends (approximately −0.2 to −0.1 ppmv decade−1) are derived from a merged zonal mean satellite data set for a latitude band around the Boulder latitude. Overall, the trend differences between the two data sets range from about 0.3 to 0.5 ppmv decade−1, depending on altitude. It has been proposed that a possible explanation for these discrepancies is a different temporal behaviour at Boulder and the zonal mean. In this work we investigate trend differences between Boulder and the zonal mean using primarily simulations from ECHAM/MESSy (European Centre for Medium-Range Weather Forecasts Hamburg/Modular Earth Submodel System) Atmospheric Chemistry (EMAC), WACCM (Whole Atmosphere Community Climate Model), CMAM (Canadian Middle Atmosphere Model) and CLaMS (Chemical Lagrangian Model of the Stratosphere). On shorter timescales we address this aspect also based on satellite observations from UARS/HALOE (Upper Atmosphere Research Satellite/Halogen Occultation Experiment), Envisat/MIPAS (Environmental Satellite/Michelson Interferometer for Passive Atmospheric Sounding) and Aura/MLS (Microwave Limb Sounder). Overall, both the simulations and observations exhibit trend differences between Boulder and the zonal mean. The differences are dependent on altitude and the time period considered. The model simulations indicate only small trend differences between Boulder and the zonal mean for the time period between the late 1980s and 2010. These are clearly not sufficient to explain the discrepancies between the trend estimates derived from the FPH observations and the merged zonal mean satellite data set. Unless the simulations underrepresent variability or the trend differences originate from smaller spatial and temporal scales than resolved by the model simulations, trends at Boulder for this time period should also be quite representative for the zonal mean and even other latitude bands. Trend differences for a decade of data are larger and need to be kept in mind when comparing results for Boulder and the zonal mean on this timescale. Beyond that, we find that the trend estimates for the time period between the late 1980s and 2010 also significantly differ among the simulations. They are larger than those derived from the merged satellite data set and smaller than the trend estimates derived from the FPH observations.

Description: Within the framework of the second SPARC (Stratosphere-troposphere Processes And their Role in Climate) water vapour assessment (WAVAS-II), profile-to-profile comparisons of stratospheric and lower mesospheric water vapour were performed by considering 33 data sets derived from satellite observations of 15 different instruments. These comparisons aimed to provide a picture of the typical biases and drifts in the observational database and to identify data-set-specific problems. The observational database typically exhibits the largest biases below 70 hPa, both in absolute and relative terms. The smallest biases are often found between 50 and 5 hPa. Typically, they range from 0.25 to 0.5 ppmv (5 % to 10 %) in this altitude region, based on the 50 % percentile over the different comparison results. Higher up, the biases increase with altitude overall but this general behaviour is accompanied by considerable variations. Characteristic values vary between 0.3 and 1 ppmv (4 % to 20 %). Obvious data-set-specific bias issues are found for a number of data sets. In our work we performed a drift analysis for data sets overlapping for a period of at least 36 months. This assessment shows a wide range of drifts among the different data sets that are statistically significant at the 2σ uncertainty level. In general, the smallest drifts are found in the altitude range between about 30 and 10 hPa. Histograms considering results from all altitudes indicate the largest occurrence for drifts between 0.05 and 0.3 ppmv decade−1. Comparisons of our drift estimates to those derived from comparisons of zonal mean time series only exhibit statistically significant differences in slightly more than 3 % of the comparisons. Hence, drift estimates from profile-to-profile and zonal mean time series comparisons are largely interchangeable. As for the biases, a number of data sets exhibit prominent drift issues. In our analyses we found that the large number of MIPAS data sets included in the assessment affects our general results as well as the bias summaries we provide for the individual data sets. This is because these data sets exhibit a relative similarity with respect to the remaining data sets, despite the fact that they are based on different measurement modes and different processors implementing different retrieval choices. Because of that, we have by default considered an aggregation of the comparison results obtained from MIPAS data sets. Results without this aggregation are provided on multiple occasions to characterise the effects due to the numerous MIPAS data sets. Among other effects, they cause a reduction of the typical biases in the observational database.

Description: Retrievals of vertical profiles of key atmospheric gases provide a critical long-term record from ground-based Fourier transform infrared (FTIR) solar absorption measurements. However, the characterization of the retrieved vertical profile structure can be difficult to validate, especially for gases with large vertical gradients and spatial–temporal variability such as water vapor. In this work, we evaluate the accuracy of the most common water vapor isotope (H216O, hereafter WV) FTIR retrievals in the lower and upper troposphere–lower stratosphere. Coincident high-quality vertically resolved WV profile measurements obtained from 2010 to 2016 with balloon-borne NOAA frost point hygrometers (FPHs) are used as reference to evaluate the performance of the retrieved profiles at two sites: Boulder (BLD), Colorado, and at the mountaintop observatory of Mauna Loa (MLO), Hawaii. For a meaningful comparison, the spatial–temporal variability has been investigated. We present results of comparisons among FTIR retrievals with unsmoothed and smoothed FPH profiles to assess WV vertical gradients. Additionally, we evaluate the quantitative impact of different a priori profiles in the retrieval of WV. An orthogonal linear regression analysis shows the best correlation among tropospheric layers using ERA-Interim (ERA-I) a priori profiles and biases are lower for unsmoothed comparisons. In Boulder, we found a negative bias of 0.02±1.9 % (r=0.95) for the 1.5–3 km layer. A larger negative bias of 11.1±3.5 % (r=0.97) was found in the lower free troposphere layer of 3–5 km attributed to rapid vertical change of WV, which is not always captured by the retrievals. The bias improves in the 5–7.5 km layer (1.0±5.3 %, r=0.94). The bias remains at about 13 % for layers above 7.5 km but below 13.5 km. At MLO the spatial mismatch is significantly larger due to the launch of the sonde being farther from the FTIR location. Nevertheless, we estimate a negative bias of 5.9±4.6 % (r=0.93) for the 3.5–5.5 km layer and 9.9±3.7 % (r=0.93) for the 5.5–7.5 km layer, and we measure positive biases of 6.2±3.6 % (r=0.95) for the 7.5–10 km layer and 12.6 % and greater values above 10 km. The agreement for the first layer is significantly better at BLD because the air masses are similar for both FTIR and FPH. Furthermore, for the first time we study the influence of different WV a priori profiles in the retrieval of selected gas profiles. Using NDACC standard retrievals we present results for hydrogen cyanide (HCN), carbon monoxide (CO), and ethane (C2H6) by taking NOAA FPH profiles as the ground truth and evaluating the impact of other WV profiles. We show that the effect is minor for C2H6 (bias

Description: Because fossil fuel derived CO2 is the only source of atmospheric CO2 that is devoid of 14C, atmospheric measurements of Δ14CO2 can be used to constrain fossil fuel emissions at local and regional scales. However, at the continental scale, atmospheric transport and other sources of variability in Δ14CO2 may influence the fossil fuel detection capability. We present a set of Δ14CO2 observations from the train-based TROICA-8 expedition across Eurasia in March–April 2004. Local perturbations in Δ14CO2 are caused by easily identifiable sources from nuclear reactors and localized pollution events. The remaining data show an increase in Δ14CO2 from Western Russia (40° E) to Eastern Siberia (120° E), consistent with depletion in 14CO2 caused by fossil fuel CO2 emissions in heavily populated Europe, and gradual dispersion of the fossil fuel plume across Northern Asia. Other tracer gas species which may be correlated with fossil fuel CO2 emissions, including carbon monoxide, sulphur hexafluoride, and perchloroethylene, were also measured and the results compared with the Δ14CO2 measurements. The sulphur hexafluoride longitudinal gradient is not significant relative to the measurement uncertainty. Carbon monoxide and perchloroethylene show large-scale trends of enriched values in Western Russia and decreasing values in Eastern Siberia, consistent with fossil fuel emissions, but exhibit significant spatial variability, especially near their primary sources in Western Russia. The clean air Δ14CO2 observations are compared with simulated spatial gradients from the TM5 atmospheric transport model. We show that the change in Δ14CO2 across the TROICA transect is due almost entirely to emissions of fossil fuel CO2, but that the magnitude of this Δ14CO2 gradient is relatively insensitive to modest uncertainties in the fossil fuel flux. In contrast, the Δ14CO2 gradient is strongly sensitive to the modeled representation of vertical mixing, suggesting that Δ14CO2 may be a useful tracer for training mixing in atmospheric transport models.

Description: Because fossil fuel derived CO2 is the only source of atmospheric CO2 that is devoid of 14C, atmospheric measurements of Δ14CO2 can be used to constrain fossil fuel emission estimates at local and regional scales. However, at the continental scale, uncertainties in atmospheric transport and other sources of variability in Δ14CO2 may influence the fossil fuel detection capability. We present a set of Δ14CO2 observations from the train-based TROICA-8 expedition across Eurasia in March–April 2004. Local perturbations in Δ14CO2 are caused by easily identifiable sources from nuclear reactors and localized pollution events. The remaining data show an increase in Δ14CO2 from Western Russia (40° E) to Eastern Siberia (120° E), consistent with depletion in 14CO2 caused by fossil fuel CO2 emissions in heavily populated Europe, and gradual dispersion of the fossil fuel plume across Northern Asia. Other trace gas species which may be correlated with fossil fuel CO2 emissions, including carbon monoxide, sulphur hexafluoride, and perchloroethylene, were also measured and the results compared with the Δ14CO2 measurements. The sulphur hexafluoride longitudinal gradient is not significant relative to the measurement uncertainty. Carbon monoxide and perchloroethylene show large-scale trends of enriched values in Western Russia and decreasing values in Eastern Siberia, consistent with fossil fuel emissions, but exhibit significant spatial variability, especially near their primary sources in Western Russia. The clean air Δ14CO2 observations are compared with simulated spatial gradients from the TM5 atmospheric transport model. We show that the change in Δ14CO2 across the TROICA transect is due almost entirely to emissions of fossil fuel CO2, but that the magnitude of this Δ14CO2 gradient is relatively insensitive to modest uncertainties in the fossil fuel flux. In contrast, the Δ14CO2 gradient is more sensitive to the modeled representation of vertical mixing, suggesting that Δ14CO2 may be a useful tracer for training mixing in atmospheric transport models.

Description: Coastal defences have long been employed to halt or slow coastal erosion. Their impact on local sediment flux and ecology has been studied in detail through field studies and numerical simulations. The non-local impact of a modified sediment flux regime on mesoscale erosion and accretion has received less attention. Morphological changes at this scale due to defended structures can be difficult to quantify or identify with field data. Engineering scale numerical models, often applied to assess the design of modern defences on local coastal erosion, tend not to cover large stretches of coast and are rarely applied to assess the impact of older structures. We extend previous work to explore the influences of sea walls on the evolution and morphological sensitivity of a pinned, soft-cliff, sandy coastline under a changing wave climate. The Holderness coast of East Yorkshire, UK, is used as a case study, represented both as a defended example with major sea walls included and a natural example where no sea defences exist. Using a mesoscale numerical coastal evolution model, stochastic wave climate data are perturbed gradually to assess the sensitivity of the coastal morphology to changing wave climate for both the defended and natural scenarios. Comparative analysis of the simulated output suggests that sea walls in the south of the region have a greater impact on sediment flux due to the increased sediment availability along this part of the coast. Multiple defended structures, including those separated by several kilometres, were found to interact with each other, producing a complex imprint on coastal morphology under a changing wave climate. Although spatially and temporally heterogeneous, sea walls generally slowed coastal recession and accumulated sediment on their up-drift side.

Description: Linking landslide size and frequency is important at both human and geological time-scales for quantifying both landslide hazards and the effectiveness of landslides in the removal of sediment from evolving landscapes. Landslide inventories are usually compiled following a particular triggering event such as an earthquake or storm, and their statistical behavior is typically characterized by an inflected power-law relationship. The occurrence of landslides is expected to be influenced by the material properties of rock and/or regolith in which failure occurs. Here we explore the statistical behavior and the controls of a secular landslide inventory (SLI) (i.e. events occurring over an indefinite time period) consisting of mapped landslide deposits and their underlying lithology (bedrock or superficial) across the United Kingdom. The magnitude-frequency distribution of this secular inventory exhibits an inflected power law relationship, well approximated by an inverse Gamma or double Pareto model. The scaling exponent for the power-law relationship is α = −1.76. The small-event rollover occurs at a significantly higher magnitude than observed in single-event landslide records, which we interpret as evidence of “landscape annealing” at these relatively short length-scales, noting the corollary that a secular dataset will tend to underestimate the frequency of small landslides. This is supported by a subset of data where a complete landslide inventory was recently mapped. Large landslides also appear to be under-represented relative to model predictions, which we interpret as a non-linear or transient landscape response as the UK emerged from the last glacial maximum and through relatively volatile conditions toward a generally more stable late Holocene climate.

Description: The characteristics of the sediment transported by rivers (e.g., sediment flux, grain size distribution – GSD –) dictate whether rivers aggrade or erode their substrate. They also condition the architecture and properties of sedimentary successions in basins. In this study, we investigate the relationship between landscape steepness and the grain size of hillslope and fluvial sediments. The study area is located within the Feather River Basin in Northern California, and studied basins are underlain exclusively by tonalite lithology. Erosion rates in the study area vary over an order of magnitude, from 〉 250 mm ka−1 in the Feather River canyon to 〈 15 mm ka−1 on an adjacent low relief plateau. We find that the coarseness of hillslope sediment increases with increasing hillslope steepness and erosion rates. We hypothesize that, in our soil samples, the measured ten-fold increase in D50 and doubling of the amount of fragments larger than 1 mm when slope increases from 0.38 to 0.83 m m−1 is due to a decrease in the residence time of rock fragments, causing particles to be exposed for shorter periods of time to processes that can reduce grain size. For slopes in excess of 0.7 m m−1, landslides and scree cones supply much coarser sediment to rivers, with D50 and D84 more than one order of magnitude larger than in soils. In the tributary basins of the Feather River, a prominent break in slope developed in response to the rapid incision of the Feather River. Downstream of the break in slope, fluvial sediment grain size increases, due to an increase in flow competence (mostly driven by channel steepening) but also by a change in sediment source and in sediment dynamics: on the plateau upstream of the break in slope, rivers transport easily mobilised fine-grained sediment derived exclusively from soils. Downstream of the break in slope, mass wasting processes supply a wide range of grain sizes that rivers entrain selectively, depending on the competence of their flow. Our results also suggest that in this study site, hillslopes respond rapidly to an increase in the rate of base-level lowering compared to rivers.

Description: Quantifying rates of erosion on cliffed coasts across a range of timescales is vital for understanding the drivers and processes of coastal change and for assessing risks posed by future cliff retreat. Historical records cover at best the last 150 years; Cosmogenic radionuclides, such as 10Be could allow us to look further into past to assess coastal change at millenial timescales. CRNs accumulate in-situ near the Earth surface and have been used extensively to quantify erosion rates, burial dates and surface exposure ages in terrestrial landscapes over the last three decades. More recently, applications in rocky coast settings have quantified the timing of mass wasting events, determined long-term-averaged rates of cliff retreat and revealed the exposure history of shore platforms. In this contribution, we developed and explored a numerical model for the accumulation of 10Be on eroding shore platforms. In a series of numerical experiments, we investigated the influence of topographic and water shielding, dynamic platform erosion processes, the presence and variation in beach cover, and heterogeneous distribution of erosion on the distribution of 10Be across shore platforms. Results demonstrate that, taking into account relative sea level change and tides, the concentration of 10Be is sensitive to rates of cliff retreat. Factors such as topographic shielding and beach cover, act to reduce 10Be concentrations on the platform, and may result in overestimation of cliff retreat rates if not accounted for. The shape of the distribution of 10Be across a shore platform can potentially reveal whether cliff retreat rates are declining or accelerating through time. Measurement of 10Be in shore platforms has great potential to allow us to quantify long-term rates of cliff retreat and platform erosion.