A semi-distributed ntegrated itrogen model for multiple source assessment in tchments (INCA): Part I — model structure and process equations

A new model has been developed for assessing multiple sources of nitrogen in catchments. The model (INCA) is process based and uses reaction kinetic equations to simulate the principal mechanisms operating. The model allows for plant uptake, surface and sub-surface pathways and can simulate up to six land uses simultaneously. The model can be applied to catchment as a semi-distributed simulation and has an inbuilt multi-reach structure for river systems. Sources of nitrogen can be from atmospheric deposition, from the terrestrial environment (e.g. agriculture, leakage from forest systems etc.), from urban areas or from direct discharges via sewage or intensive farm units. The model is a daily simulation model and can provide information in the form of time series at key sites, or as profiles down river systems or as statistical distributions. The process model is described and in a companion paper the model is applied to the River Tywi catchment in South Wales and the Great Ouse in Bedfordshire.

Morphology induced spinodal decomposition at the surface of symmetric diblock copolymer films

Atomic force microscopy is used to study the ordering dynamics of symmetric diblock copolymer films. The films order to form a lamellar structure which results in a frustration when the film thickness is incommensurate with the lamellae. By probing the morphology of incommensurate films in the early ordering stages, we discover an intermediate phase of lamellae arranged perpendicular to the film surface. This morphology is accompanied by a continuous growth in amplitude of the film surface topography with a characteristic wavelength, indicative of a spinodal process. Using selfconsistent field theory, we show that the observation of perpendicular lamellae suggests an intermediate state with parallel lamellae at the substrate and perpendicular lamellae at the free surface. The calculations confirm that the intermediate state is unstable to thickness fluctuations, thereby driving the spinodal growth of surface structures.

A case study of sea breeze blocking regulated by sea surface temperature along the English south coast

The sensitivity of sea breeze structure to sea surface temperature (SST) and coastal orography is investigated in convection-permitting Met Office Unified Model simulations of a case study along the south coast of England. Changes in SST of 1 K are shown to significantly modify the structure of the sea breeze immediately offshore. On the day of the case study, the sea breeze was partially blocked by coastal orography, particularly within Lyme Bay. The extent to which the flow is blocked depends strongly on the static stability of the marine boundary layer. In experiments with colder SST, the marine boundary layer is more stable, and the degree of blocking is more pronounced. Although a colder SST would also imply a larger land–sea temperature contrast and hence a stronger onshore wind – an effect which alone would discourage blocking – the increased static stability exerts a dominant control over whether blocking takes place. The implications of prescribing fixed SST from climatology in numerical weather prediction model forecasts of the sea breeze are discussed.

Vertical structure of stratospheric water vapour trends derived from merged satellite data

Stratospheric water vapour is a powerful greenhouse gas. The longest available record from balloon observations over Boulder, Colorado, USA shows increases in stratospheric water vapour concentrations that cannot be fully explained by observed changes in the main drivers, tropical tropopause temperatures and methane. Satellite observations could help resolve the issue, but constructing a reliable long-term data record from individual short satellite records is challenging. Here we present an approach to merge satellite data sets with the help of a chemistry–climate model nudged to observed meteorology. We use the models’ water vapour as a transfer function between data sets that overcomes issues arising from instrument drift and short overlap periods. In the lower stratosphere, our water vapour record extends back to 1988 and water vapour concentrations largely follow tropical tropopause temperatures. Lower and mid-stratospheric long-term trends are negative, and the trends from Boulder are shown not to be globally representative. In the upper stratosphere, our record extends back to 1986 and shows positive long-term trends. The altitudinal differences in the trends are explained by methane oxidation together with a strengthened lower-stratospheric and a weakened upper stratospheric circulation inferred by this analysis. Our results call into question previous estimates of surface radiative forcing based on presumed global long-term increases in water vapour concentrations in the lower stratosphere.

Surface chemistry of alanine on Cu{111}: Adsorption geometry and temperature dependence

