ALBERT

Beschreibung: The relationships between trachytes and peralkaline rhyolites (i.e. pantellerites and comendites), which occur in many continental rift systems, oceanic islands and continental intraplate settings, is unclear. To fill this gap, we have performed phase equilibrium experiments on two representative metaluminous trachytes from Pantelleria to determine both their pre-eruptive equilibration conditions (pressure, temperature, H2O content and redox state) and liquid lines of descent. Experiments were performed in the temperature range 750–950 C, pressure 0 5–1 5 kbar and fluid saturation conditions with XH2O [¼H2O/(H2OþCO2)] ranging between zero and unity. Redox conditions were fixed below the nickel–nickel oxide buffer (NNO). The results show that at 950 C and melt water contents (H2Omelt) close to saturation, trachytes are at liquidus conditions at all pressures. Clinopyroxene is the liquidus phase, being followed by iron-rich olivine and alkali feldspar. Comparison of experimental and natural phases (abundances and compositions) yields the following pre-eruptive conditions: P¼160 5 kbar, T¼925625 C, H2Omelt¼261wt %, and fO2 between NNO– 0 5 and NNO– 2. A decrease in temperature from 950 C to 750 C, as well as of H2Omelt, promotes a massive crystallization of alkali feldspar to over 80 wt %. Iron-bearing minerals show gradual iron enrichment when T and fO2 decrease, trending towards the compositions of the phenocrysts of natural pantellerites. Despite the metaluminous character of the bulk-rock compositions, residual glasses obtained after 80 wt % crystallization evolve toward comenditic compositions, owing to profuse alkali feldspar crystallization, which decreases the Al2O3 of the melt, leading to a consequent increase in the peralkalinity index [PI¼molar (Na2OþK2O)/Al2O3]. This is the first experimental demonstration that peralkaline felsic derivatives can be produced by low-pressure fractional crystallization of metaluminous mafic magmas. Our results show that the pantelleritic magmas of basalt–trachyte–rhyolite igneous suites require at least 95 wt % of parental basalt crystallization, consistent with trace element evidence. Redox conditions, through their effect on Fe–Ti oxide stabilities, control the final iron content of the evolving melt.

Beschreibung: Published

Beschreibung: 559- 588

Beschreibung: 2V. Struttura e sistema di alimentazione dei vulcani

Beschreibung: JCR Journal

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Bulletin of the American Meteorological Society 99, Suppl. S (2018): S21-S26, doi:10.1175/BAMS-D-17-0128.1.

Beschreibung: NOAA Coral Reef Conservation Program; National Science Foundation OCE 1537338, OCE 1605365, OCE 1031971

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 1555-1566, doi:10.1175/JPO-D-17-0231.1.

Beschreibung: A primary challenge in modeling flow over shallow coral reefs is accurately characterizing the bottom drag. Previous studies over continental shelves and sandy beaches suggest surface gravity waves should enhance the drag on the circulation over coral reefs. The influence of surface gravity waves on drag over four platform reefs in the Red Sea is examined using observations from 6-month deployments of current and pressure sensors burst sampling at 1Hz for 4–5min. Depth-average current fluctuations U0 within each burst are dominated by wave orbital velocities uw that account for 80%–90%of the burst variance and have a magnitude of order 10 cm s21, similar to the lower-frequency depth-average current Uavg. Previous studies have shown that the cross-reef bottom stress balances the pressure gradient over these reefs. A bottom stress estimate that neglects the waves (rCdaUavgjUavgj, where r is water density and Cda is a drag coefficient) balances the observed pressure gradient when uw is smaller than Uavg but underestimates the pressure gradient when uw is larger than Uavg (by a factor of 3–5 when uw 5 2Uavg), indicating the neglected waves enhance the bottom stress. In contrast, a bottom stress estimate that includes the waves [rCda(Uavg 1 U0)jUavg 1 U0j)] balances the observed pressure gradient independent of the relative size of uw and Uavg, indicating that this estimate accounts for the wave enhancement of the bottom stress. A parameterization proposed by Wright and Thompson provides a reasonable estimate of the total bottom stress (including the waves) given the burst-averaged current and the wave orbital velocity.

Beschreibung: The Red Sea field program was supported by Awards USA 00002 and KSA 00011 made by KAUST. S. Lentz was supported for the analysis by NSF Award OCE-1558343.

Beschreibung: 2019-01-13

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Bulletin of the American Meteorological Society 99 (2018): 513-520, doi:10.1175/BAMS-D-16-0323.1.

Beschreibung: Seasonally ice-covered marginal seas are among the most difficult regions in the Arctic to study. Physical constraints imposed by the variable presence of sea ice in all stages of growth and melt make the upper water column and air–sea ice interface especially challenging to observe. At the same time, the flow of solar energy through Alaska’s marginal seas is one of the most important regulators of their weather and climate, sea ice cover, and ecosystems. The deficiency of observing systems in these areas hampers forecast services in the region and is a major contributor to large uncertainties in modeling and related climate projections. The Arctic Heat Open Science Experiment strives to fill this observation gap with an array of innovative autonomous floats and other near-real-time weather and ocean sensing systems. These capabilities allow continuous monitoring of the seasonally evolving state of the Chukchi Sea, including its heat content. Data collected by this project are distributed in near–real time on project websites and on the Global Telecommunications System (GTS), with the objectives of (i) providing timely delivery of observations for use in weather and sea ice forecasts, for model, and for reanalysis applications and (ii) supporting ongoing research activities across disciplines. This research supports improved forecast services that protect and enhance the safety and economic viability of maritime and coastal community activities in Alaska. Data are free and open to all (see www.pmel.noaa.gov/arctic-heat/).

Beschreibung: This work was supported by NOAA Ocean and Atmospheric Research and the Joint Institute for the Study of the Atmosphere and Ocean (JISAO) under NOAA Cooperative Agreement NA15OAR4320063 and by the Innovative Technology for Arctic Exploration (ITAE) program at JISAO/PMEL. Jayne, Robbins, and Ekholm were supported by ONR (N00014-12-10110).

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Atmospheric and Oceanic Technology 35 (2018): 893-910, doi:10.1175/JTECH-D-17-0102.1.

Beschreibung: Rotary sidescan sonars are widely used to image the seabed given their high temporal and spatial resolution. This high resolution is necessary to resolve bedform dynamics and evolution; however, sidescan sonars do not directly measure bathymetry, limiting their utility. When sidescan sonars are mounted close to the seabed, bedforms may create acoustical “shadows” that render previous methods that invert the backscatter intensity to estimate bathymetry and are based on the assumption of a fully insonified seabed ineffective. This is especially true in coastal regions, where bedforms are common features whose large height relative to the water depth may significantly influence the surrounding flow. A method is described that utilizes sonar shadows to estimate bedform height and asymmetry. The method accounts for the periodic structure of bedform fields and the projection of the shadows onto adjacent bedforms. It is validated with bathymetric observations of wave-orbital ripples, with wavelengths ranging from 0.3 to 0.8 m, and tidally reversing megaripples, with wavelengths from 5 to 8 m. In both cases, bathymetric-measuring sonars were deployed in addition to a rotary sidescan sonar to provide a ground truth; however, the bathymetric sonars typically measure different and smaller areas than the rotary sidescan sonar. The shadow-based method and bathymetric-measuring sonar data produce estimates of bedform height that agree by 34.0% ± 27.2% for wave-orbital ripples and 16.6% ± 14.7% for megaripples. Errors for estimates of asymmetry are 1.9% ± 2.1% for wave-orbital ripples and 11.2% ± 9.6% for megaripples.

Beschreibung: This project is partially supported by the National Science Foundation through a Graduate Research Fellowship and a Massachusetts Institute of Technology Energy Initiative Fellowship. Additionally, funding used in developing the method was obtained from NSF Grants OCE-1634481 and OCE-1635151. Field work was funded under ONR Grants N00014-06-10329 and N00014-13-1-0364.

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 2127-2140, doi:10.1175/JPO-D-18-0035.1.

Beschreibung: Shipboard hydrographic and velocity measurements collected in summer 2014 are used to study the evolution of the freshwater coastal current in southern Greenland as it encounters Cape Farewell. The velocity structure reveals that the coastal current maintains its identity as it flows around the cape and bifurcates such that most of the flow is diverted to the outer west Greenland shelf, while a small portion remains on the inner shelf. Taking into account this inner branch, the volume transport of the coastal current is conserved, but the freshwater transport decreases on the west side of Cape Farewell. A significant amount of freshwater appears to be transported off the shelf where the outer branch flows adjacent to the shelfbreak circulation. It is argued that the offshore transposition of the coastal current is caused by the flow following the isobaths as they bend offshore because of the widening of the shelf on the west side of Cape Farewell. An analysis of the potential vorticity shows that the subsequent seaward flux of freshwater can be enhanced by instabilities of the current. This set of circumstances provides a pathway for the freshest water originating from the Arctic, as well as runoff from the Greenland ice sheet, to be fluxed into the interior Labrador Sea where it could influence convection in the basin.

Beschreibung: Funding for this project was provided by the National Science Foundation under Grant OCE-1259618.

Beschreibung: 2019-03-11

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 435-453, doi:10.1175/JPO-D-17-0122.1.

Beschreibung: Observations of surface waves, currents, and turbulence at the Columbia River mouth are used to investigate the source and vertical structure of turbulence in the surface boundary layer. Turbulent velocity data collected on board freely drifting Surface Wave Instrument Float with Tracking (SWIFT) buoys are corrected for platform motions to estimate turbulent kinetic energy (TKE) and TKE dissipation rates. Both of these quantities are correlated with wave steepness, which has previously been shown to determine wave breaking within the same dataset. Estimates of the turbulent length scale increase linearly with distance from the free surface, and roughness lengths estimated from velocity statistics scale with significant wave height. The vertical decay of turbulence is consistent with a balance between vertical diffusion and dissipation. Below a critical depth, a power-law scaling commonly applied in the literature works well to fit the data. Above this depth, an exponential scaling fits the data well. These results, which are in a surface-following reference frame, are reconciled with results from the literature in a fixed reference frame. A mapping between free-surface and mean-surface reference coordinates suggests 30% of the TKE is dissipated above the mean sea surface.

Beschreibung: Funding for this project was provided by the Office of Naval Research as part of the RIVET-II DRI, and for the DARLA group.

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 573-590, doi:10.1175/JPO-D-17-0206.1.

Beschreibung: Motivated by the proximity of the Northern Recirculation Gyre and the deep western boundary current in the North Atlantic, an idealized model is used to investigate how recirculation gyres and a deep flow along a topographic slope interact. In this two-layer quasigeostrophic model, an unstable jet imposed in the upper layer generates barotropic recirculation gyres. These are maintained by an eddy-mean balance of potential vorticity (PV) in steady state. The authors show that the topographic slope can constrain the northern recirculation gyre meridionally and that the gyre’s adjustment to the slope leads to increased eddy PV fluxes at the base of the slope. When a deep current is present along the topographic slope in the lower layer, these eddy PV fluxes stir the deep current and recirculation gyre waters. Increased proximity to the slope dampens the eddy growth rate within the unstable jet, altering the geometry of recirculation gyre forcing and leading to a decrease in overall eddy PV fluxes. These mechanisms may shape the circulation in the western North Atlantic, with potential feedbacks on the climate system.

Beschreibung: We gratefully acknowledge an AMS graduate fellowship (IALB) and U.S. National Science Foundation Grants OCE-1332667 and 1332834 (IALB and JMT).

Beschreibung: 2018-09-06

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 739-748, doi:10.1175/JPO-D-17-0089.1.

Beschreibung: McDougall and Ferrari have estimated the global deep upward diapycnal flow in the boundary layer overlying continental slopes that must balance both downward diapycnal flow in the deep interior and the formation of bottom water around Antarctica. The decrease of perimeter of isopycnal surfaces with depth and the observed decay with height above bottom of turbulent dissipation in the deep ocean play a key role in their estimate. They argue that because the perimeter of seamounts increases with depth, the net effect of mixing around seamounts is to produce net downward diapycnal flow. While this is true along much of a seamount, it is shown here that diapycnal flow of the densest water around the seamount is upward, with buoyancy being transferred from water just above. The same is true for midocean ridges, whose perimeter is constant with depth. It is argued that mixing around seamounts and especially midocean ridges contributes positively to the global deep overturning circulation, reducing the amount of turbulence demanded over the continental slopes to balance the buoyancy budget for the bottom and deep water.

Beschreibung: This work was supported by National Science Foundation Grant OCE- 1232962.

Beschreibung: 2018-09-29

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 643-646, doi:10.1175/JPO-D-17-0240.1.

Beschreibung: A simple oceanic model is presented for source–sink flow on the β plane to discuss the pathways from source to sink when transport boundary layers have large enough Reynolds numbers to be inertial in their dynamics. A representation of the flow as a Fofonoff gyre, suggested by prior work on inertial boundary layers and eddy-driven circulations in two-dimensional turbulent flows, indicates that even when the source and sink are aligned along the same western boundary the flow must intrude deep into the interior before exiting at the sink. The existence of interior pathways for the flow is thus an intrinsic property of an inertial circulation and is not dependent on particular geographical basin geometry.

Beschreibung: 2018-09-12

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 479-509, doi:10.1175/JPO-D-16-0283.1.

Beschreibung: Lateral submesoscale processes and their influence on vertical stratification at shallow salinity fronts in the central Bay of Bengal during the winter monsoon are explored using high-resolution data from a cruise in November 2013. The observations are from a radiator survey centered at a salinity-controlled density front, embedded in a zone of moderate mesoscale strain (0.15 times the Coriolis parameter) and forced by winds with a downfront orientation. Below a thin mixed layer, often ≤10 m, the analysis shows several dynamical signatures indicative of submesoscale processes: (i) negative Ertel potential vorticity (PV); (ii) low-PV anomalies with O(1–10) km lateral extent, where the vorticity estimated on isopycnals and the isopycnal thickness are tightly coupled, varying in lockstep to yield low PV; (iii) flow conditions susceptible to forced symmetric instability (FSI) or bearing the imprint of earlier FSI events; (iv) negative lateral gradients in the absolute momentum field (inertial instability); and (v) strong contribution from differential sheared advection at O(1) km scales to the growth rate of the depth-averaged stratification. The findings here show one-dimensional vertical processes alone cannot explain the vertical stratification and its lateral variability over O(1–10) km scales at the radiator survey.

Beschreibung: S. Ramachandran acknowledges support from the National Science Foundation through award OCE 1558849 and the U.S. Office of Naval Research, Grants N00014-13-1-0456 and N00014-17- 1-2355. A. Tandon acknowledges support from the U.S. Office of Naval Research, Grants N00014-13-1-0456 and N00014-17-1-2355. J. T. Farrar and R. A. Weller were supported by the U.S. Office of Naval Research, Grant N00014-13-1-0453, to collect the UCTD data and process theUCTD and shipboard meteorological data. J. Nash, J. Mackinnon, and A. F. Waterhouse acknowledge support from the U. S. Office of Naval Research, Grants N00014-13-1-0503 and N00014-14-1-0455. E. Shroyer acknowledges support from the U. S. Office of Naval Research, Grants N00014-14-10236 and N00014-15- 12634. A. Mahadevan acknowledges support fromthe U. S. Office of Naval Research, Grant N00014-13-10451. A. J. Lucas and R. Pinkel acknowledge support from the U. S. Office of Naval Research, Grant N00014-13-1-0489.

Beschreibung: 2018-08-26

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 607-623, doi:10.1175/JPO-D-17-0189.1.

Beschreibung: The roles of straining and dissipation in controlling stratification are derived analytically using a vertical salinity variance method. Stratification is produced by converting horizontal variance to vertical variance via straining, that is, differential advection of horizontal salinity gradients, and stratification is destroyed by the dissipation of vertical variance through turbulent mixing. A numerical model is applied to the Changjiang estuary in order to demonstrate the salinity variance balance and how it reveals the factors controlling stratification. The variance analysis reveals that dissipation reaches its maximum during spring tide in the Changjiang estuary, leading to the lowest stratification. Stratification increases from spring tide to neap tide because of the increasing excess of straining over dissipation. Throughout the spring–neap tidal cycle, straining is almost always larger than dissipation, indicating a net excess of production of vertical variance relative to dissipation. This excess is balanced on average by advection, which exports vertical variance out of the estuarine region into the plume. During neap tide, tidal straining shows a general tendency of destratification during the flood tide and restratification during ebb, consistent with the one-dimensional theory of tidal straining. During spring tide, however, positive straining occurs during flood because of the strong baroclinicity induced by the intensified horizontal salinity gradient. These results indicate that the salinity variance method provides a valuable approach for examining the spatial and temporal variability of stratification in estuaries and coastal environments.

Beschreibung: X. Li was supported by the China Scholarship Council. W. R. Geyer was supported by NSF Grants OCE 1736539 and OCE 1634480. J. Zhu was supported by the National Natural Science Foundation of China (41476077 and 41676083). H. Wu was supported by the National Natural Science Foundation of China (41576088 and 41776101).

Beschreibung: 2018-09-08

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 31 (2018): 6245-6261, doi:10.1175/JCLI-D-17-0513.1.

Beschreibung: Reconstructions of sea surface temperature (SST) based on instrumental observations suggest that the equatorial Pacific zonal SST gradient has increased over the twentieth century. While this increase is suggestive of the ocean dynamical thermostat mechanism of Clement et al., observations of a concurrent weakening of the zonal atmospheric (Walker) circulation are not. Here we show, using heat and momentum budget calculations on an ocean reanalysis dataset, that a seasonal weakening of the zonal atmospheric circulation is in fact consistent with cooling in the eastern equatorial Pacific (EEP) and thus an increase in the zonal SST gradient. This cooling is driven by a strengthening Equatorial Undercurrent (EUC) in response to decreased upper-ocean westward momentum associated with weakening equatorial zonal wind stress. This process can help to reconcile the seemingly contradictory twentieth-century trends in the tropical Pacific atmosphere and ocean. Moreover, it is shown that coupled general circulation models (CGCMs) do not correctly simulate this process; we identify a systematic bias in the relationship between changes in equatorial surface zonal wind stress in the EEP and EUC strength that may help to explain why observations and CGCMs have opposing trends in the zonal SST gradient over the twentieth century.

