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.

Electromagnetic atmosphere-plasma coupling: the global atmospheric electric circuit

A description is given of the global atmospheric electric circuit operating between the Earth’s surface and the ionosphere. Attention is drawn to the huge range of horizontal and vertical spatial scales, ranging from 10−9 m to 1012 m, concerned with the many important processes at work. A similarly enormous range of time scales is involved from 10−6 s to 109 s, in the physical effects and different phenomena that need to be considered. The current flowing in the global circuit is generated by disturbed weather such as thunderstorms and electrified rain/shower clouds, mostly occurring over the Earth’s land surface. The profile of electrical conductivity up through the atmosphere, determined mainly by galactic cosmic ray ionization, is a crucial parameter of the circuit. Model simulation results on the variation of the ionospheric potential, ∼250 kV positive with respect to the Earth’s potential, following lightning discharges and sprites are summarized. Experimental results comparing global circuit variations with the neutron rate recorded at Climax, Colorado, are then discussed. Within the return (load) part of the circuit in the fair weather regions remote from the generators, charge layers exist on the upper and lower edges of extensive layer clouds; new experimental evidence for these charge layers is also reviewed. Finally, some directions for future research in the subject are suggested.

Kinetic multi-layer model of gas-particle interactions in aerosols and clouds (KM-GAP): linking condensation, evaporation and chemical reactions of organics, oxidants and water

We present a novel kinetic multi-layer model for gas-particle interactions in aerosols and clouds (KM-GAP) that treats explicitly all steps of mass transport and chemical reaction of semi-volatile species partitioning between gas phase, particle surface and particle bulk. KM-GAP is based on the PRA model framework (Pöschl-Rudich-Ammann, 2007), and it includes gas phase diffusion, reversible adsorption, surface reactions, bulk diffusion and reaction, as well as condensation, evaporation and heat transfer. The size change of atmospheric particles and the temporal evolution and spatial profile of the concentration of individual chemical species can be modelled along with gas uptake and accommodation coefficients. Depending on the complexity of the investigated system, unlimited numbers of semi-volatile species, chemical reactions, and physical processes can be treated, and the model shall help to bridge gaps in the understanding and quantification of multiphase chemistry and microphysics in atmo- spheric aerosols and clouds. In this study we demonstrate how KM-GAP can be used to analyze, interpret and design experimental investigations of changes in particle size and chemical composition in response to condensation, evaporation, and chemical reaction. For the condensational growth of water droplets, our kinetic model results provide a direct link between laboratory observations and molecular dynamic simulations, confirming that the accommodation coefficient of water at 270 K is close to unity. Literature data on the evaporation of dioctyl phthalate as a function of particle size and time can be reproduced, and the model results suggest that changes in the experimental conditions like aerosol particle concentration and chamber geometry may influence the evaporation kinetics and can be optimized for eðcient probing of specific physical effects and parameters. With regard to oxidative aging of organic aerosol particles, we illustrate how the formation and evaporation of volatile reaction products like nonanal can cause a decrease in the size of oleic acid particles exposed to ozone.

Kinetic multi-layer model of gas-particle interactions in aerosols and clouds (KM-GAP): linking condensation, evaporation and chemical reactions of organics, oxidants and water

We present a novel kinetic multi-layer model for gas-particle interactions in aerosols and clouds (KMGAP) that treats explicitly all steps of mass transport and chemical reaction of semi-volatile species partitioning between gas phase, particle surface and particle bulk. KMGAP is based on the PRA model framework (P¨oschl-Rudich- Ammann, 2007), and it includes gas phase diffusion, reversible adsorption, surface reactions, bulk diffusion and reaction, as well as condensation, evaporation and heat transfer. The size change of atmospheric particles and the temporal evolution and spatial profile of the concentration of individual chemical species can be modeled along with gas uptake and accommodation coefficients. Depending on the complexity of the investigated system and the computational constraints, unlimited numbers of semi-volatile species, chemical reactions, and physical processes can be treated, and the model shall help to bridge gaps in the understanding and quantification of multiphase chemistry and microphysics in atmospheric aerosols and clouds. In this study we demonstrate how KM-GAP can be used to analyze, interpret and design experimental investigations of changes in particle size and chemical composition in response to condensation, evaporation, and chemical reaction. For the condensational growth of water droplets, our kinetic model results provide a direct link between laboratory observations and molecular dynamic simulations, confirming that the accommodation coefficient of water at 270K is close to unity (Winkler et al., 2006). Literature data on the evaporation of dioctyl phthalate as a function of particle size and time can be reproduced, and the model results suggest that changes in the experimental conditions like aerosol particle concentration and chamber geometry may influence the evaporation kinetics and can be optimized for efficient probing of specific physical effects and parameters. With regard to oxidative aging of organic aerosol particles, we illustrate how the formation and evaporation of volatile reaction products like nonanal can cause a decrease in the size of oleic acid particles exposed to ozone.

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.

A period 2 genetic variant interacts with plasma SFA to modify plasma lipid concentrations in adults with metabolic syndrome

