Steady state and dynamic modelling of nitrogen in the River Kennet: impacts of land use change since the 1930s

Steady state and dynamic models have been developed and applied to the River Kennet system. Annual nitrogen exports from the land surface to the river have been estimated based on land use from the 1930s and the 1990s. Long term modelled trends indicate that there has been a large increase in nitrogen transport into the river system driven by increased fertiliser application associated with increased cereal production, increased population and increased livestock levels. The dynamic model INCA Integrated Nitrogen in Catchments. has been applied to simulate the day-to-day transport of N from the terrestrial ecosystem to the riverine environment. This process-based model generates spatial and temporal data and reproduces the observed instream concentrations. Applying the model to current land use and 1930s land use indicates that there has been a major shift in the short term dynamics since the 1930s, with increased river and groundwater concentrations caused by both non-point source pollution from agriculture and point source discharges. �

A dynamic model of neurovascular coupling: Implications for blood vessel dilation and constriction

Neurovascular coupling in response to stimulation of the rat barrel cortex was investigated using concurrent multichannel electrophysiology and laser Doppler flowmetry. The data were used to build a linear dynamic model relating neural activity to blood flow. Local field potential time series were subject to current source density analysis, and the time series of a layer IV sink of the barrel cortex was used as the input to the model. The model output was the time series of the changes in regional cerebral blood flow (CBF). We show that this model can provide excellent fit of the CBF responses for stimulus durations of up to 16 s. The structure of the model consisted of two coupled components representing vascular dilation and constriction. The complex temporal characteristics of the CBF time series were reproduced by the relatively simple balance of these two components. We show that the impulse response obtained under the 16-s duration stimulation condition generalised to provide a good prediction to the data from the shorter duration stimulation conditions. Furthermore, by optimising three out of the total of nine model parameters, the variability in the data can be well accounted for over a wide range of stimulus conditions. By establishing linearity, classic system analysis methods can be used to generate and explore a range of equivalent model structures (e.g., feed-forward or feedback) to guide the experimental investigation of the control of vascular dilation and constriction following stimulation.

A dynamic causal model of the coupling between pulse stimulation and neural activity

We present a dynamic causal model that can explain context-dependent changes in neural responses, in the rat barrel cortex, to an electrical whisker stimulation at different frequencies. Neural responses were measured in terms of local field potentials. These were converted into current source density (CSD) data, and the time series of the CSD sink was extracted to provide a time series response train. The model structure consists of three layers (approximating the responses from the brain stem to the thalamus and then the barrel cortex), and the latter two layers contain nonlinearly coupled modules of linear second-order dynamic systems. The interaction of these modules forms a nonlinear regulatory system that determines the temporal structure of the neural response amplitude for the thalamic and cortical layers. The model is based on the measured population dynamics of neurons rather than the dynamics of a single neuron and was evaluated against CSD data from experiments with varying stimulation frequency (1–40 Hz), random pulse trains, and awake and anesthetized animals. The model parameters obtained by optimization for different physiological conditions (anesthetized or awake) were significantly different. Following Friston, Mechelli, Turner, and Price (2000), this work is part of a formal mathematical system currently being developed (Zheng et al., 2005) that links stimulation to the blood oxygen level dependent (BOLD) functional magnetic resonance imaging (fMRI) signal through neural activity and hemodynamic variables. The importance of the model described here is that it can be used to invert the hemodynamic measurements of changes in blood flow to estimate the underlying neural activity.

A model of the dynamic relationship between blood flow and volume changes during brain activation

The temporal relationship between changes in cerebral blood flow (CBF) and cerebral blood volume (CBV) is important in the biophysical modeling and interpretation of the hemodynamic response to activation, particularly in the context of magnetic resonance imaging and the blood oxygen level-dependent signal. Grubb et al. (1974) measured the steady state relationship between changes in CBV and CBF after hypercapnic challenge. The relationship CBV proportional to CBFPhi has been used extensively in the literature. Two similar models, the Balloon (Buxton et al., 1998) and the Windkessel (Mandeville et al., 1999), have been proposed to describe the temporal dynamics of changes in CBV with respect to changes in CBF. In this study, a dynamic model extending the Windkessel model by incorporating delayed compliance is presented. The extended model is better able to capture the dynamics of CBV changes after changes in CBF, particularly in the return-to-baseline stages of the response.

