Procés de destintatge del paper per flotació. Avaluació de l'eficàcia de l'eliminació de tinta

La tesi doctoral presentada té com a objectius principals l'estudi de les etapes fonamentals de desintegració i flotació en un procés de destintatge de papers vell de qualitats elevades per a poder millorar l'eficàcia d'aquestes etapes clau. Conté una revisió teòrica completa i molt actualitzada del procés de desintegració i flotació tant a nivell macroscòpic com microscòpic. La metodologia de treball en el laboratori, la posada a punt dels aparells, així com les anàlisis efectuades per a valorar la resposta del procés (anàlisi de blancor, anàlisi d'imatge i anàlisi de la concentració efectiva de tinta residual) estan descrites en el capítol de material i mètodes. La posada en marxa permet obtenir unes primeres conclusions respecte la necessitat de treballar amb una matèria primera homogènia i respecte la no significació de la temperatura de desintegració dins l'interval de treball permès al laboratori (20-50°C). L'anàlisi de les variables mecàniques de desintegració: consistència de desintegració (c), velocitat d'agitació en la desintegració (N) i temps de desintegració (t), permet de discernir que la consistència de desintegració és una variable fonamental. El valor de consistència igual al 10% marca el límit d'existència de les forces d'impacte mecànic en la suspensió fibrosa. A consistències superiors, les forces viscoses i d'acceleració dominen l'etapa de desintegració. Existeix una interacció entre la consistència i el temps de desintegració, optimitzant-se aquesta darrera variable en funció del valor de la consistència. La velocitat d'agitació és significativa només per a valors de consistència de desintegració inferiors al 10%. En aquests casos, incrementar el valor de N de 800 a 1400 rpm representa una disminució de 14 punts en el factor de destintabilitat. L'anàlisi de les variables químiques de desintegració: concentració de silicat sòdic (% Na2SiO3), peròxid d'hidrogen (% H2O2) i hidròxid sòdic (% Na2OH), proporciona resultats força significatius. El silicat sòdic presenta un efecte altament dispersant corroborat per les corbes de distribució dels diàmetres de partícula de tinta obtingudes mitjançant anàlisi d'imatges. L'hidròxid sòdic també presenta un efecte dispersant tot i que no és tant important com el del silicat sòdic. Aquests efectes dispersants són deguts principalment a l'increment de les repulsions electrostàtiques que aporten a la suspensió fibrosa aquests reactius químics fent disminuir l'eficàcia d'eliminació de l'etapa de flotació. El peròxid d'hidrogen utilitzat generalment com agent blanquejant, actua en aquests casos com a neutralitzador dels grups hidroxil provinents tant del silicat sòdic com de l'hidròxid sòdic, disminuint la repulsió electrostàtica dins la suspensió. Amb l'anàlisi de les variables hidrodinàmiques de flotació: consistència de flotació (c), velocitat d'agitació durant la flotació (N) i cabal d'aire aplicat (q), s'aconsegueix la seva optimització dins el rang de treball permès al laboratori. Valors elevats tant de la velocitat d'agitació com del cabal d'aire aplicat durant la flotació permeten eliminar majors quantitats de tinta. La consistència de flotació assoleix valors òptims depenent de les condicions de flux dins la cel·la de flotació. Les metodologies d'anàlisi emprades permeten obtenir diferents factors de destintabilitat. Entre aquests factors existeix una correlació important (determinada pels coeficients de correlació de Pearson) que permet assegurar la utilització de la blancor com a paràmetre fonamental en l'anàlisi del destintatge de papers vells, sempre i quan es complementi amb anàlisis d'imatge o bé amb anàlisi de la concentració efectiva de tinta residual. S'aconsegueixen expressions empíriques tipus exponencial que relacionen aquests factors de destintabilitat amb les variables experimentals. L' estudi de les cinètiques de flotació permet calcular les constants cinètiques (kBlancor, kERIC, kSupimp) en funció de les variables experimentals, obtenint un model empíric de flotació que relacionant-lo amb els paràmetres microscòpics que afecten realment l'eliminació de partícules de tinta, deriva en un model fonamental molt més difícil d'interpretar. Mitjançant l'estudi d'aquestes cinètiques separades per classes, també s'aconsegueix determinar que l'eficàcia d'eliminació de partícules de tinta és màxima si el seu diàmetre equivalent és superior a 50 μm. Les partícules amb diàmetres equivalents inferiors a 15 μm no són eliminades en les condicions de flotació analitzades. Es pot dir que és físicament impossible eliminar partícules de tinta de diàmetres molt diferents amb la mateixa eficiència i sota les mateixes condicions de treball. El rendiment del procés analitzat en funció de l'eliminació de sòlids per l'etapa de flotació no ha presentat relacions significatives amb cap de les variables experimentals analitzades. Únicament es pot concloure que addicionar quantitats elevades de silicat sòdic provoca una disminució tant de sòlids com de matèria inorgànica presents en les escumes de flotació.

