On Provable Algorithms for Determination of Continuous Protein Interdomain Motions from Residual Dipolar Couplings

Dynamics of biomolecules over various spatial and time scales are essential for biological functions such as molecular recognition, catalysis and signaling. However, reconstruction of biomolecular dynamics from experimental observables requires the determination of a conformational probability distribution. Unfortunately, these distributions cannot be fully constrained by the limited information from experiments, making the problem an ill-posed one in the terminology of Hadamard. The ill-posed nature of the problem comes from the fact that it has no unique solution. Multiple or even an infinite number of solutions may exist. To avoid the ill-posed nature, the problem needs to be regularized by making assumptions, which inevitably introduce biases into the result.

Here, I present two continuous probability density function approaches to solve an important inverse problem called the RDC trigonometric moment problem. By focusing on interdomain orientations we reduced the problem to determination of a distribution on the 3D rotational space from residual dipolar couplings (RDCs). We derived an analytical equation that relates alignment tensors of adjacent domains, which serves as the foundation of the two methods. In the first approach, the ill-posed nature of the problem was avoided by introducing a continuous distribution model, which enjoys a smoothness assumption. To find the optimal solution for the distribution, we also designed an efficient branch-and-bound algorithm that exploits the mathematical structure of the analytical solutions. The algorithm is guaranteed to find the distribution that best satisfies the analytical relationship. We observed good performance of the method when tested under various levels of experimental noise and when applied to two protein systems. The second approach avoids the use of any model by employing maximum entropy principles. This 'model-free' approach delivers the least biased result which presents our state of knowledge. In this approach, the solution is an exponential function of Lagrange multipliers. To determine the multipliers, a convex objective function is constructed. Consequently, the maximum entropy solution can be found easily by gradient descent methods. Both algorithms can be applied to biomolecular RDC data in general, including data from RNA and DNA molecules.

A New Method for Modeling Free Surface Flows and Fluid-structure Interaction with Ocean Applications

The computational modeling of ocean waves and ocean-faring devices poses numerous challenges. Among these are the need to stably and accurately represent both the fluid-fluid interface between water and air as well as the fluid-structure interfaces arising between solid devices and one or more fluids. As techniques are developed to stably and accurately balance the interactions between fluid and structural solvers at these boundaries, a similarly pressing challenge is the development of algorithms that are massively scalable and capable of performing large-scale three-dimensional simulations on reasonable time scales. This dissertation introduces two separate methods for approaching this problem, with the first focusing on the development of sophisticated fluid-fluid interface representations and the second focusing primarily on scalability and extensibility to higher-order methods.

We begin by introducing the narrow-band gradient-augmented level set method (GALSM) for incompressible multiphase Navier-Stokes flow. This is the first use of the high-order GALSM for a fluid flow application, and its reliability and accuracy in modeling ocean environments is tested extensively. The method demonstrates numerous advantages over the traditional level set method, among these a heightened conservation of fluid volume and the representation of subgrid structures.

Next, we present a finite-volume algorithm for solving the incompressible Euler equations in two and three dimensions in the presence of a flow-driven free surface and a dynamic rigid body. In this development, the chief concerns are efficiency, scalability, and extensibility (to higher-order and truly conservative methods). These priorities informed a number of important choices: The air phase is substituted by a pressure boundary condition in order to greatly reduce the size of the computational domain, a cut-cell finite-volume approach is chosen in order to minimize fluid volume loss and open the door to higher-order methods, and adaptive mesh refinement (AMR) is employed to focus computational effort and make large-scale 3D simulations possible. This algorithm is shown to produce robust and accurate results that are well-suited for the study of ocean waves and the development of wave energy conversion (WEC) devices.

Influence of Increased Human Presence in the Mills River Basin on Water Availability and Drought

