Hyperresolution global land surface modeling: meeting a grand challenge for monitoring Earth's terrestrial water

Monitoring Earth's terrestrial water conditions is critically important to many hydrological applications such as global food production; assessing water resources sustainability; and flood, drought, and climate change prediction. These needs have motivated the development of pilot monitoring and prediction systems for terrestrial hydrologic and vegetative states, but to date only at the rather coarse spatial resolutions (∼10–100 km) over continental to global domains. Adequately addressing critical water cycle science questions and applications requires systems that are implemented globally at much higher resolutions, on the order of 1 km, resolutions referred to as hyperresolution in the context of global land surface models. This opinion paper sets forth the needs and benefits for a system that would monitor and predict the Earth's terrestrial water, energy, and biogeochemical cycles. We discuss six major challenges in developing a system: improved representation of surface‐subsurface interactions due to fine‐scale topography and vegetation; improved representation of land‐atmospheric interactions and resulting spatial information on soil moisture and evapotranspiration; inclusion of water quality as part of the biogeochemical cycle; representation of human impacts from water management; utilizing massively parallel computer systems and recent computational advances in solving hyperresolution models that will have up to 109 unknowns; and developing the required in situ and remote sensing global data sets. We deem the development of a global hyperresolution model for monitoring the terrestrial water, energy, and biogeochemical cycles a “grand challenge” to the community, and we call upon the international hydrologic community and the hydrological science support infrastructure to endorse the effort.

Simulating dynamic crop growth with an adapted land surface model – JULES-SUCROS: model development and validation

The increasing demand for ecosystem services, in conjunction with climate change, is expected to signif- icantly alter terrestrial ecosystems. In order to evaluate the sustainability of land and water resources, there is a need for a better understanding of the relationships between crop production, land surface characteristics and the energy and water cycles. These relationships are analysed using the Joint UK Land Environment Simulator (JULES). JULES includes the full hydrological cycle and vegetation effects on the energy, water, and carbon fluxes. However, this model currently only simulates land surface processes in natural ecosystems. An adapted version of JULES for agricultural ecosystems, called JULES-SUCROS has therefore been developed. In addition to overall model improvements, JULES-SUCROS includes a dynamic crop growth structure that fully fits within and builds upon the biogeochemical modelling framework for natural vegetation. Specific agro-ecosystem features such as the development of yield-bearing organs and the phenological cycle from sowing till harvest have been included in the model. This paper describes the structure of JULES-SUCROS and evaluates the fluxes simulated with this model against FLUXNET measurements at 6 European sites. We show that JULES-SUCROS significantly improves the correlation between simulated and observed fluxes over cropland and captures well the spatial and temporal vari- ability of the growth conditions in Europe. Simulations with JULES-SUCROS highlight the importance of vegetation structure and phenology, and the impact they have on land–atmosphere interactions.

Simple model of adsorption on external surface of carbon nanotubes: a new analytical approach basing on molecular simulation data

Nitrogen adsorption on carbon nanotubes is wide- ly studied because nitrogen adsorption isotherm measurement is a standard method applied for porosity characterization. A further reason is that carbon nanotubes are potential adsorbents for separation of nitrogen from oxygen in air. The study presented here describes the results of GCMC simulations of nitrogen (three site model) adsorption on single and multi walled closed nanotubes. The results obtained are described by a new adsorption isotherm model proposed in this study. The model can be treated as the tube analogue of the GAB isotherm taking into account the lateral adsorbate-adsorbate interactions. We show that the model describes the simulated data satisfactorily. Next this new approach is applied for a description of experimental data measured on different commercially available (and characterized using HRTEM) carbon nanotubes. We show that generally a quite good fit is observed and therefore it is suggested that the observed mechanism of adsorption in the studied materials is mainly determined by adsorption on tubes separated at large distances, so the tubes behave almost independently.

Can carbon surface oxidation shift the pore size distribution curve calculated from Ar, N(2) and CO(2) adsorption isotherms? Simulation results for a realistic carbon model

