Numerical simulation of baroclinic waves with a parameterized boundary layer

A dry three-dimensional baroclinic life cycle model is used to investigate the role of turbulent fluxes of heat and momentum within the boundary layer on mid-latitude cyclones. Simulations are performed of life cycles for two basic states, both with and without turbulent fluxes. The different basic states produce cyclones with contrasting frontal and mesoscale-flow structures. The analysis focuses on the generation of potential-vorticity (PV) in the boundary layer and its subsequent transport into the free troposphere. The dynamic mechanism through which friction mitigates a barotropic vortex is that of Ekman pumping. This has often been assumed to be also the dominant mechanism for baroclinic developments. The PV framework highlights an additional, baroclinic mechanism. Positive PV is generated baroclinically due to friction to the north-east of a surface low and is transported out of the boundary layer by a cyclonic conveyor belt flow. The result is an anomaly of increased static stability in the lower troposphere which restricts the growth of the baroclinic wave. The reduced coupling between lower and upper levels can be sufficient to change the character of the upper-level evolution of the mature wave. The basic features of the baroclinic damping mechanism are robust for different frontal structures, with and without turbulent heat fluxes, and for the range of surface roughness found over the oceans.

Improved simulation of water vapour and clear-sky radiation using 24 hour forecasts from ERA40

Monthly mean water vapour and clear-sky radiation extracted from the European Centre for Medium Range Weather Forecasts 40-year reanalysis (ERA40) forecasts are assessed using satellite observations and additional reanalysis data. There is a marked improvement in the interannual variability of column-integrated water vapour (CWV) over the oceans when using the 24-hour forecasts compared with the standard 6-hour forecasts products. The spatial distribution of CWV are well simulated by the 6-hour forecasts; using the 24-hour forecasts does not degrade this simulation substantially and in many cases improves on the quality. There is also an improved simulation of clear-sky radiation from the 24-hour forecasts compared with the 6-hour forecasts based on comparison with satellite observations and empirical estimates. Further work is required to assess the quality of water vapour simulation by reanalyses over land regions. Over the oceans, it is recommended that 24-hour forecasts of CWV and clear-sky radiation are used in preference to the standard 6-hour forecast products from ERA40

Observations of the diurnal cycle of outgoing longwave radiation from the Geostationary Earth Radiation Budget instrument

The Geostationary Earth Radiation Budget instrument on Meteosat-8, located over Africa, provides unprecedented temporal sampling (~17 minutes) of the broadband emitted thermal and reflected solar radiances. We analyse the diurnal cycle of the outgoing longwave radiation (OLR) fluxes derived from the thermal radiances for July 2006. Principal component (PC) analysis separates the signals of the surface temperature response to solar heating and of the development of convective clouds. The first two PCs explain most of the OLR variations: PC1 (surface heating) explains 82.3% of the total variance and PC2 (cloud development) explains 12.8% of the variance. Convection is initiated preferentially over mountainous regions and the cloud then advects downstream in the ambient flow. Diurnal variations are much weaker over the oceans, but a coherent signal over the Gulf of Guinea suggests that the cloudiness is modulated by the diurnally varying contrast between the Gulf and the adjacent land mass.

Stratosphere-troposphere transport in a numerical simulation of midlatitude convection

The transport of stratospheric air deep into the troposphere via convection is investigated numerically using the UK Met Office Unified Model. A convective system that formed on 27 June 2004 near southeast England, in the vicinity an upper level potential vorticity anomaly and a lowered tropopause, provides the basis for analysis. Transport is diagnosed using a stratospheric tracer that can either be passed through or withheld from the model’s convective parameterization scheme. Three simulations are performed at increasingly finer resolutions, with horizontal grid lengths of 12, 4, and 1 km. In the 12 and 4 km simulations, tracer is transported deeply into the troposphere by the parameterized convection. In the 1 km simulation, for which the convective parameterization is disengaged, deep transport is still accomplished but with a much smaller magnitude. However, the 1 km simulation resolves stirring along the tropopause that does not exist in the coarser simulations. In all three simulations, the concentration of the deeply transported tracer is small, three orders of magnitude less than that of the shallow transport near the tropopause, most likely because of the efficient dilution of parcels in the lower troposphere.

