Large Space Scales Research Articles

Estimating climatological planetary boundary layer heights from radiosonde observations: Comparison of methods and uncertainty analysis

Planetary boundary layer (PBL) processes control energy, water, and pollutant exchanges between the surface and free atmosphere. However, there is no observation‐based global PBL climatology for evaluation of climate, weather, and air quality models or for characterizing PBL variability on large space and time scales. As groundwork for such a climatology, we compute PBL height by seven methods, using temperature, potential temperature, virtual potential temperature, relative humidity, specific humidity, and refractivity profiles from a 10 year, 505‐station radiosonde data set. Six methods are directly compared; they generally yield PBL height estimates that differ by several hundred meters. Relative humidity and potential temperature gradient methods consistently give higher PBL heights, whereas the parcel (or mixing height) method yields significantly lower heights that show larger and more consistent diurnal and seasonal variations (with lower nighttime and wintertime PBLs). Seasonal and diurnal patterns are sometimes associated with local climatological phenomena, such as nighttime radiation inversions, the trade inversion, and tropical convection and associated cloudiness. Surface‐based temperature inversions are a distinct type of PBL that is more common at night and in the morning than during midday and afternoon, in polar regions than in the tropics, and in winter than other seasons. PBL height estimates are sensitive to the vertical resolution of radiosonde data; standard sounding data yield higher PBL heights than high‐resolution data. Several sources of both parametric and structural uncertainty in climatological PBL height values are estimated statistically; each can introduce uncertainties of a few 100 m.

Predicting marine phytoplankton community size structure from empirical relationships with remotely sensed variables

The size composition of primary producers has a potential influence on the length of marine food chains and carbon sinking rates, thus on the proportion of primary production (PP) that is removed from the upper layers and available to higher trophic levels. While total rates of PP are widely reported, it is also necessary to account for the size composition of primary producers when developing food web models that predict consumer biomass and production. Empirical measurement of size composition over large space and time scales is not feasible, so one approach is to predict size composition from environmental variables that are measured and reported on relevant scales. Here, we describe relationships between the environment and the size composition of phytoplankton communities, using a collation of empirical measurements of size composition from sites that include polar, tropical and upwelling environments. The size composition of the phytoplankton communities can be predicted using two remotely sensed variables, chlorophyll-a concentration and sea surface temperature. Applying such relationships in combination allows prediction of the slope and location of phytoplankton size spectra and estimation of the percentage of different sized phytoplankton groups in communities.

Comparison of Different Ground-Based NDVI Measurement Methodologies to Evaluate Crop Biophysical Properties

The usage of vegetation indices such as the Normalized Difference Vegetation Index (NDVI) calculated by means of remote sensing data is widely spread for describing vegetation status on large space scale. However, a big limitation of these indices is their inadequate time resolution for agricultural purposes. This limitation could be overcome by the ground-based vegetation indices that could provide an interesting tool for integrating satellite-based value. In this work, three techniques to calculate the ground-NDVI have been evaluated for sugar beet cultivated in South Italy in all its phenological phases: the NDVI1 based on hand made reflectance measurements, the NDVI2 calculated on automatically reflectance measurements and the broadband NDVIb based on Photosynthetically Active Radiation (PAR) and global radiation measurements. The best performance was obtained by the NDVIb. Moreover, crop-microclimate-NDVI relations were investigated. In particular, the relationship between NDVI and the Leaf Area Index (LAI) was found logarithmic with a saturation of NDVI at LAI around 1.5 m2 m-2. A clear relation was found between NDVI and crop coefficient Kc experimentally determined by the ratio between actual and reference measured or modelled evapotranspirations, while the relation between NDVI and crop actual evapotranspiration was very weak and not usable for practical purposes. Lastly, no relationship between the microclimate and the NDVI was found.

Validity of Uniform Flow Hypothesis in One-Dimensional Morphodynamic Models

