A family-level macroinvertebrate biotic index for ecological assessment of lakes in Yunnan,China

Persistent economic growth in Chinese southwestern Yunnan Province is setting aquatic ecosystems in its plateau lakes under enormous pressure. While several different systems have previously been used to study these lakes, no existing methodology adequately measures both the chemical and biotic parameters of these water bodies. Here, we present a novel Biotic Monitoring Yunnan Lakes (BMYL) index that provides a general assessment tool for ecological deterioration that is caused by organic pollution. Principal Component Analysis is used to analyze the occurrence of families of macroinvertebrates and chemical properties of the lakes. In brief, families of macroinvertebrates were given a score from 10 to 1 based on sensitivity to organic enrichment and eutrophication. Sampling at each lake yielded an Average Score per Lake (ASPL) which is calculated by dividing the total BMYL by the total number of scoring families. High ASPL values characterize a biologically intact lake containing relatively large numbers of high scoring taxa, while lower ASPL values denote a polluted lake that does not support many high scoring taxa. The results of the BMYL show a notably more accurate characterization of the long-term health of concerned aquatic ecosystems than studies that use abundance levels of species or a simple analysis of chemical parameter.     

Evaluation of the 15 July 2009 Fiordland,New Zealand Tsunami in the Source Region

We undertake detailed near-field numerical modelling of the tsunami generated by the 15 July 2009 earthquake (Mw 7.8) in Fiordland, New Zealand. High resolution bathymetry and topography data at Breaksea and Dusky Sounds, and Chalky and Preservation Inlets are derived mostly from digitised New Zealand nautical charts, Shuttle Radar Topographic Mission (SRTM) 3 arc-second data, and General Bathymetric Chart of the Ocean (GEBCO) 30 s data. A combination of continuous and campaign Global Positioning System (GPS), satellite radar (ALOS/PALSAR InSAR images) and seismology data are used to constrain the seafloor deformation for the initial tsunami condition. This source model, derived independently of DART observations, provides an excellent fit to observed tsunami elevations recorded by DART buoy 55015. The model results in the near field show maximum tsunami elevations in the range 0.5–2.0 m inside the sounds and inlets with maximum flow speeds of 3.0 m/s. Along the open coast, maximum tsunami elevations reach 2.0 m. The high flow speeds through the inlets may change the inlet stratifications and water mass inside the sounds. Media reports and field reconnaissance data show some tsunami evidence at Cormorant Cove, Duck and Goose Coves, and Passage Point.     

Computing the Sensitivity Kernels for 2.5-D Seismic Waveform Inversion in Heterogeneous,Anisotropic Media

2.5-D modeling and inversion techniques are much closer to reality than the simple and traditional 2-D seismic wave modeling and inversion. The sensitivity kernels required in full waveform seismic tomographic inversion are the Fréchet derivatives of the displacement vector with respect to the independent anisotropic model parameters of the subsurface. They give the sensitivity of the seismograms to changes in the model parameters. This paper applies two methods, called ‘the perturbation method’ and ‘the matrix method’, to derive the sensitivity kernels for 2.5-D seismic waveform inversion. We show that the two methods yield the same explicit expressions for the Fréchet derivatives using a constant-block model parameterization, and are available for both the line-source (2-D) and the point-source (2.5-D) cases. The method involves two Green’s function vectors and their gradients, as well as the derivatives of the elastic modulus tensor with respect to the independent model parameters. The two Green’s function vectors are the responses of the displacement vector to the two directed unit vectors located at the source and geophone positions, respectively; they can be generally obtained by numerical methods. The gradients of the Green’s function vectors may be approximated in the same manner as the differential computations in the forward modeling. The derivatives of the elastic modulus tensor with respect to the independent model parameters can be obtained analytically, dependent on the class of medium anisotropy. Explicit expressions are given for two special cases—isotropic and tilted transversely isotropic (TTI) media. Numerical examples are given for the latter case, which involves five independent elastic moduli (or Thomsen parameters) plus one angle defining the symmetry axis.     

