Medem Channel profiles based on stand-alone model simulation and corrections with nudging

State-of-the-art process-based models have shown to be applicable to the simulation and prediction of coastal morphodynamics. On annual to decadal temporal scales, these models may show limitations in reproducing complex natural morphological evolution patterns, such as the movement of bars and tidal channels, e.g. the observed decadal migration of the Medem Channel in the Elbe Estuary, German Bight. Here a morphodynamic model is shown to simulate the hydrodynamics and sediment budgets of the domain to some extent, but fails to adequately reproduce the pronounced channel migration, due to the insufficient implementation of bank erosion processes. In order to allow for long-term simulations of the domain, a nudging method has been introduced to update the model-predicted bathymetries with observations. The model-predicted bathymetry is nudged towards true states in annual time steps. Sensitivity analysis of a user-defined correlation length scale, for the definition of the background error covariance matrix during the nudging procedure, suggests that the optimal error correlation length is similar to the grid cell size, here 80-90 m. Additionally, spatially heterogeneous correlation lengths produce more realistic channel depths than do spatially homogeneous correlation lengths. Consecutive application of the nudging method compensates for the (stand-alone) model prediction errors and corrects the channel migration pattern, with a Brier skill score of 0.78. The proposed nudging method in this study serves as an analytical approach to update model predictions towards a predefined 'true' state for the spatiotemporal interpolation of incomplete morphological data in long-term simulations.

Planktonic foraminifera in Core ML2_66, North Atlantic

Evolution of approaches and methods for reconstruction of paleoenvironmental conditions from microfossils contained in bottom sediments is assessed. Authors elaborated a new actualistic basis for such reconstructions, consisting of a database on contents of tests of planktonic foraminifers in the surface layer of Atlantic sediments and a package of mathematical tools for computer data processing. Structure of the database is described. It contains data on test contents for 29 species and varieties of planktonic foraminifers in 381 samples. A mathematical model designed for reconstructions is based on factor analysis and multidimensional spline interpolation. The model allows one to deduce Quaternary hydrological parameters (paleotemperature, paleosalinity) for standard hydrological levels down to depth of 250 m for the four seasons of the year. Reconstructions are illustrated by an example of a sedimentary core from the North Atlantic representing a period of 300 ky. During the next to last and the last maxima of continental glaciation (oxygen isotope stages 8, 6, 4, and 2), the subarctic water mass was spread here. Winter and summer surface water temperatures comprised 1-5° and 5-7°C, respectively. During interglacials and in Holocene the conditions were close to present ones: winter and summer surface water temperatures comprised 10-12 and 15-17°C, respectively. Vertical paleohydrological profiles compiled for peaks of climatostratigraphic intervals suggest that during cold intervals water stratification was stronger than during the warm ones. At depth 50 m seasonal salinity oscillations did not exceed 0.4 per mil and commonly salinity was minimum in winter and maximum in summer.

Biogenic opal and carbonate concentration in late Neogene sediments of ODP Site 114-704

