A randomized trial of the effects of the noble gases helium and argon on neuroprotection in a rodent cardiac arrest model.

BACKGROUND The noble gas xenon is considered as a neuroprotective agent, but availability of the gas is limited. Studies on neuroprotection with the abundant noble gases helium and argon demonstrated mixed results, and data regarding neuroprotection after cardiac arrest are scant. We tested the hypothesis that administration of 50% helium or 50% argon for 24 h after resuscitation from cardiac arrest improves clinical and histological outcome in our 8 min rat cardiac arrest model. METHODS Forty animals had cardiac arrest induced with intravenous potassium/esmolol and were randomized to post-resuscitation ventilation with either helium/oxygen, argon/oxygen or air/oxygen for 24 h. Eight additional animals without cardiac arrest served as reference, these animals were not randomized and not included into the statistical analysis. Primary outcome was assessment of neuronal damage in histology of the region I of hippocampus proper (CA1) from those animals surviving until day 5. Secondary outcome was evaluation of neurobehavior by daily testing of a Neurodeficit Score (NDS), the Tape Removal Test (TRT), a simple vertical pole test (VPT) and the Open Field Test (OFT). Because of the non-parametric distribution of the data, the histological assessments were compared with the Kruskal-Wallis test. Treatment effect in repeated measured assessments was estimated with a linear regression with clustered robust standard errors (SE), where normality is less important. RESULTS Twenty-nine out of 40 rats survived until day 5 with significant initial deficits in neurobehavioral, but rapid improvement within all groups randomized to cardiac arrest. There were no statistical significant differences between groups neither in the histological nor in neurobehavioral assessment. CONCLUSIONS The replacement of air with either helium or argon in a 50:50 air/oxygen mixture for 24 h did not improve histological or clinical outcome in rats subjected to 8 min of cardiac arrest.

Respiratory rhythm generation: A modeling study of the respiratory central pattern generator and the phrenic motor neuron

The respiratory central pattern generator is a collection of medullary neurons that generates the rhythm of respiration. The respiratory central pattern generator feeds phrenic motor neurons, which, in turn, drive the main muscle of respiration, the diaphragm. The purpose of this thesis is to understand the neural control of respiration through mathematical models of the respiratory central pattern generator and phrenic motor neurons. ^ We first designed and validated a Hodgkin-Huxley type model that mimics the behavior of phrenic motor neurons under a wide range of electrical and pharmacological perturbations. This model was constrained physiological data from the literature. Next, we designed and validated a model of the respiratory central pattern generator by connecting four Hodgkin-Huxley type models of medullary respiratory neurons in a mutually inhibitory network. This network was in turn driven by a simple model of an endogenously bursting neuron, which acted as the pacemaker for the respiratory central pattern generator. Finally, the respiratory central pattern generator and phrenic motor neuron models were connected and their interactions studied. ^ Our study of the models has provided a number of insights into the behavior of the respiratory central pattern generator and phrenic motor neurons. These include the suggestion of a role for the T-type and N-type calcium channels during single spikes and repetitive firing in phrenic motor neurons, as well as a better understanding of network properties underlying respiratory rhythm generation. We also utilized an existing model of lung mechanics to study the interactions between the respiratory central pattern generator and ventilation. ^

Metastatic lung cancer in a new mouse model inheriting p53 and latent K-ras mutations