Local-uniform-flow (LUF) hypothesis is a simplification of the governing equations describing river morphodynamics, which is needed to determine the evolution of the bed profile and bed-material composition in the case of large time and space scales. This paper presents a rigorous analysis of the full one-dimensional river hydrodynamic and morphodynamic mathematical model compared to its LUF approximation. The analysis establishes two criteria to assess the validity of the LUF hypothesis: (1) a criterion for rivers in equilibrium and (2) a criterion for evolving rivers (i.e., in nonequilibrium). The first criterion is based on the concept of the morphological box. Variations of the river bed longer than the box length are adequately reproduced by the LUF hypothesis, whereas only spatially averaged values are resolved within the box. The second criterion is based on the concept of an evolution window. Temporal variations represented by wave periods larger than the evolution window can be adequately reproduced by the LUF hypothesis, whereas variations with shorter periods are averaged within this window. The minimum size of morphological box and evolution window that limit the error introduced by the LUF hypothesis increases when the Froude number decreases. Further, the minimum size of the evolution window increases for decreasing sediment concentration and increasing mixing layer thickness (i.e., for larger bed forms). The LUF hypothesis is therefore best applied to small mountain rivers for which both the minimum size of the morphological box and the evolution window is relatively small, so that spatial and temporal variations can be resolved in more detail. Applications using the LUF hypothesis for large watersheds (including the lowland portion of the fluvial network) are possible, but are limited to simulations over larger spatial and temporal intervals.

Biogeography of soil microbial communities: a review and a description of the ongoing french national initiative

Microbial biogeography is the study of the distribution of microbial diversity on large scales of space and time. This science aims at understanding biodiversity regulation and its link with ecosystem biological functioning, goods and services such as maintenance of productivity, of soil and atmospheric quality, and of soil health. Although the initial concept dates from the early 20th century (Beijerinck (1913) De infusies en de ontdekking der backterien, in: Jaarboek van de Knoniklijke Akademie van Wetenschappen, Muller, Amsterdam), only recently have an increasing number of studies have investigated the biogeographical patterns of soil microbial diversity. A such delay is due to the constraints of the microbial models, the need to develop relevant molecular and bioinformatic tools to assess microbial diversity, and the non-availability of an adequate sampling strategy. Consequently, the conclusions from microbial ecology studies have rarely been generally applicable and even the fundamental power-laws differ because the taxa-area relationship and the influence of global and distal parameters on the spatial distribution of microbial communities have not been examined. In this article we define and discuss the scientific, technical and operational limits and outcomes resulting from soil microbial biogeography together with the technical and logistical feasibility. The main results are that microbial communities are not stochastically distributed on a wide scale and that biogeographical patterns are more influenced by local parameters such as soil type and land use than by distal ones, e.g. climate and geomorphology, contrary to plants and animals. We then present the European soil biological survey network, focusing on the French national initiative and the „ECOMIC-RMQS” project. The objective of the ECOMIC-RMQS project is to characterise the density and diversity of bacterial communities in all soils in the RMQS library in order to assess, for the first time, not only microbial biogeography across the whole of France but also the impact of land use on soil biodiversity (Reseau de Mesures de la Qualite des Sols = French Soil Quality Monitoring Network, 2200 soils covering all the French territory with a systematic grid of sampling). The scientific, technical and logistical outputs are examined with a view to the future prospects needed to develop this scientific domain and its applications in sustainable land use.

Resolution dependent relative dispersion statistics in a hierarchy of ocean models

In order to determine the effect of Eulerian spatial resolution on the two particle statistics of synthetic drifter trajectories, we examine a hierarchy of ocean models, starting from 2D turbulence simulations, progressing to idealized simulations of a buoyant coastal jet with ROMS, and finally to realistic HYCOM simulations of the Gulf Stream. In each case, particle dispersion at large time and space scales is found to be controlled by energetic meso-scale features of the flow that are relatively insensitive to the resolution of finer scale motions. In all cases, time-distance graphs given in terms of computed Finite Scale Lyapunov Exponents show an expected increase in the extent of exponential scaling with increasing spatial smoothing of the velocity field. The limiting value of the FSLE at small separation distances is found to scale remarkably well with the resolution of Eulerian velocity gradients as given by the average of positive Okubo–Weiss parameter values.

Discrete Event Front-tracking Simulation of a Physical Fire-spread Model

Simulation of moving interfaces such as a fire front usually requires resolution of a large-scale and detailed domain. Such computing involves the use of supercomputers to process the large amount of data and calculations. This limitation is mainly due to the fact that a large scale of space and time is usually split into nodes, cells, or matrices and the solving methods often require small time steps. In this paper we present a novel method that enables the simulation of large-scale/highresolution systems by focusing on the interface and its application to fire-spread simulation. Unlike the conventional explicit and implicit integration schemes, it is based on the discrete-event approach, which describes time advance in terms of increments of physical quantities rather than discrete time stepping. In addition, space is not split into discrete nodes or cells, but we use polygons with real coordinates. The system is described by the behavior of its interface and evolves by computing collision events of this interface in the simulation. As this simulation technique is suitable for a class of models that can explicitly provide the rate of spread, we developed a radiation-based propagation model of wild land fire. Simulations of a real large-scale fire performed by implementation of our method provide very interesting results in less than 30 s with a 3-m resolution with current personal computers.