Adsorption of l-alanine on the Cu{111} single crystal surface was investigated as a model system for interactions between small chiral modifier molecules and close-packed metal surfaces. Synchrotron-based X-ray photoelectron spectroscopy (XPS) and near-edge X-ray absorption fine structure (NEXAFS) spectroscopy are used to determine the chemical state, bond coordination and out-of-plane orientation of the molecule on the surface. Alanine adsorbs in its anionic form at room temperature, whilst at low temperature the overlayer consists of anionic and zwitterionic molecules. NEXAFS spectra exhibit a strong angular dependence of the π ⁎ resonance associated with the carboxylate group, which allows determining the tilt angle of this group with respect to the surface plane (48° ± 2°) at room temperature. Low-energy electron diffraction (LEED) shows a p(2√13x2√13)R13° superstructure with only one domain, which breaks the mirror symmetry of the substrate and, thus, induces global chirality to the surface. Temperature-programmed XPS (TP-XPS) and temperature-programmed desorption (TPD) experiments indicate that the zwitterionic form converts into the anionic species (alaninate) at 293 K. The latter desorbs/decomposes between 435 K and 445 K.

The annual-cycle modulation of meridional asymmetry in ENSO’s atmospheric response and its dependence on ENSO zonal structure

Previous studies documented that a distinct southward shift of central-Pacific low-level wind anomalies occurring during the ENSO decaying phase, is caused by an interaction between the Western Pacific annual cycle and El Niño-Southern Oscillation (ENSO) variability. The present study finds that the meridional movement of the central-Pacific wind anomalies appears only during traditional Eastern-Pacific (or EP) El Niño events rather than in Central-Pacific (CP) El Niño events in which sea surface temperature (SST) anomalies are confined to the central Pacific. The zonal structure of ENSO-related SST anomalies therefore has an important effect on meridional asymmetry in the associated atmospheric response and its modulation by the annual cycle. In contrast to EP El Niño events, the SST anomalies of CP El Niño events extend further west towards to the warm pool region with its climatological warm SSTs. In the warm pool region, relatively small SST anomalies thus are able to excite convection anomalies on both sides of the equator, even with a meridionally asymmetric SST background state. Therefore, almost meridionally symmetric precipitation and wind anomalies are observed over the central Pacific during the decaying phase of CP El Niño events. The SST anomaly pattern of La Niña events is similar to CP El Niño events with a reversed sign. Accordingly, no distinct southward displacement of the atmospheric response occurs over the central Pacific during the La Niña decaying phase. These results have important implications for ENSO climate impacts over East Asia, since the anomalous low-level anticyclone over the western North Pacific is an integral part of the annual cycle-modulated ENSO response.

On the structure and dynamics of Indian monsoon depressions

ERA-Interim reanalysis data from the past 35 years have been used with a newly-developed feature tracking algorithm to identify Indian monsoon depressions originating in or near the Bay of Bengal. These were then rotated, centralised and combined to give a fully three-dimensional 106-depression composite structure – a considerably larger sample than any previous detailed study on monsoon depressions and their structure. Many known features of depression structure are confirmed, particularly the existence of a maximum to the southwest of the centre in rainfall and other fields, and a westward axial tilt in others. Additionally, the depressions are found to have significant asymmetry due to the presence of the Himalayas; a bimodal mid-tropospheric potential vorticity core; a separation into thermally cold- (~–1.5K) and neutral- (~0K) cores near the surface with distinct properties; and that the centre has very large CAPE and very small CIN. Variability as a function of background state has also been explored, with land/coast/sea, diurnal, ENSO, active/break and Indian Ocean Dipole contrasts considered. Depressions are found to be markedly stronger during the active phase of the monsoon, as well as during La Niña. Depressions on land are shown to be more intense and more tightly constrained to the central axis. A detailed schematic diagram of a vertical cross-section through a composite depression is also presented, showing its inherent asymmetric structure.

Observations of urban boundary layer structure during a strong urban heat island event

