Resolution of sharp fronts in the presence of model error in variational data assimilation

We show that the four-dimensional variational data assimilation method (4DVar) can be interpreted as a form of Tikhonov regularization, a very familiar method for solving ill-posed inverse problems. It is known from image restoration problems that L1-norm penalty regularization recovers sharp edges in the image more accurately than Tikhonov, or L2-norm, penalty regularization. We apply this idea from stationary inverse problems to 4DVar, a dynamical inverse problem, and give examples for an L1-norm penalty approach and a mixed total variation (TV) L1–L2-norm penalty approach. For problems with model error where sharp fronts are present and the background and observation error covariances are known, the mixed TV L1–L2-norm penalty performs better than either the L1-norm method or the strong constraint 4DVar (L2-norm)method. A strength of the mixed TV L1–L2-norm regularization is that in the case where a simplified form of the background error covariance matrix is used it produces a much more accurate analysis than 4DVar. The method thus has the potential in numerical weather prediction to overcome operational problems with poorly tuned background error covariance matrices.

The diabatic contour advective semi-Lagrangian model

This article describes a novel algorithmic development extending the contour advective semi-Lagrangian model to include nonconservative effects. The Lagrangian contour representation of finescale tracer fields, such as potential vorticity, allows for conservative, nondiffusive treatment of sharp gradients allowing very high numerical Reynolds numbers. It has been widely employed in accurate geostrophic turbulence and tracer advection simulations. In the present, diabatic version of the model the constraint of conservative dynamics is overcome by including a parallel Eulerian field that absorbs the nonconservative ( diabatic) tendencies. The diabatic buildup in this Eulerian field is limited through regular, controlled transfers of this field to the contour representation. This transfer is done with a fast newly developed contouring algorithm. This model has been implemented for several idealized geometries. In this paper a single-layer doubly periodic geometry is used to demonstrate the validity of the model. The present model converges faster than the analogous semi-Lagrangian models at increased resolutions. At the same nominal spatial resolution the new model is 40 times faster than the analogous semi-Lagrangian model. Results of an orographically forced idealized storm track show nontrivial dependency of storm-track statistics on resolution and on the numerical model employed. If this result is more generally applicable, this may have important consequences for future high-resolution climate modeling.

Chemical composition observed over the mid-Atlantic and the detection of pollution signatures far from source regions

The atmospheric composition of the central North Atlantic region has been sampled using the FAAM BAe146 instrumented aircraft during the Intercontinental Transport of Ozone and Precursors (ITOP) campaign, part of the wider International Consortium for Atmospheric Research on Transport and Transformation (ICARTT). This paper presents an overview of the ITOP campaign. Between late July and early August 2004, twelve flights comprising 72 hours of measurement were made in a region from approximately 20 to 40°W and 33 to 47°N centered on Faial Island, Azores, ranging in altitude from 50 to 9000 m. The vertical profiles of O3 and CO are consistent with previous observations made in this region during 1997 and our knowledge of the seasonal cycles within the region. A cluster analysis technique is used to partition the data set into air mass types with distinct chemical signatures. Six clusters provide a suitable balance between cluster generality and specificity. The clusters are labeled as biomass burning, low level outflow, upper level outflow, moist lower troposphere, marine and upper troposphere. During this summer, boreal forest fire emissions from Alaska and northern Canada were found to provide a major perturbation of tropospheric composition in CO, PAN, organic compounds and aerosol. Anthropogenic influenced air from the continental boundary layer of the USA was clearly observed running above the marine boundary layer right across the mid-Atlantic, retaining high pollution levels in VOCs and sulfate aerosol. Upper level outflow events were found to have far lower sulfate aerosol, resulting from washout on ascent, but much higher PAN associated with the colder temperatures. Lagrangian links with flights of other aircraft over the USA and Europe show that such signatures are maintained many days downwind of emission regions. Some other features of the data set are highlighted, including the strong perturbations to many VOCs and OVOCs in this remote region.

The Convective Storm Initiation Project

