Internal solitary wave-induced flow over a corrugated bed

A combined numerical and experimental study of the propagation of an internal solitary wave (ISW) over a corrugated bed is presented, in which the amplitude and the wavelength of the corrugated bed, together with the wave amplitude and wave speed of the ISW, have been varied parametrically. Both ISWs of elevation and depression have been considered. The wave-induced currents over the corrugated bed cause flow separation at the apex of the corrugations and a sequence of lee vortices forms as a result. These vortices develop fully after the main wave has passed over the topographic feature, resulting in deformation of the overlying pycnocline and, in some instances, significant vertical mixing. It is found that the intensity of the vortex formation is dependent on both the amplitude and wavelength of the bottom topography. In the case of an ISW of depression, the generation of vertically (upward)-propagating vortices is seen to result in entrainment of fluid from a bottom boundary jet (Carr and Davies, Phys Fluids 18:016601, 2006), while, in the elevation case, a second mechanism is present to induce significant turbulent mixing in the water column. It occurs when the bottom corrugations reach into, or are very near, the pycnocline at rest. Large waves of elevation that are stable on approach to the corrugations exhibit evidence of a spatio-temporally developing shear instability as they interact with the bottom corrugation. The shear instability takes the form of billows that have a vertical extent that can reach 50% of the wave amplitude.     

Waves, circulation and vertical dependence

Longuet-Higgins and Stewart (J Fluid Mech 13:481–504, 1962; Deep-Sea Res 11:529–562, 1964) and later Phillips (1977) introduced the problem of waves incident on a beach, from deep to shallow water. From the wave energy equation and the vertically integrated continuity equation, they inferred velocities to be Stokes drift plus a return current so that the vertical integral of the combined velocities was nil. As a consequence, it can be shown that velocities of the order of Stokes drift rendered the advective term in the momentum equation negligible resulting in a simple balance between the horizontal gradients of the vertically integrated elevation and wave radiation stress terms; the latter was first derived by Longuet-Higgins and Stewart. Mellor (J Phys Oceanogr 33:1978–1989, 2003a), noting that vertically integrated continuity and momentum equations were not able to deal with three-dimensional numerical or analytical ocean models, derived a vertically dependent theory of wave–circulation interaction. It has since been partially revised and the revisions are reviewed here. The theory is comprised of the conventional, three-dimensional, continuity and momentum equations plus a vertically distributed, wave radiation stress term. When applied to the problem of waves incident on a beach with essentially zero turbulence momentum mixing, velocities are very large and the simple balance between elevation and radiation stress gradients no longer prevails. However, when turbulence mixing is reinstated, the vertically dependent radiation stresses produce vertical velocity gradients which then produce turbulent mixing; as a consequence, velocities are reduced, but are still larger by an order of magnitude compared to Stokes drift. Nevertheless, the velocity reduction is sufficient so that elevation set-down obtained from a balance between elevation gradient and radiation stress gradients is nearly coincident with that obtained by the aforementioned papers. This paper includes four appendices. The first appendix demonstrates the numerical process by which Stokes drift is excluded from the turbulence stress parameterization in the momentum equation. A second appendix determines a bottom slope criterion for the application of linear wave relations to the derivation of the wave radiation stress. The third appendix explores the possibility of generalizing results by non-dimensionalization. The final appendix applies the basic theory to a problem introduced by Bennis and Ardhuin (J Phys Oceanogr 41:2008–2012, 2011).     

Modelling sand wave migration in shallow shelf seas

Sand waves form a prominent regular pattern in the offshore seabed of sandy shallow seas. The positions of sand-wave crests and troughs slowly change in time. Sand waves are usually assumed to migrate in the direction of the residual current. This paper considers the physical mechanisms that may cause sand waves to migrate and methods to quantify the associated migration rates. We carried out a theoretical study based on the assumption that sand waves evolve as free instabilities of the system. A linear stability analysis was then performed on a 2DV morphological model describing the interaction between the vertically varying water motion and an erodible bed in a shallow sea. Here, we disrupted the basic tidal symmetry by choosing a combination of a steady current (M₀) and a sinusoidal tidal motion (M₂) as the basic flow. We allowed for two different physical mechanisms to generate the steady current: a sea surface wind stress and a pressure gradient. The results show that similar sand waves develop for both flow conditions and that these sand waves migrate slowly in the direction of the residual flow. The rates of migration and wavelengths found in this work agree with theoretical and empirical values reported in the literature.     