Beschreibung: 2019-01-11

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 1831-1848, doi:10.1175/JPO-D-18-0068.1.

Beschreibung: We present a simplified theory using reduced-gravity equations for North Atlantic Deep Water (NADW) and its variation driven by high-latitude deep-water formation. The theory approximates layer thickness on the eastern boundary with domain-averaged layer thickness and, in tandem with a mass conservation argument, retains fundamental physics for cross-equatorial flows on interannual and longer forcing time scales. Layer thickness anomalies are driven by a time-dependent northern boundary condition that imposes a southward volume flux representative of a variable source of NADW and damped by diapycnal mixing throughout the basin. Moreover, an outflowing southern boundary condition imposes a southward volume flux that generally differs from the volume flux at the northern boundary, giving rise to temporal storage of NADW within the Atlantic basin. Closed form analytic solutions for the amplitude and phase are provided when the variable source of NADW is sinusoidal. We provide a nondimensional analysis that demonstrates that solution behavior is primarily controlled by two parameters that characterize the meridional extent of the southern basin and the width of the basin relative to the equatorial deformation radius. Similar scaling applied to the time-lagged equations of Johnson and Marshall provides a clear connection to their results. Numerical simulations of reduced-gravity equations agree with analytic predictions in linear, turbulent, and diabatic regimes. The theory introduces a simple analytic framework for studying idealized buoyancy- and wind-driven cross-equatorial flows on interannual and longer time scales.

Beschreibung: This research was supported by the National Science Foundation under Grant OCE- 1634468.

Beschreibung: 2019-02-15

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 31 (2018): 7565-7581, doi:10.1175/JCLI-D-18-0108.1.

Beschreibung: There is mounting evidence that the width of the tropics has increased over the last few decades, but there are large differences in reported expansion rates. This is, likely, in part due to the wide variety of metrics that have been used to define the tropical width. Here we perform a systematic investigation into the relationship among nine metrics of the zonal-mean tropical width using preindustrial control and abrupt quadrupling of CO2 simulations from a suite of coupled climate models. It is shown that the latitudes of the edge of the Hadley cell, the midlatitude eddy-driven jet, the edge of the subtropical dry zones, and the Southern Hemisphere subtropical high covary interannually and exhibit similar long-term responses to a quadrupling of CO2. However, metrics based on the outgoing longwave radiation, the position of the subtropical jet, the break in the tropopause, and the Northern Hemisphere subtropical high have very weak covariations with the above metrics and/or respond differently to increases in CO2 and thus are not good indicators of the expansion of the Hadley cell or subtropical dry zone. The differing variability and responses to increases in CO2 among metrics highlights that care is needed when choosing metrics for studies of the width of the tropics and that it is important to make sure the metric used is appropriate for the specific phenomena and impacts being examined.

Beschreibung: DW acknowledges support from NSF AGS-1403676.

Beschreibung: 2019-02-08

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 2457-2475, doi:10.1175/JPO-D-17-0186.1.

Beschreibung: A subpolar marginal sea, like the Nordic seas, is a transition zone between the temperature-stratified subtropics (the alpha ocean) and the salinity-stratified polar regions (the beta ocean). An inflow of Atlantic Water circulates these seas as a boundary current that is cooled and freshened downstream, eventually to outflow as Deep and Polar Water. Stratification in the boundary region is dominated by a thermocline over the continental slope and a halocline over the continental shelves, separating Atlantic Water from Deep and Polar Water, respectively. A conceptual model is introduced for the circulation and water mass transformation in a subpolar marginal sea to explore the potential interaction between the alpha and beta oceans. Freshwater input into the shelf regions has a slight strengthening effect on the Atlantic inflow, but more prominently impacts the water mass composition of the outflow. This impact of freshwater, characterized by enhancing Polar Water outflow and suppressing Deep Water outflow, is strongly determined by the source location of freshwater. Concretely, perturbations in upstream freshwater sources, like the Baltic freshwater outflow into the Nordic seas, have an order of magnitude larger potential to impact water mass transports than perturbations in downstream sources like the Arctic freshwater outflow. These boundary current dynamics are directly related to the qualitative stratification in transition zones and illustrate the interaction between the alpha and beta oceans.

Beschreibung: This research was supported by the Research Council of Norway project NORTH. Support for the publication was provided by the University of Bergen. Ocean Outlook has supported a research visit for EL to Woods Hole Oceanographic Institute where much of the current work has been carried out. Support forMAS was provided by the National Science Foundation Grant OCE-1558742.

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 883-904, doi:10.1175/JPO-D-17-0084.1.

Beschreibung: The dynamics controlling the along-valley (cross shelf) flow in idealized shallow shelf valleys with small to moderate Burger number are investigated, and analytical scales of the along-valley flows are derived. This paper follows Part I, which shows that along-shelf winds in the opposite direction to coastal-trapped wave propagation (upwelling regime) force a strong up-valley flow caused by the formation of a lee wave. In contrast, along-shelf winds in the other direction (downwelling regime) do not generate a lee wave and consequently force a relatively weak net down-valley flow. The valley flows in both regimes are cyclostrophic with 0(1) Rossby number. A major difference between the two regimes is the along-shelf length scales of the along-valley flows L. In the upwelling regime Ls, depends on the valley width W, and the wavelength lambda(1w) of the coastal-trapped lee wave arrested by the along-shelf flow U-s. In the downwelling regime L depends on the inertial length scale U-s|'f and W-c. The along-valley velocity scale in the upwelling regime, given by V-u approximate to root pi H-c/H-s integral W-c lambda(1w)/2 pi L-x (1+L-x(2)/L-c(2))(-1) e(-(pi Wc)/(lambda 1w),) is based on potential vorticity (PV) conservation and lee-wave dynamics (Hs and H, are the shelf and valley depth scales, respectively, and fis the Coriolis parameter). The velocity scale in the downwelling regime, given by |v(d)| approximate to (H-s/H-s)[1 + (L-x(2)/L-x(2))](-1) fL, is based on PV conservation. The velocity scales are validated by the numerical sensitivity simulations and can be useful for observational studies of along -valley transports. The work provides a framework for investigating cross -shelf transport induced by irregular shelf bathymetry and calls for future studies of this type under realistic environmental conditions and over a broader parameter space.

Beschreibung: Both WGZ and SJL were supported by the National Science Foundation (NSF) through Grant OCE 1154575.WGZis also supported by the NSF Grant OCE 1634965 and SJL by NSF Grant OCE 1558874.

Beschreibung: 2018-10-16

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 2799-2827, doi:10.1175/JPO-D-18-0057.1.

Beschreibung: The fjords that connect Greenland’s glaciers to the ocean are gateways for importing heat to melt ice and for exporting meltwater into the ocean. The transport of heat and meltwater can be modulated by various drivers of fjord circulation, including freshwater, local winds, and shelf variability. Shelf-forced flows (also known as the intermediary circulation) are the dominant mode of variability in two major fjords of east Greenland, but we lack a dynamical understanding of the fjord’s response to shelf forcing. Building on observations from east Greenland, we use numerical simulations and analytical models to explore the dynamics of shelf-driven flows. For the parameter space of Greenlandic fjords, we find that the fjord’s response is primarily a function of three nondimensional parameters: the fjord width over the deformation radius (W/Rd), the forcing time scale over the fjord adjustment time scale, and the forcing amplitude (shelf pycnocline displacements) over the upper-layer thickness. The shelf-forced flows in both the numerical simulations and the observations can largely be explained by a simple analytical model for Kelvin waves propagating around the fjord. For fjords with W/Rd 〉 0.5 (most Greenlandic fjords), 3D dynamics are integral to understanding shelf forcing—the fjord dynamics cannot be approximated with 2D models that neglect cross-fjord structure. The volume flux exchanged between the fjord and shelf increases for narrow fjords and peaks around the resonant forcing frequency, dropping off significantly at higher- and lower-frequency forcing.

Beschreibung: This work was funded by NSF Grant OCE-1536856 and by the NOAA Climate and Global Change Postdoctoral Fellowship.

Beschreibung: © The Author(s), 2018. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Conservation Physiology 6 (2018): coy049, doi:10.1093/conphys/coy049.

Beschreibung: Male baleen whales have long been suspected to have annual cycles in testosterone, but due to difficulty in collecting endocrine samples, little direct evidence exists to confirm this hypothesis. Potential influences of stress or adrenal stress hormones (cortisol, corticosterone) on male reproduction have also been difficult to study. Baleen has recently been shown to accumulate steroid hormones during growth, such that a single baleen plate contains a continuous, multi-year retrospective record of the whale’s endocrine history. As a preliminary investigation into potential testosterone cyclicity in male whales and influences of stress, we determined patterns in immunoreactive testosterone, two glucocorticoids (cortisol and corticosterone), and stable-isotope (SI) ratios, across the full length of baleen plates from a bowhead whale (Balaena mysticetus), a North Atlantic right whale (Eubalaena glacialis) and a blue whale (Balaenoptera musculus), all adult males. Baleen was subsampled at 2 cm (bowhead, right) or 1 cm (blue) intervals and hormones were extracted from baleen powder with methanol, followed by quantification of all three hormones using enzyme immunoassays validated for baleen extract of these species. Baleen of all three males contained regularly spaced peaks in testosterone content, with number and spacing of testosterone peaks corresponding well to SI data and to species-specific estimates of annual baleen growth rate. Cortisol and corticosterone exhibited some peaks that co-occurred with testosterone peaks, while other glucocorticoid peaks occurred independent of testosterone peaks. The right whale had unusually high glucocorticoids during a period with a known entanglement in fishing gear and a possible disease episode; in the subsequent year, testosterone was unusually low. Further study of baleen testosterone patterns in male whales could help clarify conservation- and management-related questions such as age of sexual maturity, location and season of breeding, and the potential effect of anthropogenic and natural stressors on male testosterone cycles.

Beschreibung: This work was supported by (1) the Arizona Board of Regents Technology Research Initiative Fund; (2) the Center for Bioengineering Innovation at Northern Arizona University; (3) the Greenland Institute of Natural Resources; (4) the Woods Hole Oceanographic Institution Ocean Life Institute and (5) Fisheries and Ocean Canada’s (DFO) Priorities and Partnership Strategic Initiatives Fund and Oceans Protection Plan.

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 31 (2018): 4847-4863, doi:10.1175/JCLI-D-17-0802.1.

Beschreibung: The sensitivity of sea ice to the temperature of inflowing Atlantic water across the Greenland–Scotland Ridge is investigated using an eddy-resolving configuration of the Massachusetts Institute of Technology General Circulation Model with idealized topography. During the last glacial period, when climate on Greenland is known to have been extremely unstable, sea ice is thought to have covered the Nordic seas. The dramatic excursions in climate during this period, seen as large abrupt warming events on Greenland and known as Dansgaard–Oeschger (DO) events, are proposed to have been caused by a rapid retreat of Nordic seas sea ice. Here, we show that a full sea ice cover and Arctic-like stratification can exist in the Nordic seas given a sufficiently cold Atlantic inflow and corresponding low transport of heat across the Greenland–Scotland Ridge. Once sea ice is established, continued sea ice formation and melt efficiently freshens the surface ocean and makes the deeper layers more saline. This creates a strong salinity stratification in the Nordic seas, similar to today’s Arctic Ocean, with a cold fresh surface layer protecting the overlying sea ice from the warm Atlantic water below. There is a nonlinear response in Nordic seas sea ice to Atlantic water temperature with simulated large abrupt changes in sea ice given small changes in inflowing temperature. This suggests that the DO events were more likely to have occurred during periods of reduced warm Atlantic water inflow to the Nordic seas.

Beschreibung: The research was supported by the Centre for Climate Dynamics at the Bjerknes Centre for Climate Research. The research leading to these results is part of the ice2ice project funded by the European Research Council under the European Community Seventh Framework Programme (FP7/2007-2013)/ERC Grant Agreement 610055.

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 773-794, doi:10.1175/JPO-D-17-0205.1.

Beschreibung: Fourteen autonomous profiling floats, equipped with CTDs, were deployed in the deep eastern and western basins of the Gulf of Mexico over a four-year interval (July 2011–August 2015), producing a total of 706 casts. This is the first time since the early 1970s that there has been a comprehensive survey of water masses in the deep basins of the Gulf, with better vertical resolution than available from older ship-based surveys. Seven floats had 14-day cycles with parking depths of 1500 m, and the other half from the U.S. Argo program had varying cycle times. Maps of characteristic water masses, including Subtropical Underwater, Antarctic Intermediate Water (AAIW), and North Atlantic Deep Water, showed gradients from east to west, consistent with their sources being within the Loop Current (LC) and the Yucatan Channel waters. Altimeter SSH was used to characterize profiles being in LC or LC eddy water or in cold eddies. The two-layer nature of the deep Gulf shows isotherms being deeper in the warm anticyclonic LC and LC eddies and shallower in the cold cyclones. Mixed layer depths have an average seasonal signal that shows maximum depths (~60 m) in January and a minimum in June–July (~20 m). Basin-mean steric heights from 0–50-m dynamic heights and altimeter SSH show a seasonal range of ~12 cm, with significant interannual variability. The translation of LC eddies across the western basin produces a region of low homogeneous potential vorticity centered over the deepest part of the western basin.

Beschreibung: The authors were supported by the Department of the Interior, Bureau of Ocean Energy Management (BOEM), Contract M08PC20043 to Leidos, Inc., Raleigh, North Carolina.

Beschreibung: 2018-10-04

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 905-923, doi:10.1175/JPO-D-17-0133.1.

Beschreibung: Observations of turbulent kinetic energy, dissipation, and turbulent stress were collected in the middle reaches of Chesapeake Bay and were used to assess second-moment closure predictions of turbulence generated beneath breaking waves. Dissipation scaling indicates that the turbulent flow structure observed during a 10-day wind event was dominated by a three-layer response that consisted of 1) a wave transport layer, 2) a surface log layer, and 3) a tidal, bottom boundary layer limited by stable stratification. Below the wave transport layer, turbulent mixing was limited by stable stratification. Within the wave transport layer, where dissipation was balanced by a divergence in the vertical turbulent kinetic energy flux, the eddy viscosity was significantly underestimated by second-moment turbulence closure models, suggesting that breaking waves homogenized the mixed surface layer to a greater extent than the simple model of TKE diffusing away from a source at the surface. While the turbulent transport of TKE occurred largely downgradient, the intermittent downward sweeps of momentum generated by breaking waves occurred largely independent of the mean shear. The underprediction of stress in the wave transport layer by second-moment closures was likely due to the inability of the eddy viscosity model to capture the nonlocal turbulent transport of the momentum flux beneath breaking waves. Finally, the authors hypothesize that large-scale coherent turbulent eddies played a significant role in transporting momentum generated near the surface to depth.

Beschreibung: This work was supported by National Science Foundation Grants OCE-1061609 and OCE-1339032.

Beschreibung: 2018-10-19

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Atmospheric and Oceanic Technology 35 (2018): 281-297, doi:10.1175/JTECH-D-17-0076.1.

Beschreibung: The wavenumber spectrum of sea surface height (SSH) is an important indicator of the dynamics of the ocean interior. While the SSH wavenumber spectrum has been well studied at mesoscale wavelengths and longer, using both in situ oceanographic measurements and satellite altimetry, it remains largely unknown for wavelengths less than ~70 km. The Surface Water Ocean Topography (SWOT) satellite mission aims to resolve the SSH wavenumber spectrum at 15–150-km wavelengths, which is specified as one of the mission requirements. The mission calibration and validation (CalVal) requires the ground truth of a synoptic SSH field to resolve the targeted wavelengths, but no existing observational network is able to fulfill the task. A high-resolution global ocean simulation is used to conduct an observing system simulation experiment (OSSE) to identify the suitable oceanographic in situ measurements for SWOT SSH CalVal. After fixing 20 measuring locations (the minimum number for resolving 15–150-km wavelengths) along the SWOT swath, four instrument platforms were tested: pressure-sensor-equipped inverted echo sounders (PIES), underway conductivity–temperature–depth (UCTD) sensors, instrumented moorings, and underwater gliders. In the context of the OSSE, PIES was found to be an unsuitable tool for the target region and for SSH scales 15–70 km; the slowness of a single UCTD leads to significant aliasing by high-frequency motions at short wavelengths below ~30 km; an array of station-keeping gliders may meet the requirement; and an array of moorings is the most effective system among the four tested instruments for meeting the mission’s requirement. The results shown here warrant a prelaunch field campaign to further test the performance of station-keeping gliders.

