Four-dimensional variational data assimilation for high resolution nested models

Four-dimensional variational data assimilation (4D-Var) is used in environmental prediction to estimate the state of a system from measurements. When 4D-Var is applied in the context of high resolution nested models, problems may arise in the representation of spatial scales longer than the domain of the model. In this paper we study how well 4D-Var is able to estimate the whole range of spatial scales present in one-way nested models. Using a model of the one-dimensional advection–diffusion equation we show that small spatial scales that are observed can be captured by a 4D-Var assimilation, but that information in the larger scales may be degraded. We propose a modification to 4D-Var which allows a better representation of these larger scales.

Self-assembled molecular complexes and coordination polymers of CdII, hexamine, and monocarboxylates: structural analysis and theoretical studies of supramolecular interactions

Four new cadmium(II) complexes [Cd-2(bz)(4)(H2O)(4)(mu 2-hmt)]center dot Hbz center dot H2O (1), [Cd-3(bz)(6)(H2O)(6)(mu 2-hmt)(2)]center dot 6H(2)O (2), [Cd(pa)(2)(H2O)(mu(2)-hmt)](n) (3), and {[Cd-3(ac)(6)(H2O)(3)(mu(3)-hmt)(2)]center dot 6H(2)O}(n) (4) with hexamine (hmt) and monocarboxylate ions, benzoate (bz), phenylacetate (pa), or acetate (ac) have been synthesized and characterized structurally. Structure determinations reveal that 1 is dinuclear, 2 is trinuclear, 3 is a one-dimensional (1D) infinite chain, and 4 is a two-dimensional (2D) polymer with fused hexagonal rings consisting of Cd-II and hmt. All the Cd-II atoms in the four complexes (except one CdII in 2) possess seven-coordinate pentagonal bipyramidal geometry with the various chelating bidentate carboxylate groups in equatorial sites. One of the CdII ions in 2, a complex that contains two monodentate carboxylates is in a distorted octahedral environment. The bridging mode of hmt is mu 2- in complexes 1-3 but is mu 3- in complex 4. In all complexes, there are significant numbers of H-bonds, C-H/pi, and pi-pi interactions which play crucial roles in forming the supramolecular networks. The importance of the noncovalent interactions in terms of energies and geometries has been analyzed using high level ab initio calculations. The effect of the cadmium coordinated to hmt on the energetic features of the C-H/pi interaction is analyzed. Finally, the interplay between C-H/pi and pi-pi interactions observed in the crystal structure of 3 is also studied.

A new floating model level scheme for the assimilation of boundary layer top inversions: the univariate assimilation of temperature.

The assimilation of observations with a forecast is often heavily influenced by the description of the error covariances associated with the forecast. When a temperature inversion is present at the top of the boundary layer (BL), a significant part of the forecast error may be described as a vertical positional error (as opposed to amplitude error normally dealt with in data assimilation). In these cases, failing to account for positional error explicitly is shown t o r esult in an analysis for which the inversion structure is erroneously weakened and degraded. In this article, a new assimilation scheme is proposed to explicitly include the positional error associated with an inversion. This is done through the introduction of an extra control variable to allow position errors in the a priori to be treated simultaneously with the usual amplitude errors. This new scheme, referred to as the ‘floating BL scheme’, is applied to the one-dimensional (vertical) variational assimilation of temperature. The floating BL scheme is tested with a series of idealised experiments a nd with real data from radiosondes. For each idealised experiment, the floating BL scheme gives an analysis which has the inversion structure and position in agreement with the truth, and outperforms the a ssimilation which accounts only for forecast a mplitude error. When the floating BL scheme is used to assimilate a l arge sample of radiosonde data, its ability to give an analysis with an inversion height in better agreement with that observed is confirmed. However, it is found that the use of Gaussian statistics is an inappropriate description o f t he error statistics o f t he extra c ontrol variable. This problem is alleviated by incorporating a non-Gaussian description of the new control variable in the new scheme. Anticipated challenges in implementing the scheme operationally are discussed towards the end of the article.

Mixed copper, silver, and gold cyanides, (MxM′1–x)CN: tailoring chain structures to influence physical properties