The Convective Storm Initiation Project (CSIP) is an international project to understand precisely where, when, and how convective clouds form and develop into showers in the mainly maritime environment of southern England. A major aim of CSIP is to compare the results of the very high resolution Met Office weather forecasting model with detailed observations of the early stages of convective clouds and to use the newly gained understanding to improve the predictions of the model. A large array of ground-based instruments plus two instrumented aircraft, from the U.K. National Centre for Atmospheric Science (NCAS) and the German Institute for Meteorology and Climate Research (IMK), Karlsruhe, were deployed in southern England, over an area centered on the meteorological radars at Chilbolton, during the summers of 2004 and 2005. In addition to a variety of ground-based remote-sensing instruments, numerous rawin-sondes were released at one- to two-hourly intervals from six closely spaced sites. The Met Office weather radar network and Meteosat satellite imagery were used to provide context for the observations made by the instruments deployed during CSIP. This article presents an overview of the CSIP field campaign and examples from CSIP of the types of convective initiation phenomena that are typical in the United Kingdom. It shows the way in which certain kinds of observational data are able to reveal these phenomena and gives an explanation of how the analyses of data from the field campaign will be used in the development of an improved very high resolution NWP model for operational use.

Layer-by-layer electrostatic entrapment of protein molecules on superparamagnetic nanoparticle: new strategy to enhance adsorption capacity and maintain biological activity

There is a recent interest to use inorganic-based magnetic nanoparticles as a vehicle to carry biomolecules for various biophysical applications, but direct attachment of the molecules is known to alter their conformation leading to attenuation in activity. In addition, surface immobilization has been limited to monolayer coverage. It is shown that alternate depositions of negatively charged protein molecules, typically bovine serum albumin (BSA) with a positively charged aminocarbohydrate template such as glycol chitosan (GC) on magnetic iron oxide nanoparticle surface as a colloid, are carried out under pH 7.4. Circular dichroism (CD) clearly reveals that the secondary structure of the entrapped BSA sequential depositions in this manner remains totally unaltered which is in sharp contrast to previous attempts. Probing the binding properties of the entrapped BSA using small molecules (Site I and Site II drug compounds) confirms for the first time the full retention of its biological activity as compared with native BSA, which also implies the ready accessibility of the entrapped protein molecules through the porous overlayers. This work clearly suggests a new method to immobilize and store protein molecules beyond monolayer adsorption on a magnetic nanoparticle surface without much structural alteration. This may find applications in magnetic recoverable enzymes or protein delivery.

Reactive Halogens in the Marine Boundary Layer (RHaMBLe): the tropical North Atlantic experiments

The NERC UK SOLAS-funded Reactive Halogens in the Marine Boundary Layer (RHaMBLe) programme comprised three field experiments. This manuscript presents an overview of the measurements made within the two simultaneous remote experiments conducted in the tropical North Atlantic in May and June 2007. Measurements were made from two mobile and one ground-based platforms. The heavily instrumented cruise D319 on the RRS Discovery from Lisbon, Portugal to São Vicente, Cape Verde and back to Falmouth, UK was used to characterise the spatial distribution of boundary layer components likely to play a role in reactive halogen chemistry. Measurements onboard the ARSF Dornier aircraft were used to allow the observations to be interpreted in the context of their vertical distribution and to confirm the interpretation of atmospheric structure in the vicinity of the Cape Verde islands. Long-term ground-based measurements at the Cape Verde Atmospheric Observatory (CVAO) on São Vicente were supplemented by long-term measurements of reactive halogen species and characterisation of additional trace gas and aerosol species during the intensive experimental period. This paper presents a summary of the measurements made within the RHaMBLe remote experiments and discusses them in their meteorological and chemical context as determined from these three platforms and from additional meteorological analyses. Air always arrived at the CVAO from the North East with a range of air mass origins (European, Atlantic and North American continental). Trace gases were present at stable and fairly low concentrations with the exception of a slight increase in some anthropogenic components in air of North American origin, though NOx mixing ratios during this period remained below 20 pptv. Consistency with these air mass classifications is observed in the time series of soluble gas and aerosol composition measurements, with additional identification of periods of slightly elevated dust concentrations consistent with the trajectories passing over the African continent. The CVAO is shown to be broadly representative of the wider North Atlantic marine boundary layer; measurements of NO, O3 and black carbon from the ship are consistent with a clean Northern Hemisphere marine background. Aerosol composition measurements do not indicate elevated organic material associated with clean marine air. Closer to the African coast, black carbon and NO levels start to increase, indicating greater anthropogenic influence. Lower ozone in this region is possibly associated with the increased levels of measured halocarbons, associated with the nutrient rich waters of the Mauritanian upwelling. Bromide and chloride deficits in coarse mode aerosol at both the CVAO and on D319 and the continuous abundance of inorganic gaseous halogen species at CVAO indicate significant reactive cycling of halogens. Aircraft measurements of O3 and CO show that surface measurements are representative of the entire boundary layer in the vicinity both in diurnal variability and absolute levels. Above the inversion layer similar diurnal behaviour in O3 and CO is observed at lower mixing ratios in the air that had originated from south of Cape Verde, possibly from within the ITCZ. ECMWF calculations on two days indicate very different boundary layer depths and aircraft flights over the ship replicate this, giving confidence in the calculated boundary layer depth.

