SWAT use of gridded observations for simulating runoff – a Vietnam river basin study

Many research studies that focus on basin hydrology have applied the SWAT model using station data to simulate runoff. But over regions lacking robust station data, there is a problem of applying the model to study the hydrological responses. For some countries and remote areas, the rainfall data availability might be a constraint due to many different reasons such as lacking of technology, war time and financial limitation that lead to difficulty in constructing the runoff data. To overcome such a limitation, this research study uses some of the available globally gridded high resolution precipitation datasets to simulate runoff. Five popular gridded observation precipitation datasets: (1) Asian Precipitation Highly Resolved Observational Data Integration Towards the Evaluation of Water Resources (APHRODITE), (2) Tropical Rainfall Measuring Mission (TRMM), (3) Precipitation Estimation from Remote Sensing Information using Artificial Neural Network (PERSIANN), (4) Global Precipitation Climatology Project (GPCP), (5) a modified version of Global Historical Climatology Network (GHCN2) and one reanalysis dataset, National Centers for Environment Prediction/National Center for Atmospheric Research (NCEP/NCAR) are used to simulate runoff over the Dak Bla river (a small tributary of the Mekong River) in Vietnam. Wherever possible, available station data are also used for comparison. Bilinear interpolation of these gridded datasets is used to input the precipitation data at the closest grid points to the station locations. Sensitivity Analysis and Auto-calibration are performed for the SWAT model. The Nash-Sutcliffe Efficiency (NSE) and Coefficient of Determination (<i>R</i><sup>2</sup>) indices are used to benchmark the model performance. Results indicate that the APHRODITE dataset performed very well on a daily scale simulation of discharge having a good NSE of 0.54 and <i>R</i><sup>2</sup> of 0.55, when compared to the discharge simulation using station data (0.68 and 0.71). The GPCP proved to be the next best dataset that was applied to the runoff modelling, with NSE and <i>R</i><sup>2</sup> of 0.46 and 0.51, respectively. The PERSIANN and TRMM rainfall data driven runoff did not show good agreement compared to the station data as both the NSE and <i>R</i><sup>2</sup> indices showed a low value of 0.3. GHCN2 and NCEP also did not show good correlations. The varied results by using these datasets indicate that although the gauge based and satellite-gauge merged products use some ground truth data, the different interpolation techniques and merging algorithms could also be a source of uncertainties. This entails a good understanding of the response of the hydrological model to different datasets and a quantification of the uncertainties in these datasets. Such a methodology is also useful for planning on Rainfall-runoff and even reservoir/river management both at rural and urban scales

Iterative graph cuts for image segmentation with a nonlinear statistical shape prior