Recent trends in plankton and upwelling intensity off Galicia (NW Spain)

The mean intensity of the NE Atlantic upwelling system at its northern limit (Galicia, NW Spain) decreased during the last 40 years. At the same time, warming of surface waters was detected. Plankton biomass and composition are expected to reflect such changes when integrated over large time and space scales. In this study, biomass, abundance and species composition of phyto- and zooplankton were analysed to search for significant patterns of annual change and relations with upwelling intensity. Regionally integrated, mostly offshore, data were obtained from the Continuous Plankton Recorder (since 1958) whereas coastal data from Vigo and A Coruña came from the Radiales program (since 1987). No significant trends were found in phytoplankton biomass at either regional or local scales. However, there was a significant decrease in diatom abundance at regional scales and also of large species at local scales. Zooplankton abundance (mainly copepods) significantly decreased offshore but increased near the coast. Biomass of zooplankton also increased near the coast, with the fastest rates in the south. Warm-water species, like Temora stylifera, were increasingly abundant at both regional and local scales. Significant correlations between upwelling intensity and plankton suggest that climatic effects were delayed for several years. Our results indicate that the effects of large scale climatic trends on plankton communities are being effectively modulated within the pelagic ecosystem in this upwelling region.

Natural experiments in landscape evolution

AbstractA natural experiment in landscape evolution is a case study of landform development in which only one element varies significantly, and for which the driving forces, initial conditions, and/or boundary conditions are well constrained. Natural experiments provide a means of testing landscape evolution theory on the large space and time scales to which that theory applies. Natural experiments can involve either steady or transient conditions. Cases with steady conditions allow one to test predictions about the relationships among topography, erosion rates, and various attributes related to climate and material properties. Transient cases are valuable for distinguishing between models whose predictions might be similar, and therefore indistinguishable, under steady conditions. Essential ingredients of a natural experiment include minimal variation in all but one factor, good constraints on timing and/or rates, well‐characterized processes, and high quality topographic data. Other useful ingredients include information about intermediate topographic states (such as a former valley profile revealed by strath terraces), and knowledge of the time history of erosion rates. In order to deepen our understanding of the physics and chemistry of long‐term landscape evolution, there is a pressing need to identify natural experiments and develop the necessary databases to take advantage of them. Copyright © 2009 John Wiley & Sons, Ltd.

Data assimilation in a simple marine ecosystem model based on spatial biological parameterizations

Assuming that a set of constant parameters fits for marine ecosystem modeling and parameter estimation studies on large space scales is questionable since ecosystem types spanning long distances are quite different. In this study, SeaWiFS chlorophyll-a data are assimilated into a simple NPZD model by the adjoint method in a climatological physical environment provided by FOAM. To improve the assimilation results, different spatial parameterization schemes are utilized. The results show that the values of the selected sensitive parameters are spatially variable and the application of spatial parameterizations can improve the assimilation results significantly.

Human Dispersals: Mathematical Models and the Archaeological Record

The theoretical literature on human population dispersal processes at the large time and space scale is reviewed, including references to and discussions of relevant empirical data. The basic Fisher-KPP reaction-diffusion system is summarized for the single population situation, and developments relating to the Allee effect, density-dependent dispersal, time delay, advection, spatial and temporal heterogeneity, and anomalous and stratified diffusion are reviewed. Two- and three-population competitive reaction-diffusion systems of Lotka-Volterra type are also reviewed, as are dynamic approaches to carrying capacity that incorporate predator-prey instabilities, ecosystem engineering, and gene-culture coevolution.

Human Dispersals: Mathematical Models and the Archaeological Record

The theoretical literature on human population dispersal processes at the large time and space scale is reviewed, including references to and discussions of relevant empirical data. The basic Fisher-KPP reaction-diffusion system is summarized for the single population situation, and developments relating to the Allee effect, density-dependent dispersal, time delay, advection, spatial and temporal heterogeneity, and anomalous and stratified diffusion are reviewed. Two- and three-population competitive reaction-diffusion systems of Lotka-Volterra type are also reviewed, as are dynamic approaches to carrying capacity that incorporate predator-prey instabilities, ecosystem engineering, and gene-culture coevolution.