Convectively coupled equatorial waves: A new methodology for identifying wave structures in observational data

Convectively coupled equatorial waves are fundamental components of the interaction between the physics and dynamics of the tropical atmosphere. A new methodology, which isolates individual equatorial wave modes, has been developed and applied to observational data. The methodology assumes that the horizontal structures given by equatorial wave theory can be used to project upper- and lower-tropospheric data onto equatorial wave modes. The dynamical fields are first separated into eastward- and westward-moving components with a specified domain of frequency–zonal wavenumber. Each of the components for each field is then projected onto the different equatorial modes using the y structures of these modes given by the theory. The latitudinal scale yo of the modes is predetermined by data to fit the equatorial trapping in a suitable latitude belt y = ±Y. The extent to which the different dynamical fields are consistent with one another in their depiction of each equatorial wave structure determines the confidence in the reality of that structure. Comparison of the analyzed modes with the eastward- and westward-moving components in the convection field enables the identification of the dynamical structure and nature of convectively coupled equatorial waves. In a case study, the methodology is applied to two independent data sources, ECMWF Reanalysis and satellite-observed window brightness temperature (Tb) data for the summer of 1992. Various convectively coupled equatorial Kelvin, mixed Rossby–gravity, and Rossby waves have been detected. The results indicate a robust consistency between the two independent data sources. Different vertical structures for different wave modes and a significant Doppler shifting effect of the background zonal winds on wave structures are found and discussed. It is found that in addition to low-level convergence, anomalous fluxes induced by strong equatorial zonal winds associated with equatorial waves are important for inducing equatorial convection. There is evidence that equatorial convection associated with Rossby waves leads to a change in structure involving a horizontal structure similar to that of a Kelvin wave moving westward with it. The vertical structure may also be radically changed. The analysis method should make a very powerful diagnostic tool for investigating convectively coupled equatorial waves and the interaction of equatorial dynamics and physics in the real atmosphere. The results from application of the analysis method for a reanalysis dataset should provide a benchmark against which model studies can be compared.

Temperature trends derived from Stratospheric Sounding Unit radiances: the effect of increasing CO2 on the weighting function

Stratospheric Sounding Units (SSU) on the NOAA operational satellites have been the main source of near global temperature trend data above the lower stratosphere. They have been used extensively for comparison with model-derived trends. The SSU senses in the 15 micron band of CO2 and hence the weighting function is sensitive to changes in CO2 concentrations. The impact of this change in weighting function has been ignored in all recent trend analyses. We show that the apparent trends in global mean brightness temperature due to the change in weighting function vary from about -0.4 K/decade to 0.4 K/decade depending on the altitude sensed by the different SSU channels. For some channels, this apparent trend is of a similar size to the trend deduced from SSU data but ignoring the change in weighting function. In the mid-stratosphere, the revised trends are now significantly more negative and in better agreement with model-calculated trends.

Effects of sub-pixel heterogeneity on the retrieval of soil moisture from passive microwave radiometry

