Validation of non‑stationary precipitation series for site‑specific impact assessment: comparison of two statistical downscaling techniques

Statistical downscaling (SD) methods have become a popular, low-cost and accessible means of bridging the gap between the coarse spatial resolution at which climate models output climate scenarios and the finer spatial scale at which impact modellers require these scenarios, with various different SD techniques used for a wide range of applications across the world. This paper compares the Generator for Point Climate Change (GPCC) model and the Statistical DownScaling Model (SDSM)—two contrasting SD methods—in terms of their ability to generate precipitation series under non-stationary conditions across ten contrasting global climates. The mean, maximum and a selection of distribution statistics as well as the cumulative frequencies of dry and wet spells for four different temporal resolutions were compared between the models and the observed series for a validation period. Results indicate that both methods can generate daily precipitation series that generally closely mirror observed series for a wide range of non-stationary climates. However, GPCC tends to overestimate higher precipitation amounts, whilst SDSM tends to underestimate these. This infers that GPCC is more likely to overestimate the effects of precipitation on a given impact sector, whilst SDSM is likely to underestimate the effects. GPCC performs better than SDSM in reproducing wet and dry day frequency, which is a key advantage for many impact sectors. Overall, the mixed performance of the two methods illustrates the importance of users performing a thorough validation in order to determine the influence of simulated precipitation on their chosen impact sector.

HOLOGRAPHIC IMAGES OF IODINE ATOMS ON A SILVER SURFACE FROM ELECTRON-EMISSION PATTERNS

BARTON 1 has suggested that photoelectron interference patterns may be used directly as holograms to obtain atomic-resolution images of surface structures. Bulk structures have been obtained previously by this means from experimental patterns of high-energy Kikuchi(quasi-elastically scattered) and Auger electrons 2,3. Here we test the feasibility of this technique for determination of surface structures using Auger intensity patterns obtained 4,5 from iodine chemisorbed on a pseudomorphic silver monolayer on Pt{111}. By direct numerical holographic inversion, we obtain three-dimensional images which show that iodine adatoms are located in hollows of 3-fold symmetry on the surface. The images yield the site symmetry with good atomic resolution in the surface plane, but suffer from poor resolution along the Ag-I axis. We anticipate that data with better angular resolution obtained at low temperatures would improve the spatial resolution of such images.

Statistical Analysis of Small Ellerman Bomb Events

The properties of Ellerman bombs (EBs), small-scale brightenings in the Hα line wings, have proved difficult to establish because their size is close to the spatial resolution of even the most advanced telescopes. Here, we aim to infer the size and lifetime of EBs using high-resolution data of an emerging active region collected using the Interferometric BIdimensional Spectrometer (IBIS) and Rapid Oscillations of the Solar Atmosphere (ROSA) instruments as well as the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO). We develop an algorithm to track EBs through their evolution, finding that EBs can often be much smaller (around 0.3″) and shorter-lived (less than one minute) than previous estimates. A correlation between G-band magnetic bright points and EBs is also found. Combining SDO/HMI and G-band data gives a good proxy of the polarity for the vertical magnetic field. It is found that EBs often occur both over regions of opposite polarity flux and strong unipolar fields, possibly hinting at magnetic reconnection as a driver of these events.The energetics of EB events is found to follow a power-law distribution in the range of a nanoflare (1022-25 ergs).

Astrochemistry

In the last 40 years a wide range of molecules, including neutrals, cations and anions, containing up to 13 atoms—in addition to detections of C₆₀ and C₇₀ — have been found in the harsh environment of the interstellar medium. The exquisite sensitivity and very high spectral and, more recently, spatial resolution, of modern telescopes has enabled the physics of star formation to be probed through rotational line emission. In this article, I review the basic properties of interstellar clouds and the processes that initiate the chemistry and generate chemical complexity, particularly in regions of star and planet formation. Our understanding of astrochemistry has evolved over the years. Before 1990, the general consensus was that molecules were formed in binary, gas-phase, or volume, reactions, most importantly ion-neutral reactions despite the very low ionization in clouds. Since then, observations have indicated unambiguously that there is also a contribution from surface processes, particularly on the icy mantles that form around refractory grain cores in cold, dense gas. The balance between these two processes depends on particular physical conditions and can vary during the life cycle of a particular volume of interstellar cloud.The complex chemistry that occurs in space is driven mostly through interaction of thegas with cosmic ray protons, a source of ionization that enables a rich ion-neutral chemistry. In addition, I show that the interaction between the gas and the dust in cold, dense regionsalso leads to additional chemical complexity through reactions that take place in ices at onlya few tens of degrees above absolute zero. Although densities are low compared to those in terrestrial environments, the extremely long life times of interstellar clouds and their enormous sizes, enable complex molecules to be synthesised and detected. I show that in some instances, particularly in reactions involving deuterium, the rotational populations of reactants, together with spin-selection rules, can determine the detailed abundances. Although the review is mainly focused on regions associated with star formation, I also consider chemistry in other interesting astronomical regions — in the early Universe and in the envelopes formed by mass loss during the final stages of stellar evolution.

