Impacts of changing climate and land use on vegetation dynamics in a Mediterranean ecosystem: insights from paleoecology and dynamic modeling

Forests near the Mediterranean coast have been shaped by millennia of human disturbance. Consequently, ecological studies relying on modern observations or historical records may have difficulty assessing natural vegetation dynamics under current and future climate. We combined a sedimentary pollen record from Lago di Massacciucoli, Tuscany, Italy with simulations from the LandClim dynamic vegetation model to determine what vegetation preceded intense human disturbance, how past changes in vegetation relate to fire and browsing, and the potential of an extinct vegetation type under present climate. We simulated vegetation dynamics near Lago di Massaciucoli for the last 7,000 years using a local chironomid-inferred temperature reconstruction with combinations of three fire regimes (small infrequent, large infrequent, small frequent) and three browsing intensities (no browsing, light browsing, and moderate browsing), and compared model output to pollen data. Simulations with low disturbance support pollen-inferred evidence for a mixed forest dominated by Quercus ilex (a Mediterranean species) and Abies alba (a montane species). Whereas pollen data record the collapse of A. alba after 6000 cal yr bp, simulated populations expanded with declining summer temperatures during the late Holocene. Simulations with increased fire and browsing are consistent with evidence for expansion by deciduous species after A. alba collapsed. According to our combined paleo-environmental and modeling evidence, mixed Q. ilex and A. alba forests remain possible with current climate and limited disturbance, and provide a viable management objective for ecosystems near the Mediterranean coast and in regions that are expected to experience a mediterranean-type climate in the future.

MODELING THE DYNAMICS OF AIRWAY CONSTRICTION: EFFECTS OF AGONIST TRANSPORT AND BINDING

Recent advances have revealed that during exogenous airway challenge, airway diameters can not be adequately predicted by their initial diameters. Furthermore, airway diameters can also vary greatly in time on scales shorter than a breath. In order to better understand these phenomena, we developed a multiscale model which allows us to simulate aerosol challenge in the airways during ventilation. The model incorporates agonist-receptor binding kinetics to govern the temporal response of airway smooth muscle (ASM) contraction on individual airway segments, which together with airway wall mechanics, determines local airway caliber. Global agonist transport and deposition is coupled with pressure-driven flow, linking local airway constrictions with global flow dynamics. During the course of challenge, airway constriction alters the flow pattern, redistributing agonist to less constricted regions. This results in a negative feedback which may be a protective property of the normal lung. As a consequence, repetitive challenge can cause spatial constriction patterns to evolve in time, resulting in a loss of predictability of airway diameters. Additionally, the model offers new insight into several phenomena including the intra- and inter-breath dynamics of airway constriction throughout the tree structure.

Spatial and temporal variation in population genetic structure of wild Nile tilapia (Oreochromis niloticus) across Africa

Background: Reconstructing the evolutionary history of a species is challenging. It often depends not only on the past biogeographic and climatic events but also the contemporary and ecological factors, such as current connectivity and habitat heterogeneity. In fact, these factors might interact with each other and shape the current species distribution. However, to what extent the current population genetic structure reflects the past and the contemporary factors is largely unknown. Here we investigated spatio-temporal genetic structures of Nile tilapia (Oreochromis niloticus) populations, across their natural distribution in Africa. While its large biogeographic distribution can cause genetic differentiation at the paleo-biogeographic scales, its restricted dispersal capacity might induce a strong genetic structure at micro-geographic scales. Results: Using nine microsatellite loci and 350 samples from ten natural populations, we found the highest genetic differentiation among the three ichthyofaunal provinces and regions (Ethiopian, Nilotic and Sudano-Sahelian) (R(ST) = 0.38 - 0.69). This result suggests the predominant effect of paleo-geographic events at macro-geographic scale. In addition, intermediate divergences were found between rivers and lakes within the regions, presumably reflecting relatively recent interruptions of gene flow between hydrographic basins (R(ST) = 0.24 - 0.32). The lowest differentiations were observed among connected populations within a basin (R(ST) = 0.015 in the Volta basin). Comparison of temporal sample series revealed subtle changes in the gene pools in a few generations (F = 0 - 0.053). The estimated effective population sizes were 23 - 143 and the estimated migration rate was moderate (m similar to 0.094 - 0.097) in the Volta populations. Conclusions: This study revealed clear hierarchical patterns of the population genetic structuring of O. niloticus in Africa. The effects of paleo-geographic and climatic events were predominant at macro-geographic scale, and the significant effect of geographic connectivity was detected at micro-geographic scale. The estimated effective population size, the moderate level of dispersal and the rapid temporal change in genetic composition might reflect a potential effect of life history strategy on population dynamics. This hypothesis deserves further investigation. The dynamic pattern revealed at micro-geographic and temporal scales appears important from a genetic resource management as well as from a biodiversity conservation point of view.

