Computational and Experimental Modelling of Mooring Lines Dynamics for Offshore Floating Wind Turbines

El cálculo de cargas de aerogeneradores flotantes requiere herramientas de simulación en el dominio del tiempo que consideren todos los fenómenos que afectan al sistema, como la aerodinámica, la dinámica estructural, la hidrodinámica, las estrategias de control y la dinámica de las líneas de fondeo. Todos estos efectos están acoplados entre sí y se influyen mutuamente. Las herramientas integradas se utilizan para calcular las cargas extremas y de fatiga que son empleadas para dimensionar estructuralmente los diferentes componentes del aerogenerador. Por esta razón, un cálculo preciso de las cargas influye de manera importante en la optimización de los componentes y en el coste final del aerogenerador flotante. En particular, el sistema de fondeo tiene gran impacto en la dinámica global del sistema. Muchos códigos integrados para la simulación de aerogeneradores flotantes utilizan modelos simplificados que no consideran los efectos dinámicos de las líneas de fondeo. Una simulación precisa de las líneas de fondeo dentro de los modelos integrados puede resultar fundamental para obtener resultados fiables de la dinámica del sistema y de los niveles de cargas en los diferentes componentes. Sin embargo, el impacto que incluir la dinámica de los fondeos tiene en la simulación integrada y en las cargas todavía no ha sido cuantificada rigurosamente. El objetivo principal de esta investigación es el desarrollo de un modelo dinámico para la simulación de líneas de fondeo con precisión, validarlo con medidas en un tanque de ensayos e integrarlo en un código de simulación para aerogeneradores flotantes. Finalmente, esta herramienta, experimentalmente validada, es utilizada para cuantificar el impacto que un modelos dinámicos de líneas de fondeo tienen en la computación de las cargas de fatiga y extremas de aerogeneradores flotantes en comparación con un modelo cuasi-estático. Esta es una información muy útil para los futuros diseñadores a la hora de decidir qué modelo de líneas de fondeo es el adecuado, dependiendo del tipo de plataforma y de los resultados esperados. El código dinámico de líneas de fondeo desarrollado en esta investigación se basa en el método de los Elementos Finitos, utilizando en concreto un modelo ”Lumped Mass” para aumentar su eficiencia de computación. Los experimentos realizados para la validación del código se realizaron en el tanque del École Céntrale de Nantes (ECN), en Francia, y consistieron en sumergir una cadena con uno de sus extremos anclados en el fondo del tanque y excitar el extremo suspendido con movimientos armónicos de diferentes periodos. El código demostró su capacidad para predecir la tensión y los movimientos en diferentes posiciones a lo largo de la longitud de la línea con gran precisión. Los resultados indicaron la importancia de capturar la dinámica de las líneas de fondeo para la predicción de la tensión especialmente en movimientos de alta frecuencia. Finalmente, el código se utilizó en una exhaustiva evaluación del efecto que la dinámica de las líneas de fondeo tiene sobre las cargas extremas y de fatiga de diferentes conceptos de aerogeneradores flotantes. Las cargas se calcularon para tres tipologías de aerogenerador flotante (semisumergible, ”spar-buoy” y ”tension leg platform”) y se compararon con las cargas obtenidas utilizando un modelo cuasi-estático de líneas de fondeo. Se lanzaron y postprocesaron más de 20.000 casos de carga definidos por la norma IEC 61400-3 siguiendo todos los requerimientos que una entidad certificadora requeriría a un diseñador industrial de aerogeneradores flotantes. Los resultados mostraron que el impacto de la dinámica de las líneas de fondeo, tanto en las cargas de fatiga como en las extremas, se incrementa conforme se consideran elementos situados más cerca de la plataforma: las cargas en la pala y en el eje sólo son ligeramente modificadas por la dinámica de las líneas, las cargas en la base de la torre pueden cambiar significativamente dependiendo del tipo de plataforma y, finalmente, la tensión en las líneas de fondeo depende fuertemente de la dinámica de las líneas, tanto en fatiga como en extremas, en todos los conceptos de plataforma que se han evaluado. ABSTRACT The load calculation of floating offshore wind turbine requires time-domain simulation tools taking into account all the phenomena that affect the system such as aerodynamics, structural dynamics, hydrodynamics, control actions and the mooring lines dynamics. These effects present couplings and are mutually influenced. The results provided by integrated simulation tools are used to compute the fatigue and ultimate loads needed for the structural design of the different components of the wind turbine. For this reason, their accuracy has an important influence on the optimization of the components and the final cost of the floating wind turbine. In particular, the mooring system greatly affects the global dynamics of the floater. Many integrated codes for the simulation of floating wind turbines use simplified approaches that do not consider the mooring line dynamics. An accurate simulation of the mooring system within the integrated codes can be fundamental to obtain reliable results of the system dynamics and the loads. The impact of taking into account the mooring line dynamics in the integrated simulation still has not been thoroughly quantified. The main objective of this research consists on the development of an accurate dynamic model for the simulation of mooring lines, validate it against wave tank tests and then integrate it in a simulation code for floating wind turbines. This experimentally validated tool is finally used to quantify the impact that dynamic mooring models have on the computation of fatigue and ultimate loads of floating wind turbines in comparison with quasi-static tools. This information will be very useful for future designers to decide which mooring model is adequate depending on the platform type and the expected results. The dynamic mooring lines code developed in this research is based in the Finite Element Method and is oriented to the achievement of a computationally efficient code, selecting a Lumped Mass approach. The experimental tests performed for the validation of the code were carried out at the `Ecole Centrale de Nantes (ECN) wave tank in France, consisting of a chain submerged into a water basin, anchored at the bottom of the basin, where the suspension point of the chain was excited with harmonic motions of different periods. The code showed its ability to predict the tension and the motions at several positions along the length of the line with high accuracy. The results demonstrated the importance of capturing the evolution of the mooring dynamics for the prediction of the line tension, especially for the high frequency motions. Finally, the code was used for an extensive assessment of the effect of mooring dynamics on the computation of fatigue and ultimate loads for different floating wind turbines. The loads were computed for three platforms topologies (semisubmersible, spar-buoy and tension leg platform) and compared with the loads provided using a quasi-static mooring model. More than 20,000 load cases were launched and postprocessed following the IEC 61400-3 guideline and fulfilling the conditions that a certification entity would require to an offshore wind turbine designer. The results showed that the impact of mooring dynamics in both fatigue and ultimate loads increases as elements located closer to the platform are evaluated; the blade and the shaft loads are only slightly modified by the mooring dynamics in all the platform designs, the tower base loads can be significantly affected depending on the platform concept and the mooring lines tension strongly depends on the lines dynamics both in fatigue and extreme loads in all the platform concepts evaluated.

Attenuated T2 relaxation by mutual cancellation of dipole–dipole coupling and chemical shift anisotropy indicates an avenue to NMR structures of very large biological macromolecules in solution

Fast transverse relaxation of 1H, 15N, and 13C by dipole-dipole coupling (DD) and chemical shift anisotropy (CSA) modulated by rotational molecular motions has a dominant impact on the size limit for biomacromolecular structures that can be studied by NMR spectroscopy in solution. Transverse relaxation-optimized spectroscopy (TROSY) is an approach for suppression of transverse relaxation in multidimensional NMR experiments, which is based on constructive use of interference between DD coupling and CSA. For example, a TROSY-type two-dimensional 1H,15N-correlation experiment with a uniformly 15N-labeled protein in a DNA complex of molecular mass 17 kDa at a 1H frequency of 750 MHz showed that 15N relaxation during 15N chemical shift evolution and 1HN relaxation during signal acquisition both are significantly reduced by mutual compensation of the DD and CSA interactions. The reduction of the linewidths when compared with a conventional two-dimensional 1H,15N-correlation experiment was 60% and 40%, respectively, and the residual linewidths were 5 Hz for 15N and 15 Hz for 1HN at 4°C. Because the ratio of the DD and CSA relaxation rates is nearly independent of the molecular size, a similar percentagewise reduction of the overall transverse relaxation rates is expected for larger proteins. For a 15N-labeled protein of 150 kDa at 750 MHz and 20°C one predicts residual linewidths of 10 Hz for 15N and 45 Hz for 1HN, and for the corresponding uniformly 15N,2H-labeled protein the residual linewidths are predicted to be smaller than 5 Hz and 15 Hz, respectively. The TROSY principle should benefit a variety of multidimensional solution NMR experiments, especially with future use of yet somewhat higher polarizing magnetic fields than are presently available, and thus largely eliminate one of the key factors that limit work with larger molecules.

Steric and mass-induced sea level variations in the Mediterranean Sea revisited

The total sea level variation (SLV) is the combination of steric and mass␣induced SLV, whose exact shares are key to understanding the oceanic response to climate system changes. Total SLV can be observed by radar altimetry satellites such as TOPEX/POSEIDON and Jason 1/2. The steric SLV can be computed through temperature and salinity profiles from in situ measurements or from ocean general circulation models (OGCM), which can assimilate the said observations. The mass-induced SLV can be estimated from its time-variable gravity (TVG) signals. We revisit this problem in the Mediterranean Sea estimating the observed, steric, and mass-induced SLV, for the latter we analyze the latest TVG data set from the GRACE (Gravity Recovery and Climate Experiment) satellite mission launched in 2002, which is 3.5 times longer than in previous studies, with the application of a two-stage anisotropic filter to reduce the noise in high-degree and -order spherical harmonic coefficients. We confirm that the intra-annual total SLV are only produced by water mass changes, a fact explained in the literature as a result of the wind field around the Gibraltar Strait. The steric SLV estimated from the residual of “altimetry minus GRACE” agrees in phase with that estimated from OGCMs and in situ measurements, although showing a higher amplitude. The net water fluxes through both the straits of Gibraltar and Sicily have also been estimated accordingly.

Accreting magnetars: a new type of high-mass X-ray binaries?

The discovery of very slow pulsations (Pspin =5560 s) has solved the long-standing question of the nature of the compact object in the high-mass X-ray binary 4U 2206+54 but has posed new ones. According to spin evolutionary models in close binary systems, such slow pulsations require a neutron star magnetic field strength larger than the quantum critical value of 4.4 × 1013 G, suggesting the presence of a magnetar. We present the first XMM–Newton observations of 4U 2206+54 and investigate its spin evolution. We find that the observed spin-down rate agrees with the magnetar scenario. We analyse Integral Spacecraft Gamma-Ray Imager (ISGRI)/INTErnational Gamma-RAy Laboratory (INTEGRAL) observations of 4U 2206+54 to search for the previously suggested cyclotron resonance scattering feature at ∼30 keV. We do not find a clear indication of the presence of the line, although certain spectra display shallow dips, not always at 30 keV. The association of these dips with a cyclotron line is very dubious because of its apparent transient nature. We also investigate the energy spectrum of 4U 2206+54 in the energy range 0.3–10 keV with unprecedented detail and report for the first time the detection of very weak 6.5 keV fluorescence iron lines. The photoelectric absorption is consistent with the interstellar value, indicating very small amount of local matter, which would explain the weakness of the florescence lines. The lack of matter locally to the source may be the consequence of the relatively large orbital separation of the two components of the binary. The wind would be too tenuous in the vicinity of the neutron star.

Probing large-scale wind structures in Vela X–1 using off-states with INTEGRAL

Vela X–1 is the prototype of the class of wind-fed accreting pulsars in high-mass X-ray binaries hosting a supergiant donor. We have analysed in a systematic way 10 years of INTEGRAL data of Vela X–1 (22–50 keV) and we found that when outside the X-ray eclipse, the source undergoes several luminosity drops where the hard X-rays luminosity goes below ∼3 × 1035 erg s−1, becoming undetected by INTEGRAL. These drops in the X-ray flux are usually referred to as ‘off-states’ in the literature. We have investigated the distribution of these off-states along the Vela X–1 ∼ 8.9 d orbit, finding that their orbital occurrence displays an asymmetric distribution, with a higher probability to observe an off-state near the pre-eclipse than during the post-eclipse. This asymmetry can be explained by scattering of hard X-rays in a region of ionized wind, able to reduce the source hard X-ray brightness preferentially near eclipse ingress. We associate this ionized large-scale wind structure with the photoionization wake produced by the interaction of the supergiant wind with the X-ray emission from the neutron star. We emphasize that this observational result could be obtained thanks to the accumulation of a decade of INTEGRAL data, with observations covering the whole orbit several times, allowing us to detect an asymmetric pattern in the orbital distribution of off-states in Vela X–1.

An XMM-Newton view of FeKα in high-mass X-ray binaries

We present a comprehensive analysis of the whole sample of available XMM-Newton observations of high-mass X-ray binaries (HMXBs) until August 2013, focusing on the FeKα emission line. This line is key to better understanding the physical properties of the material surrounding the X-ray source within a few stellar radii (the circumstellar medium). We collected observations from 46 HMXBs and detected FeKα in 21 of them. We used the standard classification of HMXBs to divide the sample into different groups. We find that (1) different classes of HMXBs display different qualitative behaviours in the FeKα spectral region. This is visible especially in SGXBs (showing ubiquitous Fe fluorescence but not recombination Fe lines) and in γ Cass analogues (showing both fluorescent and recombination Fe lines). (2) FeKα is centred at a mean value of 6.42 keV. Considering the instrumental and fits uncertainties, this value is compatible with ionization states that are lower than Fe xviii. (3) The flux of the continuum is well correlated with the flux of the line, as expected. Eclipse observations show that the Fe fluorescence emission comes from an extended region surrounding the X-ray source. (4) We observe an inverse correlation between the X-ray luminosity and the equivalent width of FeKα (EW). This phenomenon is known as the X-ray Baldwin effect. (5) FeKα is narrow (σline< 0.15 keV), reflecting that the reprocessing material does not move at high speeds. We attempt to explain the broadness of the line in terms of three possible broadening phenomena: line blending, Compton scattering, and Doppler shifts (with velocities of the reprocessing material V ~ 1000 km s-1). (6) The equivalent hydrogen column (NH) directly correlates to the EW of FeKα, displaying clear similarities to numerical simulations. It highlights the strong link between the absorbing and the fluorescent matter. (7) The observed NH in supergiant X-ray binaries (SGXBs) is in general higher than in supergiant fast X-ray transients (SFXTs). We suggest two possible explanations: different orbital configurations or a different interaction compact object – wind. (8) Finally, we analysed the sources IGR J16320-4751 and 4U 1700-37 in more detail, covering several orbital phases. The observed variation in NH between phases is compatible with the absorption produced by the wind of their optical companions. The results clearly point to a very important contribution of the donor’s wind in the FeKα emission and the absorption when the donor is a supergiant massive star.

A VLT/FLAMES survey for massive binaries in Westerlund 1. IV. Wd1-5 – binary product and a pre-supernova companion for the magnetar CXOU J1647-45?

Context. The first soft gamma-ray repeater was discovered over three decades ago, and was subsequently identified as a magnetar, a class of highly magnetised neutron star. It has been hypothesised that these stars power some of the brightest supernovae known, and that they may form the central engines of some long duration gamma-ray bursts. However there is currently no consenus on the formation channel(s) of these objects. Aims. The presence of a magnetar in the starburst cluster Westerlund 1 implies a progenitor with a mass ≥40 M⊙, which favours its formation in a binary that was disrupted at supernova. To test this hypothesis we conducted a search for the putative pre-SN companion. Methods. This was accomplished via a radial velocity survey to identify high-velocity runaways, with subsequent non-LTE model atmosphere analysis of the resultant candidate, Wd1-5. Results. Wd1-5 closely resembles the primaries in the short-period binaries, Wd1-13 and 44, suggesting a similar evolutionary history, although it currently appears single. It is overluminous for its spectroscopic mass and we find evidence of He- and N-enrichement, O-depletion, and critically C-enrichment, a combination of properties that is difficult to explain under single star evolutionary paradigms. We infer a pre-SN history for Wd1-5 which supposes an initial close binary comprising two stars of comparable (~ 41 M⊙ + 35 M⊙) masses. Efficient mass transfer from the initially more massive component leads to the mass-gainer evolving more rapidly, initiating luminous blue variable/common envelope evolution. Reverse, wind-driven mass transfer during its subsequent WC Wolf-Rayet phase leads to the carbon pollution of Wd1-5, before a type Ibc supernova disrupts the binary system. Under the assumption of a physical association between Wd1-5 and J1647-45, the secondary is identified as the magnetar progenitor; its common envelope evolutionary phase prevents spin-down of its core prior to SN and the seed magnetic field for the magnetar forms either in this phase or during the earlier episode of mass transfer in which it was spun-up. Conclusions. Our results suggest that binarity is a key ingredient in the formation of at least a subset of magnetars by preventing spin-down via core-coupling and potentially generating a seed magnetic field. The apparent formation of a magnetar in a Type Ibc supernova is consistent with recent suggestions that superluminous Type Ibc supernovae are powered by the rapid spin-down of these objects.

Measuring the stellar wind parameters in IGR J17544-2619 and Vela X-1 constrains the accretion physics in supergiant fast X-ray transient and classical supergiant X-ray binaries

Context. Classical supergiant X-ray binaries (SGXBs) and supergiant fast X-ray transients (SFXTs) are two types of high-mass X-ray binaries (HMXBs) that present similar donors but, at the same time, show very different behavior in the X-rays. The reason for this dichotomy of wind-fed HMXBs is still a matter of debate. Among the several explanations that have been proposed, some of them invoke specific stellar wind properties of the donor stars. Only dedicated empiric analysis of the donors’ stellar wind can provide the required information to accomplish an adequate test of these theories. However, such analyses are scarce. Aims. To close this gap, we perform a comparative analysis of the optical companion in two important systems: IGR J17544-2619 (SFXT) and Vela X-1 (SGXB). We analyze the spectra of each star in detail and derive their stellar and wind properties. As a next step, we compare the wind parameters, giving us an excellent chance of recognizing key differences between donor winds in SFXTs and SGXBs. Methods. We use archival infrared, optical and ultraviolet observations, and analyze them with the non-local thermodynamic equilibrium (NLTE) Potsdam Wolf-Rayet model atmosphere code. We derive the physical properties of the stars and their stellar winds, accounting for the influence of X-rays on the stellar winds. Results. We find that the stellar parameters derived from the analysis generally agree well with the spectral types of the two donors: O9I (IGR J17544-2619) and B0.5Iae (Vela X-1). The distance to the sources have been revised and also agree well with the estimations already available in the literature. In IGR J17544-2619 we are able to narrow the uncertainty to d = 3.0 ± 0.2 kpc. From the stellar radius of the donor and its X-ray behavior, the eccentricity of IGR J17544-2619 is constrained to e< 0.25. The derived chemical abundances point to certain mixing during the lifetime of the donors. An important difference between the stellar winds of the two stars is their terminal velocities (ν∞ = 1500 km s-1 in IGR J17544-2619 and ν∞ = 700 km s-1 in Vela X-1), which have important consequences on the X-ray luminosity of these sources. Conclusions. The donors of IGR J17544-2619 and Vela X-1 have similar spectral types as well as similar parameters that physically characterize them and their spectra. In addition, the orbital parameters of the systems are similar too, with a nearly circular orbit and short orbital period. However, they show moderate differences in their stellar wind velocity and the spin period of their neutron star which has a strong impact on the X-ray luminosity of the sources. This specific combination of wind speed and pulsar spin favors an accretion regime with a persistently high luminosity in Vela X-1, while it favors an inhibiting accretion mechanism in IGR J17544-2619. Our study demonstrates that the relative wind velocity is critical in class determination for the HMXBs hosting a supergiant donor, given that it may shift the accretion mechanism from direct accretion to propeller regimes when combined with other parameters.

(Table 1) Mass-accumulation rates for the non-authigenic, inorganic, crystalline, eolian component of sediments recovered at DSDP Hole 62-463

Elevated regions in the central parts of ocean basins are excellent for study of accumulation of eolian material. The mass-accumulation rates of this sediment component appear to reflect changes in the influx of volcanic materials through the Early Cretaceous to Recent history of Deep Sea Drilling Project Site 463, on the Mid-Pacific Mountains. Four distinct episodes of eolian accumulation occurred during the Cretaceous: two periods of moderate accumulation, averaging about 0.2 to 0.3 g/cm**2/10**3 yr, 67 to 70.5 m.y. ago and 91 to 108 m.y. ago; a period of low accumulation, approximately 0.03 g/cm**2/10**3 yr, 70.5 to 90 m.y. ago; and a period of high accumulation, about 0.9 g/cm**2/10**3 yr, 109 to 117 m.y. ago (bottom of the hole). Much of the Cenozoic section is missing from Site 463. Upper Miocene to Recent sediments record an upward increase in accumulation rates, from less than 0.01 to about 0.044 g/cm**2/10**3 yr. The late Pliocene-Pleistocene peak may reflect the change to glacial-wind regimes, as well as an increase in volcanic source materials.

Eolian mass-accumulation rates and grain size at DSDP Sites 86-576 and 86-578

The mass-accumulation rate and grain size of the total eolian component of North Pacific pelagic clays at Deep Sea Drilling Project Sites 576 and 578 have been used to evaluate changes in eolian sedimentation and the intensity of atmospheric circulation that have occurred during the past 70 m.y. Eolian deposition, an indicator of source area aridity, was low in the Paleocene, Eocene, and Oligocene, apparently reflecting the humid environments of that time as well as the lack of glacial erosion products. A general increase in eoiian accumulation in the Miocene apparently reflects the relative increase in global aridity during the latter part of the Cenozoic. A dramatic increase in eolian accumulation rates in the Pliocene reflects the increased aridity and availability of glacial erosion products associated with Northern Hemisphere glaciation 2.5 m.y. ago. Eolian grain size, an indicator of wind intensity, suggests that Late Cretaceous wind strength was comparable to present-day wind strength. A sharp decrease in eolian grain size across the Paleocene/Eocene boundary is not readily interpreted, but may indicate a significant reduction in the intensity of atmospheric circulation at that time. Fine eolian grain size and low accumulation rates in the Eocene and early Oligocene are in agreement with low early Tertiary thermal gradients and less vigorous atmospheric circulation. Large increases in grain size during the Oligocene, mid-to-late Miocene, and Pliocene appear to be a response to steepening thermal gradients resulting from increasing polar isolation.

Table 2: Firn temperature, ram resistance, mass balance and submergence velo city measurements on Colle Gnifetti

Distributional patterns of glaciological parameters at the Colle Gnifetti core drilling site are described and their interrelationships are brietly discussed. Observations within a stake network established in 1980 furnish information about snow accumulation (short term balance), submergence velocity of ice tlow (long term balance), ram hardness (melt layer stratigraphy), and firn temperature. In addition, a numerical model was used to estimate local variations of available radiant energy. Melt layer formation is considerably more intensive on the south facing parts of the firn saddie where incoming radiation is high. These melt layers seem to effectively protect some of the fallen snow from wind erosion. As a result, balance ist up to one order of magnitude larger on south facing slopes. Heat applied to the surface is therefore positively correlated with balance, whereas the relation between solar radiation and firn temperature is less dear. Distributional patterns of submergence velocity confirm that the observed spatial variability of surface balance is representative for longer time periods and greatly intluences the time scale and the stratigraphy of firn and ice cores from Colle Gnifetti.

The modern band: a treatise on wind instruments, symphony band and military band,

Bk.I. pt.1. Wind instruments. - pt.2. Band instrumentation. -Bk.II. Score examples.

Formic acid mixing ratio (CIMS) and 3D wind data, plus SOSAA model runs -- Hyytiälä, Finland, Apr-Jun 2014

These are data of eddy covariance flux measurements of formic acid (HCOOH), performed by a chemical ionization mass spectrometer (CIMS) over a boreal forest canopy in Hyytiälä, Finland, in spring/summer 2014. Results from the 1-D chemical transport model runs using SOSAA (Simulate Organic vapours, Sulphuric Acid and Aerosols) are included as well. The data accompany a submission of a manuscript to Geophysical Research Letters for consideration for publication (Schobesberger et al.).

The effects of J-gate potential and interfaces on donor exchange coupling in the Kane quantum computer architecture

We calculate the electron exchange coupling for a phosphorus donor pair in silicon perturbed by a J-gate potential and the boundary effects of the silicon host geometry. In addition to the electron-electron exchange interaction we also calculate the contact hyperfine interaction between the donor nucleus and electron as a function of the varying experimental conditions. Donor separation, depth of the P nuclei below the silicon oxide layer and J-gate voltage become decisive factors in determining the strength of both the exchange coupling and hyperfine interaction-both crucial components for qubit operations in the Kane quantum computer. These calculations were performed using an anisotropic effective-mass Hamiltonian approach. The behaviour of the donor exchange coupling as a function of the parameters varied in this work provides relevant information for the experimental design of these devices.

Challenges for simultaneous nitrification, denitrification, and phosphorus removal in microbial aggregates: mass transfer limitation and nitrous oxide production

The microbial community composition and activity was investigated in aggregates from a lab-scale bioreactor, in which nitrification, denitrification and phosphorus removal occurred simultaneously. The biomass was highly enriched for polyphosphate accumulating organisms facilitating complete removal of phosphorus from the bulk liquid; however, some inorganic nitrogen still remained at the end of the reactor cycle. This was ascribed to incomplete coupling of nitrification and denitrification causing NO3- accumulation. After 2 h of aeration, denitrification was dependent on the activity of nitrifying bacteria facilitating the formation of anoxic zones in the aggregates; hence, denitrification could not occur without simultaneous nitrification towards the end of the reactor cycle. Nitrous oxide was identified as a product of denitrification, when based on stored PHA as carbon source. This observation is of critical importance to the outlook of applying PHA-driven denitrification in activated sludge processes. (c) 2004 Federation of European Microbiological Societies. Published by Elsevier B.V. All rights reserved.