Urban energy fluxes in Helsinki and their high frequency spectral corrections

Inadvertent climate modification has led to an increase in urban temperatures compared to the surrounding rural area. The main reason for the temperature rise is the altered energy portioning of input net radiation to heat storage and sensible and latent heat fluxes in addition to the anthropogenic heat flux. The heat storage flux and anthropogenic heat flux have not yet been determined for Helsinki and they are not directly measurable. To the contrary, turbulent fluxes of sensible and latent heat in addition to net radiation can be measured, and the anthropogenic heat flux together with the heat storage flux can be solved as a residual. As a result, all inaccuracies in the determination of the energy balance components propagate to the residual term and special attention must be paid to the accurate determination of the components. One cause of error in the turbulent fluxes is the fluctuation attenuation at high frequencies which can be accounted for by high frequency spectral corrections. The aim of this study is twofold: to assess the relevance of high frequency corrections to water vapor fluxes and to assess the temporal variation of the energy fluxes. Turbulent fluxes of sensible and latent heat have been measured at SMEAR III station, Helsinki, since December 2005 using the eddy covariance technique. In addition, net radiation measurements have been ongoing since July 2007. The used calculation methods in this study consist of widely accepted eddy covariance data post processing methods in addition to Fourier and wavelet analysis. The high frequency spectral correction using the traditional transfer function method is highly dependent on relative humidity and has an 11% effect on the latent heat flux. This method is based on an assumption of spectral similarity which is shown not to be valid. A new correction method using wavelet analysis is thus initialized and it seems to account for the high frequency variation deficit. Anyhow, the resulting wavelet correction remains minimal in contrast to the traditional transfer function correction. The energy fluxes exhibit a behavior characteristic for urban environments: the energy input is channeled to sensible heat as latent heat flux is restricted by water availability. The monthly mean residual of the energy balance ranges from 30 Wm-2 in summer to -35 Wm-2 in winter meaning a heat storage to the ground during summer. Furthermore, the anthropogenic heat flux is approximated to be 50 Wm-2 during winter when residential heating is important.

An assessment of observed vertical flux divergence in long-term eddy-covariance measurements over two Midwestern forest ecosystems

Vertical divergence of CO2 fluxes is observed over two Midwestern AmeriFlux forest sites. The differences in ensemble averaged hourly CO2 fluxes measured at two heights above canopy are relatively small (0.2–0.5 μmol m−2 s−1), but they are the major contributors to differences (76–256 g C m−2 or 41.8–50.6%) in estimated annual net ecosystem exchange (NEE) in 2001. A friction velocity criterion is used in these estimates but mean flow advection is not accounted for. This study examines the effects of coordinate rotation, averaging time period, sampling frequency and co-spectral correction on CO2 fluxes measured at a single height, and on vertical flux differences measured between two heights. Both the offset in measured vertical velocity and the downflow/upflow caused by supporting tower structures in upwind directions lead to systematic over- or under-estimates of fluxes measured at a single height. An offset of 1 cm s−1 and an upflow/downflow of 1° lead to 1% and 5.6% differences in momentum fluxes and nighttime sensible heat and CO2 fluxes, respectively, but only 0.5% and 2.8% differences in daytime sensible heat and CO2 fluxes. The sign and magnitude of both offset and upflow/downflow angle vary between sonic anemometers at two measurement heights. This introduces a systematic and large bias in vertical flux differences if these effects are not corrected in the coordinate rotation. A 1 h averaging time period is shown to be appropriate for the two sites. In the daytime, the absolute magnitudes of co-spectra decrease with height in the natural frequencies of 0.02–0.1 Hz but increase in the lower frequencies (<0.01 Hz). Thus, air motions in these two frequency ranges counteract each other in determining vertical flux differences, whose magnitude and sign vary with averaging time period. At night, co-spectral densities of CO2 are more positive at the higher levels of both sites in the frequency range of 0.03–0.4 Hz and this vertical increase is also shown at most frequencies lower than 0.03 Hz. Differences in co-spectral corrections at the two heights lead to a positive shift in vertical CO2 flux differences throughout the day at both sites. At night, the vertical CO2 flux differences between two measurement heights are 20–30% and 40–60% of co-spectral corrected CO2 fluxes measured at the lower levels of the two sites, respectively. Vertical differences of CO2 flux are relatively small in the daytime. Vertical differences in estimated mean vertical advection of CO2 between the two measurement heights generally do not improve the closure of the 1D (vertical) CO2 budget in the air layer between the two measurement heights. This may imply the significance of horizontal advection. However, a reliable assessment of mean advection contributions in annual NEE estimate at these two AmeriFlux sites is currently an unsolved problem.