Periods of drought and low streamflow can have profound impacts on both human and natural systems. People depend on a reliable source of water for numerous reasons including potable water supply and to produce economic value through agriculture or energy production. Aquatic ecosystems depend on water in addition to the economic benefits they provide to society through ecosystem services. Given that periods of low streamflow may become more extreme and frequent in the future, it is important to study the factors that control water availability during these times. In the absence of precipitation the slower hydrological response of groundwater systems will play an amplified role in water supply. Understanding the variability of the fraction of streamflow contribution from baseflow or groundwater during periods of drought provides insight into what future water availability may look like and how it can best be managed. The Mills River Basin in North Carolina is chosen as a case-study to test this understanding. First, obtaining a physically meaningful estimation of baseflow from USGS streamflow data via computerized hydrograph analysis techniques is carried out. Then applying a method of time series analysis including wavelet analysis can highlight signals of non-stationarity and evaluate the changes in variance required to better understand the natural variability of baseflow and low flows. In addition to natural variability, human influence must be taken into account in order to accurately assess how the combined system reacts to periods of low flow. Defining a combined demand that consists of both natural and human demand allows us to be more rigorous in assessing the level of sustainable use of a shared resource, in this case water. The analysis of baseflow variability can differ based on regional location and local hydrogeology, but it was found that baseflow varies from multiyear scales such as those associated with ENSO (3.5, 7 years) up to multi decadal time scales, but with most of the contributing variance coming from decadal or multiyear scales. It was also found that the behavior of baseflow and subsequently water availability depends a great deal on overall precipitation, the tracks of hurricanes or tropical storms and associated climate indices, as well as physiography and hydrogeology. Evaluating and utilizing the Duke Combined Hydrology Model (DCHM), reasonably accurate estimates of streamflow during periods of low flow were obtained in part due to the model’s ability to capture subsurface processes. Being able to accurately simulate streamflow levels and subsurface interactions during periods of drought can be very valuable to water suppliers, decision makers, and ultimately impact citizens. Knowledge of future droughts and periods of low flow in addition to tracking customer demand will allow for better management practices on the part of water suppliers such as knowing when they should withdraw more water during a surplus so that the level of stress on the system is minimized when there is not ample water supply.

Carbon and Nitrogen Stable Isotopic Patterns in South Florida Coastal Ecosystems: Modern and Paleoceanographic Perspectives

Long term management plans for restoration of natural flow conditions through the Everglades increase the importance of understanding potential nutrient impacts of increased freshwater delivery on coastal biogeochemistry. The present study sought to increase understanding of the coastal marine system of South Florida under modern conditions and through the anthropogenic changes in the last century, on scales ranging from individual nutrient cycle processes to seasonal patterns in organic material (OM) under varying hydrodynamic regime, to century scale analysis of sedimentary records. In all applications, carbon and nitrogen stable isotopic compositions of OM were examined as natural recorders of change and nutrient cycling in the coastal system. High spatial and temporal variability in stable isotopic compositions were observed on all time scales. During a transient phytoplankton bloom, ä15N values suggested nitrogen fixation as a nutrient source supporting enhanced productivity. Seasonally, particulate organic material (POM) from ten sites along the Florida Reef Tract and in Florida Bay demonstrated variable fluctuations dependent on hydrodynamic setting. Three separate intra-annual patterns were observed, yet statistical differences were observed between groupings of Florida Bay and Atlantic Ocean sites. The POM ä15N values ranged on a quarterly basis by 7‰, while ä13C varied by 22‰. From a sediment history perspective, four cores collected from Florida Bay further demonstrated the spatial and temporal variability of the system in isotopic composition of bulk OM over time. Source inputs of OM varied with location, with terrestrial inputs dominating proximal to Everglades freshwater discharge, seagrasses dominating in open estuary cores, and a marine mixture of phytoplankton and seagrass in a core from the boundary zone between Florida Bay and the Gulf of Mexico. Significant shifts in OM geochemistry were observed coincident with anthropogenic events of the 20th century, including railroad and road construction in the Florida Keys and Everglades, and also the extensive drainage changes in Everglades hydrology. The sediment record also preserved evidence of the major hurricanes of the last century, with excursions in geochemical composition coincident with Category 4-5 storms.

Barite accumulation rates throughout Pliocene-Pleistocene in Eastern Equatorial Pacific at ODP Site 138-849

