Extended optimal estimation techniques for sea surface temperature from the Spinning Enhanced Visible and Infra-Red Imager (SEVIRI)

Sea surface temperature (SST) can be estimated from day and night observations of the Spinning Enhanced Visible and Infra-Red Imager (SEVIRI) by optimal estimation (OE). We show that exploiting the 8.7 μm channel, in addition to the “traditional” wavelengths of 10.8 and 12.0 μm, improves OE SST retrieval statistics in validation. However, the main benefit is an improvement in the sensitivity of the SST estimate to variability in true SST. In a fair, single-pixel comparison, the 3-channel OE gives better results than the SST estimation technique presently operational within the Ocean and Sea Ice Satellite Application Facility. This operational technique is to use SST retrieval coefficients, followed by a bias-correction step informed by radiative transfer simulation. However, the operational technique has an additional “atmospheric correction smoothing”, which improves its noise performance, and hitherto had no analogue within the OE framework. Here, we propose an analogue to atmospheric correction smoothing, based on the expectation that atmospheric total column water vapour has a longer spatial correlation length scale than SST features. The approach extends the observations input to the OE to include the averaged brightness temperatures (BTs) of nearby clear-sky pixels, in addition to the BTs of the pixel for which SST is being retrieved. The retrieved quantities are then the single-pixel SST and the clear-sky total column water vapour averaged over the vicinity of the pixel. This reduces the noise in the retrieved SST significantly. The robust standard deviation of the new OE SST compared to matched drifting buoys becomes 0.39 K for all data. The smoothed OE gives SST sensitivity of 98% on average. This means that diurnal temperature variability and ocean frontal gradients are more faithfully estimated, and that the influence of the prior SST used is minimal (2%). This benefit is not available using traditional atmospheric correction smoothing.

Fireball models for flares in AE Aquarii