Lung cancer is a devastating disease with very poor prognosis. The design of better treatments for patients would be greatly aided by mouse models that closely resemble the human disease. The most common type of human lung cancer is adenocarcinoma with frequent metastasis. Unfortunately, current models for this tumor are inadequate due to the absence of metastasis. Based on the molecular findings in human lung cancer and metastatic potential of osteosarcomas in mutant p53 mouse models, I hypothesized that mice with both K-ras and p53 missense mutations might develop metastatic lung adenocarcinomas. Therefore, I incorporated both K-rasLA1 and p53RI72HΔg alleles into mouse lung cells to establish a more faithful model for human lung adenocarcinoma and for translational and mechanistic studies. Mice with both mutations ( K-rasLA1/+ p53R172HΔg/+) developed advanced lung adenocarcinomas with similar histopathology to human tumors. These lung adenocarcinomas were highly aggressive and metastasized to multiple intrathoracic and extrathoracic sites in a pattern similar to that seen in lung cancer patients. This mouse model also showed gender differences in cancer related death and developed pleural mesotheliomas in 23.2% of them. In a preclinical study, the new drug Erlotinib (Tarceva) decreased the number and size of lung lesions in this model. These data demonstrate that this mouse model most closely mimics human metastatic lung adenocarcinoma and provides an invaluable system for translational studies. ^ To screen for important genes for metastasis, gene expression profiles of primary lung adenocarcinomas and metastases were analyzed. Microarray data showed that these two groups were segregated in gene expression and had 79 highly differentially expressed genes (more than 2.5 fold changes and p<0.001). Microarray data of Bub1b, Vimentin and CCAM1 were validated in tumors by quantitative real-time PCR (QPCR). Bub1b , a mitotic checkpoint gene, was overexpressed in metastases and this correlated with more chromosomal abnormalities in metastatic cells. Vimentin, a marker of epithelial-mesenchymal transition (EMT), was also highly expressed in metastases. Interestingly, Twist, a key EMT inducer, was also highly upregulated in metastases by QPCR, and this significantly correlated with the overexpression of Vimentin in the same tumors. These data suggest EMT occurs in lung adenocarcinomas and is a key mechanism for the development of metastasis in K-ras LA1/+ p53R172HΔg/+ mice. Thus, this mouse model provides a unique system to further probe the molecular basis of metastatic lung cancer.^

Comparative studies of Hox genes in mammalian limb and vertebral pattern formation

Divergence of anterior-posterior (AP) limb pattern and differences in vertebral column morphology are the two main examples of mammalian evolution. The Hox genes (homeobox containing gene) have been implicated in driving evolution of these structures. However, regarding Hox genes, how they contribute to the generation of mammalian morphological diversities, is still unclear. Implementing comparative gene expression and phenotypic rescue studies for different mammalian Hox genes could aid in unraveling this mystery. In the first part of this thesis, the expression pattern of Hoxd13 gene, a key Hox gene in the establishment of the limb AP pattern, was examined in developing limbs of bats and mice. Bat forelimbs exhibit a pronounced asymmetric AP pattern and offer a good model to study the molecular mechanisms that contribute to the variety of mammalian limbs. The data showed that the expression domain of bat Hoxd13 was shifted prior to the asymmetric limb plate expansion, whereas its domain in mice was much more symmetric. This finding reveals a correlation between the divergence of Hoxd13 expression and the AP patterning difference in limb development. The second part of this thesis details a phenotypic rescue approach by human HOXB1-9 transgenes in mice with Hoxb1-9 deletion, The mouse mutants displayed homeosis in cervical and anterior thoracic vertebrae. The human transgenes entirely rescued the mouse mutants, suggesting that these human HOX genes have similar functions to their mouse orthologues in anterior axial skeletal patterning. The anterior expressing human HOXB transgenes such as HOXB1-3 were expressed in the mouse embryonic trunk in a similar manner as their murine orthologues. However, the anterior boundary of human HOXB9 expression domain was more posterior than that of the mouse Hoxb9 by 2-3 somites. These data provide the molecular support for the hypothesis that Hox genes are responsible for maintaining similar anterior axial skeletal architectures cervical and anterior thoracic regions, but different architectures in lumbar and posterior thoracic regions between humans and mice. ^

Body mass index trajectories and predictors among 3rd to 12th graders using growth curve mixture modeling the Child and Adolescent Trial for Cardiovascular Health (CATCH) study