Export production is an important component of the carbon cycle, modulating the climate system by transferring CO2 from the atmosphere to the deep ocean via the biological pump. Here we use barite accumulation rates to reconstruct export production in the eastern equatorial Pacific over the past 4.3 Ma. We find that export production fluctuated considerably on multiple time scales. Export production was on average higher (51 g C/m**2/yr) during the Pliocene than the Pleistocene (40 g C/m**2/yr), decreasing between 3 and 1 Ma (from more than 60 to 20 g C/m**2/yr) followed by an increase over the last million years. These trends likely reflect basin-scale changes in nutrient inventory and ocean circulation. Our record reveals decoupling between export production and temperatures on these long (million years) time scale. On orbital time scales, export production was generally higher during cold periods (glacial maxima) between 4.3 and 1.1 Ma. This could be due to stronger wind stress and higher upwelling rates during glacial periods. A shift in the timing of maximum export production to deglaciations is seen in the last ~1.1 million years. Results from this study suggest that, in the eastern equatorial Pacific, mechanisms that affect nutrient supply and/or ecosystem structure and in turn carbon export on orbital time scales differ from those operating on longer time scales and that processes linking export production and climate-modulated oceanic conditions changed about 1.1 million years ago. These observations should be accounted for in climate models to ensure better predictions of future climate change.

Middle Eocene Fe and Ca counts, Ca/Fe ratios, revised splices and astronomical age models form different Sites of IODP Exp342

Accurate age models are a tool of utmost important in paleoclimatology. Constraining the rate and pace of past climate change are at the core of paleoclimate research, as such knowledge is crucial to our understanding of the climate system. Indeed, it allows for the disentanglement of the various drivers of climate change. The scarcity of highly resolved sedimentary records from the middle Eocene (Bartonian - Lutetian Stages; 47.8 - 37.8 Ma) has led to the existence of the "Eocene astronomical time scale gap" and hindered the establishment of a comprehensive astronomical time scale (ATS) for the entire Cenozoic. Sediments from the Newfoundland Ridge drilled during Integrated Ocean Drilling Program (IODP) Expedition 342 span the Eocene gap at an unprecedented stratigraphic resolution with carbonate bearing sediments. Moreover, these sediments exhibit cyclic lithological changes that allow for an astronomical calibration of geologic time. In this study, we use the dominant obliquity imprint in XRF-derived calcium-iron ratio series (Ca/Fe) from three sites drilled during IODP Expedition 342 (U1408, U1409, U1410) to construct a floating astrochronology. We then anchor this chronology to numerical geological time by tuning 173-kyr cycles in the amplitude modulation pattern of obliquity to an astronomical solution. This study is one of the first to use the 173-kyr obliquity amplitude cycle for astrochronologic purposes, as previous studies primarily use the 405-kyr long eccentricity cycle as a tuning target to calibrate the Paleogene geologic time scale. We demonstrate that the 173-kyr cycles in obliquity's amplitude are stable between 40 and 50 Ma, which means that one can use the 173-kyr cycle for astrochronologic calibration in the Eocene. Our tuning provides new age estimates for magnetochron reversals C18n.1n - C21r and a stratigraphic framework for key sites from Expedition 342 for the Eocene. Some disagreements emerge when we compare our tuning for the interval between C19r and C20r with previous tuning attempts from the South Atlantic. We therefore present a revision of the original astronomical interpretations for the latter records, so that the various astrochronologic age models for the middle Eocene in the North- and South-Atlantic are consistent.

Benthic foraminiferal record from Santa Barbara basin