Using the virtual porous carbon model proposed by Harris et al, we study the effect of carbon surface oxidation on the pore size distribution (PSD) curve determined from simulated Ar, N(2) and CO(2) isotherms. It is assumed that surface oxidation is not destructive for the carbon skeleton, and that all pores are accessible for studied molecules (i.e., only the effect of the change of surface chemical composition is studied). The results obtained show two important things, i.e., oxidation of the carbon surface very slightly changes the absolute porosity (calculated from the geometric method of Bhattacharya and Gubbins (BG)); however, PSD curves calculated from simulated isotherms are to a greater or lesser extent affected by the presence of surface oxides. The most reliable results are obtained from Ar adsorption data. Not only is adsorption of this adsorbate practically independent from the presence of surface oxides, but, more importantly, for this molecule one can apply the slit-like model of pores as the first approach to recover the average pore diameter of a real carbon structure. For nitrogen, the effect of carbon surface chemical composition is observed due to the quadrupole moment of this molecule, and this effect shifts the PSD curves compared to Ar. The largest differences are seen for CO2, and it is clearly demonstrated that the PSD curves obtained from adsorption isotherms of this molecule contain artificial peaks and the average pore diameter is strongly influenced by the presence of electrostatic adsorbate-adsorbate as well as adsorbate-adsorbent interactions.

The importance of attractive three-point interaction in enantioselective surface chemistry: stereospecific adsorption of serine on the intrinsically chiral Cu{531} surface

Both enantiomers of serine adsorb on the intrinsically chiral Cu{531} surface in two different adsorption geometries, depending on the coverage. At saturation, substrate bonds are formed through the two oxygen atoms of the carboxylate group and the amino group (μ3 coordination), whereas at lower coverage, an additional bond is formed through the deprotonated β−OH group (μ4 coordination). The latter adsorption geometry involves substrate bonds through three side groups of the chiral center, respectively, which leads to significantly larger enantiomeric differences in adsorption geometries and energies compared to the μ3 coordination, which involves only two side groups. This relatively simple model system demonstrates, in direct comparison, that attractive interactions of three side groups with the substrate are much more effective in inducing strong enantiomeric differences in heterogeneous chiral catalyst systems than hydrogen bonds or repulsive interactions.

Influence of fermentation conditions on the surface properties and adhesion of Lactobacillus rhamnosus GG

Background: The surface properties of probiotic bacteria influence to a large extent their interactions within the gut ecosystem. There is limited amount of information on the effect of the production process on the surface properties of probiotic lactobacilli in relation to the mechanisms of their adhesion to the gastrointestinal mucosa. The aim of this work was to investigate the effect of the fermentation pH and temperature on the surface properties and adhesion ability to Caco-2 cells of the probiotic strain Lactobacillus rhamnosus GG. Results: The cells were grown at pH 5, 5.5, 6 (temperature 37 °C) and at pH 6.5 (temperature 25 °C, 30 °C and 37 °C), and their surfaces analysed by X-ray photoelectron spectrometry (XPS), Fourier transform infrared spectroscopy (FT-IR) and gel-based proteomics. The results indicated that for all the fermentation conditions, with the exception of pH 5, a higher nitrogen to carbon ratio and a lower phosphate content was observed at the surface of the bacteria, which resulted in a lower surface hydrophobicity and reduced adhesion levels to Caco-2 cells as compared to the control fermentation (pH 6.5, 37 oC). A number of adhesive proteins, which have been suggested in previous published works to take part in the adhesion of bacteria to the human gastrointestinal tract, were identified by proteomic analysis, with no significant differences between samples however. Conclusions: The temperature and the pH of the fermentation influenced the surface composition, hydrophobicity and the levels of adhesion of L. rhamnosus GG to Caco-2 cells. It was deduced from the data that a protein rich surface reduced the adhesion ability of the cells.

Ice formation and development in aged, wintertime cumulus over the UK: observations and modelling