Sensitivity of Typhoon Vamei (2001) Simulation to Planetary Boundary Layer Parameterization Using PSU/NCAR MM5

This paper investigates the sensitivity of the numerical simulations of a near equatorial Typhoon Vamei (2001) to various planetary boundary layer (PBL) parameterization schemes in the Pennsylvania State University (PSU)/National Centre for Atmospheric Research (NCAR) non-hydrostatic mesoscale model (MM5). The numerical simulations are conducted on two domains at 45 and 15 km grids nested in a one-way fashion. Four different PBL parameterization schemes including the Blackadar (BLK) scheme, the Burk–Thompson (BURKT) scheme, the NCEP Eta model scheme and the NCEP medium range forecast (MRF) model scheme are investigated. Results indicate that the intensity and propagation track of the simulated near equatorial typhoon system is not very sensitive to the different PBL treatments. The simulated minimum central pressures and the maximum surface wind speeds differ by only 5–6 hPa and 6–8 ms⁻¹, respectively. Larger variations between the simulations occur during the weakening phase of the typhoon system. While all schemes simulated the typhoon with reasonable accuracy, the ETA scheme produces the strongest storm intensity with the largest heat exchanges over the marine environment and the highest warm moisture air content in the PBL around the core of the storm.     

LES of turbulent surface shear stress and pressure-gradient-driven flow on shallow continental shelves

Turbulent shear flows on shallow continental shelves (here shallow means that the interaction with the solid, no-slip bottom is important) are of great importance because tide- and wind-driven flows on the shelf are drivers of the transfer of momentum, heat, and mass (gas) across the air–sea interface. These turbulent flows play an important role because vertical mixing and current are vectors for the transport of sediment and bioactive material on continental shelves. Understanding the dynamics of this class of flows presents complications because of the presence of a free surface and also because the flow can be driven by a pressure gradient (a tidal current), a stress at the free surface (a wind-driven current), or a combination of both. In addition, the flow can be modified by the presence of a wave field that can induce Langmuir circulation (Langmuir, Science 87:119–123, 1938). Large eddy simulation is used to quantify the effects of pressure gradient and wind shear on the distinctive structures of the turbulent flow. From these computations, an understanding of the physics governing the turbulence of pressure-driven and wind-driven flows, how they can interact in a normal or a tangential direction, and the effect of wave forcing on these flows is obtained.     

FerryBox - Application of continuous water quality observations along transects in the North Sea

In 2002, the Helmholtz-Zentrum Geesthacht, Germany, started to use FerryBox-automated monitoring systems on Ships of Opportunity to continuously record standard oceanographic, biological and chemical in situ data in the North Sea. The present study summarises the operational experience gathered since the beginning of this deployment and reflects on the potential and limits of FerryBox systems as a monitoring tool. One part relates to the instrumental performance, constancy of shipping services, and the availability and quality of the recorded in situ data. The other considers integration of the FerryBox observations in scientific applications and routine monitoring campaigns. Examples are presented that highlight the added value of the recorded data for the study of both long- and short-term variability in water mass stability, plankton communities and surface water productivity in the North Sea. Through the assessment of technical and scientific performance, it is evident that FerryBoxes have become a valuable tool in marine research that helps to fill gaps in coastal and open ocean operational observation networks.     

The May 2005 eruption of Fernandina volcano,Galápagos: The first circumferential dike intrusion observed by GPS and InSAR

