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.

An easily implemented agro-hydrological procedure with dynamic root simulation for water transfer in the crop-soil system: validation and application

Models for water transfer in the crop-soil system are key components of agro-hydrological models for irrigation, fertilizer and pesticide practices. Many of the hydrological models for water transfer in the crop-soil system are either too approximate due to oversimplified algorithms or employ complex numerical schemes. In this paper we developed a simple and sufficiently accurate algorithm which can be easily adopted in agro-hydrological models for the simulation of water dynamics. We used a dual crop coefficient approach proposed by the FAO for estimating potential evaporation and transpiration, and a dynamic model for calculating relative root length distribution on a daily basis. In a small time step of 0.001 d, we implemented algorithms separately for actual evaporation, root water uptake and soil water content redistribution by decoupling these processes. The Richards equation describing soil water movement was solved using an integration strategy over the soil layers instead of complex numerical schemes. This drastically simplified the procedures of modeling soil water and led to much shorter computer codes. The validity of the proposed model was tested against data from field experiments on two contrasting soils cropped with wheat. Good agreement was achieved between measurement and simulation of soil water content in various depths collected at intervals during crop growth. This indicates that the model is satisfactory in simulating water transfer in the crop-soil system, and therefore can reliably be adopted in agro-hydrological models. Finally we demonstrated how the developed model could be used to study the effect of changes in the environment such as lowering the groundwater table caused by the construction of a motorway on crop transpiration. (c) 2009 Elsevier B.V. All rights reserved.

Using coordinated observations in polarized white light and Faraday rotation to probe the spatial position and magnetic field of an interplanetary sheath

Coronal mass ejections (CMEs) can be continuously tracked through a large portion of the inner heliosphere by direct imaging in visible and radio wavebands. White light (WL) signatures of solar wind transients, such as CMEs, result from Thomson scattering of sunlight by free electrons and therefore depend on both viewing geometry and electron density. The Faraday rotation (FR) of radio waves from extragalactic pulsars and quasars, which arises due to the presence of such solar wind features, depends on the line-of-sight magnetic field component B ∥ and the electron density. To understand coordinated WL and FR observations of CMEs, we perform forward magnetohydrodynamic modeling of an Earth-directed shock and synthesize the signatures that would be remotely sensed at a number of widely distributed vantage points in the inner heliosphere. Removal of the background solar wind contribution reveals the shock-associated enhancements in WL and FR. While the efficiency of Thomson scattering depends on scattering angle, WL radiance I decreases with heliocentric distance r roughly according to the expression Ir –3. The sheath region downstream of the Earth-directed shock is well viewed from the L4 and L5 Lagrangian points, demonstrating the benefits of these points in terms of space weather forecasting. The spatial position of the main scattering site r sheath and the mass of plasma at that position M sheath can be inferred from the polarization of the shock-associated enhancement in WL radiance. From the FR measurements, the local B ∥sheath at r sheath can then be estimated. Simultaneous observations in polarized WL and FR can not only be used to detect CMEs, but also to diagnose their plasma and magnetic field properties.

Complex-valued B-spline neural networks for modeling and inverting Hammerstein systems

Many communication signal processing applications involve modelling and inverting complex-valued (CV) Hammerstein systems. We develops a new CV B-spline neural network approach for efficient identification of the CV Hammerstein system and effective inversion of the estimated CV Hammerstein model. Specifically, the CV nonlinear static function in the Hammerstein system is represented using the tensor product from two univariate B-spline neural networks. An efficient alternating least squares estimation method is adopted for identifying the CV linear dynamic model’s coefficients and the CV B-spline neural network’s weights, which yields the closed-form solutions for both the linear dynamic model’s coefficients and the B-spline neural network’s weights, and this estimation process is guaranteed to converge very fast to a unique minimum solution. Furthermore, an accurate inversion of the CV Hammerstein system can readily be obtained using the estimated model. In particular, the inversion of the CV nonlinear static function in the Hammerstein system can be calculated effectively using a Gaussian-Newton algorithm, which naturally incorporates the efficient De Boor algorithm with both the B-spline curve and first order derivative recursions. The effectiveness of our approach is demonstrated using the application to equalisation of Hammerstein channels.

Modeling migratory grazing of zooplankton on toxic and non-toxic phytoplankton

