Spectral deconvolution and operational use of stripping ratios in airborne radiometrics

Spectral deconvolution using stripping ratios for a set of pre-defined energy windows is the simplest means of reducing the most important part of gamma-ray spectral information. In this way, the effective interferences between the measured peaks are removed, leading, through a calibration, to clear estimates of radionuclide inventory. While laboratory measurements of stripping ratios are relatively easy to acquire, with detectors placed above small-scale calibration pads of known radionuclide concentrations, the extrapolation to measurements at altitudes where airborne survey detectors are used bring difficulties such as air-path attenuation and greater uncertainties in knowing ground level inventories. Stripping ratios are altitude dependent, and laboratory measurements using various absorbers to simulate the air-path have been used with some success. Full-scale measurements from an aircraft require a suitable location where radionuclide concentrations vary little over the field of view of the detector (which may be hundreds of metres). Monte Carlo simulations offer the potential of full-scale reproduction of gamma-ray transport and detection mechanisms. Investigations have been made to evaluate stripping ratios using experimental and Monte Carlo methods.     

Experimental and Monte Carlo simulation studies of open cylindrical radon monitoring device using CR-39 detector

There are several methods of measuring radon concentrations but nuclear track detector cylindrical dosimeters are widely employed. In this investigation, the consequence of effective volumes of the dosimeters on the registration of alpha tracks in a CR-39 detector was studied. In a series of experiments an optimum radius for a CR-39-based open cylindrical radon dosimeter was found to be about 3 cm. Monte Carlo simulation techniques have been employed to verify the experimental results. In this context, a computer code Monte Carlo simulation dosimetry (MOCSID) was developed. Monte Carlo simulation experiments gave the optimum radius of the dosimeters as 3.0 cm. The experimental results are in good agreement with those obtained by Monte Carlo design calculations. In addition to this, plate-out effects of radon progeny were also studied. It was observed that the contribution of radon progeny (218Po and 214Po) plated-out on the wall of the dosimeters increases with an increase of dosimeter radii and then decrease to 0 at a radius of about 3 cm if a point detector has been installed at the center of the dosimeter base. In the code MOCSID different types of random number generators were employed. The results of this research are very useful for designing an optimum size of radon dosimeters.     

Monte Carlo based calibration of scintillation detectors for laboratory and in situ gamma ray measurements

The calibration of scintillation detectors for gamma radiation in a well characterized setup can be transferred to other geometries using Monte Carlo simulations to account for the differences between the calibration and the other geometry. In this study a calibration facility was used that is constructed from bricks of well-known activity concentrations of ⁴⁰K and of radionuclides from the ²³⁸U- and ²³²Th-series. Transfer of the calibration was attempted to a Marinelli beaker geometry with the detector inside a lead shield and to an in situ application with the detector positioned on a sand bed. In general this resulted in good correspondence (within 5-10%) between the activity concentrations derived using the transferred calibration and activities that were derived by independent measurements. Some discrepancies were identified that were attributed to coincident summing in the natural decay series and interference of radon.     

Low level radioactivity measurements with phoswich detectors using coincident techniques and digital pulse processing analysis

A new system has been developed for the detection of low radioactivity levels of fission products and actinides using coincidence techniques. The device combines a phoswich detector for alpha/beta/gamma-ray recognition with a fast digital card for electronic pulse analysis. The phoswich can be used in a coincident mode by identifying the composed signal produced by the simultaneous detection of alpha/beta particles and X-rays/gamma particles. The technique of coincidences with phoswich detectors was proposed recently to verify the Nuclear Test Ban Treaty (NTBT) which established the necessity of monitoring low levels of gaseous fission products produced by underground nuclear explosions. With the device proposed here it is possible to identify the coincidence events and determine the energy and type of coincident particles. The sensitivity of the system has been improved by employing liquid scintillators and a high resolution low energy germanium detector. In this case it is possible to identify simultaneously by alpha/gamma coincidence transuranic nuclides present in environmental samples without necessity of performing radiochemical separation. The minimum detectable activity was estimated to be 0.01 Bq kg(-1) for 0.1 kg of soil and 1000 min counting.     

Full-spectrum analysis of natural gamma-ray spectra

