The empirical analysis of cigarette tax avoidance and illicit trade in Vietnam, 1998-2010.

Illicit trade carries the potential to magnify existing tobacco-related health care costs through increased availability of untaxed and inexpensive cigarettes. What is known with respect to the magnitude of illicit trade for Vietnam is produced primarily by the industry, and methodologies are typically opaque. Independent assessment of the illicit cigarette trade in Vietnam is vital to tobacco control policy. This paper measures the magnitude of illicit cigarette trade for Vietnam between 1998 and 2010 using two methods, discrepancies between legitimate domestic cigarette sales and domestic tobacco consumption estimated from surveys, and trade discrepancies as recorded by Vietnam and trade partners. The results indicate that Vietnam likely experienced net smuggling in during the period studied. With the inclusion of adjustments for survey respondent under-reporting, inward illicit trade likely occurred in three of the four years for which surveys were available. Discrepancies in trade records indicate that the value of smuggled cigarettes into Vietnam ranges from $100 million to $300 million between 2000 and 2010 and that these cigarettes primarily originate in Singapore, Hong Kong, Macao, Malaysia, and Australia. Notable differences in trends over time exist between the two methods, but by comparison, the industry estimates consistently place the magnitude of illicit trade at the upper bounds of what this study shows. The unavailability of annual, survey-based estimates of consumption may obscure the true, annual trend over time. Second, as surveys changed over time, estimates relying on them may be inconsistent with one another. Finally, these two methods measure different components of illicit trade, specifically consumption of illicit cigarettes regardless of origin and smuggling of cigarettes into a particular market. However, absent a gold standard, comparisons of different approaches to illicit trade measurement serve efforts to refine and improve measurement approaches and estimates.

Spatial and temporal scales of neuronal correlation in visual area V4.

The spiking activity of nearby cortical neurons is correlated on both short and long time scales. Understanding this shared variability in firing patterns is critical for appreciating the representation of sensory stimuli in ensembles of neurons, the coincident influences of neurons on common targets, and the functional implications of microcircuitry. Our knowledge about neuronal correlations, however, derives largely from experiments that used different recording methods, analysis techniques, and cortical regions. Here we studied the structure of neuronal correlation in area V4 of alert macaques using recording and analysis procedures designed to match those used previously in primary visual cortex (V1), the major input to V4. We found that the spatial and temporal properties of correlations in V4 were remarkably similar to those of V1, with two notable differences: correlated variability in V4 was approximately one-third the magnitude of that in V1 and synchrony in V4 was less temporally precise than in V1. In both areas, spontaneous activity (measured during fixation while viewing a blank screen) was approximately twice as correlated as visual-evoked activity. The results provide a foundation for understanding how the structure of neuronal correlation differs among brain regions and stages in cortical processing and suggest that it is likely governed by features of neuronal circuits that are shared across the visual cortex.

Monte Carlo Analysis and Physics Characterization of a Novel Nanoparticle Detector for Medical and Micro-dosimetry Applications

The outcomes for both (i) radiation therapy and (ii) preclinical small animal radio- biology studies are dependent on the delivery of a known quantity of radiation to a specific and intentional location. Adverse effects can result from these procedures if the dose to the target is too high or low, and can also result from an incorrect spatial distribution in which nearby normal healthy tissue can be undesirably damaged by poor radiation delivery techniques. Thus, in mice and humans alike, the spatial dose distributions from radiation sources should be well characterized in terms of the absolute dose quantity, and with pin-point accuracy. When dealing with the steep spatial dose gradients consequential to either (i) high dose rate (HDR) brachytherapy or (ii) within the small organs and tissue inhomogeneities of mice, obtaining accurate and highly precise dose results can be very challenging, considering commercially available radiation detection tools, such as ion chambers, are often too large for in-vivo use.

In this dissertation two tools are developed and applied for both clinical and preclinical radiation measurement. The first tool is a novel radiation detector for acquiring physical measurements, fabricated from an inorganic nano-crystalline scintillator that has been fixed on an optical fiber terminus. This dosimeter allows for the measurement of point doses to sub-millimeter resolution, and has the ability to be placed in-vivo in humans and small animals. Real-time data is displayed to the user to provide instant quality assurance and dose-rate information. The second tool utilizes an open source Monte Carlo particle transport code, and was applied for small animal dosimetry studies to calculate organ doses and recommend new techniques of dose prescription in mice, as well as to characterize dose to the murine bone marrow compartment with micron-scale resolution.