Multiscalarity of the Nutrient–Chlorophyll Relationship in Coastal Phytoplankton

The relationship between nitrogen (N) and phytoplankton chlorophyll a (Chl) establishes a basis for understanding eutrophication in coastal marine ecosystems. A substantial literature exists on cross-ecosystem analysis of this relationship, but there is little information on cross-scale patterns. A collection of observational records in Bedford Basin (Canada) was used to construct the N–Chl relationship at four time scales: intra-day, intra-annual, interannual, and interdecadal. Additionally, a dataset of contingent observations from 16 biogeochemical ocean provinces was used to construct the N–Chl relationship at large spatial scale. In Bedford Basin, N statistically predicts Chl at time scales that are short (intra-day, intra-annual) and long (interdecadal) but not intermediate (interannual). There is an apparent stoichiometric regularity in the dependence of Chl on N that crosses time scales. Further, an apparent similitude exists between the local pattern at long time scale and the global pattern at large space scale.

Multiscale modelling of asphaltene disaggregation

Asphaltene aggregation reduces bitumen upgrading efficiency by increasing bitumen viscosity and coke formation. Our approach to model asphaltene aggregation involves geometry optimisation by using the Harris approximation implemented in DMol3, followed by a solvation calculation by using the three-dimensional molecular theory of solvation (a.k.a. 3D-RISM) we have developed recently [A. Kovalenko, Three-dimensional RISM theory for molecular liquids and solid–liquid interfaces, in Understanding Chemical Reactivity: Molecular Theory of Solvation, F. Hirata ed., Vol. 24, Kluwer Academic Publishers, New York, NY, 2003, pp. 169–275]. From the Harris approximation, we obtain the Hirshfeld and Mulliken asphaltene atomic charges. The 3D-RISM theory allows one to model solvation at given temperature, solvent density and solvent composition. The theory predicts solvation structure and thermodynamic characteristics, such as the potential of mean force (PMF). We investigate the effect of the Hirshfeld and Mulliken charge calculation methods on the PMF values for asphaltene disaggregation in quinoline and 1-methylnaphthalene solvents at 298 and 473 K. Our PMF results predict that asphaltene disaggregation is favoured in quinoline at 473 K, whereas in 1-methylnaphthalene the asphaltene aggregate would remain undisturbed. These results are in agreement with the experiment and correlate with the molecular dynamics (MD) simulation results [T. Takanohashi, S. Sato, and R. Tanaka, Structural relaxation behaviors of three different asphaltenes using MD calculations, Petr. Sci. Technol. 22 (2004), pp. 901–904]. The statistical–mechanical 3D-RISM method probes the entire phase space and yields the solvation structure and thermodynamics at a much lower computational cost than MD, and thus gives access to solvation processes that occur on large time and space scales.

Non-linear advection–diffusion equations approximate swarming but not schooling populations

Advection–diffusion equations (ADEs) are concise and tractable mathematical descriptions of population distributions used widely to address spatial problems in applied and theoretical ecology. We assessed the potential of non-linear ADEs to approximate over very large time and space scales the spatial distributions resulting from social behaviors such as swarming and schooling, in which populations are subdivided into many groups of variable size, velocity and directional persistence. We developed a simple numerical scheme to estimate coefficients in non-linear ADEs from individual-based model (IBM) simulations. Alignment responses between neighbors within groups quantitatively and qualitatively affected how populations moved. Asocial and swarming populations, and schooling populations with weak alignment tendencies, were well approximated by non-linear ADEs. For these behaviors, numerical estimates such as ours could enhance realism and efficiency in ecosystem models of social organisms. Schooling populations with strong alignment were poorly approximated, because (in contradiction to assumptions underlying the ADE approach) effective diffusion and advection were not uniquely defined functions of local density. PDE forms other than ADEs are apparently required to approximate strongly aligning populations.

A parallel FE–FV scheme to solve fluid flow in complex geologic media