The oxygen minimum zone (OMZ) of the late Quaternary California margin experienced abrupt and dramatic changes in strength and depth in response to changes in intermediate water ventilation, ocean productivity, and climate at orbital through millennial time scales. Expansion and contraction of the OMZ is exhibited at high temporal resolution (107-126 year) by quantitative benthic foraminiferal assemblage changes in two piston cores forming a vertical profile in Santa Barbara Basin (569 m, basin floor; 481 m, near sill depth) to 34 and 24 ka, respectively. Variation in the OMZ is quantified by new benthic foraminiferal groupings and new dissolved oxygen index based on documented relations between species and water-mass oxygen concentrations. Foraminiferal-based paleoenvironmental assessments are integrated with principal component analysis, bioturbation, grain size, CaCO3, total organic carbon, and d13C to reconstruct basin oxygenation history. Fauna responded similarly between the two sites, although with somewhat different magnitude and taxonomic expression. During cool episodes (Younger Dryas and stadials), the water column was well oxygenated, most strongly near the end of the glacial episode (17-16 ka; Heinrich 1). In contrast, the OMZ was strong during warm episodes (Bølling/Allerød, interstadials, and Pre-Boreal). During the Bølling/Allerød, the OMZ shoaled to <360 m of contemporaneous sea level, its greatest vertical expansion of the last glacial cycle. Assemblages were then dominated by Bolivina tumida, reflecting high concentrations of dissolved methane in bottom waters. Short decadal intervals were so severely oxygen-depleted that no benthic foraminifera were present. The middle to late Holocene (6-0 ka) was less dysoxic than the early Holocene.

Radiocarbon, stable isotopic compositions and seawater chemical properties of Black Sea core-top sediments and CTD stations

Chloropigments and their derivative pheopigments preserved in sediments can directly be linked to photosynthesis. Their carbon and nitrogen stable isotopic compositions have been shown to be a good recorder of recent and past surface ocean environmental conditions tracing the carbon and nitrogen sources and dominant assimilation processes of the phytoplanktonic community. In this study we report results from combined compound-specific radiocarbon and stable carbon and nitrogen isotope analysis to examine the time-scales of synthesis and fate of chlorophyll-a and its degradation products pheophytin-a, pyropheophytin-a, and 132,173-cyclopheophorbide-a-enol until burial in Black Sea core-top sediments. The pigments are mainly of marine phytoplanktonic origin as implied by their stable isotopic compositions. Pigment ?15N values indicate nitrate as the major uptake substrate but 15N-depletion towards the open marine setting indicates either contribution from N2-fixation or direct uptake of ammonium from deeper waters. Radiocarbon concentrations translate into minimum and maximum pigment ages of approximately 40 to 1200 years. This implies that protective mechanisms against decomposition such as association with minerals, storage in deltaic anoxic environments, or eutrophication-induced hypoxia and light limitation are much more efficient than previously thought. Moreover, seasonal variations of nutrient source, growth period, and habitat and their associated isotopic variability are likely at least as strong as long-term trends. Combined triple isotope analysis of sedimentary chlorophyll and its primary derivatives is a powerful tool to delineate biogeochemical and diagenetic processes in the surface water and sediments, and to assess their precise time-scales.

Compilation of organic carbon distribution and sedimentology in the surface sediments on the continental margin offshore southwestern Africa

In this study we demonstrate the relevance of lateral particle transport in nepheloid layers for organic carbon (OC) accumulation and burial across high-productive continental margins. We present geochemical data from surface sediments and suspended particles in the bottom nepheloid layer (BNL) from the most productive coastal upwelling area of the modern ocean, the Benguela upwelling system offshore southwest Africa. Interpretation of depositional patterns and comparison of downslope trends in OC content, organic matter composition, and 14C age between suspended particles and surface sediments indicate that lateral particle transport is the primary mechanism controlling supply and burial of OC. We propose that effective seaward particle transport primarily along the BNL is a key process that promotes and maintains local high sedimentation rates, ultimately causing high preservation of OC in a depocenter on the upper slope offshore Namibia. As lateral transport efficiently displaces areas of enhanced OC burial from maximum production at highly productive continental margins, vertical particle flux models do not sufficiently explain the relationship between primary production and shallow-marine OC burial. On geologic time scales, the widest distribution and strongest intensity of lateral particle transport is expected during periods of rapid sea-level change. At times in the geologic past, widespread downslope lateral transport of OC thus may have been a primary driver of enhanced OC burial at deeper continental slopes and abyssal basins.

(Fig. 1D, 3E) Sea surface temperature and alkenone C37:4 for sediment core MD01-2444 during the penultimate glaciation