This dissertation examined body mass index (BMI) growth trajectories and the effects of gender, ethnicity, dietary intake, and physical activity (PA) on BMI growth trajectories among 3rd to 12th graders (9-18 years of age). Growth curve model analysis was performed using data from The Child and Adolescent Trial for Cardiovascular Health (CATCH) study. The study population included 2909 students who were followed up from grades 3-12. The main outcome was BMI at grades 3, 4, 5, 8, and 12. ^ The results revealed that BMI growth differed across two distinct developmental periods of childhood and adolescence. Rate of BMI growth was faster in middle childhood (9-11 years old or 3rd - 5th grades) than in adolescence (11-18 years old or 5th - 12th grades). Students with higher BMI at 3rd grade (baseline) had faster rates of BMI growth. Three groups of students with distinct BMI growth trajectories were identified: high, average, and low. ^ Black and Hispanic children were more likely to be in the groups with higher baseline BMI and faster rates of BMI growth over time. The effects of gender or ethnicity on BMI growth differed across the three groups. The effects of ethnicity on BMI growth were weakened as the children aged. The effects of gender on BMI growth were attenuated in the groups with a large proportion of black and Hispanic children, i.e., “high” or “average” BMI trajectory group. After controlling for gender, ethnicity, and age at baseline, in the “high BMI trajectory”, rate of yearly BMI growth in middle childhood increased 0.102 for every 500 Kcals increase (p=0.049). No significant effects of percentage of energy from total fat and saturated fat on BMI growth were found. Baseline BMI increased 0.041 for every 30 minutes increased in moderate-to-vigorous PA (MVPA) in the “low BMI trajectory”, while Baseline BMI decreased 0.345 for every 30 minutes increased in vigorous PA (VPA) in the “high BMI trajectory”. ^ Childhood overweight and obesity interventions should start at the earliest possible ages, prior to 3rd grade and continue through grade school. Interventions should focus on all children, but specifically black and Hispanic children, who are more likely to be highest at-risk. Promoting VPA earlier in childhood is important for preventing overweight and obesity among children and adolescents. Interventions should target total energy intake, rather than only percentage of energy from total fat or saturated fat. ^

Design of long-term animal bioassay for low-dose extrapolation using a multistage model

Conventional designs of animal bioassays allocate the same number of animals into control and dose groups to explore the spontaneous and induced tumor incidence rates, respectively. The purpose of such bioassays are (a) to determine whether or not the substance exhibits carcinogenic properties, and (b) if so, to estimate the human response at relatively low doses. In this study, it has been found that the optimal allocation to the experimental groups which, in some sense, minimize the error of the estimated response for low dose extrapolation is associated with the dose level and tumor risk. The number of dose levels has been investigated at the affordable experimental cost. The pattern of the administered dose, 1 MTD, 1/2 MTD, 1/4 MTD,....., etc. plus control, gives the most reasonable arrangement for the low dose extrapolation purpose. The arrangement of five dose groups may make the highest dose trivial. A four-dose design can circumvent this problem and has also one degree of freedom for testing the goodness-of-fit of the response model.^ An example using the data on liver tumors induced in mice in a lifetime study of feeding dieldrin (Walker et al., 1973) is implemented with the methodology. The results are compared with conclusions drawn from other studies. ^

RESPIRATORY COMPENSATION MECHANISMS IN THE UPPER AND LOWER RESPIRATORY TRACTS OF AN ANIMAL MODEL AFTER SUBCHRONIC INHALATION EXPOSURE TO FORMALDEHYDE: IMPLICATIONS FOR TOXICOLOGY RISK ASSESSMENT (DEPOSITION, DOSE RECEIVED, RESPONSE)

The potential for significant human populations to experience long-term inhalation of formaldehyde and reports of symptomatology due to this exposure has led to a considerable interest in the toxicologic assessment of risk from subchronic formaldehyde exposures using animal models. Since formaldehyde inhalation depresses certain respiratory parameters in addition to its other forms of toxicity, there is a potential for the alteration of the actual dose received by the exposed individual (and the resulting toxicity) due to this respiratory effect. The respiratory responses to formaldehyde inhalation and the subsequent pattern of deposition were therefore investigated in animals that had received subchronic exposure to the compound, and the potential for changes in the formaldehyde dose received due to long-term inhalation evaluated. Male Sprague-Dawley rats were exposed to either 0, 0.5, 3, or 15 ppm formaldehyde for 6 hours/day, 5 days/week for up to 6 months. The patterns of respiratory response, deposition and the compensation mechanisms involved were then determined in a series of formaldehyde test challenges to both the upper and to the lower respiratory tracts in separate groups of subchronically exposed animals and age-specific controls (four concentration groups, two time points). In both the control and pre-exposed animals, there was a characteristic recovery of respiratory parameters initially depressed by formaldehyde inhalation to at or approaching pre-exposure levels within 10 minutes of the initiation of exposure. Also, formaldehyde deposition was found to remain very high in the upper and lower tracts after long-term exposure. Therefore, there was probably little subsequent effect on the dose received by the exposed individual that was attributable to the repeated exposures. There was a diminished initial minute volume response in test challenges of both the upper and lower tracts of animals that had received at least 16 weeks of exposure to 15 ppm, with compensatory increases in tidal volume in the upper tract and respiratory rate in the lower tract. However, this dose-related effect was probably not relevant to human risk estimation because this formaldehyde dose is in excess of that experienced by human populations. ^