In this paper, a new system to measure natural gamma-radiation in situ will be presented. This system combines a high-efficiency BGO scintillation detector with full-spectrum data analysis (FSA). This technique uses the (nearly) full spectral shape and the so-called 'standard spectra' to calculate the activity concentrations of 40K, 232Th and 238U present in a geological matrix (sediment, rock, etc.). We describe the FSA and the determination of the standard spectra. Standard spectra are constructed for various geometries and a comparison in intensity and shape will be made. The performance of such a system has been compared to a more traditional system, consisting of a NaI detector in combination with the 'windows' analysis. For count rates typically encountered in field experiments, the same accuracy is obtained 10-20 times faster using the new system. This allows for shorter integration times and hence shorter measurements or a better spatial resolution. The applicability of such a system will be illustrated via an example of an airborne experiment in which the new system produced results comparable to those of much larger traditional systems. This paper will conclude with a discussion of the current status of the system and an outlook for future research.     

Monte Carlo calculation of entire in situ gamma-ray spectra

We present a method for the synthesis of entire in situ gamma-ray spectra based on Monte Carlo calculations and measured data that characterize the detector properties. The method can serve for the determination of the effective depth of 137Cs in soil based on the information contained in the low-energy part of an in situ spectrum. Effective depth is defined as the depth of a plane distribution of 137Cs beneath the surface that reproduces the fluence energy and angular distribution at 1 m above the ground of gamma rays belonging to the real 137Cs distributions. We managed to reproduce the measured in situ spectra with our method and to demonstrate that the method allows the determination of the effective depth of 137Cs with a precision of 10(-2) m. The method requires minimal experimental characterization of the detector and is not sensitive to the details of the detector model and the soil composition and density employed in the Monte Carlo calculations.     

Measurement and modelling of radioxenon plumes in the Ottawa Valley

Since 2001 a real-time radiation monitoring network of Canadian nuclear facilities and major population centres has been implemented for response to nuclear incidents including a possible terrorist attack. Unshielded NaI(Tl) spectroscopic detectors are employed to measure gamma radiation from airborne radioactivity and radioactivity deposited on the ground. These detectors are composed of a standard 3'x3' cylindrical NaI(Tl) spectrometers with data storage and integrated telemetry. Some of the detectors have been deployed in the Ottawa Valley near Chalk River Laboratories and Ottawa, which has a complex radioxenon environment due to the proximity of nuclear power reactors, and medical isotope facilities. Although not a health threat, these releases have provided an opportunity for the Canadian Meteorological Centre and the Commissariat à l'Energie Atomique to validate their meteorological models. The meteorological models of the two organizations are in good agreement on the origin and the source terms of these releases.     

Parametric uncertainty analysis on hydrodynamic coefficients in groundwater numerical models using Monte Carlo method and RPEM

Groundwater resources are the only source of water in many arid and semi-arid regions. It is important to manage these resources to have a sustainable development. However, there are many factors influencing the accuracy of the results in groundwater modeling. In this research, the uncertainty of two important groundwater model parameters (hydraulic conductivity and specific yield) were considered as the main sources of uncertainty in estimating water level in an unconfined aquifer, in Iran. For this purpose, a simple method called Rosenblueth Point Estimate Method (RPEM) was used to assess groundwater modeling parametric uncertainty, and its performance was compared with Monte Carlo method as a very complicated and time-consuming method. According to calibrated values of hydraulic conductivity and specific yield, several uncertainty intervals were considered to analyze uncertainty. The results showed that the optimum interval for hydraulic conductivity was 40% increase–30% decrease of the calibrated values in both Monte Carlo and RPEM methods. This interval for specific yield was 200% increase–90% decrease of the calibrated values. RPEM showed better performance using the evaluating indices in comparison with Monte Carlo method for both hydraulic conductivity and specific yield with 43% and 17% higher index values, respectively. These results can be used in groundwater management and future prediction of groundwater level.

     

Review of dynamical models for external dose calculations based on Monte Carlo simulations in urbanised areas