Exposure modeling of high-frequency electromagnetic fields

We developed a geospatial model that calculates ambient high-frequency electromagnetic field (HF-EMF) strengths of stationary transmission installations such as mobile phone base stations and broadcast transmitters with high spatial resolution in the order of 1 m. The model considers the location and transmission patterns of the transmitters, the three-dimensional topography, and shielding effects by buildings. The aim of the present study was to assess the suitability of the model for exposure monitoring and for epidemiological research. We modeled time-averaged HF-EMF strengths for an urban area in the city of Basel as well as for a rural area (Bubendorf). To compare modeling with measurements, we selected 20 outdoor measurement sites in Basel and 18 sites in Bubendorf. We calculated Pearson's correlation coefficients between modeling and measurements. Chance-corrected agreement was evaluated by weighted Cohen's kappa statistics for three exposure categories. Correlation between measurements and modeling of the total HF-EMF strength was 0.67 (95% confidence interval (CI): 0.33-0.86) in the city of Basel and 0.77 (95% CI: 0.46-0.91) in the rural area. In both regions, kappa coefficients between measurements and modeling were 0.63 and 0.77 for the total HF-EMF strengths and for all mobile phone frequency bands. First evaluation of our geospatial model yielded substantial agreement between modeling and measurements. However, before the model can be applied for future epidemiologic research, additional validation studies focusing on indoor values are needed to improve model validity.Journal of Exposure Science and Environmental Epidemiology (2008) 18, 183-191; doi:10.1038/sj.jes.7500575; published online 4 April 2007.

Modeling of electroosmotic and electrophoretic mobilization in capillary and microchip isoelectric focusing

Our dynamic capillary electrophoresis model which uses material specific input data for estimation of electroosmosis was applied to investigate fundamental aspects of isoelectric focusing (IEF) in capillaries or microchannels made from bare fused-silica (FS), FS coated with a sulfonated polymer, polymethylmethacrylate (PMMA) and poly(dimethylsiloxane) (PDMS). Input data were generated via determination of the electroosmotic flow (EOF) using buffers with varying pH and ionic strength. Two models are distinguished, one that neglects changes of ionic strength and one that includes the dependence between electroosmotic mobility and ionic strength. For each configuration, the models provide insight into the magnitude and dynamics of electroosmosis. The contribution of each electrophoretic zone to the net EOF is thereby visualized and the amount of EOF required for the detection of the zone structures at a particular location along the capillary, including at its end for MS detection, is predicted. For bare FS, PDMS and PMMA, simulations reveal that EOF is decreasing with time and that the entire IEF process is characterized by the asymptotic formation of a stationary steady-state zone configuration in which electrophoretic transport and electroosmotic zone displacement are opposite and of equal magnitude. The location of immobilization of the boundary between anolyte and most acidic carrier ampholyte is dependent on EOF, i.e. capillary material and anolyte. Overall time intervals for reaching this state in microchannels produced by PDMS and PMMA are predicted to be similar and about twice as long compared to uncoated FS. Additional mobilization for the detection of the entire pH gradient at the capillary end is required. Using concomitant electrophoretic mobilization with an acid as coanion in the catholyte is shown to provide sufficient additional cathodic transport for that purpose. FS capillaries dynamically double coated with polybrene and poly(vinylsulfonate) are predicted to provide sufficient electroosmotic pumping for detection of the entire IEF gradient at the cathodic column end.