Beschreibung: The authors would like to acknowledge the funding sources: the SWOT mission (JW, LF, DM); NASA Projects NNX13AE32G, NNX16AH76G, and NNX17AH54G (TF); and NNX16AH66G and NNX17AH33G (BQ). AF and MF were funded by the Keck Institute for Space Studies (which is generously supported by the W. M. Keck Foundation) through the project Science-driven Autonomous and Heterogeneous Robotic Networks: A Vision for Future Ocean Observations (http://kiss.caltech.edu/?techdev/seafloor/seafloor.html).

Beschreibung: 2018-08-07

Beschreibung: © The Author(s), 2018. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Environmental Epigenetics 4 (2018): dvy005, doi:10.1093/eep/dvy005.

Beschreibung: There is growing evidence that environmental toxicants can affect various physiological processes by altering DNA methylation patterns. However, very little is known about the impact of toxicant-induced DNA methylation changes on gene expression patterns. The objective of this study was to determine the genome-wide changes in DNA methylation concomitant with altered gene expression patterns in response to 3, 3’, 4, 4’, 5-pentachlorobiphenyl (PCB126) exposure. We used PCB126 as a model environmental chemical because the mechanism of action is well-characterized, involving activation of aryl hydrocarbon receptor, a ligand-activated transcription factor. Adult zebrafish were exposed to 10 nM PCB126 for 24 h (water-borne exposure) and brain and liver tissues were sampled at 7 days post-exposure in order to capture both primary and secondary changes in DNA methylation and gene expression. We used enhanced Reduced Representation Bisulfite Sequencing and RNAseq to quantify DNA methylation and gene expression, respectively. Enhanced reduced representation bisulfite sequencing analysis revealed 573 and 481 differentially methylated regions in the liver and brain, respectively. Most of the differentially methylated regions are located more than 10 kilobases upstream of transcriptional start sites of the nearest neighboring genes. Gene Ontology analysis of these genes showed that they belong to diverse physiological pathways including development, metabolic processes and regeneration. RNAseq results revealed differential expression of genes related to xenobiotic metabolism, oxidative stress and energy metabolism in response to polychlorinated biphenyl exposure. There was very little correlation between differentially methylated regions and differentially expressed genes suggesting that the relationship between methylation and gene expression is dynamic and complex, involving multiple layers of regulation.

Beschreibung: This work was supported by the National Institute of Health Outstanding New Environmental Scientist Award to NA (NIH R01ES024915) and Woods Hole Center for Oceans and Human Health [National Institutes of Health (NIH) grant P01ES021923 and National Science Foundation Grant OCE-1314642 to M. Hahn, J. Stegeman, NA and SK].

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 163-174, doi:10.1175/JPO-D-17-0161.1.

Beschreibung: The general problem of exchange from a shallow shelf across sharp topography to the deep ocean forced by narrow, cross-shelf wind jets is studied using quasigeostrophic theory and an idealized primitive equation numerical model. Interest is motivated by katabatic winds that emanate from narrow fjords in southeast Greenland, although similar topographically constrained wind jets are found throughout the world’s oceans. Because there is no net vorticity input by the wind, the circulation is largely confined to the region near the forcing. Circulation over the shelf is limited by bottom friction for weakly stratified flows, but stratification allows for much stronger upper-layer flows that are regulated by weak coupling to the lower layer. Over the sloping topography, the topographic beta effect limits the deep flow, while, for sufficient stratification, the upper-layer flow can cross the topography to connect the shelf to the open ocean. This can be an effective transport mechanism even for short, strong wind events because damping of the upper-layer flow is weak. A variety of transients are generated for an abrupt onset of winds, including short topography Rossby waves, long topographic Rossby waves, and inertial waves. Using parameters representative of southeast Greenland, katabatic wind events will force an offshore transport of O(0.4) Sv (1 Sv ≡ 106 m3 s−1) that, when considered for 2 days, will result in an offshore flux of O(5 × 1010) m3.

Beschreibung: MAS was supported by the National Science Foundation under Grant OCE-1533170.

Beschreibung: 2018-07-18

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 1533-1541, doi:10.1175/JPO-D-17-0267.1.

Beschreibung: Our study analyzes measurements primarily from two Floating Instrument Platform (FLIP) field programs and from the Air–Sea Interaction Tower (ASIT) site to examine the relationship between the wind and sea surface stress for contrasting conditions. The direct relationship of the surface momentum flux to U2 is found to be better posed than the relationship between and U, where U is the wind speed and is the friction velocity. Our datasets indicate that the stress magnitude often decreases significantly with height near the surface due to thin marine boundary layers and/or enhanced stress divergence close to the sea surface. Our study attempts to correct the surface stress estimated from traditional observational levels by using multiple observational levels near the surface and extrapolating to the surface. The effect of stability on the surface stress appears to be generally smaller than errors due to the stress divergence. Definite conclusions require more extensive measurements close to the sea surface.

Beschreibung: This work was supported by the U.S. Office of Naval Research through Award N00014-16-1-2600. We

Beschreibung: 2019-01-10

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 1815-1830, doi:10.1175/JPO-D-17-0275.1.

Beschreibung: Recent progress in direct numerical simulations (DNSs) of stratified turbulent flows has led to increasing attention to the validity of the constancy of the dissipation flux coefficient Γ in the Osborn’s eddy diffusivity model. Motivated by lack of observational estimates of Γ, particularly under weakly stratified deep-ocean conditions, this study estimates Γ using deep microstructure profiles collected in various regions of the North Pacific and Southern Oceans. It is shown that Γ is not constant but varies significantly with the Ozmidov/Thorpe scale ratio ROT in a fashion similar to that obtained by previous DNS studies. Efficient mixing events with Γ ~ O(1) and ROT ~ O(0.1) tend to be frequently observed in the deep ocean (i.e., weak stratification), while moderate mixing events with Γ ~ O(0.1) and ROT ~ O(1) tend to be observed in the upper ocean (i.e., strong stratification). The observed negative relationship between Γ and ROT is consistent with a simple scaling that can be derived from classical turbulence theories. In contrast, the observed results exhibit no definite relationships between Γ and the buoyancy Reynolds number Reb, although Reb has long been thought to be another key parameter that controls Γ.

Beschreibung: This study was supported by MEXT KAKENHI Grant JP15H05824 and JSPS KAKENHI Grant JP15H02131.

Beschreibung: 2019-02-15

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 31 (2018): 7751-7769, doi:10.1175/JCLI-D-18-0184.1.

Beschreibung: Decadal variability of the subsurface ocean heat content (OHC) in the Indian Ocean is investigated using a coupled climate model experiment, in which observed eastern tropical Pacific sea surface temperature (EPSST) anomalies are specified. This study intends to understand the contributions of external forcing relative to those of internal variability associated with EPSST, as well as the mechanisms by which the Pacific impacts Indian Ocean OHC. Internally generated variations associated with EPSST dominate decadal variations in the subsurface Indian Ocean. Consistent with ocean reanalyses, the coupled model reproduces a pronounced east–west dipole structure in the southern tropical Indian Ocean and discontinuities in westward-propagating signals in the central Indian Ocean around 100°E. This implies distinct mechanisms by which the Pacific impacts the eastern and western Indian Ocean on decadal time scales. Decadal variations of OHC in the eastern Indian Ocean are attributed to 1) western Pacific surface wind anomalies, which trigger oceanic Rossby waves propagating westward through the Indonesian Seas and influence Indonesian Throughflow transport, and 2) zonal wind anomalies over the central tropical Indian Ocean, which trigger eastward-propagating Kelvin waves. Decadal variations of OHC in the western Indian Ocean are linked to conditions in the Pacific via changes in the atmospheric Walker cell, which trigger anomalous wind stress curl and Ekman pumping in the central tropical Indian Ocean. Westward-propagating oceanic Rossby waves extend the influence of this anomalous Ekman pumping to the western Indian Ocean.

Beschreibung: This research was supported by the Independent Research and Development Program at WHOI to CCU, an NSF OCE PO grant (NSF OCE- 1242989) to Young-Oh Kwon, NOAA CP CVP grants (NA15OAR4310176 and NA17OAR4310255) to Hyodae Seo, and a research grant fromtheMinistry of Science and Technology of the People’s Republic of China to Tsinghua University (2017YFA0603902).

Beschreibung: 2019-02-13

Beschreibung: Author Posting. © The Authors, 2018. This article is posted here by permission of The Royal Astronomical Society for personal use, not for redistribution. The definitive version was published in Geophysical Journal International 215 (2018): 1072–1087, doi:10.1093/gji/ggy203.

Beschreibung: An earthquake rupture process can be kinematically described by rupture velocity, duration and spatial extent. These key kinematic source parameters provide important constraints on earthquake physics and rupture dynamics. In particular, core questions in earthquake science can be addressed once these properties of small earthquakes are well resolved. However, these parameters of small earthquakes are poorly understood, often limited by available data sets and methodologies. The Incorporated Research Institutions for Seismology Community Wavefield Experiment in Oklahoma deployed ∼350 three-component nodal stations within 40 km2 for a month, offering an unprecedented opportunity to test new methodologies for resolving small earthquake finite source properties in high resolution. In this study, we demonstrate the power of the nodal data set to resolve the variations in the seismic wavefield over the focal sphere due to the finite source attributes of an M2 earthquake within the array. The dense coverage allows us to tightly constrain rupture area using the second moment method even for such a small earthquake. The M2 earthquake was a strike-slip event and unilaterally propagated towards the surface at 90 per cent local S-wave speed (2.93 km s−1). The earthquake lasted ∼0.019 s and ruptured Lc ∼70 m and Wc ∼45 m. With the resolved rupture area, the stress-drop of the earthquake is estimated as 7.3 MPa for Mw 2.3. We demonstrate that the maximum and minimum bounds on rupture area are within a factor of two, much lower than typical stress-drop uncertainty, despite a suboptimal station distribution. The rupture properties suggest that there is little difference between the M2 Oklahoma earthquake and typical large earthquakes. The new three-component nodal systems have great potential for improving the resolution of studies of earthquake source properties.

Beschreibung: WF is currently supported by the Postdoctoral Scholar Program at the Woods Hole Oceanographic Institution, with funding provided by the Weston Howland Jr. Postdoctoral Scholarship. JM was partially supported by SCEC grant #17177 at Woods Hole Oceanographic Institution. This research was supported by the Southern California Earthquake Center (Contribution No. 8014). SCEC is funded by NSF Cooperative Agreement EAR-1033462 and USGS Cooperative Agreement G12AC20038.

Beschreibung: Author Posting. © The Authors, 2018. This article is posted here by permission of The Royal Astronomical Society for personal use, not for redistribution. The definitive version was published in Geophysical Journal International 215 (2018): 942–958, doi:10.1093/gji/ggy316.

Beschreibung: Surface waves recorded by global arrays have proven useful for locating tectonic earthquakes and in detecting slip events depleted in high frequency, such as glacial quakes. We develop a novel method using an aggregation of small- to continental-scale arrays to detect and locate seismic sources with Rayleigh waves at 20–50 s period. The proposed method is a hybrid approach including first dividing a large aperture aggregate array into Delaunay triangular subarrays for beamforming, and then using the resolved surface wave propagation directions and arrival times from the subarrays as data to formulate an inverse problem to locate the seismic sources and their origin times. The approach harnesses surface wave coherence and maximizes resolution of detections by combining measurements from stations spanning the whole U.S. continent. We tested the method with earthquakes, glacial quakes and landslides. The results show that the method can effectively resolve earthquakes as small as ∼M3 and exotic slip events in Greenland. We find that the resolution of the locations is non-uniform with respect to azimuth, and decays with increasing distance between the source and the array when no calibration events are available. The approach has a few advantages: the method is insensitive to seismic event type, it does not require a velocity model to locate seismic sources, and it is computationally efficient. The method can be adapted to real-time applications and can help in identifying new classes of seismic sources.

Beschreibung: WF is currently supported by the Postdoctoral Scholar Program at the Woods Hole Oceanographic Institution, with funding provided by the Weston Howland Jr. Postdoctoral Scholarship. This work was supported by National Science Foundation grant EAR-1358520 at Scripps Institution of Oceanography, UC San Diego.

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018):1941-1950, doi:10.1175/JPO-D-17-0194.1.

Beschreibung: Subglacial discharges have been observed to generate buoyant plumes along the ice face of Greenland tidewater glaciers. These plumes have been traditionally modeled using classical plume theory, and their characteristic parameters (e.g., velocity) are employed in the widely used three-equation melt parameterization. However, the applicability of plume theory for three-dimensional turbulent wall plumes is questionable because of the complex near-wall plume dynamics. In this study, corrections to the classical plume theory are introduced to account for the presence of a wall. In particular, the drag and entrainment coefficients are quantified for a three-dimensional turbulent wall plume using data from direct numerical simulations. The drag coefficient is found to be an order of magnitude larger than that for a boundary layer flow over a flat plate at a similar Reynolds number. This result suggests a significant increase in the melting estimates by the current parameterization. However, the volume flux in a wall plume is found to be one-half that of a conical plume that has 2 times the buoyancy flux. This finding suggests that the total entrainment (per unit area) of ambient water is the same and that the plume scalar characteristics (i.e., temperature and salinity) can be predicted reasonably well using classical plume theory.

Beschreibung: This work was supported by the Linné FLOW Centre at KTH and the Academy of Finland Center of Excellence Programme Grant 307331 (author Ezhova) and by VR Swedish Research Council GrantVR2014-5001 (author Brandt). Support to author Cenedese was given by NSF Project OCE-1434041.

Beschreibung: 2019-02-27

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 511-529, doi:10.1175/JPO-D-17-0140.1.

Beschreibung: The large-scale circulation of the bottom layer of the Gulf of Mexico is analyzed, with special attention to the historically least studied western basin. The analysis is based on 4 years of data collected by 158 subsurface floats parked at 1500 and 2500 m and is complemented with data collected by current meter moorings in the western basin during the same period. Three main circulation patterns stand out: a cyclonic boundary current, a cyclonic gyre in the abyssal plain, and the very high eddy kinetic energy observed in the eastern Gulf. The boundary current and the cyclonic gyre appear as distinct features, which interact in the western tip of the Yucatan shelf. The persistence and continuity of the boundary current is addressed. Although high variability is observed, the boundary flow serves as a pathway for water to travel around the western basin in approximately 2 years. An interesting discovery is the separation of the boundary current over the northwestern slope of the Yucatan shelf. The separation and retroflection of the along-slope current appears to be a persistent feature and is associated with anticyclonic eddies whose genesis mechanism remains to be understood. As the boundary flow separates, it feeds into the westward flow of the deep cyclonic gyre. The location of this gyre—named the Sigsbee Abyssal Gyre—coincides with closed geostrophic contours, so eddy–topography interaction via bottom form stresses may drive this mean flow. The contribution to the cyclonic vorticity of the gyre by modons traveling under Loop Current eddies is discussed.

Beschreibung: This work was supported by the Bureau of Ocean Energy Management (BOEM) under Contract M10PC00112 assigned to Leidos, Inc.

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 31 (2018): 4157-4174, doi:10.1175/JCLI-D-17-0654.1.

Beschreibung: Decadal variabilities in Indian Ocean subsurface ocean heat content (OHC; 50–300 m) since the 1950s are examined using ocean reanalyses. This study elaborates on how Pacific variability modulates the Indian Ocean on decadal time scales through both oceanic and atmospheric pathways. High correlations between OHC and thermocline depth variations across the entire Indian Ocean Basin suggest that OHC variability is primarily driven by thermocline fluctuations. The spatial pattern of the leading mode of decadal Indian Ocean OHC variability closely matches the regression pattern of OHC on the interdecadal Pacific oscillation (IPO), emphasizing the role of the Pacific Ocean in determining Indian Ocean OHC decadal variability. Further analyses identify different mechanisms by which the Pacific influences the eastern and western Indian Ocean. IPO-related anomalies from the Pacific propagate mainly through oceanic pathways in the Maritime Continent to impact the eastern Indian Ocean. By contrast, in the western Indian Ocean, the IPO induces wind-driven Ekman pumping in the central Indian Ocean via the atmospheric bridge, which in turn modifies conditions in the southwestern Indian Ocean via westward-propagating Rossby waves. To confirm this, a linear Rossby wave model is forced with wind stresses and eastern boundary conditions based on reanalyses. This linear model skillfully reproduces observed sea surface height anomalies and highlights both the oceanic connection in the eastern Indian Ocean and the role of wind-driven Ekman pumping in the west. These findings are also reproduced by OGCM hindcast experiments forced by interannual atmospheric boundary conditions applied only over the Pacific and Indian Oceans, respectively.

Beschreibung: This research was supported by a scholarship from the China Scholarship Council (CSC) to X. J., a research fellowship by the Alexander von Humboldt Foundation to C. C. U., an NSF OCE PO Grant (OCE- 1242989) to Y.-O. K., the ONR Young Investigator Award (N00014-15-1-2588) to H. S., and a research grant from the Ministry of Science and Technology of the People’s Republic of China to Tsinghua University (2017YFA0603902).

Beschreibung: 2018-10-30

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 31 (2018): 2771-2796, doi:10.1175/JCLI-D-17-0061.1.

Beschreibung: The Generalized Equilibrium Feedback Analysis (GEFA) is used to distinguish the influence of the Oyashio Extension (OE) and the Kuroshio Extension (KE) variability on the atmosphere from 1979 to 2014 from that of the main SST variability modes, using seasonal mean anomalies. Remote SST anomalies are associated with each single oceanic regressor, but the multivariate approach efficiently confines their SST footprints. In autumn [October–December (OND)], the OE meridional shifts are followed by a North Pacific Oscillation (NPO)-like signal. The OE influence is not investigated in winter [December–February (DJF)] because of multicollinearity, but a robust response with a strong signal over the Bering Sea is found in late winter/early spring [February–April (FMA)], a northeastward strengthening of the Aleutian low following a northward OE shift. A robust response to the KE variability is found in autumn, but not in winter and late winter when the KE SST footprint becomes increasingly small and noisy as regressors are added in GEFA. In autumn, a positive PDO is followed by a northward strengthening of the Aleutian low and a southward shift of the storm track in the central Pacific, reflecting the surface heat flux footprint in the central Pacific. In winter, the PDO shifts the maximum baroclinicity and storm track southward, the response strongly tilts westward with height in the North Pacific, and there is a negative NAO-like teleconnection. In late winter, the North Pacific NPO-like response to the PDO interferes negatively with the response to the OE and is only detected when the OE is represented in GEFA. A different PDO influence on the atmospheric circulation is found from 1958 to 1977.

