Bayesian Optimal Experimental Design Using Multilevel Monte Carlo

KAUST Repository

Ben Issaid, Chaouki; Long, Quan; Scavino, Marco; Tempone, Raul

2015-01-01

Experimental design is very important since experiments are often resource-exhaustive and time-consuming. We carry out experimental design in the Bayesian framework. To measure the amount of information, which can be extracted from the data in an experiment, we use the expected information gain as the utility function, which specifically is the expected logarithmic ratio between the posterior and prior distributions. Optimizing this utility function enables us to design experiments that yield the most informative data for our purpose. One of the major difficulties in evaluating the expected information gain is that the integral is nested and can be high dimensional. We propose using Multilevel Monte Carlo techniques to accelerate the computation of the nested high dimensional integral. The advantages are twofold. First, the Multilevel Monte Carlo can significantly reduce the cost of the nested integral for a given tolerance, by using an optimal sample distribution among different sample averages of the inner integrals. Second, the Multilevel Monte Carlo method imposes less assumptions, such as the concentration of measures, required by Laplace method. We test our Multilevel Monte Carlo technique using a numerical example on the design of sensor deployment for a Darcy flow problem governed by one dimensional Laplace equation. We also compare the performance of the Multilevel Monte Carlo, Laplace approximation and direct double loop Monte Carlo.

Bayesian Optimal Experimental Design Using Multilevel Monte Carlo

KAUST Repository

Ben Issaid, Chaouki

2015-01-07

Experimental design is very important since experiments are often resource-exhaustive and time-consuming. We carry out experimental design in the Bayesian framework. To measure the amount of information, which can be extracted from the data in an experiment, we use the expected information gain as the utility function, which specifically is the expected logarithmic ratio between the posterior and prior distributions. Optimizing this utility function enables us to design experiments that yield the most informative data for our purpose. One of the major difficulties in evaluating the expected information gain is that the integral is nested and can be high dimensional. We propose using Multilevel Monte Carlo techniques to accelerate the computation of the nested high dimensional integral. The advantages are twofold. First, the Multilevel Monte Carlo can significantly reduce the cost of the nested integral for a given tolerance, by using an optimal sample distribution among different sample averages of the inner integrals. Second, the Multilevel Monte Carlo method imposes less assumptions, such as the concentration of measures, required by Laplace method. We test our Multilevel Monte Carlo technique using a numerical example on the design of sensor deployment for a Darcy flow problem governed by one dimensional Laplace equation. We also compare the performance of the Multilevel Monte Carlo, Laplace approximation and direct double loop Monte Carlo.

Le Comte de Monte Cristo: da literatura ao cinema

OpenAIRE

Caravela, Natércia Murta Silva

2008-01-01

A presente dissertação discute o diálogo estabelecido entre literatura e cinema no tratamento da personagem principal – um homem traído que se vinga de forma cruel dos seus inimigos – na obra literária Le Comte de Monte-Cristo, de Alexandre Dumas, e nas três adaptações fílmicas escolhidas: Le Comte de Monte-Cristo de Robert Vernay (1943); The count of Monte Cristo de David Greene (1975) e The count of Monte Cristo de Kevin Reynolds (2002). O projecto centra-se na análise da ...

Physiographic and floristic gradients across topography in transitional seasonally dry evergreen forests of southeast Pará, Brazil Gradientes fisiográficos e florísticos ao longo do relevo em florestas perenifólias sazonalmente secas de transição no Sudeste do Pará, Brasil

Directory of Open Access Journals (Sweden)

James Grogan

2006-12-01