The May 2005 eruption of Fernandina volcano, Galápagos, occurred along circumferential fissures parallel to the caldera rim and fed lava flows down the steep southwestern slope of the volcano for several weeks. This was the first circumferential dike intrusion ever observed by both InSAR and GPS measurements and thus provides an opportunity to determine the subsurface geometry of these enigmatic structures that are common on Galápagos volcanoes but are rare elsewhere. Pre- and post- eruption ground deformation between 2002 and 2006 can be modeled by the inflation of two separate magma reservoirs beneath the caldera: a shallow sill at ~1 km depth and a deeper point-source at ~5 km depth, and we infer that this system also existed at the time of the 2005 eruption. The co-eruption deformation is dominated by uplift near the 2005 eruptive fissures, superimposed on a broad subsidence centered on the caldera. Modeling of the co-eruption deformation was performed by including various combinations of planar dislocations to simulate the 2005 circumferential dike intrusion. We found that a single planar dike could not match both the InSAR and GPS data. Our best-fit model includes three planar dikes connected along hinge lines to simulate a curved concave shell that is steeply dipping (~45–60°) toward the caldera at the surface and more gently dipping (~12–14°) at depth where it connects to the horizontal sub-caldera sill. The shallow sill is underlain by the deep point source. The geometry of this modeled magmatic system is consistent with the petrology of Fernandina lavas, which suggest that circumferential eruptions tap the shallowest parts of the system, whereas radial eruptions are fed from deeper levels. The recent history of eruptions at Fernandina is also consistent with the idea that circumferential and radial intrusions are sometimes in a stress-feedback relationship and alternate in time with one another.     

Shallow degassing events as a trigger for very-long-period seismicity at Kīlauea Volcano,Hawai‘i

The first eruptive activity at Kīlauea Volcano’s summit in 25 years began in March 2008 with the opening of a 35-m-wide vent in Halema‘uma‘u crater. The new activity has produced prominent very-long-period (VLP) signals corresponding with two new behaviors: episodic tremor bursts and small explosive events, both of which represent degassing events from the top of the lava column. Previous work has shown that VLP seismicity has long been present at Kīlauea’s summit, and is sourced approximately 1 km below Halema‘uma‘u. By integrating video observations, infrasound and seismic data, we show that the onset of the large VLP signals occurs within several seconds of the onset of the degassing events. This timing indicates that the VLP is caused by forces—sourced at or very near the lava free surface due to degassing—transmitted down the magma column and coupling to the surrounding rock at 1 km depth.     

Synergy and fusion of optical and synthetic aperture radar satellite data for underwater topography estimation in coastal areas

A method to obtain underwater topography for coastal areas using state-of-the-art remote sensing data and techniques worldwide is presented. The data from the new Synthetic Aperture Radar (SAR) satellite TerraSAR-X with high resolution up to 1 m are used to render the ocean waves. As bathymetry is reflected by long swell wave refraction governed by underwater structures in shallow areas, it can be derived using the dispersion relation from observed swell properties. To complete the bathymetric maps, optical satellite data of the QuickBird satellite are fused to map extreme shallow waters, e.g., in near-coast areas. The algorithms for bathymetry estimation from optical and SAR data are combined and integrated in order to cover different depth domains. Both techniques make use of different physical phenomena and mathematical treatment. The optical methods based on sunlight reflection analysis provide depths in shallow water up to 20 m in preferably calm weather conditions. The depth estimation from SAR is based on the observation of long waves and covers the areas between about 70- and 10-m water depths depending on sea state and acquisition quality. The depths in the range of 20 m up to 10 m represent the domain where the synergy of data from both sources arises. Thus, the results derived from SAR and optical sensors complement each other. In this study, a bathymetry map near Rottnest Island, Australia, is derived. QuickBird satellite optical data and radar data from TerraSAR-X have been used. The depths estimated are aligned on two different grids. The first one is a uniform rectangular mesh with a horizontal resolution of 150 m, which corresponds to an average swell wavelength observed in the 10 × 10-km SAR image acquired. The second mesh has a resolution of 150 m for depths up to 20 m (deeper domain covered by SAR-based technique) and 2.4 m resolution for the shallow domain imaged by an optical sensor. This new technique provides a platform for mapping of coastal bathymetry over a broad area on a scale that is relevant to marine planners, managers, and offshore industry.