Beschreibung: This research has received funding from the European Union 7th Framework Program (FP7 2007-2013) under Grant Agreement 308299 (NACLIM) and from NSF Grants AGS CLD 1035423 and OCE PO 1242989.

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 1367-1373, doi:10.1175/JPO-D-17-0185.1.

Beschreibung: An earlier study indicates that the side melting of icebergs subject to vertically homogeneous horizontal velocities is controlled by two distinct regimes, which depend on the melt plume behavior and produce a nonlinear dependence of side melt rate on velocity. Here, we extend this study to consider ice blocks melting in a two-layer vertically sheared flow in a laboratory setting. It is found that the use of the vertically averaged flow speed in current melt parameterizations gives an underestimate of the submarine side melt rate, in part because of the nonlinearity of the dependence of the side melt rate on flow speed but also because vertical shear in the horizontal velocity profile fundamentally changes the flow splitting around the ice block and consequently the velocity felt by the ice surface. An observational record of 90 icebergs in a Greenland fjord suggests that this effect could produce an average underestimate of iceberg side melt rates of 21%.

Beschreibung: A. F. was supported by NA14OAR4320106 from the National Oceanic and Atmospheric Administration, U.S. Department of Commerce. C. C. was supported by NSF OCE-1658079 and F. S. was supported by NSF OCE-1657601 and NSF PLR-1743693.

Beschreibung: 2018-12-12

Beschreibung: Author Posting. © The Author(s), 2018. This article is posted here by permission of Oxford University Press for personal use, not for redistribution. The definitive version was published in Geophysical Journal International 214 (2018): 2224–2235, doi:10.1093/gji/ggy201.

Beschreibung: The key kinematic earthquake source parameters: rupture velocity, duration and area, shed light on earthquake dynamics, provide direct constraints on stress drop, and have implications for seismic hazard. However, for moderate and small earthquakes, these parameters are usually poorly constrained due to limitations of the standard analysis methods. Numerical experiments by Kaneko and Shearer demonstrated that standard spectral fitting techniques can lead to roughly one order of magnitude variation in stress-drop estimates that do not reflect the actual rupture properties even for simple crack models. We utilize these models to explore an alternative approach where we estimate the rupture area directly. For the suite of models, the area averaged static stress drop is nearly constant for models with the same underlying friction law, yet corner-frequency-based stress-drop estimates vary by a factor of 5–10 even for noise-free data. Alternatively, we simulated inversions for the rupture area as parametrized by the second moments of the slip distribution. A natural estimate for the rupture area derived from the second moments is A = πLcWc, where Lc and Wc are the characteristic rupture length and width. This definition yields estimates of stress drop that vary by only 10 per cent between the models but are slightly larger than the true area averaged values. We simulate inversions for the second moments for the various models and find that the area can be estimated well when there are at least 15 available measurements of apparent duration at a variety of take-off angles. The improvement compared to azimuthally averaged corner-frequency-based approaches results from the second moments accounting for directivity and removing the assumption of a circular rupture area, both of which bias the standard approach. We also develop a new method that determines the minimum and maximum values of rupture area that are consistent with a particular data set at the 95 per cent confidence level. For the Kaneko and Shearer models with 20+ randomly distributed observations and ∼10 per cent noise levels, we find that the maximum and minimum bounds on rupture area typically vary by a factor of two and that the minimum stress drop is often more tightly constrained than the maximum.

Beschreibung: This work was supported by USGS NEHRP Award G17AP00029. The research was supported by the Southern California Earthquake Center (SCEC; Contribution No. 8013). SCEC is funded by NSF Cooperative Agreement EAR-1033462 and USGS Cooperative Agreement G12AC20038. YK was supported by both public funding from the Government of New Zealand and the Royal Society of New Zealand’s Rutherford Discovery Fellowship.

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 31 (2018): 8059-8079, doi:10.1175/JCLI-D-17-0769.1.

Beschreibung: We use the method of least squares with Lagrange multipliers to fit an ocean general circulation model to the Multiproxy Approach for the Reconstruction of the Glacial Ocean Surface (MARGO) estimate of near sea surface temperature (NSST) at the Last Glacial Maximum (LGM; circa 23–19 thousand years ago). Compared to a modern simulation, the resulting global, last-glacial ocean state estimate, which fits the MARGO data within uncertainties in a free-running coupled ocean–sea ice simulation, has global-mean NSSTs that are 2°C lower and greater sea ice extent in all seasons in both the Northern and Southern Hemispheres. Increased brine rejection by sea ice formation in the Southern Ocean contributes to a stronger abyssal stratification set principally by salinity, qualitatively consistent with pore fluid measurements. The upper cell of the glacial Atlantic overturning circulation is deeper and stronger. Dye release experiments show similar distributions of Southern Ocean source waters in the glacial and modern western Atlantic, suggesting that LGM NSST data do not require a major reorganization of abyssal water masses. Outstanding challenges in reconstructing LGM ocean conditions include reducing effects from model biases and finding computationally efficient ways to incorporate abyssal tracers in global circulation inversions. Progress will be aided by the development of coupled ocean–atmosphere–ice inverse models, by improving high-latitude model processes that connect the upper and abyssal oceans, and by the collection of additional paleoclimate observations.

Beschreibung: DEA was supported by a NSF Graduate Research Fellowship and NSF Grant OCE-1060735. OM acknowledges support from the NSF. GF was supported by NASA Award 1553749 and Simons Foundation Award 549931.

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 1969-1993, doi:10.1175/JPO-D-18-0031.1.

Beschreibung: Upstream mean semidiurnal internal tidal energy flux has been found in the Gulf Stream in hydrodynamical model simulations of the Atlantic Ocean. A major source of the energy in the simulations is the south edge of Georges Bank, where strong and resonant Gulf of Maine tidal currents are found. An explanation of the flux pattern within the Gulf Stream is that internal wave modal rays can be strongly redirected by baroclinic currents and even trapped (ducted) by current jets that feature strong velocities above the thermocline that are directed counter to the modal wavenumber vector (i.e., when the waves travel upstream). This ducting behavior is analyzed and explained here with ray-based wave propagation studies for internal wave modes with anisotropic wavenumbers, as occur in mesoscale background flow fields. Two primary analysis tools are introduced and then used to analyze the strong refraction and ducting: the generalized Jones equation governing modal properties and ray equations that are suitable for studying waves with anisotropic wavenumbers.

Beschreibung: The Woods Hole research was supported by National Science Foundation Grant OCE-1060430 and by the Office of Naval Research Grants N00014-11-1-0701 and N00014-17-1-2624. The USM research was supported by ONR Grant N00014-15-1-2288 and National Science Foundation Grant OCE-1537449.

Beschreibung: 2019-02-28

Beschreibung: Author Posting. © The Authors, 2018. This article is posted here by permission of The Royal Astronomical Society for personal use, not for redistribution. The definitive version was published in Geophysical Journal International 215 (2018): 713–735, doi:10.1093/gji/ggy313.

Beschreibung: Gas flux in volcanic conduits is often associated with long-period oscillations known as seismic tremor (Lesage et al.; Nadeau et al.). In this study, we revisit and extend the ‘magma wagging’and ‘whirling’models for seismic tremor, in order to explore the effects of gas flux on the motion of a magma column surrounded by a permeable vesicular annulus (Jellinek & Bercovici; Bercovici et al.; Liao et al.). We find that gas flux flowing through the annulus leads to a Bernoulli effect, which causes waves on the magma column to become unstable and grow. Specifically, the Bernoulli effects are associated with torques and forces acting on the magma column, increasing its angular momentum and energy. As the displacement of the magma column becomes large due to the Bernoulli effect, frictional drag on the conduit wall decelerates the motions of the column, restoring them to small amplitude. Together, the Bernoulli effect and the damping effect contribute to a self-sustained wagging-and-whirling mechanism that help explain the longevity of long-period seismic tremor.

Beschreibung: This work was supported by National Science Foundation grants EAR-1344538 and EAR-1645057

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 1375-1384, doi:10.1175/JPO-D-17-0266.1.

Beschreibung: The relationship between net mixing and the estuarine exchange flow may be quantified using a salinity variance budget. Here “mixing” is defined as the rate of destruction of volume-integrated salinity variance, and the exchange flow is quantified using the total exchange flow. These concepts are explored using an idealized 3D model estuary. It is shown that in steady state (e.g., averaging over the spring–neap cycle) the volume-integrated mixing is approximately given by Mixing ≅ SinSoutQr, where Sin and Sout are the representative salinities of in- and outflowing layers at the mouth and Qr is the river volume flux. This relationship provides an extension of the familiar Knudsen relation, in which the exchange flow is diagnosed based on knowledge of these same three quantities, quantitatively linking mixing to the exchange flow.

Beschreibung: The work was supported by the National Science Foundation through Grants OCE-1736242 to PM and OCE-1736539 to WRG and by the German Research Foundation through Grants TRR 181 and GRK 2000 to HB.

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 2703-2719, doi:10.1175/JPO-D-17-0245.1.

Beschreibung: A new set of deep float trajectory data collected in the Gulf of Mexico from 2011 to 2015 at 1500- and 2500-m depths is analyzed to describe mesoscale processes, with particular attention paid to the western Gulf. Wavelet analysis is used to identify coherent eddies in the float trajectories, leading to a census of the basinwide coherent eddy population and statistics of the eddies’ kinematic properties. The eddy census reveals a new formation region for anticyclones off the Campeche Escarpment, located northwest of the Yucatan Peninsula. These eddies appear to form locally, with no apparent direct connection to the upper layer. Once formed, the eddies drift westward along the northern edge of the Sigsbee Abyssal Gyre, located in the southwestern Gulf of Mexico over the abyssal plain. The formation mechanism and upstream sources for the Campeche Escarpment eddies are explored: the observational data suggest that eddy formation may be linked to the collision of a Loop Current eddy with the western boundary of the Gulf. Specifically, the disintegration of a deep dipole traveling under the Loop Current eddy Kraken, caused by the interaction with the northwestern continental slope, may lead to the acceleration of the abyssal gyre and the boundary current in the Bay of Campeche region.

Beschreibung: The authors were supported by the Department of the Interior, Bureau of Ocean Energy Management (BOEM), Contract M10PC00112 to Leidos, Inc., Raleigh, North Carolina.

Beschreibung: 2019-05-07

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 31 (2018): 9881-9901, doi:10.1175/JCLI-D-17-0889.1.

Beschreibung: The Atlantic meridional overturning circulation and associated poleward heat transport are balanced by northern heat loss to the atmosphere and corresponding water-mass transformation. The circulation of northward-flowing Atlantic Water at the surface and returning overflow water at depth is particularly manifested—and observed—at the Greenland–Scotland Ridge where the water masses are guided through narrow straits. There is, however, a rich variability in the exchange of water masses across the ridge on all time scales. Focusing on seasonal and interannual time scales, and particularly the gateways of the Denmark Strait and between the Faroe Islands and Shetland, we specifically assess to what extent the exchanges of water masses across the Greenland–Scotland Ridge relate to wind forcing. On seasonal time scales, the variance explained of the observed exchanges can largely be related to large-scale wind patterns, and a conceptual model shows how this wind forcing can manifest via a barotropic, cyclonic circulation. On interannual time scales, the wind stress impact is less direct as baroclinic mechanisms gain importance and observations indicate a shift in the overflows from being more barotropically to more baroclinically forced during the observation period. Overall, the observed Greenland–Scotland Ridge exchanges reflect a horizontal (cyclonic) circulation on seasonal time scales, while the interannual variability more represents an overturning circulation.

Beschreibung: This research was supported by the Research Council of Norway project NORTH (Grant 229763). Additional support for M. A. Spall was provided by National Science Foundation Grant OCE- 1558742, for T. Eldevik and S. Østerhus by the European Union’s Horizon 2020 research and innovation program project Blue-Action (Grant 727852), and for S. Østerhus by the European Framework Programs under Grant Agreement 308299 (NACLIM).

Beschreibung: Author Posting. © The Authors, 2018. This article is posted here by permission of The Royal Astronomical Society for personal use, not for redistribution. The definitive version was published in Geophysical Journal International 215 (2018): 460–473, doi:10.1093/gji/ggy152.

Beschreibung: In this work, we present a new methodology to predict grain-size distributions from geophysical data. Specifically, electric conductivity and magnetic susceptibility of seafloor sediments recovered from electromagnetic profiling data are used to predict grain-size distributions along shelf-wide survey lines. Field data from the NW Iberian shelf are investigated and reveal a strong relation between the electromagnetic properties and grain-size distribution. The here presented workflow combines unsupervised and supervised machine-learning techniques. Non-negative matrix factorization is used to determine grain-size end-members from sediment surface samples. Four end-members were found, which well represent the variety of sediments in the study area. A radial basis function network modified for prediction of compositional data is then used to estimate the abundances of these end-members from the electromagnetic properties. The end-members together with their predicted abundances are finally back transformed to grain-size distributions. A minimum spatial variation constraint is implemented in the training of the network to avoid overfitting and to respect the spatial distribution of sediment patterns. The predicted models are tested via leave-one-out cross-validation revealing high prediction accuracy with coefficients of determination (R2) between 0.76 and 0.89. The predicted grain-size distributions represent the well-known sediment facies and patterns on the NW Iberian shelf and provide new insights into their distribution, transition and dynamics. This study suggests that electromagnetic benthic profiling in combination with machine learning techniques is a powerful tool to estimate grain-size distribution of marine sediments.

Beschreibung: This work was funded through DFG Research Center/Cluster of Excellence ‘The Ocean in the Earth System’ and was part of MARUM Research Area SD

Beschreibung: Author Posting. © American Meteorological Society, 2017. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Atmospheric and Oceanic Technology 34 (2017): 2673-2682, doi:10.1175/JTECH-D-17-0054.1.

Beschreibung: Expendable bathythermographs (XBT) to profile upper-ocean temperatures from vessels in motion have been in use for some 50 years now. Developed originally for navy use, they were soon adapted by oceanographers to map out upper-ocean thermal structure and its space–-time variability from both research vessels and merchant marine vessels in regular traffic. These activities continue today. This paper describes a new technology—the Autonomous Expendable Instrument System (AXIS)—that has been developed to provide the capability to deploy XBT probes on a predefined schedule, or adaptively in response to specific events without the presence of an observer on board. AXIS is a completely self-contained system that can hold up to 12 expendable probes [XBTs, XCTDs, expendable sound velocimeter (XSV)] in any combination. A single-board Linux computer keeps track of what probes are available, takes commands from ashore via Iridium satellite on what deployment schedule to follow, and records and forwards the probe data immediately with a time stamp and the GPS position. This paper provides a brief overview of its operation, capabilities, and some examples of how it is improving coverage along two lines in the Atlantic.

Beschreibung: Initial development of AXIS mechanical design elements wasmade possible by awards from the Cecil H. and Ida M. Green Technology Innovation Fund and the Sealark Foundation to the team of Dave Fratantoni, Keith von der Heydt (WHOI), and Terry Hammar (WHOI). Construction of the first full AXIS prototype was supported by a technology grant from the National Science Foundation (OCE-0926853) and the second one through an NSF-funded (OCE-1061185) subcontract from the University of Rhode Island.

Beschreibung: 2018-06-28

Beschreibung: Author Posting. © American Meteorological Society, 2017. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 47 (2017): 879-894, doi:10.1175/JPO-D-16-0196.1.

Beschreibung: Models show that surface cooling over a sloping continental shelf should give rise to baroclinic instability and thus tend toward gravitationally stable density stratification. Less is known about how alongshore winds affect this process, so the role of surface momentum input is treated here by means of a sequence of idealized, primitive equation numerical model calculations. The effects of cooling rate, wind amplitude and direction, bottom slope, bottom friction, and rotation rate are all considered. All model runs lead to instability and an eddy field. While instability is not strongly affected by upwelling-favorable alongshore winds, wind-driven downwelling substantially reduces eddy kinetic energy, largely because the downwelling circulation plays a similar role to baroclinic instability by flattening isotherms and so reducing available potential energy. Not surprisingly, cross-shelf winds appear to have little effect. Analysis of the model runs leads to quantitative relations for the wind effect on eddy kinetic energy for the equilibrium density stratification (which increases as the cooling rate increases) and for eddy length scale.

Beschreibung: This research was supported by the National Science Foundation Physical Oceanography Program through Grant OCE-1433953.

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Bulletin of the American Meteorological Society 99 (2018): 61-82, doi:10.1175/BAMS-D-16-0254.1.

Beschreibung: Accurate estimation of the climate sensitivity requires a better understanding of the nexus between polar marine ecosystem responses to warming, changes in sea ice extent, and emissions of marine biogenic aerosol (MBA). Sea ice brine channels contain very high concentrations of MBA precursors that, once ventilated, have the potential to alter cloud microphysical properties, such as cloud droplet number, and the regional radiative energy balance. In contrast to temperate latitudes, where the pelagic phytoplankton are major sources of MBAs, the seasonal sea ice dynamic plays a key role in determining MBA concentration in both the Arctic and Antarctic. We review the current knowledge of MBA sources and the link between ice melt and emissions of aerosol precursors in the polar oceans. We illustrate the processes by examining decadal-scale time series in various satellite-derived parameters such as aerosol optical depth (AOD), sea ice extent, and phytoplankton biomass in the sea ice zones of both hemispheres. The sharpest gradients in aerosol indicators occur during the spring period of ice melt. In sea ice–covered waters, the peak in AOD occurs well before the annual maximum in biomass in both hemispheres. The results provide strong evidence that suggests seasonal changes in sea ice and ocean biology are key drivers of the polar aerosol cycle. The positive trend in annual-mean Antarctic sea ice extent is now almost one-third of the magnitude of the annual-mean decrease in Arctic sea ice, suggesting the potential for different patterns of aerosol emissions in the future.