In situ high resolution aircraft measurements of cloud microphysical properties were made in coordination with ground based remote sensing observations of a line of small cumulus clouds, using Radar and Lidar, as part of the Aerosol Properties, PRocesses And InfluenceS on the Earth's climate (APPRAISE) project. A narrow but extensive line (~100 km long) of shallow convective clouds over the southern UK was studied. Cloud top temperatures were observed to be higher than −8 °C, but the clouds were seen to consist of supercooled droplets and varying concentrations of ice particles. No ice particles were observed to be falling into the cloud tops from above. Current parameterisations of ice nuclei (IN) numbers predict too few particles will be active as ice nuclei to account for ice particle concentrations at the observed, near cloud top, temperatures (−7.5 °C). The role of mineral dust particles, consistent with concentrations observed near the surface, acting as high temperature IN is considered important in this case. It was found that very high concentrations of ice particles (up to 100 L−1) could be produced by secondary ice particle production providing the observed small amount of primary ice (about 0.01 L−1) was present to initiate it. This emphasises the need to understand primary ice formation in slightly supercooled clouds. It is shown using simple calculations that the Hallett-Mossop process (HM) is the likely source of the secondary ice. Model simulations of the case study were performed with the Aerosol Cloud and Precipitation Interactions Model (ACPIM). These parcel model investigations confirmed the HM process to be a very important mechanism for producing the observed high ice concentrations. A key step in generating the high concentrations was the process of collision and coalescence of rain drops, which once formed fell rapidly through the cloud, collecting ice particles which caused them to freeze and form instant large riming particles. The broadening of the droplet size-distribution by collision-coalescence was, therefore, a vital step in this process as this was required to generate the large number of ice crystals observed in the time available. Simulations were also performed with the WRF (Weather, Research and Forecasting) model. The results showed that while HM does act to increase the mass and number concentration of ice particles in these model simulations it was not found to be critical for the formation of precipitation. However, the WRF simulations produced a cloud top that was too cold and this, combined with the assumption of continual replenishing of ice nuclei removed by ice crystal formation, resulted in too many ice crystals forming by primary nucleation compared to the observations and parcel modelling.

The role of protein hydrophobicity in thionin–phospholipid interactions: a comparison of α1 and α2-purothionin adsorbed anionic phospholipid monolayers