Towards medium-term (order of months) morphodynamic modelling of the Teign estuary, UK

The main objective of this paper is to address the principal mechanisms involved in the medium-term (order of months to years) morphodynamic evolution of estuaries through the application of a process-based numerical modelling. The Teign estuary (Teignmouth, UK) is the selected site. The system is forced by the macrotidal semi-diurnal tide in the English Channel and is perturbed to a minor extent by high river discharge events (freshets). Although waves have a definite influence on the adjacent coastal area, Wells (Teignmouth Quay Development Environmental Statement: Changes to Physical Processes. Report R.984c:140. ABP Marine Environmental Research Ltd., Southampton, 2002b) suggested that swell waves do not enter the estuary. Hence, wave effects are neglected in this study, as only tides and the river discharge are taken into account. The sediment grain size is highly variable, but mainly sandy. Within the frame of the COAST3D project (), four bathymetric surveys of the adjacent coastal area were carried out at a nearly weekly intervals. The outer estuary and the adjacent coastal area were also surveyed every 6 months as part of the COASTVIEW project (). Based on these data and on continuously measured parameters, such as water level, waves, wind and river discharge, numerical modelling of the morphodynamic processes can be tested. To replicate the morphological changes in the medium-term within a feasible simulation time, forcing conditions are reduced through the use of an input reduction method (called ensemble technique). In this study, simulations are based on the coupling between Telemac-2D and its non-cohesive sediment transport module, Sisyphe (version 5.3 for both modules). Three different sediment transport formulae were tested: (1) Engelund and Hansen (A monograph on sediment transport in alluvial streams, 3rd edn. Technological University of Denmark, Copenhagen, 1967) including the modifications proposed by Chollet and Cunge (J Hydraul Eng 17(1):1–13, 1979); (2) Bijker (Mechanics of sediment transport by the combination of waves and current. In: Design and reliability of coastal structures. 23rd international conference on Coastal Engineering, pp 147–173, 1968) and (3) Soulsby (Dynamics of Marine Sands. A manual for practical applications. HR Wallingford, Wallingford, p 142, 1997) modified version of van Rijn [J Hydraul Eng 110(10):1431–1456, 1984a, J Hydraul Eng 110(11):1613–1641, 1984b] formulation. Both a qualitative (i.e. visual comparison) and a quantitative tool [Brier Skill Score (BSS); described in Sutherland et al. in Coast Eng 51:917–939, 2004b] are applied to assess the similarity of simulations when compared to model predictions and observations. Tests confirmed the reliability and time efficiency of the ensemble technique, since it reproduced very well the results of a reference run, a computation based on the observed boundary conditions. For the spring-neap cycle modelled, the BSS was of 0.91 (a perfect modelling would have a BSS of 1), with a reduction in the simulation time on the order of 80%. For the 6-month-period simulation, results were also excellent: BSS=0.92 and a computer time reduction of 85%. In principle, this method has the advantage of being applied to any process-based numerical model.     

Wave-induced topographic formstress in baroclinic channel flow