Linear viscoelasticity from molecular dynamics simulation of entangled polymers

The linear viscoelastic (LVE) spectrum is one of the primary fingerprints of polymer solutions and melts, carrying information about most relaxation processes in the system. Many single chain theories and models start with predicting the LVE spectrum to validate their assumptions. However, until now, no reliable linear stress relaxation data were available from simulations of multichain systems. In this work, we propose a new efficient way to calculate a wide variety of correlation functions and mean-square displacements during simulations without significant additional CPU cost. Using this method, we calculate stress−stress autocorrelation functions for a simple bead−spring model of polymer melt for a wide range of chain lengths, densities, temperatures, and chain stiffnesses. The obtained stress−stress autocorrelation functions were compared with the single chain slip−spring model in order to obtain entanglement related parameters, such as the plateau modulus or the molecular weight between entanglements. Then, the dependence of the plateau modulus on the packing length is discussed. We have also identified three different contributions to the stress relaxation:  bond length relaxation, colloidal and polymeric. Their dependence on the density and the temperature is demonstrated for short unentangled systems without inertia.

A double-Gerdien instrument for simultaneous bipolar air conductivity measurements on balloon platforms

A bipolar air conductivity instrument is described for use with a standard disposable meteorological radiosonde package. It is intended to provide electrical measurements at cloud boundaries, where the ratio of the bipolar air conductivities is affected by the presence of charged particles. The sensors are two identical Gerdien-type electrodes, which, through a voltage decay method, measure positive and negative air conductivities simultaneously. Voltage decay provides a thermally stable approach and a novel low current leakage electrometer switch is described which initiates the decay sequence. The radiosonde supplies power and telemetry, as well as measuring simultaneous meteorological data. A test flight using a tethered balloon determined positive (σ+) and negative (σ−) conductivities of σ+ = 2.77±0.2 fS m−1 and σ− = 2.82±0.2 fS m−1, respectively, at 400 m aloft, with σ+/σ− = 0.98±0.04.

New insights on growth mechanisms of protein clusters at surfaces: an AFM and simulation study

Despite its relevance to a wide range of technological and fundamental areas, a quantitative understanding of protein surface clustering dynamics is often lacking. In inorganic crystal growth, surface clustering of adatoms is well described by diffusion-aggregation models. In such models, the statistical properties of the aggregate arrays often reveal the molecular scale aggregation processes. We investigate the potential of these theories to reveal hitherto hidden facets of protein clustering by carrying out concomitant observations of lysozyme adsorption onto mica surfaces, using atomic force microscopy. and Monte Carlo simulations of cluster nucleation and growth. We find that lysozyme clusters diffuse across the substrate at a rate that varies inversely with size. This result suggests which molecular scale mechanisms are responsible for the mobility of the proteins on the substrate. In addition the surface diffusion coefficient of the monomer can also be extracted from the comparison between experiments and simulations. While concentrating on a model system of lysozyme-on-mica, this 'proof of concept' study successfully demonstrates the potential of our approach to understand and influence more biomedically applicable protein-substrate couples.

Simulation of the Earth's radiation budget by the European Centre for Medium-Range Weather Forecasts 40-year reanalysis (ERA40)