Field data-based simulations of geologic systems require much computational time because of their mathematical complexity and the often desired large scales in space and time. To conduct accurate simulations in an acceptable time period, methods to reduce runtime are required. A parallelization approach is attractive because fast multi-processor clusters are nowadays readily available. Here we report on our recent efforts to parallelize our multiphysics code CSMP ++ (Complex System Modelling Platform). In particular, we describe a parallel finite element–finite volume method for multi-phase fluid flow in heterogeneous porous media. We take a domain partitioning approach where the finite element mesh is partitioned into sub-domains, assigning each of them to a single processor. For each sub-domain a local finite volume mesh is constructed. We can now solve advection–dispersion type equations taking an operator splitting approach: Pressure diffusion is calculated with an implicit finite element method and advection with an implicit or explicit finite volume scheme. We have tested the accuracy, robustness and computational speedup of our new parallel scheme on a Linux cluster by means of three geologic applications. All tests give excellent computational speedup with increasing number of up to 32 processors. These results broaden the range of possible simulations in terms of spatial and temporal scale and resolution as well as numerical accuracy up to two orders of magnitude.

Energy flow in plates: Analysis of bending waves and exact quadratic formulation for one-dimensional fields

Thickness-averaged energy densities and energy flow are analyzed in plates, accounting for the hysteretic damping of the material and the external exciting load distribution, without any more approximation than those of the kinematic models. Bending waves, which are governed by bi-harmonic equations, are studied. Complex kinetic- and strain-energy densities, structural intensity and, its divergence and curl are written using some quadratic variables. An exact quadratic formulation, which is based on energy densities and another quadratic variable, is obtained for one-dimensional bending fields, on account of the simplifications then occurring. The fundamental solutions for quadratic variables in an unloaded portion of plate are analyzed; only some of their components, those with the larger space scale, are often considered in energy models of the literature. Then the quadratic formulation for one-dimensional bending waves is applied to a semi-infinite plate excited by a concentrated force. This proves the equivalence of the quadratic and the displacement formulations for modeling energy fields in plates, and the limitation of the quadratic formulation that needs more intensive computations than the displacement formulation for numerical applications.

Large Scale Dynamics of the Persistent Turning Walker Model of Fish Behavior

This paper considers a new model of individual displacement, based on fish motion, the so-called Persistent Turning Walker (PTW) model, which involves an Ornstein-Uhlenbeck process on the curvature of the particle trajectory. The goal is to show that its large time and space scale dynamics is of diffusive type, and to provide an analytic expression of the diffusion coefficient. Two methods are investigated. In the first one, we compute the large time asymptotics of the variance of the individual stochastic trajectories. The second method is based on a diffusion approximation of the kinetic formulation of these stochastic trajectories. The kinetic model is a Fokker-Planck type equation posed in an extended phase-space involving the curvature among the kinetic variables. We show that both methods lead to the same value of the diffusion constant. We present some numerical simulations to illustrate the theoretical results.

Numerical simulation of a stratigraphic model

A multi-lithology diffusive stratigraphic model is considered, which simulates at large scales in space and time the infill of sedimentary basins governed by the interaction between tectonics displacements, eustatic variations, sediment supply, and sediment transport laws. The model accounts for the mass conservation of each sediment lithology resulting in a mixed parabolic, hyperbolic system of partial differential equations (PDEs) for the lithology concentrations and the sediment thickness. It also takes into account a limit on the rock alteration velocity modeled as a unilaterality constraint. To obtain a robust, fast, and accurate simulation, fully and semi-implicit finite volume discre tization schemes are derived for which the existence of stable solutions is proved. Then, the set of nonlinear equations is solved using a Newton algorithm adapted to the unilaterality constraint, and preconditioning strategies are defined for the solution of the linear system at each Newton iteration. They are based on an algebraic approximate decoupling of the sediment thickness and the concentration variables as well as on a proper preconditioning of each variable. These algorithms are studied and compared in terms of robustness, scalability, and efficiency on two real basin test cases.

Transfers through argillaceous rocks over large space and time scales: Insights given by water stable isotopes

In order to characterize the large-scale transport properties of the Toarcian/Domerian formation in the Tournemire Underground Research Laboratory, the isotopic composition ( δ 18O and δ 2H) of pore water was determined on samples from a 250 m borehole drilled from the tunnel ground down to the lower aquifer. The isotope data were obtained from core samples by diffusive vapour exchange. The comparison with previous profiles obtained from vacuum distillation shows a clear discrepancy for all the stratigraphic levels. This demonstrates that artefacts, induced by distillation and already highlighted on only one stratigraphic level, occurred all over the shale sequence. Moreover, diffusion parameters were derived from both radial diffusion experiments and petrophysical determinations. To evaluate the large-scale transport properties of the formation, we performed a series of diffusive model calculations and compared them to experimental data. In agreement with the hydrogeological history, the best simulations considered an activation of the embedding aquifers 11–13 Ma ago, corresponding to the beginning of the diffusive process through the argillaceous formations.