Large-scale zonal flow driven across submarine topography establishes standing Rossby waves. In the presence of stratification, the wave pattern can be represented by barotropic and baroclinic Rossby waves of mixed planetary topographic nature, which are locked to the topography. In the balance of momentum, the wave pattern manifests itself as topographic formstress. This wave-induced formstress has the net effect of braking the flow and reducing the zonal transport. Locally, it may lead to acceleration, and the parts induced by the barotropic and baroclinic waves may have opposing effects. This flow regime occurs in the circumpolar flow around Antarctica. The different roles that the wave-induced formstress plays in homogeneous and stratified flows through a zonal channel are analyzed with the BARBI (BARotropic-Baroclinic-Interaction ocean model, Olbers and Eden, J Phys Oceanogr 33:2719–2737, 2003) model. It is used in complete form and in a low-order version to clarify the different regimes. It is shown that the barotropic formstress arises by topographic locking due to viscous friction and the baroclinic one due to eddy-induced density advection. For the sinusoidal topography used in this study, the transport obeys a law in which friction and wave-induced formstress act as additive resistances, and windstress, the effect of Ekman pumping on the density stratification, and the buoyancy forcing (diapycnal mixing of the stratified water column) of the potential energy stored in the stratification act as additive forcing functions. The dependence of the resistance on the system parameters (lateral viscosity ε, lateral diffusivity κ of eddy density advection, Rossby radius λ, and topography height δ) as well as the dependence of transport on the forcing functions are determined. While the current intensity in a channel with homogeneous density decreases from the viscous flat bottom case in an inverse quadratic law ~δ ^–2 with increasing topography height and always depends on ε, a stratified system runs into a saturated state in which the transport becomes independent of δ and ε and is determined by the density diffusivity κ rather than the viscosity: κ/λ ² acts as a vertical eddy viscosity, and the transport is λ ²/κ times the applied forcing. Critical values for the topographic heights in these regimes are identified.     

Field study on flow structures within aquatic vegetation under combined currents and small-scale waves

Field measurements were conducted to study the influence of aquatic vegetation on flow structures in floodplains under combined currents and wind-driven waves. Wave and turbulent velocities were decomposed from the time series of instantaneous velocity and analysed separately. In the present study, the wind waves were small, leading to the ratios of wave excursion (E_w) to stem spacing (S) for all cases tested here were less than 0.5. This caused the vertical distributions of time-averaged velocity (U_horiz) and turbulent kinetic energy (TKE) impacted by vegetation similar with the vegetated flow structures under pure current conditions. For emergent vegetation, U_horiz and TKE distributed uniformly through the entire water column or increased slightly from bed to water surface. Similar distributions were present in the lower part of submerged vegetation. In the upper part of submerged vegetation, U_horiz and TKE increased rapidly toward water surface and TKE reached its maximum near the top of vegetation. The measured wave orbital velocity (U_w) fitted linear wave theory well through the entire water depth for both the emergent and submerged cases, so that with small E_w/S the wave velocity was not attenuated within vegetation and U_w within canopy can be predicted by the linear wave theory under combined currents and waves. However, wind-driven waves made the turbulence generated near the top of canopy penetrate a deeper depth into vegetation than predictions under pure current conditions.     

Variations of vertical velocity in the deep oceans simulated by a 1/10° OGCM

This paper analyzes variations of vertical velocity w simulated by the 1/10° Ocean General Circulation Model for the Earth Simulator (OFES). Strong w-variability is found in the deep oceans. When w is WKBJ-normalized, the standard deviation averaged over the Southern Ocean increases with depth and is larger than 8 × 10^− 3 cm/s throughout the water column below 1,500 m. Evidences are presented that link this w-variability to internal waves generated by quasi-steady currents over topography. The aliasing errors in lag-3-day correlations suggest a bottom generation of near-inertial waves. A scale analysis indicates that vertically propagating waves that can be resolved by the OFES model are waves with frequencies of the order of inertial frequency and wavelengths comparable to the order of the grid size. The vertical energy flux associated with these waves is substantial. When integrated globally, the vertical energy flux is upward in the upper 4 km and reaches maximum values of about 0.8 TW at about 1 to 2 km depth. Thus, the w-variability in the 1/10° OFES integration points not only to a strong bottom generation of near-inertial internal waves in the deep Southern Ocean but also to the possibility that the power provided by internal waves generated by non-tidal currents over topography can be comparable to the power provided by internal waves generated by tidal flows over topography.     

Sensitivity analysis of a morphodynamic modelling system applied to a coastal lagoon inlet

This work investigates the recent morphological changes at the inlet of a complex coastal system (Ria de Aveiro lagoon, Portugal). This study was carried out using bathymetric data analysis and numerical simulations obtained with the 2DH morphodynamic modelling system MORSYS2D. The present simulations considered only tidal forcing, and a sensitivity analysis was performed by tuning the formula used to compute the sediment transports. A non-uniform sediment grain size distribution for the Ria de Aveiro inlet is considered in the numerical simulations, based on surveys performed in this area. The model results are analysed to assess if they resemble the observed trends of erosion and deposition, as calculated from bathymetric data. A quantitative analysis of the differences between the bathymetric changes obtained through surveys and the numerical results over a period of 3 years considering different sediment transport formulations shows that the formulations of Ackers and White (1973) and Engelund and Hansen (1967) are the ones that best describe the morphodynamic changes driven by tides in the Ria de Aveiro inlet.     