Annual textural banding in Holocene estuarine silts, Severn Estuary Levels (SW Britain): patterns, cause and implications

A high-resolution textural study has been made of laminated and banded estuarine silts exposed intertidally at representative localities and horizons in the Holocene deposits of the Severn Estuary Levels. The laminae, on a submillimetre to millimetre scale, are sharp-based, graded couplets formed of a lower silty part overlain by a finer-textured clayey element. The centimetre- to decimetre-scale banding is formed of laminae in alternating, gradually intergrading sets of relatively coarse and relative fine-grained examples. At outcrop in the field, the banding is recognizable because the coarse sets prove to be recessive to varying degrees under the influence of weathering and current action. Independent evidence at two localities points toward an annual origin for the banding; at a third it arose during part of what appears to have been a relatively short period. Quantified physical arguments suggest that the textural banding is a response of suspended fine sediment to marked seasonal changes in sea temperature and windiness. The banded silts occur in four distinct stratigraphical contexts and record high deposition rates (order 0.01-0.1 m/yr). Because physical factors determine their textures, the silts potentially afford insights in all contexts into aspects of changing Holocene climatic conditions. In one context, the thickness of the bands points to high (order 0.01-0.1 m/yr) but comparatively short-lived (order 10s-100s yrs) rates of relative water-level rise. In the others, however, the banding has no implications for sea-level behaviour, and simply records gross environmental disequilibrium, for example, the recovery of mudflats/marshes after an erosional episode. Similarly, because on account of their rapid accumulation the banded silts preserve animal and human tracks and trackways especially well, they provide an archive of animal and human behaviour in the area during the Holocene.

Incorporating Legacy Soil pH Databases into Digital Soil Maps

Soil data and reliable soil maps are imperative for environmental management. conservation and policy. Data from historical point surveys, e.g. experiment site data and farmers fields can serve this purpose. However, legacy soil information is not necessarily collected for spatial analysis and mapping such that the data may not have immediately useful geo-references. Methods are required to utilise these historical soil databases so that we can produce quantitative maps of soil propel-ties to assess spatial and temporal trends but also to assess where future sampling is required. This paper discusses two such databases: the Representative Soil Sampling Scheme which has monitored the agricultural soil in England and Wales from 1969 to 2003 (between 400 and 900 bulked soil samples were taken annually from different agricultural fields); and the former State Chemistry Laboratory, Victoria, Australia where between 1973 and 1994 approximately 80,000 soil samples were submitted for analysis by farmers. Previous statistical analyses have been performed using administrative regions (with sharp boundaries) for both databases, which are largely unrelated to natural features. For a more detailed spatial analysis that call be linked to climate and terrain attributes, gradual variation of these soil properties should be described. Geostatistical techniques such as ordinary kriging are suited to this. This paper describes the format of the databases and initial approaches as to how they can be used for digital soil mapping. For this paper we have selected soil pH to illustrate the analyses for both databases.