Physical and chemical matrix effects in soil carbon quantification using laser-induced breakdown spectroscopy

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Factor analysis used to peak position optimization and spectral interference correction in inductively coupled plasma atomic emission spectrometry

Target transformation factor analysis was used to correct spectral interference in inductively coupled plasma atomic emission spectrometry (ICP-BES) for the determination of rare earth impurities in high purity thulium oxide. Data matrix was constructed with pure and mixture vectors and background vector. A method based on an error evaluation function was proposed to optimize the peak position, so the influence of the peak position shift in spectral scans on the determination was eliminated or reduced. Satisfactory results were obtained using factor analysis and the proposed peak position optimization method.

APPLICATION OF FACTOR-ANALYSIS TO CORRECTION OF SPECTRAL OVERLAP INTERFERENCE IN INDUCTIVELY COUPLED PLASMA ATOMIC EMISSION-SPECTROMETRY

Correction of spectral overlap interference in inductively coupled plasma atomic emission spectrometry by factor analysis is attempted. For the spectral overlap of two known lines, a data matrix can be composed from one or two pure spectra and a spectrum of the mixture. The data matrix is decomposed into a spectra matrix and a concentration matrix by target transformation factor analysis. The component concentration of interest in a binary mixture is obtained from the concentration matrix and interference from the other component is eliminated. This method is applied to correcting spectral interference of yttrium on the determination of copper and aluminium: satisfactory results are obtained. This method may also be applied to correcting spectral overlap interference for more than two lines. Like other methods of correcting spectral interferences, factor analysis can only be used for additive spectral overlap. Results obtained from measurements on copper/yttrium mixtures with different white noise added show that random errors in measurement data do not significantly affect the results of the correction method.

3D refraction correction and extraction of clinical parameters from spectral domain optical coherence tomography of the cornea.

Capable of three-dimensional imaging of the cornea with micrometer-scale resolution, spectral domain-optical coherence tomography (SDOCT) offers potential advantages over Placido ring and Scheimpflug photography based systems for accurate extraction of quantitative keratometric parameters. In this work, an SDOCT scanning protocol and motion correction algorithm were implemented to minimize the effects of patient motion during data acquisition. Procedures are described for correction of image data artifacts resulting from 3D refraction of SDOCT light in the cornea and from non-idealities of the scanning system geometry performed as a pre-requisite for accurate parameter extraction. Zernike polynomial 3D reconstruction and a recursive half searching algorithm (RHSA) were implemented to extract clinical keratometric parameters including anterior and posterior radii of curvature, central cornea optical power, central corneal thickness, and thickness maps of the cornea. Accuracy and repeatability of the extracted parameters obtained using a commercial 859nm SDOCT retinal imaging system with a corneal adapter were assessed using a rigid gas permeable (RGP) contact lens as a phantom target. Extraction of these parameters was performed in vivo in 3 patients and compared to commercial Placido topography and Scheimpflug photography systems. The repeatability of SDOCT central corneal power measured in vivo was 0.18 Diopters, and the difference observed between the systems averaged 0.1 Diopters between SDOCT and Scheimpflug photography, and 0.6 Diopters between SDOCT and Placido topography.

Towards Multi-Sensor Spectral Alignment Through Post Measurement Calibration Correction

Semiconductor manufactures are increasing reliant on optical emission spectroscopy (OES) to source information on plasma characteristics and process change. However, nonlinearities in the response of OES sensors and errors in their calibration lead to discrepancies in observed wavelength detector response. This paper presents a technique for the retrospective spectral calibration of multiple OES sensors. Underlying methodology is given, and alignment performance is evaluated using OES recordings from a semiconductor plasma process. The paper concludes with a discussion of results and suggests avenues for future work.

Mathematical equation correction to spectral and transport interferences in high-resolution continuum source flame atomic absorption spectrometry: determination of lead in phosphoric acid