A discontinuous finite element baroclinic marine model on unstructured prismatic meshes

We describe the space discretization of a three-dimensional baroclinic finite element model, based upon a discontinuous Galerkin method, while the companion paper (Comblen et al. 2010a) describes the discretization in time. We solve the hydrostatic Boussinesq equations governing marine flows on a mesh made up of triangles extruded from the surface toward the seabed to obtain prismatic three-dimensional elements. Diffusion is implemented using the symmetric interior penalty method. The tracer equation is consistent with the continuity equation. A Lax–Friedrichs flux is used to take into account internal wave propagation. By way of illustration, a flow exhibiting internal waves in the lee of an isolated seamount on the sphere is simulated. This enables us to show the advantages of using an unstructured mesh, where the resolution is higher in areas where the flow varies rapidly in space, the mesh being coarser far from the region of interest. The solution exhibits the expected wave structure. Linear and quadratic shape functions are used, and the extension to higher-order discretization is straightforward.     

On the reynolds stress and lateral eddy viscosity due to transient oceanic Rossby waves

Summary From numerical solutions of a wind-driven homogeneous ocean model, anegative lateral eddy viscosity of the order 10⁴ cm² sec^–1 is inferred from the large-scale time-dependent currents in the interior of an enclosed shallow basin. The transient Rossby waves in this region produce a systematic convergence of eddy momentum at the latitude of the maximum average eastward current, and thus effect a transfer of zonal momentum from the large-scale eddies to the mean flow. In this sense they are analogous to the Rossby waves in the atmospheric general circulation, and it is speculated that such waves may help to maintain the mean zonal ocean currents. Although this negative viscosity induced by the large-scale transients is relatively small compared with the prescribed lateral viscosity of 10⁸ cm² sec^–1 and should be given a quite different physical interpretation, it is evidently an important viscous effect for the mean flow in the interior of the basin. The prescribed viscosity, on the other hand, is effective in controlling the model's simulated sub-grid scale dissipation, which occurs almost entirely in the nearby steady boundary currents.     

An automatic scheme for baseline correction of strong-motion records in coseismic deformation determination

Coseismic deformation can be determined from strong-motion records of large earthquakes. Iwan et al. (Bull Seismol Soc Am 75:1225–1246, 1985) showed that baseline corrections are often required to obtain reliable coseismic deformation because baseline offsets lead to unrealistic permanent displacements. Boore (Bull Seismol Soc Am 91:1199–1211, 2001) demonstrated that different choices of time points for baseline correction can yield realistically looking displacements, but with variable amplitudes. The baseline correction procedure of Wu and Wu (J Seismol 11:159–170, 2007) improved upon Iwan et al. (Bull Seismol Soc Am 75:1225–1246, 1985) and achieved stable results. However, their time points for baseline correction were chosen by a recursive process with an artificial criterion. In this study, we follow the procedure of Wu and Wu (J Seismol 11:159–170, 2007) but use the ratio of energy distribution in accelerograms as the criterion to determine the time points of baseline correction automatically, thus avoiding the manual choice of time points and speeding up the estimation of coseismic deformation. We use the 1999 Chi-Chi earthquake in central Taiwan and the 2003 Chengkung and 2006 Taitung earthquakes in eastern Taiwan to illustrate this new approach. Comparison between the results from this and previous studies shows that our new procedure is suitable for quick and reliable determination of coseismic deformation from strong-motion records.     

Can zonally symmetric inertial waves drive an oscillating zonal mean flow?