The plant defence proteins α1- and α2-purothionin (Pth) are type 1 thionins from common wheat (Triticum aestivum). These highly homologous proteins possess characteristics common amongst antimicrobial peptides and proteins, that is, cationic charge, amphiphilicity and hydrophobicity. Both α1- and α2-Pth possess the same net charge, but differ in relative hydrophobicity as determined by C18 reversed phase HPLC. Brewster angle microscopy, X-ray and neutron reflectometry, external reflection FTIR and associated surface pressure measurements demonstrated that α1 and α2-Pth interact strongly with condensed phase 1,2-dipalmitoyl-sn-glycero-3-phospho-(1'-rac-glycerol) (DPPG) monolayers at the air/liquid interface. Both thionins disrupted the in-plane structure of the anionic phospholipid monolayer, removing lipid during this process and both penetrated the lipid monolayer in addition to adsorbing as a single protein layer to the lipid head-group. However, analysis of the interfacial structures revealed that the α2-Pth showed faster disruption of the lipid film and removed more phospholipid (12%) from the interface than α1-Pth. Correlating the protein properties and lipid binding activity suggests that hydrophobicity plays a key role in the membrane lipid removal activity of thionins.

Estimating root–soil contact from 3D X-ray microtomographs

Adequate contact with the soil is essential for water and nutrient adsorption by plant roots, but the determination of root–soil contact is a challenging task because it is difficult to visualize roots in situ and quantify their interactions with the soil at the scale of micrometres. A method to determine root–soil contact using X-ray microtomography was developed. Contact areas were determined from 3D volumetric images using segmentation and iso-surface determination tools. The accuracy of the method was tested with physical model systems of contact between two objects (phantoms). Volumes, surface areas and contact areas calculated from the measured phantoms were compared with those estimated from image analysis. The volume was accurate to within 0.3%, the surface area to within 2–4%, and the contact area to within 2.5%. Maize and lupin roots were grown in soil (<2 mm) and vermiculite at matric potentials of −0.03 and −1.6 MPa and in aggregate fractions of 4–2, 2–1, 1–0.5 and < 0.5 mm at a matric potential of −0.03 MPa. The contact of the roots with their growth medium was determined from 3D volumetric images. Macroporosity (>70 µm) of the soil sieved to different aggregate fractions was calculated from binarized data. Root-soil contact was greater in soil than in vermiculite and increased with decreasing aggregate or particle size. The differences in root–soil contact could not be explained solely by the decrease in porosity with decreasing aggregate size but may also result from changes in particle and aggregate packing around the root.

Analysis of the seasonal cycle within the first international urban land-surface model comparison

A number of urban land-surface models have been developed in recent years to satisfy the growing requirements for urban weather and climate interactions and prediction. These models vary considerably in their complexity and the processes that they represent. Although the models have been evaluated, the observational datasets have typically been of short duration and so are not suitable to assess the performance over the seasonal cycle. The First International Urban Land-Surface Model comparison used an observational dataset that spanned a period greater than a year, which enables an analysis over the seasonal cycle, whilst the variety of models that took part in the comparison allows the analysis to include a full range of model complexity. The results show that, in general, urban models do capture the seasonal cycle for each of the surface fluxes, but have larger errors in the summer months than in the winter. The net all-wave radiation has the smallest errors at all times of the year but with a negative bias. The latent heat flux and the net storage heat flux are also underestimated, whereas the sensible heat flux generally has a positive bias throughout the seasonal cycle. A representation of vegetation is a necessary, but not sufficient, condition for modelling the latent heat flux and associated sensible heat flux at all times of the year. Models that include a temporal variation in anthropogenic heat flux show some increased skill in the sensible heat flux at night during the winter, although their daytime values are consistently overestimated at all times of the year. Models that use the net all-wave radiation to determine the net storage heat flux have the best agreement with observed values of this flux during the daytime in summer, but perform worse during the winter months. The latter could result from a bias of summer periods in the observational datasets used to derive the relations with net all-wave radiation. Apart from these models, all of the other model categories considered in the analysis result in a mean net storage heat flux that is close to zero throughout the seasonal cycle, which is not seen in the observations. Models with a simple treatment of the physical processes generally perform at least as well as models with greater complexity.

Sensitivity of direct radiative forcing by mineral dust to particle characteristics

Airborne dust affects the Earth's energy balance — an impact that is measured in terms of the implied change in net radiation (or radiative forcing, in W m-2) at the top of the atmosphere. There remains considerable uncertainty in the magnitude and sign of direct forcing by airborne dust under current climate. Much of this uncertainty stems from simplified assumptions about mineral dust-particle size, composition and shape, which are applied in remote sensing retrievals of dust characteristics and dust-cycle models. Improved estimates of direct radiative forcing by dust will require improved characterization of the spatial variability in particle characteristics to provide reliable information dust optical properties. This includes constraints on: (1) particle-size distribution, including discrimination of particle subpopulations and quantification of the amount of dust in the sub-10 µm to <0.1 µm mass fraction; (2) particle composition, specifically the abundance of iron oxides, and whether particles consist of single or multi-mineral grains; (3) particle shape, including degree of sphericity and surface roughness, as a function of size and mineralogy; and (4) the degree to which dust particles are aggregated together. The use of techniques that measure the size, composition and shape of individual particles will provide a better basis for optical modelling.

Key conclusions of the first international urban land surface model comparison project

he first international urban land surface model comparison was designed to identify three aspects of the urban surface-atmosphere interactions: (1) the dominant physical processes, (2) the level of complexity required to model these, and 3) the parameter requirements for such a model. Offline simulations from 32 land surface schemes, with varying complexity, contributed to the comparison. Model results were analysed within a framework of physical classifications and over four stages. The results show that the following are important urban processes; (i) multiple reflections of shortwave radiation within street canyons, (ii) reduction in the amount of visible sky from within the canyon, which impacts on the net long-wave radiation, iii) the contrast in surface temperatures between building roofs and street canyons, and (iv) evaporation from vegetation. Models that use an appropriate bulk albedo based on multiple solar reflections, represent building roof surfaces separately from street canyons and include a representation of vegetation demonstrate more skill, but require parameter information on the albedo, height of the buildings relative to the width of the streets (height to width ratio), the fraction of building roofs compared to street canyons from a plan view (plan area fraction) and the fraction of the surface that is vegetated. These results, whilst based on a single site and less than 18 months of data, have implications for the future design of urban land surface models, the data that need to be measured in urban observational campaigns, and what needs to be included in initiatives for regional and global parameter databases.