Um método de correção de interferência espectral e de transporte é proposto, e foi aplicado para minimizar interferências por moléculas de PO produzidas em chama ar-acetileno e de transporte causada pela variação da concentração de ácido fosfórico. Átomos de Pb e moléculas de PO absorvem a 217,0005 nm, então Atotal217,0005 nm = A Pb217,0005 nm + A PO217,0005 nm. Monitorando o comprimento de onda alternativo de PO em 217,0458 nm, é possível calcular a contribuição relativa de PO na absorbância total a 217,0005 nm: A Pb217,0005 nm = Atotal217,0005 nm - A PO217,0005 nm = Atotal217,0005 nm - k (A PO217,0458 nm). O fator de correção k é a razão entre os coeficientes angulares de duas curvas analíticas para P obtidas a 217,0005 e 217,0458 nm (k = b217,0005 nm/b217,0458 nm). Fixando-se a taxa de aspiração da amostra em 5,0 ml min-1, e integrando-se a absorbância no comprimento de onda a 3 pixels, curvas analíticas para Pb (0,1 - 1,0 mg L-1) foram obtidas com coeficientes de correlação típicos > 0,9990. As correlações lineares entre absorbância e concentração de P nos comprimentos de onda 217,0005 e 217,0458 foram > 0,998. O limite de detecção de Pb foi 10 µg L-1. O método de correção proposto forneceu desvios padrão relativos (n=12) de 2,0 a 6,0%, ligeiramente menores que os obtidos sem correção (1,4-4,3%). As recuperações de Pb adicionado às amostras de ácido fosfórico variaram de 97,5 a 100% (com correção pelo método proposto) e de 105 a 230% (sem correção).

A simple correction for multiple testing for single-nucleotide polymorphisms in linkage disequilibrium with each other

In this report, we describe a simple correction for multiple testing of single-nucleotide polymorphisms (SNPs) in linkage disequilibrium (LD) with each other, on the basis of the spectral decomposition (SpD) of matrices of pairwise LD between SNPs. This method provides a useful alternative to more computationally intensive permutation tests. Additionally, output from SNPSpD includes eigenvalues, principal-component coefficients, and factor "loadings" after varimax rotation, enabling the selection of a subset of SNPs that optimize the information in a genomic region.

Correction of Ocular Artifacts in EEG Recordings using Empirical Mode Decomposition

EEG recordings are often contaminated with ocular artifacts such as eye blinks and eye movements. These artifacts may obscure underlying brain activity in the electroencephalogram (EEG) data and make the analysis of the data difficult. In this paper, we explore the use of empirical mode decomposition (EMD) based filtering technique to correct the eye blinks and eye movementartifacts in single channel EEG data. In this method, the single channel EEG data containing ocular artifact is segmented such that the artifact in each of the segment is considered as some type of slowly varying trend in the dataand the EMD is used to remove the trend. The filtering is done using partial reconstruction from components of the decomposition. The method is completely data dependent and hence adaptive and nonlinear. Experimental results are provided to check the applicability of the method on real EEG data and the results are quantified using power spectral density (PSD) as a measure. The method has given fairlygood results and does not make use of any preknowledge of artifacts or the EEG data used.

Autoregressive spectral array for graphical display of EEG data

A graphical display of the frequency content of background,electroencephalogram (EEG) activity is obtained by calculating the spectral estimates using autocorrelation autoregressive method and the classical Fourier transform method, Display of spectral content of consecutive data segments is made using hidden-line suppression technique so as to get a spectral array, The autoregressive spectral array (ASA) is found to be sensitive to baseline drift, Following baseline correction the autoregressive technique is found to be superior to the Fourier method of compressed spectral array (CSA) in detecting the transitions in the frequencies of the signal. The smoothed ASA gives a better picture of transitions and changes in the background activity, The ASA can be made to adapt to specific changes of dominant frequencies while eliminating unnecessary peaks in the spectrum. The utility,of the ASA for background EEG analysis is discussed,

Accelerogram processing using reliability bounds and optimal correction methods

This study addresses the problem of obtaining reliable velocities and displacements from accelerograms, a concern which often arises in earthquake engineering. A closed-form acceleration expression with random parameters is developed to test any strong-motion accelerogram processing method. Integration of this analytical time history yields the exact velocities, displacements and Fourier spectra. Noise and truncation can also be added. A two-step testing procedure is proposed and the original Volume II routine is used as an illustration. The main sources of error are identified and discussed. Although these errors may be reduced, it is impossible to extract the true time histories from an analog or digital accelerogram because of the uncertain noise level and missing data. Based on these uncertainties, a probabilistic approach is proposed as a new accelerogram processing method. A most probable record is presented as well as a reliability interval which reflects the level of error-uncertainty introduced by the recording and digitization process. The data is processed in the frequency domain, under assumptions governing either the initial value or the temporal mean of the time histories. This new processing approach is tested on synthetic records. It induces little error and the digitization noise is adequately bounded. Filtering is intended to be kept to a minimum and two optimal error-reduction methods are proposed. The "noise filters" reduce the noise level at each harmonic of the spectrum as a function of the signal-to-noise ratio. However, the correction at low frequencies is not sufficient to significantly reduce the drifts in the integrated time histories. The "spectral substitution method" uses optimization techniques to fit spectral models of near-field, far-field or structural motions to the amplitude spectrum of the measured data. The extremes of the spectrum of the recorded data where noise and error prevail are then partly altered, but not removed, and statistical criteria provide the choice of the appropriate cutoff frequencies. This correction method has been applied to existing strong-motion far-field, near-field and structural data with promising results. Since this correction method maintains the whole frequency range of the record, it should prove to be very useful in studying the long-period dynamics of local geology and structures.