Full Text Available Seasonally dry evergreen forests in southeast Pará, Brazil are transitional between taller closed forests of the interior Amazon Basin and woodland savannas (cerrados of Brazil's south-central plains. We describe abiotic and biotic gradients in this region near the frontier town of Redenção where forest structure and composition grade subtly across barely undulating topography. Annual precipitation averaged 1859 mm between 1995-2001, with nearly zero rainfall during the dry season months of June August. Annual vertical migrations of deep-soil water caused by seasonal rainfall underlie edaphic and floristic differences between high- and low-ground terrain. Low-ground soils are hydromorphic, shaped by perching water tables during the wet season, pale gray, brown, or white in color, with coarse texture, low moisture retention during the dry season, and relatively high macro-nutrient status in the surface horizons. Forest canopies on low ground are highly irregular, especially along seasonal streams, while overstory community composition differs demonstrably from that on high ground. High-ground soils are dystrophic, well-drained through the wet season, brown or red-yellow in color, with finer texture, higher moisture retention, and low macro-nutrient status in the surface horizons compared to low-ground soils. Forest canopies are, on average, taller, more regular, and more closed on high ground. Low-ground areas can be envisioned as energy and nutrient sinks, where, because of hydrologic cycles, canopy disturbance likely occurs more frequently than at high-ground positions if not necessarily at larger scales.As florestas perenifólias sazonalmente secas no sudeste do Pará, Brasil, são áreas de transição entre as florestas fechadas mais altas do interior da Bacia Amazônica e o cerrado das planícies da região Sul-Central do Brasil. Descrevemos os gradientes abióticos e bióticos nessa região próxima da cidade de Redenção, onde a

Present status of transport code development based on Monte Carlo method

International Nuclear Information System (INIS)

Nakagawa, Masayuki

1985-01-01

The present status of development in Monte Carlo code is briefly reviewed. The main items are the followings; Application fields, Methods used in Monte Carlo code (geometry spectification, nuclear data, estimator and variance reduction technique) and unfinished works, Typical Monte Carlo codes and Merits of continuous energy Monte Carlo code. (author)

Fragilidade ambiental na bacia hidrográfica do Alto Parnaíba

OpenAIRE

Alves de Melo, Nivaneide

2007-01-01

Esta pesquisa desenvolveu estudos integrados sobre a bacia hidrográfica do Alto Parnaíba, no Piauí, considerando a atuação das atividades humanas sobre o ambiente natural, a fim de determinar o grau de alteração desse ambiente, a partir da presença antrópica neste local, além de propor ações para restabelecimento de uma situação de equilíbrio ambiental. Pois, quando os eventos sobre a paisagem são de origem antrópica e de orientação econômica, os impactos poderão causar dano...

Estudio farmacoepidemiológico de uso de antiinflamatorios no esteroideos en pacientes de alto riesgo cardiovascular

Directory of Open Access Journals (Sweden)

Jorge Enrique Machado-Alba

Full Text Available Con el objetivo de determinar la frecuencia de uso prolongado de antiinflamatorios no esteroideos (AINES en pacientes colombianos de alto riesgo cardiovascular (ARC se desarrolló un estudio retrospectivo en el cual se identificaron pacientes de ARC que usaron AINES por más de cinco meses continuos entre enero de 2011 y marzo de 2013. Se identificó a los pacientes que recibían crónicamente nitratos, digitálicos, y clopidogrel y ácido acetil salicílico (ASA, quienes fueron identificados como de ARC. Se realizó un análisis de frecuencias de uso según la comedicación recibida. Se encontró uso concomitante de AINES en el 0,35% de los consumidores de nitratos (tiempo promedio: 9,5 ± 4,4 meses, en el 0,36% de los consumidores de clopidogrel y ASA (tiempo promedio: 9,3 ± 3,4 meses, y en el 0,4% de los consumidores de digitálicos (10,2 ± 4,6 meses. Se concluye que existe una baja proporción de uso de AINES de manera crónica en pacientes de alto riesgo cardiovascular.

Vazões máximas e mínimas para bacias hidrográficas da região alto Rio Grande, MG Maximum and minimum discharges for Alto Rio Grande region basins, Minas Gerais state, Brazil

Directory of Open Access Journals (Sweden)

Carlos Rogério de Mello

2010-04-01

Full Text Available Vazões máximas são grandezas hidrológicas aplicadas a projetos de obras hidráulicas e vazões mínimas são utilizadas para a avaliação das disponibilidades hídricas em bacias hidrográficas e comportamento do escoamento subterrâneo. Neste estudo, objetivou-se à construção de intervalos de confiança estatísticos para vazões máximas e mínimas diárias anuais e sua relação com as características fisiográficas das 6 maiores bacias hidrográficas da região Alto Rio Grande à montante da represa da UHE-Camargos/CEMIG. As distribuições de probabilidades Gumbel e Gama foram aplicadas, respectivamente, para séries históricas de vazões máximas e mínimas, utilizando os estimadores de Máxima Verossimilhança. Os intervalos de confiança constituem-se em uma importante ferramenta para o melhor entendimento e estimativa das vazões, sendo influenciado pelas características geológicas das bacias. Com base nos mesmos, verificou-se que a região Alto Rio Grande possui duas áreas distintas: a primeira, abrangendo as bacias Aiuruoca, Carvalhos e Bom Jardim, que apresentaram as maiores vazões máximas e mínimas, significando potencialidade para cheias mais significativas e maiores disponibilidades hídricas; a segunda, associada às bacias F. Laranjeiras, Madre de Deus e Andrelândia, que apresentaram as menores disponibilidades hídricas.Maximum discharges are applied to hydraulic structure design and minimum discharges are used to characterize water availability in hydrographic basins and subterranean flow. This study is aimed at estimating the confidence statistical intervals for maximum and minimum annual discharges and their relationship wih the physical characteristics of basins in the Alto Rio Grande Region, State of Minas Gerais. The study was developed for the six (6 greatest Alto Rio Grande Region basins at upstream of the UHE-Camargos/CEMIG reservoir. Gumbel and Gama probability distribution models were applied to the

Successful vectorization - reactor physics Monte Carlo code

International Nuclear Information System (INIS)

Martin, W.R.

1989-01-01

Most particle transport Monte Carlo codes in use today are based on the ''history-based'' algorithm, wherein one particle history at a time is simulated. Unfortunately, the ''history-based'' approach (present in all Monte Carlo codes until recent years) is inherently scalar and cannot be vectorized. In particular, the history-based algorithm cannot take advantage of vector architectures, which characterize the largest and fastest computers at the current time, vector supercomputers such as the Cray X/MP or IBM 3090/600. However, substantial progress has been made in recent years in developing and implementing a vectorized Monte Carlo algorithm. This algorithm follows portions of many particle histories at the same time and forms the basis for all successful vectorized Monte Carlo codes that are in use today. This paper describes the basic vectorized algorithm along with descriptions of several variations that have been developed by different researchers for specific applications. These applications have been mainly in the areas of neutron transport in nuclear reactor and shielding analysis and photon transport in fusion plasmas. The relative merits of the various approach schemes will be discussed and the present status of known vectorization efforts will be summarized along with available timing results, including results from the successful vectorization of 3-D general geometry, continuous energy Monte Carlo. (orig.)

Bayesian phylogeny analysis via stochastic approximation Monte Carlo

KAUST Repository

Cheon, Sooyoung; Liang, Faming

2009-01-01

in simulating from the posterior distribution of phylogenetic trees, rendering the inference ineffective. In this paper, we apply an advanced Monte Carlo algorithm, the stochastic approximation Monte Carlo algorithm, to Bayesian phylogeny analysis. Our method

Feminicidios en la frontera chilena: el caso de Alto Hospicio

Directory of Open Access Journals (Sweden)

Ainhoa Vásquez Mejías

2016-01-01

Full Text Available Entre 1999 y el 2001, varias adolescentes desaparecieron en el norte de Chile y, posteriormente, fueron encontradas muertas. La novela Alto Hospicio (2008, de Rodrigo Ramos Bañados, recrea estos feminicidios desde la visión del cómplice del único inculpado por estos crímenes. En el presente artículo se plantea que el escritor utiliza la noción de frontera como aparato crítico, con el fin de desestabilizar barreras espaciales, simbólicas y textuales que se agrupan en cuatro ejes: Bolivia-Chile, santas-putas, racionalidad-locura, ficción-realidad. Así, la novela cuestiona y deconstruye los límites para hablarnos de estos feminicidios desde un punto indeterminado, que funciona como espejo de la incertidumbre e irresolución de este caso en la justicia chilena.

Reflections on early Monte Carlo calculations

International Nuclear Information System (INIS)

Spanier, J.

1992-01-01

Monte Carlo methods for solving various particle transport problems developed in parallel with the evolution of increasingly sophisticated computer programs implementing diffusion theory and low-order moments calculations. In these early years, Monte Carlo calculations and high-order approximations to the transport equation were seen as too expensive to use routinely for nuclear design but served as invaluable aids and supplements to design with less expensive tools. The earliest Monte Carlo programs were quite literal; i.e., neutron and other particle random walk histories were simulated by sampling from the probability laws inherent in the physical system without distoration. Use of such analogue sampling schemes resulted in a good deal of time being spent in examining the possibility of lowering the statistical uncertainties in the sample estimates by replacing simple, and intuitively obvious, random variables by those with identical means but lower variances

Forests

Science.gov (United States)

Louis R. Iverson; Mark W. Schwartz

1994-01-01

Originally diminished by development, forests are coming back: forest biomass is accumulating. Forests are repositories for many threatened species. Even with increased standing timber, however, biodiversity is threatened by increased forest fragmentation and by exotic species.

Reconstruction of Monte Carlo replicas from Hessian parton distributions

Energy Technology Data Exchange (ETDEWEB)

Hou, Tie-Jiun [Department of Physics, Southern Methodist University,Dallas, TX 75275-0181 (United States); Gao, Jun [INPAC, Shanghai Key Laboratory for Particle Physics and Cosmology,Department of Physics and Astronomy, Shanghai Jiao-Tong University, Shanghai 200240 (China); High Energy Physics Division, Argonne National Laboratory,Argonne, Illinois, 60439 (United States); Huston, Joey [Department of Physics and Astronomy, Michigan State University,East Lansing, MI 48824 (United States); Nadolsky, Pavel [Department of Physics, Southern Methodist University,Dallas, TX 75275-0181 (United States); Schmidt, Carl; Stump, Daniel [Department of Physics and Astronomy, Michigan State University,East Lansing, MI 48824 (United States); Wang, Bo-Ting; Xie, Ke Ping [Department of Physics, Southern Methodist University,Dallas, TX 75275-0181 (United States); Dulat, Sayipjamal [Department of Physics and Astronomy, Michigan State University,East Lansing, MI 48824 (United States); School of Physics Science and Technology, Xinjiang University,Urumqi, Xinjiang 830046 (China); Center for Theoretical Physics, Xinjiang University,Urumqi, Xinjiang 830046 (China); Pumplin, Jon; Yuan, C.P. [Department of Physics and Astronomy, Michigan State University,East Lansing, MI 48824 (United States)

2017-03-20

We explore connections between two common methods for quantifying the uncertainty in parton distribution functions (PDFs), based on the Hessian error matrix and Monte-Carlo sampling. CT14 parton distributions in the Hessian representation are converted into Monte-Carlo replicas by a numerical method that reproduces important properties of CT14 Hessian PDFs: the asymmetry of CT14 uncertainties and positivity of individual parton distributions. The ensembles of CT14 Monte-Carlo replicas constructed this way at NNLO and NLO are suitable for various collider applications, such as cross section reweighting. Master formulas for computation of asymmetric standard deviations in the Monte-Carlo representation are derived. A correction is proposed to address a bias in asymmetric uncertainties introduced by the Taylor series approximation. A numerical program is made available for conversion of Hessian PDFs into Monte-Carlo replicas according to normal, log-normal, and Watt-Thorne sampling procedures.

Sampling from a polytope and hard-disk Monte Carlo

International Nuclear Information System (INIS)

Kapfer, Sebastian C; Krauth, Werner

2013-01-01

The hard-disk problem, the statics and the dynamics of equal two-dimensional hard spheres in a periodic box, has had a profound influence on statistical and computational physics. Markov-chain Monte Carlo and molecular dynamics were first discussed for this model. Here we reformulate hard-disk Monte Carlo algorithms in terms of another classic problem, namely the sampling from a polytope. Local Markov-chain Monte Carlo, as proposed by Metropolis et al. in 1953, appears as a sequence of random walks in high-dimensional polytopes, while the moves of the more powerful event-chain algorithm correspond to molecular dynamics evolution. We determine the convergence properties of Monte Carlo methods in a special invariant polytope associated with hard-disk configurations, and the implications for convergence of hard-disk sampling. Finally, we discuss parallelization strategies for event-chain Monte Carlo and present results for a multicore implementation

Problems in radiation shielding calculations with Monte Carlo methods

International Nuclear Information System (INIS)

Ueki, Kohtaro

1985-01-01

The Monte Carlo method is a very useful tool for solving a large class of radiation transport problem. In contrast with deterministic method, geometric complexity is a much less significant problem for Monte Carlo calculations. However, the accuracy of Monte Carlo calculations is of course, limited by statistical error of the quantities to be estimated. In this report, we point out some typical problems to solve a large shielding system including radiation streaming. The Monte Carlo coupling technique was developed to settle such a shielding problem accurately. However, the variance of the Monte Carlo results using the coupling technique of which detectors were located outside the radiation streaming, was still not enough. So as to bring on more accurate results for the detectors located outside the streaming and also for a multi-legged-duct streaming problem, a practicable way of ''Prism Scattering technique'' is proposed in the study. (author)

Cluster monte carlo method for nuclear criticality safety calculation

International Nuclear Information System (INIS)

Pei Lucheng

1984-01-01

One of the most important applications of the Monte Carlo method is the calculation of the nuclear criticality safety. The fair source game problem was presented at almost the same time as the Monte Carlo method was applied to calculating the nuclear criticality safety. The source iteration cost may be reduced as much as possible or no need for any source iteration. This kind of problems all belongs to the fair source game prolems, among which, the optimal source game is without any source iteration. Although the single neutron Monte Carlo method solved the problem without the source iteration, there is still quite an apparent shortcoming in it, that is, it solves the problem without the source iteration only in the asymptotic sense. In this work, a new Monte Carlo method called the cluster Monte Carlo method is given to solve the problem further

Wielandt acceleration for MCNP5 Monte Carlo eigenvalue calculations

International Nuclear Information System (INIS)

Brown, F.

2007-01-01

Monte Carlo criticality calculations use the power iteration method to determine the eigenvalue (k eff ) and eigenfunction (fission source distribution) of the fundamental mode. A recently proposed method for accelerating convergence of the Monte Carlo power iteration using Wielandt's method has been implemented in a test version of MCNP5. The method is shown to provide dramatic improvements in convergence rates and to greatly reduce the possibility of false convergence assessment. The method is effective and efficient, improving the Monte Carlo figure-of-merit for many problems. In addition, the method should eliminate most of the underprediction bias in confidence intervals for Monte Carlo criticality calculations. (authors)

Monte Carlo shielding analyses using an automated biasing procedure

International Nuclear Information System (INIS)

Tang, J.S.; Hoffman, T.J.

1988-01-01

A systematic and automated approach for biasing Monte Carlo shielding calculations is described. In particular, adjoint fluxes from a one-dimensional discrete ordinates calculation are used to generate biasing parameters for a Monte Carlo calculation. The entire procedure of adjoint calculation, biasing parameters generation, and Monte Carlo calculation has been automated. The automated biasing procedure has been applied to several realistic deep-penetration shipping cask problems. The results obtained for neutron and gamma-ray transport indicate that with the automated biasing procedure Monte Carlo shielding calculations of spent-fuel casks can be easily performed with minimum effort and that accurate results can be obtained at reasonable computing cost

Applications of the Monte Carlo method in radiation protection

International Nuclear Information System (INIS)

Kulkarni, R.N.; Prasad, M.A.

1999-01-01

This paper gives a brief introduction to the application of the Monte Carlo method in radiation protection. It may be noted that an exhaustive review has not been attempted. The special advantage of the Monte Carlo method has been first brought out. The fundamentals of the Monte Carlo method have next been explained in brief, with special reference to two applications in radiation protection. Some sample current applications have been reported in the end in brief as examples. They are, medical radiation physics, microdosimetry, calculations of thermoluminescence intensity and probabilistic safety analysis. The limitations of the Monte Carlo method have also been mentioned in passing. (author)

Pore-scale uncertainty quantification with multilevel Monte Carlo

KAUST Repository

Icardi, Matteo; Hoel, Haakon; Long, Quan; Tempone, Raul

2014-01-01

. Since there are no generic ways to parametrize the randomness in the porescale structures, Monte Carlo techniques are the most accessible to compute statistics. We propose a multilevel Monte Carlo (MLMC) technique to reduce the computational cost

Forest structure in low-diversity tropical forests: a study of Hawaiian wet and dry forests.

Science.gov (United States)

Ostertag, Rebecca; Inman-Narahari, Faith; Cordell, Susan; Giardina, Christian P; Sack, Lawren

2014-01-01

The potential influence of diversity on ecosystem structure and function remains a topic of significant debate, especially for tropical forests where diversity can range widely. We used Center for Tropical Forest Science (CTFS) methodology to establish forest dynamics plots in montane wet forest and lowland dry forest on Hawai'i Island. We compared the species diversity, tree density, basal area, biomass, and size class distributions between the two forest types. We then examined these variables across tropical forests within the CTFS network. Consistent with other island forests, the Hawai'i forests were characterized by low species richness and very high relative dominance. The two Hawai'i forests were floristically distinct, yet similar in species richness (15 vs. 21 species) and stem density (3078 vs. 3486/ha). While these forests were selected for their low invasive species cover relative to surrounding forests, both forests averaged 5->50% invasive species cover; ongoing removal will be necessary to reduce or prevent competitive impacts, especially from woody species. The montane wet forest had much larger trees, resulting in eightfold higher basal area and above-ground biomass. Across the CTFS network, the Hawaiian montane wet forest was similar to other tropical forests with respect to diameter distributions, density, and aboveground biomass, while the Hawai'i lowland dry forest was similar in density to tropical forests with much higher diversity. These findings suggest that forest structural variables can be similar across tropical forests independently of species richness. The inclusion of low-diversity Pacific Island forests in the CTFS network provides an ∼80-fold range in species richness (15-1182 species), six-fold variation in mean annual rainfall (835-5272 mm yr(-1)) and 1.8-fold variation in mean annual temperature (16.0-28.4°C). Thus, the Hawaiian forest plots expand the global forest plot network to enable testing of ecological theory for

Forest structure in low-diversity tropical forests: a study of Hawaiian wet and dry forests.

Directory of Open Access Journals (Sweden)

Rebecca Ostertag

Full Text Available The potential influence of diversity on ecosystem structure and function remains a topic of significant debate, especially for tropical forests where diversity can range widely. We used Center for Tropical Forest Science (CTFS methodology to establish forest dynamics plots in montane wet forest and lowland dry forest on Hawai'i Island. We compared the species diversity, tree density, basal area, biomass, and size class distributions between the two forest types. We then examined these variables across tropical forests within the CTFS network. Consistent with other island forests, the Hawai'i forests were characterized by low species richness and very high relative dominance. The two Hawai'i forests were floristically distinct, yet similar in species richness (15 vs. 21 species and stem density (3078 vs. 3486/ha. While these forests were selected for their low invasive species cover relative to surrounding forests, both forests averaged 5->50% invasive species cover; ongoing removal will be necessary to reduce or prevent competitive impacts, especially from woody species. The montane wet forest had much larger trees, resulting in eightfold higher basal area and above-ground biomass. Across the CTFS network, the Hawaiian montane wet forest was similar to other tropical forests with respect to diameter distributions, density, and aboveground biomass, while the Hawai'i lowland dry forest was similar in density to tropical forests with much higher diversity. These findings suggest that forest structural variables can be similar across tropical forests independently of species richness. The inclusion of low-diversity Pacific Island forests in the CTFS network provides an ∼80-fold range in species richness (15-1182 species, six-fold variation in mean annual rainfall (835-5272 mm yr(-1 and 1.8-fold variation in mean annual temperature (16.0-28.4°C. Thus, the Hawaiian forest plots expand the global forest plot network to enable testing of ecological

CRECIMIENTO DEL MAÍZ EN VERTISOLES CON ALTO ALUMINIO EN LA BAIXADA MARANHENSE PRE-AMAZONIA, BRASIL

Directory of Open Access Journals (Sweden)

Alessandro Costa da Silva

2014-01-01

Full Text Available Crecimiento del maíz en vertisoles con alto aluminio en la Baixada Maranhense pre-Amazonia, Brasil. El objetivo de este trabajo fue evaluar el crecimiento del maíz en suelos con alto contenido de aluminio. Se midió el efecto del Al3+ en raíces y la cantidad de materia seca (raíz, hoja y tallo de maíz. Se efectuó la caracterización físico-química de cuatro muestras de suelo con alto aluminio colectadas del horizonte Ap, en tres municipios de la región conocida como Baixada Maranhense (Pre-Amazonia, Brasil: Santa Rita (SR, Arari (AR y Vitoria do Mearim (VM y un testigo colectado en el municipio de São Luís, Área del Núcleo de Tecnología Rural (T. El estudio, ejecutado en 2009, se llevó a cabo en invernadero y se utilizó 2 dm3 de suelo por maceta. Asimismo las muestras fueron divididas en muestras con y sin fertilización. La variación en la longitud de la raíz y de materia seca de las hojas difirió significativamente entre tratados con y sin fertilizante, excepto en la muestra de la localidad T. La producción de materia seca de raíz, tallo y hoja fue mayor en todos los suelos cuando se fertilizó. El suelo testigo también superó a todos los demás en cuanto a producción de materia seca en la raíz, posiblemente como resultado de una menor cantidad de Al3+ (1,2 cmolc/dm3 en comparación con los suelos SR, AR y VM (6,8; 8,0 y 7,0 cmolc/dm3 respectivamente. Se concluye la fertilización reduce el efecto detrimental del aluminio en la producción de maíz en la Baixada Maranhense.

Current and future applications of Monte Carlo

International Nuclear Information System (INIS)

Zaidi, H.

2003-01-01

Full text: The use of radionuclides in medicine has a long history and encompasses a large area of applications including diagnosis and radiation treatment of cancer patients using either external or radionuclide radiotherapy. The 'Monte Carlo method'describes a very broad area of science, in which many processes, physical systems, and phenomena are simulated by statistical methods employing random numbers. The general idea of Monte Carlo analysis is to create a model, which is as similar as possible to the real physical system of interest, and to create interactions within that system based on known probabilities of occurrence, with random sampling of the probability density functions (pdfs). As the number of individual events (called 'histories') is increased, the quality of the reported average behavior of the system improves, meaning that the statistical uncertainty decreases. The use of the Monte Carlo method to simulate radiation transport has become the most accurate means of predicting absorbed dose distributions and other quantities of interest in the radiation treatment of cancer patients using either external or radionuclide radiotherapy. The same trend has occurred for the estimation of the absorbed dose in diagnostic procedures using radionuclides as well as the assessment of image quality and quantitative accuracy of radionuclide imaging. As a consequence of this generalized use, many questions are being raised primarily about the need and potential of Monte Carlo techniques, but also about how accurate it really is, what would it take to apply it clinically and make it available widely to the nuclear medicine community at large. Many of these questions will be answered when Monte Carlo techniques are implemented and used for more routine calculations and for in-depth investigations. In this paper, the conceptual role of the Monte Carlo method is briefly introduced and followed by a survey of its different applications in diagnostic and therapeutic

PANIC: A General-purpose Panoramic Near-infrared Camera for the Calar Alto Observatory

Science.gov (United States)

Cárdenas Vázquez, M.-C.; Dorner, B.; Huber, A.; Sánchez-Blanco, E.; Alter, M.; Rodríguez Gómez, J. F.; Bizenberger, P.; Naranjo, V.; Ibáñez Mengual, J.-M.; Panduro, J.; García Segura, A. J.; Mall, U.; Fernández, M.; Laun, W.; Ferro Rodríguez, I. M.; Helmling, J.; Terrón, V.; Meisenheimer, K.; Fried, J. W.; Mathar, R. J.; Baumeister, H.; Rohloff, R.-R.; Storz, C.; Verdes-Montenegro, L.; Bouy, H.; Ubierna, M.; Fopp, P.; Funke, B.

2018-02-01

PANIC7 is the new PAnoramic Near-Infrared Camera for Calar Alto and is a project jointly developed by the MPIA in Heidelberg, Germany, and the IAA in Granada, Spain, for the German-Spanish Astronomical Center at Calar Alto Observatory (CAHA; Almería, Spain). This new instrument works with the 2.2 m and 3.5 m CAHA telescopes covering a field of view of 30 × 30 arcmin and 15 × 15 arcmin, respectively, with a sampling of 4096 × 4096 pixels. It is designed for the spectral bands from Z to K S , and can also be equipped with narrowband filters. The instrument was delivered to the observatory in 2014 October and was commissioned at both telescopes between 2014 November and 2015 June. Science verification at the 2.2 m telescope was carried out during the second semester of 2015 and the instrument is now at full operation. We describe the design, assembly, integration, and verification process, the final laboratory tests and the PANIC instrument performance. We also present first-light data obtained during the commissioning and preliminary results of the scientific verification. The final optical model and the theoretical performance of the camera were updated according to the as-built data. The laboratory tests were made with a star simulator. Finally, the commissioning phase was done at both telescopes to validate the camera real performance on sky. The final laboratory test confirmed the expected camera performances, complying with the scientific requirements. The commissioning phase on sky has been accomplished.

Nuevos hallazgos de arte rupestre en el alto Duero

Directory of Open Access Journals (Sweden)

Teógenes ORTEGO

2009-11-01

Full Text Available La riqueza del importante foco de pinturas rupestres prehistóricas, descubiertas en agosto de 1951 en el monte Valonsadero, de Soria, va aumentando en extensión y variedad a medida que el tiempo nos permite nuevas exploraciones.Uno de los covachos con muestras de este arte, últimamente descubierto, se encuentra en el lugar denominado "Cuerda del Torilejo' de cuyo estudio anticipamos esta breve información.

FLORÍSTICA E ESTRUTURA EM ÁREA DE FLORESTA ESTACIONAL DECIDUAL NA REGIÃO DO ALTO URUGUAI, RS / FLORISTIC AND STRUCTURE IN AREA OF SEASONAL DECIDUOUS FOREST IN THE REGION OF THE ALTO URUGUAI, RS

Directory of Open Access Journals (Sweden)

F. R. Lambrecht

2016-06-01

Full Text Available           O trabalho teve como objetivos realizar a análise florística das estruturas horizontal e vertical em um fragmento de Floresta Estacional Decidual (FED na região do Alto Uruguai, em Frederico Westphalen, RS, Brasil. Foi realizado um censo em uma área amostral de 10.000 m², subdividida em 100 subparcelas de 10x10 m. Realizou-se a mensuração e a identificação de todos os indivíduos de espécies arbóreas com circunferência a altura do peito (CAP≥ 31,4cm. Com os dados obtidos, foram realizados cálculos de densidade (D, frequência (F, dominância (Do, índice de valor de importância (IVI, índice de valor de cobertura (IVC e a estratificação dos indivíduos. O maior número de espécies concentrou-se na família Fabaceae (10 e o maior número de indivíduos na família Meliaceae (147 ou 25,3%. Trichilia claussenii C. DC. Nectandra megapotamica (Spreng. Mez, Achatocarpus praecox Griseb. e Eugenia rostrifolia D. Legrand representaram 38,28% da densidade relativa e Holocalyx balansae Micheli, Trichilia claussenii, Nectandra megapotamica, Cordia americana (L. Gottshling & J.E.Mill., Phytolacca dioica  L., Achatocarpus praecox, e Eugenia rostrifolia representaram 50,5% da dominância relativa. Considerando a estrutura vertical, as sete espécies mais frequentes apresentaram uma maior quantidade de indivíduos no estrato médio. Conclui-se que a área florestal encontra-se em estágio intermediário de sucessão.

Quantum statistical Monte Carlo methods and applications to spin systems

International Nuclear Information System (INIS)

Suzuki, M.

1986-01-01

A short review is given concerning the quantum statistical Monte Carlo method based on the equivalence theorem that d-dimensional quantum systems are mapped onto (d+1)-dimensional classical systems. The convergence property of this approximate tansformation is discussed in detail. Some applications of this general appoach to quantum spin systems are reviewed. A new Monte Carlo method, ''thermo field Monte Carlo method,'' is presented, which is an extension of the projection Monte Carlo method at zero temperature to that at finite temperatures

Sentidos de Vitória/Derrota para os Pais Segundo Atletas do Alto Rendimento

OpenAIRE

Amblard, Isabela; Cruz, Fatima Leite

2015-01-01

Este estudo compreendeu as representações sociais da vitória/derrota para os pais segundo atletas-adolescentes do esporte de alto rendimento, na cidade do Recife. A adolescência é compreendida a partir da Psicologia Social-histórica, e o embasamento teórico-metodológico da Teoria das Representações Sociais abordou os sujeitos em diferentes contextos socioculturais, lugares de pertencimento, experiências, crenças, saberes e sentimentos compartilhados. Adotou-se a perspectiva pluri metodológica...

SPQR: a Monte Carlo reactor kinetics code

International Nuclear Information System (INIS)

Cramer, S.N.; Dodds, H.L.

1980-02-01

The SPQR Monte Carlo code has been developed to analyze fast reactor core accident problems where conventional methods are considered inadequate. The code is based on the adiabatic approximation of the quasi-static method. This initial version contains no automatic material motion or feedback. An existing Monte Carlo code is used to calculate the shape functions and the integral quantities needed in the kinetics module. Several sample problems have been devised and analyzed. Due to the large statistical uncertainty associated with the calculation of reactivity in accident simulations, the results, especially at later times, differ greatly from deterministic methods. It was also found that in large uncoupled systems, the Monte Carlo method has difficulty in handling asymmetric perturbations

Optix: A Monte Carlo scintillation light transport code

Energy Technology Data Exchange (ETDEWEB)

Safari, M.J., E-mail: mjsafari@aut.ac.ir [Department of Energy Engineering and Physics, Amir Kabir University of Technology, PO Box 15875-4413, Tehran (Iran, Islamic Republic of); Afarideh, H. [Department of Energy Engineering and Physics, Amir Kabir University of Technology, PO Box 15875-4413, Tehran (Iran, Islamic Republic of); Ghal-Eh, N. [School of Physics, Damghan University, PO Box 36716-41167, Damghan (Iran, Islamic Republic of); Davani, F. Abbasi [Nuclear Engineering Department, Shahid Beheshti University, PO Box 1983963113, Tehran (Iran, Islamic Republic of)

2014-02-11

The paper reports on the capabilities of Monte Carlo scintillation light transport code Optix, which is an extended version of previously introduced code Optics. Optix provides the user a variety of both numerical and graphical outputs with a very simple and user-friendly input structure. A benchmarking strategy has been adopted based on the comparison with experimental results, semi-analytical solutions, and other Monte Carlo simulation codes to verify various aspects of the developed code. Besides, some extensive comparisons have been made against the tracking abilities of general-purpose MCNPX and FLUKA codes. The presented benchmark results for the Optix code exhibit promising agreements. -- Highlights: • Monte Carlo simulation of scintillation light transport in 3D geometry. • Evaluation of angular distribution of detected photons. • Benchmark studies to check the accuracy of Monte Carlo simulations.

Bayesian phylogeny analysis via stochastic approximation Monte Carlo

KAUST Repository

Cheon, Sooyoung

2009-11-01

Monte Carlo methods have received much attention in the recent literature of phylogeny analysis. However, the conventional Markov chain Monte Carlo algorithms, such as the Metropolis-Hastings algorithm, tend to get trapped in a local mode in simulating from the posterior distribution of phylogenetic trees, rendering the inference ineffective. In this paper, we apply an advanced Monte Carlo algorithm, the stochastic approximation Monte Carlo algorithm, to Bayesian phylogeny analysis. Our method is compared with two popular Bayesian phylogeny software, BAMBE and MrBayes, on simulated and real datasets. The numerical results indicate that our method outperforms BAMBE and MrBayes. Among the three methods, SAMC produces the consensus trees which have the highest similarity to the true trees, and the model parameter estimates which have the smallest mean square errors, but costs the least CPU time. © 2009 Elsevier Inc. All rights reserved.

Mejoras tecnológicas en el proceso de inyección de carbón pulverizado en el horno alto

Directory of Open Access Journals (Sweden)

Babich, A.

1996-04-01

Full Text Available Blast furnace operation technology with pulverized coal injection (PCI has been carried out in order to replace the maximum amount of metallurgical coke in the burden. PC burning has been studied and methods and designs for the intensification of its combustion have been developed within the raceway. Recommendations for optimizing coal grinding have been developed. Problems of PC distribution inside the circumference of the furnace were investigated. Compensating technology for changes in reduction, heat exchange and other processes under PCI have been developed.

Se estudia la tecnología de operación del horno alto con inyección de carbón pulverizado (ICP con el fin de sustituir el máximo posible de coque siderúrgico en la carga y se analizan aspectos de la combustión del carbón pulverizado (CP y los métodos y dispositivos elaborados para intensificar su combustión en la zona de toberas del horno alto. Se ofrecen recomendaciones para optimizar la molienda del carbón, se estudia la problemática de la distribución del CP a lo largo de la periferia del horno y se describe una tecnología elaborada para compensar las perturbaciones producidas por la inyección de CP en el horno alto sobre diferentes parámetros, tales como grado de reducción directa, intercambio térmico y otras.

Forest ownership dynamics of southern forests

Science.gov (United States)

Brett J. Butler; David N. Wear

2013-01-01

Key FindingsPrivate landowners hold 86 percent of the forest area in the South; two-thirds of this area is owned by families or individuals.Fifty-nine percent of family forest owners own between 1 and 9 acres of forest land, but 60 percent of family-owned forests are in holdings of 100 acres or more.Two-...

Present status and future prospects of neutronics Monte Carlo

International Nuclear Information System (INIS)

Gelbard, E.M.

1990-01-01

It is fair to say that the Monte Carlo method, over the last decade, has grown steadily more important as a neutronics computational tool. Apparently this has happened for assorted reasons. Thus, for example, as the power of computers has increased, the cost of the method has dropped, steadily becoming less and less of an obstacle to its use. In addition, more and more sophisticated input processors have now made it feasible to model extremely complicated systems routinely with really remarkable fidelity. Finally, as we demand greater and greater precision in reactor calculations, Monte Carlo is often found to be the only method accurate enough for use in benchmarking. Cross section uncertainties are now almost the only inherent limitations in our Monte Carlo capabilities. For this reason Monte Carlo has come to occupy a special position, interposed between experiment and other computational techniques. More and more often deterministic methods are tested by comparison with Monte Carlo, and cross sections are tested by comparing Monte Carlo with experiment. In this way one can distinguish very clearly between errors due to flaws in our numerical methods, and those due to deficiencies in cross section files. The special role of Monte Carlo as a benchmarking tool, often the only available benchmarking tool, makes it crucially important that this method should be polished to perfection. Problems relating to Eigenvalue calculations, variance reduction and the use of advanced computers are reviewed in this paper. (author)

Diffusion Monte Carlo approach versus adiabatic computation for local Hamiltonians

Science.gov (United States)

Bringewatt, Jacob; Dorland, William; Jordan, Stephen P.; Mink, Alan

2018-02-01

Most research regarding quantum adiabatic optimization has focused on stoquastic Hamiltonians, whose ground states can be expressed with only real non-negative amplitudes and thus for whom destructive interference is not manifest. This raises the question of whether classical Monte Carlo algorithms can efficiently simulate quantum adiabatic optimization with stoquastic Hamiltonians. Recent results have given counterexamples in which path-integral and diffusion Monte Carlo fail to do so. However, most adiabatic optimization algorithms, such as for solving MAX-k -SAT problems, use k -local Hamiltonians, whereas our previous counterexample for diffusion Monte Carlo involved n -body interactions. Here we present a 6-local counterexample which demonstrates that even for these local Hamiltonians there are cases where diffusion Monte Carlo cannot efficiently simulate quantum adiabatic optimization. Furthermore, we perform empirical testing of diffusion Monte Carlo on a standard well-studied class of permutation-symmetric tunneling problems and similarly find large advantages for quantum optimization over diffusion Monte Carlo.

Frequency domain Monte Carlo simulation method for cross power spectral density driven by periodically pulsed spallation neutron source using complex-valued weight Monte Carlo

International Nuclear Information System (INIS)

Yamamoto, Toshihiro

2014-01-01

Highlights: • The cross power spectral density in ADS has correlated and uncorrelated components. • A frequency domain Monte Carlo method to calculate the uncorrelated one is developed. • The method solves the Fourier transformed transport equation. • The method uses complex-valued weights to solve the equation. • The new method reproduces well the CPSDs calculated with time domain MC method. - Abstract: In an accelerator driven system (ADS), pulsed spallation neutrons are injected at a constant frequency. The cross power spectral density (CPSD), which can be used for monitoring the subcriticality of the ADS, is composed of the correlated and uncorrelated components. The uncorrelated component is described by a series of the Dirac delta functions that occur at the integer multiples of the pulse repetition frequency. In the present paper, a Monte Carlo method to solve the Fourier transformed neutron transport equation with a periodically pulsed neutron source term has been developed to obtain the CPSD in ADSs. Since the Fourier transformed flux is a complex-valued quantity, the Monte Carlo method introduces complex-valued weights to solve the Fourier transformed equation. The Monte Carlo algorithm used in this paper is similar to the one that was developed by the author of this paper to calculate the neutron noise caused by cross section perturbations. The newly-developed Monte Carlo algorithm is benchmarked to the conventional time domain Monte Carlo simulation technique. The CPSDs are obtained both with the newly-developed frequency domain Monte Carlo method and the conventional time domain Monte Carlo method for a one-dimensional infinite slab. The CPSDs obtained with the frequency domain Monte Carlo method agree well with those with the time domain method. The higher order mode effects on the CPSD in an ADS with a periodically pulsed neutron source are discussed

Liquefaction Hazard Maps for Three Earthquake Scenarios for the Communities of San Jose, Campbell, Cupertino, Los Altos, Los Gatos, Milpitas, Mountain View, Palo Alto, Santa Clara, Saratoga, and Sunnyvale, Northern Santa Clara County, California

Science.gov (United States)

Holzer, Thomas L.; Noce, Thomas E.; Bennett, Michael J.

2008-01-01

Maps showing the probability of surface manifestations of liquefaction in the northern Santa Clara Valley were prepared with liquefaction probability curves. The area includes the communities of San Jose, Campbell, Cupertino, Los Altos, Los Gatos Milpitas, Mountain View, Palo Alto, Santa Clara, Saratoga, and Sunnyvale. The probability curves were based on complementary cumulative frequency distributions of the liquefaction potential index (LPI) for surficial geologic units in the study area. LPI values were computed with extensive cone penetration test soundings. Maps were developed for three earthquake scenarios, an M7.8 on the San Andreas Fault comparable to the 1906 event, an M6.7 on the Hayward Fault comparable to the 1868 event, and an M6.9 on the Calaveras Fault. Ground motions were estimated with the Boore and Atkinson (2008) attenuation relation. Liquefaction is predicted for all three events in young Holocene levee deposits along the major creeks. Liquefaction probabilities are highest for the M7.8 earthquake, ranging from 0.33 to 0.37 if a 1.5-m deep water table is assumed, and 0.10 to 0.14 if a 5-m deep water table is assumed. Liquefaction probabilities of the other surficial geologic units are less than 0.05. Probabilities for the scenario earthquakes are generally consistent with observations during historical earthquakes.

Neutron point-flux calculation by Monte Carlo

International Nuclear Information System (INIS)

Eichhorn, M.

1986-04-01

A survey of the usual methods for estimating flux at a point is given. The associated variance-reducing techniques in direct Monte Carlo games are explained. The multigroup Monte Carlo codes MC for critical systems and PUNKT for point source-point detector-systems are represented, and problems in applying the codes to practical tasks are discussed. (author)

Technical and scientific report of the Alto project

International Nuclear Information System (INIS)

Essabaa, S.; Gardes, D.; Grialou, D.; Ibrahim, F.; Le Scornet, J.C.

2002-01-01

The Alto project means the installation of an electron linear accelerator inside the experimental area of the tandem accelerator of the nuclear physics institute of Orsay (IPNO, France). This linear accelerator comes from CERN where it was operating as a pre-injector for LEP. This equipment will allow IPNO'teams to perform fast kinetics studies in a domain different from that of ELYSE accelerator. The time resolution will not be as high as that of ELYSE (picosecond) but will be sufficient (microsecond) to produce free radicals in aqueous and gaseous media. The main expectations of this installation can be classified according 3 axis: 1) basic research (mainly the study of nuclear matter through photo-fission, 2) research and development of accelerators (by providing a test bench for new high frequency systems and superconducting components), and 3) applied research for industry concerning: biochemistry under irradiation, radiation sensibility, DNA breaking, food and drug sterilization and behaviour of electronic components under irradiation. This rapport details the research program that could be achieved with this equipment, describes its contributions in terms of economic development, cooperation with industry, student training, and specifies the needed investment and the operating and maintenance costs. (A.C.)

Research on perturbation based Monte Carlo reactor criticality search

International Nuclear Information System (INIS)

Li Zeguang; Wang Kan; Li Yangliu; Deng Jingkang

2013-01-01

Criticality search is a very important aspect in reactor physics analysis. Due to the advantages of Monte Carlo method and the development of computer technologies, Monte Carlo criticality search is becoming more and more necessary and feasible. Traditional Monte Carlo criticality search method is suffered from large amount of individual criticality runs and uncertainty and fluctuation of Monte Carlo results. A new Monte Carlo criticality search method based on perturbation calculation is put forward in this paper to overcome the disadvantages of traditional method. By using only one criticality run to get initial k_e_f_f and differential coefficients of concerned parameter, the polynomial estimator of k_e_f_f changing function is solved to get the critical value of concerned parameter. The feasibility of this method was tested. The results show that the accuracy and efficiency of perturbation based criticality search method are quite inspiring and the method overcomes the disadvantages of traditional one. (authors)

Shell model the Monte Carlo way

International Nuclear Information System (INIS)

Ormand, W.E.

1995-01-01

The formalism for the auxiliary-field Monte Carlo approach to the nuclear shell model is presented. The method is based on a linearization of the two-body part of the Hamiltonian in an imaginary-time propagator using the Hubbard-Stratonovich transformation. The foundation of the method, as applied to the nuclear many-body problem, is discussed. Topics presented in detail include: (1) the density-density formulation of the method, (2) computation of the overlaps, (3) the sign of the Monte Carlo weight function, (4) techniques for performing Monte Carlo sampling, and (5) the reconstruction of response functions from an imaginary-time auto-correlation function using MaxEnt techniques. Results obtained using schematic interactions, which have no sign problem, are presented to demonstrate the feasibility of the method, while an extrapolation method for realistic Hamiltonians is presented. In addition, applications at finite temperature are outlined

Shell model the Monte Carlo way

Energy Technology Data Exchange (ETDEWEB)

Ormand, W.E.

1995-03-01

The formalism for the auxiliary-field Monte Carlo approach to the nuclear shell model is presented. The method is based on a linearization of the two-body part of the Hamiltonian in an imaginary-time propagator using the Hubbard-Stratonovich transformation. The foundation of the method, as applied to the nuclear many-body problem, is discussed. Topics presented in detail include: (1) the density-density formulation of the method, (2) computation of the overlaps, (3) the sign of the Monte Carlo weight function, (4) techniques for performing Monte Carlo sampling, and (5) the reconstruction of response functions from an imaginary-time auto-correlation function using MaxEnt techniques. Results obtained using schematic interactions, which have no sign problem, are presented to demonstrate the feasibility of the method, while an extrapolation method for realistic Hamiltonians is presented. In addition, applications at finite temperature are outlined.

Monte Carlo learning/biasing experiment with intelligent random numbers

International Nuclear Information System (INIS)

Booth, T.E.

1985-01-01

A Monte Carlo learning and biasing technique is described that does its learning and biasing in the random number space rather than the physical phase-space. The technique is probably applicable to all linear Monte Carlo problems, but no proof is provided here. Instead, the technique is illustrated with a simple Monte Carlo transport problem. Problems encountered, problems solved, and speculations about future progress are discussed. 12 refs

Temperature variance study in Monte-Carlo photon transport theory

International Nuclear Information System (INIS)

Giorla, J.

1985-10-01

We study different Monte-Carlo methods for solving radiative transfer problems, and particularly Fleck's Monte-Carlo method. We first give the different time-discretization schemes and the corresponding stability criteria. Then we write the temperature variance as a function of the variances of temperature and absorbed energy at the previous time step. Finally we obtain some stability criteria for the Monte-Carlo method in the stationary case [fr

Monte Carlo applications to radiation shielding problems

International Nuclear Information System (INIS)

Subbaiah, K.V.

2009-01-01

Monte Carlo methods are a class of computational algorithms that rely on repeated random sampling of physical and mathematical systems to compute their results. However, basic concepts of MC are both simple and straightforward and can be learned by using a personal computer. Uses of Monte Carlo methods require large amounts of random numbers, and it was their use that spurred the development of pseudorandom number generators, which were far quicker to use than the tables of random numbers which had been previously used for statistical sampling. In Monte Carlo simulation of radiation transport, the history (track) of a particle is viewed as a random sequence of free flights that end with an interaction event where the particle changes its direction of movement, loses energy and, occasionally, produces secondary particles. The Monte Carlo simulation of a given experimental arrangement (e.g., an electron beam, coming from an accelerator and impinging on a water phantom) consists of the numerical generation of random histories. To simulate these histories we need an interaction model, i.e., a set of differential cross sections (DCS) for the relevant interaction mechanisms. The DCSs determine the probability distribution functions (pdf) of the random variables that characterize a track; 1) free path between successive interaction events, 2) type of interaction taking place and 3) energy loss and angular deflection in a particular event (and initial state of emitted secondary particles, if any). Once these pdfs are known, random histories can be generated by using appropriate sampling methods. If the number of generated histories is large enough, quantitative information on the transport process may be obtained by simply averaging over the simulated histories. The Monte Carlo method yields the same information as the solution of the Boltzmann transport equation, with the same interaction model, but is easier to implement. In particular, the simulation of radiation

Forests

International Nuclear Information System (INIS)

Melin, J.

1997-01-01

Forests have the capacity to trap and retain radionuclides for a substantial period of time. The dynamic behaviour of nutrients, pollution and radionuclides in forests is complex. The rotation period of a forest stand in the Nordic countries is about 100 years, whilst the time for decomposition of organic material in a forest environment can be several hundred years. This means that any countermeasure applied in the forest environment must have an effect for several decades, or be reapplied continuously for long periods of time. To mitigate the detrimental effect of a contaminated forest environment on man, and to minimise the economic loss in trade of contaminated forest products, it is necessary to understand the mechanisms of transfer of radionuclides through the forest environment. It must also be stressed that any countermeasure applied in the forest environment must be evaluated with respect to long, as well as short term, negative effects, before any decision about remedial action is taken. Of the radionuclides studied in forests in the past, radiocaesium has been the main contributor to dose to man. In this document, only radiocaesium will be discussed since data on the impact of other radionuclides on man are too scarce for a proper evaluation. (EG)

Forests and Forest Cover - MDC_NaturalForestCommunity

Data.gov (United States)

NSGIC Local Govt | GIS Inventory — A point feature class of NFCs - Natural Forest Communities. Natural Forest Community shall mean all stands of trees (including their associated understory) which...

Randomized quasi-Monte Carlo simulation of fast-ion thermalization

Science.gov (United States)

Höök, L. J.; Johnson, T.; Hellsten, T.

2012-01-01

This work investigates the applicability of the randomized quasi-Monte Carlo method for simulation of fast-ion thermalization processes in fusion plasmas, e.g. for simulation of neutral beam injection and radio frequency heating. In contrast to the standard Monte Carlo method, the quasi-Monte Carlo method uses deterministic numbers instead of pseudo-random numbers and has a statistical weak convergence close to {O}(N^{-1}) , where N is the number of markers. We have compared different quasi-Monte Carlo methods for a neutral beam injection scenario, which is solved by many realizations of the associated stochastic differential equation, discretized with the Euler-Maruyama scheme. The statistical convergence of the methods is measured for time steps up to 214.

A Monte Carlo algorithm for the Vavilov distribution

International Nuclear Information System (INIS)

Yi, Chul-Young; Han, Hyon-Soo

1999-01-01

Using the convolution property of the inverse Laplace transform, an improved Monte Carlo algorithm for the Vavilov energy-loss straggling distribution of the charged particle is developed, which is relatively simple and gives enough accuracy to be used for most Monte Carlo applications

Adaptive Multilevel Monte Carlo Simulation

KAUST Repository

Hoel, H

2011-08-23

This work generalizes a multilevel forward Euler Monte Carlo method introduced in Michael B. Giles. (Michael Giles. Oper. Res. 56(3):607–617, 2008.) for the approximation of expected values depending on the solution to an Itô stochastic differential equation. The work (Michael Giles. Oper. Res. 56(3):607– 617, 2008.) proposed and analyzed a forward Euler multilevelMonte Carlo method based on a hierarchy of uniform time discretizations and control variates to reduce the computational effort required by a standard, single level, Forward Euler Monte Carlo method. This work introduces an adaptive hierarchy of non uniform time discretizations, generated by an adaptive algorithmintroduced in (AnnaDzougoutov et al. Raùl Tempone. Adaptive Monte Carlo algorithms for stopped diffusion. In Multiscale methods in science and engineering, volume 44 of Lect. Notes Comput. Sci. Eng., pages 59–88. Springer, Berlin, 2005; Kyoung-Sook Moon et al. Stoch. Anal. Appl. 23(3):511–558, 2005; Kyoung-Sook Moon et al. An adaptive algorithm for ordinary, stochastic and partial differential equations. In Recent advances in adaptive computation, volume 383 of Contemp. Math., pages 325–343. Amer. Math. Soc., Providence, RI, 2005.). This form of the adaptive algorithm generates stochastic, path dependent, time steps and is based on a posteriori error expansions first developed in (Anders Szepessy et al. Comm. Pure Appl. Math. 54(10):1169– 1214, 2001). Our numerical results for a stopped diffusion problem, exhibit savings in the computational cost to achieve an accuracy of ϑ(TOL),from(TOL−3), from using a single level version of the adaptive algorithm to ϑ(((TOL−1)log(TOL))2).

USOS DE LA TIERRA, CONFLICTO Y PODER POLÍTICO EN LA FRONTERA AGRÍCOLA DEL ALTO SINÚ 1940-1970

Directory of Open Access Journals (Sweden)

Ruth Esther Gutierrez

2016-08-01

Full Text Available El presente artículo presenta un recorrido histórico, a partir del cual aborda las perspectivas historiográficas desde donde se ha estudiado la manera cómo la subregión del alto Sinú se convirtió en una frontera agrícola, que a mediados del siglo XX y más específicamente entre 1940 y 1970, despertó un profundo interés en las élites económicas y políticas nacionales y regionales, quienes abiertamente iniciaron un proceso de acaparamiento de las tierras e instauración de sistemas laborales que incluían la sujeción de la mano de obra a las corrientes bipartidistas liberales y conservadoras. Así mismo, se presentan las formas contestatarias que adoptaron los campesinos del alto Sinú frente al monopolio de la tierra ostentado por gamonales y funcionarios públicos.

Monte Carlo computation in the applied research of nuclear technology

International Nuclear Information System (INIS)

Xu Shuyan; Liu Baojie; Li Qin

2007-01-01

This article briefly introduces Monte Carlo Methods and their properties. It narrates the Monte Carlo methods with emphasis in their applications to several domains of nuclear technology. Monte Carlo simulation methods and several commonly used computer software to implement them are also introduced. The proposed methods are demonstrated by a real example. (authors)

Nested Sampling with Constrained Hamiltonian Monte Carlo

OpenAIRE

Betancourt, M. J.

2010-01-01

Nested sampling is a powerful approach to Bayesian inference ultimately limited by the computationally demanding task of sampling from a heavily constrained probability distribution. An effective algorithm in its own right, Hamiltonian Monte Carlo is readily adapted to efficiently sample from any smooth, constrained distribution. Utilizing this constrained Hamiltonian Monte Carlo, I introduce a general implementation of the nested sampling algorithm.

Combating Forest Corruption: the Forest Integrity Network

NARCIS (Netherlands)

Gupta, A.; Siebert, U.

2004-01-01

This article describes the strategies and activities of the Forest Integrity Network. One of the most important underlying causes of forest degradation is corruption and related illegal logging. The Forest Integrity Network is a timely new initiative to combat forest corruption. Its approach is to

Projeto : Plano de promoção para os roteiros municipais do Alto Alentejo

OpenAIRE

Dimas, Alicia

2012-01-01

Dissertação de Mestrado em Marketing e Promoção Turística apresentada à ESTM - Escola Superior de Turismo e Tecnologia do Mar do Instituto Politécnico de Leiria O presente projeto pretende delinear as linhas mestras de um plano de promoção para o produto Roteiros Municipais. O produto, cujo objetivo é auxiliar a atividade turística na região do Alto Alentejo (Portugal), foi colocado recentemente no mercado sem qualquer estratégia de marketing associada – este projeto pretende preencher ess...

Prevalencia de trastornos psicológicos en deportistas españoles de alto rendimiento

OpenAIRE

Marí Cortés, Josep; Pérez Recio, Guillermo

1997-01-01

El artículo compara la prevalencia de psicopatologías en una muestra de 1376 deportistas españoles de alto rendimiento con la incidencia de las mismas alteraciones psicológicas en la población general. Para explicar la prevalencia hallada en la muestra de deportistas, se sugieren y analizan tres hipótesis (la autoselección, el efecto psicológico beneficioso la práctica de actividad física y la generalización, a los conflictos personales, de las pautas de enfrentamiento que han resultado previ...

Statistics of Monte Carlo methods used in radiation transport calculation

International Nuclear Information System (INIS)

Datta, D.

2009-01-01

Radiation transport calculation can be carried out by using either deterministic or statistical methods. Radiation transport calculation based on statistical methods is basic theme of the Monte Carlo methods. The aim of this lecture is to describe the fundamental statistics required to build the foundations of Monte Carlo technique for radiation transport calculation. Lecture note is organized in the following way. Section (1) will describe the introduction of Basic Monte Carlo and its classification towards the respective field. Section (2) will describe the random sampling methods, a key component of Monte Carlo radiation transport calculation, Section (3) will provide the statistical uncertainty of Monte Carlo estimates, Section (4) will describe in brief the importance of variance reduction techniques while sampling particles such as photon, or neutron in the process of radiation transport

Multiple histogram method and static Monte Carlo sampling

NARCIS (Netherlands)

Inda, M.A.; Frenkel, D.

2004-01-01

We describe an approach to use multiple-histogram methods in combination with static, biased Monte Carlo simulations. To illustrate this, we computed the force-extension curve of an athermal polymer from multiple histograms constructed in a series of static Rosenbluth Monte Carlo simulations. From

Shell model Monte Carlo methods

International Nuclear Information System (INIS)

Koonin, S.E.

1996-01-01

We review quantum Monte Carlo methods for dealing with large shell model problems. These methods reduce the imaginary-time many-body evolution operator to a coherent superposition of one-body evolutions in fluctuating one-body fields; resultant path integral is evaluated stochastically. We first discuss the motivation, formalism, and implementation of such Shell Model Monte Carlo methods. There then follows a sampler of results and insights obtained from a number of applications. These include the ground state and thermal properties of pf-shell nuclei, thermal behavior of γ-soft nuclei, and calculation of double beta-decay matrix elements. Finally, prospects for further progress in such calculations are discussed. 87 refs

Discrete Diffusion Monte Carlo for Electron Thermal Transport

Science.gov (United States)

Chenhall, Jeffrey; Cao, Duc; Wollaeger, Ryan; Moses, Gregory

2014-10-01

The iSNB (implicit Schurtz Nicolai Busquet electron thermal transport method of Cao et al. is adapted to a Discrete Diffusion Monte Carlo (DDMC) solution method for eventual inclusion in a hybrid IMC-DDMC (Implicit Monte Carlo) method. The hybrid method will combine the efficiency of a diffusion method in short mean free path regions with the accuracy of a transport method in long mean free path regions. The Monte Carlo nature of the approach allows the algorithm to be massively parallelized. Work to date on the iSNB-DDMC method will be presented. This work was supported by Sandia National Laboratory - Albuquerque.

Monte Carlo techniques in diagnostic and therapeutic nuclear medicine

International Nuclear Information System (INIS)

Zaidi, H.

2002-01-01

Monte Carlo techniques have become one of the most popular tools in different areas of medical radiation physics following the development and subsequent implementation of powerful computing systems for clinical use. In particular, they have been extensively applied to simulate processes involving random behaviour and to quantify physical parameters that are difficult or even impossible to calculate analytically or to determine by experimental measurements. The use of the Monte Carlo method to simulate radiation transport turned out to be the most accurate means of predicting absorbed dose distributions and other quantities of interest in the radiation treatment of cancer patients using either external or radionuclide radiotherapy. The same trend has occurred for the estimation of the absorbed dose in diagnostic procedures using radionuclides. There is broad consensus in accepting that the earliest Monte Carlo calculations in medical radiation physics were made in the area of nuclear medicine, where the technique was used for dosimetry modelling and computations. Formalism and data based on Monte Carlo calculations, developed by the Medical Internal Radiation Dose (MIRD) committee of the Society of Nuclear Medicine, were published in a series of supplements to the Journal of Nuclear Medicine, the first one being released in 1968. Some of these pamphlets made extensive use of Monte Carlo calculations to derive specific absorbed fractions for electron and photon sources uniformly distributed in organs of mathematical phantoms. Interest in Monte Carlo-based dose calculations with β-emitters has been revived with the application of radiolabelled monoclonal antibodies to radioimmunotherapy. As a consequence of this generalized use, many questions are being raised primarily about the need and potential of Monte Carlo techniques, but also about how accurate it really is, what would it take to apply it clinically and make it available widely to the medical physics

Monte Carlo strategies in scientific computing

CERN Document Server

Liu, Jun S

2008-01-01

This paperback edition is a reprint of the 2001 Springer edition This book provides a self-contained and up-to-date treatment of the Monte Carlo method and develops a common framework under which various Monte Carlo techniques can be "standardized" and compared Given the interdisciplinary nature of the topics and a moderate prerequisite for the reader, this book should be of interest to a broad audience of quantitative researchers such as computational biologists, computer scientists, econometricians, engineers, probabilists, and statisticians It can also be used as the textbook for a graduate-level course on Monte Carlo methods Many problems discussed in the alter chapters can be potential thesis topics for masters’ or PhD students in statistics or computer science departments Jun Liu is Professor of Statistics at Harvard University, with a courtesy Professor appointment at Harvard Biostatistics Department Professor Liu was the recipient of the 2002 COPSS Presidents' Award, the most prestigious one for sta...

Off-diagonal expansion quantum Monte Carlo.

Science.gov (United States)

Albash, Tameem; Wagenbreth, Gene; Hen, Itay

2017-12-01

We propose a Monte Carlo algorithm designed to simulate quantum as well as classical systems at equilibrium, bridging the algorithmic gap between quantum and classical thermal simulation algorithms. The method is based on a decomposition of the quantum partition function that can be viewed as a series expansion about its classical part. We argue that the algorithm not only provides a theoretical advancement in the field of quantum Monte Carlo simulations, but is optimally suited to tackle quantum many-body systems that exhibit a range of behaviors from "fully quantum" to "fully classical," in contrast to many existing methods. We demonstrate the advantages, sometimes by orders of magnitude, of the technique by comparing it against existing state-of-the-art schemes such as path integral quantum Monte Carlo and stochastic series expansion. We also illustrate how our method allows for the unification of quantum and classical thermal parallel tempering techniques into a single algorithm and discuss its practical significance.

Dynamic bounds coupled with Monte Carlo simulations

Energy Technology Data Exchange (ETDEWEB)

Rajabalinejad, M., E-mail: M.Rajabalinejad@tudelft.n [Faculty of Civil Engineering, Delft University of Technology, Delft (Netherlands); Meester, L.E. [Delft Institute of Applied Mathematics, Delft University of Technology, Delft (Netherlands); Gelder, P.H.A.J.M. van; Vrijling, J.K. [Faculty of Civil Engineering, Delft University of Technology, Delft (Netherlands)

2011-02-15

For the reliability analysis of engineering structures a variety of methods is known, of which Monte Carlo (MC) simulation is widely considered to be among the most robust and most generally applicable. To reduce simulation cost of the MC method, variance reduction methods are applied. This paper describes a method to reduce the simulation cost even further, while retaining the accuracy of Monte Carlo, by taking into account widely present monotonicity. For models exhibiting monotonic (decreasing or increasing) behavior, dynamic bounds (DB) are defined, which in a coupled Monte Carlo simulation are updated dynamically, resulting in a failure probability estimate, as well as a strict (non-probabilistic) upper and lower bounds. Accurate results are obtained at a much lower cost than an equivalent ordinary Monte Carlo simulation. In a two-dimensional and a four-dimensional numerical example, the cost reduction factors are 130 and 9, respectively, where the relative error is smaller than 5%. At higher accuracy levels, this factor increases, though this effect is expected to be smaller with increasing dimension. To show the application of DB method to real world problems, it is applied to a complex finite element model of a flood wall in New Orleans.

Coded aperture optimization using Monte Carlo simulations

International Nuclear Information System (INIS)

Martineau, A.; Rocchisani, J.M.; Moretti, J.L.

2010-01-01

Coded apertures using Uniformly Redundant Arrays (URA) have been unsuccessfully evaluated for two-dimensional and three-dimensional imaging in Nuclear Medicine. The images reconstructed from coded projections contain artifacts and suffer from poor spatial resolution in the longitudinal direction. We introduce a Maximum-Likelihood Expectation-Maximization (MLEM) algorithm for three-dimensional coded aperture imaging which uses a projection matrix calculated by Monte Carlo simulations. The aim of the algorithm is to reduce artifacts and improve the three-dimensional spatial resolution in the reconstructed images. Firstly, we present the validation of GATE (Geant4 Application for Emission Tomography) for Monte Carlo simulations of a coded mask installed on a clinical gamma camera. The coded mask modelling was validated by comparison between experimental and simulated data in terms of energy spectra, sensitivity and spatial resolution. In the second part of the study, we use the validated model to calculate the projection matrix with Monte Carlo simulations. A three-dimensional thyroid phantom study was performed to compare the performance of the three-dimensional MLEM reconstruction with conventional correlation method. The results indicate that the artifacts are reduced and three-dimensional spatial resolution is improved with the Monte Carlo-based MLEM reconstruction.

Recovery of Forest Canopy Parameters by Inversion of Multispectral LiDAR Data

Directory of Open Access Journals (Sweden)

Andrew Wallace

2012-02-01

Full Text Available We describe the use of Bayesian inference techniques, notably Markov chain Monte Carlo (MCMC and reversible jump MCMC (RJMCMC methods, to recover forest structural and biochemical parameters from multispectral LiDAR (Light Detection and Ranging data. We use a variable dimension, multi-layered model to represent a forest canopy or tree, and discuss the recovery of structure and depth profiles that relate to photochemical properties. We first demonstrate how simple vegetation indices such as the Normalized Differential Vegetation Index (NDVI, which relates to canopy biomass and light absorption, and Photochemical Reflectance Index (PRI which is a measure of vegetation light use efficiency, can be measured from multispectral data. We further describe and demonstrate our layered approach on single wavelength real data, and on simulated multispectral data derived from real, rather than simulated, data sets. This evaluation shows successful recovery of a subset of parameters, as the complete recovery problem is ill-posed with the available data. We conclude that the approach has promise, and suggest future developments to address the current difficulties in parameter inversion.

Randomized quasi-Monte Carlo simulation of fast-ion thermalization

International Nuclear Information System (INIS)

Höök, L J; Johnson, T; Hellsten, T

2012-01-01

This work investigates the applicability of the randomized quasi-Monte Carlo method for simulation of fast-ion thermalization processes in fusion plasmas, e.g. for simulation of neutral beam injection and radio frequency heating. In contrast to the standard Monte Carlo method, the quasi-Monte Carlo method uses deterministic numbers instead of pseudo-random numbers and has a statistical weak convergence close to O(N -1 ), where N is the number of markers. We have compared different quasi-Monte Carlo methods for a neutral beam injection scenario, which is solved by many realizations of the associated stochastic differential equation, discretized with the Euler-Maruyama scheme. The statistical convergence of the methods is measured for time steps up to 2 14 . (paper)

Usefulness of the Monte Carlo method in reliability calculations

International Nuclear Information System (INIS)

Lanore, J.M.; Kalli, H.

1977-01-01

Three examples of reliability Monte Carlo programs developed in the LEP (Laboratory for Radiation Shielding Studies in the Nuclear Research Center at Saclay) are presented. First, an uncertainty analysis is given for a simplified spray system; a Monte Carlo program PATREC-MC has been written to solve the problem with the system components given in the fault tree representation. The second program MONARC 2 has been written to solve the problem of complex systems reliability by the Monte Carlo simulation, here again the system (a residual heat removal system) is in the fault tree representation. Third, the Monte Carlo program MONARC was used instead of the Markov diagram to solve the simulation problem of an electric power supply including two nets and two stand-by diesels

Combinatorial nuclear level density by a Monte Carlo method

International Nuclear Information System (INIS)

Cerf, N.

1994-01-01

We present a new combinatorial method for the calculation of the nuclear level density. It is based on a Monte Carlo technique, in order to avoid a direct counting procedure which is generally impracticable for high-A nuclei. The Monte Carlo simulation, making use of the Metropolis sampling scheme, allows a computationally fast estimate of the level density for many fermion systems in large shell model spaces. We emphasize the advantages of this Monte Carlo approach, particularly concerning the prediction of the spin and parity distributions of the excited states,and compare our results with those derived from a traditional combinatorial or a statistical method. Such a Monte Carlo technique seems very promising to determine accurate level densities in a large energy range for nuclear reaction calculations

Monte Carlo variance reduction approaches for non-Boltzmann tallies

International Nuclear Information System (INIS)

Booth, T.E.

1992-12-01

Quantities that depend on the collective effects of groups of particles cannot be obtained from the standard Boltzmann transport equation. Monte Carlo estimates of these quantities are called non-Boltzmann tallies and have become increasingly important recently. Standard Monte Carlo variance reduction techniques were designed for tallies based on individual particles rather than groups of particles. Experience with non-Boltzmann tallies and analog Monte Carlo has demonstrated the severe limitations of analog Monte Carlo for many non-Boltzmann tallies. In fact, many calculations absolutely require variance reduction methods to achieve practical computation times. Three different approaches to variance reduction for non-Boltzmann tallies are described and shown to be unbiased. The advantages and disadvantages of each of the approaches are discussed

Variational Monte Carlo Technique

Indian Academy of Sciences (India)

Home; Journals; Resonance – Journal of Science Education; Volume 19; Issue 8. Variational Monte Carlo Technique: Ground State Energies of Quantum Mechanical Systems. Sukanta Deb. General Article Volume 19 Issue 8 August 2014 pp 713-739 ...

La propiedad rústica en la provincia de Jaén: el solar forestal público

Directory of Open Access Journals (Sweden)

José Domingo Sánchez Martínez

2000-01-01

Full Text Available Se explica el origen, la evolución y la tipología de los montes sujetos a propiedad y/o gestión pública en la provincia de Jaén, diferenciando tres grandes ámbitos forestales que se desarrollan en las periferias montañosas meridional, oriental y septentrional. Se relaciona la sucesión de diferentes fases de la política forestal española del último siglo y medio con el tratamiento concedido en cada momento a las cuestiones patrimoniales.

A model of Alto Lazio boiling water reactor using the LEGO code balance of plant simulation

International Nuclear Information System (INIS)

Spelta, S.; Garbossa, G.B.

1989-01-01

An extensive effort has been made at the Italian National Electricity Board (ENEL) to construct and validate a LEGO model capable of simulating the operational transients of the Alto Lazio Nuclear Station, a two twin-units site with BWR/6 class reactors, rated at 2894 MWt and with Mark III containment. LEGO is a modular package developed at the Research and development Department of the Italian National Electricity Board (CRA-ENEL) for computer aided modeling of fossil-fired and nuclear steam power plants. In this paper a system analysis model capable of describing steady-state and transient performance of the Balance of Plant (BOP) of the Alto Lazio Power Station is presented. This is one of two companion papers devoted to the description of the overall plant model including both the Nuclear Steam Supply System (NSSS) and the BOP. In the paper, after a brief summary of the main LEGO characteristics, a description of the BOP lay-out is presented. The overall model, which has been set-up, including control systems and automation, is very detailed and consists of almost 2000 differential or algebraic equations. After a brief description of the mathematical model, two significant transients obtained using the overall model are presented and discussed

Índice hiperbárico en la prevención primaria de las complicaciones hipertensivas del embarazo de alto riesgo

Directory of Open Access Journals (Sweden)

Alfonso Otero González

2015-11-01

Resultados: En las gestantes remitidas a consulta de Nefrología (38,2%, diagnosticadas de alto riesgo de PE y tratadas con AAS 100 mg nocturno (desde la semana 17 se redujo la incidencia de episodios de PE un 96,94%.

Discrete diffusion Monte Carlo for frequency-dependent radiative transfer

International Nuclear Information System (INIS)

Densmore, Jeffery D.; Thompson, Kelly G.; Urbatsch, Todd J.

2011-01-01

Discrete Diffusion Monte Carlo (DDMC) is a technique for increasing the efficiency of Implicit Monte Carlo radiative-transfer simulations. In this paper, we develop an extension of DDMC for frequency-dependent radiative transfer. We base our new DDMC method on a frequency integrated diffusion equation for frequencies below a specified threshold. Above this threshold we employ standard Monte Carlo. With a frequency-dependent test problem, we confirm the increased efficiency of our new DDMC technique. (author)

Variations in dung beetles assemblages (Coleoptera: Scarabaeidae within two rain forest habitats in French Guiana Variaciones de las agregaciones de escarabajos coprófagos (Coleoptera: Scarabaeidae en dos hábitats del bosque tropical en la Guayana Francesa

Directory of Open Access Journals (Sweden)

François Feer

2013-06-01

Full Text Available The structure of dung beetle communities inhabiting tropical forests are known to be sensitive to many kinds of environmental changes such as microclimate related to vegetation structure. I examined Scarabaeinae assemblages in two sites of undisturbed high forest and two sites of low forest forming a transitional zone with the open habitat of an inselberg in French Guiana. Sampling was made with pitfall and flight interception traps during 2003 and 2004. The driest and warmest conditions characterized the low forest sites. Across two years we obtained 2 927 individuals from 61 species with pitfall traps and 1 431 individuals from 85 species with flight interception traps. Greater species richness and abundance characterized all sites sampled with pitfall traps during 2003 more than 2004. In 2003 no differences were detected among sites by rarefaction analyses. In 2004 the species richest high forest site was significantly different from one of the low forest sites. For both years Clench model asymptotes for species richness were greater in high forest than in low forest sites. For both years, mean per-trap species richness, abundance and biomass among high forest sites were similar and higher than in low forest sites, especially where the lowest humidity and the highest temperature were recorded. Within the two low forest sites, species richness and abundance recorded during the second year, decreased with distance to edge. Different dominant roller species characterized the pitfall samples in one site of low forest and in other sites. Small variations in microclimatic conditions correlated to canopy height and openness likely affected dung beetle assemblages but soil depth and the presence of large mammals providing dung resource may also play a significant role.Es conocido que la estructura de las comunidades de escarabajos coprófagos que habitan los bosques tropicales es sensible a diferentes tipos de cambios ambientales tales como el

Modified Monte Carlo procedure for particle transport problems

International Nuclear Information System (INIS)

Matthes, W.

1978-01-01

The simulation of photon transport in the atmosphere with the Monte Carlo method forms part of the EURASEP-programme. The specifications for the problems posed for a solution were such, that the direct application of the analogue Monte Carlo method was not feasible. For this reason the standard Monte Carlo procedure was modified in the sense that additional properly weighted branchings at each collision and transport process in a photon history were introduced. This modified Monte Carlo procedure leads to a clear and logical separation of the essential parts of a problem and offers a large flexibility for variance reducing techniques. More complex problems, as foreseen in the EURASEP-programme (e.g. clouds in the atmosphere, rough ocean-surface and chlorophyl-distribution in the ocean) can be handled by recoding some subroutines. This collision- and transport-splitting procedure can of course be performed differently in different space- and energy regions. It is applied here only for a homogeneous problem

Uncertainty analysis in Monte Carlo criticality computations

International Nuclear Information System (INIS)

Qi Ao

2011-01-01

Highlights: ► Two types of uncertainty methods for k eff Monte Carlo computations are examined. ► Sampling method has the least restrictions on perturbation but computing resources. ► Analytical method is limited to small perturbation on material properties. ► Practicality relies on efficiency, multiparameter applicability and data availability. - Abstract: Uncertainty analysis is imperative for nuclear criticality risk assessments when using Monte Carlo neutron transport methods to predict the effective neutron multiplication factor (k eff ) for fissionable material systems. For the validation of Monte Carlo codes for criticality computations against benchmark experiments, code accuracy and precision are measured by both the computational bias and uncertainty in the bias. The uncertainty in the bias accounts for known or quantified experimental, computational and model uncertainties. For the application of Monte Carlo codes for criticality analysis of fissionable material systems, an administrative margin of subcriticality must be imposed to provide additional assurance of subcriticality for any unknown or unquantified uncertainties. Because of a substantial impact of the administrative margin of subcriticality on economics and safety of nuclear fuel cycle operations, recently increasing interests in reducing the administrative margin of subcriticality make the uncertainty analysis in criticality safety computations more risk-significant. This paper provides an overview of two most popular k eff uncertainty analysis methods for Monte Carlo criticality computations: (1) sampling-based methods, and (2) analytical methods. Examples are given to demonstrate their usage in the k eff uncertainty analysis due to uncertainties in both neutronic and non-neutronic parameters of fissionable material systems.

Forest resources of the Nez Perce National Forest

Science.gov (United States)

Michele Disney

2010-01-01

As part of a National Forest System cooperative inventory, the Interior West Forest Inventory and Analysis (IWFIA) Program of the USDA Forest Service conducted a forest resource inventory on the Nez Perce National Forest using a nationally standardized mapped-plot design (for more details see the section "Inventory methods"). This report presents highlights...

TNO Photometry and Spectroscopy at ESO and Calar Alto

Science.gov (United States)

Boehnhardt, H.; Sekiguchi, T.; Vair, M.; Hainaut, O.; Delahodde, C.; West, R. M.; Tozzi, G. P.; Barrera, L.; Birkle, K.; Watanabe, J.; Meech, K.

New photometry and spectroscopy of Transneptunian objects (TNO) has been obtained at ESO (VLT+FORS1, NTT+SOFI) and the Calar Alto (3.5m+MOSCA) observatory. BVRI photometry of more than 10 objects confirms the general colour-colour distribution of TNOs found previously. Quasi-simultaneous spectroscopy in the visible wavelength range of 5 TNOs did not reveal any spectral signature apart from the spetral gradients which are in agreement with the broadband colours. JHK filter photometry of 3 objects indicates that the reddening may only occur in the near-IR at least in some cases. Using new observations from the ESO VLT the lightcurve, colours and spectrum of 1996TO66 are investigated: the rotation period of 6.25h is confirmed, also the change in the lightcurve between 1997 and 1998 which indicates an exceptional behaviour in this object (temporary cometary activity ?). The 1999 photometry and spectroscopy in the visible revealed solar colours, no reddening and no spectral features. V-R colour changes over the rotation phase are not found. This works is done in colaboration with:

HTLV-I en población de alto riesgo sexual de Pisco, Ica, Perú.

OpenAIRE

GARRIDO, Patricia; ANICAMA, Rolando; GOTUZZO, Eduardo; CHAUCA, Gloria; WATTS, Douglas

2013-01-01

Objetivo: Se estudiaron 141 personas con alto riesgo sexual en la ciudad de Pisco para detectar infección por HTLV-I. Material y Métodos: Se encuestaron y se tomaron muestras de sangre a 141 personas que involucró a trabajadoras sexuales (32), varones homosexuales (54), y varones bisexuales(55). Resultados: Tres de treintidós (10.4%) trabajadoras sexuales fueron positivas; uno de cincuenticuatro (1.9%) de varones homosexuales y ninguno de 55 bisexuales. Hubo una elevada frecuencia de parejas,...

Trastorno de espectro autista de alto funcionamiento y su relación con la familia

OpenAIRE

González Villafáñez, María Felisa

2014-01-01

El objetivo de este trabajo consiste en analizar las dificultades que tienen los adolescentes con Trastorno del Espectro Autista de Alto Funcionamiento (TEA-AF), en sus relaciones sociales y familiares. Este análisis ha requerido la revisión de las diversas definiciones y teorías explicativas del citado trastorno. A su vez también se describen las terapias conocidas hasta hoy para crear un programa de apoyo a los familiares (padres, hermanos, abuelos…) que se relacionan con dicho adolescente....

IRMANDADE DO ROSÁRIO DOS PRETOS DE SANTA EFIGÊNIA DO ALTO DA CRUZ: PROPOSTA DE GESTÃO DE UM ARQUIVO EM OURO PRETO (MG).

OpenAIRE

Silva, Laurimar Gomes da

2009-01-01

A proposta deste projeto visa atender a uma demanda da Irmandade do Rosário dos Pretos de Santa Efigênia do Alto da Cruz na cidade de Ouro Preto (MG), a fim de consolidar parcerias no sentido de preservar e gerir o arquivo. Fundada em 1719, compondo uma das vinte irmandades encontradas na cidade, a Irmandade do Rosário dos Pretos de Santa Efigênia do Alto da Cruz se particulariza por ser a responsável pela administração da igreja barroca de Santa Efigênia, que no imaginário ...

An Overview of the Monte Carlo Application ToolKit (MCATK)

Energy Technology Data Exchange (ETDEWEB)

Trahan, Travis John [Los Alamos National Lab. (LANL), Los Alamos, NM (United States)

2016-01-07

MCATK is a C++ component-based Monte Carlo neutron-gamma transport software library designed to build specialized applications and designed to provide new functionality in existing general-purpose Monte Carlo codes like MCNP; it was developed with Agile software engineering methodologies under the motivation to reduce costs. The characteristics of MCATK can be summarized as follows: MCATK physics – continuous energy neutron-gamma transport with multi-temperature treatment, static eigenvalue (k and α) algorithms, time-dependent algorithm, fission chain algorithms; MCATK geometry – mesh geometries, solid body geometries. MCATK provides verified, unit-tested Monte Carlo components, flexibility in Monte Carlo applications development, and numerous tools such as geometry and cross section plotters. Recent work has involved deterministic and Monte Carlo analysis of stochastic systems. Static and dynamic analysis is discussed, and the results of a dynamic test problem are given.

General and biochemical properties of forest soils on serpentinized areas of Galicia; Propiedades generales y bioquimicas de suelos forestales en areas serpentinizadas de Galicia

Energy Technology Data Exchange (ETDEWEB)

Pena, Wagner; Briceno, Jorge A., E-mail: wpenac@cia.ucr.ac.c [Univ. de Costa Rica, Centro de Investigaciones Agronomicas (Costa Rica); Leiros de la Pena, Maria del Carmen [Univ. de Santiago de Compostela, Dept. Edafologia y Quimica Agricola (Spain)

2005-07-15

The use agricultural and forest soils is normal in the territory of the Autonomous Community of Galicia for submitting the physical, chemical and agroclimatic optimum. The region has areas by the characteristics of the environment has been impossible to use agricultural and forestry. Serpentinite soils developed on frequently drew attention to the high degree of infertility, caused by nutrient deficiency and excess toxic metal content. The results of general physical and chemical properties have been analyzed of the total metal content and biochemical parameters of 28 soils on serpentinite in A Coruna Galicia, Spain. The infertility of soil organic matter, nutrients and low microbial activity and biochemistry have been caused by high level of pH and nutritional imbalance due to the high content of Mg{sup 2+}. OM (organic matter) content, total N, soluble P and exchangeable K{sup +}, were lower due to the underdevelopment of vegetation and the potential toxic effects that have metal on the biological activity in the middle. (author) [Spanish] El uso agricola y forestal de los suelos es normal en el territorio de la Comunidad Autonoma de Galicia por haber presentado las condiciones fisicas, quimicas y agroclimaticas optimas. La region tiene zonas que por las caracteristicas del medio ha sido imposible el uso agricola y forestal. Los suelos desarrollados sobre serpentinita es frecuente que hayan llamado la atencion por el elevado grado de infertilidad, provocado por la deficiencia de nutrimentos y por el exceso contenido en metales toxicos. Los resultados de las propiedades generales fisicas y quimicas han sido analizados del contenido total de metales y de parametros bioquimicos de 28 suelos sobre serpentinita en A Coruna, Galicia, Espana. La infertilidad de los suelos de materia organica, nutrimentos y baja actividad microbiana y bioquimica han sido causadas por el alto nivel de pH y el desbalance nutricional por el alto contenido de Mg{sup 2+}. El contenido de MO

Biodiversity and distribution of arbuscular mycorrhizal fungi in Araucaria angustifolia forest Biodiversidade e distribuição de fungos micorrízicos arbusculares em floresta de Araucaria angustifolia

Directory of Open Access Journals (Sweden)

Milene Moreira

2007-08-01

Full Text Available Araucaria angustifolia (Bert. O. Ktze. is an endangered Brazilian coniferous tree that has been almost exterminated in the native areas because of uncontrolled wood exploitation. This tree has been shown to be highly dependent on arbuscular mycorrhizal fungi (AMF and, therefore, AMF may be essential for forest sustainability and biological diversity. Root colonization, density and diversity of AMF spores were assessed in two Araucaria forest stands at the State Park of Alto Ribeira (PETAR, at two sampling dates: May and October. A comparison was made between a mature native stand composed of Araucaria trees mixed into a variety of tropical trees and shrubs, without any sign of anthropogenic interference (FN and an Araucaria stand planted in 1987 (R, which has been used as a pasture. Assessments included percent root colonization, AMF spore numbers and species richness, Simpson's dominance index (Is, and Shannon's diversity index (H. Mycorrhizal root colonization did not differ between ecosystems in May. In October, however, the native stand (FN presented a higher colonization than the planted forest (R, and the root colonization was more intense than in May. When considering both sampling periods and forests, 27 species of AM fungi, with higher numbers of spores in FN than in R were found. Canonical discriminant analysis (CDA indicated Shannon's diversity index as the ecological attribute that contributed the most to distinguish between forest ecosystems, with higher value of H in FN in relation to R. CDA showed to be a useful tool for the study of ecological attributes.O pinheiro brasileiro, Araucaria angustifolia (Bert. O. Ktze é uma espécie arbórea ameaçada de extinção, quase exterminada pela exploração descontrolada de madeira. Essa espécie demonstra alta dependência de fungos micorrízicos arbusculares (FMA, de grande importância para a estabilidade do ecossistema e manutenção da biodiversidade. Avaliou-se a coloniza

Efficiency and accuracy of Monte Carlo (importance) sampling

NARCIS (Netherlands)

Waarts, P.H.

2003-01-01

Monte Carlo Analysis is often regarded as the most simple and accurate reliability method. Be-sides it is the most transparent method. The only problem is the accuracy in correlation with the efficiency. Monte Carlo gets less efficient or less accurate when very low probabilities are to be computed

Suppression of the initial transient in Monte Carlo criticality simulations; Suppression du regime transitoire initial des simulations Monte-Carlo de criticite

Energy Technology Data Exchange (ETDEWEB)

Richet, Y

2006-12-15

Criticality Monte Carlo calculations aim at estimating the effective multiplication factor (k-effective) for a fissile system through iterations simulating neutrons propagation (making a Markov chain). Arbitrary initialization of the neutron population can deeply bias the k-effective estimation, defined as the mean of the k-effective computed at each iteration. A simplified model of this cycle k-effective sequence is built, based on characteristics of industrial criticality Monte Carlo calculations. Statistical tests, inspired by Brownian bridge properties, are designed to discriminate stationarity of the cycle k-effective sequence. The initial detected transient is, then, suppressed in order to improve the estimation of the system k-effective. The different versions of this methodology are detailed and compared, firstly on a plan of numerical tests fitted on criticality Monte Carlo calculations, and, secondly on real criticality calculations. Eventually, the best methodologies observed in these tests are selected and allow to improve industrial Monte Carlo criticality calculations. (author)

Coal exploration in the Alto San Jorge area, Cordoba Department. Exploracion de carbones en el Ato San Jorge, Departamento de Cordoba

Energy Technology Data Exchange (ETDEWEB)

Ospina, L H; Oquendo, G G [Geominas Ltda, Medellin (Colombia)

1989-01-01

A Mining Feasibility Study in the Area of Alto San Jorge, Department of Cordoba, Colombia, was commissioned by CARBOCOL S.A. to the Consortium Geominas-NACI. An area of 800 Ka2 was explored to define surface mining possibilities within two subareas referred to as Alto San Jorge and San Pedro Ure. Rocks of Cretaceous, Tertiary and Quaternary age crop out in the zone. In the subarea Alto San Jorge the principal structure is a syncline with a south-north direction. The San Pedro Ure subarea is formed by undulations with flanks of low dip, the most important being the San Antonio Syncline because it contains the mining block. The geological study of the surface demonstrated the existence of coal in the Oligocene Cienaga de Oro Formation and the Niocene Cerrito Formation, with potential resources of 6.3 billion tons. The subsequent exploration of the subsoil, with 20.618 m of drilling, permitted determination of demonstrated reserves in the order of 2.9 billion tons within two areas. In the sector selected for the mine plan, in the area of San Pedro-Puerto Libertador, 7.791 m of drilling was accomplished to define a demonstrated reserve of 515 million tons of coal down to a depth of 200. The combustible type coal has 5.000 cal/g. Complete mining schedules were developed at the prefeasibility level for two surface mines with productions of 1.5 MMTY and 4 MMTY. 9 figs., 3 tabs., 28 refs.

Monte Carlo criticality analysis for dissolvers with neutron poison

International Nuclear Information System (INIS)

Yu, Deshun; Dong, Xiufang; Pu, Fuxiang.

1987-01-01

Criticality analysis for dissolvers with neutron poison is given on the basis of Monte Carlo method. In Monte Carlo calculations of thermal neutron group parameters for fuel pieces, neutron transport length is determined in terms of maximum cross section approach. A set of related effective multiplication factors (K eff ) are calculated by Monte Carlo method for the three cases. Related numerical results are quite useful for the design and operation of this kind of dissolver in the criticality safety analysis. (author)

Fifty-year spatiotemporal analysis of landscape changes in the Mont Saint-Hilaire UNESCO Biosphere Reserve (Quebec, Canada).

Science.gov (United States)

Béliveau, Marc; Germain, Daniel; Ianăş, Ana-Neli

2017-05-01

Diachronic analysis with a GIS-based classification of land-use changes based on aerial photographs, orthophotos, topographic maps, geotechnical reports, urban plans, and using landscape metrics has permitted insight into the driving forces responsible for landscape fragmentation in the Mont Saint-Hilaire (MSH) Biosphere Reserve over the period 1958-2015. Although the occurrence of exogenous factors, such as extreme weather and fires, can have a significant influence on the fragmentation of the territory in time and space, the accelerated development of the built environment (+470%) is nevertheless found to be primarily responsible for landscape fragmentation and the loss of areas formerly occupied by orchards, agriculture, and woodlands. The landscape metrics used corroborate these results, with a simplification of the shape of polygons, and once again reveal the difficulties of harmonizing different land uses. MSH has become somewhat of a forest island in a sea of residential development and agriculture. To counter this isolation of fragmented habitat components, forest corridors have been proposed and developed for the Biosphere Reserve and particularly for the core area. Two corridors, to the north and south, are used to connect the protected area and other wooded areas at the regional scale, in order to promote genetic exchange between populations of various species. In that regard, the forest buffer zone around the hill continues to play a key role and has great ecological value for species and ecological preservation and conservation. However, appropriate management and landscape preservation actions should recognize and focus on landscape composition and the associated geographical configuration.

Improvements for Monte Carlo burnup calculation

Energy Technology Data Exchange (ETDEWEB)

Shenglong, Q.; Dong, Y.; Danrong, S.; Wei, L., E-mail: qiangshenglong@tsinghua.org.cn, E-mail: d.yao@npic.ac.cn, E-mail: songdr@npic.ac.cn, E-mail: luwei@npic.ac.cn [Nuclear Power Inst. of China, Cheng Du, Si Chuan (China)

2015-07-01

Monte Carlo burnup calculation is development trend of reactor physics, there would be a lot of work to be done for engineering applications. Based on Monte Carlo burnup code MOI, non-fuel burnup calculation methods and critical search suggestions will be mentioned in this paper. For non-fuel burnup, mixed burnup mode will improve the accuracy of burnup calculation and efficiency. For critical search of control rod position, a new method called ABN based on ABA which used by MC21 will be proposed for the first time in this paper. (author)

Monte Carlo dose distributions for radiosurgery

International Nuclear Information System (INIS)

Perucha, M.; Leal, A.; Rincon, M.; Carrasco, E.

2001-01-01

The precision of Radiosurgery Treatment planning systems is limited by the approximations of their algorithms and by their dosimetrical input data. This fact is especially important in small fields. However, the Monte Carlo methods is an accurate alternative as it considers every aspect of particle transport. In this work an acoustic neurinoma is studied by comparing the dose distribution of both a planning system and Monte Carlo. Relative shifts have been measured and furthermore, Dose-Volume Histograms have been calculated for target and adjacent organs at risk. (orig.)

Quasi Monte Carlo methods for optimization models of the energy industry with pricing and load processes; Quasi-Monte Carlo Methoden fuer Optimierungsmodelle der Energiewirtschaft mit Preis- und Last-Prozessen

Energy Technology Data Exchange (ETDEWEB)

Leoevey, H.; Roemisch, W. [Humboldt-Univ., Berlin (Germany)

2015-07-01

We discuss progress in quasi Monte Carlo methods for numerical calculation integrals or expected values and justify why these methods are more efficient than the classic Monte Carlo methods. Quasi Monte Carlo methods are found to be particularly efficient if the integrands have a low effective dimension. That's why We also discuss the concept of effective dimension and prove on the example of a stochastic Optimization model of the energy industry that such models can posses a low effective dimension. Modern quasi Monte Carlo methods are therefore for such models very promising. [German] Wir diskutieren Fortschritte bei Quasi-Monte Carlo Methoden zur numerischen Berechnung von Integralen bzw. Erwartungswerten und begruenden warum diese Methoden effizienter sind als die klassischen Monte Carlo Methoden. Quasi-Monte Carlo Methoden erweisen sich als besonders effizient, falls die Integranden eine geringe effektive Dimension besitzen. Deshalb diskutieren wir auch den Begriff effektive Dimension und weisen am Beispiel eines stochastischen Optimierungsmodell aus der Energiewirtschaft nach, dass solche Modelle eine niedrige effektive Dimension besitzen koennen. Moderne Quasi-Monte Carlo Methoden sind deshalb fuer solche Modelle sehr erfolgversprechend.

Shell model Monte Carlo methods

International Nuclear Information System (INIS)

Koonin, S.E.; Dean, D.J.; Langanke, K.

1997-01-01

We review quantum Monte Carlo methods for dealing with large shell model problems. These methods reduce the imaginary-time many-body evolution operator to a coherent superposition of one-body evolutions in fluctuating one-body fields; the resultant path integral is evaluated stochastically. We first discuss the motivation, formalism, and implementation of such Shell Model Monte Carlo (SMMC) methods. There then follows a sampler of results and insights obtained from a number of applications. These include the ground state and thermal properties of pf-shell nuclei, the thermal and rotational behavior of rare-earth and γ-soft nuclei, and the calculation of double beta-decay matrix elements. Finally, prospects for further progress in such calculations are discussed. (orig.)

Monte Carlo Methods in ICF

Science.gov (United States)

Zimmerman, George B.

Monte Carlo methods appropriate to simulate the transport of x-rays, neutrons, ions and electrons in Inertial Confinement Fusion targets are described and analyzed. The Implicit Monte Carlo method of x-ray transport handles symmetry within indirect drive ICF hohlraums well, but can be improved 50X in efficiency by angular biasing the x-rays towards the fuel capsule. Accurate simulation of thermonuclear burn and burn diagnostics involves detailed particle source spectra, charged particle ranges, inflight reaction kinematics, corrections for bulk and thermal Doppler effects and variance reduction to obtain adequate statistics for rare events. It is found that the effects of angular Coulomb scattering must be included in models of charged particle transport through heterogeneous materials.

Monte Carlo methods in ICF

International Nuclear Information System (INIS)

Zimmerman, George B.

1997-01-01

Monte Carlo methods appropriate to simulate the transport of x-rays, neutrons, ions and electrons in Inertial Confinement Fusion targets are described and analyzed. The Implicit Monte Carlo method of x-ray transport handles symmetry within indirect drive ICF hohlraums well, but can be improved 50X in efficiency by angular biasing the x-rays towards the fuel capsule. Accurate simulation of thermonuclear burn and burn diagnostics involves detailed particle source spectra, charged particle ranges, inflight reaction kinematics, corrections for bulk and thermal Doppler effects and variance reduction to obtain adequate statistics for rare events. It is found that the effects of angular Coulomb scattering must be included in models of charged particle transport through heterogeneous materials

Improved safety features in the design of Alto Lazio NPP

International Nuclear Information System (INIS)

Bava, G.; Cianciolo, T.; Del Nero, G.

1988-01-01

The ALTO LAZIO Nuclear Power Plant, two 1000Mwe units, is a BWR 6/MARK III located about 100 km north of Rome, on the Tyrrhenian Sea Coasts. The construction of the plant started in 1978, but it has recently been stopped by a Government decision following a national referendum, when the units were about 70% completed. This paper is mainly intended to illustrate the major safety features which have been implemented as result of specific requirements issued by the safety authority (ENEA DISP) during the construction permit stage or the subsequent licensing process. One of the tools used to identify the need for design modifications has been a comprehensive reliability analysis of safety system: in the paper the methods used and the major results obtained by this study are briefly presented. Also, the approach used in the investigation of severe accidents and major applications in the area of plant design and emergency procedures are briefly discussed; furthermore the trend toward a simpler mitigation concept is described

BREM5 electroweak Monte Carlo

International Nuclear Information System (INIS)

Kennedy, D.C. II.

1987-01-01

This is an update on the progress of the BREMMUS Monte Carlo simulator, particularly in its current incarnation, BREM5. The present report is intended only as a follow-up to the Mark II/Granlibakken proceedings, and those proceedings should be consulted for a complete description of the capabilities and goals of the BREMMUS program. The new BREM5 program improves on the previous version of BREMMUS, BREM2, in a number of important ways. In BREM2, the internal loop (oblique) corrections were not treated in consistent fashion, a deficiency that led to renormalization scheme-dependence; i.e., physical results, such as cross sections, were dependent on the method used to eliminate infinities from the theory. Of course, this problem cannot be tolerated in a Monte Carlo designed for experimental use. BREM5 incorporates a new way of treating the oblique corrections, as explained in the Granlibakken proceedings, that guarantees renormalization scheme-independence and dramatically simplifies the organization and calculation of radiative corrections. This technique is to be presented in full detail in a forthcoming paper. BREM5 is, at this point, the only Monte Carlo to contain the entire set of one-loop corrections to electroweak four-fermion processes and renormalization scheme-independence. 3 figures

Forest rights

DEFF Research Database (Denmark)

Balooni, Kulbhushan; Lund, Jens Friis

2014-01-01

One of the proposed strategies for implementation of reducing emissions from deforestation and forest degradation plus (REDD+) is to incentivize conservation of forests managed by communities under decentralized forest management. Yet, we argue that this is a challenging road to REDD+ because...... conservation of forests under existing decentralized management arrangements toward a push for extending the coverage of forests under decentralized management, making forest rights the hard currency of REDD+....

PEPSI: a Monte Carlo generator for polarized leptoproduction

International Nuclear Information System (INIS)

Mankiewicz, L.

1992-01-01

We describe PEPSI (Polarized Electron Proton Scattering Interactions) a Monte Carlo program for the polarized deep inelastic leptoproduction mediated by electromagnetic interaction. The code is a modification of the LEPTO 4.3 Lund Monte Carlo for unpolarized scattering and requires the standard polarization-independent JETSET routines to perform fragmentation into final hadrons. (orig.)

Importance estimation in Monte Carlo modelling of neutron and photon transport

International Nuclear Information System (INIS)

Mickael, M.W.

1992-01-01

The estimation of neutron and photon importance in a three-dimensional geometry is achieved using a coupled Monte Carlo and diffusion theory calculation. The parameters required for the solution of the multigroup adjoint diffusion equation are estimated from an analog Monte Carlo simulation of the system under investigation. The solution of the adjoint diffusion equation is then used as an estimate of the particle importance in the actual simulation. This approach provides an automated and efficient variance reduction method for Monte Carlo simulations. The technique has been successfully applied to Monte Carlo simulation of neutron and coupled neutron-photon transport in the nuclear well-logging field. The results show that the importance maps obtained in a few minutes of computer time using this technique are in good agreement with Monte Carlo generated importance maps that require prohibitive computing times. The application of this method to Monte Carlo modelling of the response of neutron porosity and pulsed neutron instruments has resulted in major reductions in computation time. (Author)

Forest structure in low diversity tropical forests: a study of Hawaiian wet and dry forests

Science.gov (United States)

R. Ostertag; F. Inman-Narahari; S. Cordell; C.P. Giardina; L. Sack

2014-01-01

The potential influence of diversity on ecosystem structure and function remains a topic of significant debate, especially for tropical forests where diversity can range widely. We used Center for Tropical Forest Science (CTFS) methodology to establish forest dynamics plots in montane wet forest and lowland dry forest on Hawai‘i Island. We compared the species...

Avaliação da conduta conservadora na lesão intraepitelial cervical de alto grau

Directory of Open Access Journals (Sweden)

Nelson Shozo Uchimura

2012-06-01

Full Text Available OBJETIVO: Analisar a associação entre a conduta conservadora em lesão intraepitelial cervical de alto grau com o índice de recidiva da neoplasia e faixa etária. MÉTODOS: Estudo transversal e retrospectivo realizado com 509 mulheres (15-76 anos atendidas no período de 1996 a 2006, com colpocitologia oncótica alterada, em um serviço público de referência em Maringá, PR. Os dados foram coletados dos prontuários médicos e estudadas as variáveis diagnóstico definitivo, tipos de tratamento, ocorrência da lesão e recidivas, analisados por meio de testes de associação de qui-quadrado de Pearson e teste exato de Fisher. RESULTADOS: A lesão intraepitelial cervical de alto grau ocorreu em 168 casos; destes, 31 mulheres foram submetidas à amputação cônica, 104 a cirurgias de alta frequência, nove histerectomizadas e 24 receberam conduta conservadora. Dentre as mulheres com lesão de alto grau e tratadas de forma conservadora, oito (33,3% recidivaram, enquanto dentre as submetidas à conduta não conservadora dez (6,9% recidivaram, sendo essa diferença estatisticamente significante (p = 0,0009, RP = 4,8 (IC95% 2,11;10,93. Para aquelas que fizeram o seguimento clínico-citológico, três (30,0% e, dentre as cauterizadas, cinco (35,7% recidivaram no prazo de três anos, sem diferença significante (p = 0,5611. A recidiva abaixo e acima de 30 anos ocorreu, respectivamente, em sete (13,8% e 11 (12,2% mulheres (p = 0,9955. CONCLUSÕES: A idade da mulher não influencia o prognóstico de recidiva. O tratamento conservador deve ser indicado como conduta de exceção, dada a alta taxa de recidiva, e o seguimento deve ser rigoroso, com acompanhamento citológico e colposcópico de até três anos, período em que ocorre a maioria das recidivas.

Does participatory forest management promote sustainable forest utilisation in Tanzania?

DEFF Research Database (Denmark)

Treue, Thorsten; Ngaga, Y.M.; Meilby, Henrik

2014-01-01

Over the past 20 years, Participatory Forest Management (PFM) has become a dominant forest management strategy in Tanzania, covering more than 4.1 million hectares. Sustainable forest use and supply of wood products to local people are major aims of PFM. This paper assesses the sustainability...... of forest utilisation under PFM, using estimates of forest condition and extraction rates based on forest inventories and 480 household surveys from 12 forests; seven under Community Based Forest Management (CBFM), three under Joint Forest Management (JFM) and two under government management (non......-PFM). Extraction of products is intense in forests close to Dar es Salaam, regardless of management regime. Further from Dar es Salaam, harvesting levels in forests under PFM are, with one prominent exception, broadly sustainable. Using GIS data from 116 wards, it is shown that half of the PFM forests in Tanzania...

Iterative acceleration methods for Monte Carlo and deterministic criticality calculations

Energy Technology Data Exchange (ETDEWEB)

Urbatsch, T.J.

1995-11-01

If you have ever given up on a nuclear criticality calculation and terminated it because it took so long to converge, you might find this thesis of interest. The author develops three methods for improving the fission source convergence in nuclear criticality calculations for physical systems with high dominance ratios for which convergence is slow. The Fission Matrix Acceleration Method and the Fission Diffusion Synthetic Acceleration (FDSA) Method are acceleration methods that speed fission source convergence for both Monte Carlo and deterministic methods. The third method is a hybrid Monte Carlo method that also converges for difficult problems where the unaccelerated Monte Carlo method fails. The author tested the feasibility of all three methods in a test bed consisting of idealized problems. He has successfully accelerated fission source convergence in both deterministic and Monte Carlo criticality calculations. By filtering statistical noise, he has incorporated deterministic attributes into the Monte Carlo calculations in order to speed their source convergence. He has used both the fission matrix and a diffusion approximation to perform unbiased accelerations. The Fission Matrix Acceleration method has been implemented in the production code MCNP and successfully applied to a real problem. When the unaccelerated calculations are unable to converge to the correct solution, they cannot be accelerated in an unbiased fashion. A Hybrid Monte Carlo method weds Monte Carlo and a modified diffusion calculation to overcome these deficiencies. The Hybrid method additionally possesses reduced statistical errors.

Iterative acceleration methods for Monte Carlo and deterministic criticality calculations

International Nuclear Information System (INIS)

Urbatsch, T.J.

1995-11-01

If you have ever given up on a nuclear criticality calculation and terminated it because it took so long to converge, you might find this thesis of interest. The author develops three methods for improving the fission source convergence in nuclear criticality calculations for physical systems with high dominance ratios for which convergence is slow. The Fission Matrix Acceleration Method and the Fission Diffusion Synthetic Acceleration (FDSA) Method are acceleration methods that speed fission source convergence for both Monte Carlo and deterministic methods. The third method is a hybrid Monte Carlo method that also converges for difficult problems where the unaccelerated Monte Carlo method fails. The author tested the feasibility of all three methods in a test bed consisting of idealized problems. He has successfully accelerated fission source convergence in both deterministic and Monte Carlo criticality calculations. By filtering statistical noise, he has incorporated deterministic attributes into the Monte Carlo calculations in order to speed their source convergence. He has used both the fission matrix and a diffusion approximation to perform unbiased accelerations. The Fission Matrix Acceleration method has been implemented in the production code MCNP and successfully applied to a real problem. When the unaccelerated calculations are unable to converge to the correct solution, they cannot be accelerated in an unbiased fashion. A Hybrid Monte Carlo method weds Monte Carlo and a modified diffusion calculation to overcome these deficiencies. The Hybrid method additionally possesses reduced statistical errors

Environmental and ecological architects: Guidelines for the Chilean temperate rainforest management derived from the monito del monte (Dromiciops gliroides conservation Arquitectos ambientales y ecológicos: Pautas para la gestión ambiental del bosque templado lluvioso de Chile derivadas de la conservación del monito del monte (Dromiciops gliroides

Directory of Open Access Journals (Sweden)

FRANCISCO E FONTÚRBEL

2011-06-01

Full Text Available Land use change is one of the main biodiversity threats. Due to this change, natural habitats such as the South American temperate rainforests have been rapidly degraded, fragmented, and lost. Consequently, the management and conservation of the remaining forest is a priority and having an appropriate environmental policy is mandatory for this purpose. Conservation actions in the temperate rainforest have been addressed from an individual species perspective, giving less attention to the ecosystem level conservation. Moreover, conservation-related information has not reflected yet on environmental policy development. We used the case study of the monito del monte (Dromiciops gliroides to illustrate how our current ecological and conservation knowledge of a species could be used to generate a new environmental policy. Dromiciops gliroides is a forest-dependent species with an important ecological role and quite unique evolutionary status. In order to guarantee the persistence of D. gliroides, we propose two habitat management components to be incorporated in management plans: structure and connectivity. Structure refers to spatial arrangement and key structural elements that determine habitat quality and connectivity refers to functional connectivity at the landscape level. The conservation of the monito del monte might also contribute to the conservation of many other forest-dependent species. By conserving such species it will be possible to conserve the ecological interactions and the eco-evolutionary processes, which ultimately determine the conservation of the temperate rainforest.El cambio de uso de la tierra es una de las principales amenazas a la biodiversidad. Hábitats naturales como el bosque templado lluvioso de Sudamérica han sido degradados, fragmentados y eliminados rápidamente. En consecuencia, es prioritario el manejo y conservación de los remanentes, para lo que es necesario contar con una política ambiental acorde. Las

Selection of sweet potato clones for the region Alto Vale do Jequitinhonha Seleção de clones de batata-doce para a região do Alto Vale do Jequitinhonha

Directory of Open Access Journals (Sweden)

Valter C de Andrade Júnior

2009-09-01

Full Text Available An experiment was carried out from December 2005 to July 2006, in the Universidade Federal dos Vales do Jequitinhonha e Mucuri (UFVJM, in Diamantina, Minas Gerais State, Brazil, aiming at selecting sweet potato clones for the Alto Vale do Jequitinhonha. We evaluated nine clones from the UFVJM germplasm bank, using cultivars Brazlândia Branca, Brazlândia Roxa, and Princesa as controls. The experimental design was blocks at random, with four replications. Plants were harvested seven months after transplanting. We assessed the fresh mass yield of vines and roots, as well as root shape and resistance to soil insects. Genotypes did not differ from each other for the fresh mass yield of vines (ranging from 3.81 to 11.76 t ha-1. The total yield of roots ranged from 22.0 to 45.4 t ha-1 and clones BD-06, BD-113-TO, BD-15, BD-38, BD-25, BD-61, and cultivar Princesa had statistically the highest figures. However, only clone BD-06 significantly overcame the control cultivars Brazlândia Branca and Brazlândia Roxa. Clone BD-06 had also the highest commercial yield of roots (38.58 t ha-1, statically similar to most of the other clones and cultivar Princesa (25.87 t ha-1, but superior to cultivars Brazlândia Branca and Brazlândia Roxa. Most of the clones tested, including clone BD-06, produced good shaped roots and were resistant to soil insects. Considering our results, clone BD-06 stood out as a good option for growing sweet potato in the Upper Valley of Jequitinhonha.Com o objetivo de selecionar clones de batata-doce para a região do alto Vale do Jequitinhonha, conduziu-se um experimento de dezembro de 2005 a julho de 2006, no CampusJK da Universidade Federal dos Vales do Jequitinhonha e Mucuri, UFVJM, município de Diamantina-MG. Foram avaliados nove clones de batata-doce pertencentes ao banco de germoplasma da UFVJM, juntamente com as cultivares Brazlândia Branca, Brazlândia Roxa e Princesa, utilizadas como testemunhas. O delineamento experimental

Estudo de Caso de Sistemas Meteorológicos em Altos Níveis Sobre a América do Sul: Análise Sinótica

Directory of Open Access Journals (Sweden)

André Becker Nunes

2017-01-01

Full Text Available O Vórtice ciclônico de altos níveis (VCAN sobre o nordeste brasileiro (NeB e a Alta da Bolívia (AB são os típicos sistemas de pressão em altos níveis (250 hPa da América do Sul, especialmente nas estações mais quentes. Somados a eles, também pode-se considerar o VCAN extratropical, de característica transiente e encontrado o ano todo. Neste trabalho, é apresentada uma análise sinótica, baseada na literatura especializada, de um caso observado em 01 de fevereiro de 2005 em que os três sistemas ocorreram simultaneamente. Algumas características encontradas na literatura foram observadas nos sistemas deste caso. O campo de temperatura mais associado ao sistema extratropical foi o de 500 hPa, e não o de 250 hPa. Os movimentos verticais mais relevantes foram encontrados somente no sistema extratropical: intensos movimentos ascendentes à sotavento do eixo do cavado em altos níveis. Não se encontrou um sistema em superfície especificamente associado ao VCAN do NeB, enquanto que se observou a baixa continental (ou Baixa do Chaco associada à Alta da Bolívia e um intenso ciclone extratropical associado ao VCAN extratropical. Tal ciclone em superfície acarretou em vento sul e consequente resfriamento da camada, o que pôde ser observado pela amplitude do cavado em altos níveis, o que é raro para esta época do ano.

Urban Forests

Science.gov (United States)

David Nowak

2016-01-01

Urban forests (and trees) constitute the second forest resource considered in this report. We specifically emphasize the fact that agricultural and urban forests exist on a continuum defined by their relationship (and interrelationship) with a given landscape. These two forest types generally serve different purposes, however. Whereas agricultural forests are...

Study on random number generator in Monte Carlo code

International Nuclear Information System (INIS)

Oya, Kentaro; Kitada, Takanori; Tanaka, Shinichi

2011-01-01

The Monte Carlo code uses a sequence of pseudo-random numbers with a random number generator (RNG) to simulate particle histories. A pseudo-random number has its own period depending on its generation method and the period is desired to be long enough not to exceed the period during one Monte Carlo calculation to ensure the correctness especially for a standard deviation of results. The linear congruential generator (LCG) is widely used as Monte Carlo RNG and the period of LCG is not so long by considering the increasing rate of simulation histories in a Monte Carlo calculation according to the remarkable enhancement of computer performance. Recently, many kinds of RNG have been developed and some of their features are better than those of LCG. In this study, we investigate the appropriate RNG in a Monte Carlo code as an alternative to LCG especially for the case of enormous histories. It is found that xorshift has desirable features compared with LCG, and xorshift has a larger period, a comparable speed to generate random numbers, a better randomness, and good applicability to parallel calculation. (author)

Rb-Sr geochronology from Barro Alto Complex, Goias: metamorphism evidence of high degree and continental collision around 1300 Ma ago in Central Brazil

International Nuclear Information System (INIS)

Fuck, R.A.; Neves, B.B.B.; Cordani, U.G.; Kawashita, K.

1989-01-01

Rb-Sr geochronologic investigation carried out on rocks from the Barro Alto Complex, Goias, yielded iso chronic ages of 1266 +- 17 Ma, for felsic rocks from the granulite belt and 1330 +- 67 Ma, for gneisses belonging to the Juscelandia Sequence. Rb-Sr isotope measurements suggest that Barro Alto rocks have undergone an important metamorphic event during middle Proterozoic times, around 1300 Ma ago. During that event, volcanic and sedimentary rocks of Juscelandia Sequence, as well as the underlying gabbros-anorthosite layered complex, underwent deformation and recrystallization under amphibolite facies conditions. Deformation and metamorphism took place during the collision of two continental blocks, which resulted in a southeastward directed thrust complex, allowing the exposure of granulite slices from the middle-lower crust of the overthrusted block. (author)

Combinatorial geometry domain decomposition strategies for Monte Carlo simulations

Energy Technology Data Exchange (ETDEWEB)

Li, G.; Zhang, B.; Deng, L.; Mo, Z.; Liu, Z.; Shangguan, D.; Ma, Y.; Li, S.; Hu, Z. [Institute of Applied Physics and Computational Mathematics, Beijing, 100094 (China)

2013-07-01

Analysis and modeling of nuclear reactors can lead to memory overload for a single core processor when it comes to refined modeling. A method to solve this problem is called 'domain decomposition'. In the current work, domain decomposition algorithms for a combinatorial geometry Monte Carlo transport code are developed on the JCOGIN (J Combinatorial Geometry Monte Carlo transport INfrastructure). Tree-based decomposition and asynchronous communication of particle information between domains are described in the paper. Combination of domain decomposition and domain replication (particle parallelism) is demonstrated and compared with that of MERCURY code. A full-core reactor model is simulated to verify the domain decomposition algorithms using the Monte Carlo particle transport code JMCT (J Monte Carlo Transport Code), which has being developed on the JCOGIN infrastructure. Besides, influences of the domain decomposition algorithms to tally variances are discussed. (authors)

Combinatorial geometry domain decomposition strategies for Monte Carlo simulations

International Nuclear Information System (INIS)

Li, G.; Zhang, B.; Deng, L.; Mo, Z.; Liu, Z.; Shangguan, D.; Ma, Y.; Li, S.; Hu, Z.

2013-01-01

Analysis and modeling of nuclear reactors can lead to memory overload for a single core processor when it comes to refined modeling. A method to solve this problem is called 'domain decomposition'. In the current work, domain decomposition algorithms for a combinatorial geometry Monte Carlo transport code are developed on the JCOGIN (J Combinatorial Geometry Monte Carlo transport INfrastructure). Tree-based decomposition and asynchronous communication of particle information between domains are described in the paper. Combination of domain decomposition and domain replication (particle parallelism) is demonstrated and compared with that of MERCURY code. A full-core reactor model is simulated to verify the domain decomposition algorithms using the Monte Carlo particle transport code JMCT (J Monte Carlo Transport Code), which has being developed on the JCOGIN infrastructure. Besides, influences of the domain decomposition algorithms to tally variances are discussed. (authors)

An analysis on the costs of Belo Monte hydroelectric power plant; Uma analise sobre os custos da hidreletrica Belo Monte

Energy Technology Data Exchange (ETDEWEB)

Silva, Marcos Vinicius Miranda da [Universidade de Sao Paulo (PIPGE/USP), SP (Brazil). Programa Interunidades de Pos-Graduacao em Energia], e-mail: energiapara@yahoo.com.br

2008-07-01

The Belo Monte hydropower plant's low generation cost is among the arguments used by Centrais Eletricas do Norte do Brazil (ELETRONORTE), a Brazilian state electric utility, to make possible its construction. This paper shows that the generation cost presented by ELETROBRAS is very low in relation to the world pattern of cost and probably unrealistic. It also shows that the generation cost cannot be used separately to determine the Belo Monte dam's economic feasibility. There is the need to include other costs, such as: socio environmental degradation and control, financial compensation for using the hydraulic resources, transmission and thermal backup stations, beyond, evidently, generation cost for assuring the credibility of the Belo Monte hydropower plant's economic analysis. (author)

Monte Carlo method applied to medical physics

International Nuclear Information System (INIS)

Oliveira, C.; Goncalves, I.F.; Chaves, A.; Lopes, M.C.; Teixeira, N.; Matos, B.; Goncalves, I.C.; Ramalho, A.; Salgado, J.

2000-01-01

The main application of the Monte Carlo method to medical physics is dose calculation. This paper shows some results of two dose calculation studies and two other different applications: optimisation of neutron field for Boron Neutron Capture Therapy and optimization of a filter for a beam tube for several purposes. The time necessary for Monte Carlo calculations - the highest boundary for its intensive utilisation - is being over-passed with faster and cheaper computers. (author)

A radiating shock evaluated using Implicit Monte Carlo Diffusion

International Nuclear Information System (INIS)

Cleveland, M.; Gentile, N.

2013-01-01

Implicit Monte Carlo [1] (IMC) has been shown to be very expensive when used to evaluate a radiation field in opaque media. Implicit Monte Carlo Diffusion (IMD) [2], which evaluates a spatial discretized diffusion equation using a Monte Carlo algorithm, can be used to reduce the cost of evaluating the radiation field in opaque media [2]. This work couples IMD to the hydrodynamics equations to evaluate opaque diffusive radiating shocks. The Lowrie semi-analytic diffusive radiating shock benchmark[a] is used to verify our implementation of the coupled system of equations. (authors)

The Monte Carlo method the method of statistical trials

CERN Document Server

Shreider, YuA

1966-01-01

The Monte Carlo Method: The Method of Statistical Trials is a systematic account of the fundamental concepts and techniques of the Monte Carlo method, together with its range of applications. Some of these applications include the computation of definite integrals, neutron physics, and in the investigation of servicing processes. This volume is comprised of seven chapters and begins with an overview of the basic features of the Monte Carlo method and typical examples of its application to simple problems in computational mathematics. The next chapter examines the computation of multi-dimensio

Contribution to the knowledge of the fungal biodiversity of Ordesa and Monte Perdido National Park II

Directory of Open Access Journals (Sweden)

F. Pancorbo

2017-11-01

Full Text Available This paper continues the taxonomic revision of the species collected during the campaign of 2015 in the National Park of Ordesa and Monte Perdido. 409 taxa are added to the previous inventory, some of them from alpine-subalpine ecology, of which 76% correspond to phylum Basidiomycota and 22.2% to Ascomycota. They are presented in the form of a check-list, followed by 20 selected taxonomic descriptions of interesting, infrequent species, and those though to be new in the peninsular territory. Among the species determined, six species were included in the proposals for the inventory of protected species of the Iberian Peninsula and/or Aragon. This paper presents a first approximation, as a platform for later evaluations, of the beech conservation degree in some forests from the Park through the occurrence of indicator species.

Phylobetadiversity among forest types in the Brazilian Atlantic Forest complex.

Science.gov (United States)

Duarte, Leandro Da Silva; Bergamin, Rodrigo Scarton; Marcilio-Silva, Vinícius; Seger, Guilherme Dubal Dos Santos; Marques, Márcia Cristina Mendes

2014-01-01

Phylobetadiversity is defined as the phylogenetic resemblance between communities or biomes. Analyzing phylobetadiversity patterns among different vegetation physiognomies within a single biome is crucial to understand the historical affinities between them. Based on the widely accepted idea that different forest physiognomies within the Southern Brazilian Atlantic Forest constitute different facies of a single biome, we hypothesize that more recent phylogenetic nodes should drive phylobetadiversity gradients between the different forest types within the Atlantic Forest, as the phylogenetic divergence among those forest types is biogeographically recent. We compiled information from 206 checklists describing the occurrence of shrub/tree species across three different forest physiognomies within the Southern Brazilian Atlantic Forest (Dense, Mixed and Seasonal forests). We analyzed intra-site phylogenetic structure (phylogenetic diversity, net relatedness index and nearest taxon index) and phylobetadiversity between plots located at different forest types, using five different methods differing in sensitivity to either basal or terminal nodes (phylogenetic fuzzy weighting, COMDIST, COMDISTNT, UniFrac and Rao's H). Mixed forests showed higher phylogenetic diversity and overdispersion than the other forest types. Furthermore, all forest types differed from each other in relation phylobetadiversity patterns, particularly when phylobetadiversity methods more sensitive to terminal nodes were employed. Mixed forests tended to show higher phylogenetic differentiation to Dense and Seasonal forests than these latter from each other. The higher phylogenetic diversity and phylobetadiversity levels found in Mixed forests when compared to the others likely result from the biogeographical origin of several taxa occurring in these forests. On one hand, Mixed forests shelter several temperate taxa, like the conifers Araucaria and Podocarpus. On the other hand, tropical groups, like

Lípidos séricos en escolares y adolescentes sanos chilenos de estrato socioeconómico alto

Directory of Open Access Journals (Sweden)

Zacarías S. José, Dr.

2012-11-01

Conclusiones: En niños y adolescentes chilenos sanos de estratos socioeconómicos altos, hay una alta proporción con colesterol total y colesterol LDL en zonas de riesgo o aumentadas. Las concentraciones de colesterol de adolescentes mujeres son significativamente mayores que las de varones. No se observó una asociación de lípidos séricos con perímetro abdominal o con IMC.

Automated Monte Carlo biasing for photon-generated electrons near surfaces.

Energy Technology Data Exchange (ETDEWEB)

Franke, Brian Claude; Crawford, Martin James; Kensek, Ronald Patrick

2009-09-01

This report describes efforts to automate the biasing of coupled electron-photon Monte Carlo particle transport calculations. The approach was based on weight-windows biasing. Weight-window settings were determined using adjoint-flux Monte Carlo calculations. A variety of algorithms were investigated for adaptivity of the Monte Carlo tallies. Tree data structures were used to investigate spatial partitioning. Functional-expansion tallies were used to investigate higher-order spatial representations.

Applicability of quasi-Monte Carlo for lattice systems

International Nuclear Information System (INIS)

Ammon, Andreas; Deutsches Elektronen-Synchrotron; Hartung, Tobias; Jansen, Karl; Leovey, Hernan; Griewank, Andreas; Mueller-Preussker, Michael

2013-11-01

This project investigates the applicability of quasi-Monte Carlo methods to Euclidean lattice systems in order to improve the asymptotic error scaling of observables for such theories. The error of an observable calculated by averaging over random observations generated from ordinary Monte Carlo simulations scales like N -1/2 , where N is the number of observations. By means of quasi-Monte Carlo methods it is possible to improve this scaling for certain problems to N -1 , or even further if the problems are regular enough. We adapted and applied this approach to simple systems like the quantum harmonic and anharmonic oscillator and verified an improved error scaling of all investigated observables in both cases.

Applicability of quasi-Monte Carlo for lattice systems

Energy Technology Data Exchange (ETDEWEB)

Ammon, Andreas [Berlin Humboldt-Univ. (Germany). Dept. of Physics; Deutsches Elektronen-Synchrotron (DESY), Zeuthen (Germany). John von Neumann-Inst. fuer Computing NIC; Hartung, Tobias [King' s College London (United Kingdom). Dept. of Mathematics; Jansen, Karl [Deutsches Elektronen-Synchrotron (DESY), Zeuthen (Germany). John von Neumann-Inst. fuer Computing NIC; Leovey, Hernan; Griewank, Andreas [Berlin Humboldt-Univ. (Germany). Dept. of Mathematics; Mueller-Preussker, Michael [Berlin Humboldt-Univ. (Germany). Dept. of Physics

2013-11-15

This project investigates the applicability of quasi-Monte Carlo methods to Euclidean lattice systems in order to improve the asymptotic error scaling of observables for such theories. The error of an observable calculated by averaging over random observations generated from ordinary Monte Carlo simulations scales like N{sup -1/2}, where N is the number of observations. By means of quasi-Monte Carlo methods it is possible to improve this scaling for certain problems to N{sup -1}, or even further if the problems are regular enough. We adapted and applied this approach to simple systems like the quantum harmonic and anharmonic oscillator and verified an improved error scaling of all investigated observables in both cases.

Uniform distribution and quasi-Monte Carlo methods discrepancy, integration and applications

CERN Document Server

Kritzer, Peter; Pillichshammer, Friedrich; Winterhof, Arne

2014-01-01

The survey articles in this book focus on number theoretic point constructions, uniform distribution theory, and quasi-Monte Carlo methods. As deterministic versions of the Monte Carlo method, quasi-Monte Carlo rules enjoy increasing popularity, with many fruitful applications in mathematical practice, as for example in finance, computer graphics, and biology.

Clinical implementation of full Monte Carlo dose calculation in proton beam therapy

International Nuclear Information System (INIS)

Paganetti, Harald; Jiang, Hongyu; Parodi, Katia; Slopsema, Roelf; Engelsman, Martijn

2008-01-01

The goal of this work was to facilitate the clinical use of Monte Carlo proton dose calculation to support routine treatment planning and delivery. The Monte Carlo code Geant4 was used to simulate the treatment head setup, including a time-dependent simulation of modulator wheels (for broad beam modulation) and magnetic field settings (for beam scanning). Any patient-field-specific setup can be modeled according to the treatment control system of the facility. The code was benchmarked against phantom measurements. Using a simulation of the ionization chamber reading in the treatment head allows the Monte Carlo dose to be specified in absolute units (Gy per ionization chamber reading). Next, the capability of reading CT data information was implemented into the Monte Carlo code to model patient anatomy. To allow time-efficient dose calculation, the standard Geant4 tracking algorithm was modified. Finally, a software link of the Monte Carlo dose engine to the patient database and the commercial planning system was established to allow data exchange, thus completing the implementation of the proton Monte Carlo dose calculation engine ('DoC++'). Monte Carlo re-calculated plans are a valuable tool to revisit decisions in the planning process. Identification of clinically significant differences between Monte Carlo and pencil-beam-based dose calculations may also drive improvements of current pencil-beam methods. As an example, four patients (29 fields in total) with tumors in the head and neck regions were analyzed. Differences between the pencil-beam algorithm and Monte Carlo were identified in particular near the end of range, both due to dose degradation and overall differences in range prediction due to bony anatomy in the beam path. Further, the Monte Carlo reports dose-to-tissue as compared to dose-to-water by the planning system. Our implementation is tailored to a specific Monte Carlo code and the treatment planning system XiO (Computerized Medical Systems Inc

Forest tenure and sustainable forest management

Science.gov (United States)

J.P. Siry; K. McGinley; F.W. Cubbage; P. Bettinger

2015-01-01

We reviewed the principles and key literature related to forest tenure and sustainable forest management, and then examined the status of sustainable forestry and land ownership at the aggregate national level for major forested countries. The institutional design principles suggested by Ostrom are well accepted for applications to public, communal, and private lands....

Structure of neutron rich nuclei of Germanium and Gallium beyond N equals 50 at Alto

International Nuclear Information System (INIS)

Lebois, M.

2008-09-01

The gamma rays following the beta decay of the following very neutron-rich isotopes: 82,83,84 Ga produced by photo-fission, have been studied at the newly built ISOL facility in Orsay: ALTO. In ALTO the interaction of an electron beam with U 238 target generates a continuous spectra of Bremsstrahlung gamma radiation that triggers U 238 fission. The fission fragments are then ionized, extracted and mass-separated. The analysis of the data has shown the existence of an isomer in 31 84 Ga 53 and has enabled us to confirm known results on 32 83 Ge 51 energy levels including the gamma transition between the 1/2+ state at 247,7 KeV and the fundamental state. We have also proposed the first energy level scheme for 33 84 As 51 . In order to understand the structure of the nucleus we have used the Thankappan and True model that gives a description of the coupling between the pair-pair core (half-magical) and the single nucleon. This model applied to the N=51 chain ( 38 89 Sr 51 , 36 87 Kr 51 , 34 85 Se 51 , 32 83 Ge 51 and 30 81 Zn 51 ) has allowed us to see the main features of odd isotope structure. We have also confirmed previous results concerning the nature of the states in the following decay 31 83 Ga 52 → 32 83 Ge 51

Exponential convergence on a continuous Monte Carlo transport problem

International Nuclear Information System (INIS)

Booth, T.E.

1997-01-01

For more than a decade, it has been known that exponential convergence on discrete transport problems was possible using adaptive Monte Carlo techniques. An adaptive Monte Carlo method that empirically produces exponential convergence on a simple continuous transport problem is described

A flexible coupling scheme for Monte Carlo and thermal-hydraulics codes

Energy Technology Data Exchange (ETDEWEB)

Hoogenboom, J. Eduard, E-mail: J.E.Hoogenboom@tudelft.nl [Delft University of Technology (Netherlands); Ivanov, Aleksandar; Sanchez, Victor, E-mail: Aleksandar.Ivanov@kit.edu, E-mail: Victor.Sanchez@kit.edu [Karlsruhe Institute of Technology, Institute of Neutron Physics and Reactor Technology, Eggenstein-Leopoldshafen (Germany); Diop, Cheikh, E-mail: Cheikh.Diop@cea.fr [CEA/DEN/DANS/DM2S/SERMA, Commissariat a l' Energie Atomique, Gif-sur-Yvette (France)

2011-07-01

A coupling scheme between a Monte Carlo code and a thermal-hydraulics code is being developed within the European NURISP project for comprehensive and validated reactor analysis. The scheme is flexible as it allows different Monte Carlo codes and different thermal-hydraulics codes to be used. At present the MCNP and TRIPOLI4 Monte Carlo codes can be used and the FLICA4 and SubChanFlow thermal-hydraulics codes. For all these codes only an original executable is necessary. A Python script drives the iterations between Monte Carlo and thermal-hydraulics calculations. It also calls a conversion program to merge a master input file for the Monte Carlo code with the appropriate temperature and coolant density data from the thermal-hydraulics calculation. Likewise it calls another conversion program to merge a master input file for the thermal-hydraulics code with the power distribution data from the Monte Carlo calculation. Special attention is given to the neutron cross section data for the various required temperatures in the Monte Carlo calculation. Results are shown for an infinite lattice of PWR fuel pin cells and a 3 x 3 fuel BWR pin cell cluster. Various possibilities for further improvement and optimization of the coupling system are discussed. (author)

A flexible coupling scheme for Monte Carlo and thermal-hydraulics codes

International Nuclear Information System (INIS)

Hoogenboom, J. Eduard; Ivanov, Aleksandar; Sanchez, Victor; Diop, Cheikh

2011-01-01

A coupling scheme between a Monte Carlo code and a thermal-hydraulics code is being developed within the European NURISP project for comprehensive and validated reactor analysis. The scheme is flexible as it allows different Monte Carlo codes and different thermal-hydraulics codes to be used. At present the MCNP and TRIPOLI4 Monte Carlo codes can be used and the FLICA4 and SubChanFlow thermal-hydraulics codes. For all these codes only an original executable is necessary. A Python script drives the iterations between Monte Carlo and thermal-hydraulics calculations. It also calls a conversion program to merge a master input file for the Monte Carlo code with the appropriate temperature and coolant density data from the thermal-hydraulics calculation. Likewise it calls another conversion program to merge a master input file for the thermal-hydraulics code with the power distribution data from the Monte Carlo calculation. Special attention is given to the neutron cross section data for the various required temperatures in the Monte Carlo calculation. Results are shown for an infinite lattice of PWR fuel pin cells and a 3 x 3 fuel BWR pin cell cluster. Various possibilities for further improvement and optimization of the coupling system are discussed. (author)

Análisis de la operación del horno alto mediante la exergía térmica

Directory of Open Access Journals (Sweden)

Babich, A.

1998-05-01

Full Text Available In order to know the blast furnace thermic state, the thermic exergy can be used, which considers the temperatures level of the process. The magnitude analysis of thermic exergy for blast furnace operation using pulverized coal, natural gas and oxigen injection by tuyeres, shows that natural gas injection is less effective than pulverized coal and enriched blast with oxygen simultaneous injection.

Para conocer el estado térmico del horno alto se puede utilizar la exergía térmica, que considera el nivel de temperaturas del proceso. El análisis de la magnitud de la exergía térmica, para la operación del horno alto, utilizando inyecciones por toberas de carbón pulverizado, de gas natural y viento enriquecido en oxígeno, en ambos casos muestra que la inyección de gas natural es menos efectiva que la inyección simultánea de carbón pulverizado y viento enriquecido con oxígeno.

Isotopic depletion with Monte Carlo

International Nuclear Information System (INIS)

Martin, W.R.; Rathkopf, J.A.

1996-06-01

This work considers a method to deplete isotopes during a time- dependent Monte Carlo simulation of an evolving system. The method is based on explicitly combining a conventional estimator for the scalar flux with the analytical solutions to the isotopic depletion equations. There are no auxiliary calculations; the method is an integral part of the Monte Carlo calculation. The method eliminates negative densities and reduces the variance in the estimates for the isotope densities, compared to existing methods. Moreover, existing methods are shown to be special cases of the general method described in this work, as they can be derived by combining a high variance estimator for the scalar flux with a low-order approximation to the analytical solution to the depletion equation

Multilevel sequential Monte-Carlo samplers

KAUST Repository

Jasra, Ajay

2016-01-01

Multilevel Monte-Carlo methods provide a powerful computational technique for reducing the computational cost of estimating expectations for a given computational effort. They are particularly relevant for computational problems when approximate distributions are determined via a resolution parameter h, with h=0 giving the theoretical exact distribution (e.g. SDEs or inverse problems with PDEs). The method provides a benefit by coupling samples from successive resolutions, and estimating differences of successive expectations. We develop a methodology that brings Sequential Monte-Carlo (SMC) algorithms within the framework of the Multilevel idea, as SMC provides a natural set-up for coupling samples over different resolutions. We prove that the new algorithm indeed preserves the benefits of the multilevel principle, even if samples at all resolutions are now correlated.

Monte Carlo methods in ICF

International Nuclear Information System (INIS)

Zimmerman, G.B.

1997-01-01

Monte Carlo methods appropriate to simulate the transport of x-rays, neutrons, ions and electrons in Inertial Confinement Fusion targets are described and analyzed. The Implicit Monte Carlo method of x-ray transport handles symmetry within indirect drive ICF hohlraums well, but can be improved 50X in efficiency by angular biasing the x-rays towards the fuel capsule. Accurate simulation of thermonuclear burn and burn diagnostics involves detailed particle source spectra, charged particle ranges, inflight reaction kinematics, corrections for bulk and thermal Doppler effects and variance reduction to obtain adequate statistics for rare events. It is found that the effects of angular Coulomb scattering must be included in models of charged particle transport through heterogeneous materials. copyright 1997 American Institute of Physics

Multilevel sequential Monte-Carlo samplers

KAUST Repository

Jasra, Ajay

2016-01-05

Multilevel Monte-Carlo methods provide a powerful computational technique for reducing the computational cost of estimating expectations for a given computational effort. They are particularly relevant for computational problems when approximate distributions are determined via a resolution parameter h, with h=0 giving the theoretical exact distribution (e.g. SDEs or inverse problems with PDEs). The method provides a benefit by coupling samples from successive resolutions, and estimating differences of successive expectations. We develop a methodology that brings Sequential Monte-Carlo (SMC) algorithms within the framework of the Multilevel idea, as SMC provides a natural set-up for coupling samples over different resolutions. We prove that the new algorithm indeed preserves the benefits of the multilevel principle, even if samples at all resolutions are now correlated.

Adesão ao Tratamento em Gestação de Alto Risco

Directory of Open Access Journals (Sweden)

Fabíola Langaro

Full Text Available No Brasil, aproximadamente 15% das gestações são de alto risco, sendo os diagnósticos de diabetes gestacional e hipertensão as causas mais frequentes dessa condição. Para a efetividade dos tratamentos indicados nesses casos, bem como alcance de desfechos favoráveis, é essencial a adesão da paciente às recomendações fornecidas pelas equipes de saúde, o que torna fundamental a avaliação desse indicador. A pesquisa apresentada avaliou a adesão de participantes de um programa de atendimento interdisciplinar a gestantes de alto risco desenvolvido em um hospital geral privado por meio de três instrumentos: Questionário Morisky-Green de adesão ao uso de medicação, Questionário de Adesão ao Tratamento e Pergunta descritiva de avaliação qualitativa. Durante oito meses, 83 gestantes responderam aos questionários, que possibilitaram identificar níveis de adesão ao longo do tempo de participação no programa, situações de não adesão mais frequentes, bem como variáveis que interferem nas medidas de adesão ao uso de medicação e ao tratamento global. Considera-se que os índices levantados são passíveis de uso como indicadores clínicos e gerenciais, possibilitando a comprovação da efetividade e melhoria constante das rotinas de saúde. Tais mecanismos são úteis para definição de protocolos, discussão de casos clínicos e mesmo para feedbackàs próprias pacientes.

Parallel MCNP Monte Carlo transport calculations with MPI

International Nuclear Information System (INIS)

Wagner, J.C.; Haghighat, A.

1996-01-01

The steady increase in computational performance has made Monte Carlo calculations for large/complex systems possible. However, in order to make these calculations practical, order of magnitude increases in performance are necessary. The Monte Carlo method is inherently parallel (particles are simulated independently) and thus has the potential for near-linear speedup with respect to the number of processors. Further, the ever-increasing accessibility of parallel computers, such as workstation clusters, facilitates the practical use of parallel Monte Carlo. Recognizing the nature of the Monte Carlo method and the trends in available computing, the code developers at Los Alamos National Laboratory implemented the message-passing general-purpose Monte Carlo radiation transport code MCNP (version 4A). The PVM package was chosen by the MCNP code developers because it supports a variety of communication networks, several UNIX platforms, and heterogeneous computer systems. This PVM version of MCNP has been shown to produce speedups that approach the number of processors and thus, is a very useful tool for transport analysis. Due to software incompatibilities on the local IBM SP2, PVM has not been available, and thus it is not possible to take advantage of this useful tool. Hence, it became necessary to implement an alternative message-passing library package into MCNP. Because the message-passing interface (MPI) is supported on the local system, takes advantage of the high-speed communication switches in the SP2, and is considered to be the emerging standard, it was selected

Monte Carlo systems used for treatment planning and dose verification

Energy Technology Data Exchange (ETDEWEB)

Brualla, Lorenzo [Universitaetsklinikum Essen, NCTeam, Strahlenklinik, Essen (Germany); Rodriguez, Miguel [Centro Medico Paitilla, Balboa (Panama); Lallena, Antonio M. [Universidad de Granada, Departamento de Fisica Atomica, Molecular y Nuclear, Granada (Spain)

2017-04-15

General-purpose radiation transport Monte Carlo codes have been used for estimation of the absorbed dose distribution in external photon and electron beam radiotherapy patients since several decades. Results obtained with these codes are usually more accurate than those provided by treatment planning systems based on non-stochastic methods. Traditionally, absorbed dose computations based on general-purpose Monte Carlo codes have been used only for research, owing to the difficulties associated with setting up a simulation and the long computation time required. To take advantage of radiation transport Monte Carlo codes applied to routine clinical practice, researchers and private companies have developed treatment planning and dose verification systems that are partly or fully based on fast Monte Carlo algorithms. This review presents a comprehensive list of the currently existing Monte Carlo systems that can be used to calculate or verify an external photon and electron beam radiotherapy treatment plan. Particular attention is given to those systems that are distributed, either freely or commercially, and that do not require programming tasks from the end user. These systems are compared in terms of features and the simulation time required to compute a set of benchmark calculations. (orig.) [German] Seit mehreren Jahrzehnten werden allgemein anwendbare Monte-Carlo-Codes zur Simulation des Strahlungstransports benutzt, um die Verteilung der absorbierten Dosis in der perkutanen Strahlentherapie mit Photonen und Elektronen zu evaluieren. Die damit erzielten Ergebnisse sind meist akkurater als solche, die mit nichtstochastischen Methoden herkoemmlicher Bestrahlungsplanungssysteme erzielt werden koennen. Wegen des damit verbundenen Arbeitsaufwands und der langen Dauer der Berechnungen wurden Monte-Carlo-Simulationen von Dosisverteilungen in der konventionellen Strahlentherapie in der Vergangenheit im Wesentlichen in der Forschung eingesetzt. Im Bemuehen, Monte

Avances en la regulación y control del estado térmico del horno alto

Directory of Open Access Journals (Sweden)

Babich, A.

1997-04-01

Full Text Available An automated system of direction (SAD for the control of the thermal state of the hearth of the blast furnace is used, taking into account the interrelation existing between the pig iron temperature and its silicon content. By the magnitude of the thermal exergy of the process, which determines the heat work capacity the variations in the heating of the hearth can be valued and carried out the necessary control and regulation actions, such as the variations in the pulverized coal injection and in the sinter/coke relation in the melting bed of the blast furnace.

Se utiliza un sistema automatizado de dirección (SAD para el control del estado térmico del crisol del horno alto, teniendo en cuenta la relación que existe entre la temperatura del arrabio y su contenido de silicio. Por la magnitud de la exergía térmica del proceso, que determina la capacidad de trabajo del calor, se pueden valorar las variaciones de calentamiento del crisol y realizar las acciones necesarias de control y regulación, como son las variaciones en la inyección de carbón pulverizado y en la relación sinterizado/coque en el lecho de fusión del horno alto.

Multilevel Monte Carlo in Approximate Bayesian Computation

KAUST Repository

Jasra, Ajay

2017-02-13

In the following article we consider approximate Bayesian computation (ABC) inference. We introduce a method for numerically approximating ABC posteriors using the multilevel Monte Carlo (MLMC). A sequential Monte Carlo version of the approach is developed and it is shown under some assumptions that for a given level of mean square error, this method for ABC has a lower cost than i.i.d. sampling from the most accurate ABC approximation. Several numerical examples are given.

Monte Carlo simulation of Markov unreliability models

International Nuclear Information System (INIS)

Lewis, E.E.; Boehm, F.

1984-01-01

A Monte Carlo method is formulated for the evaluation of the unrealibility of complex systems with known component failure and repair rates. The formulation is in terms of a Markov process allowing dependences between components to be modeled and computational efficiencies to be achieved in the Monte Carlo simulation. Two variance reduction techniques, forced transition and failure biasing, are employed to increase computational efficiency of the random walk procedure. For an example problem these result in improved computational efficiency by more than three orders of magnitudes over analog Monte Carlo. The method is generalized to treat problems with distributed failure and repair rate data, and a batching technique is introduced and shown to result in substantial increases in computational efficiency for an example problem. A method for separating the variance due to the data uncertainty from that due to the finite number of random walks is presented. (orig.)

A residual Monte Carlo method for discrete thermal radiative diffusion

International Nuclear Information System (INIS)

Evans, T.M.; Urbatsch, T.J.; Lichtenstein, H.; Morel, J.E.

2003-01-01

Residual Monte Carlo methods reduce statistical error at a rate of exp(-bN), where b is a positive constant and N is the number of particle histories. Contrast this convergence rate with 1/√N, which is the rate of statistical error reduction for conventional Monte Carlo methods. Thus, residual Monte Carlo methods hold great promise for increased efficiency relative to conventional Monte Carlo methods. Previous research has shown that the application of residual Monte Carlo methods to the solution of continuum equations, such as the radiation transport equation, is problematic for all but the simplest of cases. However, the residual method readily applies to discrete systems as long as those systems are monotone, i.e., they produce positive solutions given positive sources. We develop a residual Monte Carlo method for solving a discrete 1D non-linear thermal radiative equilibrium diffusion equation, and we compare its performance with that of the discrete conventional Monte Carlo method upon which it is based. We find that the residual method provides efficiency gains of many orders of magnitude. Part of the residual gain is due to the fact that we begin each timestep with an initial guess equal to the solution from the previous timestep. Moreover, fully consistent non-linear solutions can be obtained in a reasonable amount of time because of the effective lack of statistical noise. We conclude that the residual approach has great potential and that further research into such methods should be pursued for more general discrete and continuum systems

O USO DO BIOFERTILIZANTE EM SOLOS AGRÍCOLAS DO CERRADO DA REGIÃO DO ALTO PARANAÍBA (MG

Directory of Open Access Journals (Sweden)

Rodrigo Cavalcanti do Nascimento

2010-12-01

Full Text Available Since the implementation of the integrated project in pork meat production, the region of Alto Paranaiba (MGfaced a new sector restructuring, which resulted in the installation of several agribusiness companies. However,the region is experiencing a strong environmental impact caused by the accumulation of pig slurry and themisuse of technological mechanisms that attempt to correct and fertilize soils without proper accountability.Aware of the various environmental concerns on a global scale that have been highlighted in international environmentalconferences (e.g., Kyoto Protocol, transnational companies, in partnership with pork meat producersfrom the region, began to invest in anaerobic biodigesters, which allow to reduce the organic load by 96%, providingelectrical energy (biogas and producing biofertilizer. From the results of research done by Embrapa, theapplication of biofertilizer on the coffee and corn plantation has brought gains in productivity and plague elimination.Analyzing the soil after application, has found harmful high chemical concentrations to the environment.Therefore, this article aims to analyze the positive and negative aspects. This article aims to analyze the use ofbiofertilizers in agricultural areas of the Alto Paranaiba (MG region, specifically in the municipalities of Patosde Minas and Patrocinio.

Spatial and temporal changes in land cover and its relation to the wildfire risk in Montes de Zuera (NE-Spain)

Science.gov (United States)

Iñiguez-Remón, E.; Badía-Villas, D.; Ibarra-Benlloch, P.

2012-04-01

Land cover changes in a territory can be as result of the diverse human activities, and also, by the dynamic natural ecosystem. The analysis of these changes constitutes a fundamental indicator in improving the knowledge towards a potential sustainable development (OSE, 2006). This paper analyzes the evolution of the land cover at detailed scale 1:5000, along the years 1957, 1984, and 2005, in the Montes de Zuera (UTM: 665000-680000 East /4640000-4650000 North), nearby Zaragoza city (NE-Spain). Despite they are located in the midst of the semiarid Central Ebro Basin, the area have an uncommon, extensive and dense forest cover. For these reasons, Montes de Zuera form part of the "Natura 2000 Network" of Aragon. Because this region has a high incidence of wildfires and, as a measure to prevent them, land cover changes have been analyzed in this study. To get this objective, a methodology is developed based on GIS and Remote Sensing tools, and this makes it possible to establish the relationship between the dynamic land cover and "fire hazard". According to FAO (1986), fire hazard is the measure of risk of wildfire, which is explained by the presence of combustible materials available to burn (Chuvieco et al, 2004). It is understood that the areas with major fire hazard present a high vulnerability to suffer vast wildfires. To develop this work, the type of vegetation, its distribution and the state of maturity of the forest masses are analyzed. Moreover, the wildfires that occurred during the last five decades have been located, and the forest masses and the processes (anthropic and natural) associated to them have been characterized in the regions affected by vast wildfires. Moreover, the cartography of fire hazard is obtained, and this is useful tool to fre prevention and to promote the sustainable management of the forest masses. Two regions with high rate of fire hazards are differentiated (region 1 and region 2). The Aleppo pine cover occupies the majority of

Stochastic approximation Monte Carlo importance sampling for approximating exact conditional probabilities

KAUST Repository

Cheon, Sooyoung

2013-02-16

Importance sampling and Markov chain Monte Carlo methods have been used in exact inference for contingency tables for a long time, however, their performances are not always very satisfactory. In this paper, we propose a stochastic approximation Monte Carlo importance sampling (SAMCIS) method for tackling this problem. SAMCIS is a combination of adaptive Markov chain Monte Carlo and importance sampling, which employs the stochastic approximation Monte Carlo algorithm (Liang et al., J. Am. Stat. Assoc., 102(477):305-320, 2007) to draw samples from an enlarged reference set with a known Markov basis. Compared to the existing importance sampling and Markov chain Monte Carlo methods, SAMCIS has a few advantages, such as fast convergence, ergodicity, and the ability to achieve a desired proportion of valid tables. The numerical results indicate that SAMCIS can outperform the existing importance sampling and Markov chain Monte Carlo methods: It can produce much more accurate estimates in much shorter CPU time than the existing methods, especially for the tables with high degrees of freedom. © 2013 Springer Science+Business Media New York.

Stochastic approximation Monte Carlo importance sampling for approximating exact conditional probabilities

KAUST Repository

Cheon, Sooyoung; Liang, Faming; Chen, Yuguo; Yu, Kai

2013-01-01

Importance sampling and Markov chain Monte Carlo methods have been used in exact inference for contingency tables for a long time, however, their performances are not always very satisfactory. In this paper, we propose a stochastic approximation Monte Carlo importance sampling (SAMCIS) method for tackling this problem. SAMCIS is a combination of adaptive Markov chain Monte Carlo and importance sampling, which employs the stochastic approximation Monte Carlo algorithm (Liang et al., J. Am. Stat. Assoc., 102(477):305-320, 2007) to draw samples from an enlarged reference set with a known Markov basis. Compared to the existing importance sampling and Markov chain Monte Carlo methods, SAMCIS has a few advantages, such as fast convergence, ergodicity, and the ability to achieve a desired proportion of valid tables. The numerical results indicate that SAMCIS can outperform the existing importance sampling and Markov chain Monte Carlo methods: It can produce much more accurate estimates in much shorter CPU time than the existing methods, especially for the tables with high degrees of freedom. © 2013 Springer Science+Business Media New York.

Estimating forest conversion rates with annual forest inventory data

Science.gov (United States)

Paul C. Van Deusen; Francis A. Roesch

2009-01-01

The rate of land-use conversion from forest to nonforest or natural forest to forest plantation is of interest for forest certification purposes and also as part of the process of assessing forest sustainability. Conversion rates can be estimated from remeasured inventory plots in general, but the emphasis here is on annual inventory data. A new estimator is proposed...

Contributon Monte Carlo

International Nuclear Information System (INIS)

Dubi, A.; Gerstl, S.A.W.

1979-05-01

The contributon Monte Carlo method is based on a new recipe to calculate target responses by means of volume integral of the contributon current in a region between the source and the detector. A comprehensive description of the method, its implementation in the general-purpose MCNP code, and results of the method for realistic nonhomogeneous, energy-dependent problems are presented. 23 figures, 10 tables

Bayesian Monte Carlo method

International Nuclear Information System (INIS)

Rajabalinejad, M.

2010-01-01

To reduce cost of Monte Carlo (MC) simulations for time-consuming processes, Bayesian Monte Carlo (BMC) is introduced in this paper. The BMC method reduces number of realizations in MC according to the desired accuracy level. BMC also provides a possibility of considering more priors. In other words, different priors can be integrated into one model by using BMC to further reduce cost of simulations. This study suggests speeding up the simulation process by considering the logical dependence of neighboring points as prior information. This information is used in the BMC method to produce a predictive tool through the simulation process. The general methodology and algorithm of BMC method are presented in this paper. The BMC method is applied to the simplified break water model as well as the finite element model of 17th Street Canal in New Orleans, and the results are compared with the MC and Dynamic Bounds methods.

Closed-shell variational quantum Monte Carlo simulation for the ...

African Journals Online (AJOL)

Closed-shell variational quantum Monte Carlo simulation for the electric dipole moment calculation of hydrazine molecule using casino-code. ... Nigeria Journal of Pure and Applied Physics ... The variational quantum Monte Carlo (VQMC) technique used in this work employed the restricted Hartree-Fock (RHF) scheme.

New Approaches and Applications for Monte Carlo Perturbation Theory

Energy Technology Data Exchange (ETDEWEB)

Aufiero, Manuele; Bidaud, Adrien; Kotlyar, Dan; Leppänen, Jaakko; Palmiotti, Giuseppe; Salvatores, Massimo; Sen, Sonat; Shwageraus, Eugene; Fratoni, Massimiliano

2017-02-01

This paper presents some of the recent and new advancements in the extension of Monte Carlo Perturbation Theory methodologies and application. In particular, the discussed problems involve Brunup calculation, perturbation calculation based on continuous energy functions, and Monte Carlo Perturbation Theory in loosely coupled systems.

Recommender engine for continuous-time quantum Monte Carlo methods

Science.gov (United States)

Huang, Li; Yang, Yi-feng; Wang, Lei

2017-03-01

Recommender systems play an essential role in the modern business world. They recommend favorable items such as books, movies, and search queries to users based on their past preferences. Applying similar ideas and techniques to Monte Carlo simulations of physical systems boosts their efficiency without sacrificing accuracy. Exploiting the quantum to classical mapping inherent in the continuous-time quantum Monte Carlo methods, we construct a classical molecular gas model to reproduce the quantum distributions. We then utilize powerful molecular simulation techniques to propose efficient quantum Monte Carlo updates. The recommender engine approach provides a general way to speed up the quantum impurity solvers.

Rapid Monte Carlo Simulation of Gravitational Wave Galaxies

Science.gov (United States)

Breivik, Katelyn; Larson, Shane L.

2015-01-01

With the detection of gravitational waves on the horizon, astrophysical catalogs produced by gravitational wave observatories can be used to characterize the populations of sources and validate different galactic population models. Efforts to simulate gravitational wave catalogs and source populations generally focus on population synthesis models that require extensive time and computational power to produce a single simulated galaxy. Monte Carlo simulations of gravitational wave source populations can also be used to generate observation catalogs from the gravitational wave source population. Monte Carlo simulations have the advantes of flexibility and speed, enabling rapid galactic realizations as a function of galactic binary parameters with less time and compuational resources required. We present a Monte Carlo method for rapid galactic simulations of gravitational wave binary populations.

Acceleration of monte Carlo solution by conjugate gradient method

International Nuclear Information System (INIS)

Toshihisa, Yamamoto

2005-01-01

The conjugate gradient method (CG) was applied to accelerate Monte Carlo solutions in fixed source problems. The equilibrium model based formulation enables to use CG scheme as well as initial guess to maximize computational performance. This method is available to arbitrary geometry provided that the neutron source distribution in each subregion can be regarded as flat. Even if it is not the case, the method can still be used as a powerful tool to provide an initial guess very close to the converged solution. The major difference of Monte Carlo CG to deterministic CG is that residual error is estimated using Monte Carlo sampling, thus statistical error exists in the residual. This leads to a flow diagram specific to Monte Carlo-CG. Three pre-conditioners were proposed for CG scheme and the performance was compared with a simple 1-D slab heterogeneous test problem. One of them, Sparse-M option, showed an excellent performance in convergence. The performance per unit cost was improved by four times in the test problem. Although direct estimation of efficiency of the method is impossible mainly because of the strong problem-dependence of the optimized pre-conditioner in CG, the method seems to have efficient potential as a fast solution algorithm for Monte Carlo calculations. (author)

PERHITUNGAN VaR PORTOFOLIO SAHAM MENGGUNAKAN DATA HISTORIS DAN DATA SIMULASI MONTE CARLO

Directory of Open Access Journals (Sweden)

WAYAN ARTHINI

2012-09-01

Full Text Available Value at Risk (VaR is the maximum potential loss on a portfolio based on the probability at a certain time.  In this research, portfolio VaR values calculated from historical data and Monte Carlo simulation data. Historical data is processed so as to obtain stock returns, variance, correlation coefficient, and variance-covariance matrix, then the method of Markowitz sought proportion of each stock fund, and portfolio risk and return portfolio. The data was then simulated by Monte Carlo simulation, Exact Monte Carlo Simulation and Expected Monte Carlo Simulation. Exact Monte Carlo simulation have same returns and standard deviation  with historical data, while the Expected Monte Carlo Simulation satistic calculation similar to historical data. The results of this research is the portfolio VaR  with time horizon T=1, T=10, T=22 and the confidence level of 95 %, values obtained VaR between historical data and Monte Carlo simulation data with the method exact and expected. Value of VaR from both Monte Carlo simulation is greater than VaR historical data.

Textiles para turistas: tejedoras y comerciantes en los Altos de Chiapas

Directory of Open Access Journals (Sweden)

Eugenia Bayona Escat

2013-01-01

Full Text Available El presente trabajo indaga en el papel que ocupan las mujeres indígenas como productoras y vendedoras de piezas textiles, consideradas artesanales y de tradición indígena, con alta demanda en el mercado turístico de la región de los Altos de Chiapas, México. La presencia del turismo en la zona ha provocado la consolidación de un complejo mercado local de textiles en el que se entrelazan tejedoras, distribuidoras y comerciantes, cadenas mercantiles y múltiples espacios de venta. A través de este protagonismo femenino y de la comercialización de lo “indígena”, el trabajo indaga en un contexto local que permite a las mujeres indígenas tomar nuevas posiciones sociales a la vez que integra aspectos globales sobre mercancías e imaginarios sociales.

Monte Carlo methods for the reliability analysis of Markov systems

International Nuclear Information System (INIS)

Buslik, A.J.

1985-01-01

This paper presents Monte Carlo methods for the reliability analysis of Markov systems. Markov models are useful in treating dependencies between components. The present paper shows how the adjoint Monte Carlo method for the continuous time Markov process can be derived from the method for the discrete-time Markov process by a limiting process. The straightforward extensions to the treatment of mean unavailability (over a time interval) are given. System unavailabilities can also be estimated; this is done by making the system failed states absorbing, and not permitting repair from them. A forward Monte Carlo method is presented in which the weighting functions are related to the adjoint function. In particular, if the exact adjoint function is known then weighting factors can be constructed such that the exact answer can be obtained with a single Monte Carlo trial. Of course, if the exact adjoint function is known, there is no need to perform the Monte Carlo calculation. However, the formulation is useful since it gives insight into choices of the weight factors which will reduce the variance of the estimator

A general transform for variance reduction in Monte Carlo simulations

International Nuclear Information System (INIS)

Becker, T.L.; Larsen, E.W.

2011-01-01

This paper describes a general transform to reduce the variance of the Monte Carlo estimate of some desired solution, such as flux or biological dose. This transform implicitly includes many standard variance reduction techniques, including source biasing, collision biasing, the exponential transform for path-length stretching, and weight windows. Rather than optimizing each of these techniques separately or choosing semi-empirical biasing parameters based on the experience of a seasoned Monte Carlo practitioner, this General Transform unites all these variance techniques to achieve one objective: a distribution of Monte Carlo particles that attempts to optimize the desired solution. Specifically, this transform allows Monte Carlo particles to be distributed according to the user's specification by using information obtained from a computationally inexpensive deterministic simulation of the problem. For this reason, we consider the General Transform to be a hybrid Monte Carlo/Deterministic method. The numerical results con rm that the General Transform distributes particles according to the user-specified distribution and generally provide reasonable results for shielding applications. (author)

Forest Classification Based on Forest texture in Northwest Yunnan Province

Science.gov (United States)

Wang, Jinliang; Gao, Yan; Wang, Xiaohua; Fu, Lei

2014-03-01

Forest texture is an intrinsic characteristic and an important visual feature of a forest ecological system. Full utilization of forest texture will be a great help in increasing the accuracy of forest classification based on remote sensed data. Taking Shangri-La as a study area, forest classification has been based on the texture. The results show that: (1) From the texture abundance, texture boundary, entropy as well as visual interpretation, the combination of Grayscale-gradient co-occurrence matrix and wavelet transformation is much better than either one of both ways of forest texture information extraction; (2) During the forest texture information extraction, the size of the texture-suitable window determined by the semi-variogram method depends on the forest type (evergreen broadleaf forest is 3×3, deciduous broadleaf forest is 5×5, etc.). (3)While classifying forest based on forest texture information, the texture factor assembly differs among forests: Variance Heterogeneity and Correlation should be selected when the window is between 3×3 and 5×5 Mean, Correlation, and Entropy should be used when the window in the range of 7×7 to 19×19 and Correlation, Second Moment, and Variance should be used when the range is larger than 21×21.

Forest Classification Based on Forest texture in Northwest Yunnan Province

International Nuclear Information System (INIS)

Wang, Jinliang; Gao, Yan; Fu, Lei; Wang, Xiaohua

2014-01-01

Forest texture is an intrinsic characteristic and an important visual feature of a forest ecological system. Full utilization of forest texture will be a great help in increasing the accuracy of forest classification based on remote sensed data. Taking Shangri-La as a study area, forest classification has been based on the texture. The results show that: (1) From the texture abundance, texture boundary, entropy as well as visual interpretation, the combination of Grayscale-gradient co-occurrence matrix and wavelet transformation is much better than either one of both ways of forest texture information extraction; (2) During the forest texture information extraction, the size of the texture-suitable window determined by the semi-variogram method depends on the forest type (evergreen broadleaf forest is 3×3, deciduous broadleaf forest is 5×5, etc.). (3)While classifying forest based on forest texture information, the texture factor assembly differs among forests: Variance Heterogeneity and Correlation should be selected when the window is between 3×3 and 5×5; Mean, Correlation, and Entropy should be used when the window in the range of 7×7 to 19×19; and Correlation, Second Moment, and Variance should be used when the range is larger than 21×21

A Monte Carlo approach to combating delayed completion of ...

African Journals Online (AJOL)

The objective of this paper is to unveil the relevance of Monte Carlo critical path analysis in resolving problem of delays in scheduled completion of development projects. Commencing with deterministic network scheduling, Monte Carlo critical path analysis was advanced by assigning probability distributions to task times.

Simulation of Forest Cover Dynamics for Eastern Eurasian Boreal Forests

Science.gov (United States)

Shugart, H. H.; Yan, X.; Zhang, N.; Isaev, A. S.; Shuman, J. K.

2006-12-01

We are developing and testing a boreal zone forest dynamics model capable of simulating the forest cover dynamics of the Eurasian boreal forest, a major biospheric ecosystem with potentially large roles in the planetary carbon cycle and in the feedback between terrestrial surface and the atmosphere. In appreciating the role of this region in the coupling between atmosphere and terrestrial surface, on must understand the interactions between CO2 source/sink relationships (associated with growing or clearing forests) and the albedo effects (from changes in terrestrial surface cover). There is some evidence that in the Eurasian Boreal zone, the Carbon budget effects from forest change may oppose the albedo changes. This creates complex feedbacks between surface and atmosphere and motivates the need for a forest dynamics model that simultaneous represents forest vegetation and carbon storage and release. A forest dynamics model applied to Eastern Eurasia, FAREAST, has been tested using three types of information: 1. Direct species composition comparisons between simulated and observed mature forests at the same locations; 2. Forest type comparisons between simulated and observed forests along altitudinal gradients of several different mountains; 3. Comparison with forest stands in different succession stages of simulated forests. Model comparisons with independent data indicate the FAREAST model is capable of representing many of the broad features of the forests of Northeastern China. After model validation in the Northeast China region, model applications were developed for the forests of the Russian Far East. Continental-scale forest cover can be simulated to a relatively realistic degree using a forest gap model with standard representations of individual-plant processes. It appears that such a model, validated relatively locally in this case, in Northeastern China, can then be applied over a much larger region and under conditions of climatic change.

Perturbation based Monte Carlo criticality search in density, enrichment and concentration

International Nuclear Information System (INIS)

Li, Zeguang; Wang, Kan; Deng, Jingkang

2015-01-01

Highlights: • A new perturbation based Monte Carlo criticality search method is proposed. • The method could get accurate results with only one individual criticality run. • The method is used to solve density, enrichment and concentration search problems. • Results show the feasibility and good performances of this method. • The relationship between results’ accuracy and perturbation order is discussed. - Abstract: Criticality search is a very important aspect in reactor physics analysis. Due to the advantages of Monte Carlo method and the development of computer technologies, Monte Carlo criticality search is becoming more and more necessary and feasible. Existing Monte Carlo criticality search methods need large amount of individual criticality runs and may have unstable results because of the uncertainties of criticality results. In this paper, a new perturbation based Monte Carlo criticality search method is proposed and discussed. This method only needs one individual criticality calculation with perturbation tallies to estimate k eff changing function using initial k eff and differential coefficients results, and solves polynomial equations to get the criticality search results. The new perturbation based Monte Carlo criticality search method is implemented in the Monte Carlo code RMC, and criticality search problems in density, enrichment and concentration are taken out. Results show that this method is quite inspiring in accuracy and efficiency, and has advantages compared with other criticality search methods

European mixed forests

DEFF Research Database (Denmark)

Bravo-Oviedo, Andres; Pretzsch, Hans; Ammer, Christian

2014-01-01

Aim of study: We aim at (i) developing a reference definition of mixed forests in order to harmonize comparative research in mixed forests and (ii) review the research perspectives in mixed forests. Area of study: The definition is developed in Europe but can be tested worldwide. Material...... and Methods: Review of existent definitions of mixed forests based and literature review encompassing dynamics, management and economic valuation of mixed forests. Main results: A mixed forest is defined as a forest unit, excluding linear formations, where at least two tree species coexist at any...... density in mixed forests, (iii) conversion of monocultures to mixed-species forest and (iv) economic valuation of ecosystem services provided by mixed forests. Research highlights: The definition is considered a high-level one which encompasses previous attempts to define mixed forests. Current fields...

PENENTUAN HARGA OPSI BELI TIPE ASIA DENGAN METODE MONTE CARLO-CONTROL VARIATE

Directory of Open Access Journals (Sweden)

NI NYOMAN AYU ARTANADI

2017-01-01

Full Text Available Option is a contract between the writer and the holder which entitles the holder to buy or sell an underlying asset at the maturity date for a specified price known as an exercise price. Asian option is a type of financial derivatives which the payoff taking the average value over the time series of the asset price. The aim of the study is to present the Monte Carlo-Control Variate as an extension of Standard Monte Carlo applied on the calculation of the Asian option price. Standard Monte Carlo simulations 10.000.000 generate standard error 0.06 and the option price convergent at Rp.160.00 while Monte Carlo-Control Variate simulations 100.000 generate standard error 0.01 and the option price convergent at Rp.152.00. This shows the Monte Carlo-Control Variate achieve faster option price toward convergent of the Monte Carlo Standar.

Monte Carlo numerical study of lattice field theories

International Nuclear Information System (INIS)

Gan Cheekwan; Kim Seyong; Ohta, Shigemi

1997-01-01

The authors are interested in the exact first-principle calculations of quantum field theories which are indeed exact ones. For quantum chromodynamics (QCD) at low energy scale, a nonperturbation method is needed, and the only known such method is the lattice method. The path integral can be evaluated by putting a system on a finite 4-dimensional volume and discretizing space time continuum into finite points, lattice. The continuum limit is taken by making the lattice infinitely fine. For evaluating such a finite-dimensional integral, the Monte Carlo numerical estimation of the path integral can be obtained. The calculation of light hadron mass in quenched lattice QCD with staggered quarks, 3-dimensional Thirring model calculation and the development of self-test Monte Carlo method have been carried out by using the RIKEN supercomputer. The motivation of this study, lattice QCD formulation, continuum limit, Monte Carlo update, hadron propagator, light hadron mass, auto-correlation and source size dependence are described on lattice QCD. The phase structure of the 3-dimensional Thirring model for a small 8 3 lattice has been mapped. The discussion on self-test Monte Carlo method is described again. (K.I.)

Continuous energy Monte Carlo method based lattice homogeinzation

International Nuclear Information System (INIS)

Li Mancang; Yao Dong; Wang Kan

2014-01-01

Based on the Monte Carlo code MCNP, the continuous energy Monte Carlo multi-group constants generation code MCMC has been developed. The track length scheme has been used as the foundation of cross section generation. The scattering matrix and Legendre components require special techniques, and the scattering event method has been proposed to solve this problem. Three methods have been developed to calculate the diffusion coefficients for diffusion reactor core codes and the Legendre method has been applied in MCMC. To the satisfaction of the equivalence theory, the general equivalence theory (GET) and the superhomogenization method (SPH) have been applied to the Monte Carlo method based group constants. The super equivalence method (SPE) has been proposed to improve the equivalence. GET, SPH and SPE have been implemented into MCMC. The numerical results showed that generating the homogenization multi-group constants via Monte Carlo method overcomes the difficulties in geometry and treats energy in continuum, thus provides more accuracy parameters. Besides, the same code and data library can be used for a wide range of applications due to the versatility. The MCMC scheme can be seen as a potential alternative to the widely used deterministic lattice codes. (authors)

Sentidos de Vitória/Derrota para os Pais Segundo Atletas do Alto Rendimento

Directory of Open Access Journals (Sweden)

Isabela Amblard

Full Text Available Este estudo compreendeu as representações sociais da vitória/derrota para os pais segundo atletas-adolescentes do esporte de alto rendimento, na cidade do Recife. A adolescência é compreendida a partir da Psicologia Social-histórica, e o embasamento teórico-metodológico da Teoria das Representações Sociais abordou os sujeitos em diferentes contextos socioculturais, lugares de pertencimento, experiências, crenças, saberes e sentimentos compartilhados. Adotou-se a perspectiva pluri metodológica com variados recursos de coleta e análise progressiva de dados. Participaram 101 atletas-adolescentes do esporte de alto rendimento, nas modalidades natação e vôlei. Os instrumentos utilizados foram questionários de associação livre e entrevistas semidirigidas. Na análise dos dados, o software EVOC e as técnicas de análise temática de conteúdo de Bardin. Identificamos as representações de vitória para os pais, nas dimensões: pessoal e motivacional, e afetivo-emocional e, nas representações sociais da derrota para os pais, além destas, a dimensão técnica. O sentido de vitória para os pais apareceu ampliado: o contexto esportivo, a escolarização, e o crescimento pessoal e profissional na vida. A derrota para os pais é representada como oportunidade de aprendizado e superação para o atleta-adolescente, porém, eles mostraram que necessitam do apoio afetivo de sua família para lidar com a autoculpabilização nas situações de fracasso.

Implications of Monte Carlo Statistical Errors in Criticality Safety Assessments

International Nuclear Information System (INIS)

Pevey, Ronald E.

2005-01-01

Most criticality safety calculations are performed using Monte Carlo techniques because of Monte Carlo's ability to handle complex three-dimensional geometries. For Monte Carlo calculations, the more histories sampled, the lower the standard deviation of the resulting estimates. The common intuition is, therefore, that the more histories, the better; as a result, analysts tend to run Monte Carlo analyses as long as possible (or at least to a minimum acceptable uncertainty). For Monte Carlo criticality safety analyses, however, the optimization situation is complicated by the fact that procedures usually require that an extra margin of safety be added because of the statistical uncertainty of the Monte Carlo calculations. This additional safety margin affects the impact of the choice of the calculational standard deviation, both on production and on safety. This paper shows that, under the assumptions of normally distributed benchmarking calculational errors and exact compliance with the upper subcritical limit (USL), the standard deviation that optimizes production is zero, but there is a non-zero value of the calculational standard deviation that minimizes the risk of inadvertently labeling a supercritical configuration as subcritical. Furthermore, this value is shown to be a simple function of the typical benchmarking step outcomes--the bias, the standard deviation of the bias, the upper subcritical limit, and the number of standard deviations added to calculated k-effectives before comparison to the USL

Biased Monte Carlo optimization: the basic approach

International Nuclear Information System (INIS)

Campioni, Luca; Scardovelli, Ruben; Vestrucci, Paolo

2005-01-01

It is well-known that the Monte Carlo method is very successful in tackling several kinds of system simulations. It often happens that one has to deal with rare events, and the use of a variance reduction technique is almost mandatory, in order to have Monte Carlo efficient applications. The main issue associated with variance reduction techniques is related to the choice of the value of the biasing parameter. Actually, this task is typically left to the experience of the Monte Carlo user, who has to make many attempts before achieving an advantageous biasing. A valuable result is provided: a methodology and a practical rule addressed to establish an a priori guidance for the choice of the optimal value of the biasing parameter. This result, which has been obtained for a single component system, has the notable property of being valid for any multicomponent system. In particular, in this paper, the exponential and the uniform biases of exponentially distributed phenomena are investigated thoroughly

PRODUCCIÓN DE HOJARASCA FINA EN BOSQUES ALTO ANDINOS DE ANTIOQUIA, COLOMBIA FINE LITTER PRODUCTION IN HIGH ANDEAN FORESTS FROM ANTIOQUIA, COLOMBIA

Directory of Open Access Journals (Sweden)

Claudia Marcela Zapata Duque

2007-06-01

Full Text Available En bosques montanos naturales de Quercus humboldtii y reforestados (Pinus patula y Cupressus lusitanica de la región de Piedras Blancas, Antioquia (Colombia, fue evaluada la producción de hojarasca por un periodo de 2 años. Se utilizaron trampas de hojarasca con el fin de recoger el material desprendido del dosel para su posterior separación en fracciones y pesado respectivo. El promedio de caída de hojarasca anual para Q. humboldtii, P. patula y C. lusitanica fue de 7877,20; 8362,47 y 3725,97 kg ha-1año-1 respectivamente; siendo la fracción foliar la que mayor participación tuvo en la producción total. Mediante análisis de regresión lineal múltiple se ajustaron modelos de producción de hojarasca según fracciones por cobertura en función de diferentes variables hidrológicas, tales como la intensidad y la cantidad de lluvia del período simultáneo a la recolección de la hojarasca o inmediatamente anterior.Litter production was measured over two years in a montane oak forest (Quercus humboldtii and in pine (Pinus patula and cypress (Cupressus lusitanica plantations in Piedras Blancas, Antioquia ( Colombia . Litter traps were used in order to collect litterfall to be subsequently separated into fractions and weighed. Annual mean litterfall for Q. humboldtii, P. patula and C. lusitanica was of 7877,20; 8362,47 and 3725,97 kg ha-1year-1 respectively; being the leaf fraction of highest participation in total production. Multiple linear regression models were used to fit litter production for each fraction and forest cover as a function of different hydrological variables such as intensity and quantity of precipitation, both during the period when the leaf litter was collected and immediately preceding one.

Forest crimes as a threat to sustainable forest management

Directory of Open Access Journals (Sweden)

S. Özden

2016-08-01

Full Text Available From ancient times to the present day, forest public relations has been an issue on the agenda. This relationship’s purpose was initially needed for shelter and nutrition; however today this process has changed with urbanization, overpopulation and understanding the new functions of forests. When land ownership became a tool of production, offenses occurred in order to convert forestlands to agricultural lands. So the vast majority of the world’s forests have been lost for this reason. Today, deforestation is occurring in tropical countries that are expecting to gain agricultural area. The purpose of this study is to investigate the relationship between urbanization and the qualitative and quantitative characteristics of forest crimes, which are a major obstacle for sustainable forestry. Although forests cover about 27 % of Turkey’s territory, the forests are losing viability; the status of wood raw material per unit area and the total area of the country in the ratio of productive forests are becoming critical in Turkey. Turkey’s rugged terrain and factors such as human interventions, fires, deforestation for agriculture, illegal cuttings, or improper grazing reduce existing forests or cause deterioration of their structure. In the past, deforestation, as a result of human interventions in Turkey, was done by forest villagers who live in rural areas. The forest crimes depend on various socio-economic reasons and have many adverse effects on the sustainability of forest and forest existence. In developed countries, illegal interventions such as opening, grazing, cutting, occupation, use, settlement, or hunting crimes have been largely eliminated because of the absence of cadastral problems, the existence of more responsive people to protect the environment and forests and a rural population, which has a higher standard of living. In the last 20 years, there has been both a dramatic decrease in the population living in rural areas and a

Self-learning Monte Carlo (dynamical biasing)

International Nuclear Information System (INIS)

Matthes, W.

1981-01-01

In many applications the histories of a normal Monte Carlo game rarely reach the target region. An approximate knowledge of the importance (with respect to the target) may be used to guide the particles more frequently into the target region. A Monte Carlo method is presented in which each history contributes to update the importance field such that eventually most target histories are sampled. It is a self-learning method in the sense that the procedure itself: (a) learns which histories are important (reach the target) and increases their probability; (b) reduces the probabilities of unimportant histories; (c) concentrates gradually on the more important target histories. (U.K.)

Tromboembolismo pulmonar masivo de alto riesgo asociado a foramen oval permeable High-risk massive pulmonary thromboembolism associated with patent foramen ovale

Directory of Open Access Journals (Sweden)

Antonio Miranda

2012-04-01

Full Text Available La alta mortalidad de los pacientes con tromboembolismo pulmonar masivo de alto riesgo amerita un enfoque terapéutico enérgico e invasivo que incluya la embolectomía pulmonar quirúrgica en aquellos pacientes con contraindicación para trombolisis o trombolisis fallida. Describimos un caso de tromboembolismo pulmonar masivo de alto riesgo que recibió tratamiento quirúrgico en vez de trombolisis debido a que al momento del diagnóstico presentaba un trombo móvil a través de un foramen oval permeable con altísima posibilidad de embolismo paradójico arterial.High mortality rate associated with massive pulmonary embolism requires an aggressive invasive approach including surgical pulmonary embolectomy when thrombolytic therapy has failed or is contraindicated. We describe a case of high-risk massive pulmonary embolism who underwent surgical treatment due to the presence of a mobile intracardiac clot in a patent foramen ovale, and the possible risk of paradoxical arterial embolism.

RNA folding kinetics using Monte Carlo and Gillespie algorithms.

Science.gov (United States)

Clote, Peter; Bayegan, Amir H

2018-04-01

RNA secondary structure folding kinetics is known to be important for the biological function of certain processes, such as the hok/sok system in E. coli. Although linear algebra provides an exact computational solution of secondary structure folding kinetics with respect to the Turner energy model for tiny ([Formula: see text]20 nt) RNA sequences, the folding kinetics for larger sequences can only be approximated by binning structures into macrostates in a coarse-grained model, or by repeatedly simulating secondary structure folding with either the Monte Carlo algorithm or the Gillespie algorithm. Here we investigate the relation between the Monte Carlo algorithm and the Gillespie algorithm. We prove that asymptotically, the expected time for a K-step trajectory of the Monte Carlo algorithm is equal to [Formula: see text] times that of the Gillespie algorithm, where [Formula: see text] denotes the Boltzmann expected network degree. If the network is regular (i.e. every node has the same degree), then the mean first passage time (MFPT) computed by the Monte Carlo algorithm is equal to MFPT computed by the Gillespie algorithm multiplied by [Formula: see text]; however, this is not true for non-regular networks. In particular, RNA secondary structure folding kinetics, as computed by the Monte Carlo algorithm, is not equal to the folding kinetics, as computed by the Gillespie algorithm, although the mean first passage times are roughly correlated. Simulation software for RNA secondary structure folding according to the Monte Carlo and Gillespie algorithms is publicly available, as is our software to compute the expected degree of the network of secondary structures of a given RNA sequence-see http://bioinformatics.bc.edu/clote/RNAexpNumNbors .

MONTE: the next generation of mission design and navigation software

Science.gov (United States)

Evans, Scott; Taber, William; Drain, Theodore; Smith, Jonathon; Wu, Hsi-Cheng; Guevara, Michelle; Sunseri, Richard; Evans, James

2018-03-01

The Mission analysis, Operations and Navigation Toolkit Environment (MONTE) (Sunseri et al. in NASA Tech Briefs 36(9), 2012) is an astrodynamic toolkit produced by the Mission Design and Navigation Software Group at the Jet Propulsion Laboratory. It provides a single integrated environment for all phases of deep space and Earth orbiting missions. Capabilities include: trajectory optimization and analysis, operational orbit determination, flight path control, and 2D/3D visualization. MONTE is presented to the user as an importable Python language module. This allows a simple but powerful user interface via CLUI or script. In addition, the Python interface allows MONTE to be used seamlessly with other canonical scientific programming tools such as SciPy, NumPy, and Matplotlib. MONTE is the prime operational orbit determination software for all JPL navigated missions.

A NEW MONTE CARLO METHOD FOR TIME-DEPENDENT NEUTRINO RADIATION TRANSPORT

International Nuclear Information System (INIS)

Abdikamalov, Ernazar; Ott, Christian D.; O'Connor, Evan; Burrows, Adam; Dolence, Joshua C.; Löffler, Frank; Schnetter, Erik

2012-01-01

Monte Carlo approaches to radiation transport have several attractive properties such as simplicity of implementation, high accuracy, and good parallel scaling. Moreover, Monte Carlo methods can handle complicated geometries and are relatively easy to extend to multiple spatial dimensions, which makes them potentially interesting in modeling complex multi-dimensional astrophysical phenomena such as core-collapse supernovae. The aim of this paper is to explore Monte Carlo methods for modeling neutrino transport in core-collapse supernovae. We generalize the Implicit Monte Carlo photon transport scheme of Fleck and Cummings and gray discrete-diffusion scheme of Densmore et al. to energy-, time-, and velocity-dependent neutrino transport. Using our 1D spherically-symmetric implementation, we show that, similar to the photon transport case, the implicit scheme enables significantly larger timesteps compared with explicit time discretization, without sacrificing accuracy, while the discrete-diffusion method leads to significant speed-ups at high optical depth. Our results suggest that a combination of spectral, velocity-dependent, Implicit Monte Carlo and discrete-diffusion Monte Carlo methods represents a robust approach for use in neutrino transport calculations in core-collapse supernovae. Our velocity-dependent scheme can easily be adapted to photon transport.

A NEW MONTE CARLO METHOD FOR TIME-DEPENDENT NEUTRINO RADIATION TRANSPORT

Energy Technology Data Exchange (ETDEWEB)

Abdikamalov, Ernazar; Ott, Christian D.; O' Connor, Evan [TAPIR, California Institute of Technology, MC 350-17, 1200 E California Blvd., Pasadena, CA 91125 (United States); Burrows, Adam; Dolence, Joshua C. [Department of Astrophysical Sciences, Princeton University, Peyton Hall, Ivy Lane, Princeton, NJ 08544 (United States); Loeffler, Frank; Schnetter, Erik, E-mail: abdik@tapir.caltech.edu [Center for Computation and Technology, Louisiana State University, 216 Johnston Hall, Baton Rouge, LA 70803 (United States)

2012-08-20

Monte Carlo approaches to radiation transport have several attractive properties such as simplicity of implementation, high accuracy, and good parallel scaling. Moreover, Monte Carlo methods can handle complicated geometries and are relatively easy to extend to multiple spatial dimensions, which makes them potentially interesting in modeling complex multi-dimensional astrophysical phenomena such as core-collapse supernovae. The aim of this paper is to explore Monte Carlo methods for modeling neutrino transport in core-collapse supernovae. We generalize the Implicit Monte Carlo photon transport scheme of Fleck and Cummings and gray discrete-diffusion scheme of Densmore et al. to energy-, time-, and velocity-dependent neutrino transport. Using our 1D spherically-symmetric implementation, we show that, similar to the photon transport case, the implicit scheme enables significantly larger timesteps compared with explicit time discretization, without sacrificing accuracy, while the discrete-diffusion method leads to significant speed-ups at high optical depth. Our results suggest that a combination of spectral, velocity-dependent, Implicit Monte Carlo and discrete-diffusion Monte Carlo methods represents a robust approach for use in neutrino transport calculations in core-collapse supernovae. Our velocity-dependent scheme can easily be adapted to photon transport.

Structural rearrangements detected in newly-formed hexaploid

Indian Academy of Sciences (India)

Institute of Biotechnology and Bioengineering (IBB), Centre of Genomics and Biotechnology (CGB), University of Tras-os-Montes and Alto Douro (UTAD), P.O. Box 1013, 5001-801 VilaReal, Portugal; Instituto de Agricultura Sostenible (IAS), Consejo Superior de Investigaciones Científicas (CSIC), Apartado 4084, 14080 ...

Therapeutic Applications of Monte Carlo Calculations in Nuclear Medicine

International Nuclear Information System (INIS)

Coulot, J

2003-01-01

Monte Carlo techniques are involved in many applications in medical physics, and the field of nuclear medicine has seen a great development in the past ten years due to their wider use. Thus, it is of great interest to look at the state of the art in this domain, when improving computer performances allow one to obtain improved results in a dramatically reduced time. The goal of this book is to make, in 15 chapters, an exhaustive review of the use of Monte Carlo techniques in nuclear medicine, also giving key features which are not necessary directly related to the Monte Carlo method, but mandatory for its practical application. As the book deals with therapeutic' nuclear medicine, it focuses on internal dosimetry. After a general introduction on Monte Carlo techniques and their applications in nuclear medicine (dosimetry, imaging and radiation protection), the authors give an overview of internal dosimetry methods (formalism, mathematical phantoms, quantities of interest). Then, some of the more widely used Monte Carlo codes are described, as well as some treatment planning softwares. Some original techniques are also mentioned, such as dosimetry for boron neutron capture synovectomy. It is generally well written, clearly presented, and very well documented. Each chapter gives an overview of each subject, and it is up to the reader to investigate it further using the extensive bibliography provided. Each topic is discussed from a practical point of view, which is of great help for non-experienced readers. For instance, the chapter about mathematical aspects of Monte Carlo particle transport is very clear and helps one to apprehend the philosophy of the method, which is often a difficulty with a more theoretical approach. Each chapter is put in the general (clinical) context, and this allows the reader to keep in mind the intrinsic limitation of each technique involved in dosimetry (for instance activity quantitation). Nevertheless, there are some minor remarks to

Grain-boundary melting: A Monte Carlo study

DEFF Research Database (Denmark)

Besold, Gerhard; Mouritsen, Ole G.

1994-01-01

Grain-boundary melting in a lattice-gas model of a bicrystal is studied by Monte Carlo simulation using the grand canonical ensemble. Well below the bulk melting temperature T(m), a disordered liquidlike layer gradually emerges at the grain boundary. Complete interfacial wetting can be observed...... when the temperature approaches T(m) from below. Monte Carlo data over an extended temperature range indicate a logarithmic divergence w(T) approximately - ln(T(m)-T) of the width of the disordered layer w, in agreement with mean-field theory....

Modelling mixed forest growth : a review of models for forest management

NARCIS (Netherlands)

Porte, A.; Bartelink, H.H.

2002-01-01

Most forests today are multi-specific and heterogeneous forests (`mixed forests'). However, forest modelling has been focusing on mono-specific stands for a long time, only recently have models been developed for mixed forests. Previous reviews of mixed forest modelling were restricted to certain

Analysis of error in Monte Carlo transport calculations

International Nuclear Information System (INIS)

Booth, T.E.

1979-01-01

The Monte Carlo method for neutron transport calculations suffers, in part, because of the inherent statistical errors associated with the method. Without an estimate of these errors in advance of the calculation, it is difficult to decide what estimator and biasing scheme to use. Recently, integral equations have been derived that, when solved, predicted errors in Monte Carlo calculations in nonmultiplying media. The present work allows error prediction in nonanalog Monte Carlo calculations of multiplying systems, even when supercritical. Nonanalog techniques such as biased kernels, particle splitting, and Russian Roulette are incorporated. Equations derived here allow prediction of how much a specific variance reduction technique reduces the number of histories required, to be weighed against the change in time required for calculation of each history. 1 figure, 1 table

PyMercury: Interactive Python for the Mercury Monte Carlo Particle Transport Code

International Nuclear Information System (INIS)

Iandola, F.N.; O'Brien, M.J.; Procassini, R.J.

2010-01-01

Monte Carlo particle transport applications are often written in low-level languages (C/C++) for optimal performance on clusters and supercomputers. However, this development approach often sacrifices straightforward usability and testing in the interest of fast application performance. To improve usability, some high-performance computing applications employ mixed-language programming with high-level and low-level languages. In this study, we consider the benefits of incorporating an interactive Python interface into a Monte Carlo application. With PyMercury, a new Python extension to the Mercury general-purpose Monte Carlo particle transport code, we improve application usability without diminishing performance. In two case studies, we illustrate how PyMercury improves usability and simplifies testing and validation in a Monte Carlo application. In short, PyMercury demonstrates the value of interactive Python for Monte Carlo particle transport applications. In the future, we expect interactive Python to play an increasingly significant role in Monte Carlo usage and testing.

Neutron flux calculation by means of Monte Carlo methods

International Nuclear Information System (INIS)

Barz, H.U.; Eichhorn, M.

1988-01-01

In this report a survey of modern neutron flux calculation procedures by means of Monte Carlo methods is given. Due to the progress in the development of variance reduction techniques and the improvements of computational techniques this method is of increasing importance. The basic ideas in application of Monte Carlo methods are briefly outlined. In more detail various possibilities of non-analog games and estimation procedures are presented, problems in the field of optimizing the variance reduction techniques are discussed. In the last part some important international Monte Carlo codes and own codes of the authors are listed and special applications are described. (author)

Forest Stakeholder Participation in Improving Game Habitat in Swedish Forests

Directory of Open Access Journals (Sweden)

Eugene E. Ezebilo

2012-07-01

Full Text Available Although in Sweden the simultaneous use of forests for timber production and game hunting are both of socioeconomic importance it often leads to conflicting interests. This study examines forest stakeholder participation in improving game habitat to increase hunting opportunities as well as redistribute game activities in forests to help reduce browsing damage in valuable forest stands. The data for the study were collected from a nationwide survey that involved randomly selected hunters and forest owners in Sweden. An ordered logit model was used to account for possible factors influencing the respondents’ participation in improving game habitat. The results showed that on average, forest owning hunters were more involved in improving game habitat than non-hunting forest owners. The involvement of non-forest owning hunters was intermediate between the former two groups. The respondents’ participation in improving game habitat were mainly influenced by factors such as the quantity of game meat obtained, stakeholder group, forests on hunting grounds, the extent of risk posed by game browsing damage to the economy of forest owners, importance of bagging game during hunting, and number of hunting days. The findings will help in designing a more sustainable forest management strategy that integrates timber production and game hunting in forests.

Effects of Deforestation and Forest Degradation on Forest Carbon Stocks in Collaborative Forests, Nepal

Directory of Open Access Journals (Sweden)

Ram Asheshwar MANDAL

2012-12-01

Full Text Available There are some key drivers that favor deforestation and forest degradation. Consequently, levels of carbon stock are affected in different parts of same forest types. But the problem lies in exploring the extent of the effects on level of carbon stocking. This paper highlights the variations in levels of carbon stocks in three different collaborative forests of same forest type i.e. tropical sal (Shorea robusta forest in Mahottari district of the central Terai in Nepal. Three collaborative forests namely Gadhanta-Bardibas Collaborative Forest (CFM, Tuteshwarnath CFM and Banke- Maraha CFM were selected for research site. Interview and workshops were organized with the key informants that include staffs, members and representatives of CFMs to collect the socio-economic data and stratified random sampling was applied to collect the bio-physical data to calculate the carbon stocks. Analysis was carried out using statistical tools. It was found five major drivers namely grazing, fire, logging, growth of invasive species and encroachment. It was found highest carbon 269.36 ton per ha in Gadhanta- Bardibash CFM. The findings showed that the levels of carbon stocks in the three studied CFMs are different depending on how the drivers of deforestation and forest degradation influence over them.

Estimating national forest carbon stocks and dynamics: combining models and remotely sensed information

Science.gov (United States)

Smallman, Thomas Luke; Exbrayat, Jean-François; Bloom, Anthony; Williams, Mathew

2017-04-01

Forests are a critical component of the global carbon cycle, storing significant amounts of carbon, split between living biomass and dead organic matter. The carbon budget of forests is the most uncertain component of the global carbon cycle - it is currently impossible to quantify accurately the carbon source/sink strength of forest biomes due to their heterogeneity and complex dynamics. It has been a major challenge to generate robust carbon budgets across landscapes due to data scarcity. Models have been used for estimating carbon budgets, but outputs have lacked an assessment of uncertainty, making a robust assessment of their reliability and accuracy challenging. Here a Metropolis Hastings - Markov Chain Monte Carlo (MH-MCMC) data assimilation framework has been used to combine remotely sensed leaf area index (MODIS), biomass (where available) and deforestation estimates, in addition to forest planting information from the UK's national forest inventory, an estimate of soil carbon from the Harmonized World Database (HWSD) and plant trait information with a process model (DALEC) to produce a constrained analysis with a robust estimate of uncertainty of the UK forestry carbon budget between 2000 and 2010. Our analysis estimates the mean annual UK forest carbon sink at -3.9 MgC ha-1 yr-1 with a 95 % confidence interval between -4.0 and -3.1 MgC ha-1yr-1. The UK national forest inventory (NFI) estimates the mean UK forest carbon sink to be between -1.4 and -5.5 MgC ha-1 yr-1. The analysis estimate for total forest biomass stock in 2010 is estimated at 229 (177/232) TgC, while the NFI an estimated total forest biomass carbon stock of 216 TgC. Leaf carbon area (LCA) is a key plant trait which we are able to estimate using our analysis. Comparison of median estimates for (LCA) retrieved from the analysis and a UK land cover map show higher and lower values for LCA are estimated areas dominated by needle leaf and broad leaf forests forest respectively, consistent with

Markov Chain Monte Carlo

Indian Academy of Sciences (India)

Home; Journals; Resonance – Journal of Science Education; Volume 7; Issue 3. Markov Chain Monte Carlo - Examples. Arnab Chakraborty. General Article Volume 7 Issue 3 March 2002 pp 25-34. Fulltext. Click here to view fulltext PDF. Permanent link: https://www.ias.ac.in/article/fulltext/reso/007/03/0025-0034. Keywords.

Forest soils

Science.gov (United States)

Charles H. (Hobie) Perry; Michael C. Amacher

2009-01-01

Productive soils are the foundation of sustainable forests throughout the United States. Forest soils are generally subjected to fewer disturbances than agricultural soils, particularly those that are tilled, so forest soils tend to have better preserved A-horizons than agricultural soils. Another major contrast between forest and agricultural soils is the addition of...

Caracterización de las gestantes de alto riesgo obstétrico (ARO en el departamento de Sucre (Colombia, 2015

Directory of Open Access Journals (Sweden)

Judith Martínez Royert

2016-01-01

Full Text Available Objetivo: Caracterizar las gestantes de Alto Riesgo Obstétrico ( ARO que acuden a una IPS pública en el departamento de Sucre, Colombia (periodo enero, febrero y marzo de 2015. Material y métodos: Estudio cuantitativo, descriptivo. La muestra la conformaron 123 gestantes ARO . Se utilizó un instrumento elaborado por las investigadoras; se sometió a validez de constructo y contenido, y análisis de consistencia interna mediante el alfa de Cronbach. Resultados: El 13,18 % de las gestantes eran menores de 18 años; 38,2% procedentes de la capital y 19,5 % de la región del San Jorge; 66 % no manifestaron antecedentes patológicos; 13,8 % presentaron complicaciones de amenaza de aborto o de parto pretérmino; 37 % eran nulípara; 20.3 % tenían cesárea anterior; 22.8 % sufrieron abortos; 54.5 % manifestaron tensión emocional y mal humor; 82.9% no programaron el embarazo; 24 % con periodo intergenésico de 1 año; 55.3 % (68 gestantes se encontraban entre la semana 30 y 40 de gestación al momento de participar en el estudio. Conclusiones: La subregión de la Sabana y San Jorge fueron las que presentaron mayor número de gestantes de alto riesgo. Entre las patologías preexistentes más frecuentes se encontró anemias y migrañas, así como las del sistema endocrino y respiratorio. Esta investigación servirá como referente para proporcionar conocimiento respecto al perfil de las gestantes de alto riesgo en Sucre, para que los profesionales involucrados en su atención desempeñen un rol que permita contribuir al control y prevención de las complicaciones en ellas y en la reducción significativa de la mortalidad materna.

Monte Carlo studies of high-transverse-energy hadronic interactions

International Nuclear Information System (INIS)

Corcoran, M.D.

1985-01-01

A four-jet Monte Carlo calculation has been used to simulate hadron-hadron interactions which deposit high transverse energy into a large-solid-angle calorimeter and limited solid-angle regions of the calorimeter. The calculation uses first-order QCD cross sections to generate two scattered jets and also produces beam and target jets. Field-Feynman fragmentation has been used in the hadronization. The sensitivity of the results to a few features of the Monte Carlo program has been studied. The results are found to be very sensitive to the method used to ensure overall energy conservation after the fragmentation of the four jets is complete. Results are also sensitive to the minimum momentum transfer in the QCD subprocesses and to the distribution of p/sub T/ to the jet axis and the multiplicities in the fragmentation. With reasonable choices of these features of the Monte Carlo program, good agreement with data at Fermilab/CERN SPS energies is obtained, comparable to the agreement achieved with more sophisticated parton-shower models. With other choices, however, the calculation gives qualitatively different results which are in strong disagreement with the data. These results have important implications for extracting physics conclusions from Monte Carlo calculations. It is not possible to test the validity of a particular model or distinguish between different models unless the Monte Carlo results are unambiguous and different models exhibit clearly different behavior

Análisis económico y ambiental de la agricultura alternativa. Caso: Alto Ricaurte - Boyacá

OpenAIRE

González A., Carolina; Gutiérrez B., Patricia

1992-01-01

El trabajo: Análisis económico y ambiental de la agricultura alternativa, caso: Alto Ricaurte Boyacá, aborda un tema de trascendencia como es la propuesta de un modelo de agricultura alternativa como base para un desarrollo sustentable, lo cual es, precisamente una de las principales preocupaciones de la comunidad internacional que motivó la convocatoria de la llamada Cumbre de la tierra que por estos días se reúne en Río de Janeiro.

Dependence of hydropower energy generation on forests in the Amazon Basin at local and regional scales.

Science.gov (United States)

Stickler, Claudia M; Coe, Michael T; Costa, Marcos H; Nepstad, Daniel C; McGrath, David G; Dias, Livia C P; Rodrigues, Hermann O; Soares-Filho, Britaldo S

2013-06-04

Tropical rainforest regions have large hydropower generation potential that figures prominently in many nations' energy growth strategies. Feasibility studies of hydropower plants typically ignore the effect of future deforestation or assume that deforestation will have a positive effect on river discharge and energy generation resulting from declines in evapotranspiration (ET) associated with forest conversion. Forest loss can also reduce river discharge, however, by inhibiting rainfall. We used land use, hydrological, and climate models to examine the local "direct" effects (through changes in ET within the watershed) and the potential regional "indirect" effects (through changes in rainfall) of deforestation on river discharge and energy generation potential for the Belo Monte energy complex, one of the world's largest hydropower plants that is currently under construction on the Xingu River in the eastern Amazon. In the absence of indirect effects of deforestation, simulated deforestation of 20% and 40% within the Xingu River basin increased discharge by 4-8% and 10-12%, with similar increases in energy generation. When indirect effects were considered, deforestation of the Amazon region inhibited rainfall within the Xingu Basin, counterbalancing declines in ET and decreasing discharge by 6-36%. Under business-as-usual projections of forest loss for 2050 (40%), simulated power generation declined to only 25% of maximum plant output and 60% of the industry's own projections. Like other energy sources, hydropower plants present large social and environmental costs. Their reliability as energy sources, however, must take into account their dependence on forests.

Monte Carlo methods and applications in nuclear physics

International Nuclear Information System (INIS)

Carlson, J.

1990-01-01

Monte Carlo methods for studying few- and many-body quantum systems are introduced, with special emphasis given to their applications in nuclear physics. Variational and Green's function Monte Carlo methods are presented in some detail. The status of calculations of light nuclei is reviewed, including discussions of the three-nucleon-interaction, charge and magnetic form factors, the coulomb sum rule, and studies of low-energy radiative transitions. 58 refs., 12 figs

Forest Management

Science.gov (United States)

S. Hummel; K. L. O' Hara

2008-01-01

Global variation in forests and in human cultures means that a single method for managing forests is not possible. However, forest management everywhere shares some common principles because it is rooted in physical and biological sciences like chemistry and genetics. Ecological forest management is an approach that combines an understanding of universal processes with...

Forest hydrology

Science.gov (United States)

Ge Sun; Devendra Amatya; Steve McNulty

2016-01-01

Forest hydrology studies the distribution, storage, movement, and quality of water and the hydrological processes in forest-dominated ecosystems. Forest hydrological science is regarded as the foundation of modern integrated water¬shed management. This chapter provides an overview of the history of forest hydrology and basic principles of this unique branch of...

Monte Carlo and analytic simulations in nanoparticle-enhanced radiation therapy

Directory of Open Access Journals (Sweden)

Paro AD

2016-09-01

Full Text Available Autumn D Paro,1 Mainul Hossain,2 Thomas J Webster,1,3,4 Ming Su1,4 1Department of Chemical Engineering, Northeastern University, Boston, MA, USA; 2NanoScience Technology Center and School of Electrical Engineering and Computer Science, University of Central Florida, Orlando, Florida, USA; 3Excellence for Advanced Materials Research, King Abdulaziz University, Jeddah, Saudi Arabia; 4Wenzhou Institute of Biomaterials and Engineering, Chinese Academy of Science, Wenzhou Medical University, Zhejiang, People’s Republic of China Abstract: Analytical and Monte Carlo simulations have been used to predict dose enhancement factors in nanoparticle-enhanced X-ray radiation therapy. Both simulations predict an increase in dose enhancement in the presence of nanoparticles, but the two methods predict different levels of enhancement over the studied energy, nanoparticle materials, and concentration regime for several reasons. The Monte Carlo simulation calculates energy deposited by electrons and photons, while the analytical one only calculates energy deposited by source photons and photoelectrons; the Monte Carlo simulation accounts for electron–hole recombination, while the analytical one does not; and the Monte Carlo simulation randomly samples photon or electron path and accounts for particle interactions, while the analytical simulation assumes a linear trajectory. This study demonstrates that the Monte Carlo simulation will be a better choice to evaluate dose enhancement with nanoparticles in radiation therapy. Keywords: nanoparticle, dose enhancement, Monte Carlo simulation, analytical simulation, radiation therapy, tumor cell, X-ray 

Texas' forests, 2008

Science.gov (United States)

James W. Bentley; Consuelo Brandeis; Jason A. Cooper; Christopher M. Oswalt; Sonja N. Oswalt; KaDonna Randolph

2014-01-01

This bulletin describes forest resources of the State of Texas at the time of the 2008 forest inventory. This bulletin addresses forest area, volume, growth, removals, mortality, forest health, timber product output, and the economy of the forest sector.

Changing markets - Medicinal plants in the markets of La Paz and El Alto, Bolivia.

Science.gov (United States)

Bussmann, Rainer W; Paniagua Zambrana, Narel Y; Moya Huanca, Laura Araseli; Hart, Robbie

2016-12-04

Given the importance of local markets as a source of medicinal plants for both healers and the population, literature on market flows and the value of the plant material traded is rather scarce. This stands in contrast to wealth of available information for other components of Bolivian ethnobotany. The present study attempts to remedy this situation by providing a detailed inventory of medicinal plant markets in the La Paz-El Alto metropolitan area, hypothesizing that both species composition, and medicinal applications, have changed considerably over time. From October 2013-October 2015 semi-structured interviews were conducted with 39 plant vendors between October 2013 and October 2015 in the Mercado Rodriguez, Mercado Calle Santa Cruz, Mercado Cohoni, Mercado Cota Cota, and Mercado Seguencoma and Mercado El Alto in order to elucidate more details on plant usage and provenance. The results of the present study were then compared to previous inventories of medicinal plants in La Paz and El Alto studies to elucidate changes over time and impact of interview techniques. In this study we encountered 163 plant species belonging to 127 genera and 58 families. In addition, 17 species could not be identified. This species richness is considerably higher than that reported in previous studies (2005, 129 species of 55 families; 2015, 94 identified species). While the overall distribution of illness categories is in line with older reports the number of species used per application, as well as the applications per species, were much higher in the present study. Overall, informant consensus was relatively low, which might be explained by the large number of new species that have entered the local pharmacopoeia in the last decade, although some species might simply have been missed by previous studies. In course of the present study it became apparent that even well known species might often be replaced by other apparently similar but botanically unrelated species due to

Oncogenes E6-E7 de los Papilomavirus Humanos de alto riesgo detectados por PCR en Biopsias de pene incluidas en parafina

Directory of Open Access Journals (Sweden)

I Guerrero

1999-01-01

Full Text Available La alta prevalencia del papilomavirus humano (PVH, referida a nivel mundial, en lesiones genitales de ambos sexos, el rol del varón como reservorio pasivo del virus, y el incremento de la mortalidad por cáncer genital en la mujer en nuestro país, motiva la detección y correlación de los oncogenes de los PVH de alto riesgo con la neoplasia de pene. Informamos de diez casos de biopsias de carcinoma escamoso de pene, incluidos en parafina, los cuales fueron investigados para la presencia de los oncogenes E6-E7 de PVH de alto riesgo, utilizando cebadores tipo específico para los PVH -16 y 18, mediante la reacción en cadena de la polimerasa (PCR. El 40% de los casos mostró un producto de amplificación ADN E6-E7 de los PVH estudiados, correspondiendo el 75% de ellos a detección simple por PVH-18 y el 25% presentó detección mixta ADN E6 - E7 del PVH-16 y 18 simultáneamente. El producto de amplificación fue sometido a comprobación por análisis de restricción específico. La prevalencia obtenida de los oncogenes E6-E7 de los PVH de alto riesgo, usando un método tan sensible como la PCR, apoya el rol de estos virus en el proceso de carcinogénesis de la neoplasia de pene.

High-efficiency wavefunction updates for large scale Quantum Monte Carlo

Science.gov (United States)

Kent, Paul; McDaniel, Tyler; Li, Ying Wai; D'Azevedo, Ed

Within ab intio Quantum Monte Carlo (QMC) simulations, the leading numerical cost for large systems is the computation of the values of the Slater determinants in the trial wavefunctions. The evaluation of each Monte Carlo move requires finding the determinant of a dense matrix, which is traditionally iteratively evaluated using a rank-1 Sherman-Morrison updating scheme to avoid repeated explicit calculation of the inverse. For calculations with thousands of electrons, this operation dominates the execution profile. We propose a novel rank- k delayed update scheme. This strategy enables probability evaluation for multiple successive Monte Carlo moves, with application of accepted moves to the matrices delayed until after a predetermined number of moves, k. Accepted events grouped in this manner are then applied to the matrices en bloc with enhanced arithmetic intensity and computational efficiency. This procedure does not change the underlying Monte Carlo sampling or the sampling efficiency. For large systems and algorithms such as diffusion Monte Carlo where the acceptance ratio is high, order of magnitude speedups can be obtained on both multi-core CPU and on GPUs, making this algorithm highly advantageous for current petascale and future exascale computations.

La fiesta xita: patrimonio biocultural mazahua de San Pedro el Alto, México

Directory of Open Access Journals (Sweden)

Alicia Y. Vásquez González

2016-01-01

Full Text Available Ante la crisis ambiental y de la modernidad, se propone conservar la riqueza cultural y biológica de la humanidad como patrimonio. Así, se valora y revalora el patrimonio vivo, como es la relación entre culturas y biodiversidad, que se manifiestan en paisajes y territorios de los pueblos originarios. Se propone la fiesta xita (viejos como patrimonio biocultural de la milpa mazahua. El estudio etnográfico se realizó en San Pedro el Alto, México. La fiesta xita es parte del patrimonio biocultural mazahua, mediante la cual se entretejen agrobiodiversidad y cultura, lo que permite que los mazahuas mantengan y desarrollen sus actividades sociales y culturales que han sido heredadas de generación en generación.

Estudio preliminar de la fauna Amphibia del valle de Sibundoy, Alto Putumayo, Colombia

Directory of Open Access Journals (Sweden)

Jonh Jairo Mueses-Cisneros

2004-07-01

(Valle de Sibundoy y por localidades muestreadas, se hace una comparación con las faunas realizadas en el transecto Napo en Ecuador y en el transecto Montañita-Alto Gabinete en Caquetá y se presenta una clave taxonómica y una serie de fotografías para facilitar su  reconocimiento. El estudio además aporta cinco especies nuevas para la ciencia, amplía el rango altitudinal de otras ocho y presenta tres más para ser incluidas al listado general de las especies de anfibios de Colombia, las cuales habían sido ya anteriormente reportadas para el país por otros autores, pero no habían sido tenidas en cuenta en el último listado.

The Monte Carlo Simulation Method for System Reliability and Risk Analysis

CERN Document Server

Zio, Enrico

2013-01-01

Monte Carlo simulation is one of the best tools for performing realistic analysis of complex systems as it allows most of the limiting assumptions on system behavior to be relaxed. The Monte Carlo Simulation Method for System Reliability and Risk Analysis comprehensively illustrates the Monte Carlo simulation method and its application to reliability and system engineering. Readers are given a sound understanding of the fundamentals of Monte Carlo sampling and simulation and its application for realistic system modeling.   Whilst many of the topics rely on a high-level understanding of calculus, probability and statistics, simple academic examples will be provided in support to the explanation of the theoretical foundations to facilitate comprehension of the subject matter. Case studies will be introduced to provide the practical value of the most advanced techniques.   This detailed approach makes The Monte Carlo Simulation Method for System Reliability and Risk Analysis a key reference for senior undergra...

Niveles subjetivos de estrés-recuperación en deportistas Costarricenses de alto rendimiento

OpenAIRE

Sánchez Ureña, M.Sc. Braulio; Ureña Bonilla, Pedro; Calleja González, Julio

2014-01-01

El objetivo del presente estudio fue determinar los niveles estrés-recuperación en deportistas costarricenses de alto rendimiento. Metodología: participaron un total de 239 sujetos, con una media de edad de 25,4 ± 4,8 años y con un promedio de 6,4 ± 4,6 años de entrenamiento en la primera división de su disciplina deportiva. Se utilizó la versión española del Cuestionario de Estrés-Recuperación para deportistas Resultados: los valores promedio para las variables relacionadas a la recuperación...

Invernaderos sostenibles para la producción de hortalizas en zonas alto andinas del Cusco

OpenAIRE

Zanabria Pacheco, Pedro Fortunato

2015-01-01

La investigación tiene como objetivo determinar la importancia en el uso de los invernaderos sostenibles para la producción de hortalizas en las zonas alto andinas del Curso. Para ello, se han diseñado y construido dos prototipos de invernaderos a 3330 msnm, variando la orientación y forma del techo, permitiendo la producción de hortalizas en la estación fría. A la vez está investigación servirá como base para aplicar el proyecto sobre la tecnología de producción sostenible de hortalizas bajo...

Estudio florístico y aportaciones a la conservación del alto Cabriel (Cuenca)

OpenAIRE

Mayoral García-Berlanga, Olga

2011-01-01

El alto Cabriel es un área montañosa inmersa en plena Serranía de Cuenca (Sistema Ibérico meridional), con altitudes entre los 910 y los 1.840 m. La zona de estudio com-prende 93.311 hectáreas -repartidas en 20 cuadrículas UTM de 10 km de lado- y 16 muni-cipios. Casi la mitad de la longitud del río Cabriel, afluente izquierdo del Júcar, discurre por la zona de estudio. Dominan las litologías calcáreas, destacando dos afloramientos de areniscas del Bundsandstein. Desde un punto de vista bio...

Monte Carlo simulation applied to alpha spectrometry

International Nuclear Information System (INIS)

Baccouche, S.; Gharbi, F.; Trabelsi, A.

2007-01-01

Alpha particle spectrometry is a widely-used analytical method, in particular when we deal with pure alpha emitting radionuclides. Monte Carlo simulation is an adequate tool to investigate the influence of various phenomena on this analytical method. We performed an investigation of those phenomena using the simulation code GEANT of CERN. The results concerning the geometrical detection efficiency in different measurement geometries agree with analytical calculations. This work confirms that Monte Carlo simulation of solid angle of detection is a very useful tool to determine with very good accuracy the detection efficiency.

Monte Carlo simulation of neutron scattering instruments

International Nuclear Information System (INIS)

Seeger, P.A.

1995-01-01

A library of Monte Carlo subroutines has been developed for the purpose of design of neutron scattering instruments. Using small-angle scattering as an example, the philosophy and structure of the library are described and the programs are used to compare instruments at continuous wave (CW) and long-pulse spallation source (LPSS) neutron facilities. The Monte Carlo results give a count-rate gain of a factor between 2 and 4 using time-of-flight analysis. This is comparable to scaling arguments based on the ratio of wavelength bandwidth to resolution width

Simulation of transport equations with Monte Carlo

International Nuclear Information System (INIS)

Matthes, W.

1975-09-01

The main purpose of the report is to explain the relation between the transport equation and the Monte Carlo game used for its solution. The introduction of artificial particles carrying a weight provides one with high flexibility in constructing many different games for the solution of the same equation. This flexibility opens a way to construct a Monte Carlo game for the solution of the adjoint transport equation. Emphasis is laid mostly on giving a clear understanding of what to do and not on the details of how to do a specific game

Microcanonical Monte Carlo

International Nuclear Information System (INIS)

Creutz, M.

1986-01-01

The author discusses a recently developed algorithm for simulating statistical systems. The procedure interpolates between molecular dynamics methods and canonical Monte Carlo. The primary advantages are extremely fast simulations of discrete systems such as the Ising model and a relative insensitivity to random number quality. A variation of the algorithm gives rise to a deterministic dynamics for Ising spins. This model may be useful for high speed simulation of non-equilibrium phenomena

A contribution Monte Carlo method

International Nuclear Information System (INIS)

Aboughantous, C.H.

1994-01-01

A Contribution Monte Carlo method is developed and successfully applied to a sample deep-penetration shielding problem. The random walk is simulated in most of its parts as in conventional Monte Carlo methods. The probability density functions (pdf's) are expressed in terms of spherical harmonics and are continuous functions in direction cosine and azimuthal angle variables as well as in position coordinates; the energy is discretized in the multigroup approximation. The transport pdf is an unusual exponential kernel strongly dependent on the incident and emergent directions and energies and on the position of the collision site. The method produces the same results obtained with the deterministic method with a very small standard deviation, with as little as 1,000 Contribution particles in both analog and nonabsorption biasing modes and with only a few minutes CPU time

Illinois' Forests 2005

Science.gov (United States)

Susan J. Crocker; Gary J. Brand; Brett J. Butler; David E. Haugen; Dick C. Little; Dacia M. Meneguzzo; Charles H. Perry; Ronald J. Piva; Barry T. Wilson; Christopher W. Woodall

2009-01-01

The first full, annualized inventory of Illinois' forests reports more than 4.5 million acres of forest land with an average of 459 trees per acre. Forest land is dominated by oak/hickory forest types, which occupy 65 percent of total forest land area. Seventy-two percent of forest land consists of sawtimber, 20 percent contains poletimber, and 8 percent contains...

Minnesota's Forests 2008

Science.gov (United States)

Patrick D. Miles; David Heinzen; Manfred E. Mielke; Christopher W. Woodall; Brett J. Butler; Ron J. Piva; Dacia M. Meneguzzo; Charles H. Perry; Dale D. Gormanson; Charles J. Barnett

2011-01-01

The second full annual inventory of Minnesota's forests reports 17 million acres of forest land with an average volume of more than 1,000 cubic feet per acre. Forest land is dominated by the aspen forest type, which occupies nearly 30 percent of the total forest land area. Twenty-eight percent of forest land consists of sawtimber, 35 percent poletimber, 35 percent...

A tale of two "forests": random forest machine learning AIDS tropical forest carbon mapping.

Science.gov (United States)

Mascaro, Joseph; Asner, Gregory P; Knapp, David E; Kennedy-Bowdoin, Ty; Martin, Roberta E; Anderson, Christopher; Higgins, Mark; Chadwick, K Dana

2014-01-01

Accurate and spatially-explicit maps of tropical forest carbon stocks are needed to implement carbon offset mechanisms such as REDD+ (Reduced Deforestation and Degradation Plus). The Random Forest machine learning algorithm may aid carbon mapping applications using remotely-sensed data. However, Random Forest has never been compared to traditional and potentially more reliable techniques such as regionally stratified sampling and upscaling, and it has rarely been employed with spatial data. Here, we evaluated the performance of Random Forest in upscaling airborne LiDAR (Light Detection and Ranging)-based carbon estimates compared to the stratification approach over a 16-million hectare focal area of the Western Amazon. We considered two runs of Random Forest, both with and without spatial contextual modeling by including--in the latter case--x, and y position directly in the model. In each case, we set aside 8 million hectares (i.e., half of the focal area) for validation; this rigorous test of Random Forest went above and beyond the internal validation normally compiled by the algorithm (i.e., called "out-of-bag"), which proved insufficient for this spatial application. In this heterogeneous region of Northern Peru, the model with spatial context was the best preforming run of Random Forest, and explained 59% of LiDAR-based carbon estimates within the validation area, compared to 37% for stratification or 43% by Random Forest without spatial context. With the 60% improvement in explained variation, RMSE against validation LiDAR samples improved from 33 to 26 Mg C ha(-1) when using Random Forest with spatial context. Our results suggest that spatial context should be considered when using Random Forest, and that doing so may result in substantially improved carbon stock modeling for purposes of climate change mitigation.

A tale of two "forests": random forest machine learning AIDS tropical forest carbon mapping.

Directory of Open Access Journals (Sweden)

Joseph Mascaro

Full Text Available Accurate and spatially-explicit maps of tropical forest carbon stocks are needed to implement carbon offset mechanisms such as REDD+ (Reduced Deforestation and Degradation Plus. The Random Forest machine learning algorithm may aid carbon mapping applications using remotely-sensed data. However, Random Forest has never been compared to traditional and potentially more reliable techniques such as regionally stratified sampling and upscaling, and it has rarely been employed with spatial data. Here, we evaluated the performance of Random Forest in upscaling airborne LiDAR (Light Detection and Ranging-based carbon estimates compared to the stratification approach over a 16-million hectare focal area of the Western Amazon. We considered two runs of Random Forest, both with and without spatial contextual modeling by including--in the latter case--x, and y position directly in the model. In each case, we set aside 8 million hectares (i.e., half of the focal area for validation; this rigorous test of Random Forest went above and beyond the internal validation normally compiled by the algorithm (i.e., called "out-of-bag", which proved insufficient for this spatial application. In this heterogeneous region of Northern Peru, the model with spatial context was the best preforming run of Random Forest, and explained 59% of LiDAR-based carbon estimates within the validation area, compared to 37% for stratification or 43% by Random Forest without spatial context. With the 60% improvement in explained variation, RMSE against validation LiDAR samples improved from 33 to 26 Mg C ha(-1 when using Random Forest with spatial context. Our results suggest that spatial context should be considered when using Random Forest, and that doing so may result in substantially improved carbon stock modeling for purposes of climate change mitigation.

Forest and water relationships: hydrologic implications of forestation campaigns in China

Science.gov (United States)

Ge Sun; Guoyi Zhou; Zhiqiang Zhang; Xiaohua Wei; Steven G. McNulty; James Vose

2005-01-01

Reforestation and afforestation (referred to forestation thereafter) campaigns in the past two decades have resulted in great increases in both forest land area and forest ecosystem productivity in China. Although the ecological benefits of forests are well accepted, the hydrologic consequences of man-made forests by forestation are unclear. Debate and confusion on the...

El poblamiento humano antiguo en el valle alto del Lozoya (Madrid

Directory of Open Access Journals (Sweden)

Belén Márquez

2008-01-01

Full Text Available La localización en superficie en el término municipal de Lozoya (Madrid de un conjunto de piezas de industria lítica achelense permite retrotraer el primer poblamiento de esta zona al Pleistoceno medio, cuando hasta la fecha las evidencias más antiguas en el valle alto del río Lozoya se remontaban a inicios del Pleistoceno superior. El nuevo hallazgo sitúa al río Lozoya en relación cultural con el poblamiento humano antiguo de los ríos Jarama y Manzanares.The discovery of a surface sample of Acheulean lithic industry at the Lozoya district (Madrid takes back to the Middle Pleistocene the first human peopling of this area. Until now themost ancient evidences at the Upper Lozoya Valley went back to the initial Upper Pleistocene. The new findings place the Lozoya River in cultural relation with the ancient human population of Jarama and Manzanares Rivers.

Potential of the Russian forests and forest industries

Energy Technology Data Exchange (ETDEWEB)

Anttonen, T.; Petrov, A.P. [eds.

1997-12-31

The publication contains the proceedings of the seminar `Potential of the Russian Forests and Forest Industries` held in Moscow, May 14-16, 1997. The seminar was one step along the road to spread knowledge and become acquainted with forestry and forest industries in northern Europe and Russia. The seminar proceedings contain a lot of fresh information concerning forestry and forest industries in Russia. Both have undergone many changes and reforms during the last few years

Ecological consequences of forest elephant declines for Afrotropical forests.

Science.gov (United States)

Poulsen, John R; Rosin, Cooper; Meier, Amelia; Mills, Emily; Nuñez, Chase L; Koerner, Sally E; Blanchard, Emily; Callejas, Jennifer; Moore, Sarah; Sowers, Mark

2017-10-27

Poaching is rapidly extirpating African forest elephants (Loxodonta cyclotis) from most of their historical range, leaving vast areas of elephant-free tropical forest. Elephants are ecological engineers that create and maintain forest habitat; thus, their loss will have large consequences for the composition and structure of Afrotropical forests. Through a comprehensive literature review, we evaluated the roles of forest elephants in seed dispersal, nutrient recycling, and herbivory and physical damage to predict the cascading ecological effects of their population declines. Loss of seed dispersal by elephants will favor tree species dispersed abiotically and by smaller dispersal agents, and tree species composition will depend on the downstream effects of changes in elephant nutrient cycling and browsing. Loss of trampling and herbivory of seedlings and saplings will result in high tree density with release from browsing pressures. Diminished seed dispersal by elephants and high stem density are likely to reduce the recruitment of large trees and thus increase homogeneity of forest structure and decrease carbon stocks. The loss of ecological services by forest elephants likely means Central African forests will be more like Neotropical forests, from which megafauna were extirpated thousands of years ago. Without intervention, as much as 96% of Central African forests will have modified species composition and structure as elephants are compressed into remaining protected areas. Stopping elephant poaching is an urgent first step to mitigating these effects, but long-term conservation will require land-use planning that incorporates elephant habitat into forested landscapes that are being rapidly transformed by industrial agriculture and logging. © 2017 Society for Conservation Biology.

The impact of Monte Carlo simulation: a scientometric analysis of scholarly literature

CERN Document Server

Pia, Maria Grazia; Bell, Zane W; Dressendorfer, Paul V

2010-01-01

A scientometric analysis of Monte Carlo simulation and Monte Carlo codes has been performed over a set of representative scholarly journals related to radiation physics. The results of this study are reported and discussed. They document and quantitatively appraise the role of Monte Carlo methods and codes in scientific research and engineering applications.

Exact Monte Carlo for molecules

International Nuclear Information System (INIS)

Lester, W.A. Jr.; Reynolds, P.J.

1985-03-01

A brief summary of the fixed-node quantum Monte Carlo method is presented. Results obtained for binding energies, the classical barrier height for H + H 2 , and the singlet-triplet splitting in methylene are presented and discussed. 17 refs

Evaluación de la exposición a selenio en Los Altos de Jalisco, México Evaluation of the exposure to selenium in Los Altos de Jalisco, México

Directory of Open Access Journals (Sweden)

Roberto Hurtado-Jiménez

2007-08-01

Full Text Available OBJETIVO: Evaluar la exposición a selenio (Se vía agua potable en los habitantes de Los Altos de Jalisco. MATERIAL Y MÉTODOS: Se determinó la concentración de Se en 125 pozos y se estimaron los niveles de exposición a Se en bebés, niños y adultos. RESULTADOS: La dosis de exposición y la ingestión de Se vía agua potable variaron en los siguientes rangos: a bebés: 1.3-6.7 µg/kg/d y 12.6-67.2 µg/d; b niños: 0.8-4.5 µg/kg/d y 16.8-89.6 µg/d; c adultos: 0.6-3.0 µg/kg/d y 33.6-179.2 µg/d. CONCLUSIONES: En este caso, la exposición a Se representa un riesgo potencial para la salud de la población, ya que en la mayoría de los casos es mayor que la recomendada por organismos internacionales de salud. Sin embargo, no es tan alta como para esperar la ocurrencia de selenosis.OBJECTIVE: To evaluate the exposure to selenium in drinking water in Los Altos de Jalisco (Jalisco State Heights. MATERIALS AND METHODS: The concentration of selenium was determined in 125 water wells, and the exposure doses to selenium were estimated for babies, children and adults. RESULTS: The estimated values of the exposure doses to selenium and total intake of selenium were in the following ranges, respectively: (a babies: 1.3-6.7 µg/kg/d and 12.6-67.2 µg/d; (b children: 0.8-4.5 µg/kg/d and 16.8-89.6 µg/d, (c adults: 0.6-3.0 µg/kg/d and 33.6-179.2 µg/d. CONCLUSIONS: The estimated exposure levels to selenium were higher than those recommended as optimum by international health organizations, representing a potential health risk. Nevertheless, estimated values are not high enough to produce selenosis.

No-compromise reptation quantum Monte Carlo

International Nuclear Information System (INIS)

Yuen, W K; Farrar, Thomas J; Rothstein, Stuart M

2007-01-01

Since its publication, the reptation quantum Monte Carlo algorithm of Baroni and Moroni (1999 Phys. Rev. Lett. 82 4745) has been applied to several important problems in physics, but its mathematical foundations are not well understood. We show that their algorithm is not of typical Metropolis-Hastings type, and we specify conditions required for the generated Markov chain to be stationary and to converge to the intended distribution. The time-step bias may add up, and in many applications it is only the middle of a reptile that is the most important. Therefore, we propose an alternative, 'no-compromise reptation quantum Monte Carlo' to stabilize the middle of the reptile. (fast track communication)

Exploring cluster Monte Carlo updates with Boltzmann machines.

Science.gov (United States)

Wang, Lei

2017-11-01

Boltzmann machines are physics informed generative models with broad applications in machine learning. They model the probability distribution of an input data set with latent variables and generate new samples accordingly. Applying the Boltzmann machines back to physics, they are ideal recommender systems to accelerate the Monte Carlo simulation of physical systems due to their flexibility and effectiveness. More intriguingly, we show that the generative sampling of the Boltzmann machines can even give different cluster Monte Carlo algorithms. The latent representation of the Boltzmann machines can be designed to mediate complex interactions and identify clusters of the physical system. We demonstrate these findings with concrete examples of the classical Ising model with and without four-spin plaquette interactions. In the future, automatic searches in the algorithm space parametrized by Boltzmann machines may discover more innovative Monte Carlo updates.

Exploring cluster Monte Carlo updates with Boltzmann machines

Science.gov (United States)

Wang, Lei

2017-11-01

Boltzmann machines are physics informed generative models with broad applications in machine learning. They model the probability distribution of an input data set with latent variables and generate new samples accordingly. Applying the Boltzmann machines back to physics, they are ideal recommender systems to accelerate the Monte Carlo simulation of physical systems due to their flexibility and effectiveness. More intriguingly, we show that the generative sampling of the Boltzmann machines can even give different cluster Monte Carlo algorithms. The latent representation of the Boltzmann machines can be designed to mediate complex interactions and identify clusters of the physical system. We demonstrate these findings with concrete examples of the classical Ising model with and without four-spin plaquette interactions. In the future, automatic searches in the algorithm space parametrized by Boltzmann machines may discover more innovative Monte Carlo updates.

Monte Carlo simulation of continuous-space crystal growth

International Nuclear Information System (INIS)

Dodson, B.W.; Taylor, P.A.

1986-01-01

We describe a method, based on Monte Carlo techniques, of simulating the atomic growth of crystals without the discrete lattice space assumed by conventional Monte Carlo growth simulations. Since no lattice space is assumed, problems involving epitaxial growth, heteroepitaxy, phonon-driven mechanisms, surface reconstruction, and many other phenomena incompatible with the lattice-space approximation can be studied. Also, use of the Monte Carlo method circumvents to some extent the extreme limitations on simulated timescale inherent in crystal-growth techniques which might be proposed using molecular dynamics. The implementation of the new method is illustrated by studying the growth of strained-layer superlattice (SLS) interfaces in two-dimensional Lennard-Jones atomic systems. Despite the extreme simplicity of such systems, the qualitative features of SLS growth seen here are similar to those observed experimentally in real semiconductor systems

Effect of error propagation of nuclide number densities on Monte Carlo burn-up calculations

International Nuclear Information System (INIS)

Tohjoh, Masayuki; Endo, Tomohiro; Watanabe, Masato; Yamamoto, Akio

2006-01-01

As a result of improvements in computer technology, the continuous energy Monte Carlo burn-up calculation has received attention as a good candidate for an assembly calculation method. However, the results of Monte Carlo calculations contain the statistical errors. The results of Monte Carlo burn-up calculations, in particular, include propagated statistical errors through the variance of the nuclide number densities. Therefore, if statistical error alone is evaluated, the errors in Monte Carlo burn-up calculations may be underestimated. To make clear this effect of error propagation on Monte Carlo burn-up calculations, we here proposed an equation that can predict the variance of nuclide number densities after burn-up calculations, and we verified this equation using enormous numbers of the Monte Carlo burn-up calculations by changing only the initial random numbers. We also verified the effect of the number of burn-up calculation points on Monte Carlo burn-up calculations. From these verifications, we estimated the errors in Monte Carlo burn-up calculations including both statistical and propagated errors. Finally, we made clear the effects of error propagation on Monte Carlo burn-up calculations by comparing statistical errors alone versus both statistical and propagated errors. The results revealed that the effects of error propagation on the Monte Carlo burn-up calculations of 8 x 8 BWR fuel assembly are low up to 60 GWd/t

Indiana's Forests 2008

Science.gov (United States)

Christopher W. Woodall; Mark N. Webb; Barry T. Wilson; Jeff Settle; Ron J. Piva; Charles H. Perry; Dacia M. Meneguzzo; Susan J. Crocker; Brett J. Butler; Mark Hansen; Mark Hatfield; Gary Brand; Charles. Barnett

2011-01-01

The second full annual inventory of Indiana's forests reports more than 4.75 million acres of forest land with an average volume of more than 2,000 cubic feet per acre. Forest land is dominated by the white oak/red oak/hickory forest type, which occupies nearly a third of the total forest land area. Seventy-six percent of forest land consists of sawtimber, 16...

Mixed-Forest Species Establishment in a Monodominant Forest in Central Africa: Implications for Tropical Forest Invasibility

Science.gov (United States)

Peh, Kelvin S.-H.; Sonké, Bonaventure; Séné, Olivier; Djuikouo, Marie-Noël K.; Nguembou, Charlemagne K.; Taedoumg, Hermann; Begne, Serge K.; Lewis, Simon L.

2014-01-01

Background Traits of non-dominant mixed-forest tree species and their synergies for successful co-occurrence in monodominant Gilbertiodendron dewevrei forest have not yet been investigated. Here we compared the tree species diversity of the monodominant forest with its adjacent mixed forest and then determined which fitness proxies and life history traits of the mixed-forest tree species were most associated with successful co-existence in the monodominant forest. Methodology/Principal Findings We sampled all trees (diameter in breast height [dbh]≥10 cm) within 6×1 ha topographically homogenous areas of intact central African forest in SE Cameroon, three independent patches of G. dewevrei-dominated forest and three adjacent areas (450–800 m apart). Monodominant G. dewevrei forest had lower sample-controlled species richness, species density and population density than its adjacent mixed forest in terms of stems with dbh≥10 cm. Analysis of a suite of population-level characteristics, such as relative abundance and geographical distribution, and traits such as wood density, height, diameter at breast height, fruit/seed dispersal mechanism and light requirement–revealed after controlling for phylogeny, species that co-occur with G. dewevrei tend to have higher abundance in adjacent mixed forest, higher wood density and a lower light requirement. Conclusions/Significance Our results suggest that certain traits (wood density and light requirement) and population-level characteristics (relative abundance) may increase the invasibility of a tree species into a tropical closed-canopy system. Such knowledge may assist in the pre-emptive identification of invasive tree species. PMID:24844914

Mixed-forest species establishment in a monodominant forest in central Africa: implications for tropical forest invasibility.

Directory of Open Access Journals (Sweden)

Kelvin S-H Peh

Full Text Available BACKGROUND: Traits of non-dominant mixed-forest tree species and their synergies for successful co-occurrence in monodominant Gilbertiodendron dewevrei forest have not yet been investigated. Here we compared the tree species diversity of the monodominant forest with its adjacent mixed forest and then determined which fitness proxies and life history traits of the mixed-forest tree species were most associated with successful co-existence in the monodominant forest. METHODOLOGY/PRINCIPAL FINDINGS: We sampled all trees (diameter in breast height [dbh]≥10 cm within 6×1 ha topographically homogenous areas of intact central African forest in SE Cameroon, three independent patches of G. dewevrei-dominated forest and three adjacent areas (450-800 m apart. Monodominant G. dewevrei forest had lower sample-controlled species richness, species density and population density than its adjacent mixed forest in terms of stems with dbh≥10 cm. Analysis of a suite of population-level characteristics, such as relative abundance and geographical distribution, and traits such as wood density, height, diameter at breast height, fruit/seed dispersal mechanism and light requirement-revealed after controlling for phylogeny, species that co-occur with G. dewevrei tend to have higher abundance in adjacent mixed forest, higher wood density and a lower light requirement. CONCLUSIONS/SIGNIFICANCE: Our results suggest that certain traits (wood density and light requirement and population-level characteristics (relative abundance may increase the invasibility of a tree species into a tropical closed-canopy system. Such knowledge may assist in the pre-emptive identification of invasive tree species.

Monte Carlo simulation of neutron counters for safeguards applications

International Nuclear Information System (INIS)

Looman, Marc; Peerani, Paolo; Tagziria, Hamid

2009-01-01

MCNP-PTA is a new Monte Carlo code for the simulation of neutron counters for nuclear safeguards applications developed at the Joint Research Centre (JRC) in Ispra (Italy). After some preliminary considerations outlining the general aspects involved in the computational modelling of neutron counters, this paper describes the specific details and approximations which make up the basis of the model implemented in the code. One of the major improvements allowed by the use of Monte Carlo simulation is a considerable reduction in both the experimental work and in the reference materials required for the calibration of the instruments. This new approach to the calibration of counters using Monte Carlo simulation techniques is also discussed.

Using a Numerical Model to Assess the Geomorphic Impacts of Forest Management Scenarios on Streams

Science.gov (United States)

Davidson, S. L.; Eaton, B. C.

2014-12-01

In-stream large wood governs the morphology of many small to intermediate streams, while riparian vegetation influences bank strength and channel pattern. Forest management practices such as harvesting and fire suppression therefore dramatically influence channel processes and associated aquatic habitat. The primary objective of this research is to compare the impacts of three common forest scenarios - natural fire disturbance, forest harvesting with a riparian buffer, and fire suppression - on the volume of in-channel wood and the complexity of aquatic habitat in channels at a range of scales. Each scenario is explored through Monte Carlo simulations run over a period of 1000 years using a numerical reach scale channel simulator (RSCS), with variations in tree toppling rate and forest density used to represent each forest management trajectory. The habitat complexity associated with each scenario is assessed based on the area of the bed occupied by pools and spawning sized sediment, the availability of wood cover, and the probability of avulsion. Within the fire scenario, we also use the model to separately investigate the effects of root decay and recovery on equilibrium channel geometry by varying the rooting depth and associated bank strength through time. The results show that wood loading and habitat complexity are influenced by the timing and magnitude of wood recruitment, as well as channel scale. The forest harvesting scenario produces the lowest wood loads and habitat complexity so long as the buffer width is less than the average mature tree height. The natural fire cycle produces the greatest wood loading and habitat complexity, but also the greatest variability because these streams experience significant periods without wood recruitment as forests regenerate. In reaches that experience recurrent fires, width increases in the post-fire period as roots decay, at times producing a change in channel pattern when a threshold width to depth ratio is

Monte Carlo methods and applications in nuclear physics

Energy Technology Data Exchange (ETDEWEB)

Carlson, J.

1990-01-01

Monte Carlo methods for studying few- and many-body quantum systems are introduced, with special emphasis given to their applications in nuclear physics. Variational and Green's function Monte Carlo methods are presented in some detail. The status of calculations of light nuclei is reviewed, including discussions of the three-nucleon-interaction, charge and magnetic form factors, the coulomb sum rule, and studies of low-energy radiative transitions. 58 refs., 12 figs.

Forest owners' perceptions of ecotourism: Integrating community values and forest conservation.

Science.gov (United States)

Rodríguez-Piñeros, Sandra; Mayett-Moreno, Yesica

2015-03-01

The use of forest land for ecotourism has been well accepted due to its ability to provide income to local people and to conserve the forest. Preparing the forest with infrastructure to attract and educate visitors has been reported of importance. This study applied Q methodology in a small rural community of the State of Puebla, Mexico, to reveal forest owners' perceptions to build infrastructure in their forest as part of their ecotourism project. It also discloses forest owners' underlying motives to use their forest for ecotourism. Ecotourism is perceived as a complementary activity to farming that would allow women to be involved in community development. Low impact infrastructure is desired due to forest owners' perception to preserve the forest for the overall community well-being.

Forest inventory: role in accountability for sustainable forest management

Science.gov (United States)

Lloyd C. Irland

2007-01-01

Forest inventory can play several roles in accountability for sustainable forest management. A first dimension is accountability for national performance. The new field of Criteria and Indicators is an expression of this need. A more familiar role for the U.S. Department of Agriculture Forest Service Forest Inventory and Analysis (FIA) program is for assessment and...

Research on Monte Carlo improved quasi-static method for reactor space-time dynamics

International Nuclear Information System (INIS)

Xu Qi; Wang Kan; Li Shirui; Yu Ganglin

2013-01-01

With large time steps, improved quasi-static (IQS) method can improve the calculation speed for reactor dynamic simulations. The Monte Carlo IQS method was proposed in this paper, combining the advantages of both the IQS method and MC method. Thus, the Monte Carlo IQS method is beneficial for solving space-time dynamics problems of new concept reactors. Based on the theory of IQS, Monte Carlo algorithms for calculating adjoint neutron flux, reactor kinetic parameters and shape function were designed and realized. A simple Monte Carlo IQS code and a corresponding diffusion IQS code were developed, which were used for verification of the Monte Carlo IQS method. (authors)

Kontrola tačnosti rezultata u simulacijama Monte Karlo / Accuracy control in Monte Carlo simulations

Directory of Open Access Journals (Sweden)

Nebojša V. Nikolić

2010-04-01

Full Text Available U radu je demonstrirana primena metode automatizovanog ponavljanja nezavisnih simulacionih eksperimenata sa prikupljanjem statistike slučajnih procesa, u dostizanju i kontroli tačnosti simulacionih rezultata u simulaciji sistema masovnog opsluživanja Monte Karlo. Metoda se zasniva na primeni osnovnih stavova i teorema matematičke statistike i teorije verovatnoće. Tačnost simulacionih rezultata dovedena je u direktnu vezu sa brojem ponavljanja simulacionih eksperimenata. / The paper presents an application of the Automated Independent Replication with Gathering Statistics of the Stochastic Processes Method in achieving and controlling the accuracy of simulation results in the Monte Carlo queuing simulations. The method is based on the application of the basic theorems of the theory of probability and mathematical statistics. The accuracy of the simulation results is linked with a number of independent replications of simulation experiments.

Final Report: 06-LW-013, Nuclear Physics the Monte Carlo Way

International Nuclear Information System (INIS)

Ormand, W.E.

2009-01-01

This is document reports the progress and accomplishments achieved in 2006-2007 with LDRD funding under the proposal 06-LW-013, 'Nuclear Physics the Monte Carlo Way'. The project was a theoretical study to explore a novel approach to dealing with a persistent problem in Monte Carlo approaches to quantum many-body systems. The goal was to implement a solution to the notorious 'sign-problem', which if successful, would permit, for the first time, exact solutions to quantum many-body systems that cannot be addressed with other methods. In this document, we outline the progress and accomplishments achieved during FY2006-2007 with LDRD funding in the proposal 06-LW-013, 'Nuclear Physics the Monte Carlo Way'. This project was funded under the Lab Wide LDRD competition at Lawrence Livermore National Laboratory. The primary objective of this project was to test the feasibility of implementing a novel approach to solving the generic quantum many-body problem, which is one of the most important problems being addressed in theoretical physics today. Instead of traditional methods based matrix diagonalization, this proposal focused a Monte Carlo method. The principal difficulty with Monte Carlo methods, is the so-called 'sign problem'. The sign problem, which will discussed in some detail later, is endemic to Monte Carlo approaches to the quantum many-body problem, and is the principal reason that they have not been completely successful in the past. Here, we outline our research in the 'shifted-contour method' applied the Auxiliary Field Monte Carlo (AFMC) method

Lattice gauge theories and Monte Carlo simulations

International Nuclear Information System (INIS)

Rebbi, C.

1981-11-01

After some preliminary considerations, the discussion of quantum gauge theories on a Euclidean lattice takes up the definition of Euclidean quantum theory and treatment of the continuum limit; analogy is made with statistical mechanics. Perturbative methods can produce useful results for strong or weak coupling. In the attempts to investigate the properties of the systems for intermediate coupling, numerical methods known as Monte Carlo simulations have proved valuable. The bulk of this paper illustrates the basic ideas underlying the Monte Carlo numerical techniques and the major results achieved with them according to the following program: Monte Carlo simulations (general theory, practical considerations), phase structure of Abelian and non-Abelian models, the observables (coefficient of the linear term in the potential between two static sources at large separation, mass of the lowest excited state with the quantum numbers of the vacuum (the so-called glueball), the potential between two static sources at very small distance, the critical temperature at which sources become deconfined), gauge fields coupled to basonic matter (Higgs) fields, and systems with fermions

Time step length versus efficiency of Monte Carlo burnup calculations

International Nuclear Information System (INIS)

Dufek, Jan; Valtavirta, Ville

2014-01-01

Highlights: • Time step length largely affects efficiency of MC burnup calculations. • Efficiency of MC burnup calculations improves with decreasing time step length. • Results were obtained from SIE-based Monte Carlo burnup calculations. - Abstract: We demonstrate that efficiency of Monte Carlo burnup calculations can be largely affected by the selected time step length. This study employs the stochastic implicit Euler based coupling scheme for Monte Carlo burnup calculations that performs a number of inner iteration steps within each time step. In a series of calculations, we vary the time step length and the number of inner iteration steps; the results suggest that Monte Carlo burnup calculations get more efficient as the time step length is reduced. More time steps must be simulated as they get shorter; however, this is more than compensated by the decrease in computing cost per time step needed for achieving a certain accuracy

Evaluation of some energy options for the replacement of wood in cooking food in the Region Alto Patia

International Nuclear Information System (INIS)

Santacruz, Carlos A; Macia, Andres F

2006-01-01

In the present article an environmental evaluation appears financier and social of the different power solutions for the substitution from the wood in the food baking. First begins to describe the region of the Alto Patia, soon the objective population defines, this is village the San Juanito, in addition and they define the evaluation criteria the propose technologies that are: Solar, Electric, coal, gas and finally it is selected most suitable for this region.

Uncertainty in Forest Net Present Value Estimations

Directory of Open Access Journals (Sweden)

Ilona Pietilä

2010-09-01

Full Text Available Uncertainty related to inventory data, growth models and timber price fluctuation was investigated in the assessment of forest property net present value (NPV. The degree of uncertainty associated with inventory data was obtained from previous area-based airborne laser scanning (ALS inventory studies. The study was performed, applying the Monte Carlo simulation, using stand-level growth and yield projection models and three alternative rates of interest (3, 4 and 5%. Timber price fluctuation was portrayed with geometric mean-reverting (GMR price models. The analysis was conducted for four alternative forest properties having varying compartment structures: (A a property having an even development class distribution, (B sapling stands, (C young thinning stands, and (D mature stands. Simulations resulted in predicted yield value (predicted NPV distributions at both stand and property levels. Our results showed that ALS inventory errors were the most prominent source of uncertainty, leading to a 5.1–7.5% relative deviation of property-level NPV when an interest rate of 3% was applied. Interestingly, ALS inventory led to significant biases at the property level, ranging from 8.9% to 14.1% (3% interest rate. ALS inventory-based bias was the most significant in mature stand properties. Errors related to the growth predictions led to a relative standard deviation in NPV, varying from 1.5% to 4.1%. Growth model-related uncertainty was most significant in sapling stand properties. Timber price fluctuation caused the relative standard deviations ranged from 3.4% to 6.4% (3% interest rate. The combined relative variation caused by inventory errors, growth model errors and timber price fluctuation varied, depending on the property type and applied rates of interest, from 6.4% to 12.6%. By applying the methodology described here, one may take into account the effects of various uncertainty factors in the prediction of forest yield value and to supply the

Tropical Forest Gain and Interactions amongst Agents of Forest Change

Directory of Open Access Journals (Sweden)

Sean Sloan

2016-02-01

Full Text Available The tropical deforestation literature advocates multi-agent enquiry in recognition that key dynamics arise from inter-agent interactions. Studies of tropical forest-cover gain have lagged in this respect. This article explores the roles and key aspects of interactions shaping natural forest regeneration and active reforestation in Eastern Panama since 1990. It employs household surveys of agricultural landholders, interviews with community forest-restoration organisations, archival analysis of plantation reforestation interests, satellite image analysis of forest-cover change, and the consideration of State reforestation policies. Forest-cover gain reflected a convergence of interests and land-use trends amongst agents. Low social and economic costs of sustained interaction and organisation enabled extensive forest-cover gain, but low transaction costs did not. Corporate plantation reforestation rose to the fore of regional forest-cover gain via opportunistic land sales by ranchers and economic subsidies indicative of a State preference for autonomous, self-organising forest-cover gain. This reforestation follows a recent history of neoliberal frontier development in which State-backed loggers and ranchers similarly displaced agriculturalists. Community institutions, long neglected by the State, struggled to coordinate landholders and so effected far less forest-cover gain. National and international commitments to tropical forest restoration risk being similarly characterised as ineffective by a predominance of industrial plantation reforestation without greater State support for community forest management.