Ocean temperature response to idealized Gleissberg and de Vries solar cycles in a comprehensive climate model

The ~90-year Gleissberg and ~200-year de Vries cycles have been identified as two distinctive quasi-periodic components of Holocene solar activity. Evidence exists for the impact of such multi-decadal to centennial-scale variability in total solar irradiance (TSI) on climate, but concerning the ocean, this evidence is mainly restricted to the surface response. Here we use a comprehensive global climate model to study the impact of idealized solar forcing, representing the Gleissberg and de Vries cycles, on global ocean potential temperature at different depth levels, after a recent proxy record indicates a signal of TSI anomalies in the northeastern Atlantic at mid-depth. Potential impacts of TSI anomalies on deeper oceanic levels are climatically relevant due to their possible effect on ocean circulation by altering water mass characteristics. Simulated solar anomalies are shown to penetrate the ocean down to at least deep-water levels. Despite the fact that the two forcing periods differ only by a factor of ~2, the spatial pattern of response is significantly distinctive between the experiments, suggesting different mechanisms for solar signal propagation. These are related to advection by North Atlantic Deep Water flow (200-year forcing), and barotropic adjustment in the South Atlantic in response to a latitudinal shift of the westerly wind belt (90-year forcing).

Contents and accumulation rates of the rare earth and some other elements in metalliferous sediments from DSDP Hole 92-598

The sediments recovered during DSDP Leg 92 (Site 598) include a complete 16 m.y. record of hydrothermal sedimentation along the western flank of the East Pacific Rise at 19°S. Fifty samples from this sediment column were analyzed to test the hypothesis that the REE composition of the hydrothermal component is primarily acquired via scavenging from seawater. Site 598 provides an ideal sample suite for this purpose: the sediments are lithologically "simple," primarily consisting of a mixture of hydrothermal materials and biogenous carbonates; the composition of the hydrothermal component is essentially constant through space and time; and the sediments have undergone minimal diagenetic alteration. The following observations suggest the above-stated hypothesis is true. The Ce anomaly as well as key indices of light and heavy REE behavior all show that the REE pattern of hydrothermal sediments approaches that of seawater with increasing paleodistance from the rise crest. Moreover, shale-normalized REE patterns are similar to that of seawater, varying only in absolute REE content: the REE content increases with distance from the paleo-rise crest and exhibits a pronounced increase in sediments deposited below the paleolysocline. Based on significant correlative relationships between paleodistance from the rise crest and both the concentration and mass accumulation rates (MARs) of REEs and Fe, we conclude the REEs in the hydrothermal component are derived from the interaction of seawater and Fe in the hydrothermal plume.

Concentration of the rare-earth elements in metalliferous sediments from DSDP Leg 92 holes

DSDP Leg 92 drilled at four sites along an east-west transect at 19°S on the western flank of the East Pacific Rise (EPR), in an area where sediments are essentially a mixture of hydrothermal and biogenic components, with only a minimal contribution of clastic material. Rare-earth element (REE) data on the metalliferous (non-carbonate) fraction of samples ranging in age from ~2 to ~27 Ma indicate the existence of two distinct groups of patterns corresponding to two broad age groups, one <=8 Ma, the other >=10 Ma. Within each group, REE patterns have characteristics which are near-uniform, despite large variations in total REE abundances. Sediments of the younger group are enriched in light REE (LREE) relative to deep bottom waters influenced by the hydrothermal plume extending west from the EPR at 19°S. Sediments of the older groups show further relative LREE enrichment and/or heavy REE (HREE) depletion. Surficial sediments deposited beneath the lysocline have high Sum REE concentrations resulting from slow accumulation rates, and patterns resembling older sediments due to early diagenetic effects. A correlation between the mass accumulation rates (MAR) of Sum REE and Fe + Mn suggests that ferromanganese particulate matter supplied by the hydrothermal plume scavenges REE; during this process the LREE are preferentially removed from plume seawater. The MAR of Fe + Mn shows a general decrease with age above basement, whereas Sum REE concentrations in the metalliferous component increase with age above basement. This supports the Ruhlin and Owen model wherein limited scavenging of REE, due to rapid burial of sediment near the palaeo-axis, leads to low concentrations (but high MAR-values) for the REE. Following deposition and burial of the hydrothermal component, further relative flattening of the REE pattern takes place, probably the result of diagenetic reactions over several million years. Phase partitioning data indicate that the proportion of REE residing in more poorly crystalline phases tends to increase with age (from ~45% to 90% of Sum REE). This suggests that as initial ferromanganese precipitates undergo diagenetic recrystallization, REE are transferred to the poorly crystalline phases, and/or are scavenged from pore waters by these phases. Because of the various modifications to REE patterns apparently produced both in the water column and post-depositional settings, the REE patterns of metalliferous sediments will not reflect fine-scale REE variations in associated oceanic water masses.