A simple formulation relating the L-band microwave brightness temperature detected by a passive microwave radiometer to the near surface soil moisture was developed using MICRO-SWEAT, a coupled microwave emission model and soil-vegetation-atmosphere-transfer (SVAT) scheme. This simple model provides an ideal tool with which to explore the impact of sub-pixel heterogeneity on the retrieval of soil moisture from microwave brightness temperatures. In the case of a bare soil pixel, the relationship between apparent emissivity and surface soil moisture is approximately linear, with the clay content of the soil influencing just the intercept of this relationship. It is shown that there are no errors in the retrieved soil moisture from a bare soil pixel that is heterogeneous in soil moisture and texture. However, in the case of a vegetated pixel, the slope of the relationship between apparent emissivity and surface soil moisture decreases with increasing vegetation. Therefore for a pixel that is heterogeneous in vegetation and soil moisture, errors can be introduced into the retrieved soil moisture. Generally, under moderate conditions, the retrieved soil moisture is within 3% of the actual soil moisture. Examples illustrating this discussion use data collected during the Southern Great Plains '97 Experiment (SGP97).

Moving towards a more mechanistic approach in the determination of soil heat flux from remote measurements - II. Diurnal shape of soil heat flux

For vegetated surfaces, calculation of soil heat flux, G, with the Exact or Analytical method requires a harmonic analysis of below-canopy soil surface temperature, to obtain the shape of the diurnal course of G. When determining G with remote sensing methods, only composite (vegetation plus soil) radiometric brightness temperature is available. This paper presents a simple equation that relates the sum of the harmonic terms derived for the composite radiometric surface temperature to that of belowcanopy soil surface temperature. The thermal inertia, Gamma(,) for which a simple equation has been presented in a companion paper, paper I, is used to set the magnitude of G. To assess the success of the method proposed in this paper for the estimation of the diurnal shape of G, a comparison was made between 'remote' and in situ calculated values from described field sites. This indicated that the proposed method was suitable for the estimation of the shape of G for a variety of vegetation types and densities. The approach outlined in paper I, to obtain Gamma, was then combined with the estimated harmonic terms to predict estimates of G, which were compared to values predicted by empirical remote methods found in the literature. This indicated that the method proposed in the combination of papers I and II gave reliable estimates of G, which, in comparison to the other methods, resulted in more realistic predictions for vegetated surfaces. This set of equations can also be used for bare and sparsely vegetated soils, making it a universally applicable method. (C) 2007 Elsevier B.V. All rights reserved.

A sensitivity analysis of soil moisture retrieval from the t-w microwave emission model

The potential of the τ-ω model for retrieving the volumetric moisture content of bare and vegetated soil from dual polarisation passive microwave data acquired at single and multiple angles is tested. Measurement error and several additional sources of uncertainty will affect the theoretical retrieval accuracy. These include uncertainty in the soil temperature, the vegetation structure and consequently its microwave singlescattering albedo, and uncertainty in soil microwave emissivity based on its roughness. To test the effects of these uncertainties for simple homogeneous scenes, we attempt to retrieve soil moisture from a number of simulated microwave brightness temperature datasets generated using the τ-ω model. The uncertainties for each influence are estimated and applied to curves generated for typical scenarios, and an inverse model used to retrieve the soil moisture content, vegetation optical depth and soil temperature. The effect of each influence on the theoretical soil moisture retrieval limit is explored, the likelihood of each sensor configuration meeting user requirements is assessed, and the most effective means of improving moisture retrieval indicated.

Owls and rabbits: predation against substandard individuals of an easy prey