The radiation budget simulated by the European Centre for Medium-Range Weather Forecasts (ECMWF) 40-year reanalysis (ERA40) is evaluated for the period 1979–2001 using independent satellite data and additional model data. This provides information on the quality of the radiation products and indirect evaluation of other aspects of the climate produced by ERA40. The climatology of clear-sky outgoing longwave radiation (OLR) is well captured by ERA40. Underestimations of about 10 W m−2 in clear-sky OLR over tropical convective regions by ERA40 compared to satellite data are substantially reduced when the satellite sampling is taken into account. The climatology of column-integrated water vapor is well simulated by ERA40 compared to satellite data over the ocean, indicating that the simulation of downward clear-sky longwave fluxes at the surface is likely to be good. Clear-sky absorbed solar radiation (ASR) and clear-sky OLR are overestimated by ERA40 over north Africa and high-latitude land regions. The observed interannual changes in low-latitude means are not well reproduced. Using ERA40 to analyze trends and climate feedbacks globally is therefore not recommended. The all-sky radiation budget is poorly simulated by ERA40. OLR is overestimated by around 10 W m−2 over much of the globe. ASR is underestimated by around 30 W m−2 over tropical ocean regions. Away from marine stratocumulus regions, where cloud fraction is underestimated by ERA40, the poor radiation simulation by ERA40 appears to be related to inaccurate radiative properties of cloud rather than inaccurate cloud distributions.

MJO Simulation Diagnostics

The Madden–Julian oscillation (MJO) interacts with and influences a wide range of weather and climate phenomena (e.g., monsoons, ENSO, tropical storms, midlatitude weather), and represents an important, and as yet unexploited, source of predictability at the subseasonal time scale. Despite the important role of the MJO in climate and weather systems, current global circulation models (GCMs) exhibit considerable shortcomings in representing this phenomenon. These shortcomings have been documented in a number of multimodel comparison studies over the last decade. However, diagnosis of model performance has been challenging, and model progress has been difficult to track, because of the lack of a coherent and standardized set of MJO diagnostics. One of the chief objectives of the U.S. Climate Variability and Predictability (CLIVAR) MJO Working Group is the development of observation-based diagnostics for objectively evaluating global model simulations of the MJO in a consistent framework. Motivation for this activity is reviewed, and the intent and justification for a set of diagnostics is provided, along with specification for their calculation, and illustrations of their application. The diagnostics range from relatively simple analyses of variance and correlation to more sophisticated space–time spectral and empirical orthogonal function analyses. These diagnostic techniques are used to detect MJO signals, to construct composite life cycles, to identify associations of MJO activity with the mean state, and to describe interannual variability of the MJO.

Laboratory simulation of terrestrial meteorite weathering using the Bensour (LL6) ordinary chondrite

Laboratory dissolution experiments using the LL6 ordinary chondrite Bensour demonstrate that meteoritic minerals readily react with distilled water at low temperatures, liberating ions into solution and forming reaction products. Three experiments were performed, all for 68 days and at atmospheric fO(2) but using a range of water/rock ratios and different ternperatures. Experiments I and 2 were batch experiments and undertaken at room temperature, whereas in experiment 3, condensed boiling water was dripped onto meteorite subsamples within a Soxhlet extractor. Solutions from experiment 1 were chemically analyzed at the end of the experiment, whereas aliquots were extracted from experiments 2 and 3 for analysis at regular intervals. In all three experiments, a very significant proportion of the Na, Cl, and K within the Bensour subsamples entered solution, demonstrating that chlorapatite and feldspar were especially susceptible to dissolution. Concentrations of Mg, Al, Si, Ca, and Fe in solution were strongly affected by the precipitation of reaction products and Mg and Ca may also have been removed by sorption. Calculations predict saturation of experimental solutions with respect to Al hydroxides, Fe oxides, and Fe (oxy)hydroxides, which would have frequently been accompanied by hydrous aluminosilicates. Some reaction products were identified and include silica, a Mg-rich silicate, Fe oxides, and Fe (oxy)hydroxides. The implications of these results are that even very short periods of subaerial exposure of ordinary chondrites will lead to dissolution of primary minerals and crystallization of weathering products that are likely to include aluminosilicates and silicates, Mg-Ca carbonates, and sulfates in addition to the ubiquitous Fe oxides and (oxy)hydroxides.