Supernovae and radio transients in M82

We present optical and near-infrared (NIR) photometry and NIR spectroscopy of SN 2004am, the only optically detected supernova (SN) in M82. These demonstrate that SN 2004am was a highly reddened Type II-P SN similar to the low-luminosity Type II-P events such as SNe 1997D and 2005cs. We show that SN 2004am was located coincident with the obscured super star cluster M82-L, and from the cluster age infer a progenitor mass of 12{^{+ 7}_{- 3}} M⊙. In addition to this, we present a high spatial resolution Gemini-North Telescope K-band adaptive optics image of the site of SN 2008iz and a second transient of uncertain nature, both detected so far only at radio wavelengths. Using image subtraction techniques together with archival data from the Hubble Space Telescope, we are able to recover a NIR transient source coincident with both objects. We find the likely extinction towards SN 2008iz to be not more than AV ˜ 10. The nature of the second transient remains elusive and we regard an extremely bright microquasar in M82 as the most plausible scenario.

Binaural Reproduction of Finite Difference Simulations Using Spherical Array Processing

Due to its efficiency and simplicity, the finite-difference time-domain method is becoming a popular choice for solving wideband, transient problems in various fields of acoustics. So far, the issue of extracting a binaural response from finite difference simulations has only been discussed in the context of embedding a listener geometry in the grid. In this paper, we propose and study a method for binaural response rendering based on a spatial decomposition of the sound field. The finite difference grid is locally sampled using a volumetric array of receivers, from which a plane wave density function is computed and integrated with free-field head related transfer functions, in the spherical harmonics domain. The volumetric array is studied in terms of numerical robustness and spatial aliasing. Analytic formulas that predict the performance of the array are developed, facilitating spatial resolution analysis and numerical binaural response analysis for a number of finite difference schemes. Particular emphasis is placed on the effects of numerical dispersion on array processing and on the resulting binaural responses. Our method is compared to a binaural simulation based on the image method. Results indicate good spatial and temporal agreement between the two methods.

Evolution and Enabling Capabilities of Spatially Resolved Techniques for the Characterization of Heterogeneously Catalyzed Reactions

The development and optimization of catalysts and catalytic processes requires knowledge of reaction kinetics and mechanisms. In traditional catalyst kinetic characterization, the gas composition is known at the inlet, and the exit flow is measured to determine changes in concentration. As such, the progression of the chemistry within the catalyst is not known. Technological advances in electromagnetic and physical probes have made visualizing the evolution of the chemistry within catalyst samples a reality, as part of a methodology commonly known as spatial resolution. Herein, we discuss and evaluate the development of spatially resolved techniques, including the evolutions and achievements of this growing area of catalytic research. The impact of such techniques is discussed in terms of the invasiveness of physical probes on catalytic systems, as well as how experimentally obtained spatial profiles can be used in conjunction with kinetic modelling. Furthermore, some aims and aspirations for further evolution of spatially resolved techniques are considered.

Laser-driven x-ray and neutron source development for industrial applications of plasma accelerators