In the present paper zonal mean flow excitation by inertial waves is studied in analogy to mean flow excitation by gravity waves that plays an important role for the quasi-biennial oscillation in the equatorial atmosphere. In geophysical flows that are stratified and rotating, pure gravity and inertial waves correspond to the two limiting cases: gravity waves neglect rotation, inertial waves neglect stratification. The former are more relevant for fluids like the atmosphere, where stratification is dominant, the latter for the deep oceans or planet cores, where rotation dominates. In the present study a hierarchy of simple analytical and numerical models of zonally symmetric inertial wave-mean flow interactions is considered and the results are compared with data from a laboratory experiment. The main findings can be summarised as follows: (i) when the waves are decoupled from the mean flow they just drive a retrograde (eastward) zonal mean flow, independent of the sign of the meridional phase speed; (ii) when coupling is present and the zonal mean flow is assumed to be steady, the waves can drive vertically alternating jets, but still, in contrast to the gravity wave case, the structure is independent of the sign of the meridional phase speed; (iii) when coupling is present and time-dependent zonal mean flows are considered the waves can drive vertically and temporarily oscillating mean flows. The comparison with laboratory data from a rotating annulus experiment shows a qualitative agreement. It appears that the experiment captures the basic elements of the inertial wave mean flow coupling. The results might be relevant to understand how the Equatorial Deep Jets can be maintained against dissipation, a process currently discussed controversially.     

Selection and ranking of ground motion models for seismic hazard analysis in the Pyrenees

The issue addressed in this paper is the objective selection of appropriate ground motion models for seismic hazard assessment in the Pyrenees. The method of Scherbaum et al. (2004a) is applied in order to rank eight published ground motion models relevant to intraplate or to low deformation rate contexts. This method is based on a transparent and data-driven process which quantifies the model fit and also measures how well the underlying model assumptions are met. The method is applied to 15 accelerometric records obtained in the Pyrenees for events of local magnitude between 4.8 and 5.1, corresponding to moment magnitudes ranging from 3.7 to 3.9. Only stations at rock sites are considered. A total of 720 spectral amplitudes are used to rank the selected ground motion models. Some control parameters of these models, such as magnitude and distance definitions, may vary from one model to the other. It is thus important to correct the selected models for their difference with respect to the magnitude and distance definitions used for the Pyrenean data. Our analysis shows that, with these corrections, some of the ground motion models successfully fit the data. These are the Lussou et al. (2001) and the Berge-Thierry et al. (2003) models. According to the selected ground motion models, a possible scenario of a magnitude 6 event is proposed; it predicts response spectra accelerations of 0.08–0.1 g at 1 Hz at a hypocentral distance of 10 km.     

Moment inequality and complete convergence of moving average processes under asymptotically linear negative quadrant dependence assumptions

Let {Y, Y _i , −∞ < i < ∞} be a doubly infinite sequence of identically distributed and asymptotically linear negative quadrant dependence random variables, {a _i , −∞ < i < ∞} an absolutely summable sequence of real numbers. We are inspired by Wang et al. (Econometric Theory 18:119–139, 2002) and Salvadori (Stoch Environ Res Risk Assess 17:116–140, 2003). And Salvadori (Stoch Environ Res Risk Assess 17:116–140, 2003) have obtained Linear combinations of order statistics to estimate the quantiles of generalized pareto and extreme values distributions. In this paper, we prove the complete convergence of under some suitable conditions. The results obtained improve and generalize the results of Li et al. (1992) and Zhang (1996). The results obtained extend those for negative associated sequences and ρ^*-mixing sequences. CIC Number O211, AMS (2000) Subject Classification 60F15, 60G50 Research supported by National Natural Science Foundation of China     

Harmonization check of M w within the central, northern, and northwestern European earthquake catalogue (CENEC)