Hardware design changes were implemented to reduce the overall fiber diameter to <0.9 mm for the nano-crystalline scintillator based fiber optic detector (NanoFOD) system. Lower limits of device sensitivity were found to be approximately 0.05 cGy/s. Herein, this detector was demonstrated to perform quality assurance of clinical 192Ir HDR brachytherapy procedures, providing comparable dose measurements as thermo-luminescent dosimeters and accuracy within 20% of the treatment planning software (TPS) for 27 treatments conducted, with an inter-quartile range ratio to the TPS dose value of (1.02-0.94=0.08). After removing contaminant signals (Cerenkov and diode background), calibration of the detector enabled accurate dose measurements for vaginal applicator brachytherapy procedures. For 192Ir use, energy response changed by a factor of 2.25 over the SDD values of 3 to 9 cm; however a cap made of 0.2 mm thickness silver reduced energy dependence to a factor of 1.25 over the same SDD range, but had the consequence of reducing overall sensitivity by 33%.

For preclinical measurements, dose accuracy of the NanoFOD was within 1.3% of MOSFET measured dose values in a cylindrical mouse phantom at 225 kV for x-ray irradiation at angles of 0, 90, 180, and 270˝. The NanoFOD exhibited small changes in angular sensitivity, with a coefficient of variation (COV) of 3.6% at 120 kV and 1% at 225 kV. When the NanoFOD was placed alongside a MOSFET in the liver of a sacrificed mouse and treatment was delivered at 225 kV with 0.3 mm Cu filter, the dose difference was only 1.09% with use of the 4x4 cm collimator, and -0.03% with no collimation. Additionally, the NanoFOD utilized a scintillator of 11 µm thickness to measure small x-ray fields for microbeam radiation therapy (MRT) applications, and achieved 2.7% dose accuracy of the microbeam peak in comparison to radiochromic film. Modest differences between the full-width at half maximum measured lateral dimension of the MRT system were observed between the NanoFOD (420 µm) and radiochromic film (320 µm), but these differences have been explained mostly as an artifact due to the geometry used and volumetric effects in the scintillator material. Characterization of the energy dependence for the yttrium-oxide based scintillator material was performed in the range of 40-320 kV (2 mm Al filtration), and the maximum device sensitivity was achieved at 100 kV. Tissue maximum ratio data measurements were carried out on a small animal x-ray irradiator system at 320 kV and demonstrated an average difference of 0.9% as compared to a MOSFET dosimeter in the range of 2.5 to 33 cm depth in tissue equivalent plastic blocks. Irradiation of the NanoFOD fiber and scintillator material on a 137Cs gamma irradiator to 1600 Gy did not produce any measurable change in light output, suggesting that the NanoFOD system may be re-used without the need for replacement or recalibration over its lifetime.

For small animal irradiator systems, researchers can deliver a given dose to a target organ by controlling exposure time. Currently, researchers calculate this exposure time by dividing the total dose that they wish to deliver by a single provided dose rate value. This method is independent of the target organ. Studies conducted here used Monte Carlo particle transport codes to justify a new method of dose prescription in mice, that considers organ specific doses. Monte Carlo simulations were performed in the Geant4 Application for Tomographic Emission (GATE) toolkit using a MOBY mouse whole-body phantom. The non-homogeneous phantom was comprised of 256x256x800 voxels of size 0.145x0.145x0.145 mm3. Differences of up to 20-30% in dose to soft-tissue target organs was demonstrated, and methods for alleviating these errors were suggested during whole body radiation of mice by utilizing organ specific and x-ray tube filter specific dose rates for all irradiations.