A preliminary composite depth section was generated for Site 704 by splicing Holes 704A and 704B together over the interval 0-350 mbsf (0-9 m.y.). High-resolution carbonate and opal data from the cores were correlated with the calcium and silicon signals from the GST logging run in Hole 704B to identify missing and disturbed intervals in the cores. Paleomagnetic and biostratigraphic age boundaries were then transferred to the composite depth records to obtain an age model, and sedimentation rates were calculated by linear interpolation between datums. Algorithms relating measured dry-bulk density to carbonate content and depth were generated to produce predicted values of density for every sample. Accumulation rates of bulk, carbonate, opal, and terrigenous sediment components were then computed to generate a record of sediment deposition on the Meteor Rise that has a resolution of better than 200,000 yr for the period from 8.6 to 1.0 m.y. From 8.6 to 2.5 m.y., bulk-accumulation rates on the Meteor Rise averaged less than 2 g/cm**2/1000 yr and were dominated by carbonate deposition. The first significant opal deposition (6.0 m.y.) punctuated a brief (less than 0.6 Ma) approach of the Polar Front Zone (PFZ) northward that heralded a period of increasing severity of periodic carbonate dissolution events (terrigenous maxima) that abruptly terminated at 4.8 m.y. (base of the Thvera Subchron), synchronous with the reflooding of the Mediterranean after the Messinian salinity crisis. From 4.8 to 2.5 m.y., carbonate again dominated deposition, and the PFZ was far south except during brief northward excursions bracketing 4.2-3.9, 3.3-2.9, and 2.8-2.7 m.y. At 2.5 m.y., all components of bulk-accumulation rates increased dramatically (up to 15 g/cm2/1000 yr), and by 2.4 m.y., a pattern of alternating, high-amplitude carbonate and opal cyclicity marked the initiation of rapid glacial to interglaci·l swings in the position of the PFZ, synchronous with the "onset" of major Northern Hemisphere glaciation. Both mass-accumulation rates and the amplitude of the cycles decreased by about 2 m.y., but opal accumulation rates remained high up through the base of the Jaramillo (0.98 m.y.). From 1.9 to 1 m.y., the record is characterized by moderate amplitude fluctuations in carbonate and opal. This record of opal accumulation rates is interpreted as a long-term "Polar Front Indicator" that monitors the advance and retreat of the opal-rich PFZ northward (southward) toward (away from) the Meteor Rise in the subantarctic sector of the South Atlantic Ocean. The timing of PFZ migrations in the subantarctic South Atlantic Ocean is remarkably similar to Pliocene-Pleistocene climate records deduced from benthic oxygen isotope records in the North Atlantic Ocean (Raymo et al., 1989, doi:10.1029/PA004i004p00413; Ruddiman et al., 1989, doi:10.1029/PA004i004p00353). These include northward migrations during "cold" intervals containing strong glacial isotope stages (2.4-2.3, 2.1-2.0, 1.95-1.55, 1.45-1.30 m.y. and at about 1.13 and 1.09 m.y.) and southward migrations during "warm" intervals containing weak glacial and/or strong interglacial stages (2.45-2.40, 2.30-2.10, 2.00-1.95, 1.52-1.45, 1.30-1.18, 1.11, and 1.06-0.93 m.y.). Although our preliminary composite record is not continuous (some stages are obviously missing), there is hope that future work will identify these missing intervals in the as yet incomplete Hole 704B and will extend this high-resolution Southern Hemisphere climate record back to 8.6 m.y.

(Table 2) Age controls used to derive a time scale for DSDP Site 73-522

The stratotype section for the base of the Miocene is at a reversed (below) to normal (above) magnetic transition that is claimed to represent magnetic chron C6Cn.2n (o). Deep Sea Drilling Project (DSDP) Site 522 is the only location we are aware of that unambiguously records the three normal events of C6Cn. We have quantitatively determined the range of the short-lived nannofossil Sphenolithus delphix and the lower limit of S. disbelemnos in DSDP Holes 522 and 522A in order to calibrate their precise relationship to the magnetostratigraphy and to confirm the completeness of the record at this site. Astronomical tuning of Ocean Drilling Program (ODP) Sites 926, 928, and 929 shows that S. disbelemnos appears at 22.67 Ma and that the entire range of S. delphix is from about 22.98 Ma to 23.24 Ma. Using these ages, linear interpolation in DSDP Site 522 suggests that the age of C6Cn.2n (o) and of the Oligocene-Miocene boundary is 22.92+/-0.04 Ma. Our value, conservatively expressed as 22.9+/-0.1 Ma, is 0.9 m.y. younger than the currently accepted age of the Oligocene-Miocene boundary and of C6Cn.2n (o), which was assigned an age of 23.8 Ma, based on an estimate of 23.8+/-1 Ma for the Oligocene-Miocene boundary. The bulk-sediment carbon isotope data from DSDP Site 522 is correlated to the record from benthic foraminifera at ODP Site 929 to refine the calibration of magnetic reversals from C6Cn.1n (o) to C7n.2n (o) at DSDP Site 522 on the astronomical time scale.

Multimodal Imaging of Cotton Wool Spots in Branch Retinal Vein Occlusion.

PURPOSE: To describe and follow cotton wool spots (CWS) in branch retinal vein occlusion (BRVO) using multimodal imaging. METHODS: In this prospective cohort study including 24 patients with new-onset BRVO, CWS were described and analyzed in color fundus photography (CF), spectral domain optical coherence tomography (SD-OCT), infrared (IR) and fluorescein angiography (FA) every 3 months for 3 years. The CWS area on SD-OCT and CF was evaluated using OCT-Tool-Kit software: CWS were marked in each single OCT B-scan and the software calculated the area by interpolation. RESULTS: 29 central CWS lesions were found. 100% of these CWS were visible on SD-OCT, 100% on FA and 86.2% on IR imaging, but only 65.5% on CF imaging. CWS were visible for 12.4 ± 7.5 months on SD-OCT, for 4.4 ± 3 months and 4.3 ± 3.4 months on CF and on IR, respectively, and for 17.5 ± 7.1 months on FA. The evaluated CWS area on SD-OCT was larger than on CF (0.26 ± 0.17 mm(2) vs. 0.13 ± 0.1 mm(2), p < 0.0001). The CWS area on SD-OCT and surrounding pathology such as intraretinal cysts, avascular zones and intraretinal hemorrhage were predictive for how long CWS remained visible (r(2) = 0.497, p < 0.002). CONCLUSIONS: The lifetime and presentation of CWS in BRVO seem comparable to other diseases. SD-OCT shows a higher sensitivity for detecting CWS compared to CF. The duration of visibility of CWS varies among different image modalities and depends on the surrounding pathology and the CWS size.