Binary mixed-metal variants of the one-dimensional MCN compounds (M = Cu, Ag, and Au) have been prepared and characterized using powder X-ray diffraction, vibrational spectroscopy, and total neutron diffraction. A solid solution with the AgCN structure exists in the (CuxAg1–x)CN system over the range (0 ≤ x ≤ 1). Line phases with compositions (Cu1/2Au1/2)CN, (Cu7/12Au5/12)CN, (Cu2/3Au1/3)CN, and (Ag1/2Au1/2)CN, all of which have the AuCN structure, are found in the gold-containing systems. Infrared and Raman spectroscopies show that complete ordering of the type [M–C≡N–M′–N≡C−]n occurs only in (Cu1/2Au1/2)CN and (Ag1/2Au1/2)CN. The sense of the cyanide bonding was determined by total neutron diffraction to be [Ag–NC–Au–CN−]n in (Ag1/2Au1/2)CN and [Cu–NC–Au–CN−]n in (Cu1/2Au1/2)CN. In contrast, in (Cu0.50Ag0.50)CN, metal ordering is incomplete, and strict alternation of metals does not occur. However, there is a distinct preference (85%) for the N end of the cyanide ligand to be bonded to copper and for Ag–CN–Cu links to predominate. Contrary to expectation, aurophilic bonding does not appear to be the controlling factor which leads to (Cu1/2Au1/2)CN and (Ag1/2Au1/2)CN adopting the AuCN structure. The diffuse reflectance, photoluminescence, and 1-D negative thermal expansion (NTE) behaviors of all three systems are reported and compared with those of the parent cyanide compounds. The photophysical properties are strongly influenced both by the composition of the individual chains and by how such chains pack together. The NTE behavior is also controlled by structure type: the gold-containing mixed-metal cyanides with the AuCN structure show the smallest contraction along the chain length on heating.

Structural variations in self-assembled coordination complexes of Zn(II) with hexamethylenetetramine and isomeric 2-, 3- and 4-nitrobenzoates

Three new zinc(II)-hexamethylenetetramine (hmt) complexes [Zn-2(4-nbz)(4)(mu(2)-hmt)(OH2)(hmt)] (1). [Zn-2(2-nbz)(4)(mu(2)-hmt)(2)](n) (2) and [Zn-3(3-nbz)(4)(mu(2)-hmt)(mu(2)-OH)(mu(3)-OH)](n) (3) with three isomeric nitrobenzoate, [4-nbz = 4-nitrobenzoate, 2-nbz = 2-nitrobenzoate and 3-nbz = 3-nitrobenzoate] have been synthesized and structurally characterized by X-ray crystallography. Their identities have also been established by elemental analysis: IR, NMR, UV-Vis and mass spectral studies. 1 is a dinuclear complex formed by bridging hmt with mu(2) coordinating mode. The geometry around the Zn centers in 1 is distorted tetrahedral. Paddle-wheel centrosymmetric Zn-2(2-nbz)(4) units of complex 2 are interconnected by mu(2)-hmt forming a one-dimensional chain with square-pyramidal geometries around the Zn centers. Compound 3 contains a mu(2)/mu(3)-hydroxido and mu(2)-hmt bridged 1D chain. In this complex, varied geometries around the Zn centers are observed viz, tetrahedral, square pyramidal and trigonal bipyramidal. Various weak forces, i.e. lone pair-pi, pi-pi and CH-pi interactions, play a key role in stabilizing the observed structures for complexes 1,2 and 3. This series of complexes demonstrates that although the nitro group does not coordinate to the metal center, its presence at the 2-, 3- or 4-position of the phenyl ring has a striking effect on the dimensionality as well as the structure of the resulted coordination polymers, probably due to the participation of the nitro group in 1.p.center dot center dot center dot pi and/or C-H center dot center dot center dot pi interactions.

A model of the three-dimensional evolution of Arctic melt ponds on first-year and multiyear sea ice