ON THE SIGNIFICANCE OF WAVELENGTH POSITIONING ACCURACY IN MULTICOMPONENT ANALYSIS TECHNIQUES FOR THE CORRECTION OF LINE INTERFERENCES IN INDUCTIVELY COUPLED PLASMA ATOMIC EMISSION-SPECTROMETRY

The present paper reports some definite evidence for the significance of wavelength positioning accuracy in multicomponent analysis techniques for the correction of line interferences in inductively coupled plasma atomic emission spectrometry (ICP-AES). Using scanning spectrometers commercially available today, a large relative error, DELTA(A) may occur in the estimated analyte concentration, owing to wavelength positioning errors, unless a procedure for data processing can eliminate the problem of optical instability. The emphasis is on the effect of the positioning error (deltalambda) in a model scan, which is evaluated theoretically and determined experimentally. A quantitative relation between DELTA(A) and deltalambda, the peak distance, and the effective widths of the analysis and interfering lines is established under the assumption of Gaussian line profiles. The agreement between calculated and experimental DELTA(A) is also illustrated. The DELTA(A) originating from deltalambda is independent of the net analyte/interferent signal ratio; this contrasts with the situation for the positioning error (dlambda) in a sample scan, where DELTA(A) decreases with an increase in the ratio. Compared with dlambda, the effect of deltalambda is generally less significant.

WAVELENGTH POSITIONING ERRORS AND TRUE DETECTION LIMIT IN RELATION TO MULTICOMPONENT ANALYSIS TECHNIQUES FOR THE CORRECTION OF LINE INTERFERENCES IN INDUCTIVELY COUPLED PLASMA ATOMIC EMISSION-SPECTROMETRY

The present paper deals with the evaluation of the relative error (DELTA(A)) in estimated analyte concentrations originating from the wavelength positioning error in a sample scan when multicomponent analysis (MCA) techniques are used for correcting line interferences in inductively coupled plasma atomic emission spectrometry. In the theoretical part, a quantitative relation of DELTA(A) with the extent of line overlap, bandwidth and the magnitude of the positioning error is developed under the assumption of Gaussian line profiles. The measurements of eleven samples covering various typical line interferences showed that the calculated DELTA(A) generally agrees well with the experimental one. An expression of the true detection limit associated with MCA techniques was thus formulated. With MCA techniques, the determination of the analyte and interferent concentrations depend on each other while with conventional correction techniques, such as the three-point method, the estimate of interfering signals is independent of the analyte signals. Therefore. a given positioning error results in a larger DELTA(A) and hence a higher true detection limit in the case of MCA techniques than that in the case of conventional correction methods. although the latter could be a reasonable approximation of the former when the peak distance expressed in the effective width of the interfering line is larger than 0.4. In the light of the effect of wavelength positioning errors, MCA techniques have no advantages over conventional correction methods unless the former can bring an essential reduction ot the positioning error.

Coherent wavelength division multiplexing: a novel transmission format for high spectral density optical communication networks

In this thesis a novel transmission format, named Coherent Wavelength Division Multiplexing (CoWDM) for use in high information spectral density optical communication networks is proposed and studied. In chapter I a historical view of fibre optic communication systems as well as an overview of state of the art technology is presented to provide an introduction to the subject area. We see that, in general the aim of modern optical communication system designers is to provide high bandwidth services while reducing the overall cost per transmitted bit of information. In the remainder of the thesis a range of investigations, both of a theoretical and experimental nature are carried out using the CoWDM transmission format. These investigations are designed to consider features of CoWDM such as its dispersion tolerance, compatibility with forward error correction and suitability for use in currently installed long haul networks amongst others. A high bit rate optical test bed constructed at the Tyndall National Institute facilitated most of the experimental work outlined in this thesis and a collaboration with France Telecom enabled long haul transmission experiments using the CoWDM format to be carried out. An amount of research was also carried out on ancillary topics such as optical comb generation, forward error correction and phase stabilisation techniques. The aim of these investigations is to verify the suitability of CoWDM as a cost effective solution for use in both current and future high bit rate optical communication networks