On the barrier properties of the cornea: a microscopy study of the penetration of fluorescently labeled nanoparticles, polymers, and sodium fluorescein

Overcoming the natural defensive barrier functions of the eye remains one of the greatest challenges of ocular drug delivery. Cornea is a chemical and mechanical barrier preventing the passage of any foreign bodies including drugs into the eye, but the factors limiting penetration of permeants and nanoparticulate drug delivery systems through the cornea are still not fully understood. In this study, we investigate these barrier properties of the cornea using thiolated and PEGylated (750 and 5000 Da) nanoparticles, sodium fluorescein, and two linear polymers (dextran and polyethylene glycol). Experiments used intact bovine cornea in addition to bovine cornea de-epithelialized or tissues pretreated with cyclodextrin. It was shown that corneal epithelium is the major barrier for permeation; pretreatment of the cornea with β-cyclodextrin provides higher permeation of low molecular weight compounds, such as sodium fluorescein, but does not enhance penetration of nanoparticles and larger molecules. Studying penetration of thiolated and PEGylated (750 and 5000 Da) nanoparticles into the de-epithelialized ocular tissue revealed that interactions between corneal surface and thiol groups of nanoparticles were more significant determinants of penetration than particle size (for the sizes used here). PEGylation with polyethylene glycol of a higher molecular weight (5000 Da) allows penetration of nanoparticles into the stroma, which proceeds gradually, after an initial 1 h lag phase.

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.

Bdellovibrio bacteriovorus HD100 guards against Pseudomonas tolaasii brown-blotch lesions on the surface of post-harvest Agaricus bisporus supermarket mushrooms

BACKGROUND: Pseudomonas tolaasii is a problematic pathogen of cultured mushrooms, forming dark brown 'blotches' on mushroom surfaces and causing spoilage during crop growth and post-harvest . Treating P. tolaasii infection is difficult, as other, commensal bacterial species such as Pseudomonas putida are necessary for mushroom growth, so treatments must be relatively specific. RESULTS: We have found that P. tolaasii is susceptible to predation in vitro by the δ-proteobacterium Bdellovibrio bacteriovorus. This effect also occurred in funga, where B. bacteriovorus was administered to post-harvest mushroom caps before and after administration of the P. tolaasii pathogen. A significant, visible improvement in blotch appearance, after incubation, was observed on administration of Bdellovibrio. A significant reduction in viable P. tolaasii cell numbers, recovered from the mushroom tissue, was detected. This was accompanied by a more marked reduction in blotch severity on Bdellovibrio administration. We found that there was in some cases an accompanying overgrowth of presumed-commensal, non-Pseudomonas bacteria on post-harvest mushroom caps after Bdellovibrio-treatment. These bacteria were identified (by 16SrRNA gene sequencing) as Enterobacter species, which were seemingly resistant to predation. We visualised predatory interactions occuring between B. bacteriovorus and P. tolaasii on the post-harvest mushroom cap surface by Scanning Electron Microscopy, seeing predatory invasion of P. tolaasii by B. bacteriovorus in funga. This anti-P. tolaasii effect worked well in post-harvest supermarket mushrooms, thus Bdellovibrio was not affected by any pre-treatment of mushrooms for commercial/consumer purposes. CONCLUSIONS: The soil-dwelling B. bacteriovorus HD100 preys upon and kills P. tolaasii, on mushroom surfaces, and could therefore be applied to prevent spoilage in post-harvest situations where mushrooms are stored and packaged for sale.