Past glacials can be thought of as natural experiments in which variations in boundary conditions influenced the character of climate change. However, beyond the last glacial, an integrated view of orbital- and millennial-scale changes and their relation to the record of glaciation has been lacking. Here, we present a detailed record of variations in the land-ocean system from the Portuguese margin during the penultimate glacial and place it within the framework of ice-volume changes, with particular reference to European ice-sheet dynamics. The interaction of orbital- and millennial-scale variability divides the glacial into an early part with warmer and wetter overall conditions and prominent climate oscillations, a transitional mid-part, and a late part with more subdued changes as the system entered a maximum glacial state. The most extreme event occurred in the mid-part and was associated with melting of the extensive European ice sheet and maximum discharge from the Fleuve Manche river. This led to disruption of the meridional overturning circulation, but not a major activation of the bipolar seesaw. In addition to stadial duration, magnitude of freshwater forcing, and background climate, the evidence also points to the influence of the location of freshwater discharges on the extent of interhemispheric heat transport.

Nd isotopes an geochemistry of bulk deep sea sediments of the Atlantic Ocean

Nd isotopes preserved in fossil fish teeth and ferromanganese crusts have become a common tool for tracking variations in water mass composition and circulation through time. Studies of Nd isotopes extracted from Pleistocene to Holocene bulk sediments using hydroxylamine hydrochloride (HH) solution yield high resolution records of Nd isotopes that can be interpreted in terms of deep water circulation, but concerns about diagenesis and potential contamination of the seawater signal limit application of this technique to geologically young samples. In this study we demonstrate that Nd extracted from the > 63 µm, decarbonated fraction of older Ocean Drilling Program (ODP) sediments using a 0.02 M HH solution produces Nd isotopic ratios that are within error of values from cleaned fossil fish teeth collected from the same samples, indicating that the HH-extractions are robust recorders of deep sea Nd isotopes. This excellent correlation was achieved for 94 paired fish teeth and HH-extraction samples ranging in age from the Miocene to Cretaceous, distributed throughout the north, tropical and south Atlantic, and composed of a range of lithologies including carbonate-rich oozes/chalks and black shales. The strong Nd signal recovered from Cretaceous anoxic black shale sequences is unlikely to be associated with ferromanganese oxide coatings, but may be derived from abundant phosphatic fish teeth and debris or organic matter in these samples. In contrast to the deep water Nd isotopic signal, Sr isotopes from HH-extractions are often offset from seawater values, suggesting that evaluation of Sr isotopes is a conservative test for the integrity of Nd isotopes in the HH fraction. However, rare earth elements (REE) from the HH-extractions and fish teeth produce distinctive middle REE bulge patterns that may prove useful for evaluating whether the Nd isotopic signal represents uncontaminated seawater. Alternatively, a few paired HH-extraction and cleaned fish teeth samples from each site of interest can be used to verify the seawater composition of the HH-extractions. The similarity between isotopic values for the HH-extraction and fish teeth illustrates that the extensive cleaning protocol applied to fish teeth samples is not necessary in typical, carbonate-rich, deep sea sediments.

Geochemistry of sediment samples from ODP Leg 135 holes

The results of inductively coupled argon plasma (ICAP) chemical analyses carried out on some 300 core samples from Ocean Drilling Program Sites 834, 835, 838, and 839 are presented. These sites were drilled during Leg 135 in the Lau Basin. The data are compared with total gamma (SGR) wireline logs at Sites 834 and 835. Pliocene (Piacenzian) nannofossil Zone CN12, which has been identified at Sites 834 and 835, is examined in detail using spectral analyses on core and wireline logs. The potassium and calcium concentrations from the core material were used to calculate an objective depth-to-geological time stretching function, which improved the stratigraphic correlation between sites. The integrated use of chemical analyses, wireline-log data and paleomagnetic results improved confidence in the correlations obtained. Although no significant sedimentation periodicities were obtained from the two sites, a common concentration of energy between 30 and 60 k.y. was recorded.

(Table 1) Age determination of North Atlantic sediment cores