Monte Carlo analysis was used on 1 µm resolution CT images of a mouse femur and a mouse vertebra to calculate the dose gradients within the bone marrow (BM) compartment of mice based on different radiation beam qualities relevant to x-ray and isotope type irradiators. Results and findings indicated that soft x-ray beams (160 kV at 0.62 mm Cu HVL and 320 kV at 1 mm Cu HVL) lead to substantially higher dose to BM within close proximity to mineral bone (within about 60 µm) as compared to hard x-ray beams (320 kV at 4 mm Cu HVL) and isotope based gamma irradiators (137Cs). The average dose increases to the BM in the vertebra for these four aforementioned radiation beam qualities were found to be 31%, 17%, 8%, and 1%, respectively. Both in-vitro and in-vivo experimental studies confirmed these simulation results, demonstrating that the 320 kV, 1 mm Cu HVL beam caused statistically significant increased killing to the BM cells at 6 Gy dose levels in comparison to both the 320 kV, 4 mm Cu HVL and the 662 keV, 137Cs beams.

Semi-Parametric Analysis of Extreme Forest Fires

Forest fires can cause extensive damage to natural resources and properties. They can also destroy wildlife habitat, affect the forest ecosystem and threaten human lives. In this paper incidences of extreme wildland fires are modelled by a point process model which incorporates time-trend. A model based on a generalised Pareto distribution is used to model data on acres of wildland burnt by extreme fire in the US since 1825. A semi-parametric smoothing approach, which is very useful in exploratory analysis of changes in extremes, is illustrated with the maximum likelihood method to estimate model parameters.

Spatial variability of air pollution in the vicinity of a permanent monitoring station in central Paris

The strong spatial and temporal variability of traffic-related air pollution detected at roadside locations in a number of European cities has raised the question of how representative the site and time period of air quality measurements actually can be. To address this question, a 7-month sampling campaign was carried out on a major road axis (Avenue Leclerc) leading to a very busy intersection (Place Basch) in central Paris, covering the surroundings of a permanent air quality monitoring station. This station has recorded the highest CO and NOx concentrations during recent years in the region of Paris. Diffusive BTX samplers as well as a mobile monitoring unit equipped with real-time CO, NOx and O3 analysers and meteorological instruments were used to reveal the small-scale pollution gradients and their temporal trends near the permanent monitoring station. The diffusive measurements provided 7-day averages of benzene, toluene, xylene and other hydrocarbons at different heights above the ground and distances from the kerb covering summer and winter periods. Relevant traffic and meteorological data were also obtained on an hourly basis. Furthermore, three semiempirical dispersion models (STREET-SRI, OSPM and AEOLIUS) were tested for an asymmetric canyon location in Av. Leclerc. The analysis of this comprehensive data set has helped to assess the representativeness of air quality monitoring information.

Fracture mechanics analysis of cracks in solder joint intermetallic compounds

The trend towards miniaturization of electronic products leads to the need for very small sized solder joints. Therefore, there is a higher reliability risk that too large a fraction of solder joints will transform into Intermetallic Compounds (IMCs) at the solder interface. In this paper, fracture mechanics study of the IMC layer for SnPb and Pb-free solder joints was carried out using finite element numerical computer modelling method. It is assumed that only one crack is present in the IMC layer. Linear Elastic Fracture Mechanics (LEFM) approach is used for parametric study of the Stress Intensity Factors (SIF, KI and KII), at the predefined crack in the IMC layer of solder butt joint tensile sample. Contrary to intuition, it is revealed that a thicker IMC layer in fact increases the reliability of solder joint for a cracked IMC. Value of KI and KII are found to decrease with the location of the crack further away from the solder interfaces while other parameters are constant. Solder thickness and strain rate were also found to have a significant influence on the SIF values. It has been found that soft solder matrix generates non-uniform plastic deformation across the solder-IMC interface near the crack tip that is responsible to obtain higher KI and KII.

Local examination of skin diffusion using FTIR spectroscopic imaging and multivariate target factor analysis

In the context of trans-dermal drug delivery it is very important to have mechanistic insight into the barrier function of the skin's stratum corneum and the diffusion mechanisms of topically applied drugs. Currently spectroscopic imaging techniques are evolving which enable a spatial examination of various types of samples in a dynamic way. ATR-FTIR imaging opens up the possibility to monitor spatial diffusion profiles across the stratum corneum of a skin sample. Multivariate data analyses methods based on factor analysis are able to provide insight into the large amount of spectroscopically complex and highly overlapping signals generated. Multivariate target factor analysis was used for spectral resolution and local diffusion profiles with time through stratum corneum. A model drug, 4-cyanophenol in polyethylene glycol 600 and water was studied. Results indicate that the average diffusion profiles between spatially different locations show similar profiles despite the heterogeneous nature of the biological sample and the challenging experimental set-up.