Dynamic hydrological modeling in drylands with TRMM based rainfall

This paper introduces and evaluates DryMOD, a dynamic water balance model of the key hydrological process in drylands that is based on free, public-domain datasets. The rainfall model of DryMOD makes optimal use of spatially disaggregated Tropical Rainfall Measuring Mission (TRMM) datasets to simulate hourly rainfall intensities at a spatial resolution of 1-km. Regional-scale applications of the model in seasonal catchments in Tunisia and Senegal characterize runoff and soil moisture distribution and dynamics in response to varying rainfall data inputs and soil properties. The results highlight the need for hourly-based rainfall simulation and for correcting TRMM 3B42 rainfall intensities for the fractional cover of rainfall (FCR). Without FCR correction and disaggregation to 1 km, TRMM 3B42 based rainfall intensities are too low to generate surface runoff and to induce substantial changes to soil moisture storage. The outcomes from the sensitivity analysis show that topsoil porosity is the most important soil property for simulation of runoff and soil moisture. Thus, we demonstrate the benefit of hydrological investigations at a scale, for which reliable information on soil profile characteristics exists and which is sufficiently fine to account for the heterogeneities of these. Where such information is available, application of DryMOD can assist in the spatial and temporal planning of water harvesting according to runoff-generating areas and the runoff ratio, as well as in the optimization of agricultural activities based on realistic representation of soil moisture conditions.

Calibration of structure in a distributed forecasting model for a semiarid flash flood: dynamic surface storage and channel roughness

Flash floods pose a significant danger for life and property. Unfortunately, in arid and semiarid environment the runoff generation shows a complex non-linear behavior with a strong spatial and temporal non-uniformity. As a result, the predictions made by physically-based simulations in semiarid areas are subject to great uncertainty, and a failure in the predictive behavior of existing models is common. Thus better descriptions of physical processes at the watershed scale need to be incorporated into the hydrological model structures. For example, terrain relief has been systematically considered static in flood modelling at the watershed scale. Here, we show that the integrated effect of small distributed relief variations originated through concurrent hydrological processes within a storm event was significant on the watershed scale hydrograph. We model these observations by introducing dynamic formulations of two relief-related parameters at diverse scales: maximum depression storage, and roughness coefficient in channels. In the final (a posteriori) model structure these parameters are allowed to be both time-constant or time-varying. The case under study is a convective storm in a semiarid Mediterranean watershed with ephemeral channels and high agricultural pressures (the Rambla del Albujón watershed; 556 km 2 ), which showed a complex multi-peak response. First, to obtain quasi-sensible simulations in the (a priori) model with time-constant relief-related parameters, a spatially distributed parameterization was strictly required. Second, a generalized likelihood uncertainty estimation (GLUE) inference applied to the improved model structure, and conditioned to observed nested hydrographs, showed that accounting for dynamic relief-related parameters led to improved simulations. The discussion is finally broadened by considering the use of the calibrated model both to analyze the sensitivity of the watershed to storm motion and to attempt the flood forecasting of a stratiform event with highly different behavior.

Dynamic auroral structure in the vicinity of the polar cusp - multipoint observations during southward and northward IMF