Stable isotope, trace metal, alkenone paleothermometry, and radiocarbon methods have been applied to sediment cores in the western subpolar North Atlantic between Hudson Strait and Cape Hatteras to reveal the history of climate in that region over the past ~11 kyr. We focus on cores from the Laurentian Fan, which is known to have rapid and continuous accumulation of hemipelagic sediment. Although results among our various proxy data are not always in agreement, the weight of the evidence (alkenone sea surface temperature (SST), d18O and abundance of Globigerinoides ruber) indicates a continual cooling of surface waters over Laurentian Fan, from about 18°C in the early Holocene to about 8°C today. Alternatively, Mg/Ca data on planktonic foraminifera indicate no systematic change in Holocene SST. The inferred long-term decrease in SST was probably driven by decreasing seasonality of Northern Hemisphere insolation. Two series of proxy data show the gradual cooling was interrupted by a two-step cold pulse that began 8500 years ago, and lasted about 700 years. Although this event is associated with the final deglaciation of Hudson Bay, there is no d18O minimum anywhere in the Labrador Sea, yet there is some evidence for it as far south as Cape Hatteras. Finally, although the 8200 year B.P. event has been implicated in decreasing North Atlantic ventilation, and hence widespread temperature depression on land and at sea, we find inconsistent evidence for a change at that time in deep ocean nutrient content at ~4 km water depth.

Stable isotopes, radionuclides, and calculated sea surface temperature of two sediment cores in the Sicily Channel

The Eastern Mediterranean Transient (EMT) occurred in the Aegean Sea from 1988 to 1995 and is the most significant intermediate-to-deep Mediterranean overturning perturbation reported by instrumental records. The EMT was likely caused by accumulation of high salinity waters in the Levantine and enhanced heat loss in the Aegean Sea, coupled with surface water freshening in the Sicily Channel. It is still unknown whether similar transients occurred in the past and, if so, what their forcing processes were. In this study, sediments from the Sicily Channel document surface water freshening (SCFR) at 1910±12, 1812±18, 1725±25 and 1580±30 CE. A regional ocean hindcast links SCFR to enhanced deep-water production and in turn to strengthened Mediterranean thermohaline circulation. Independent evidence collected in the Aegean Sea supports this reconstruction, showing that enhanced bottom water ventilation in the Eastern Mediterranean was associated with each SCFR event. Comparison between the records and multi-decadal atmospheric circulation patterns and climatic external forcings indicates that Mediterranean circulation destabilisation occurs during positive North Atlantic Oscillation (NAO) and negative Atlantic Multidecadal Oscillation (AMO) phases, reduced solar activity and strong tropical volcanic eruptions. They may have recurrently produced favourable deep-water formation conditions, both increasing salinity and reducing temperature on multi-decadal time scales.

Stable isotope and Mg/Ca record of IODP Site 339-U1385

Although millennial-scale climate variability (<10 ka) has been well studied during the last glacial cycles, little is known about this important aspect of climate in the early Pleistocene, prior to the Middle Pleistocene Transition. Here we present an early Pleistocene climate record at centennial resolution for two representative glacials (marine isotope stages (MIS) 37-41 from approximately 1235 to 1320 ka) during the "41 ka world" at Integrated Ocean Drilling Program Site U1385 (the "Shackleton Site") on the southwest Iberian margin. Millennial-scale climate variability was suppressed during interglacial periods (MIS 37, MIS 39, and MIS 41) and activated during glacial inceptions when benthic d18O exceeded 3.2 per mil. Millennial variability during glacials MIS 38 and MIS 40 closely resembled Dansgaard-Oeschger events from the last glacial (MIS 3) in amplitude, shape, and pacing. The phasing of oxygen and carbon isotope variability is consistent with an active oceanic thermal bipolar see-saw between the Northern and Southern Hemispheres during most of the prominent stadials. Surface cooling was associated with systematic decreases in benthic carbon isotopes, indicating concomitant changes in the meridional overturning circulation. A comparison to other North Atlantic records of ice rafting during the early Pleistocene suggests that freshwater forcing, as proposed for the late Pleistocene, was involved in triggering or amplifying perturbations of the North Atlantic circulation that elicited a bipolar see-saw response. Our findings support similarities in the operation of the climate system occurring on millennial time scales before and after the Middle Pleistocene Transition despite the increases in global ice volume and duration of the glacial cycles.