After an accidental release of radionuclides to the inhabited environment the external gamma irradiation from deposited radioactivity contributes significantly to the radiation exposure of the population for extended periods. For evaluating this exposure pathway, three main model requirements are needed: (i) to calculate the air kerma value per photon emitted per unit source area, based on Monte Carlo (MC) simulations; (ii) to describe the distribution and dynamics of radionuclides on the diverse urban surfaces; and (iii) to combine all these elements in a relevant urban model to calculate the resulting doses according to the actual scenario. This paper provides an overview about the different approaches to calculate photon transport in urban areas and about several dose calculation codes published. Two types of Monte Carlo simulations are presented using the global and the local approaches of photon transport. Moreover, two different philosophies of the dose calculation, the "location factor method" and a combination of relative contamination of surfaces with air kerma values are described. The main features of six codes (ECOSYS, EDEM2M, EXPURT, PARATI, TEMAS, URGENT) are highlighted together with a short model-model features intercomparison.     

Evaluation of Monte Carlo-based calibrations of HPGe detectors for in situ gamma-ray spectrometry

The aim of this work was to evaluate the use of Monte Carlo-based calibrations for in situ gamma-ray spectrometry. We have performed in situ measurements at five different sites in Sweden using HPGe detectors to determine ground deposition activity levels of ¹³⁷Cs from the 1986 Chernobyl accident. Monte Carlo-calculated efficiency calibration factors were compared with corresponding values calculated using a more traditional semi-empirical method. In addition, results for the activity ground deposition were also compared with activity densities found in soil samples. In order to facilitate meaningful comparisons between the different types of results, the combined standard uncertainty of in situ measurements was assessed for both calibration methods. Good agreement, both between the two calibration methods, and between in situ measurements and soil samples, was found at all five sites. Uncertainties in in situ measurements for the given measurement conditions, about 20 years after the fallout occurred, were found to be in the range 15–20% (with a coverage factor k = 1, i.e. with a confidence interval of about 68%).     

Recent advances in aerial gamma-ray surveying

Aerial gamma-ray surveying uses NaI(Tl) detectors mounted in small aircraft to measure gamma radiation, emitted from the earth's surface. The data are collected as gamma-ray spectra, typically with 1 s counting times, from which are derived K, U and Th concentrations in the ground. Applications of aerial surveying include geological mapping for mineral exploration, soil mapping for agriculture, pollution studies and location of lost sources. Recent advances in applying statistical methods to the spectral data have resulted in large reductions in the noise levels in the surveys. Some of the methods available to do this include noise adjusted singular value decomposition (NASVD) [Proceedings of Exploration 97: Fourth Decennial International Conference on Mineral Exploration (1997) 753] and maximum noise fraction (MNF) and enhanced MNF (eMNF) [Explor. Geophys. 31 (2000) 73]. These methods, in general, apply normalization for variance to the spectra, use a principal component method to obtain the "significant" components of the data and reconstruct cleaned spectra, which are then processed in a standard manner to get radionuclide concentrations. However, they differ in the detail of the application and thus give slightly different results. In this paper, the application of noise reduction methods to various synthetic surveys is used to examine the strengths and weaknesses of the methods. In tests where there are high correlations between U and Th, the eMNF method performs best although the results are improved by prior clustering of the data by the Th/U ratio. If the data show no correlations, then the effectiveness of all the noise removal methods is reduced. If a data set is small (<1500 spectra), then MNF appears to be the better method. Consideration of the various tests suggests an optimum process whereby spectra are sorted into groups by the Th/U ratio of areas identified in a standard processing and then cleaned by eMNF or MNF, depending on the number of spectra in each group.     

Quantitative modeling of freshwater stress in the nine water basins of Tanzania

Freshwater scarcity is a global issue of environmental concern that threatens agricultural production and human health. In this study, we established freshwater stress indices (WSIs) for the nine water basins of Tanzania by using the quantity of freshwater available and various water uses. The relationship between water availability and different water uses, including environmental water requirements, was analyzed, with uncertainty and sensitivity analysis performed by a Monte Carlo simulation technique. Extreme WSI values close to 1.00 were obtained in the Rufiji, Pangani, and Wami-ruvu basins, Internal drainage, and Lake Rukwa, while low and moderate WSI values ranging from 0.03 to 0.84 were found in Lake Victoria and the Ruvuma, Tanganyika, and Nyasa basins. This study adds further knowledge on the level of freshwater scarcity, relationships between water availability and different water uses, and suggests policy options to reduce freshwater scarcity at the basin level for sustainable water supply.     

Impacts of climate change on food security and agriculture sector in Malaysia