The low- and high-latitude boundary layers of the earth's magnetosphere [low-latitude boundary layer (LLBL) and mantle] play important roles in transferring momentum and energy from the solar wind to the magnetosphere-ionosphere system. Particle precipitation, field-aligned current, auroral emission, ionospheric ion drift and ground magnetic perturbations are among the low-altitude parameters that show signatures of various plasma processes in the LLBL and the magnetopause current layer. Magnetic merging events, Kelvin-Helmholtz waves, and pressure pulses excited by the variable solar wind/magnetosheath plasma are examples of boundary phenomena that may be coupled to the ionosphere via field-aligned currents. Optical auroral observation, by photometry and all-sky TV cameras, is a unique technique for investigating the spatial and temporal structure of the electron precipitation associated with such phenomena. However, the distinction between the different boundary layer plasma populations cannot in general be unambiguously determined by optics alone. Additional information, such as satellite observations of particle boundaries and field-aligned currents, is needed in order to identify the plasma source(s) and the magnetosphere-ionosphere coupling mode(s). Two categories of auroral activity/structure in the vicinity of the polar cusp are discussed in this paper, based on combined ground and satellite data. In one case, the quasi-periodic sequence of auroral events at the polar cap boundary involves accelerated electrons (< 1 keV) moving poleward (< 1 km s-1) and azimuthally along the persistent cusp/cleft arc poleward boundary with velocities (< 4 km s-1), comparable to the local ionospheric ion drift during periods of southward IMF. A critical question is whether or not the optical events signify a corresponding plasma flow across the open/closed field line boundary in such cases. Near-simultaneous observations of magnetopause flux transfer events (FTEs) and such optical/ion drift events are reported. The reverse pattern of motion of discrete auroral forms is observed during positive interplanetary magnetic field (IMF) B(Z), i.e. equatorward motion into the cusp/cleft background arc from the poleward edge. Combined satellite and ground-based information for the latter cases indicate a source mechanism, poleward of the cusp at the high-latitude magnetopause or plasma mantle, giving rise to strong momentum transfer and electron precipitation structures within a approximately 200 km-wide latitudinal zone at the cusp/cleft poleward boundary. The striking similarities of auroral electrodynamics in the cleft/mantle region during northward and southward IMF indicate that a qualitatively similar solar wind-magnetosphere coupling mode is operating. It is suggested that, in both cases, the discrete auroral forms represent temporal/spatial structure of larger-scale convection over the polar magnetosphere.

Memory and burstiness in dynamic networks

A discrete-time random process is described, which can generate bursty sequences of events. A Bernoulli process, where the probability of an event occurring at time t is given by a fixed probability x, is modified to include a memory effect where the event probability is increased proportionally to the number of events that occurred within a given amount of time preceding t. For small values of x the interevent time distribution follows a power law with exponent −2−x. We consider a dynamic network where each node forms, and breaks connections according to this process. The value of x for each node depends on the fitness distribution, \rho(x), from which it is drawn; we find exact solutions for the expectation of the degree distribution for a variety of possible fitness distributions, and for both cases where the memory effect either is, or is not present. This work can potentially lead to methods to uncover hidden fitness distributions from fast changing, temporal network data, such as online social communications and fMRI scans.

Time-lag in biological responses to landscape changes in a highly dynamic Atlantic forest region

Time-lagged responses of biological variables to landscape modifications are widely recognized, but rarely considered in ecological studies. In order to test for the existence of time-lags in the response of trees, small mammals, birds and frogs to changes in fragment area and connectivity, we studied a fragmented and highly dynamic landscape in the Atlantic forest region. We also investigated the biological correlates associated with differential responses among taxonomic groups. Species richness and abundance for four taxonomic groups were measured in 21 secondary forest fragments during the same period (2000-2002), following a standardized protocol. Data analyses were based on power regressions and model selection procedures. The model inputs included present (2000) and past (1962, 1981) fragment areas and connectivity, as well as observed changes in these parameters. Although past landscape structure was particularly relevant for trees, all taxonomic groups (except small mammals) were affected by landscape dynamics, exhibiting a time-lagged response. Furthermore, fragment area was more important for species groups with lower dispersal capacity, while species with higher dispersal ability had stronger responses to connectivity measures. Although these secondary forest fragments still maintain a large fraction of their original biodiversity, the delay in biological response combined with high rates of deforestation and fast forest regeneration imply in a reduction in the average age of the forest. This also indicates that future species losses are likely, especially those that are more strictly-forest dwellers. Conservation actions should be implemented to reduce species extinction, to maintain old-growth forests and to favour the regeneration process. Our results demonstrate that landscape history can strongly affect the present distribution pattern of species in fragmented landscapes, and should be considered in conservation planning. (C) 2009 Elsevier Ltd. All rights reserved.

Wavelet-based dynamic time warping