Spatial, seasonal and long-term fluctuations of plankton in relation to hydroclimatic features in the English Channel, Celtic Sea and Bay of Biscay

Plankton collected by the Continuous Plankton Recorder (CPR) survey were investigated for the English Channel, Celtic Sea and Bay of Biscay from 1979 to 1995. The main goal was to study the relationship between climate and plankton and to understand the factors influencing it. In order to take into account the spatial and temporal structure of biological data, a three-mode principal component analysis (PCA) was developed. It not only identified 5 zones characterised by their similar biological composition and by the seasonal and inter-annual evolution of the plankton, it also made species associations based on their location and year-to-year change. The studied species have stronger year-to-year fluctuations in abundance over the English Channel and Celtic Sea than the species offshore in the Bay of Biscay. The changes in abundance of plankton in the English Channel are negatively related to inter-annual changes of climatic conditions from December to March (North Atlantic Oscillation [NAO] index and air temperature). Thus, the negative relationship shown by Fromentin and Planque (1996; Mar Ecol Prog Ser 134:111-118) between year-to-year changes of Calanus finmarchicus abundance in the northern North Atlantic and North Sea and NAO was also found for the most abundant copepods in the Channel. However, the hypothesis proposed to explain the plankton/NAO relationship is different for this region and a new hypothesis is proposed. In the Celtic Sea, a relationship between the planktonic assemblage and the air temperature was detected, but it is weaker than for the English Channel. No relationship was found for the Bay of Biscay. Thus, the local physical environment and the biological composition of these zones appear to modify the relationship between winter climatic conditions and the year-to-year fluctuations of the studied planktonic species. This shows, therefore, that the relationship between climate and plankton is difficult to generalise.

Environmental selection and resource allocation determine spatial patterns in picophytoplankton cell size

Here we describe a new trait-based model for cellular resource allocation that we use to investigate the relative importance of different drivers for small cell size in phytoplankton. Using the model, we show that increased investment in nonscalable structural components with decreasing cell size leads to a trade-off between cell size, nutrient and light affinity, and growth rate. Within the most extreme nutrient-limited, stratified environments, resource competition theory then predicts a trend toward larger minimum cell size with increasing depth. We demonstrate that this explains observed trends using a marine ecosystem model that represents selection and adaptation of a diverse community defined by traits for cell size and subcellular resource allocation. This framework for linking cellular physiology to environmental selection can be used to investigate the adaptive response of the marine microbial community to environmental conditions and the adaptive value of variations in cellular physiology.

DISCOVERY 2010: Spatial and temporal variability in a dynamic polar ecosystem

The Scotia Sea has been a focus of biological- and physical oceanographic study since the Discovery expeditions in the early 1900s. It is a physically energetic region with some of the highest levels of productivity in the Southern Ocean. It is also a region within which there have been greater than average levels of change in upper water column temperature. We describe the results of three cruises transecting the central Scotia Sea from south to north in consecutive years and covering spring, summer and autumn periods. We also report on some community level syntheses using both current-day and historical data from this region. A wide range of parameters were measured during the field campaigns, covering the physical oceanography of the region, air–sea CO2 fluxes, macro- and micronutrient concentrations, the composition and biomass of the nano-, micro- and mesoplankton communities, and the distribution and biomass of Antarctic krill and mesopelagic fish. Process studies examined the effect of iron-stress on the physiology of primary producers, reproduction and egestion in Antarctic krill and the transfer of stable isotopes between trophic layers, from primary consumers up to birds and seals. Community level syntheses included an examination of the biomass-spectra, food-web modelling, spatial analysis of multiple trophic layers and historical species distributions. The spatial analyses in particular identified two distinct community types: a northern warmer water community and a southern cold community, their boundary being broadly consistent with the position of the Southern Antarctic Circumpolar Current Front (SACCF). Temperature and ice cover appeared to be the dominant, over-riding factors in driving this pattern. Extensive phytoplankton blooms were a major feature of the surveys, and were persistent in areas such as South Georgia. In situ and bioassay measurements emphasised the important role of iron inputs as facilitators of these blooms. Based on seasonal DIC deficits, the South Georgia bloom was found to contain the strongest seasonal carbon uptake in the ice-free zone of the Southern Ocean. The surveys also encountered low-production, iron-limited regions, a situation more typical of the wider Southern Ocean. The response of primary and secondary consumers to spatial and temporal heterogeneity in production was complex. Many of the life-cycles of small pelagic organisms showed a close coupling to the seasonal cycle of food availability. For instance, Antarctic krill showed a dependence on early, non-ice-associated blooms to facilitate early reproduction. Strategies to buffer against environmental variability were also examined, such as the prevalence of multiyear life-cycles and variability in energy storage levels. Such traits were seen to influence the way in which Scotia Sea communities were structured, with biomass levels in the larger size classes being higher than in other ocean regions. Seasonal development also altered trophic function, with the trophic level of higher predators increasing through the course of the year as additional predator-prey interactions emerged in the lower trophic levels. Finally, our studies re-emphasised the role that the simple phytoplankton-krill-higher predator food chain plays in this Southern Ocean region, particularly south of the SACCF. To the north, alternative food chains, such as those involving copepods, macrozooplankton and mesopelagic fish, were increasingly important. Continued ocean warming in this region is likely to increase the prevalence of such alternative such food chains with Antarctic krill predicted to move southwards.