We examine the flaring behaviour of the cataclysmic variable AE Aqr in the context of the `magnetic propeller' model for this system. The flares are thought to arise from collisions between high-density regions in the material expelled from the system after interaction with the rapidly rotating magnetosphere of the white dwarf. We calculate the first quantitative models for the flaring and calculate the time-dependent emergent optical spectra from the resulting hot, expanding ball of gas. We compare the results under different assumptions to observations and derive values for the mass, length-scale and temperature of the material involved in the flare. We see that the fits suggest that the secondary star in this system has Population II composition.

On the Nusselt number for frazil ice growth—a correction to ”Frazil evolution in channels“ by Lars Hammar and Hung-Tao Shen

The growth (melt) rate of frazil ice is governed by heat transfer away from (towards) the ice crystal, which can be represented by the Nusselt number. We discuss choices for the Nusselt number and turbulent length scale appropriate for frazil ice and note an inaccuracy in the study ”Frazil evolution in channels“ by Lars Hammar and Hung-Tao Shen, which has also led to potentially significant errors in several other papers. We correct this error and suggest an appropriate strategy for determining the Nusselt number applicable to frazil ice growth and melting.

Influence of mass diffusion in sea ice dynamical models

The mixing of floes of different thickness caused by repeated deformation of the ice cover is modeled as diffusion, and the mass balance equation for sea ice accounting for mass diffusion is developed. The effect of deformational diffusion on the ice thickness balance is shown to reach 1% of the divergence effect, which describes ridging and lead formation. This means that with the same accuracy the mass balance equation can be written in terms of mean velocity rather than mean mass-weighted velocity, which one should correctly use for a multicomponent fluid such as sea ice with components identified by floe thickness. Mixing (diffusion) of sea ice also occurs because of turbulent variations in wind and ocean drags that are unresolved in models. Estimates of the importance of turbulent mass diffusion on the dynamic redistribution of ice thickness are determined using empirical data for the turbulent diffusivity. For long-time-scale prediction (≫5 days), where unresolved atmospheric motion may have a length scale on the order of the Arctic basin and the time scale is larger than the synoptic time scale of atmospheric events, turbulent mass diffusion can exceed 10% of the divergence effect. However, for short-time-scale prediction, for example, 5 days, the unresolved scales are on the order of 100 km, and turbulent diffusion is about 0.1% of the divergence effect. Because inertial effects are small in the dynamics of the sea ice pack, diffusive momentum transfer can be disregarded.

Estimating correlated observation error statistics using an ensemble transform Kalman filter

For certain observing types, such as those that are remotely sensed, the observation errors are correlated and these correlations are state- and time-dependent. In this work, we develop a method for diagnosing and incorporating spatially correlated and time-dependent observation error in an ensemble data assimilation system. The method combines an ensemble transform Kalman filter with a method that uses statistical averages of background and analysis innovations to provide an estimate of the observation error covariance matrix. To evaluate the performance of the method, we perform identical twin experiments using the Lorenz ’96 and Kuramoto-Sivashinsky models. Using our approach, a good approximation to the true observation error covariance can be recovered in cases where the initial estimate of the error covariance is incorrect. Spatial observation error covariances where the length scale of the true covariance changes slowly in time can also be captured. We find that using the estimated correlated observation error in the assimilation improves the analysis.

Structure of the planetary boundary layer over southeast England: modeling and measurements

The Weather Research and Forecasting model was applied to analyze variations in the planetary boundary layer (PBL) structure over Southeast England including central and suburban London. The parameterizations and predictive skills of two nonlocal mixing PBL schemes, YSU and ACM2, and two local mixing PBL schemes, MYJ and MYNN2, were evaluated over a variety of stability conditions, with model predictions at a 3 km grid spacing. The PBL height predictions, which are critical for scaling turbulence and diffusion in meteorological and air quality models, show significant intra-scheme variance (> 20%), and the reasons are presented. ACM2 diagnoses the PBL height thermodynamically using the bulk Richardson number method, which leads to a good agreement with the lidar data for both unstable and stable conditions. The modeled vertical profiles in the PBL, such as wind speed, turbulent kinetic energy (TKE), and heat flux, exhibit large spreads across the PBL schemes. The TKE predicted by MYJ were found to be too small and show much less diurnal variation as compared with observations over London. MYNN2 produces better TKE predictions at low levels than MYJ, but its turbulent length scale increases with height in the upper part of the strongly convective PBL, where it should decrease. The local PBL schemes considerably underestimate the entrainment heat fluxes for convective cases. The nonlocal PBL schemes exhibit stronger mixing in the mean wind fields under convective conditions than the local PBL schemes and agree better with large-eddy simulation (LES) studies.

Improving a global model from the boundary layer: total turbulent energy and the neutral limit Prandtl number

Model intercomparisons have identified important deficits in the representation of the stable boundary layer by turbulence parametrizations used in current weather and climate models. However, detrimental impacts of more realistic schemes on the large-scale flow have hindered progress in this area. Here we implement a total turbulent energy scheme into the climate model ECHAM6. The total turbulent energy scheme considers the effects of Earth’s rotation and static stability on the turbulence length scale. In contrast to the previously used turbulence scheme, the TTE scheme also implicitly represents entrainment flux in a dry convective boundary layer. Reducing the previously exaggerated surface drag in stable boundary layers indeed causes an increase in southern hemispheric zonal winds and large-scale pressure gradients beyond observed values. These biases can be largely removed by increasing the parametrized orographic drag. Reducing the neutral limit turbulent Prandtl number warms and moistens low-latitude boundary layers and acts to reduce longstanding radiation biases in the stratocumulus regions, the Southern Ocean and the equatorial cold tongue that are common to many climate models.

Theoretical insight into diagnosing observation error correlations using observation-minus-background and observation-minus-analysis statistics

To improve the quantity and impact of observations used in data assimilation it is necessary to take into account the full, potentially correlated, observation error statistics. A number of methods for estimating correlated observation errors exist, but a popular method is a diagnostic that makes use of statistical averages of observation-minus-background and observation-minus-analysis residuals. The accuracy of the results it yields is unknown as the diagnostic is sensitive to the difference between the exact background and exact observation error covariances and those that are chosen for use within the assimilation. It has often been stated in the literature that the results using this diagnostic are only valid when the background and observation error correlation length scales are well separated. Here we develop new theory relating to the diagnostic. For observations on a 1D periodic domain we are able to the show the effect of changes in the assumed error statistics used in the assimilation on the estimated observation error covariance matrix. We also provide bounds for the estimated observation error variance and eigenvalues of the estimated observation error correlation matrix. We demonstrate that it is still possible to obtain useful results from the diagnostic when the background and observation error length scales are similar. In general, our results suggest that when correlated observation errors are treated as uncorrelated in the assimilation, the diagnostic will underestimate the correlation length scale. We support our theoretical results with simple illustrative examples. These results have potential use for interpreting the derived covariances estimated using an operational system.

Control of nanomaterial self-assembly in ultrasonically levitated droplets

We demonstrate that acoustic trapping can be used to levitate and manipulate droplets of soft matter, in particular, lyotropic mesophases formed from selfassembly of different surfactants and lipids, which can be analyzed in a contact-less manner by X-ray scattering in a controlled gas-phase environment. On the macroscopic length scale, the dimensions and the orientation of the particle are shaped by the ultrasonic field, while on the microscopic length scale the nanostructure can be controlled by varying the humidity of the atmosphere around the droplet. We demonstrate levitation and in situ phase transitions of micellar, hexagonal, bicontinuous cubic, and lamellar phases. The technique opens up a wide range of new experimental approaches of fundamental importance for environmental, biological, and chemical research.

Simulação computacional da interação fluido-estrutura em bombas de cavidades progressivas

The pumping through progressing cavities system has been more and more employed in the petroleum industry. This occurs because of its capacity of elevation of highly viscous oils or fluids with great concentration of sand or other solid particles. A Progressing Cavity Pump (PCP) consists, basically, of a rotor - a metallic device similar to an eccentric screw, and a stator - a steel tube internally covered by a double helix, which may be rigid or deformable/elastomeric. In general, it is submitted to a combination of well pressure with the pressure generated by the pumping process itself. In elastomeric PCPs, this combined effort compresses the stator and generates, or enlarges, the clearance existing between the rotor and the stator, thus reducing the closing effect between their cavities. Such opening of the sealing region produces what is known as fluid slip or slippage, reducing the efficiency of the PCP pumping system. Therefore, this research aims to develop a transient three-dimensional computational model that, based on single-lobe PCP kinematics, is able to simulate the fluid-structure interaction that occurs in the interior of metallic and elastomeric PCPs. The main goal is to evaluate the dynamic characteristics of PCP s efficiency based on detailed and instantaneous information of velocity, pressure and deformation fields in their interior. To reach these goals (development and use of the model), it was also necessary the development of a methodology for generation of dynamic, mobile and deformable, computational meshes representing fluid and structural regions of a PCP. This additional intermediary step has been characterized as the biggest challenge for the elaboration and running of the computational model due to the complex kinematic and critical geometry of this type of pump (different helix angles between rotor and stator as well as large length scale aspect ratios). The processes of dynamic generation of meshes and of simultaneous evaluation of the deformations suffered by the elastomer are fulfilled through subroutines written in Fortan 90 language that dynamically interact with the CFX/ANSYS fluid dynamic software. Since a structural elastic linear model is employed to evaluate elastomer deformations, it is not necessary to use any CAE package for structural analysis. However, an initial proposal for dynamic simulation using hyperelastic models through ANSYS software is also presented in this research. Validation of the results produced with the present methodology (mesh generation, flow simulation in metallic PCPs and simulation of fluid-structure interaction in elastomeric PCPs) is obtained through comparison with experimental results reported by the literature. It is expected that the development and application of such a computational model may provide better details of the dynamics of the flow within metallic and elastomeric PCPs, so that better control systems may be implemented in the artificial elevation area by PCP

Teleparallel equivalent of non-Abelian Kaluza-Klein theory

Based on the equivalence between a gauge theory for the translation group and general relativity, a teleparallel version of the non-Abelian Kaluza-Klein theory is constructed. In this theory, only the fiber-space turns out to be higher dimensional, spacetime being kept always four dimensional. The resulting model is a gauge theory that unifies, in the Kaluza-Klein sense, gravitational and gauge fields. In contrast with the ordinary Kaluza-Klein models, this theory defines a natural length scale for the compact submanifold of the fiber space, which is shown to be of the order of the Planck length.

Structural modifications in silica sonogels prepared with additions of poly(vinyl alcohol)

Silica wet gels were prepared from acid sonohydrolysis of tetraethoxysilane (TEOS) and additions of poly(vinyl alcohol) (PVA)-water solution. Aerogels were obtained from supercritical CO(2) extraction. The samples were studied by thermal gravimetric (TG) analysis, small-angle X-ray scattering (SAXS), and nitrogen adsorption. The structure of wet gels can be described as a mass fractal with dimension D equal to 2.0 on the whole length scale experimentally probed by SAXS, from similar to 0.3 to similar to 15 nm. Pure and low-PVA-addition wet gels exhibit an upper cutoff accounting for a finite characteristic length xi of the mass fractal structure. Additions , of PVA increase without modifying D, which was attributed to a steric effect of the polymer in the structure. The pore volume fraction of the aerogels diminishes typically about 11% with respect to that of the wet gels, although nitrogen adsorption could be underestimating some porosity. The pore size distribution of the aerogels is shifted toward the mesopore region with the additions of PVA, in a straight relationship with the increase of xi in the wet gels. The thermal stability of the pore size distribution of the aerogels was studied up to 1000 degrees C.

Structural features of silicas prepared in n-heptane/water/ethanol/sodium dodecylsulfate microemulsions

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

From sol to aerogel: a study of the nanostructural characteristics of TEOS derived sonogels

The nanostructural characteristics of acid-catalyzed sonogels are studied along the aging process at 60 degreesC in saturated conditions and after the CO, supercritical extraction (aerogel). The structural evolution was studied by means of small-angle X-ray scattering (SAXS) and UV-Visible absorption techniques. The sonogel exhibits a mass fractal structure in a length scale between zeta - 1/q(0) similar to 5.3 and a(1) similar to 1/q(m) similar to 0.22 nm, as the length scale probed by SAXS. The apparent mass fractal dimension lightly increases from 2.0 for fresh gel until 2.2 for 14 days aging in wet conditions. The UV absorption also increases with the aging time in wet conditions. Both observations are consistent with the syneresis process accompanying the polycondensation progress during aging in saturated conditions. For long aging times, the wet sonogels show a light transition from a mass to a surface fractal. in a very small interval of the length scale, developing an extremely rough surface with fractal dimension D-S similar to 2.9, the fractal characteristics of the sonogels practically do not change with the alcohol exchange. With the CO2 supercritical extraction (aerogel). The interval in the length scale in which the surface fractal is defined increases, while the surface fractal dimension diminishes to D-S similar to 2.5. The mass fractal characteristics are less apparent in the aerogels. (C) 2001 Elsevier B.V. B.V. All rights reserved.

Mass fractal characteristics of wet sonogels as determined by small-angle x-ray scattering and differential scanning calorimetry

Low density silica sonogels were prepared from acid sonohydrolysis of tetraethoxysilane. Wet gels were studied by small-angle x-ray scattering (SAXS) and differential scanning calorimetry (DSC). The DSC tests were carried out under a heating rate of 2 degrees C/min from -120 degrees C up to 30 degrees C. Aerogels were obtained by CO(2) supercritical extraction and characterized by nitrogen adsorption and SAXS. The DSC thermogram displays two distinct endothermic peaks. The first, a broad peak extending from about -80 degrees C up to practically 0 degrees C, was associated to the melting of ice nanocrystals with a crystal size distribution with pore diameter ranging from 1 or 2 nm up to about 60 nm, as estimated from Thomson's equation. The second, a sharp peak with onset temperature close to 0 degrees C, was attributed to the melting of macroscopic crystals. The DSC incremental nanopore volume distribution is in reasonable agreement with the incremental pore volume distribution of the aerogel as determined from nitrogen adsorption. No macroporosity was detected by nitrogen adsorption, probably because the adsorption method applies stress on the sample during measurement, leading to a underestimation of pore volume, or because often positive curvature of the solid surface is in aerogels, making the nitrogen condensation more difficult. According to the SAXS results, the solid network of the wet gels behaves as a mass fractal structure with mass fractal dimension D=2.20 +/- 0.01 in a characteristic length scale below xi=7.9 +/- 0.1 nm. The mass fractal characteristics of the wet gels have also been probed from DSC data by means of an earlier applied modeling for generation of a mass fractal from the incremental pore volume distribution curves. The results are shown to be in interesting agreement with the results from SAXS.