Age model for ODP Leg 175 sites

Upwelling along the western coast of Africa south of the equator may be partitioned into three major areas, each having its own dynamics and history: (1) the eastern equatorial region, comprising the Congo Fan and the area of Mid-Angola; (2) the Namibia upwelling system, extending from the Walvis Ridge to Lüderitz; and (3) the Cape Province region, where upwelling is subdued. The highest nutrient contents in thermocline waters are in the northern region, the lowest in the southern one. Wind effects are at a maximum near the southern end of the Namibia upwelling system, and maximum productivity occurs near Walvis Bay, where the product between upwelling rate and nutrient content of upwelled waters is at a maximum. In the Congo/Angola region, opal tends to follow organic carbon quite closely in the Quaternary record. However, organic carbon has a strong precessional component, while opal does not. Despite relatively low opal content, sediments off Angola show the same patterns as those off the Congo; thus, they are part of the same regime. The spectrum shows nonlinear interference patterns between high- and low-latitude forcing, presumably tied to thermocline fertility and wind. On Walvis Ridge, as in the Congo-Angola region, the organic matter record behaves normally; that is, supply is high during glacial periods. In contrast, interglacial periods are favorable for opal deposition. The pattern suggests reduction in silicate content of the thermocline during glacial periods. The reversed phase (opal abundant during interglacials) persists during the entire Pleistocene and can be demonstrated deep into the Pliocene, not just on Walvis Ridge but all the way to the Oranje River and off the Cape Province. From comparison with other regions, it appears that silicate is diminished in the global thermocline, on average, whenever winds become strong enough to substantially shorten the residence time of silicate in upper waters (Walvis Hypothesis, solving the Walvis Paradox of reversed phase in opal deposition). The central discovery during Leg 175 was the documentation of a late Pliocene opal maximum for the entire Namibia upwelling system (early Matuyama Diatom Maximum [MDM]). The maximum is centered on the period between the end of the Gauss Chron and the beginning of the Olduvai Chron. A rather sharp increase in both organic matter deposition and opal deposition occurs near 3 Ma in the middle of the Gauss Chron, in association with a series of major cooling steps. As concerns organic matter, high production persists at least to 1 Ma, when there are large changes in variability, heralding subsequent pulsed production periods. From 3 to 2 Ma, organic matter and opal deposition run more or less parallel, but after 2 Ma opal goes out of phase with organic matter. Apparently, this is the point when silicate becomes limiting to opal production. Thus, the MDM conundrum is solved by linking planetary cooling to increased mixing and upwelling (ramping up to the MDM) and a general removal of silicate from the upper ocean through excess precipitation over global supply (ramping down from the MDM). The hypothesis concerning the origin of the Namibia opal acme or MDM is fundamentally the same as the Walvis Hypothesis, stating that glacial conditions result in removal of silicate from the thermocline (and quite likely from the ocean as a whole, given enough time). The Namibia opal acme, and other opal maxima in the latest Neogene in other regions of the ocean, marks the interval when a cooling ocean selectively removes the abundant silicate inherited from a warm ocean. When the excess silicate is removed, the process ceases. According to the data gathered during Leg 175, major upwelling started in the late part of the late Miocene. Presumably, this process contributed to the drawing down of carbon dioxide from the atmosphere, helping to prepare the way for Northern Hemisphere glaciation.

Stable and strontium isotope ratios and age determination of ODP Site 130-803