Pulsed beams of energetic x-rays and neutrons from intense laser interactions with solid foils are promising for applications where bright, small emission area sources, capable of multi-modal delivery are ideal. Possible end users of laser-driven multi-modal sources are those requiring advanced non-destructive inspection techniques in industry sectors of high value commerce such as aerospace, nuclear and advanced manufacturing. We report on experimental work that demonstrates multi-modal operation of high power laser-solid interactions for neutron and x-ray beam generation. Measurements and Monte Carlo radiation transport simulations show that neutron yield is increased by a factor ∼2 when a 1 mm copper foil is placed behind a 2 mm lithium foil, compared to using a 2 cm block of lithium only. We explore x-ray generation with a 10 picosecond drive pulse in order to tailor the spectral content for radiography with medium density alloy metals. The impact of using >1 ps pulse duration on laser-accelerated electron beam generation and transport is discussed alongside the optimisation of subsequent bremsstrahlung emission in thin, high atomic number target foils. X-ray spectra are deconvolved from spectrometer measurements and simulation data generated using the GEANT4 Monte Carlo code. We also demonstrate the unique capability of laser-driven x-rays in being able to deliver single pulse high spatial resolution projection imaging of thick metallic objects. Active detector radiographic imaging of industrially relevant sample objects with a 10 ps drive pulse is presented for the first time, demonstrating that features of 200 μm size are resolved when projected at high magnification.

ALMA-resolved salt emission traces the chemical footprint and inner wind morphology of VY Canis Majoris

(abreviated) We aim to study the inner-wind structure (R<250 Rstar) of the well-known red supergiant VY CMa. We analyse high spatial resolution (~0".24x0".13) ALMA Science Verification (SV) data in band 7 in which four thermal emission lines of gaseous sodium chloride (NaCl) are present at high signal-to-noise ratio. For the first time, the NaCl emission in the inner wind region of VY CMa is spatially resolved. The ALMA observations reveal the contribution of up to four different spatial regions. The NaCl emission pattern is different compared to the dust continuum and TiO2 emission already analysed from the ALMA SV data. The emission can be reconciled with an axisymmetric geometry, where the lower density polar/rotation axis has a position angle of ~50 degrees measured from north to east. However, this picture can not capture the full morphological diversity, and discrete mass ejection events need to be invoked to explain localized higher-density regions. The velocity traced by the gaseous NaCl line profiles is significantly lower than the average wind terminal velocity, and much slower than some of the fastest mass ejections, signalling a wide range of characteristic speeds for the mass loss. Gaseous NaCl is detected far beyond the main dust condensation region. Realising the refractory nature of this metal halide, this hints at a chemical process preventing all NaCl from condensing onto dust grains. We show that in the case of the ratio of the surface binding temperature to the grain temperature being ~50, only some 10% of NaCl remains in gaseous form, while for lower values of this ratio thermal desorption efficiently evaporates NaCl. Photodesorption by stellar photons seems not to be a viable explanation for the detection of gaseous NaCl at 220 Rstar from the central star, and instead, we propose shock-induced sputtering driven by localized mass ejection events as alternative.

Self-referencing thermometry in the nanoscale