During winter the ocean surface in polar regions freezes over to form sea ice. In the summer the upper layers of sea ice and snow melts producing meltwater that accumulates in Arctic melt ponds on the surface of sea ice. An accurate estimate of the fraction of the sea ice surface covered in melt ponds is essential for a realistic estimate of the albedo for global climate models. We present a melt-pond–sea-ice model that simulates the three-dimensional evolution of melt ponds on an Arctic sea ice surface. The advancements of this model compared to previous models are the inclusion of snow topography; meltwater transport rates are calculated from hydraulic gradients and ice permeability; and the incorporation of a detailed one-dimensional, thermodynamic radiative balance. Results of model runs simulating first-year and multiyear sea ice are presented. Model results show good agreement with observations, with duration of pond coverage, pond area, and ice ablation comparing well for both the first-year ice and multiyear ice cases. We investigate the sensitivity of the melt pond cover to changes in ice topography, snow topography, and vertical ice permeability. Snow was found to have an important impact mainly at the start of the melt season, whereas initial ice topography strongly controlled pond size and pond fraction throughout the melt season. A reduction in ice permeability allowed surface flooding of relatively flat, first-year ice but had little impact on the pond coverage of rougher, multiyear ice. We discuss our results, including model shortcomings and areas of experimental uncertainty.

Groundwater residence time distributions in peatlands: implications for peat decomposition and accumulation

Peat soils consist of poorly decomposed plant detritus, preserved by low decay rates, and deep peat deposits are globally significant stores in the carbon cycle. High water tables and low soil temperatures are commonly held to be the primary reasons for low peat decay rates. However, recent studies suggest a thermodynamic limit to peat decay, whereby the slow turnover of peat soil pore water may lead to high concentrations of phenols and dissolved inorganic carbon. In sufficient concentrations, these chemicals may slow or even halt microbial respiration, providing a negative feedback to peat decay. We document the analysis of a simple, one-dimensional theoretical model of peatland pore water residence time distributions (RTDs). The model suggests that broader, thicker peatlands may be more resilient to rapid decay caused by climate change because of slow pore water turnover in deep layers. Even shallow peat deposits may also be resilient to rapid decay if rainfall rates are low. However, the model suggests that even thick peatlands may be vulnerable to rapid decay under prolonged high rainfall rates, which may act to flush pore water with fresh rainwater. We also used the model to illustrate a particular limitation of the diplotelmic (i.e., acrotelm and catotelm) model of peatland structure. Model peatlands of contrasting hydraulic structure exhibited identical water tables but contrasting RTDs. These scenarios would be treated identically by diplotelmic models, although the thermodynamic limit suggests contrasting decay regimes. We therefore conclude that the diplotelmic model be discarded in favor of model schemes that consider continuous variation in peat properties and processes.

Estimating interchannel observation-error correlations for IASI radiance data in the Met Office system

The optimal utilisation of hyper-spectral satellite observations in numerical weather prediction is often inhibited by incorrectly assuming independent interchannel observation errors. However, in order to represent these observation-error covariance structures, an accurate knowledge of the true variances and correlations is needed. This structure is likely to vary with observation type and assimilation system. The work in this article presents the initial results for the estimation of IASI interchannel observation-error correlations when the data are processed in the Met Office one-dimensional (1D-Var) and four-dimensional (4D-Var) variational assimilation systems. The method used to calculate the observation errors is a post-analysis diagnostic which utilises the background and analysis departures from the two systems. The results show significant differences in the source and structure of the observation errors when processed in the two different assimilation systems, but also highlight some common features. When the observations are processed in 1D-Var, the diagnosed error variances are approximately half the size of the error variances used in the current operational system and are very close in size to the instrument noise, suggesting that this is the main source of error. The errors contain no consistent correlations, with the exception of a handful of spectrally close channels. When the observations are processed in 4D-Var, we again find that the observation errors are being overestimated operationally, but the overestimation is significantly larger for many channels. In contrast to 1D-Var, the diagnosed error variances are often larger than the instrument noise in 4D-Var. It is postulated that horizontal errors of representation, not seen in 1D-Var, are a significant contributor to the overall error here. Finally, observation errors diagnosed from 4D-Var are found to contain strong, consistent correlation structures for channels sensitive to water vapour and surface properties.

Computer simulations of the structure of colloidal ferrofluids