Beschreibung: Matrai and Middleton were supported by the National Science Foundation (PLR-1417517). Jones was partially funded by a grant from the Australian Research Council (DP150101649).

Beschreibung: Author Posting. © American Meteorological Society, 2017. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 47 (2017): 855-866, doi:10.1175/JPO-D-16-0194.1.

Beschreibung: Mesoscale eddies shape the Beaufort Gyre response to Ekman pumping, but their transient dynamics are poorly understood. Climate models commonly use the Gent–McWilliams (GM) parameterization, taking the eddy streamfunction to be proportional to an isopycnal slope s and an eddy diffusivity K. This local-in-time parameterization leads to exponential equilibration of currents. Here, an idealized, eddy-resolving Beaufort Gyre model is used to demonstrate that carries a finite memory of past ocean states, violating a key GM assumption. As a consequence, an equilibrating gyre follows a spiral sink trajectory implying the existence of a damped mode of variability—the eddy memory (EM) mode. The EM mode manifests during the spinup as a 15% overshoot in isopycnal slope (2000 km3 freshwater content overshoot) and cannot be explained by the GM parameterization. An improved parameterization is developed, such that is proportional to an effective isopycnal slope , carrying a finite memory γ of past slopes. Introducing eddy memory explains the model results and brings to light an oscillation with a period ≈ 50 yr, where the eddy diffusion time scale TE ~ 10 yr and γ ≈ 6 yr are diagnosed from the eddy-resolving model. The EM mode increases the Ekman-driven gyre variance by γ/TE ≈ 50% ± 15%, a fraction that stays relatively constant despite both time scales decreasing with increased mean forcing. This study suggests that the EM mode is a general property of rotating turbulent flows and highlights the need for better observational constraints on transient eddy field characteristics.

Beschreibung: GEM acknowledges the Stanback Postdoctoral Fellowship Fund at Caltech and the Howland Postdoctoral Program Fund at WHOI. MAS was supported by NSF Grants PLR-1415489 and OCE- 1232389. AFT acknowledges support from NSF OCE- 1235488.

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 2209-2219, doi:10.1175/JPO-D-18-0070.1.

Beschreibung: Published observations of subinertial ocean current variability show that the vertical structure is often well described by a vertical mode that has a node of horizontal velocity at the bottom rather than the traditional node of vertical velocity. The theory of forced and free linear Rossby waves in a continuously stratified ocean with a sloping bottom and bottom friction is treated here to see if frictional effects can plausibly contribute to this phenomenon. For parameter values representative of the mesoscale, bottom dissipation by itself appears to be too weak to be an explanation, although caution is required because the present approach uses a linear model to address a nonlinear phenomenon. One novel outcome is the emergence of a short-wave, bottom-trapped, strongly damped mode that is present even with a flat bottom.

Beschreibung: Partial funding for this article is provided by the National Science Foundation Physical Oceanography section through Award OCE-1433953.

Beschreibung: 2019-03-17

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 30 (2018): 1297-1314, doi:10.1175/JCLI-D-17-0286.1.

Beschreibung: The variance of a jet’s position in latitude is found to be related to its average speed: when a jet becomes stronger, its variability in latitude decreases. This relationship is shown to hold for observed midlatitude jets around the world and also across a hierarchy of numerical models. North Atlantic jet variability is shown to be modulated on decadal time scales, with decades of a strong, steady jet being interspersed with decades of a weak, variable jet. These modulations are also related to variations in the basinwide occurrence of high-impact blocking events. A picture emerges of complex multidecadal jet variability in which recent decades do not appear unusual. An underlying barotropic mechanism is proposed to explain this behavior, related to the change in refractive properties of a jet as it strengthens, and the subsequent effect on the distribution of Rossby wave breaking.

Beschreibung: We would like to acknowledge funding from NERC and the Research Council of Norway project jetSTREAM under Grants NE/ L01047X/1 (IMPETUS) and 231716, respectively, for a contribution to the work presented here. EAB is supported in part by the NSF Climate and Large-Scale Dynamics Program under Grant 1545675. Y-OK was supported by the NSF Climate and Large-Scale Dynamics Program under Grant 1355339. KW was supported by the Joint UK BEIS/Defra Met Office Hadley Centre Climate Programme (GA01101). RL was supported by the Met Office and the National Centre for Atmospheric Science.

Beschreibung: 2018-07-29

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 31 (2018): 4309-4327, doi:10.1175/JCLI-D-17-0407.1.

Beschreibung: Multidecadal hydroclimate variability has been expressed as “megadroughts” (dry periods more severe and prolonged than observed over the twentieth century) and corresponding “megapluvial” wet periods in many regions around the world. The risk of such events is strongly affected by modes of coupled atmosphere–ocean variability and by external impacts on climate. Accurately assessing the mechanisms for these interactions is difficult, since it requires large ensembles of millennial simulations as well as long proxy time series. Here, the Community Earth System Model (CESM) Last Millennium Ensemble is used to examine statistical associations among megaevents, coupled climate modes, and forcing from major volcanic eruptions. El Niño–Southern Oscillation (ENSO) strongly affects hydroclimate extremes: larger ENSO amplitude reduces megadrought risk and persistence in the southwestern United States, the Sahel, monsoon Asia, and Australia, with corresponding increases in Mexico and the Amazon. The Atlantic multidecadal oscillation (AMO) also alters megadrought risk, primarily in the Caribbean and the Amazon. Volcanic influences are felt primarily through enhancing AMO amplitude, as well as alterations in the structure of both ENSO and AMO teleconnections, which lead to differing manifestations of megadrought. These results indicate that characterizing hydroclimate variability requires an improved understanding of both volcanic climate impacts and variations in ENSO/AMO teleconnections.

Beschreibung: This work is supported by NSF EaSM Grants AGS-1243125 and NCAR-1243107 to The University of Arizona.

Beschreibung: 2018-11-03

Beschreibung: Author Posting. © American Meteorological Society, 2018. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 31 (2018): 8627-8643, doi:10.1175/JCLI-D-18-0010.1.

Beschreibung: Drought has severe consequences for humans and their environment, yet we have a limited understanding of the drivers of drought across the full range of time scales on which it occurs. Here, the atmosphere and ocean conditions that drive this continuum of drought variability in southwestern North America (SWNA) are studied using the latest observationally based products, paleoclimate reconstructions, and state-of-the-art Earth system model simulations of the last millennium. A novel application of the self-organizing maps (SOM) methodology allows for a visualization of the continuum of climate states coinciding with thousands of droughts of varying lengths in last millennium simulations from the Community Earth System Model (CESM), the Goddard Institute for Space Studies Model E2-R (GISS E2-R), and eight other members from phase 5 of the Coupled Model Intercomparison Project (CMIP5). It is found that most droughts are associated with a cool Pacific decadal oscillation (PDO) pattern, but persistent droughts can coincide with a variety of ocean–atmosphere states, including time periods showing a warm PDO or weak ocean–atmosphere anomalies. Many CMIP5 models simulate similar SWNA teleconnection patterns, but the SOM analysis demonstrates that models simulate different continuums of ocean–atmosphere states coinciding with droughts of different lengths, suggesting fundamental differences in their drought dynamics. These findings have important implications for our understanding and simulation of the drivers of persistent drought, and for their potential predictability.

Beschreibung: The National Science Foundation EaSM2 Grant (AGS1243125) supported this work.

Beschreibung: Author Posting. © Author(s), 2017. This article is posted here by permission of Oxford University Press for personal use, not for redistribution. The definitive version was published in Geophysical Journal International 212 (2018): 1429–1449, doi:10.1093/gji/ggx488.

Beschreibung: We conducted detailed analyses of a global array of trenches, revealing systematic intra- and intertrench variations in plate bending characteristics. The intratrench variations of the Manila and Mariana Trenches were analysed in detail as end-member cases of the relatively young (16–36 Ma) and old (140–160 Ma) subducting plates, respectively. Meanwhile, the intertrench variability was investigated for a global array of additional trenches including the Philippine, Kuril, Japan, Izu-Bonin, Aleutian, Tonga-Kermadec, Middle America, Peru, Chile, Sumatra and Java Trenches. Results of the analysis show that the trench relief (W0) and width (X0) of all systems are controlled primarily by the faulting-reduced elastic thickness near the trench axis (Tme) and affected only slightly by the initial unfaulted thickness (TMe) of the incoming plate. The reduction in Te has caused significant deepening and narrowing of trench valleys. For the cases of relatively young or old plates, the plate age could be a dominant factor in controlling the trench bending shape, regardless the variations in axial loadings. Our calculations also show that the axial loading and stresses of old subducting plates can vary significantly along the trench axis. In contrast, the young subducting plates show much smaller values and variations in axial loading and stresses.

Beschreibung: This work was supported by Chinese Academy of Sciences Grants (Y4SL021001, QYZDY-SSW-DQC005, YZ201325 and YZ201534), National Natural Science Foundation of China Grants (91628301, U1606401, 41376063 and 41706056) and HKSAR Research Grant Council Grants (24601515, 14313816).

Beschreibung: Author Posting. © American Meteorological Society, 2017. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 47 (2017): 2999-3013, doi:10.1175/JPO-D-17-0129.1.

Beschreibung: Initial results are presented from a yearlong, high-resolution (~2 km) numerical simulation covering the east Greenland shelf and the Iceland and Irminger Seas. The model hydrography and circulation in the vicinity of Denmark Strait show good agreement with available observational datasets. This study focuses on the variability of the Denmark Strait overflow (DSO) by detecting and characterizing boluses and pulses, which are the two dominant mesoscale features in the strait. The authors estimate that the yearly mean southward volume flux of the DSO is about 30% greater in the presence of boluses and pulses. On average, boluses (pulses) are 57.1 (27.5) h long, occur every 3.2 (5.5) days, and are more frequent during the summer (winter). Boluses (pulses) increase (decrease) the overflow cross-sectional area, and temperatures around the overflow interface are colder (warmer) by about 2.6°C (1.8°C). The lateral extent of the boluses is much greater than that of the pulses. In both cases the along-strait equatorward flow of dense water is enhanced but more so for pulses. The sea surface height (SSH) rises by 4–10 cm during boluses and by up to 5 cm during pulses. The SSH anomaly contours form a bowl (dome) during boluses (pulses), and the two features cross the strait with a slightly different orientation. The cross streamflow changes direction; boluses (pulses) are associated with veering (backing) of the horizontal current. The model indicates that boluses and pulses play a major role in controlling the variability of the DSO transport into the Irminger Sea.

Beschreibung: This work was supported by the NSF Grants OCE-1433448, OCE-1633124, and OCE- 1259618 and the Institute for Data Intensive Engineering and Science (IDIES) seed grant funding.

Beschreibung: 2018-06-13

Beschreibung: Author Posting. © American Meteorological Society, 2017. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 47 (2017): 2927-2947, doi:10.1175/JPO-D-17-0083.1.

Beschreibung: Motivated by observations in Hudson shelf valley showing stronger onshore than offshore flows, this study investigates wind-driven flows in idealized shallow shelf valleys. This first part of a two-part sequence focuses on the mechanism of the asymmetrical flow response in a valley to along-shelf winds of opposite directions. Model simulations show that (i) when the wind is in the opposite direction to coastal-trapped wave (CTW) phase propagation, the shelf flow turns onshore in the valley and generates strong up-valley transport and a standing meander on the upstream side (in the sense of CTW phase propagation) of the valley, and (ii) when the wind is in the same direction as CTW phase propagation, the flow forms a symmetric onshore detour pattern over the valley with negligible down-valley transport. Comparison of the modeled upstream meanders in the first scenario with CTW characteristics confirms that the up-valley flow results from CTWs being arrested by the wind-driven shelf flow establishing lee waves. The valley bathymetry generates an initial excessive onshore pressure gradient force that drives the up-valley flow and induces CTW lee waves that sustain the up-valley flow. When the wind-driven shelf flow aligns with CTW phase propagation, the initial disturbance generated in the valley propagates away, allowing the valley flow to adjust to roughly follow isobaths. Because of the similarity in the physical setup, this mechanism of arrested CTWs generating stronger onshore than offshore flow is expected to be applicable to the flow response in slope canyons to along-isobath background flows of opposite directions.

Beschreibung: WGZ and SJL were supported by the National Science Foundation through GrantOCE1154575.WGZ is also supported by the NSF Grant OCE 1634965 and SJL by NSF Grant OCE 1558874.

Beschreibung: 2018-06-08

Beschreibung: Author Posting. © American Meteorological Society, 2017. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 48 (2018): 29-44, doi:10.1175/JPO-D-17-0016.1.

Beschreibung: The cospectrum of the horizontal and vertical turbulent velocity fluctuations, an essential tool for understanding measurements of the turbulent Reynolds shear stress, often departs in the ocean from the shape that has been established in the atmospheric surface layer. Here, we test the hypothesis that this departure is caused by advection of standard boundary layer turbulence by the random oscillatory velocities produced by surface gravity waves. The test is based on a model with two elements. The first is a representation of the spatial structure of the turbulence, guided by rapid distortion theory, and consistent with the one-dimensional cospectra that have been measured in the atmosphere. The second model element is a map of the spatial structure of the turbulence to the temporal fluctuations measured at fixed sensors, assuming advection of frozen turbulence by the velocities associated with surface waves. The model is adapted to removal of the wave velocities from the turbulent fluctuations using spatial filtering. The model is tested against previously published laboratory measurements under wave-free conditions and two new sets of measurements near the seafloor in the coastal ocean in the presence of waves. Although quantitative discrepancies exist, the model captures the dominant features of the laboratory and field measurements, suggesting that the underlying model physics are sound.

Beschreibung: This research was supported by National Science Foundation Ocean Sciences Division Award 1356060 and the U.S. Geological Survey Coastal and Marine Geology Program.

Beschreibung: Author Posting. © American Meteorological Society, 2017. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 30 (2017): 9871-9895, doi:10.1175/JCLI-D-17-0009.1.

Beschreibung: Two large ensembles (LEs) of historical climate simulations are used to compare how various statistical methods estimate the sea surface temperature (SST) changes due to anthropogenic and other external forcing, and how their removal affects the internally generated Atlantic multidecadal oscillation (AMO), Pacific decadal oscillation (PDO), and the SST footprint of the Atlantic meridional overturning circulation (AMOC). Removing the forced SST signal by subtracting the global mean SST (GM) or a linear regression on it (REGR) leads to large errors in the Pacific. Multidimensional ensemble empirical mode decomposition (MEEMD) and quadratic detrending only efficiently remove the forced SST signal in one LE, and cannot separate the short-term response to volcanic eruptions from natural SST variations. Removing a linear trend works poorly. Two methods based on linear inverse modeling (LIM), one where the leading LIM mode represents the forced signal and another using an optimal perturbation filter (LIMopt), perform consistently well. However, the first two LIM modes are sometimes needed to represent the forced signal, so the more robust LIMopt is recommended. In both LEs, the natural AMO variability seems largely driven by the AMOC in the subpolar North Atlantic, but not in the subtropics and tropics, and the scatter in the AMOC–AMO correlation is large between individual ensemble members. In three observational SST reconstructions for 1900–2015, linear and quadratic detrending, MEEMD, and GM yield somewhat different AMO behavior, and REGR yields smaller PDO amplitudes. Based on LIMopt, only about 30% of the AMO variability is internally generated, as opposed to more than 90% for the PDO. The natural SST variability contribution to global warming hiatus is discussed.

Beschreibung: Support from the NOAA Climate Program Office Climate Variability and Predictability program (NA13OAR4310139), NSF EaSM2 (OCE-84298900), the European Community Horizon 2020 Framework under Grant Agreement 727852 (Blue-Action), and the ANR MORDICUS project (ANR-13-SENV-0002-02) is gratefully acknowledged.

Beschreibung: 2018-05-16

Beschreibung: When cumulonimbus clouds aggregate, developing into a single entity with precipitation covering a horizontal scale of hundreds of kilometers, they are called mesoscale convective systems (MCSs). They account for much of Earth’s precipitation, generate severe weather events and flooding, produce prodigious cirriform anvil clouds, and affect the evolution of the larger-scale circulation. Understanding the inner workings of MCSs has resulted from developments in observational technology and modeling. Time–space conversion of ordinary surface and upper-air observations provided early insight into MCSs, but deeper understanding has followed field campaigns using increasingly sophisticated radars, better aircraft instrumentation, and an ever-widening range of satellite instruments, especially satellite-borne radars. High-resolution modeling and theoretical insights have shown that aggregated cumulonimbus clouds induce a mesoscale circulation consisting of air overturning on a scale larger than the scale of individual convective up- and downdrafts. These layers can be kilometers deep and decoupled from the boundary layer in elevated MCSs. Cooling in the lower troposphere and heating aloft characterize the stratiform regions of MCSs. As a result, long-lived MCSs with large stratiform regions have a top-heavy heating profile that generates potential vorticity in midlevels, thus influencing the larger-scale circulation within which the MCSs occur. Global satellite data show MCSs varying in structure, depending on the prevailing large-scale circulation and topography. These patterns are likely to change with global warming. In addition, environmental pollution affects MCS structure and dynamics subtly. Feedbacks of MCSs therefore need to be included or parameterized in climate models.