Migratory grazing of zooplankton between non-toxic phytoplankton (NTP) and toxic phytoplankton (TPP) is a realistic phenomena unexplored so far. The present article is a first step in this direction. A mathematical model of NTP–TPP-zooplankton with constant and variable zooplankton migration is proposed and analyzed. The asymptotic dynamics of the model system around the biologically feasible equilibria is explored through local stability analysis. The dynamics of the proposed system is explored and displayed for different combination of migratory parameters and toxin inhibition parameters. Our analysis suggests that the migratory grazing of zooplankton has a significant role in determining the dynamic stability and oscillation of phytoplankton zooplankton systems.

How well can we measure the ocean’s mean dynamic topography from space?

Recent gravity missions have produced a dramatic improvement in our ability to measure the ocean’s mean dynamic topography (MDT) from space. To fully exploit this oceanic observation, however, we must quantify its error. To establish a baseline, we first assess the error budget for an MDT calculated using a 3rd generation GOCE geoid and the CLS01 mean sea surface (MSS). With these products, we can resolve MDT spatial scales down to 250 km with an accuracy of 1.7 cm, with the MSS and geoid making similar contributions to the total error. For spatial scales within the range 133–250 km the error is 3.0 cm, with the geoid making the greatest contribution. For the smallest resolvable spatial scales (80–133 km) the total error is 16.4 cm, with geoid error accounting for almost all of this. Relative to this baseline, the most recent versions of the geoid and MSS fields reduce the long and short-wavelength errors by 0.9 and 3.2 cm, respectively, but they have little impact in the medium-wavelength band. The newer MSS is responsible for most of the long-wavelength improvement, while for the short-wavelength component it is the geoid. We find that while the formal geoid errors have reasonable global mean values they fail capture the regional variations in error magnitude, which depend on the steepness of the sea floor topography.

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.

Modeling the observed proton aurora and ionospheric convection responses to changes in the IMF clock angle: 1. Persistence of cusp proton aurora

We employ a numerical model of cusp ion precipitation and proton aurora emission to fit variations of the peak Doppler-shifted Lyman-a intensity observed on 26 November 2000 by the SI-12 channel of the FUV instrument on the IMAGE satellite. The major features of this event appeared in response to two brief swings of the interplanetary magnetic field (IMF) toward a southward orientation. We reproduce the observed spatial distributions of this emission on newly opened field lines by combining the proton emission model with a model of the response of ionospheric convection. The simulations are based on the observed variations of the solar wind proton temperature and concentration and the interplanetary magnetic field clock angle. They also allow for the efficiency, sampling rate, integration time and spatial resolution of the FUV instrument. The good match (correlation coefficient 0.91, significant at the 98% level) between observed and modeled variations confirms the time constant (about 4 min) for the rise and decay of the proton emissions predicted by the model for southward IMF conditions. The implications for the detection of pulsed magnetopause reconnection using proton aurora are discussed for a range of interplanetary conditions.

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.

Probabilistic downscaling of remote sensing data with applications for multi-scale biogeochemical flux modeling

Upscaling ecological information to larger scales in space and downscaling remote sensing observations or model simulations to finer scales remain grand challenges in Earth system science. Downscaling often involves inferring subgrid information from coarse-scale data, and such ill-posed problems are classically addressed using regularization. Here, we apply two-dimensional Tikhonov Regularization (2DTR) to simulate subgrid surface patterns for ecological applications. Specifically, we test the ability of 2DTR to simulate the spatial statistics of high-resolution (4 m) remote sensing observations of the normalized difference vegetation index (NDVI) in a tundra landscape. We find that the 2DTR approach as applied here can capture the major mode of spatial variability of the high-resolution information, but not multiple modes of spatial variability, and that the Lagrange multiplier (γ) used to impose the condition of smoothness across space is related to the range of the experimental semivariogram. We used observed and 2DTR-simulated maps of NDVI to estimate landscape-level leaf area index (LAI) and gross primary productivity (GPP). NDVI maps simulated using a γ value that approximates the range of observed NDVI result in a landscape-level GPP estimate that differs by ca 2% from those created using observed NDVI. Following findings that GPP per unit LAI is lower near vegetation patch edges, we simulated vegetation patch edges using multiple approaches and found that simulated GPP declined by up to 12% as a result. 2DTR can generate random landscapes rapidly and can be applied to disaggregate ecological information and compare of spatial observations against simulated landscapes.