Prediction of spatial orientation and morphology of calcaneonavicular coalitions

BACKGROUND: Calcaneonavicular coalitions (CNC) have been reported to be associated with anatomical aberrations of either the calcaneus and/or navicular bones. These morphological abnormalities may complicate accurate surgical resection. Three-dimensional analysis of spatial orientation and morphological characteristics may help in preoperative planning of resection. MATERIALS AND METHODS: Sixteen feet with a diagnosis of CNC were evaluated by means of 3-D CT modeling. Three angles were defined that were expressed in relation to one reproducible landmark (lateral border of the calcaneus): the dorsoplantar inclination, anteroposterior inclination, and socket angle. The depth and width of the coalitions were measured and calculated to obtain the estimated contact surface. Three-dimensional reconstructions of the calcanei served to evaluate the presence, distortion or absence of the anterior calcaneal facet and presence of a navicular beak. The interrater correlations were assessed in order to obtain values for the accuracy of the measurement methods. Sixteen normal feet were used as controls for comparison of the socket angle; anatomy of the anterior calcaneal facet and navicular beak as well. RESULTS: The dorsoplantar inclination angle averaged 50 degrees (+/-17), the anteroposterior inclination angle 64 degrees (+/-15), and the pathologic socket angle 98 degrees (+/-11). The average contact area was 156 mm(2). Ninety-four percent of all patients in the CNC group revealed a plantar navicular beak. In 50% of those patients the anterior calcaneal facet was replaced by the navicular portion and in 44% the facet was totally missing. In contrast, the socket angle in the control group averaged 77 degrees (+/-18), which was found to be statistically different than the CNC group (p = 0.0004). Only 25% of the patients in the control group had a plantar navicular beak. High, statistically significant interrater correlations were found for all measured angles. CONCLUSION: Computer-aided CT analysis and reconstructions help to determine the spatial orientations of CNC in space and provide useful information in order to anticipate morphological abnormalities of the calcaneus and navicular.

Modeling of human P450 oxidoreductase structure by in silico mutagenesis and MD simulation

P450 oxidoreductase (POR) is the obligate electron donor for microsomal cytochrome P450s and mutations in POR cause several metabolic disorders. We have modeled the structure of human P450 oxidoreductase by in silico amino acid replacements in the rat POR crystal structure. The rat POR has 94% homology with human POR and 38 amino acids were replaced to make its sequence identical to human POR. Several rounds of molecular dynamic simulations refined the model and removed structural clashes from side chain alterations of replaced amino acids. This approach has the advantage of keeping the cofactor contacts and structural features of the core enzyme intact which could not be achieved by homology based approaches. The final model from our approach was of high quality and compared well with experimentally determined structures of other PORs. This model will be used for analyzing the structural implications of mutations and polymorphisms in human POR.

Modeling Framework for the Establishment of the Apical-Basal Embryonic Axis in Plants

The apical-basal axis of the early plant embryo determines the body plan of the adult organism. To establish a polarized embryonic axis, plants evolved a unique mechanism that involves directional, cell-to-cell transport of the growth regulator auxin. Auxin transport relies on PIN auxin transporters 1], whose polar subcellular localization determines the flow directionality. PIN-mediated auxin transport mediates the spatial and temporal activity of the auxin response machinery 2-7] that contributes to embryo patterning processes, including establishment of the apical (shoot) and basal (root) embryo poles 8]. However, little is known of upstream mechanisms guiding the (re)polarization of auxin fluxes during embryogenesis 9]. Here, we developed a model of plant embryogenesis that correctly generates emergent cell polarities and auxin-mediated sequential initiation of apical-basal axis of plant embryo. The model relies on two precisely localized auxin sources and a feedback between auxin and the polar, subcellular PIN transporter localization. Simulations reproduced PIN polarity and auxin distribution, as well as previously unknown polarization events during early embryogenesis. The spectrum of validated model predictions suggests that our model corresponds to a minimal mechanistic framework for initiation and orientation of the apical-basal axis to guide both embryonic and postembryonic plant development.