Comparing copepod time-series in the north of Spain: Spatial autocorrelation of community composition

Four time-series of copepod species biomass in the north of Spain were contrasted to demonstrate spatial autocorrelation of local communities and their responses to short-term local and regional variability in oceanographic conditions. The series represented coastal and oceanic environments along a marked gradient of influence of seasonal upwelling from Galicia to the Mar Cantábrico (S Bay of Biscay), and each one included at least 10 years of continuous data collected at monthly frequency. Community composition (i.e. species number and diversity) was very consistent through the region, but local variations in the presence of new species and the relative proportions of common species allowed for the characterisation of the response to the environment at each site. Small-sized species were more frequent near the coast. A few species, however, captured the main patterns of variability in all series. Calanus helgolandicus and Acartia (mainly Acartia clausi) were generally the main contributors to total biomass, while other species as Paracalanus parvus and Clausocalanus spp. were important only at some locations. Most copepod indices were positively correlated with upwelling, either considering the whole community (biomass, species richness and diversity) or individual species, but only in the coastal series analysed since 1991. Copepods in the nearby ocean, however, showed negative correlations with upwelling in the period 1960–1986. The effects of upwelling may have been modulated by local factors, as showed by the increases in biomass, number of species and diversity in associations with increases in sea surface temperature in Galicia, while in the Mar Cantábrico only the warming-tolerant species increased and those typical of upwelling decreased. Density stratification of the water column was associated with decreases in total copepod biomass in Galicia, while it favoured the increase in species richness in the Mar Cantábrico. Nearly all significant responses of copepods to environmental variability were delayed by up to 5 months, showing the importance of considering time-lags in the analysis of temporal responses of zooplankton.

Long-Term Retrospective Analysis of Mackerel Spawning in the North Sea: A New Time Series and Modeling Approach to CPR Data

We present a unique view of mackerel (Scomber scombrus) in the North Sea based on a new time series of larvae caught by the Continuous Plankton Recorder (CPR) survey from 1948-2005, covering the period both before and after the collapse of the North Sea stock. Hydrographic backtrack modelling suggested that the effect of advection is very limited between spawning and larvae capture in the CPR survey. Using a statistical technique not previously applied to CPR data, we then generated a larval index that accounts for both catchability as well as spatial and temporal autocorrelation. The resulting time series documents the significant decrease of spawning from before 1970 to recent depleted levels. Spatial distributions of the larvae, and thus the spawning area, showed a shift from early to recent decades, suggesting that the central North Sea is no longer as important as the areas further west and south. These results provide a consistent and unique perspective on the dynamics of mackerel in this region and can potentially resolve many of the unresolved questions about this stock

Spatial Heterogeneity and Genetic Variation in the Copepod Neocalanus Cristatus Along Two Transects in the North Pacific Sampled By the Continuous Plankton Recorder