Catchment-scale modelling of flow and nutrient transport in the Chalk unsaturated zone

The unsaturated zone exerts a major control on the delivery of nutrients to Chalk streams, yet flow and transport processes in this complex, dual-porosity medium have remained controversial. A major challenge arises in characterising these processes, both at the detailed mechanistic level and at an appropriate level for inclusion within catchment-scale models for nutrient management. The lowland catchment research (LOCAR) programme in the UK has provided a unique set of comprehensively instrumented groundwater-dominated catchments. Of these, the Pang and Lambourn, tributaries of the Thames near Reading, have been a particular focus for research into subsurface processes and surface water-groundwater interactions. Data from LOCAR and other sources, along with a new dual permeability numerical model of the Chalk, have been used to explore the relative roles of matrix and fracture flow within the unsaturated zone and resolve conflicting hypotheses of response. From the improved understanding gained through these explorations, a parsimonious conceptualisation of the general response of flow and transport within the Chalk unsaturated zone was formulated. This paper summarises the modelling and data findings of these explorations, and describes the integration of the new simplified unsaturated zone representation with a catchment-scale model of nutrients (INCA), resulting in a new model for catchment-scale flow and transport within Chalk systems: INCA-Chalk. This model is applied to the Lambourn, and results, including hindcast and forecast simulations, are presented. These clearly illustrate the decadal time-scales that need to be considered in the context of nutrient management and the EU Water Framework Directive. (C) 2007 Elsevier B.V. All rights reserved.

Diet and diversity at later medieval Fishergate: The isotopic evidence

We present the results of stable carbon and nitrogen isotope analysis of bone collagen for 155 individuals buried at the Later Medieval (13th to early 16th century AD) Gilbertine priory of St. Andrew, Fishergate in the city of York (UK). The data show significant variation in the consumption of marine foods between males and females as well as between individuals buried in different areas of the priory. Specifically, individuals from the crossing of the church and the cloister garth had consumed significantly less marine protein than those from other locations. Isotope data for four individuals diagnosed with diffuse idiopathic skeletal hyperostosis (DISH) are consistent with a diet rich in animal protein. We also observe that isotopic signals of individuals with perimortem sharp force trauma are unusual in the context of the Fishergate dataset. We discuss possible explanations for these patterns and suggest that there may have been a specialist hospital or a local tradition of burying victims of violent conflict at the priory. The results demonstrate how the integration of archaeological, osteological, and isotopic data can provide novel information about Medieval burial and society.

Long-term change, interannual and intra-seasonal variability in climate and glacier mass balance in the Central Greater Caucasus, Russia

Long-term trends, interannual and intra-seasonal variability in the mass-balance record from Djankuat glacier, central Greater Caucasus, Russia, are related to local climate change, synoptic and large-scale anomalies in atmospheric circulation. A clear warming signal emerged in the central Greater Caucasus in the early 1990s, leading to a strong increase in ablation. In the absence of a compensating change in winter accumulation, the net mass balance of Djankuat has declined. The highest value of seasonal ablation on record was registered in the summer of 2000. At the beginning of the 21st century these trends reversed. Ablation was below average even in the summer of 2003, which was unusually warm in western Europe. Precipitation and winter accumulation were high, allowing for a partial recovery of net mass balance. The interannual variability in the components of mass balance is weakly related to the North Atlantic Oscillation (NAO) and the Scandinavian teleconnection patterns, but there is a clear link with the large-scale circulation anomalies represented by the Rossby pattern. Five synoptic categories have been identified for the ablation season of 2005, revealing a strong separation between components of radiation budget, air temperature and daily melt. Air temperature is the main control over melt. The highest values of daily ablation are related to the strongly positive NAO which forces high net radiation, and to the warm and moist advection from the Black Sea.

The formation of large-scale current sheets within magnetic clouds