The past ecology of Abies alba provides new perspectives on future responses of silver fir forests to global warming

Paleoecology can provide valuable insights into the ecology of species that complement observation and experiment-based assessments of climate impact dynamics. New paleoecological records (e.g., pollen, macrofossils) from the Italian Peninsula suggest a much wider climatic niche of the important European tree species Abies alba (silver fir) than observed in its present spatial range. To explore this discrepancy between current and past distribution of the species, we analyzed climatic data (temperature, precipitation, frost, humidity, sunshine) and vegetation-independent paleoclimatic reconstructions (e.g., lake levels, chironomids) and use global coupled carbon-cycle climate (NCAR CSM1.4) and dynamic vegetation (LandClim) modeling. The combined evidence suggests that during the mid-Holocene (6000 years ago), prior to humanization of vegetation, A. alba formed forests under conditions that exceeded the modern (1961-1990) upper temperature limit of the species by 5-7°C (July means). Annual precipitation during this natural period was comparable to today (>700-800 mm), with drier summers and wetter winters. In the meso-Mediterranean to sub-Mediterranean forests A. alba co-occurred with thermophilous taxa such as Quercus ilex, Q. pubescens, Olea europaea, Phillyrea, Arbutus, Cistus, Tilia, Ulmus, Acer, Hedera helix, Ilex aquifolium, Taxus, and Vitis. Results from the last interglacial (ca. 130 000-115 000 BP), when human impact was negligible, corroborate the Holocene evidence. Thermophilous Mediterranean A. alba stands became extinct during the last 5000 years when land-use pressure and specifically excessive anthropogenic fire and browsing disturbance increased. Our results imply that the ecology of this key European tree species is not yet well understood. On the basis of the reconstructed realized climatic niche of the species, we anticipate that the future geographic range of A. alba may not contract regardless of migration success, even if climate should become significantly warmer than today with summer temperatures increasing by up to 5-7°C, as long as precipitation does not fall below 700-800 mm/yr, and anthropogenic disturbance (e.g., fire, browsing) does not become excessive. Our finding contradicts recent studies that projected range contractions under global-warming scenarios, but did not factor how millennia of human impacts reduced the realized climatic niche of A. alba.

Spatial gradients of temperature, accumulation and delta18O-ice in Greenland over a series of Dansgaard-Oeschger events

Air and water stable isotope measurements from four Greenland deep ice cores (GRIP, GISP2, NGRIP and NEEM) are investigated over a series of Dansgaard–Oeschger events (DO 8, 9 and 10), which are representative of glacial millennial scale variability. Combined with firn modeling, air isotope data allow us to quantify abrupt temperature increases for each drill site (1σ = 0.6 °C for NEEM, GRIP and GISP2, 1.5 °C for NGRIP). Our data show that the magnitude of stadial–interstadial temperature increase is up to 2 °C larger in central and North Greenland than in northwest Greenland: i.e., for DO 8, a magnitude of +8.8 °C is inferred, which is significantly smaller than the +11.1 °C inferred at GISP2. The same spatial pattern is seen for accumulation increases. This pattern is coherent with climate simulations in response to reduced sea-ice extent in the Nordic seas. The temporal water isotope (δ18O)–temperature relationship varies between 0.3 and 0.6 (±0.08) ‰ °C−1 and is systematically larger at NEEM, possibly due to limited changes in precipitation seasonality compared to GISP2, GRIP or NGRIP. The gas age−ice age difference of warming events represented in water and air isotopes can only be modeled when assuming a 26% (NGRIP) to 40% (GRIP) lower accumulation than that derived from a Dansgaard–Johnsen ice flow model.