Aplicación de Métodos Meshless al Análisis de Problemas de Autovalores

La presente Tesis Doctoral aborda la aplicación de métodos meshless, o métodos sin malla, a problemas de autovalores, fundamentalmente vibraciones libres y pandeo. En particular, el estudio se centra en aspectos tales como los procedimientos para la resolución numérica del problema de autovalores con estos métodos, el coste computacional y la viabilidad de la utilización de matrices de masa o matrices de rigidez geométrica no consistentes. Además, se acomete en detalle el análisis del error, con el objetivo de determinar sus principales fuentes y obtener claves que permitan la aceleración de la convergencia. Aunque en la actualidad existe una amplia variedad de métodos meshless en apariencia independientes entre sí, se han analizado las diferentes relaciones entre ellos, deduciéndose que el método Element-Free Galerkin Method [Método Galerkin Sin Elementos] (EFGM) es representativo de un amplio grupo de los mismos. Por ello se ha empleado como referencia en este análisis. Muchas de las fuentes de error de un método sin malla provienen de su algoritmo de interpolación o aproximación. En el caso del EFGM ese algoritmo es conocido como Moving Least Squares [Mínimos Cuadrados Móviles] (MLS), caso particular del Generalized Moving Least Squares [Mínimos Cuadrados Móviles Generalizados] (GMLS). La formulación de estos algoritmos indica que la precisión de los mismos se basa en los siguientes factores: orden de la base polinómica p(x), características de la función de peso w(x) y forma y tamaño del soporte de definición de esa función. Se ha analizado la contribución individual de cada factor mediante su reducción a un único parámetro cuantificable, así como las interacciones entre ellos tanto en distribuciones regulares de nodos como en irregulares. El estudio se extiende a una serie de problemas estructurales uni y bidimensionales de referencia, y tiene en cuenta el error no sólo en el cálculo de autovalores (frecuencias propias o carga de pandeo, según el caso), sino también en términos de autovectores. This Doctoral Thesis deals with the application of meshless methods to eigenvalue problems, particularly free vibrations and buckling. The analysis is focused on aspects such as the numerical solving of the problem, computational cost and the feasibility of the use of non-consistent mass or geometric stiffness matrices. Furthermore, the analysis of the error is also considered, with the aim of identifying its main sources and obtaining the key factors that enable a faster convergence of a given problem. Although currently a wide variety of apparently independent meshless methods can be found in the literature, the relationships among them have been analyzed. The outcome of this assessment is that all those methods can be grouped in only a limited amount of categories, and that the Element-Free Galerkin Method (EFGM) is representative of the most important one. Therefore, the EFGM has been selected as a reference for the numerical analyses. Many of the error sources of a meshless method are contributed by its interpolation/approximation algorithm. In the EFGM, such algorithm is known as Moving Least Squares (MLS), a particular case of the Generalized Moving Least Squares (GMLS). The accuracy of the MLS is based on the following factors: order of the polynomial basis p(x), features of the weight function w(x), and shape and size of the support domain of this weight function. The individual contribution of each of these factors, along with the interactions among them, has been studied in both regular and irregular arrangement of nodes, by means of a reduction of each contribution to a one single quantifiable parameter. This assessment is applied to a range of both one- and two-dimensional benchmarking cases, and includes not only the error in terms of eigenvalues (natural frequencies or buckling load), but also of eigenvectors

Analysis of frame-compatible subsampling structures for efficient 3DTV broadcast