Large data sets covering large areas and time spans and composed of many different independent sources raise the question of the obtained degree of harmonization. The present study is an analysis of the harmonization with respect to the moment magnitude M _w within the earthquake catalogue for central, northern, and northwestern Europe (CENEC). The CENEC earthquake catalogue (Grünthal et al., J Seismol, 2009) contains parameters for over 8,000 events in the time period 1000–2004 with magnitude M _w ≥ 3.5. Only about 2% of the data used for CENEC have original M _w magnitudes derived directly from digital data. Some of the local catalogues and data files providing data give M _w, but calculated by the respective agency from other magnitude measures or intensity. About 60% of the local data give strength measures other than M _w, and these have to be transformed by us using available formulae or new regressions based on original M _w data. Although all events are thus unified to M _w magnitude, inhomogeneity in the M _w obtained from over 40 local catalogues and data files and 50 special studies is inevitable. Two different approaches have been followed to investigate the compatibility of the different M _w sets throughout CENEC. The first harmonization check is performed using M _w from moment tensor solutions from SMTS and Pondrelli et al. (Phys Earth Planet Inter 130:71–101, 2002; Phys Earth Planet Inter 164:90–112, 2007). The method to derive the SMTS is described, e.g., by Braunmiller et al. (Tectonophysics 356:5–22, 2002) and Bernardi et al. (Geophys J Int 157:703–716, 2004), and the data are available in greater extent since 1997. One check is made against the M _w given in national catalogues and another against the M _w derived by applying different empirical relations developed for CENEC. The second harmonization check concerns the vast majority of data in CENEC related to earthquakes prior to 1997 or where no moment tensor based M _w exists. In this case, an empirical relation for the M _w dependence on epicentral intensity (I ₀) and focal depth (h) was derived for 41 master events, i.e., earthquakes, located all over central Europe, with high-quality data. To include also the data lacking h, the corresponding depth-independent relation for these 41 events was also derived. These equations are compared with the different sets of data from which CENEC has been composed, and the goodness of fit is demonstrated for each set. The vast majority of the events are very well or reasonably consistent with the respective relation so that the data can be said to be harmonized with respect to M _w, but there are exceptions, which are discussed in detail.     

Summer surface layer thermal response to surface gravity waves in the Yellow Sea

The Princeton Ocean Model (POM) with generalized coordinate system (POMgcs) is used to study the summer surface-layer thermal response to surface gravity waves in the Yellow Sea (YS). The parameterization schemes of wave breaking developed by Mellor and Blumberg (J Phys Oceanogr 34:693–698, 2004) and Kantha and Clayson (Ocean Model 6:101–124, 2004), respectively, and Stokes production developed by Kantha and Clayson (Ocean Model 6:101–124, 2004) are both included in the Mellor–Yamada turbulence closure model Mellor and Yamada (Rev Geophys 20:851–875, 1982) of POMgcs. Numerical results show that surface gravity waves impact the depth of surface mixed layer of temperature in the YS in summer. The surface mixed layer in the YS cannot be reproduced well and has a visible difference from the observation if the parameterization schemes are not included. A diagnostic analysis of turbulent kinetic energy suggests that both Stokes production and wave breaking play key roles in enhancing the turbulent mixing near the sea surface in the YS. Stokes production seems to have a greater impact throughout the upper mixed layer in the YS in summer than that of wave breaking. In addition, a diagnostic analysis of the momentum balance shows that Coriolis–Stokes forcing has a significant effect on the momentum budget in the upper layer in the YS, and surface gravity waves are able to reduce the velocity of mean flow near the surface and make the mean flow near the surface more homogeneous vertically in the YS.     

Measurement of long waves at the harbor of Marina di Carrara,Italy

This paper presents the analysis of measurements of long waves at the harbor of Marina di Carrara, Italy. Nine pressure gauges (eight in the inner harbor and one outside the harbor) were used to record continuously the water surface elevation, for about 3 years; incoming short waves were measured using a directional buoy at about 500 m offshore the breakwaters. The analysis is carried out splitting the energy of the long waves into two bands, i.e., very long waves (VLW; f < 0.003 Hz) and long waves (LW; 0.003 < f < 0.030 Hz); it is found that LW energy is strongly correlated with the energy of the incoming short waves while no correlation (or very poor) is found for VLW. During energetic sea states, surface elevation spectra of VLW and LW, normalized against the wave energy of each frequency band, appear to be autosimilar; this suggests that the spectra of the incoming long waves have a mostly constant shape. Simple numerical computations of harbor resonance, carried out using a mild-slope equation linearized model, indicate that the amplification of the height of the incoming LW generated by short waves (calculated using empirical formulae from previous researches by Melito et al. 2007) inside the harbor can be reasonably predicted using such simple approach, assuming that the LW spectra are rectangular. The shape of the amplification diagram calculated by the numerical model shows reasonable (but not perfect) agreement with the shape of the normalized measured spectra inside the harbor.     