We present a macrogeographic study of spatial heterogeneity in an important subarctic Pacific copepod and describe the first genetic analysis of population structure using Continuous Plankton Recorder (CPR) samples. Samples of Neocalanus cristatus were collected at a constant depth of similar to 7 m from two CPR tow-routes, (i) an east-west similar to 6500-km transect from Vancouver Island, Canada to Hokkaido Island, Japan, and (ii) a north-south transect of similar to 2250 km from Anchorage, Alaska to Tacoma, Washington. Analysis of these samples revealed three features of the biology of N. cristatus. First, N. cristatus undergoes small-scale diel vertical migration that is larger among stages CV- adult (3-6 times more abundant at 7 m at night), than stages CI-CIV (only 2-4 times higher at night). Secondly, while there were no regions where N. cristatus did not appear, each transect sampled a few large-scale macrogeographic patches. Thirdly, an analysis of molecular variation, using a partial sequence of the N. cristatus cytochrome oxidase I gene, revealed that 7.3% (P < 0.0001) of the total genetic variation among N. cristatus sampled from macrogeographic patches by the CPR could be explained by spatial heterogeneity. We suggest that spatial heterogeneity at macrogeographic scales may be important in plankton evolution.

Genetic analysis of Continuous Plankton Recorder (CPR) samples

Genetic analysis of Continuous Plankton Recorder (CPR) samples is enabling greater taxonomic resolution and the study of plankton population structure. Here, we present some results from the genetic analysis of CPR samples collected in the North Sea and north-eastern Atlantic that reveal the impacts of climate on benthic-pelagic coupling and the food web. We show that pronounced changes in the North Sea meroplankton are related to an increased abundance and spatial distribution of the larvae of the benthic echinoderm, Echinocardium cordatum. Key stages of reproduction in E. cordatum, gametogenesis and spawning, are influenced by winter and spring sea temperature (January-May). A stepwise increase in sea temperature after 1987, which has created warmer conditions earlier in the year, together with increased summer phytoplankton, may benefit the reproduction and survival of this benthic species. Competition between the larvae of E. cordatum and other holozooplanlcton taxa may now be altering the trophodynamics of the summer pelagic ecosystem. In the north-eastern Atlantic the genetic analysis of fish larvae sampled by the CPR has revealed an unprecedented increase in the abundance of juvenile snake pipefish, Entelurus aequoreiis since 2002. We argue that increased sea surface temperatures in winter and spring when the eggs of E. aqueoreus, which are brooded by the male, are developing and the young larvae are growing in the plankton are a likely cause. The increased abundance of this species in Atlantic and adjacent European seas already appears to be influencing the marine food web.

Microbial spatial variability: An example from the Celtic Sea

In July 2004, dominant populations of microbial ultraplankton (<5 μm), in the surface of the Celtic Sea (between UK and Eire), were repeatedly mapped using flow cytometry, at 1.5 km resolution over a region of diameter 100 km. The numerically dominant representatives of all basic functional types were enumerated including one group of phototrophic bacteria (Syn), two groups of phytoplankton (PP, NP), three groups of heterotrophic bacterioplankton (HB) and the regionally dominant group of heterotrophic protists (HP). The distributions of all organisms showed strong spatial variability with little relation to variability in physical fields such as salinity and temperature. Furthermore, there was little agreement between distributions of different organisms. The only linear correlation consistently explaining more than 50% of the variance between any pairing of the organism groups enumerated is between two different groups of HB. Specifically, no linear, or non-linear, relationship is found between any pairings of SYB, PP or HB groups with their protist predators HP. Looking for multiple dependencies, factor analysis reveals three groupings: Syn, PP and low nucleic acid content HB (LNA); high nucleic acid content HB (HNA); HP and NP. Even the manner in which the spatial variability of Syn, PP and HB abundance varies as a function of lengthscale (represented by a semivariogram) differs significantly from that for HP. In summary, although all microbial planktonic groups enumerated are present and numerically dominant throughout the region studied, at face value the relationships between them seem weak. Nevertheless, the behaviour of a simple, illustrative ecological model, with strongly interacting phototrophs and heterotrophs, with stochastic forcing, is shown to be consistent with the observed poor correlations and differences in how spatial variability varies with lengthscale. Thus, our study suggests that a comparison of microbial abundances alone may not discern strong underlying trophic interactions. Specific knowledge of these processes, in particular grazing, will be required to explain the causes of the observed microbial spatial variability and its resulting consequences for the functioning of the ecosystem.