Shape-based regularization has proven to be a useful method for delineating objects within noisy images where one has prior knowledge of the shape of the targeted object. When a collection of possible shapes is available, the specification of a shape prior using kernel density estimation is a natural technique. Unfortunately, energy functionals arising from kernel density estimation are of a form that makes them impossible to directly minimize using efficient optimization algorithms such as graph cuts. Our main contribution is to show how one may recast the energy functional into a form that is minimizable iteratively and efficiently using graph cuts.Comment: Revision submitted to JMIV (02/24/13

Commensurate and Incommensurate Vortex States in Superconductors with Periodic Pinning Arrays

As a function of applied field, we find a rich variety of ordered and partially-ordered vortex lattice configurations in systems with square or triangular arrays of pinning sites. We present formulas that predict the matching fields at which commensurate vortex configurations occur and the vortex lattice orientation with respect to the pinning lattice. Our results are in excellent agreement with recent imaging experiments on square pinning arrays [K. Harada et al., Science 274, 1167 (1996)].Comment: 9 pages, 3 figures. Accepted to Physical Review

The Link Between Shocks, Turbulence and Magnetic Reconnection in Collisionless Plasmas

Global hybrid (electron fluid, kinetic ions) and fully kinetic simulations of the magnetosphere have been used to show surprising interconnection between shocks, turbulence and magnetic reconnection. In particular collisionless shocks with their reflected ions that can get upstream before retransmission can generate previously unforeseen phenomena in the post shocked flows: (i) formation of reconnecting current sheets and magnetic islands with sizes up to tens of ion inertial length. (ii) Generation of large scale low frequency electromagnetic waves that are compressed and amplified as they cross the shock. These \u27wavefronts\u27 maintain their integrity for tens of ion cyclotron times but eventually disrupt and dissipate their energy. (iii) Rippling of the shock front, which can in turn lead to formation of fast collimated jets extending to hundreds of ion inertial lengths downstream of the shock. The jets, which have high dynamical pressure, \u27stir\u27 the downstream region, creating large scale disturbances such as vortices, sunward flows, and can trigger flux ropes along the magnetopause. This phenomenology closes the loop between shocks, turbulence and magnetic reconnection in ways previously unrealized. These interconnections appear generic for the collisionless plasmas typical of space, and are expected even at planar shocks, although they will also occur at curved shocks as occur at planets or around ejecta

Role of Esrrg in the Fibrate-Mediated Regulation of Lipid Metabolism Genes in Human ApoA-I Transgenic Mice

We have used a new ApoA-I transgenic mouse model to identify by global gene expression profiling, candidate genes that affect lipid and lipoprotein metabolism in response to fenofibrate treatment. Multilevel bioinformatical analysis and stringent selection criteria (2-fold change, 0% false discovery rate) identified 267 significantly changed genes involved in several molecular pathways. The fenofibrate-treated group did not have significantly altered levels of hepatic human APOA-I mRNA and plasma ApoA-I compared with the control group. However, the treatment increased cholesterol levels to 1.95-fold mainly due to the increase in high-density lipoprotein (HDL) cholesterol. The observed changes in HDL are associated with the upregulation of genes involved in phospholipid biosynthesis and lipid hydrolysis, as well as phospholipid transfer protein. Significant upregulation was observed in genes involved in fatty acid transport and β-oxidation, but not in those of fatty acid and cholesterol biosynthesis, Krebs cycle and gluconeogenesis. Fenofibrate changed significantly the expression of seven transcription factors. The estrogen receptor-related gamma gene was upregulated 2.36-fold and had a significant positive correlation with genes of lipid and lipoprotein metabolism and mitochondrial functions, indicating an important role of this orphan receptor in mediating the fenofibrate-induced activation of a specific subset of its target genes.National Institutes of Health (HL48739 and HL68216); European Union (LSHM-CT-2006-0376331, LSHG-CT-2006-037277); the Biomedical Research Foundation of the Academy of Athens; the Hellenic Cardiological Society; the John F Kostopoulos Foundatio

Smooth analysis of the condition number and the least singular value

Let \a be a complex random variable with mean zero and bounded variance. Let $N_{n}$ be the random matrix of size $n$ whose entries are iid copies of \a and $M$ be a fixed matrix of the same size. The goal of this paper is to give a general estimate for the condition number and least singular value of the matrix $M + N_{n}$, generalizing an earlier result of Spielman and Teng for the case when \a is gaussian. Our investigation reveals an interesting fact that the "core" matrix $M$ does play a role on tail bounds for the least singular value of $M+N_{n}$. This does not occur in Spielman-Teng studies when \a is gaussian. Consequently, our general estimate involves the norm $\|M\|$. In the special case when $\|M\|$ is relatively small, this estimate is nearly optimal and extends or refines existing results.Comment: 20 pages. An erratum to the published version has been adde

Molecular markers for discriminating Streptococcus pyogenes and S. dysgalactiae subspecies equisimilis

Given the increasing aetiological importance of Streptococcus dysgalactiae subspecies equisimilis in diseases which are primarily attributed to S. pyogenes, molecular markers are essential to distinguish these species and delineate their epidemiology more precisely. Many clinical microbiology laboratories rely on agglutination reactivity and biochemical tests to distinguish them. These methods have limitations which are particularly exacerbated when isolates with mixed properties are encountered. In order to provide additional distinguishing parameters that could be used to unequivocally discriminate these two common pathogens, we assess here three molecular targets: the speB gene, intergenic region upstream of the scpG gene (IRSG) and virPCR. Of these, the former two respectively gave positive and negative results for S. pyogenes, and negative and positive results for S. dysgalactiae subsp. equisimilis. Thus, a concerted use of these nucleic acid-based methods is particularly helpful in epidemiological surveillance to accurately assess the relative contribution of these species to streptococcal infections and diseases

Elongation mechanism of the ion shaping of embedded gold nanoparticles under swift heavy ion irradiation

The elongation process under swift heavy ion irradiation (74 MeV Kr ions) of gold NPs, with a diameter in the range 10-30 nm, and embedded in a silica matrix has been investigated by combining experiment and simulation techniques: three-dimensional thermal spike (3DTS), molecular dynamics (MD) and a phenomenological simulation code specially developed for this study. 3DTS simulations evidence the formation of a track in the host matrix and the melting of the NP after the passage of the impinging ion. MD simulations demonstrate that melted NPs have enough time to expand after each ion impact. Our phenomenological simulation relies on the expansion of the melted NP, which flows in the track in silica with modified (lower) density, followed by its recrystallization upon cooling. Finally, the elongation of the spherical NP into a cylindrical one, with a length proportional to its initial size and a width close to the diameter of the track, is the result of the superposition of the independent effects of each expansion/recrystallization process occurring for each ion impact. In agreement with experiment, the simulation shows the gradual elongation of spherical NPs in the ion-beam direction until their widths saturate in the steady state and reach a value close to the track diameter. Moreover, the simulations indicate that the expansion of the gold NP is incomplete at each ion impact.Peer reviewe

L^{2}-restriction bounds for eigenfunctions along curves in the quantum completely integrable case

We show that for a quantum completely integrable system in two dimensions,the $L^{2}$-normalized joint eigenfunctions of the commuting semiclassical pseudodifferential operators satisfy restriction bounds ofthe form $\int_{\gamma} |\phi_{j}^{\hbar}|^2 ds = {\mathcal O}(|\log \hbar|)$ for generic curves $\gamma$ on the surface. We also prove that the maximal restriction bounds of Burq-Gerard-Tzvetkov are always attained for certain exceptional subsequences of eigenfunctions.Comment: Correct some typos and added some more detail in section