The interactions among the multiple factors regulating predator-prey relationships make predation a more complex process than previously thought. The degree to which substandard individuals are captured disproportionately seems to be better a function of the difficulty of prey capture than of the hunting techniques (coursing vs. ambushing predators). That is, when the capture and killing of a prey species is easy, substandard individuals will be predated in proportion to their occurrence in the prey population. In the present study, we made use of eagle owls Bubo bubo and their main prey, the rabbit Oryctolagus cuniculus: (a) the brightness of the white tails of rabbits seems to be correlated with the physical condition of individuals, (b) by using the tails of predated rabbits as an index of individual condition, we found that eagle owls seem to prefer substandard individuals (characterized by duller tails), and (c) by using information from continuous radiotracking of 14 individuals, we suggest that the difficulty of rabbit capture could be low. Although the relative benefits of preying on substandard individuals should considerably decrease when a predator is attacking an easy prey, we hypothesise that the eagle owl preference for substandard individuals could be due to the easy detection of poor individuals by a visual cue, the brightness of the rabbit tail. Several elements allow us to believe that this form of visual communication between a prey and one of its main predators could be more widespread than previously thought. In fact: (a) visual signalling plays a relevant role in intraspecific communication in eagle owls and, consequently, visual signals could also play a role in interspecific interactions, and (b) empirical studies showed that signals may inform the predator that it has been perceived, or that the prey is in a sufficiently healthy state to elude the predator.

Remote sensing of cloud properties using ground-based measurements of zenith radiance

We have conducted the first extensive field test of two new methods to retrieve optical properties for overhead clouds that range from patchy to overcast. The methods use measurements of zenith radiance at 673 and 870 nm wavelengths and require the presence of green vegetation in the surrounding area. The test was conducted at the Atmospheric Radiation Measurement Program Oklahoma site during September–November 2004. These methods work because at 673 nm (red) and 870 nm (near infrared (NIR)), clouds have nearly identical optical properties, while vegetated surfaces reflect quite differently. The first method, dubbed REDvsNIR, retrieves not only cloud optical depth τ but also radiative cloud fraction. Because of the 1-s time resolution of our radiance measurements, we are able for the first time to capture changes in cloud optical properties at the natural timescale of cloud evolution. We compared values of τ retrieved by REDvsNIR to those retrieved from downward shortwave fluxes and from microwave brightness temperatures. The flux method generally underestimates τ relative to the REDvsNIR method. Even for overcast but inhomogeneous clouds, differences between REDvsNIR and the flux method can be as large as 50%. In addition, REDvsNIR agreed to better than 15% with the microwave method for both overcast and broken clouds. The second method, dubbed COUPLED, retrieves τ by combining zenith radiances with fluxes. While extra information from fluxes was expected to improve retrievals, this is not always the case. In general, however, the COUPLED and REDvsNIR methods retrieve τ to within 15% of each other.

The effect of extratropical cyclones on satellite-retrieved aerosol properties over ocean

Extratropical cyclones may have a signicant effect on column aerosol properties over ocean. European Centre for Medium Range Weather Forecasts (ECMWF) derived storm-centric composites of MODerate resolution Imaging Spectroradiometer (MODIS) and Advanced Along-Track Scanning Radiometer (AATSR) aerosol optical depth and aerosol size parameters are produced for the North Atlantic and the South Atlantic oceans. It is found that retrieved aerosol optical depth and aerosol size both increase near the center of the composite extratropical cyclones. Using composites of ECMWF ERA-Interim reanalysis data, it is demonstrated that wind speed is a considerably more likely explanatory variable than relative humidity for the aerosol observations. A comparison of composites for both MODIS and AATSR, which uses a wind speed dependent sea-surface brightness model in the aerosol retrieval, suggests that although surface brightness eects may contribute towards some of the observations, wind speed dependent emission of sea salt also appears to make a signicant contribution to the observed aerosol properties.

Robotic movement elicits automatic imitation

Recent behavioural and neuroimaging studies have found that observation of human movement, but not of robotic movement, gives rise to visuomotor priming. This implies that the 'mirror neuron' or 'action observation–execution matching' system in the premotor and parietal cortices is entirely unresponsive to robotic movement. The present study investigated this hypothesis using an 'automatic imitation' stimulus–response compatibility procedure. Participants were required to perform a prespecified movement (e.g. opening their hand) on presentation of a human or robotic hand in the terminal posture of a compatible movement (opened) or an incompatible movement (closed). Both the human and the robotic stimuli elicited automatic imitation; the prespecified action was initiated faster when it was cued by the compatible movement stimulus than when it was cued by the incompatible movement stimulus. However, even when the human and robotic stimuli were of comparable size, colour and brightness, the human hand had a stronger effect on performance. These results suggest that effector shape is sufficient to allow the action observation–matching system to distinguish human from robotic movement. They also indicate, as one would expect if this system develops through learning, that to varying degrees both human and robotic action can be 'simulated' by the premotor and parietal cortices.