Stratigraphic information from strontium, oxygen, and carbon isotopic ratios has been integrated with diatom and planktonic foraminifer datums to refine the Oligocene to early Miocene chemostratigraphy of Site 803. The Sr isotope results are based on analyses of mixed species of planktonic foraminifer and bulk carbonate samples. 87Sr/86Sr ratios of bulk carbonate samples are, in most cases, less radiogenic than contemporaneous seawater. Estimated sediment ages based on planktonic foraminifer 87Sr/86Sr ratios, using the Sr-isotope-age relation determined by Hess and others in 1989, are in moderately good agreement with the biostratigraphic ages. Chronological resolution is significantly enhanced with the correlation of oxygen and carbon isotope records to those of the standard Oligocene section tied to the Geomagnetic Polarity Time Scale at Site 522. Ages revised by this method and other published ages of planktonic foraminifer datums are used to revise the Oligocene stratigraphy of Site 77 to correlate the stable isotope records of Sites 77 and 803. Comparison of the Cibicidoides stable isotope records of Sites 77 and 574 with paleodepths below 2500 m in the central equatorial Pacific, and Site 803 at about 2000-m paleodepth in the Ontong Java Plateau reveals inversions in the vertical d18O gradient at several times during the Oligocene and in the early Miocene. The shallower water site had significantly-higher d18O values than the deeper water sites after the earliest Oligocene 18O enrichment and before 34.5 Ma, in the late Oligocene from 27.5 to at least 25 Ma, and in the early Miocene from 22.5 to 20.5 Ma. It is not possible to ascertain if the d18O inversion persisted during the Oligocene/Miocene transition because the deeper sites have hiatuses spanning this interval. We interpret this pattern to reflect that waters at about 2000 m depth were cold and may have formed from mixing with colder waters originating in northern or southern high-latitude regions. The deeper water appear to have been warmer and may have been a mixture with warm saline waters from mid- or low-latitude regions. No apparent vertical d13C gradient is present during the Oligocene, suggesting that the age difference of these water masses was small.

Age model, stable isotopes, major element proportions and endmember unmixing results from four western tropical Atlantic cores

We investigate changes in the delivery and oceanic transport of Amazon sediments related to terrestrial climate variations over the last 250 ka. We present high-resolution geochemical records from four marine sediment cores located between 5 and 12° N along the northern South American margin. The Amazon River is the sole source of terrigenous material for sites at 5 and 9° N, while the core at 12° N receives a mixture of Amazon and Orinoco detrital particles. Using an endmember unmixing model, we estimated the relative proportions of Amazon Andean material ("%-Andes", at 5 and 9° N) and of Amazon material ("%-Amazon", at 12° N) within the terrigenous fraction. The %-Andes and %-Amazon records exhibit significant precessional variations over the last 250 ka that are more pronounced during interglacials in comparison to glacial periods. High %-Andes values observed during periods of high austral summer insolation reflect the increased delivery of suspended sediments by Andean tributaries and enhanced Amazonian precipitation, in agreement with western Amazonian speleothem records. Increased Amazonian rainfall reflects the intensification of the South American monsoon in response to enhanced land-ocean thermal gradient and moisture convergence. However, low %-Amazon values obtained at 12° N during the same periods seem to contradict the increased delivery of Amazon sediments. We propose that reorganizations in surface ocean currents modulate the northwestward transport of Amazon material. In agreement with published records, the seasonal North Brazil Current retroflection is intensified (or prolonged in duration) during cold substages of the last 250 ka (which correspond to intervals of high DJF or low JJA insolation) and deflects eastward the Amazon sediment and freshwater plume.

Results of a Finite-Element Sea-Ice Ocean Model (FESOM) set-up to study the interannual to decadal variability in the deep-water formation rates

The characteristics of a global set-up of the Finite-Element Sea-Ice Ocean Model under forcing of the period 1958-2004 are presented. The model set-up is designed to study the variability in the deep-water mass formation areas and was therefore regionally better resolved in the deep-water formation areas in the Labrador Sea, Greenland Sea, Weddell Sea and Ross Sea. The sea-ice model reproduces realistic sea-ice distributions and variabilities in the sea-ice extent of both hemispheres as well as sea-ice transport that compares well with observational data. Based on a comparison between model and ocean weather ship data in the North Atlantic, we observe that the vertical structure is well captured in areas with a high resolution. In our model set-up, we are able to simulate decadal ocean variability including several salinity anomaly events and corresponding fingerprint in the vertical hydrography. The ocean state of the model set-up features pronounced variability in the Atlantic Meridional Overturning Circulation as well as the associated mixed layer depth pattern in the North Atlantic deep-water formation areas.

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.