Na ´ultima d´ecada emergiu uma linha de investiga¸c˜ao muito activa em term´ometros n˜ao invasivos e precisos que possam determinar temperatura `a escala nanom´etrica. Esta investiga¸c˜ao foi fortemente estimulada pelas numerosas solicita¸c˜oes da nanotecnologia e da biomedicina, por exemplo. Uma das abordagens mais promissoras prop˜oe o uso de i˜oes trivalentes de lantan´ıdeos que apresenta propriedades fotoluminescentes que dependem da temperatura. Neste trabalho demonstra-se que esta t´ecnica combina as vantagens de te um limite de detec¸c˜ao de 0.5 graus com sensibilidade at´e 4.5 % · K−1. Este term´ometro molecular pode ser processado em filmes finos ou nanopart´ıculas, abrindo os campos de aplica¸c˜ao a diferentes utiliza¸c˜oes. As nanopart´ıculas de s´ılica produzidas s˜ao caracterizadas na presen¸ca e na ausˆencia de i˜oes lantan´ıdeos. Sem o metal, as nanopart´ıculas de APTES/TEOS demonstram ser luminescentes sob excita¸c˜ao UV sem necessidade de utilizar qualquer tratamento t´ermico. O rendimento quˆantico de emiss˜ao depende apenas da propor¸c˜ao dos silanos e pode atingir o valor de 0.15 ± 0.02. A co-dopagem destas nanopart´ıculas com Eu3+ e Tb3+ permite obter sondas com resposta raciom´etrica, com a possibilidade de ajustar a gama de temperaturas de opera¸c˜ao e a sensibilidade, via desenho inteligente da matriz de suporte e dos ligandos de β-dicetona que est˜ao coordenados ao i˜ao met´alico. Quando processados como filmes, este term´ometro permite o mapeamento de temperaturas com resolu¸c˜ao espacial 1.8 μm. A racionaliza¸c˜ao da dependˆencia de temperatura ´e uma ferramenta ´util para desenvolver term´ometros que operam em gamas de temperatura espec´ıficos (e.g. gama de temperatura fisiol´ogica, 290-340 K) com sensibilidade acima de 0.5 % · K−1. A combina¸c˜ao de esfor¸cos de um grande n´umero de diversas disciplinas ir´a previsivelmente permitir o surgimento de term´ometros moleculares novos e sofisticados, preenchendo os principais requisitos das nanociencias.

Estudo de otimização de imagem e dose em mamografia digital

O objetivo deste trabalho foi efetuar um estudo de otimização dos parâmetros de exposição do mamógrafo digital da marca GE, modelo Senographe DS instalado no Instituto Português de Oncologia de Coimbra (IPOC) baseado numa análise contraste detalhe e avaliar o impacto do resultado obtido na melhoria da capacidade de detecção clínica das estruturas da mama. O mamógrafo em estudo dispõe de um sistema de controle automático da exposição designado por sistema AOP (Automatic Optimization of Parameters) que foi otimizado pelo fabricante a partir da razão contraste ruído. O sistema AOP é usado na prática clínica na quase totalidade dos exames de mamografia realizados no IPOC. Tendo em conta o tipo de estrutura que se pretende visualizar na mamografia, nomeadamente estruturas de baixo contraste como as massas e estruturas de dimensão submilimétrica como as microcalcificações, a análise contraste detalhe poderia constituir uma abordagem mais adequada para o estudo de otimização dos parâmetros de exposição uma vez que permitiria uma avaliação conjunta do contraste, da resolução espacial e do ruído da imagem. Numa primeira fase do trabalho foi efetuada a caracterização da prática clínica realizada no mamógrafo em estudo em termos de espessura de mama comprimida “típica”, dos parâmetros técnicos de exposição e das opções de processamento das imagens aplicadas pelo sistema AOP (combinação alvo/filtro, tensão aplicada na ampola - kVp e produto corrente-tempo da ampola - mAs). Numa segunda fase foi realizado um estudo de otimização da qualidade da imagem versus dose na perspectiva dos parâmetros físicos. Para tal foi efetuada uma análise contrastedetalhe no objeto simulador de mama CDMAM e usada uma figura de mérito definida a partir do IQFinv (inverted image quality figure) e da dose glandular média. Os resultados apontaram para uma diferença entre o ponto ótimo resultante do estudo de otimização e o ponto associado aos parâmetros de exposição escolhidos pelo sistema AOP, designadamente no caso da mama pequena. Sendo a qualidade da imagem na perspectiva clínica um conceito mais complexo cujo resultado da apreciação de uma imagem de boa qualidade deve ter em conta as diferenças entre observadores, foi efetuado na última parte deste trabalho um estudo do impacto clínico da proposta de otimização da qualidade de imagem. A partir das imagens realizadas com o objeto simulador antropomórfico TOR MAM simulando uma mama pequena, seis médicos(as) radiologistas com mais de 5 anos de experiência em mamografia avaliaram as imagens “otimizadas” obtidas utilizando-se os parâmetros técnicos de exposição resultantes do estudo de otimização e a imagem resultante da escolha do sistema AOP. A análise estatística das avaliações feitas pelos médicos indica que a imagem “otimizada” da mama pequena apresenta uma melhor visualização das microcalcificações sem perda da qualidade da imagem na deteção de fibras e de massas em comparação com a imagem “standard”. Este trabalho permitiu introduzir uma nova definição da figura de mérito para o estudo de otimização da qualidade da imagem versus dose em mamografia. Permitiu também estabelecer uma metodologia consistente que pode facilmente ser aplicada a qualquer outro mamógrafo, contribuindo para a área da otimização em mamografia digital que é uma das áreas mais relevantes no que toca à proteção radiológica do paciente.

Development of a single photon counting computed tomography system using MPGDs

The development of computed tomography systems with energy resolving detectors is a current challenge in medical physics and biomedical engineering. A computed tomography system of this kind allows getting complementary informations relatively to conventional systems, that can help the medical diagnosis, being of great interest in medicine. The work described in this thesis is related to the development of a computed tomography system using micropattern gaseous detectors, which allow storing, simultaneously, information about the interaction position and the energy of each single photon that interacts with the detector. This kind of detectors has other advantages concerning the cost and characteristics of operation when compared with solid state detectors. Tomographic acquisitions were performed using a MicroHole & Strip Plate based detector, which allowed reconstructing cross-sectional images using energy windows, applying the energy weighting technique and performing multi-slice and tri-dimensional reconstructions. The contrast-to-noise ratio was improved by 31% by applying the energy weighting technique, comparing with the corresponding image obtained with the current medical systems. A prototype of a computed tomography with flexibility to change the detector was developed, making it possible to apply different detectors based on Thick-COBRA. Several images acquired with these detectors are presented and demonstrate their applicability in X-ray imaging. When operating in NeCH4, the detector allowed a charge gain of 8 104, an energy resolution of 20% (full width at half maximum at 8 keV), a count rate of 1 106 Hz/mm2, a very stable operation (gain fluctuations below 5%) and a spacial resolution of 1.2 mm for an energy photon of 3.6 keV. Operating the detector in pure Kr allowed increasing the detection efficiency and achieving a charge gain of 2 104, an energy resolution of 32% (full width at half maximum at 22 keV), a count rate of 1 105 Hz/mm2, very stable operation and a spatial resolution of 500 m. The software already existing in the group was improved and tools to correct geometric misalignments of the system were also developed. The reconstructions obtained after geometrical correction are free of artefacts due to the referred misalignments.

Combining sampling, profiling, remote sensing and modelling to efficiently monitoring estuarine and coastal water quality

Monitoring of coastal and estuarine water quality has been traditionally performed by sampling with subsequent laboratory analysis. This has the disadvantages of low spatial and temporal resolution and high cost. In the last decades two alternative techniques have emerged to overcome this drawback: profiling and remote sensing. Profiling using multi-parameter sensors is now in a commercial stage. It can be used, tied to a boat, to obtain a quick “picture” of the system. The spatial resolution thus increases from single points to a line coincident with the boat track. The temporal resolution however remains unchanged since campaigns and resources involved are basically the same. The need for laboratory analysis was reduced but not eliminated because parameters like nutrients, microbiology or metals are still difficult to obtain with sensors and validation measurements are still needed. In the last years the improvement in satellite resolution has enabled its use for coastal and estuarine water monitoring. Although spatial coverage and resolution of satellite images in the present is already suitable to coastal and estuarine monitoring, temporal resolution is naturally limited to satellite passages and cloud cover. With this panorama the best approach to water monitoring is to integrate and combine data from all these sources. The natural tools to perform this integration are numerical models. Models benefit from the different sources of data to obtain a better calibration. After calibration they can be used to extend spatially and temporally the methods resolution. In Algarve (South of Portugal) a monitoring effort using this approach is being undertaken. The monitoring effort comprises five different locations including coastal waters, estuaries and coastal lagoons. The objective is to establish the base line situation to evaluate the impact of Waste Water Treatment Plants design and retrofitting. The field campaigns include monthly synoptic profiling, using an YSI 6600 multi-parameter system, laboratory analysis and fixed stations. The remote sensing uses ENVISAT\MERIS Level 2 Full Resolution data. This data is combined and used with the MOHID modelling system to obtain an integrate description of the systems. The results show the limitations of each method and the ability of the modelling system to integrate the results and to produce a comprehensive picture of the system.

Fast segmentation of 3D data using an octree

The algorithm developed uses an octree pyramid in which noise is reduced at the expense of the spatial resolution. At a certain level an unsupervised clustering without spatial connectivity constraints is applied. After the classification, isolated voxels and insignificant regions are removed by assigning them to their neighbours. The spatial resolution is then increased by the downprojection of the regions, level by level. At each level the uncertainty of the boundary voxels is minimised by a dynamic selection and classification of these, using an adaptive 3D filtering. The algorithm is tested using different data sets, including NMR data.

Partial volume correction using bimodal contrast agents in PET-MRI

Tese de mestrado integrado em Engenharia Biomédica e Biofísica, apresentada à Universidade de Lisboa, através da Faculdade de Ciências, 2014