Beschreibung: The central change in understanding of the ocean circulation during the past 100 years has been its emergence as an intensely time-dependent, effectively turbulent and wave-dominated, flow. Early technologies for making the difficult observations were adequate only to depict large-scale, quasi-steady flows. With the electronic revolution of the past 50+ years, the emergence of geophysical fluid dynamics, the strongly inhomogeneous time-dependent nature of oceanic circulation physics finally emerged. Mesoscale (balanced), submesoscale oceanic eddies at 100-km horizontal scales and shorter, and internal waves are now known to be central to much of the behavior of the system. Ocean circulation is now recognized to involve both eddies and larger-scale flows with dominant elements and their interactions varying among the classical gyres, the boundary current regions, the Southern Ocean, and the tropics.

Beschreibung: Although atmospheric observing systems were already an important part of meteorology before the American Meteorological Society was established in 1919, the past 100 years have seen a steady increase in their numbers and types. Examples of how observing systems were developed and how they have enabled major scientific discoveries are presented. These examples include observing systems associated with the boundary layer, the upper air, clouds and precipitation, and solar and terrestrial radiation. Widely used specialized observing systems such as radar, lidar, and research aircraft are discussed, and examples of applications to weather forecasting and climate are given. Examples drawn from specific types of chemical measurements, such as ozone and carbon dioxide, are included. Sources of information on observing systems, including other chapters of this monograph, are also discussed. The past 100 years has been characterized by synergism between societal needs for weather observations and the needs of fundamental meteorological research into atmospheric processes. In the latter half of the period, observing system improvements have been driven by the increasing demands for higher-resolution data for numerical models, the need for long-term measurements, and for more global coverage. This has resulted in a growing demand for data access and for integrating data from an increasingly wide variety of observing system types and networks. These trends will likely continue.

Beschreibung: The polar regions present several unique challenges to meteorology, including remoteness and a harsh environment. We summarize the evolution of polar meteorology in both hemispheres, beginning with measurements made during early expeditions and concluding with the recent decades in which polar meteorology has been central to global challenges such as the ozone hole, weather prediction, and climate change. Whereas the 1800s and early 1900s provided data from expeditions and only a few subarctic stations, the past 100 years have seen great advances in the observational network and corresponding understanding of the meteorology of the polar regions. For example, a persistent view in the early twentieth century was of an Arctic Ocean dominated by a permanent high pressure cell, a glacial anticyclone. With increased observations, by the 1950s it became apparent that, while anticyclones are a common feature of the Arctic circulation, cyclones are frequent and may be found anywhere in the Arctic. Technology has benefited polar meteorology through advances in instrumentation, especially autonomously operated instruments. Moreover, satellite remote sensing and computer models revolutionized polar meteorology. We highlight the four International Polar Years and several high-latitude field programs of recent decades. We also note outstanding challenges, which include understanding of the role of the Arctic in variations of midlatitude weather and climate, the ability to model surface energy exchanges over a changing Arctic Ocean, assessments of ongoing and future trends in extreme events in polar regions, and the role of internal variability in multiyear-to-decadal variations of polar climate.

Beschreibung: The field of atmospheric science has been enhanced by its long-standing collaboration with entities with specific needs. This chapter and the two subsequent ones describe how applications have worked to advance the science at the same time that the science has served the needs of society. This chapter briefly reviews the synergy between the applications and advancing the science. It specifically describes progress in weather modification, aviation weather, and applications for security. Each of these applications has resulted in enhanced understanding of the physics and dynamics of the atmosphere, new and improved observing equipment, better models, and a push for greater computing power.

Beschreibung: Over the past century, the atmospheric and related sciences have seen incredible advances in our understanding of Earth’s environment and our ability to monitor and predict its behavior. These advances have had a profound impact on society and have been integrated into every aspect of daily life. The American Meteorological Society (AMS) has been instrumental in supporting these advances throughout its first 100 years of existence as a scientific and professional society serving the community of professionals in the atmospheric and related oceanic and hydrologic sciences. AMS has provided opportunities for researchers and practitioners to share their scientific findings and build fruitful collaborations to further the science and its application. Through strategic initiatives at key points in its history, AMS has pushed the science forward—highlighting areas ripe for development, creating frameworks for interdisciplinary interactions, and providing innovative approaches to the dissemination of research results. As a society made up of the scientific community and led by many of the most prominent scientists of their time, AMS has been able to respond to, and often anticipate, the needs of its community.

Beschreibung: A century ago, meteorologists regarded tropical cyclones as shallow vortices, extending upward only a few kilometers into the troposphere, and nothing was known about their physics save that convection was somehow involved. As recently as 1938, a major hurricane struck the densely populated northeastern United States with no warning whatsoever, killing hundreds. In the time since the American Meteorological Society was founded, however, tropical cyclone research blossomed into an endeavor of great breadth and depth, encompassing fields ranging from atmospheric and oceanic dynamics to biogeochemistry, and the precision and scope of forecasts and warnings have achieved a level of success that would have been regarded as impossible only a few decades ago. This chapter attempts to document the extraordinary progress in tropical cyclone research over the last century and to suggest some avenues for productive research over the next one.

Beschreibung: 〈span〉〈div〉Summary〈/div〉The strength of the lithosphere plays a key role in the formation and evolution of tectonic plate boundaries. Localized lithospheric deformation associated with plate tectonics requires a mechanism for weakening across the entire width of the lithosphere, including the strongest cold ductile region. We explore the microphysics of weakening of lithospheric materials, and in particular the coupled evolution of mineral grain size and intragranular defects and their control on lithospheric strength. We propose a model for the interaction between grain-boundaries and dislocation density to reduce the net free energy of grains during dynamic recrystallization (DRX). The driving forces for DRX arise from heterogeneity in dislocation density and grain boundary curvature. Our model shows that grain growth driven by variation in grain boundary curvature can be impeded by variation in dislocation density; this occurs because as the grains grow, to minimize their surface energy, their dislocation density and associated internal energy may increase and offset the driving forces for grain growth. The correlation between grain size and dislocation density can for example arise because the dislocation accumulation in smaller grains is suppressed due to the large stress that is needed to bend and elongate a short dislocation (as dictated by the small grain size), while the larger grains can have long dislocations and reach a steady state dislocation density dictated by the applied stress. In a lithospheric setting, slower grain growth means that it would require less mechanical work to establish weak localized shear zones through grain damage, and retard the healing of previously damaged zones. Furthermore, the competition of two different time-scales - that of grain growth and the dislocation kinetics - can lead to oscillating behavior over 1 to 10 years as the grain size and dislocation density advance towards their steady states. These oscillations are likely to have an effect on the rheology of lithospheric rocks, e.g. their strengthening and weakening through time, and have a potential application to geological processes such as postseismic creep in ductile shear zones.〈/span〉

Beschreibung: 〈span〉〈div〉Summary〈/div〉Locating and monitoring passive seismic sources provides us important information for studying subsurface rock deformation, injected fluid migration, regional stress conditions as well as fault rupture mechanism. In this paper, we present a novel passive-source monitoring approach using vector-based elastic time-reversal imaging. By solving the elastic wave equation using observed multicomponent records as boundary conditions, we first compute back-propagated elastic wavefields in the subsurface. Then, we separate the extrapolated wavefields into compressional (P-wave) and shear (S-wave) modes using the vector Helmholtz decomposition. A zero-lag cross-correlation imaging condition is applied to the separated pure-mode vector wavefields to produce passive-source images. We compare imaging results using three implementations, i.e., dot-product, energy and power. Numerical experiments demonstrate that the power imaging condition gives us the highest resolution and is less sensitive to the presence of random noises. To capture the propagation of microseismic fracture and earthquake rupture, we modify the traditional zero-lag cross-correlation imaging condition by summing the multiplication of the separated P- and S-wavefields within local time windows, which enables us to capture the temporal and spatial evolution of earthquake rupture. 2D and 3D numerical examples demonstrate that the proposed method is capable of accurately locating point sources, as well as delineating dynamic propagation of hydraulic fracture and earthquake rupture.〈/span〉

Beschreibung: 〈span〉〈div〉Summary〈/div〉Magnitudes of differential stress in the lithosphere, especially in the crust, are still disputed. Earthquake-based stress drop estimates indicate median values 〈 10 MPa whereas the lateral variation of gravitational potential energy per unit area (〈span〉GPE〈/span〉) across significant relief indicates stress magnitudes of ca. 100 MPa in average across a 100 km thick lithosphere between the Indian lowland and the Tibetan plateau. These standard 〈span〉GPE〈/span〉-based stress estimates correspond to membrane stresses, because they are associated with a deformation that is uniform with depth. We show here with new analytical results that lateral variations in 〈span〉GPE〈/span〉 can also cause bending moments and related bending stresses of several hundreds of MPa. Furthermore, we perform two-dimensional thermo-mechanical numerical simulations (1) to evaluate estimates for membrane and bending stresses based on 〈span〉GPE〈/span〉 variations, (2) to quantify minimum crustal stress magnitudes that are required to maintain the topographic relief between Indian lowland and Tibetan plateau for ca. 10 Ma and (3) to quantify the corresponding relative contribution of crustal strength to the total lithospheric strength. The numerical model includes viscoelastoplastic deformation, gravity and heat transfer. The model configuration is based on density fields from the CRUST1.0 data set and from a geophysically and petrologically constrained density model based on 〈span〉in situ〈/span〉 field campaigns. The numerical results indicate that values of differential stress in the upper crust must be 〉 ca. 180 MPa, corresponding to a friction angle of ca. 10°, to maintain the topographic relief between lowland and plateau for 〉 10 Ma. The relative contribution of crustal strength to total lithospheric strength varies considerably laterally. In the region between lowland and plateau and inside the plateau the depth-integrated crustal strength is approximately equal to the depth-integrated strength of the mantle lithosphere. Simple analytical formulae predicting the lateral variation of depth-integrated stresses agree with numerically calculated stress fields, which show both the accuracy of the numerical results and the applicability of simple, rheology-independent, analytical predictions to highly variable, rheology-dependent, stress fields. Our results indicate that (1) crustal strength can be locally equal to mantle lithosphere strength and that (2) crustal stresses must be at least one order of magnitude larger than median stress drops in order to support the plateau relief over a duration of ca. 10 Ma.〈/span〉

Beschreibung: 〈span〉〈div〉Summary〈/div〉Scanning magnetometers are increasingly used to characterize the magnetization of mineral grains in rock samples. Up-scaling this measurement technique to large numbers of individual particles is hampered by the intrinsic non-uniqueness of potential-field inversion. Here it is shown that this problem can be circumvented by adding tomographic information that determines the location of the possible field sources. Standard potential theory is used to prove a uniqueness theorem which completely characterizes the mathematical background of the corresponding source-localized inversion. It exactly resolves under which conditions a potential field measurement on a surface can be uniquely decomposed into signals from the different source regions. The intrinsic non-uniqueness of potential field inversion prevents that the source distribution inside the tomographically outlined regions can be recovered, but the potential field of each region is uniquely defined. For scanning magnetometers in rock magnetism, this result implies that magnetic dipole vectors of large numbers of individual magnetic particles can be reliably reconstructed from surface scans of the magnetic field, if the particle positions are independently determined. This provides an incentive to improve scanning methods for future paleomagnetic applications.〈/span〉

Beschreibung: 〈span〉〈div〉Summary〈/div〉The most common earthquake forecasting models assume that the magnitude of the next earthquake is independent from the past. This feature severely limits the capability to forecast large earthquakes with high probabilities. Here we investigate empirically on the magnitude-independence assumption, exploring if: i) background and triggered earthquakes have the same frequency-magnitude distribution, ii) variations of seismicity in the space-time-magnitude domain encode some information on the future earthquakes size. For this purpose, and to verify the stability of the findings, we consider seismic catalogues covering different space-time-magnitude windows, such as the Alto Tiberina Near Fault Observatory (TABOO), the California and Japanese seismic catalogues. Our approach is inspired by the nearest-neighbour method proposed by Baiesi & Paczuski (2004) and elaborated by Zaliapin et al. (2008) to distinguish between triggered and background earthquakes. Here we implement the same metric-based correlation to identify the precursory seismicity of any triggered earthquake; this allows us to analyse, for each triggered earthquake, the space-time-magnitude distribution of the seismicity that likely contributed to its occurrence. Our results show that the magnitude-independence assumption holds reasonably well in all catalogues, with a remarkable exception that is consistent with a previous independent study; this departure from the magnitude-independence assumption shows that larger events tend to nucleate at a higher distance from the ongoing sequence. We also notice that the reliability of this assumption may depend on the spatial scale considered; it holds for seismic catalogues of large areas, but we identify possible departures in small areas, reflecting different ways to release locally seismic energy. Finally, we come across an important issue that may lead to misleading results in similar studies, i.e., if a seismic catalogue appears overall complete above a fixed magnitude threshold, it may still yield spurious signals into the analysis. Specifically, we show that some significant departures from the magnitude-independence assumption do not survive when considering spatiotemporal variations of the magnitude of completeness.〈/span〉

Beschreibung: 〈span〉〈div〉Summary〈/div〉Vp/Vs models provide important complementary information to Vp and Vs models, relevant to lithology, rock damage, partial melting, water saturation, etc. However, seismic tomography using body-wave traveltime data from local or regional earthquakes does not constrain Vp/Vs well due to the different resolution of Vp and Vs models, with the Vp models usually better constrained than Vs. Since surface-wave data are most sensitive to Vs, which leads to complementary strengths with respect to body-wave data, we jointly invert body- and surface-wave data to better resolve the Vp/Vs models. In order to show the robustness of our joint inversion method, we compare the results to other approaches, including dividing Vp by Vs models and Vp/Vs parameterization with body-wave or both body- and surface-wave data, using synthetic data and real data from the southern California plate boundary region. We confirm that Vp/Vs models from body-wave inversion obtained by division tend to include artifacts, even when both Vp and Vs models seem reasonable. The joint inversion with Vp/Vs parameterization is found to improve the Vp/Vs ratio model significantly and the resultant Vp/Vs model shows more geologically consistent features, such as high Vp/Vs along fault traces at shallow depths likely indicating fault-related damage. The Vp/Vs model also exhibits contrasts at intermediate depths along the Peninsular Range compositional boundary, and high Vp/Vs in the lower crust near the Salton Sea region correlated with high heat flow and may indicate partial melting. The improved Vp/Vs as well as individual Vp and Vs models are useful for earthquake relocation, high-resolution Moho depth imaging, and interpretation of other data and tectonic evolution in the region.〈/span〉

Beschreibung: 〈span〉〈div〉Summary〈/div〉Interpretation of surface fault scarps and palaeoseismic trenches is a key component of estimating fault slip rates, earthquake recurrence rates and maximum magnitudes for hazard assessments. Often these analyses rely on the assumption that successive earthquakes all breached the surface and that the ruptures are recorded topographically, or by the deposits exposed in a trench. The M〈sub〉〈span〉w〈/span〉〈/sub〉7.2 1992 Suusamyr earthquake, Kyrgyzstan, is an apparently problematic case for such analyses because its ruptures show significant displacement but are only mapped as having broken the surface along small, disparate portions of the fault. Here we present the results of surveys conducted along the Suusamyr Fault to establish whether that is the case. Two sets of ruptures were identified following the earthquake. They are unusually short for their displacement and are separated by a 25 km gap. Using satellite imagery, high-resolution digital elevation models and palaeoseismic trenching we first reassess the distribution of the 1992 ruptures and then reconstruct the Holocene earthquake record to establish the extent to which the 1992 earthquake is representative of the rupture behaviour of this fault. We find evidence for at least two prehistoric surface rupturing earthquakes in the Holocene: one ∼3 ka and one 〉8 ka that, along with the modern event, gives recurrence intervals of ∼3 kyr and ∼5 kyr. Within spatial gaps between segments of the 1992 ruptures there are clear prehistoric surface ruptures and the ruptures in each prehistoric earthquake were discontinuous. We conclude that there is significant variability in the surface rupture pattern of successive earthquakes on the Suusamyr Fault, with implications for the completeness of palaeoseismic records obtained from thrust scarps.〈/span〉

Beschreibung: 〈span〉〈div〉Summary〈/div〉We provide a six-component (6-C) polarization model for 〈span〉P〈/span〉-, 〈span〉SV〈/span〉-, 〈span〉SH〈/span〉-, Rayleigh-, and Love-waves both inside an elastic medium as well as at the free surface. It is shown that single-station 6-C data comprised of three components of rotational motion and three components of translational motion provide the opportunity to unambiguously identify the wave type, propagation direction, and local 〈span〉P〈/span〉- and 〈span〉S〈/span〉-wave velocities at the receiver location by use of polarization analysis. To extract such information by conventional processing of three-component (3-C) translational data would require large and dense receiver arrays. The additional rotational components allow the extension of the rank of the coherency matrix used for polarization analysis. This enables us to accurately determine the wave type and wave parameters (propagation direction and velocity) of seismic phases, even if more than one wave is present in the analysis time window. This is not possible with standard, pure-translational 3-C recordings. In order to identify modes of vibration and to extract the accompanying wave parameters, we adapt the multiple signal classification algorithm (MUSIC). Due to the strong nonlinearity of the MUSIC estimator function, it can be used to detect the presence of specific wave types within the analysis time window at very high resolution. We show how the extracted wavefield properties can be used, in a fully automated way, to separate the wavefield into its different wave modes using only a single 6-C recording station. As an example, we apply the method to remove surface wave energy while preserving the underlying reflection signal and to suppress energy originating from undesired directions, such as side-scattered waves.〈/span〉