The evolution of the television market is led by 3DTV technology, and this tendency can accelerate during the next years according to expert forecasts. However, 3DTV delivery by broadcast networks is not currently developed enough, and acts as a bottleneck for the complete deployment of the technology. Thus, increasing interest is dedicated to ste-reo 3DTV formats compatible with current HDTV video equipment and infrastructure, as they may greatly encourage 3D acceptance. In this paper, different subsampling schemes for HDTV compatible transmission of both progressive and interlaced stereo 3DTV are studied and compared. The frequency characteristics and preserved frequency content of each scheme are analyzed, and a simple interpolation filter is specially designed. Finally, the advantages and disadvantages of the different schemes and filters are evaluated through quality testing on several progressive and interlaced video sequences.

Kernel bandwidth estimation for moving object detection in non-stabilized cameras

The evolution of the television market is led by 3DTV technology, and this tendency can accelerate during the next years according to expert forecasts. However, 3DTV delivery by broadcast networks is not currently developed enough, and acts as a bottleneck for the complete deployment of the technology. Thus, increasing interest is dedicated to ste-reo 3DTV formats compatible with current HDTV video equipment and infrastructure, as they may greatly encourage 3D acceptance. In this paper, different subsampling schemes for HDTV compatible transmission of both progressive and interlaced stereo 3DTV are studied and compared. The frequency characteristics and preserved frequency content of each scheme are analyzed, and a simple interpolation filter is specially designed. Finally, the advantages and disadvantages of the different schemes and filters are evaluated through quality testing on several progressive and interlaced video sequences.

Formation energy of vacancies in FeCr alloys: Dependence on Cr concentration

A modified version of the concentration-dependent model (CDM) potential (A. Caro et al., Phys. Rev. Lett. 95 (2005) 075702) [1] has been developed to study defects in Fe–Cr for different Crconcentrations. A comparison between this new potential and DFT results for a variety of point defect configurations is performed in order to test its reliability for radiation damage studies. The effect of Crconcentration on the vacancyformationenergy in Fe–Cr alloys is analyzed in detail. This study shows a linear dependence of the vacancyformationenergy on Crconcentration for values above 6% of Cr. However, the formationenergy deviates from the linear interpolation in the region below 6% Crconcentration. In order to understand this behavior, the influence of the relative positions between Cr atoms and vacant sites on the vacancyformationenergy has been studied.

Improving runoff estimates from regional climate models: a performance analysis in Spain

An important step to assess water availability is to have monthly time series representative of the current situation. In this context, a simple methodology is presented for application in large-scale studies in regions where a properly calibrated hydrologic model is not available, using the output variables simulated by regional climate models (RCMs) of the European project PRUDENCE under current climate conditions (period 1961–1990). The methodology compares different interpolation methods and alternatives to generate annual times series that minimise the bias with respect to observed values. The objective is to identify the best alternative to obtain bias-corrected, monthly runoff time series from the output of RCM simulations. This study uses information from 338 basins in Spain that cover the entire mainland territory and whose observed values of natural runoff have been estimated by the distributed hydrological model SIMPA. Four interpolation methods for downscaling runoff to the basin scale from 10 RCMs are compared with emphasis on the ability of each method to reproduce the observed behaviour of this variable. The alternatives consider the use of the direct runoff of the RCMs and the mean annual runoff calculated using five functional forms of the aridity index, defined as the ratio between potential evapotranspiration and precipitation. In addition, the comparison with respect to the global runoff reference of the UNH/GRDC dataset is evaluated, as a contrast of the “best estimator” of current runoff on a large scale. Results show that the bias is minimised using the direct original interpolation method and the best alternative for bias correction of the monthly direct runoff time series of RCMs is the UNH/GRDC dataset, although the formula proposed by Schreiber (1904) also gives good results

Theoretical Analysis of the No-Slip Boundary Condition Enforcement in SPH Methods

The aim of the present work is to provide an in-depth analysis of the most representative mirroring techniques used in SPH to enforce boundary conditions (BC) along solid profiles. We specifically refer to dummy particles, ghost particles, and Takeda et al. [Prog. Theor. Phys. 92 (1994), 939] boundary integrals. The analysis has been carried out by studying the convergence of the first- and second-order differential operators as the smoothing length (that is, the characteristic length on which relies the SPH interpolation) decreases. These differential operators are of fundamental importance for the computation of the viscous drag and the viscous/diffusive terms in the momentum and energy equations. It has been proved that close to the boundaries some of the mirroring techniques leads to intrinsic inaccuracies in the convergence of the differential operators. A consistent formulation has been derived starting from Takeda et al. boundary integrals (see the above reference). This original formulation allows implementing no-slip boundary conditions consistently in many practical applications as viscous flows and diffusion problems.