Volcanic forcing for climate modeling: a new microphysics-based data set covering years 1600–present

As the understanding and representation of the impacts of volcanic eruptions on climate have improved in the last decades, uncertainties in the stratospheric aerosol forcing from large eruptions are now linked not only to visible optical depth estimates on a global scale but also to details on the size, latitude and altitude distributions of the stratospheric aerosols. Based on our understanding of these uncertainties, we propose a new model-based approach to generating a volcanic forcing for general circulation model (GCM) and chemistry–climate model (CCM) simulations. This new volcanic forcing, covering the 1600–present period, uses an aerosol microphysical model to provide a realistic, physically consistent treatment of the stratospheric sulfate aerosols. Twenty-six eruptions were modeled individually using the latest available ice cores aerosol mass estimates and historical data on the latitude and date of eruptions. The evolution of aerosol spatial and size distribution after the sulfur dioxide discharge are hence characterized for each volcanic eruption. Large variations are seen in hemispheric partitioning and size distributions in relation to location/date of eruptions and injected SO2 masses. Results for recent eruptions show reasonable agreement with observations. By providing these new estimates of spatial distributions of shortwave and long-wave radiative perturbations, this volcanic forcing may help to better constrain the climate model responses to volcanic eruptions in the 1600–present period. The final data set consists of 3-D values (with constant longitude) of spectrally resolved extinction coefficients, single scattering albedos and asymmetry factors calculated for different wavelength bands upon request. Surface area densities for heterogeneous chemistry are also provided.

Modeling of ore-forming and geoenvironmental systems: Roles of fluid flow and chemical reaction processes

Ore-forming and geoenviromental systems commonly involve coupled fluid flowand chemical reaction processes. The advanced numerical methods and computational modeling have become indispensable tools for simulating such processes in recent years. This enables many hitherto unsolvable geoscience problems to be addressed using numerical methods and computational modeling approaches. For example, computational modeling has been successfully used to solve ore-forming and mine site contamination/remediation problems, in which fluid flow and geochemical processes play important roles in the controlling dynamic mechanisms. The main purpose of this paper is to present a generalized overview of: (1) the various classes and models associated with fluid flow/chemically reacting systems in order to highlight possible opportunities and developments for the future; (2) some more general issues that need attention in the development of computational models and codes for simulating ore-forming and geoenviromental systems; (3) the related progresses achieved on the geochemical modeling over the past 50 years or so; (4) the general methodology for modeling of oreforming and geoenvironmental systems; and (5) the future development directions associated with modeling of ore-forming and geoenviromental systems.

Coherence and phase locking in the scalp EEG and between LORETA model sources, and microstates as putative mechanisms of brain temporo-spatial functional organization

Brain electric mechanisms of temporary, functional binding between brain regions are studied using computation of scalp EEG coherence and phase locking, sensitive to time differences of few milliseconds. However, such results if computed from scalp data are ambiguous since electric sources are spatially oriented. Non-ambiguous results can be obtained using calculated time series of strength of intracerebral model sources. This is illustrated applying LORETA modeling to EEG during resting and meditation. During meditation, time series of LORETA model sources revealed a tendency to decreased left-right intracerebral coherence in the delta band, and to increased anterior-posterior intracerebral coherence in the theta band. An alternate conceptualization of functional binding is based on the observation that brain electric activity is discontinuous, i.e., that it occurs in chunks of up to about 100 ms duration that are detectable as quasi-stable scalp field configurations of brain electric activity, called microstates. Their functional significance is illustrated in spontaneous and event-related paradigms, where microstates associated with imagery- versus abstract-type mentation, or while reading positive versus negative emotion words showed clearly different regions of cortical activation in LORETA tomography. These data support the concept that complete brain functions of higher order such as a momentary thought might be incorporated in temporal chunks of processing in the range of tens to about 100 ms as quasi-stable brain states; during these time windows, subprocesses would be accepted as members of the ongoing chunk of processing.