The structure of a ferrofluid under the influence of an external magnetic field is expected to become anisotropic due to the alignment of the dipoles into the direction of the external field, and subsequently to the formation of particle chains due to the attractive head to tail orientations of the ferrofluid particles. Knowledge about the structure of a colloidal ferrofluid can be inferred from scattering data via the measurement of structure factors. We have used molecular-dynamics simulations to investigate the structure of both monodispersed and polydispersed ferrofluids. The results for the isotropic structure factor for monodispersed samples are similar to previous data by Camp and Patey that were obtained using an alternative Monte Carlo simulation technique, but in a different parameter region. Here we look in addition at bidispersed samples and compute the anisotropic structure factor by projecting the q vector onto the XY and XZ planes separately, when the magnetic field was applied along the z axis. We observe that the XY- plane structure factor as well as the pair distribution functions are quite different from those obtained for the XZ plane. Further, the two- dimensional structure factor patterns are investigated for both monodispersed and bidispersed samples under different conditions. In addition, we look at the scaling exponents of structure factors. Our results should be of value to interpret scattering data on ferrofluids obtained under the influence of an external field.

Joint multiple dictionary learning for tensor sparse coding

Traditional dictionary learning algorithms are used for finding a sparse representation on high dimensional data by transforming samples into a one-dimensional (1D) vector. This 1D model loses the inherent spatial structure property of data. An alternative solution is to employ Tensor Decomposition for dictionary learning on their original structural form —a tensor— by learning multiple dictionaries along each mode and the corresponding sparse representation in respect to the Kronecker product of these dictionaries. To learn tensor dictionaries along each mode, all the existing methods update each dictionary iteratively in an alternating manner. Because atoms from each mode dictionary jointly make contributions to the sparsity of tensor, existing works ignore atoms correlations between different mode dictionaries by treating each mode dictionary independently. In this paper, we propose a joint multiple dictionary learning method for tensor sparse coding, which explores atom correlations for sparse representation and updates multiple atoms from each mode dictionary simultaneously. In this algorithm, the Frequent-Pattern Tree (FP-tree) mining algorithm is employed to exploit frequent atom patterns in the sparse representation. Inspired by the idea of K-SVD, we develop a new dictionary update method that jointly updates elements in each pattern. Experimental results demonstrate our method outperforms other tensor based dictionary learning algorithms.

Polymorphism and optical properties in [NH4][InSe2]

The solvothermal synthesis and characterization of two indium selenides with stoichiometry [NH4][InSe2] is described. Yellow [NH4][InSe2] (1), which exhibits a layered structure, was initially prepared in an aqueous solution of trans-1,4-diaminocyclohexane, and subsequently using a concentrated ammonia solution. A red polymorph of one-dimensional character, [NH4][InSe2] (2), was obtained using 3,5-dimethylpyridine as solvent. [NH4][InSe2] (1) crystallizes in the non-centrosymmetric space group Cc (a=11.5147(6), b=11.3242(6), c=15.9969(9) Å and β=100.354(3)°). The structural motif of the layers is the In4Se10 adamantane unit, composed of four corner-linked InSe4 tetrahedra. These units are linked by their corners, forming [InSe2]− layers which are stacked back to back along the c-direction, and interspaced by [NH4]+cations. The one-dimensional polymorph, (2), crystallizes in the tetragonal space group, I4/mcm (a=8.2519(16), c=6.9059 (14) Å). This structure contains infinite chains of edge-sharing InSe4 tetrahedra separated by [NH4]+ cations.

A systematic method of parameterisation estimation using data assimilation

In numerical weather prediction, parameterisations are used to simulate missing physics in the model. These can be due to a lack of scientific understanding or a lack of computing power available to address all the known physical processes. Parameterisations are sources of large uncertainty in a model as parameter values used in these parameterisations cannot be measured directly and hence are often not well known; and the parameterisations themselves are also approximations of the processes present in the true atmosphere. Whilst there are many efficient and effective methods for combined state/parameter estimation in data assimilation (DA), such as state augmentation, these are not effective at estimating the structure of parameterisations. A new method of parameterisation estimation is proposed that uses sequential DA methods to estimate errors in the numerical models at each space-time point for each model equation. These errors are then fitted to pre-determined functional forms of missing physics or parameterisations that are based upon prior information. We applied the method to a one-dimensional advection model with additive model error, and it is shown that the method can accurately estimate parameterisations, with consistent error estimates. Furthermore, it is shown how the method depends on the quality of the DA results. The results indicate that this new method is a powerful tool in systematic model improvement.