It has long been known that the urban surface energy balance is different to that of a rural surface, and that heating of the urban surface after sunset gives rise to the Urban Heat Island (UHI). Less well known is how flow and turbulence structure above the urban surface are changed during different phases of the urban boundary layer (UBL). This paper presents new observations above both an urban and rural surface and investigates how much UBL structure deviates from classical behaviour. A 5-day, low wind, cloudless, high pressure period over London, UK, was chosen for analysis, during which there was a strong UHI. Boundary layer evolution for both sites was determined by the diurnal cycle in sensible heat flux, with an extended decay period of approximately 4 h for the convective UBL. This is referred to as the “Urban Convective Island” as the surrounding rural area was already stable at this time. Mixing height magnitude depended on the combination of regional temperature profiles and surface temperature. Given the daytime UHI intensity of 1.5∘C, combined with multiple inversions in the temperature profile, urban and rural mixing heights underwent opposite trends over the period, resulting in a factor of three height difference by the fifth day. Nocturnal jets undergoing inertial oscillations were observed aloft in the urban wind profile as soon as the rural boundary layer became stable: clear jet maxima over the urban surface only emerged once the UBL had become stable. This was due to mixing during the Urban Convective Island reducing shear. Analysis of turbulent moments (variance, skewness and kurtosis) showed “upside-down” boundary layer characteristics on some mornings during initial rapid growth of the convective UBL. During the “Urban Convective Island” phase, turbulence structure still resembled a classical convective boundary layer but with some influence from shear aloft, depending on jet strength. These results demonstrate that appropriate choice of Doppler lidar scan patterns can give detailed profiles of UBL flow. Insights drawn from the observations have implications for accuracy of boundary conditions when simulating urban flow and dispersion, as the UBL is clearly the result of processes driven not only by local surface conditions but also regional atmospheric structure.

Infrared and millimetre-wave scintillometry in the suburban environment – Part 1: Structure parameters

Scintillometry, a form of ground-based remote sensing, provides the capability to estimate surface heat fluxes over scales of a few hundred metres to kilometres. Measurements are spatial averages, making this technique particularly valuable over areas with moderate heterogeneity such as mixed agricultural or urban environments. In this study, we present the structure parameters of temperature and humidity, which can be related to the sensible and latent heat fluxes through similarity theory, for a suburban area in the UK. The fluxes are provided in the second paper of this two-part series. A millimetre-wave scintillometer was combined with an infrared scintillometer along a 5.5 km path over northern Swindon. The pairing of these two wavelengths offers sensitivity to both temperature and humidity fluctuations, and the correlation between wavelengths is also used to retrieve the path-averaged temperature–humidity correlation. Comparison is made with structure parameters calculated from an eddy covariance station located close to the centre of the scintillometer path. The performance of the measurement techniques under different conditions is discussed. Similar behaviour is seen between the two data sets at sub-daily timescales. For the two summer-to-winter periods presented here, similar evolution is displayed across the seasons. A higher vegetation fraction within the scintillometer source area is consistent with the lower Bowen ratio observed (midday Bowen ratio < 1) compared with more built-up areas around the eddy covariance station. The energy partitioning is further explored in the companion paper.

Improved simulation of fire–vegetation interactions in the Land surface Processes and eXchanges dynamic global vegetation model (LPX-Mv1)