Dynamic Time Warping (DTW), a pattern matching technique traditionally used for restricted vocabulary speech recognition, is based on a temporal alignment of the input signal with the template models. The principal drawback of DTW is its high computational cost as the lengths of the signals increase. This paper shows extended results over our previously published conference paper, which introduces an optimized version of the DTW I hat is based on the Discrete Wavelet Transform (DWT). (C) 2008 Elsevier B.V. All rights reserved.

Dynamic Changes in the Mental Rotation Network Revealed by Pattern Recognition Analysis of fMRI Data

We investigated the temporal dynamics and changes in connectivity in the mental rotation network through the application of spatio-temporal support vector machines (SVMs). The spatio-temporal SVM [Mourao-Miranda, J., Friston, K. J., et al. (2007). Dynamic discrimination analysis: A spatial-temporal SVM. Neuroimage, 36, 88-99] is a pattern recognition approach that is suitable for investigating dynamic changes in the brain network during a complex mental task. It does not require a model describing each component of the task and the precise shape of the BOLD impulse response. By defining a time window including a cognitive event, one can use spatio-temporal fMRI observations from two cognitive states to train the SVM. During the training, the SVM finds the discriminating pattern between the two states and produces a discriminating weight vector encompassing both voxels and time (i.e., spatio-temporal maps). We showed that by applying spatio-temporal SVM to an event-related mental rotation experiment, it is possible to discriminate between different degrees of angular disparity (0 degrees vs. 20 degrees, 0 degrees vs. 60 degrees, and 0 degrees vs. 100 degrees), and the discrimination accuracy is correlated with the difference in angular disparity between the conditions. For the comparison with highest accuracy (08 vs. 1008), we evaluated how the most discriminating areas (visual regions, parietal regions, supplementary, and premotor areas) change their behavior over time. The frontal premotor regions became highly discriminating earlier than the superior parietal cortex. There seems to be a parcellation of the parietal regions with an earlier discrimination of the inferior parietal lobe in the mental rotation in relation to the superior parietal. The SVM also identified a network of regions that had a decrease in BOLD responses during the 100 degrees condition in relation to the 0 degrees condition (posterior cingulate, frontal, and superior temporal gyrus). This network was also highly discriminating between the two conditions. In addition, we investigated changes in functional connectivity between the most discriminating areas identified by the spatio-temporal SVM. We observed an increase in functional connectivity between almost all areas activated during the 100 degrees condition (bilateral inferior and superior parietal lobe, bilateral premotor area, and SMA) but not between the areas that showed a decrease in BOLD response during the 100 degrees condition.

Influência de parâmetros moleculares em funções de correlação temporal na dinâmica de solvatação mecânica

No presente trabalho descrevemos nossos resultados relativos à investigação da dinâmica de solvatação mecânica por meio de simulações por dinâmica molecular, respeitando o regime da resposta linear, em sistemas-modelo de argônio líquido com um soluto monoatômico ou diatômico dissolvido. Estudamos sistematicamente a influência dos parâmetros moleculares dos solutos (tamanho, polarizabilidade) e da densidade frente a vários modelos de solvatação. Funções de Correlação Temporal da Energia de Solvatação foram calculadas com relação à correlações de n-corpos (n = 2; 3) distinguindo interações repulsivas e atrativas para ambos os sistemas líquidos. Também obtivemos segundas derivadas temporais dessas funções referindo-se à parcelas translacionais, rotacionais e roto-translacionais na solução do diatômico. Encontramos que funções de correlação temporal coletivas podem ser razoavelmente bem aproximadas por correlações binárias a densidades baixas e, a densidades altas, correlações ternárias tornam-se mais importantes produzindo um descorrelacionamento mais rápido das funções coletivas devido a efeitos de cancelamento parciais. As funções de correlação para interações repulsivas e atrativas exibem comportamentos dinâmicos independentes do modelo de solvatação devido a fatores de escalonamento linear que afetam apenas as amplitudes das dessas funções de correlação temporal. Em geral, os sistemas com grau de liberdade rotacional apresentam tempos de correlação mais curtos para a dinâmica coletiva e tempos de correlação mais longos para as funções binárias e ternárias. Finalmente, esse estudo mostra que os sistemas contendo o diatômico relaxam-se predominantemente por mecanismos translacionais binários em modelos de solvatação envolvendo alterações apenas na polarizabilidade do soluto, e por mecanismos rotacionais atrativos binários em modelos envolvendo alterações no comprimento de ligação.