Beschreibung: 〈span〉〈div〉Summary〈/div〉The Sentinel-1 mission comprises two synthetic aperture radar satellites, each with a 12 day orbital repeat, orbiting 6 days apart within a narrow tube. The mission design promises the ability to respond quickly to earthquakes with InSAR, and to facilitate production of interferograms with good interferometric correlation globally. We report on our efforts to study global seismicity using Sentinel-1 Interferometric Wide-Swath data between April 2015 and December 2016. We select 35 potentially detectable terrestrial earthquakes in the range 5.5 ≤ 〈span〉Mw〈/span〉 ≤ 7.8 on the basis of their locations, depths and magnitudes, and process the first post-event interferogram with the shortest possible time-span for each using the ISCE software. We evaluate each interferogram for earthquake deformation signals by visual inspection. We can identify deformation signals attributable to earthquakes in 18 of these interferograms (51%); a further six interferograms (17%) have ambiguous interferometric phase affected by tropospheric noise. 11 events (31%) could not be identified from their interferograms. The majority of these failed detections were due to interferogram decorrelation, particularly apparent for earthquakes that occurred between 15°N and 15°S, where climate conditions promote dense vegetation. The majority of the ambiguous interferograms are affected by tropospheric noise, suggesting that techniques to mitigate such noise could improve detection performance. The largest event we do not detect with Sentinel-1 data is a 〈span〉Mw〈/span〉7.0 earthquake that occurred in Vanuatu in April 2016; we also fail to detect the 2016 〈span〉Mw〈/span〉6.2 Kurayoshi earthquake in one out of two possible 24-day interferograms. We propose these as upper and lower estimates on the magnitude of completeness for earthquakes studied with Sentinel-1 data; to raise the magnitude of completeness we suggest that more frequent (e.g. six day) recurrence may be necessary in low latitude areas.〈/span〉

Beschreibung: 〈span〉〈div〉Summary〈/div〉The correct estimation of site-specific attenuation is crucial for the assessment of seismic hazard. Downhole instruments provide in this context valuable information to constrain attenuation directly from data. In this study, we apply an interferometric approach to this problem by deconvolving seismic motions recorded at depth with those recorded at the surface. In doing so, incident and surface-reflected waves can be separated. We apply this technique not only to earthquake data but also to recordings of ambient vibrations. We compute the transfer function between incident and surface-reflected waves in order to infer frequency dependent quality factors for S-waves. The method is applied to a 87 m deep borehole sensor and a co-located surface instrument situated at a hard-rock site in West Bohemia/Vogtland, Germany. We show that the described method provides comparable attenuation estimates using either earthquake data or ambient noise for frequencies between 5-15 Hz. Moreover, a single hour of noise recordings seems to be sufficient to yield stable deconvolution traces and quality factors, thus, offering a fast and easy way to derive attenuation estimates from borehole recordings even in low to mid seismicity regions.〈/span〉

Beschreibung: 〈span〉〈div〉Summary〈/div〉Injection of CO〈sub〉2〈/sub〉 into tight reservoirs produces both gravity change and ground deformation, which provide great opportunities for more accurate coupled inverse modelling. In this study, we incorporate signals generated from several synthetic models to estimate the CO〈sub〉2〈/sub〉 distribution in the reservoir. A relationship is found that connects density variations to volumetric changes associated with injected CO〈sub〉2〈/sub〉, taking advantage of a common set of model parameters for both gravitational and geo-mechanical inverse modelling. This is achieved by assuming that the injected CO〈sub〉2〈/sub〉 increases pressure in the reservoir, which in turn generates extra porosity that is then filled in by the CO〈sub〉2〈/sub〉 mass in the generated space. Tikhonov regularization, supported by the Generalized Cross Validation (GCV) technique for finding the optimized model, is used to solve the ill-posed inverse problems. The results indicate that with a combination of gravity and ground deformation monitoring, the uncertainty and ambiguity in gravimetric modelling due to high levels of noise is mitigated by implementing highly accurate ground deformation measurements, which normally have a higher signal to noise ratio.〈/span〉

Beschreibung: 〈span〉〈div〉Summary〈/div〉As the number of seismic sensors grows, it is becoming increasingly difficult for analysts to pick seismic phases manually and comprehensively, yet such efforts are fundamental to earthquake monitoring. Despite years of improvements in automatic phase picking, it is difficult to match the performance of experienced analysts. A more subtle issue is that different seismic analysts may pick phases differently, which can introduce bias into earthquake locations. We present a deep-neural-network-based arrival-time picking method called ”PhaseNet” that picks the arrival times of both P and S waves. Deep neural networks have recently made rapid progress in feature learning, and with sufficient training, have achieved super-human performance in many applications. PhaseNet uses three-component seismic waveforms as input and generates probability distributions of P arrivals, S arrivals, and noise as output. We engineer PhaseNet such that peaks in the probability distributions provide accurate arrival times for both P and S waves. PhaseNet is trained on the prodigious available data set provided by analyst-labeled P and S arrival times from the Northern California Earthquake Data Center. The dataset we use contains more than seven hundred thousand waveform samples extracted from over thirty years of earthquake recordings. We demonstrate that PhaseNet achieves much higher picking accuracy and recall rate than existing methods when applied to the waveforms of known earthquakes, which has the potential to increase the number of S-wave observations dramatically over what is currently available. This will enable both improved locations and improved shear wave velocity models.〈/span〉

Beschreibung: 〈span〉〈div〉Summary〈/div〉The low frequency earthquakes (LFEs) that constitute tectonic tremor are often inferred to be slow: to have durations of 0.2 to 0.5 s, a factor of 10 to 100 longer than those of typical 〈span〉MW〈/span〉 1-2 earthquakes. Here we examine LFEs near Parkfield, CA in order to assess several proposed explanations for LFEs’ long durations. We determine LFE rupture areas and location distributions using a new approach, similar to directivity analysis, where we examine how signals coming from various locations within LFEs’ finite rupture extents create differences in the apparent source time functions recorded at various stations. We use synthetic ruptures to determine how much the LFE signals recorded at each station would be modified by spatial variations of the source-station travel time within the rupture area given various possible rupture diameters, and then compare those synthetics with the data. Our synthetics show that the methodology can identify inter-station variations created by heterogeneous slip distributions or complex rupture edges, and thus lets us estimate LFE rupture extents for unilateral or bilateral ruptures. To obtain robust estimates of the sources’ similarity across stations, we stack signals from thousands of LFEs, using an empirical Green’s function approach to isolate the LFEs’ apparent source time functions from the path effects. Our analysis of LFEs in Parkfield implies that LFEs’ apparent source time functions are similar across stations at frequencies up to 8 to 16 Hz, depending on the family. The inter-station coherence observed at these relatively high frequencies, or short wavelengths (down to 0.2 to 0.5 km), suggest that LFEs in each of the 7 families examined occur on asperities. They are clustered in patches with sub-1-km diameters. The individual LFEs’ rupture diameters are estimated to be smaller than 1.1 km for all families, and smaller than 0.5 km and 1 km for the two shallowest families, which were previously found to have 0.2-s durations. Coupling the diameters with the durations suggests that it is possible to model these 〈span〉MW〈/span〉 1-2 LFEs with earthquake-like rupture speeds: around 70% of the shear wave speed. However, that rupture speed matches the data only at the edge of our uncertainty estimates for the family with highest coherence. The data for that family are better matched if LFEs have rupture velocities smaller than 40% of the shear wave speed, or if LFEs have different rupture dynamics. They could have long rise times, contain composite sub-ruptures, or have slip distributions that persist from event to event.〈/span〉

Beschreibung: 〈span〉〈div〉Summary〈/div〉On June 24, 2015, a 230,000 cubic metre landslide slid into the triangle bayou at the intersection of the Yangtze and Daning Rivers and generated a river tsunami that ran up 6.2 metres on the opposite shoreline at Wushan town. The slope failure and resulting waves killed two people and damaged many shipping facilities. Based on field surveys and eyewitness observations, we apply the ‘Tsunami Squares’ method to model the Hongyanzi landslide and its generated waves. Landslide simulations indicate a maximum impact velocity of ∼16 m/s that matches well with an eyewitness video. The computed post-slide mass stopped on the near riverbed with a shape fitting the observed geological profile. Tsunami simulations reveal a large region of wave impacts that coincide with the observed runup heights. The successful reproduction of the dynamics of this landslide-generated river tsunami emphasizes the capacity and efficiency of Tsunami Squares modeling in emergency reaction and risk assessment.〈/span〉

Beschreibung: 〈span〉〈div〉Summary〈/div〉Improvement of global 3D Earth density and velocity models is based in part on measurements of Earth’s normal mode eigenfrequencies and splitting function coefficients. Despite many methods developed inconsistency in measurements still exists and it is difficult to understand which results are more precise, that is which methods introduce less systematic biases in the measurements. Therefore, the main goal of this study is to test the performances of typically used techniques in low-frequency normal mode studies: the optimal sequence estimation stacking method and the autoregressive method in the frequency domain, where validation of the estimates is performed with the phasor walkout method. Motivations for their utilization are their easy and fast implementation and their accurate performances when it comes to eigenfrequency estimates. For this purpose, we first perform the analysis with synthetic seismograms in order to evaluate how the station distributions and noise levels impact the estimates of eigenfrequencies and structure coefficients. Synthetic seismograms are calculated for a 3D realistic Earth model, which includes Earth’s rotation as well as ellipticity and other lateral heterogeneities. They were computed by means of normal mode summation and a perturbation theory for modes up to 1 mHz. The three methods above are also applied to long-period seismometer and superconducting gravimeter data recorded after six earthquakes of magnitude greater than 8.3. Finally, our study shows that the optimal sequence estimation is sensitive to the station distribution under the noise influence, while the autoregressive method for frequency estimation gives us reasonable estimates within the estimated error bars. Moreover, we present new estimates of eigenfrequencies and Q-factors for 〈sub〉0〈/sub〉S〈sub〉2, 0〈/sub〉S〈sub〉3, 2〈/sub〉S〈sub〉1〈/sub〉 and 〈sub〉3〈/sub〉S〈sub〉1〈/sub〉 multiplets. A new value for the 〈span〉c〈/span〉〈sub〉20〈/sub〉 structure coefficient of 〈sub〉0〈/sub〉S〈sub〉2〈/sub〉 multiplet −0.7233 ± 0.0623 μHz is obtained.〈/span〉

Beschreibung: 〈span〉〈div〉Summary〈/div〉We propose a novel approach to compute the gravity field due to density anomaly in both outside and inside of the solid Earth with high accuracy and efficiency. The high accuracy comes from the direct employment of the analytic gravitation solution between any point on a two-dimensional (2D) plane in the horizontal direction and individual mass cubes. The high computational efficiency comes from two aspects: 1, the application of the highly efficient 2D discrete convolution algorithm; and 2, a newly developed algorithm for the optimized computation of the weight coefficient matrix. Numerical examples for applying to compare with analytical solutions demonstrated its excellent accuracy. Comparison with other state-of-the-art gravity modeling algorithm has proved that this algorithm has superior performance in both accuracy and efficiency. Application to analyze real topography demonstrated the practicality. This algorithm will be an attractive candidate for carrying out the forward modeling step in geophysical inversion problems with the claimed and proved advantages.〈/span〉

Beschreibung: 〈span〉〈div〉Summary〈/div〉We assess the stress field that drives lithospheric deformation, focusing on the active New Zealand plate boundary zone between the Australian and Pacific plates. Here there is a rich database for the horizontal velocity field and crustal structure, used to derive gravitational potential energy (GPE). We solve the stress balance equations, in the context of a thin sheet model of a viscously deforming lithosphere, for characteristic deviatoric stresses and viscosities, defined as the integral of these with depth (divided by the layer thickness), using the stress method of Flesch 〈span〉et al.〈/span〉 (2001), where the input parameters are the fields of strain rate and GPE that apply to the sheet. Synthetic tests show that the stress method is able to resolve the stress field to high (20 per cent) levels of noise in the input strain rates, and the mean stress is a very robust feature of the inversions, regardless of noise levels. We invert for the stress and viscosity fields in New Zealand, calculating the field of GPE from topography/bathymetry and crustal data, and the field of strain rates from either a long term (multi millennial) velocity field inferred from the rate and pattern of Quaternary faulting, or a short term (decadal) velocity field directly observed with decadal Global Positioning System (GPS) measurements. In addition, we consider the effect of shear stresses on the subducted plate interface along the Hikurangi Margin (5–15 MPa), or regionally. We explore the effect of GPE on the inversion results by calculating these for a range of deforming layers (ie thin sheet with 35–150 km thickness), in effect sampling the lithospheric strength in the crust and mantle. The results show that the derived stress magnitudes (square root second invariant of the stress deviator) are in the range 0–35 MPa, with mean values of 13 ± 1 MPa for all models, comparable to typical earthquake stress drops. Gravitationally induced stresses account for approximately half of the full deviatoric stress. Effective characteristic viscosities are 0.5 - 5 × 10〈sup〉21〈/sup〉 Pas in the deforming zone, with an approximate inverse relation between strain rate and viscosity, most likely controlled by thermal structure and/or lithology.〈/span〉

Beschreibung: 〈span〉〈div〉Summary〈/div〉I present a source-independent fracture imaging method to use passive seismic data for mapping subwavelength natural fractures. Unlike conventional source-dependent imaging that often adopts reflection-type seismic imaging with known source that is not available in passive seismic surveys, the proposed fracture imaging approach relies on the transmission and diffraction data without the need for source information. I assume that passive seismic data can be decomposed into two types of data: primary transmission wave data and diffraction (coda) wave data. The imaging formula states that primary waves should coincide with coda waves at scatterer points at the time of scattering. Instead of generating source wavefields in the conventional imaging method, the proposed method only need to propagate transmission wave data and diffraction wave data from the receiver arrays and apply an imaging condition to produce an image of fractures. This imaging procedure can be used for processing P-wave or S-wave. In synthetic examples, I evaluate the proposed method in several aspects: inaccurate source location, inaccurate velocity model, sparse receivers and irregular receiver spacing, elastic data, and joint surface and borehole acquisitions. I found that the proposed approach performed well (or even better) comparable to source-dependent fracture imaging when assuming exact source information is known. With perturbed source locations with random shifts (e.g., estimated source location with errors), however, fractures were missing in the source-dependent fracture imaging results but the proposed approach was not influenced. In the presence of velocity errors and sparse and irregular receiver spacing, the proposed method produces better fracture images than the source-dependent imaging results.〈/span〉

Beschreibung: 〈span〉〈div〉SUMMARY〈/div〉We introduce a formalism for estimating local spatial averages of the core-mantle boundary (CMB) radial magnetic field and its time derivatives, based on magnetic field observations collected by low-Earth-orbit satellites. This provides a useful alternative to conventional core field modelling based on global spherical harmonic basis functions, where noise in the polar regions maps into all harmonics, and model regularization and spectral truncation are required. A powerful perspective offered by the proposed technique is formal appraisal of the spatial resolution and variance of the resulting field averages. We use the Green’s functions for the Neumann boundary value problem to link the satellite observations to the radial magnetic field on the CMB and estimate field averages using a modified Backus-Gilbert inversion approach. Our approach builds on the Subtractive Optimally Localized Averages (SOLA) method developed in helioseismology, that seeks averaging kernels as close as possible to a chosen target kernel. We are able to account for both internal and external field sources and can easily incorporate data error covariance information, for example describing along-track serial error correlation. As a proof of concept we present a global map collecting local estimates of the radial main field (MF) constructed on a grid at the CMB with one degree spacing in latitude and longitude, derived from one month of three component vector magnetic field data collected by the 〈span〉Swarm〈/span〉 satellite trio, using data from dark and geomagnetically quiet times. Using sums and differences of the field components taken along track and in the east-west direction we obtain estimates with spatial resolution kernel widths varying between 18 and 54 degrees depending on the latitude, and a standard deviation of approximately 10μT (i.e. 5% of the mean CMB field amplitude). The morphology of our CMB radial field map agrees well with results from conventional spherical harmonic field models. In a second application, we determine local estimates of the average rate of change, or secular variation (SV), of the radial field at the CMB, initially considering two year time windows, and performing the analysis on data collected by either the 〈span〉Swarm〈/span〉 or CHAMP satellites. We obtain stable local estimates of the SV at the CMB, and present maps of estimates with averaging kernel widths of approximately 42, 33 and 30 degrees on the equator, with corresponding standard derivations of 0.25, 2.5 and 5 μT/yr. By subtracting SV estimates constructed at different epochs we are able to calculate the local aggregated secular acceleration (SA) and to study its time changes. Differencing SV estimates 2 years apart, and considering an averaging kernel width of 42 degrees on the equator, we obtain SA maps very similar to those found in the CHAOS-6-x7 field model truncated at SH degree 10. Using our approach we are able to directly control the width of the spatial averaging kernel and the length of the time window, enabling us to directly study the robustness of the inferred SA. Pushing to higher resolution in time, considering one year differences of SV estimates constructed using one year windows, we are able to track the evolution of coherent SA structures in time-longitude plots at the equator. At 25° W in mid 2007 we find a distinctive SA ’cross-over’ event, with strong, oppositely signed and adjacent, SA features rapidly changing sign within a year. Our method is well suited for studying such spatio-temporally localized SA events at high resolution; there will be further opportunities for such investigations as the time series of data provided by the 〈span〉Swarm〈/span〉 mission lengthens.〈/span〉