The Land surface Processes and eXchanges (LPX) model is a fire-enabled dynamic global vegetation model that performs well globally but has problems representing fire regimes and vegetative mix in savannas. Here we focus on improving the fire module. To improve the representation of ignitions, we introduced a reatment of lightning that allows the fraction of ground strikes to vary spatially and seasonally, realistically partitions strike distribution between wet and dry days, and varies the number of dry days with strikes. Fuel availability and moisture content were improved by implementing decomposition rates specific to individual plant functional types and litter classes, and litter drying rates driven by atmospheric water content. To improve water extraction by grasses, we use realistic plant-specific treatments of deep roots. To improve fire responses, we introduced adaptive bark thickness and post-fire resprouting for tropical and temperate broadleaf trees. All improvements are based on extensive analyses of relevant observational data sets. We test model performance for Australia, first evaluating parameterisations separately and then measuring overall behaviour against standard benchmarks. Changes to the lightning parameterisation produce a more realistic simulation of fires in southeastern and central Australia. Implementation of PFT-specific decomposition rates enhances performance in central Australia. Changes in fuel drying improve fire in northern Australia, while changes in rooting depth produce a more realistic simulation of fuel availability and structure in central and northern Australia. The introduction of adaptive bark thickness and resprouting produces more realistic fire regimes in Australian savannas. We also show that the model simulates biomass recovery rates consistent with observations from several different regions of the world characterised by resprouting vegetation. The new model (LPX-Mv1) produces an improved simulation of observed vegetation composition and mean annual burnt area, by 33 and 18% respectively compared to LPX.

Turbulent structure and scaling of the inertial subrange in a stratocumulus-topped boundary layer observed by a Doppler lidar

The turbulent structure of a stratocumulus-topped marine boundary layer over a 2-day period is observed with a Doppler lidar at Mace Head in Ireland. Using profiles of vertical velocity statistics, the bulk of the mixing is identified as cloud driven. This is supported by the pertinent feature of negative vertical velocity skewness in the sub-cloud layer which extends, on occasion, almost to the surface. Both coupled and decoupled turbulence characteristics are observed. The length and timescales related to the cloud-driven mixing are investigated and shown to provide additional information about the structure and the source of the mixing inside the boundary layer. They are also shown to place constraints on the length of the sampling periods used to derive products, such as the turbulent dissipation rate, from lidar measurements. For this, the maximum wavelengths that belong to the inertial subrange are studied through spectral analysis of the vertical velocity. The maximum wavelength of the inertial subrange in the cloud-driven layer scales relatively well with the corresponding layer depth during pronounced decoupled structure identified from the vertical velocity skewness. However, on many occasions, combining the analysis of the inertial subrange and vertical velocity statistics suggests higher decoupling height than expected from the skewness profiles. Our results show that investigation of the length scales related to the inertial subrange significantly complements the analysis of the vertical velocity statistics and enables a more confident interpretation of complex boundary layer structures using measurements from a Doppler lidar.

Electronic structure of Pd multi-layers on Re(0001): the role of charge transfer

Understanding the origin of the properties of metal-supported metal thin films is important for the rational design of bimetallic catalysts and other applications, but it is generally difficult to separate effects related to strain from those arising from interface interactions. Here we use density functional (DFT) theory to examine the structure and electronic behavior of few-layer palladium films on the rhenium (0001) surface, where there is negligible interfacial strain and therefore other effects can be isolated. Our DFT calculations predict stacking sequences and interlayer separations in excellent agreement with quantitative low-energy electron diffraction experiments. By theoretically simulating the Pd core-level X-ray photoemission spectra (XPS) of the films, we are able to interpret and assign the basic features of both low-resolution and high-resolution XPS measurements. The core levels at the interface shift to more negative energies, rigidly following the shifts in the same direction of the valence d-band center. We demonstrate that the valence band shift at the interface is caused by charge transfer from Re to Pd, which occurs mainly to valence states of hybridized s-p character rather than to the Pd d-band. Since the d-band filling is roughly constant, there is a correlation between the d-band center shift and its bandwidth. The resulting effect of this charge transfer on the valence d-band is thus analogous to the application of a lateral compressive strain on the adlayers. Our analysis suggests that charge transfer should be considered when describing the origin of core and valence band shifts in other metal / metal adlayer systems.

Surface chemistry of alanine on Ni{111}