Genetic variants of Period 2 (PER2), a circadian clock gene, have been linked to metabolic syndrome (MetS). However, it is still unknown whether these genetic variants interact with the various types of plasma fatty acids. This study investigated whether common single nucleotide polymorphisms (SNPs) in the PER2 locus (rs934945 and rs2304672) interact with various classes of plasma fatty acids to modulate plasma lipid metabolism in 381 participants with MetS in the European LIPGENE study. Interestingly, the rs2304672 SNP interacted with plasma total SFA concentrations to affect fasting plasma TG, TG-rich lipoprotein (TRL-TG), total cholesterol, apoC-II, apoB, and apoB-48 concentrations (P-interaction < 0.001–0.046). Carriers of the minor allele (GC+GG) with the highest SFA concentration (>median) had a higher plasma TG concentration (P = 0.001) and higher TRL-TG (P < 0.001) than the CC genotype. In addition, participants carrying the minor G allele for rs2304672 SNP and with a higher SFA concentration (>median) had higher plasma concentrations of apo C-II (P < 0.001), apo C-III (P = 0.009), and apoB-48 (P = 0.028) compared with the homozygotes for the major allele (CC). In summary, the rs2304672 polymorphism in the PER2 gene locus may influence lipid metabolism by interacting with the plasma total SFA concentration in participants with MetS. The understanding of these gene-nutrient interactions could help to provide a better knowledge of the pathogenesis in MetS.

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.

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.

Gene-nutrient interactions on the phosphoenolpyruvate carboxykinase influence insulin sensitivity in metabolic syndrome subjects

Genetic background may interact with habitual dietary fat composition, and affect development of the metabolic syndrome (MetS). The phosphoenolpyruvate carboxykinase gene (PCK1) plays a significant role regulating glucose metabolism, and fatty acids are key metabolic regulators, which interact with transcription factors and influence glucose metabolism. We explored genetic variability at the PCK1 gene locus in relation to degree of insulin resistance and plasma fatty acid levels in MetS subjects. Moreover, we analyzed the PCK1 gene expression in the adipose tissue of a subgroup of MetS subjects according to the PCK1 genetic variants.

The bile acid receptor TGR5 does not interact with β-arrestins or traffic to endosomes but transmits sustained signals from plasma membrane rafts.

TGR5 is a G protein-coupled receptor that mediates bile acid (BA) effects on energy balance, inflammation, digestion and sensation. The mechanisms and spatiotemporal control of TGR5 signaling are poorly understood. We investigated TGR5 signaling and trafficking in transfected HEK293 cells and colonocytes (NCM460) that endogenously express TGR5. BAs (deoxycholic acid, DCA, taurolithocholic acid, TLCA) and the selective agonists oleanolic acid (OA) and 3-(2-chlorophenyl)-N-(4-chlorophenyl)-N, 5-dimethylisoxazole-4-carboxamide (CCDC) stimulated cAMP formation but did not induce TGR5 endocytosis or recruitment of β-arrestins, assessed by confocal microscopy. DCA, TLCA and OA did not stimulate TGR5 association with β-arrestin 1/2 or G protein-coupled receptor kinase (GRK) 2/5/6, determined by bioluminescence resonance energy transfer. CCDC stimulated a low level of TGR5 interaction with β-arrestin2 and GRK2. DCA induced cAMP formation at the plasma membrane and cytosol, determined using exchange factor directly regulated by cAMP (Epac2)-based reporters, but cAMP signals did not desensitize. AG1478, an inhibitor of epidermal growth factor receptor (EGFR) tyrosine kinase, the metalloprotease inhibitor batimastat, and methyl-β-cyclodextrin and filipin, which block lipid raft formation, prevented DCA stimulation of extracellular signal regulated kinase (ERK1/2). BRET analysis revealed TGR5 and EGFR interactions that were blocked by disruption of lipid rafts. DCA stimulated TGR5 redistribution to plasma membrane microdomains, localized by immunogold electron microscopy. Thus, TGR5 does not interact with β-arrestins, desensitize or traffic to endosomes. TGR5 signals from plasma membrane rafts that facilitate EGFR interaction and transactivation. An understanding of the spatiotemporal control of TGR5 signaling provides insights into the actions of BAs and therapeutic TGR5 agonists/antagonists.

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.

PPARGC1A sequence variation and cardiovascular risk-factor levels: a study of the main genetic effects and gene x environment interactions in children from the European Youth Heart Study.

AIMS/HYPOTHESIS: The PPARGC1A gene coactivates multiple nuclear transcription factors involved in cellular energy metabolism and vascular stasis. In the present study, we genotyped 35 tagging polymorphisms to capture all common PPARGC1A nucleotide sequence variations and tested for association with metabolic and cardiovascular traits in 2,101 Danish and Estonian boys and girls from the European Youth Heart Study, a multicentre school-based cross-sectional cohort study. METHODS: Fasting plasma glucose concentrations, anthropometric variables and blood pressure were measured. Habitual physical activity and aerobic fitness were objectively assessed using uniaxial accelerometry and a maximal aerobic exercise stress test on a bicycle ergometer, respectively. RESULTS: In adjusted models, nominally significant associations were observed for BMI (rs10018239, p = 0.039), waist circumference (rs7656250, p = 0.012; rs8192678 [Gly482Ser], p = 0.015; rs3755863, p = 0.02; rs10018239, beta = -0.01 cm per minor allele copy, p = 0.043), systolic blood pressure (rs2970869, p = 0.018) and fasting glucose concentrations (rs11724368, p = 0.045). Stronger associations were observed for aerobic fitness (rs7656250, p = 0.005; rs13117172, p = 0.008) and fasting glucose concentrations (rs7657071, p = 0.002). None remained significant after correcting for the number of statistical comparisons. We proceeded by testing for gene x physical activity interactions for the polymorphisms that showed nominal evidence of association in the main effect models. None of these tests was statistically significant. CONCLUSIONS/INTERPRETATION: Variants at PPARGC1A may influence several metabolic traits in this European paediatric cohort. However, variation at PPARGC1A is unlikely to have a major impact on cardiovascular or metabolic health in these children.