Beschreibung: 〈span〉〈div〉Summary〈/div〉The evolution of the ratio between P- and S-wave velocity (〈span〉V〈/span〉〈sub〉P〈/sub〉/〈span〉V〈/span〉〈sub〉S〈/sub〉) with increasing fluid-saturated porosity is computed for isotropic rocks containing spheroidal pores. The ratio 〈span〉V〈/span〉〈sub〉P〈/sub〉/〈span〉V〈/span〉〈sub〉S〈/sub〉 is shown to either decrease or increase with increasing porosity, depending on the aspect ratio α of the pores, fluid to solid bulk modulus ratio ζ, and Poisson’s ratio ν〈sub〉0〈/sub〉 of the solid constituents of the rock. A critical initial Poisson’s ratio ν〈sub〉0, crit〈/sub〉 is computed, separating cases where 〈span〉V〈/span〉〈sub〉P〈/sub〉/〈span〉V〈/span〉〈sub〉S〈/sub〉 increases (if ν〈sub〉0〈/sub〉 〈 ν〈sub〉0, crit〈/sub〉) or 〈span〉decreases〈/span〉 (if ν〈sub〉0〈/sub〉 〉 ν〈sub〉0, crit〈/sub〉) with increasing porosity. For thin cracks and highly compressible fluids, ν〈sub〉0, crit〈/sub〉 is approximated by 0.157 ζ/α, whereas for spherical pores ν〈sub〉0, crit〈/sub〉 is given by 0.2 + 0.8ζ. When ν〈sub〉0〈/sub〉 is close to ν〈sub〉0, crit〈/sub〉, the evolution of 〈span〉V〈/span〉〈sub〉P〈/sub〉/〈span〉V〈/span〉〈sub〉S〈/sub〉 with increasing fluid-saturated porosity is near neutral and depends on subtle changes in pore shape and fluid properties. This regime is found to be relevant to partially dehydrated serpentinites in subduction zones (porosity of aspect ratio near 0.1 and ζ in the range 0.01–0.1), and makes detection of these rocks and possibly elevated fluid pressures difficult from 〈span〉V〈/span〉〈sub〉P〈/sub〉/〈span〉V〈/span〉〈sub〉S〈/sub〉 only.〈/span〉

Beschreibung: 〈span〉〈div〉Summary〈/div〉Secondary microseismic sources emit seismic waves over long time spans. Reoccurring signals with similar slowness and frequency therefore arrive at seismic arrays. Blind source separation techniques can be used to identify and isolate such reoccurring signals from other signals and from diffuse seismic noise. Along these lines, we use non-negative matrix factorization as blind source separation technique to decompose continuous seismic array records. We model the recorded energy as a mixture of a few components with static slowness-frequency and time dependent amplitudes. Components and amplitudes are fitted to optimally explain the recorded seismic energy over time. These components represent secondary microseismic signals with quasi-static slowness-frequency vector and fluctuating amplitude. Each fitted component reveals the geographical origin (through the slowness-frequency vector) and time evolution of an active secondary microseism with high precision because it is separated from other signals and diffuse seismic noise. Furthermore, relative travel times can be automatically extracted for the signals that correspond to a specific component that can potentially be used in tomographic studies. We show two examples of seismic signals that were extracted with this technique, one focusing on P-waves from the typhoons Goni and Atsani, and another showing secondary microseism PKP signals from typhoon Glenda.〈/span〉

Beschreibung: 〈span〉〈div〉Summary〈/div〉Locating microseismic events is essential for many areas of seismology including volcano and earthquake monitoring and reservoir engineering. Due to the large number of microseismic events in these settings, an automated seismic location method is required to perform real time seismic monitoring. The measurement environment requires a precise and noise-resistant event location method for seismic monitoring. In this paper, we apply Multichannel Coherency Migration (MCM) to automatically locate microseismic events of induced and volcano-tectonic seismicity using sparse and irregular monitoring arrays. Compared to other migration-based methods, in spite of the often sparse and irregular distribution of the monitoring arrays, the MCM can show better location performance and obtain more consistent location results with the catalogue obtained by manual picking. Our MCM method successfully locates many triggered volcano-tectonic events with local magnitude smaller that 0, which demonstrates its applicability on locating very small earthquakes. Our synthetic event location example at a carbon capture and storage site shows that continuous and coherent drilling noise in industrial settings will pose great challenges for source imaging. However, automatic quality control techniques including filtering in the frequency domain and weighting are used to automatically select high quality data, and can thus effectively reduce the effects of continuous drilling noise and improve source imaging quality. The location performance of the MCM method for synthetic and real microseismic datasets demonstrates that the MCM method can perform as a reliable and automatic seismic waveform analysis tool to locate microseismic events.〈/span〉

Beschreibung: 〈span〉〈div〉Summary〈/div〉The unification of local height systems has been a classical geodetic problem for a long time, the main challenges of which are the estimation of offsets between different height systems and the correction of tilts along the levelling lines. It has been proposed to address these challenges with clock networks. The latest generation of optical clocks as well as the dedicated frequency links, e.g., optical fibres, are now approaching to deliver the comparison of frequencies at the level of 1.0 × 10〈sup〉−18〈/sup〉. It corresponds to an accuracy of about 1.0 cm in height difference. Clock networks can thus serve as a powerful tool to connect local height systems. To verify the idea, we carried out simulations using the EUVN/2000 (European Unified Vertical Network) as a priori input. Four local height systems were simulated from the EUVN/2000 by introducing individual offsets and tilts, and were re-unified by using measurements in clock networks. The results demonstrate the great potential of clock networks for height system unification. In case that the offsets between different height systems and tilts along national levelling lines in both longitudinal and latitudinal directions are considered, three or four clocks measurements for each local region are sufficient for the unification. These clocks are to be interconnected and should be properly arranged so that they can sense the levelling tilts where necessary. Our results also indicate that even clocks with one magnitude poorer accuracy than the desired ones can still unify the height systems to some extent, but it may cause a shift for the re-unified system.〈/span〉

Beschreibung: 〈span〉〈div〉Summary〈/div〉We investigate slip-distribution models of the 2011 Tohoku earthquake, with a particular focus on diffracted tsunamis and uplift-induced waves along the back-arc region of the Japanese Island Arc.The 2011 Tohoku earthquake produced a large amplitude tsunami that diffracted around Kyushu Island before reaching Korea. At the same time, this earthquake co-seismically induced short-period small-amplitude sea waves in the East Sea. We performed tsunami simulations using seven fault models of the Tohoku earthquake to examine whether the models can accurately reproduce the observed waveforms in the open sea of the western Pacific Ocean, the South Sea of Korea, and the coast of the East Sea. For each fault model, we investigate tsunami features due to geomorphological characteristics of the Korean Peninsula in the Korea offshore. To determine which slip distribution model shows a good performance in the tsunami simulations, we set three criteria; the delay time between observations and synthetic waveforms, the normalized mean residual, and the normalized RMS misfit. Depending on the study region, all models show varying degrees of accuracy. The fault dimensions and the amount of slip have a larger effect on the RMS misfit then the slip distribution patterns of the fault models for observations along the Korean coast and the western coast of Japan.〈/span〉

Beschreibung: 〈span〉〈div〉SUMMARY〈/div〉The study of radiation patterns in poroelastic media allows us to visually explore the possibility of reconstructing crucial model properties describing reservoir rocks, and to examine the coupling effects between different parameters during full-waveform inversion (FWI). In this paper, we derive analytical formulae for the radiation patterns of single parameter perturbations in fluid-saturated porous media by deriving scattered wavefields based on plane-wave theory and the far-field approximation. We illustrate these scattered wavefields via their radiation patterns expressed as a function of the angle between the incident and scattered waves. To simplify the algebra, we consider poroelastic waves at seismic (low) frequencies, where the fast compressional wave and shear wave are propagating modes but the slow compressional wave is severely dispersive. To verify our derivation of the analytical radiation patterns, we also compute them numerically by perturbing one parameter at a single point, keeping the other parameters fixed at their background values. We find that all the analytical radiation patterns match the wavefronts of the numerically computed scattered wavefields, well indicating that our derivations are correct. Parameters such as the solid density, fluid density, viscosity of the fluid, and intrinsic permeability, have similar radiation patterns and thus show strong coupling effects. Therefore, we anticipate difficulties in recovering these parameters in a multi-parameter FWI procedure. In an attempt to mitigate these trade-offs, we analyze different parameterizations which result in different radiation patterns. As the patterns that we observe are similar to those of the elastic case, we anticipate that parameter separation might be easier when inverting for velocities than for moduli.〈/span〉

Beschreibung: 〈span〉〈div〉Abstract〈/div〉In the Sea of Marmara, areas of gas seepage or cold seeps are tightly related to the faults system and understanding the spatial and temporal dynamics in gas-related processes is crucial for geohazard mitigation. Although acoustic surveys proved to be efficient in detecting and locating cold seeps, temporal variability or trends in the gas-related processes are still poorly understood. Two arrays of 10 ocean bottom seismometers were deployed in the western part of the Sea of Marmara in 2011 and 2014, respectively. In addition to the local seismic events, the instruments recorded a large number of short duration events and long-lasting tremors. Short duration events are impulsive signals with duration 〈 1 s, amplitude well above the noise level and a frequency spectrum with one or two narrow peaks. They are not correlated from one site to another, suggesting a very local source. Tremors consist of sequences of clustered impulsive signals lasting for minutes to more than an hour with a multi-peak frequency spectrum. Based on evidence of known seepage and by analogy with volcanic and hydrothermal models, we suggest that short duration events and tremors are associated with gas migration and seepage. There is a relationship between tremors associated with gas emission and the local seismicity, although not systematic. Rather than triggering gas migration out of the seabed, locally strong earthquakes act as catalysts when gas is already present or gas emission is already initiated.〈/span〉

Beschreibung: 〈span〉〈div〉SUMMARY〈/div〉Long-period Rayleigh wave Horizontal to Vertical amplitude (H/V) ratios at a station provide information about local earth structure that is complementary to phase velocity. However, a number of studies (Ferreira and Woodhouse, 2007; Lin 〈span〉et al.〈/span〉, 2012) have observed that significant scatter appears in these measurements making it difficult to use H/V ratio measurements to resolve earth structure. Some of the scatter in these measurements has been attributed to local geological structure while some has remained unaccounted for. Most Global Seismographic Network (GSN) stations contain two nearby high-quality broadband seismometers (e.g. in the same vault, but on different piers or in different boreholes). For each broadband sensor in the IRIS/USGS component of the GSN, we estimate H/V ratios of fundamental mode Rayleigh waves using M 〉 6.5 earthquakes from 2001 to 2018 (around 19,000 measurements). We compute these ratios at a number of discrete periods (25, 50, 75, 100, and 150 s) and find that for well-isolated Rayleigh waves (windows where the correlation coefficients between radial and the phase-shifted vertical components are greater than 0.9) significant scatter in H/V ratios occurs between co-located sensors (greater than 25 per cent at 100 s period). This suggests the scatter in H/V ratio measurements can be at least partially attributed to extremely local phenomena such as sensor emplacement in the vault. We also find that H/V ratios can vary as a function of event back-azimuth, indicating that care must be taken when computing average ratios for a station, as a large number of events from a given region could bias H/V ratio measurements at a station.〈/span〉

Beschreibung: 〈span〉〈div〉SUMMARY〈/div〉We analyse noise characteristics of spatial and temporal correlation of 260 continuous GPS sites from Crustal Movement Observation Network of China (CMONOC). These data sets were mainly collected between 2010 and 2016, with an average of 6 yr of position time-series. In the functional analysis, a clear regional dependence of seasonal movements has been observed and other significant periodic signals are detected nearby the GPS draconitic period and its harmonics. The distribution of these periodic signals shows a spatial correlation, along with non-negligible local inconsistencies. In the stochastic analysis, impacts of the periodicities on the noise assessment have been investigated and Maximum likelihood estimation is used to study noise properties for deseasonalized residual time-series having the seasonal signals removed and filtered residual time-series having other periodic signals removed further. We demonstrate that for both solutions, the flicker noise is thought to be the dominant time-correlated noise and velocity uncertainties may be underestimated 8–10 times if assuming a pure white noise. Ignoring the periodicities could bias the estimation of noise amplitude, spectral index and velocity uncertainty. After removing the periodic signals, the median flicker noise magnitude shows an average of ∼10 per cent reductions and the noise process shifts closer towards white noise. A correlation between the index variations and RMS variations has been observed, indicating that the index varies more significantly for sites with more periodic signals removed. Besides, the spectral index in the vertical component has a better spatial correlation than that in the horizontal and the spatial distribution of the index of deseasonalized solutions seems to correlate well with the amplitudes of seasonal signals, probably implying common sources of spatial variations of these characteristics. Furthermore, analyses of intersite correlations indicate that the correlation induced from the periodic signals displays a similar pattern to the deseasonalized solutions, confirming that the period pattern is spatially correlated and can induce time-series correlations. Finally, the stochastic processes of the common mode noise are predominantly featured by spatially correlated flicker and white noise over a wide range, consistent with previous results.〈/span〉

Beschreibung: 〈span〉〈div〉SUMMARY〈/div〉As temperature increases with depth and the creep resistance of rock decreases exponentially, a high-viscosity sub-lithospheric layer, just beneath the ‘elastic’ lithosphere is expected to exist. Depending on the temperature profile, a low-viscosity asthenosphere may also exist if the temperature deeper down gets high enough. Since the temperature profile is expected to change laterally – especially from below the oceans to cratonic areas underneath continents, rock properties of the lithosphere, high-viscosity sub-lithosphere and low-viscosity asthenosphere are expected to change laterally. Our aim is to constrain sub-lithospheric properties (depth, thickness and viscosity), lateral lithospheric thickness variations and asthenospheric properties using observed GIA data. A Coupled Laplace-Finite Element Method is used to compute gravitationally self-consistent sea level with time-dependent coastline and rotational feedback in addition to changes in deformation, gravity and the state of stress. We start with the VM5a-ICE-6G_C model combination and then modify the lithospheric, sub-lithospheric and asthenospheric properties (including lateral thickness variation) while keeping the mantle viscosities the same as VM5a. Through this study, we confirm that the sub-lithospheric and asthenospheric properties can significantly affect the predicted global relative sea level (RSL), present-day gravity rate-of-change (g-dot) and uplift rate (u-dot) in Laurentia and Fennoscandia. In addition, incorporating the elastic lithosphere with lateral thickness variation, sub-lithosphere and asthenosphere can improve the fit to global RSL, but the predicted peak values of g-dot and u-dot in Laurentia may decrease slightly but not significant enough to affect the fit to the observed data. Our results prefer an elastic lithosphere that has maximum thickness of 140 km under continental cratons but reduces to 60 km underneath the oceans. The results preferred depth of the asthenospheric bottom is around 190–200 km with asthenospheric viscosity around 10〈sup〉20〈/sup〉Pa s. Finally, we show that the best laterally heterogeneous mantle model we found in previous publication when combined with the lithosphere with lateral thickness variaion gives the best fit to global RSL and peak g-dot and u-dot in Laurentia simultaneously.〈/span〉

Beschreibung: 〈span〉〈div〉SUMMARY〈/div〉We implement an analytical model based on flexural deflection of a thin elastic disc to investigate the magnitude of lithospheric decompression caused by deglaciation at upper crustal magmatic reservoirs. Considering a published numerical climate model describing the space–time evolution of deglaciation after the Last Glacial Maximum (LGM) along the Southern Volcanic Zone (SVZ) of the Andes, we demonstrate that changes in pressure at upper crustal levels (〈10 km depth) at the scale of several hundred years are of the order of 10–100 MPa. Total decompression and decompression rate (300–150 kPa yr〈sup〉−1〈/sup〉) are 1–2 orders of magnitude larger than values previously estimated by other authors who assume that glacial loads are supported by an elastic half-space, that is, of infinite elastic thickness. The large decompression caused by flexural unbending of an elastic plate of finite thickness as assumed here can easily surpass the tensile strength of rocks (5–20 MPa), creating adequate conditions for failure of the reservoir walls, dike propagation inside and outside the reservoir and the eventual collapse of the reservoir accompanying an explosive eruption. We apply our results to the analysis of post-glacial eruptions of SVZ volcanoes, which erupted large volumes (〉10 km〈sup〉3〈/sup〉) of mafic ignimbrites hundreds to thousands of years after deglaciation onset. We show that this time lag is necessary to achieve a decompression of several tens of megapascals at depths of several kilometres that are consistent with the location of magmatic reservoirs as estimated by independent petrologic, seismic and/or geodetic studies. Moreover the northward increase of this time lag is in agreement with a smaller size of the Andean ice cap in the north than in the south during the LGM. For wet, volatile-rich magmas typical of subduction zones, the effect of large decompression at upper crustal reservoirs caused by flexural unbending of the lithosphere after deglaciation could play a major role in promoting large explosive eruptions through devolatization of the magma, during past deglaciation events as demonstrated here for the LGM along the SVZ and current accelerated ice retreat caused by climate change over large segments of subduction-related arcs at higher latitudes.〈/span〉

Beschreibung: 〈span〉〈div〉SUMMARY〈/div〉The unification of local height systems has been a classical geodetic problem for a long time, the main challenges of which are the estimation of offsets between different height systems and the correction of tilts along the levelling lines. It has been proposed to address these challenges with clock networks. The latest generation of optical clocks as well as the dedicated frequency links, for example optical fibres, are now approaching to deliver the comparison of frequencies at the level of 1.0 × 10〈sup〉−18〈/sup〉. It corresponds to an accuracy of about 1.0 cm in height difference. Clock networks can thus serve as a powerful tool to connect local height systems. To verify the idea, we carried out simulations using the EUVN/2000 (European Unified Vertical Network) as 〈span〉apriori〈/span〉 input. Four local height systems were simulated from the EUVN/2000 by introducing individual offsets and tilts, and were reunified by using measurements in clock networks. The results demonstrate the great potential of clock networks for height system unification. In case that the offsets between different height systems and tilts along national levelling lines in both longitudinal and latitudinal directions are considered, three or four clocks measurements for each local region are sufficient for the unification. These clocks are to be interconnected and should be properly arranged so that they can sense the levelling tilts where necessary. Our results also indicate that even clocks with one magnitude poorer accuracy than the desired ones can still unify the height systems to some extent, but it may cause a shift for the reunified system.〈/span〉