The adsorption of L-alanine on Ni{111} has been studied as a 10 model of enantioselective heterogeneous catalysts. Synchrotron-based X-ray 11 photoelectron spectroscopy and near-edge X-ray absorption fine structure 12 (NEXAFS) spectroscopy were used to determine the chemical state, bond 13 coordination, and out-of-plane orientation of the molecule on the surface. 14 Alanine adsorbs in anionic and zwitterionic forms between 250 and ≈320 K. 15 NEXAFS spectra exhibit a strong angular dependence of the π* resonance 16 associated with the carboxylate group, which is compatible with two distinct 17 orientations with respect to the surface corresponding to the bidentate and 18 tridentate binding modes. Desorption and decomposition begin together at 19 ≈300 K, with decomposition occurring in a multistep process up to ≈450 K. Comparison with previous studies of amino acid 20 adsorption on metal surfaces shows that this is among the lowest decomposition temperatures found so far and lower than typical 21 temperatures used for hydrogenation reactions where modified Ni catalysts are used.

The atomic structure of low-index surfaces of the intermetallic compound InPd

The intermetallic compound InPd (CsCl type of crystal structure with a broad compositional range) is considered as a candidate catalyst for the steam reforming of methanol. Single crystals of this phase have been grown to study the structure of its three low-index surfaces under ultra-high vacuum conditions, using low energy electron diffraction (LEED), X-ray photoemission spectroscopy (XPS), and scanning tunneling microscopy (STM). During surface preparation, preferential sputtering leads to a depletion of In within the top few layers for all three surfaces. The near-surface regions remain slightly Pd-rich until annealing to ∼580 K. A transition occurs between 580 and 660 K where In segregates towards the surface and the near-surface regions become slightly In-rich above ∼660 K. This transition is accompanied by a sharpening of LEED patterns and formation of flat step-terrace morphology, as observed by STM. Several superstructures have been identified for the different surfaces associated with this process. Annealing to higher temperatures (≥750 K) leads to faceting via thermal etching as shown for the (110) surface, with a bulk In composition close to the In-rich limit of the existence domain of the cubic phase. The Pd-rich InPd(111) is found to be consistent with a Pd-terminated bulk truncation model as shown by dynamical LEED analysis while, after annealing at higher temperature, the In-rich InPd(111) is consistent with an In-terminated bulk truncation, in agreement with density functional theory (DFT) calculations of the relative surface energies. More complex surface structures are observed for the (100) surface. Additionally, individual grains of a polycrystalline sample are characterized by micro-spot XPS and LEED as well as low-energy electron microscopy. Results from both individual grains and “global” measurements are interpreted based on comparison to our single crystals findings, DFT calculations and previous literature.

Using Distributed Temperature Sensing to monitor field scale dynamics of ground surface temperature and related substrate heat flux

We present one of the first studies of the use of Distributed Temperature Sensing (DTS) along fibre-optic cables to purposely monitor spatial and temporal variations in ground surface temperature (GST) and soil temperature, and provide an estimate of the heat flux at the base of the canopy layer and in the soil. Our field site was at a groundwater-fed wet meadow in the Netherlands covered by a canopy layer (between 0-0.5 m thickness) consisting of grass and sedges. At this site, we ran a single cable across the surface in parallel 40 m sections spaced by 2 m, to create a 40×40 m monitoring field for GST. We also buried a short length (≈10 m) of cable to depth of 0.1±0.02 m to measure soil temperature. We monitored the temperature along the entire cable continuously over a two-day period and captured the diurnal course of GST, and how it was affected by rainfall and canopy structure. The diurnal GST range, as observed by the DTS system, varied between 20.94 and 35.08◦C; precipitation events acted to suppress the range of GST. The spatial distribution of GST correlated with canopy vegetation height during both day and night. Using estimates of thermal inertia, combined with a harmonic analysis of GST and soil temperature, substrate and soil-heat fluxes were determined. Our observations demonstrate how the use of DTS shows great promise in better characterising area-average substrate/soil heat flux, their spatiotemporal variability, and how this variability is affected by canopy structure. The DTS system is able to provide a much richer data set than could be obtained from point temperature sensors. Furthermore, substrate heat fluxes derived from GST measurements may be able to provide improved closure of the land surface energy balance in micrometeorological field studies. This will enhance our understanding of how hydrometeorological processes interact with near-surface heat fluxes.