This study attempts to analyze the short- and long-run impacts of the probable change in rainfall and temperature simultaneously on food availability and access to food issues, as the two dimensions of food security, in Malaysia. It uses an integrated method comprising of a stochastic method and a computable general equilibrium model using the latest (2010) input–output table published in 2015. The stochastic method, which relates to the Monte Carlo simulation, provides the probable changes in rainfall and temperature patterns and their probability of occurrence based on historical data of rainfall and temperature and crop productivity. It was found that, simultaneous variation of rainfall and temperature, in both the short- and long-run, contracts the economic performance of Malaysia. Findings also show the negative impact of rainfall–temperature variability, in both time periods, on food availability and access to food due to a reduction in the supply of agricultural products, a commodity inflation pressure and a reduction in household income. Moreover, results suggest that the climate variability shocks lead to a reduction in the consumption and welfare of all household groups, particularly in rural areas.     

A whole-plant mathematical model for the phytoextraction of lead (Pb) by maize

Phytoextraction is a technology that uses plants to remove heavy metals from contaminated soils. Although it is economically attractive compared to other methods, little attention has been paid to the mathematical modeling of the mechanisms involved. In this work, we simulate the phytoextraction of Pb using a mechanistic system dynamics modelling approach and the physiology model of maize (Zea mays) as a model system as it is a good Pb accumulator and translocator. Simulation results showed that precipitation is the most important mechanism related to the uptake of Pb from the ground. The most important model parameters have been identified through sensitivity analysis.     

Effect of a global warming-induced increase in typhoon intensity on urban productivity in Taiwan

A number of scientists have recently conducted research that shows that tropical cyclone intensity is likely to increase in the future due to the warming effect of greenhouse gases on surface sea temperatures. The aim of this paper is to establish what would be the likely decrease in the productivity of urban workers due to an increase in tropical cyclone-related downtime. The methodology used simulates future tropical cyclones by magnifying the intensity of historical tropical cyclones between the years 1978 and 2008. It then uses a Monte Carlo simulation to obtain the expected number of hours that a certain area can expect to be affected by winds of a given strength. It shows how annual downtime from tropical cyclones could increase from 1.5% nowadays to up to 2.2% by 2085, an increase of almost 50%. This decrease in productivity could result in a loss of up to 0.7% of the annual Taiwanese GDP by 2085.     

A Monte Carlo study of radon detection in cylindrical diffusion chambers

The functioning of radon diffusion chambers was studied using the Monte Carlo code RAMMX developed here. The alpha particles from radon are assumed randomly produced in the volume of the cylinder, and those from the progeny are assumed to originate randomly at the cylindrical surface. The energy spectrum, the distribution of incident angles, and the distribution of path lengths of the alpha particles on the detector were obtained. These quantities vary depending on input parameters such as initial alpha particle energy, radius and depth of the diffusion chamber, detector size and atmospheric pressure. The calculated energy spectrum for both ²²²Rn and ²²⁰Rn was compared with experiment, permitting the identification of each peak and its origin, and a better understanding of radon monitoring. Three aspects not considered in previous calculations are progeny alphas coming from surfaces of the monitor, taking into account the atmospheric pressure, and including the isotope ²²⁰Rn.     

Location, identification, and size distribution of depleted uranium grains in reservoir sediments

The location, nature, and size distribution of uranium-rich grains in sediment layers can be identified by sunbursts of etched particle tracks if each sample is pressed against a track detector, next irradiated with thermal neutrons, and the detectors then chemically etched to reveal fission tracks. The total track abundance from the sample is a measure of the 235U content; hence, if the bulk uranium (mostly 238U) has been measured, the two sets of results give the depletion or enrichment of the uranium. Sunbursts of tracks mark the locations of low-abundance, high-uranium grains allowing them to be singled out for further study.     

Imprecision of dose predictions for radionuclides released to the environment: An application of a Monte Carlo simulation technique

An evaluation of the imprecision in dose predictions has been performed using current dose assessment models and present knowledge of the variability or uncertainty in model parameter values. The propagation of parameter uncertainties is demonstrated using a Monte Carlo technique for elemental ¹³¹I transported via the pasture-cow-milk-child pathway. The results indicate that when site-specific information is not available the imprecision inherent in the predictions for this pathway is potentially large. Generally, the 99th percentile in thyroid dose for children was predicted to be approximately an order of magnitude greater than the median value. The potential consequences of the imprecision in dose for radiation protection purposes are discussed.