Cartografia temporal e análise geoambiental da dinâmica da foz do Rio Piranhas-Açu, regiao de Macau-RN, com base em imagens Landast 5-MT

This work embraces the application of Landsat 5-TM digital images, comprising August 2 1989 and September 22 1998, for temporal mapping and geoenvironmental analysis of the dynamic of Piranhas-Açu river mouth, situated in the Macau (RN) region. After treatment using several digital processing techniques (e.g. colour composition in RGB, ratio of bands, principal component analysis, index methods, among others), it was possible to generate several image products and multitemporal maps of the coastal morphodynamics of the studied area. Using the image products it was possible the identification and characterization of the principal elements of interest (vegetation, soil, geology and water) in the surface of the studied area, associating the spectral characteristics of these elements to that presented by the image products resulting of the digital processing. Thus, it was possible to define different types of soils: Amd, AQd6, SK1 and LVe4; vegetation grouping: open arboreal-shrubby caatinga, closed arborealshrubby caatinga, closed arboreal caatinga, mangrove vegetation, dune vegetation and areas predominately constituted by juremas; geological units: quaternary units beach sediments, sand banks, dune flats, barrier island, mobile dunes, fixed dunes, alluvium, tidal and inundation flats, and sandy facies of the Potengi Formation; tertiary-quaternary units Barreiras Formation grouped to the clayey facies of the Potengi Formation, Macau Formation grouped to the sediments of the Tibau Formation; Cretaceous units Jandaíra Formation; moreover it was to identify the sea/land limit, shallow submersed areas and suspended sediments. The multitemporal maps of the coastal morphodynamics allowed the identification and a semi-quantitative evoluation of regions which were submitted to erosive and constructive processes in the last decade. This semi-quantitative evoluation in association with an geoenvironmental characterization of the studied area are important data to the elaboration of actions that may minimize the possible/probable impacts caused by the implantation of the Polo Gas/Sal and to the monitoring of areas explorated by the petroleum and salt industries

Variação temporal da resistividade elétrica em contaminação por gasolina

The definition of models and base parameters for application of investigation tools in environment with high complexity are a basic premise to any sciences. The geophysics is a science with solid theoretical base and applications in diversified areas of the sciences geological, astronomical, meteorological, among many others. Its application in environmental studies is relatively recent and needs further research. To understand the behavior of resident contaminants in a dynamic and complex environment as the geological, it requests studies in scale laboratorial, under control of factors seasonality variable. This work simulates a leak of gasoline in soil, under conditions and in laboratory scale, with the objective of monitoring the temporary behavior of the hydrocarbon under the optics of variation of the parameter physical electric resistivity. The results indicate increase of the resistivity in recent periods the contamination, followed for stability in the values and finally fall and return tendency to the natural conditions.

Spatial and temporal variability in soil CO2-C emissions and relation to soil temperature at King George Island, maritime Antarctica

Few studies have examined the effects of temperature on spatial and temporal trends in soil CO2-C emissions in Antarctica. In this work, we present in situ measurements of CO2-C emissions and assess their relation with soil temperature, using dynamic chambers. We found an exponential relation between CO2 emissions and soil temperature, with the value of Q10 being close to 2.1. Mean emission rates were as low as 0.026 and 0.072 g of CO2-C m-2 h-1 for bare soil and soil covered with moss, respectively, and as high as 0.162 g of CO2-C m-2 h-1 for soil covered with grass, Deschampsia antarctica Desv. (Poaceae). A spatial variability analysis conducted using a 60-point grid, for an area with mosses (Sannionia uncianata) and D. antarctica, yielded a spherical semivariogram model for CO2-C emissions with a range of 1 m. The results suggest that soil temperature is a controlling factor on temporal variations in soil CO2-C emissions, although spatial variations appear to be more strongly related to the distribution of vegetation types. © 2010 Elsevier B.V. and NIPR.