A simple column model to explore anticipated problems in variational assimilation of satellite observations

We investigate a simplified form of variational data assimilation in a fully nonlinear framework with the aim of extracting dynamical development information from a sequence of observations over time. Information on the vertical wind profile, w(z ), and profiles of temperature, T (z , t), and total water content, qt (z , t), as functions of height, z , and time, t, are converted to brightness temperatures at a single horizontal location by defining a two-dimensional (vertical and time) variational assimilation testbed. The profiles of T and qt are updated using a vertical advection scheme. A basic cloud scheme is used to obtain the fractional cloud amount and, when combined with the temperature field, this information is converted into a brightness temperature, using a simple radiative transfer scheme. It is shown that our model exhibits realistic behaviour with regard to the prediction of cloud, but the effects of nonlinearity become non-negligible in the variational data assimilation algorithm. A careful analysis of the application of the data assimilation scheme to this nonlinear problem is presented, the salient difficulties are highlighted, and suggestions for further developments are discussed.

The diurnal cycle in the tropics

A global archive of high-resolution (3-hourly, 0.58 latitude–longitude grid) window (11–12 mm) brightness temperature (Tb) data from multiple satellites is being developed by the European Union Cloud Archive User Service (CLAUS) project. It has been used to construct a climatology of the diurnal cycle in convection, cloudiness, and surface temperature for all regions of the Tropics. An example of the application of the climatology to the evaluation of the climate version of the U.K. Met. Office Unified Model (UM), version HadAM3, is presented. The characteristics of the diurnal cycle described by the CLAUS data agree with previous observational studies, demonstrating the universality of the characteristics of the diurnal cycle for land versus ocean, clear sky versus convective regimes. It is shown that oceanic deep convection tends to reach its maximum in the early morning. Continental convection generally peaks in the evening, although there are interesting regional variations, indicative of the effects of complex land–sea and mountain–valley breezes, as well as the life cycle of mesoscale convective systems. A striking result from the analysis of the CLAUS data has been the extent to which the strong diurnal signal over land is spread out over the adjacent oceans, probably through gravity waves of varying depths. These coherent signals can be seen for several hundred kilometers and in some instances, such as over the Bay of Bengal, can lead to substantial diurnal variations in convection and precipitation. The example of the use of the CLAUS data in the evaluation of the Met. Office UM has demonstrated that the model has considerable difficulty in capturing the observed phase of the diurnal cycle in convection, which suggests some fundamental difficulties in the model’s physical parameterizations. Analysis of the diurnal cycle represents a powerful tool for identifying and correcting model deficiencies.

The effect of nonlinearity on the variational assimilation of satellite observations using a simple column model

Cloud imagery is not currently used in numerical weather prediction (NWP) to extract the type of dynamical information that experienced forecasters have extracted subjectively for many years. For example, rapidly developing mid-latitude cyclones have characteristic signatures in the cloud imagery that are most fully appreciated from a sequence of images rather than from a single image. The Met Office is currently developing a technique to extract dynamical development information from satellite imagery using their full incremental 4D-Var (four-dimensional variational data assimilation) system. We investigate a simplified form of this technique in a fully nonlinear framework. We convert information on the vertical wind field, w(z), and profiles of temperature, T(z, t), and total water content, qt (z, t), as functions of height, z, and time, t, to a single brightness temperature by defining a 2D (vertical and time) variational assimilation testbed. The profiles of w, T and qt are updated using a simple vertical advection scheme. We define a basic cloud scheme to obtain the fractional cloud amount and, when combined with the temperature field, we convert this information into a brightness temperature, having developed a simple radiative transfer scheme. With the exception of some matrix inversion routines, all our code is developed from scratch. Throughout the development process we test all aspects of our 2D assimilation system, and then run identical twin experiments to try and recover information on the vertical velocity, from a sequence of observations of brightness temperature. This thesis contains a comprehensive description of our nonlinear models and assimilation system, and the first experimental results.