Magnetic clouds are a class of interplanetary coronal mass ejections (CME) predominantly characterised by a smooth rotation in the magnetic field direction, indicative of a magnetic flux rope structure. Many magnetic clouds, however, also contain sharp discontinuities within the smoothly varying magnetic field, suggestive of narrow current sheets. In this study we present observations and modelling of magnetic clouds with strong current sheet signatures close to the centre of the apparent flux rope structure. Using an analytical magnetic flux rope model, we demonstrate how such current sheets can form as a result of a cloud’s kinematic propagation from the Sun to the Earth, without any external forces or influences. This model is shown to match observations of four particular magnetic clouds remarkably well. The model predicts that current sheet intensity increases for increasing CME angular extent and decreasing CME radial expansion speed. Assuming such current sheets facilitate magnetic reconnection, the process of current sheet formation could ultimately lead a single flux rope becoming fragmented into multiple flux ropes. This change in topology has consequences for magnetic clouds as barriers to energetic particle propagation.

Magnetic cloud distortion resulting from propagation through a structured solar wind: Models and observations

Numerical simulations of magnetic clouds (MCs) propagating through a structured solar wind suggest that MC-associated magnetic flux ropes are highly distorted by inhomogeneities in the ambient medium. In particular, a solar wind configuration of fast wind from high latitudes and slow wind at low latitudes, common at periods close to solar minimum, should distort the cross section of magnetic clouds into concave-outward structures. This phenomenon has been reported in observations of shock front orientations, but not in the body of magnetic clouds. In this study an analytical magnetic cloud model based upon a kinematically distorted flux rope is modified to simulate propagation through a structured medium. This new model is then used to identify specific time series signatures of the resulting concave-outward flux ropes. In situ observations of three well studied magnetic clouds are examined with comparison to the model, but the expected concave-outward signatures are not present. Indeed, the observations are better described by the convex-outward flux rope model. This may be due to a sharp latitudinal transition from fast to slow wind, resulting in a globally concave-outward flux rope, but with convex-outward signatures on a local scale.

Intercomparison of aircraft instruments on board the C-130 and Falcon-20 over southern Germany during EXPORT2000

In the summer 2000 EXPORT aircraft campaign (European eXport of Precursors and Ozone by long-Range Transport), two comprehensively instrumented research aircraft measuring a variety of chemical species flew wing tip to wing tip for a period of one and a quarter hours. During this interval a comparison was undertaken of the measurements of nitrogen oxide (NO), odd nitrogen species (NOy), carbon monoxide (CO) and ozone (O3). The comparison was performed at two different flight levels, which provided a 10-fold variation in the concentrations of both NO (10 to 1000 parts per trillion by volume (pptv)) and NOy (200 to over 2500 pptv). Large peaks of NO and NOy observed from the Falcon 20, which were at first thought to be from the exhaust of the C-130, were also detected on the 4 channel NOxy instrument aboard the C-130. These peaks were a good indication that both aircraft were in the same air mass and that the Falcon 20 was not in the exhaust plume of the C-130. Correlations and statistical analysis are presented between the instruments used on the two separate aircraft platforms. These were found to be in good agreement giving a high degree of correlation for the ambient air studied. Any deviations from the correlations are accounted for in the estimated inaccuracies of the instruments. These results help to establish that the instruments aboard the separate aircraft are reliably able to measure the corresponding chemical species in the range of conditions sampled and that data collected by both aircraft can be co-ordinated for purposes of interpretation.

Monte Carlo phase diagram for diblock copolymer melts

The phase diagram for diblock copolymer melts is evaluated from lattice-based Monte Carlo simulations using parallel tempering, improving upon earlier simulations that used sequential temperature scans. This new approach locates the order-disorder transition (ODT) far more accurately by the occurrence of a sharp spike in the heat capacity. The present study also performs a more thorough investigation of finite-size effects, which reveals that the gyroid (G) morphology spontaneously forms in place of the perforated-lamellar (PL) phase identified in the earlier study. Nevertheless, there still remains a small region where the PL phase appears to be stable. Interestingly, the lamellar (L) phase next to this region exhibits a small population of transient perforations, which may explain previous scattering experiments suggesting a modulated-lamellar (ML) phase.