New light on a dark subject: comment

In their article, “New light on a dark subject: On the use of fluorescence data to deduce redox states of natural organic matter,” Macalady and Walton-Day (2009) subjected natural organic matter (NOM) samples to oxidation, reduction, and photochemical transformation. Fluorescence spectra were obtained on samples, which were diluted “to bring maximum uvvisible absorbance values below 1.0.” The spectra were fit to the Cory and McKnight (2005) parallel factor analysis (PARAFAC) model, and consistent variation in the redox state of quinone-like moieties was not detected. Based on these results they concluded that fitting fluorescence spectra to the Cory and McKnight (2005) PARAFAC model to obtain information about the redox state of quinone-like moieties in NOM is problematic. Recognizing that collection and correction of fluorescence spectra requires consideration of many factors, we investigated the potential for inner-filter effects to obscure the ability of fluorescence spectroscopy to quantify the redox state of quinone-like moieties. We collected fluorescence spectra on Pony Lake and Suwannee River fulvic acid standards that were diluted to cover a range of absorbance wavelengths, and fit these spectra to the Cory and McKnight (2005) PARAFAC model. Our results suggest that, in order for the commonly used inner-filter correction to effectively remove inner-filter effects, samples should be diluted such that the absorbance at 254 nm is less than 0.3 prior to the collection of fluorescence spectra. This finding indicates that inner-filter effects may have obscured changes in the redox signature of fluorescence spectra of the highly absorbing samples studied by Macalady and Walton-Day (2009).     

Approximate recovery of coseismic deformation from Taiwan strong-motion records

Since 1990, digital strong-motion accelerographs and global positioning system (GPS) instruments have been widely deployed in the Taiwan region (Shin et al. 2003; Yu et al. 2001). The 1999 Chi-Chi, Mw 7.6 earthquake and the 2003 Chengkung, Mw 6.8 earthquake were well recorded by both digital accelerographs and GPS instruments. These data offer a good opportunity to determine coseismic displacements from strong-motion records and to compare the results with those derived from GPS measurements. As noted by Boore (2001), a double integration of the acceleration data often leads to unreasonable results, and baseline corrections are therefore required in most cases before the integration. Based on the works of Iwan et al. (1985) and Boore (2001), we developed an improved method for baseline correction and validated it using an extensive set of data from shake-table tests of a known “step” displacement on 249 accelerographs. Our baseline correction method recovered about 97% of the actual displacement from the shake-table data. We then applied this baseline correction method to compute coseismic displacements from the strong-motion data of the Chi-Chi and Chengkung earthquakes. Our results agree favorably with the coseismic displacements determined by the GPS measurements at nearby sites. The ratio of seismic to geodetic displacement varies from 0.78 to 1.41, with an average of about 1.05.     

Observations of estuarine circulation and solitary internal waves in a highly energetic tidal channel

Despite vigorous tidal and wind mixing, observations in an estuarine tidal inlet in the Wadden Sea show that during part of the tidal cycle, vertical stratification and internal waves may still develop. Acoustic Doppler current profiler (ADCP) and conductivity, temperature, depth observations, collected over the past 6 years at 13 h anchor stations (ASs), reveal that these occur especially during slack tide, when there is little wind and large freshwater discharge from nearby Lake IJssel. Measurements with a moored ADCP show that in the same tidal phase, strong cross-channel circulation develops, which may suddenly reverse circulation sense due to passing density fronts. In the vertically stratified phase that follows after the front passage, propagating mode-one solitary internal waves are observed. These are resonantly generated during decelerating tidal ebb currents when the (shear) flow passes a transcritical regime (Froude number equal to 1). A combination of photographs (including one from the International Space Station), bathymetric data, and ASs data leads to the discovery of yet another source of internal waves in this area, produced during slackening tide by propagating lee waves that develop over a deep trench. We suggest that both the cross-channel circulation as well as the (solitary) internal waves may locally be of importance for the (re)distribution and transport of sediments and nutrients and may influence tidally averaged transports.