Effects of Thomson-scattering geometry on white-light imaging of an interplanetary shock: synthetic observations from forward magnetohydrodynamic modelling

Stereoscopic white-light imaging of a large portion of the inner heliosphere has been used to track interplanetary coronal mass ejections. At large elongations from the Sun, the white-light brightness depends on both the local electron density and the efficiency of the Thomson-scattering process. To quantify the effects of the Thomson-scattering geometry, we study an interplanetary shock using forward magnetohydrodynamic simulation and synthetic white-light imaging. Identifiable as an inclined streak of enhanced brightness in a time–elongation map, the travelling shock can be readily imaged by an observer located within a wide range of longitudes in the ecliptic. Different parts of the shock front contribute to the imaged brightness pattern viewed by observers at different longitudes. Moreover, even for an observer located at a fixed longitude, a different part of the shock front will contribute to the imaged brightness at any given time. The observed brightness within each imaging pixel results from a weighted integral along its corresponding ray-path. It is possible to infer the longitudinal location of the shock from the brightness pattern in an optical sky map, based on the east–west asymmetry in its brightness and degree of polarisation. Therefore, measurement of the interplanetary polarised brightness could significantly reduce the ambiguity in performing three-dimensional reconstruction of local electron density from white-light imaging.

Tracing troposphere-to-stratosphere transport above a mid-latitude deep convective system

Within the project SPURT (trace gas measurements in the tropopause region) a variety of trace gases have been measured in situ in order to investigate the role of dynamical and chemical processes in the extra-tropical tropopause region. In this paper we report on a flight on 10 November 2001 leading from Hohn, Germany (52�N) to Faro, Portugal (37�N) through a strongly developed deep stratospheric intrusion. This streamer was associated with a large convective system over the western Mediterranean with potentially significant troposphere-to-stratosphere transport. Along major parts of the flight we measured unexpectedly high NOy mixing ratios. Also H2O mixing ratios were significantly higher than stratospheric background levels confirming the extraordinary chemical signature of the probed air masses in the interior of the streamer. Backward trajectories encompassing the streamer enable to analyze the origin and physical characteristics of the air masses and to trace troposphere-to-stratosphere transport. Near the western flank of the streamer features caused by long range transport, such as tropospheric filaments characterized by sudden drops in the O3 and NOy mixing ratios and enhanced CO and H2O can be reconstructed in great detail using the reverse domain filling technique. These filaments indicate a high potential for subsequent mixing with the stratospheric air. At the south-western edge of the streamer a strong gradient in the NOy and the O3 mixing ratios coincides very well with a sharp gradient in potential vorticity in the ECMWF fields. In contrast, in the interior of the streamer the observed highly elevated NOy and H2O mixing ratios up to a potential temperature level of 365K and potential vorticity values of maximum 10 PVU cannot be explained in terms of resolved troposphere-to-stratosphere transport along the backward trajectories. Also mesoscale simulations with a High Resolution Model reveal no direct evidence for convective H2O injection up to this level. Elevated H2O mixing ratios in the ECMWF and HRM are seen only up to about tropopause height at 340 hPa and 270 hPa, respectively, well below flight altitude of about 200 hPa. However, forward tracing of the convective influence as identified by satellite brightness temperature measurements and counts of lightning strokes shows that during this part of the flight the aircraft was closely following the border of an air mass which was heavily impacted by convective activity over Spain and Algeria. This is evidence that deep convection at mid-latitudes may have a large impact on the tracer distribution of the lowermost stratosphere reaching well above the thunderstorms anvils as claimed by recent studies using cloud-resolving models.