Contribution Of Waves Research Articles

The Contribution of High-Frequency Wind-Generated Surface Waves to the Stokes Drift

AbstractThe effects of nonbreaking surface waves on upper-ocean dynamics enter the wave-averaged primitive equations through the Stokes drift. Through the resulting upper-ocean dynamics, Stokes drift is a catalyst for the fluxes of heat and trace gases between the atmosphere and ocean. However, estimates of the Stokes drift rely crucially on properly resolving the wave spectrum. In this paper, using state-of-the-art spatial measurements (in situ and airborne remote sensing) from a number of different field campaigns, with environmental conditions ranging from 2 to 13 m s−1wind speed and significant wave height of up to 4 m, we characterize the properties of the surface wave field across the equilibrium and saturation ranges and provide a simple parameterization of the transition between the two regimes that can easily be implemented in numerical wave models. We quantify the error associated with instrument measurement limitations, or incomplete numerical parameterizations, and propose forms for the continuation of these spectra to properly estimate the Stokes drift. Depending on the instrument and the sea state, predictions of surface Stokes drift may be underestimated by more than 50%.

Three-body structure of B19 : Finite-range effects in two-neutron halo nuclei

The structure and $B(E1)$ transition strength of $^{19}$B are investigated in a $^{17}\text{B}+n+n$ model, triggered by a recent experiment showing that $^{19}$B exhibits a well pronounced two-neutron halo structure. Preliminary analysis of the experimental data was performed by employing contact $n$-$n$ interactions, which are known to underestimate the $s$-wave content in other halo nuclei such as $^{11}$Li. In the present work, the three-body hyperspherical formalism with finite-range two-body interactions is used to describe $^{19}$B. In particular, two different finite-range $n$-$n$ interactions will be used, as well as a simple central Gaussian potential whose range is progressively reduced. The purpose is to determine the main properties of the nucleus and investigate how they change when using contact-like $n$-$n$ potentials. Special attention is also paid to the dependence on the prescription used to account for three-body effects, i.e., a three-body force or a density-dependent $n$-$n$ potential. We have found that the three-body model plus finite-range potentials provide a description of $^{19}$B consistent with the experimental data. The results are essentially independent of the short-distance details of the two-body potentials, giving rise to an $(s_{1/2})^2$ content of about 55%, clearly larger than the initial estimates. Very little dependence has been found as well on the prescription used for the three-body effects. The total computed $B(E1)$ strength is compatible with the experimental result, although we slightly overestimate the data around the low-energy peak of the $dB(E1)/d\varepsilon$ distribution. Finally, we show that a reduction of the $n$-$n$ interaction range produces a significant reduction of the $s$-wave contribution, which then should be expected in calculations using contact interactions.

Measuring ocean total surface current velocity with the KuROS and KaRADOC airborne near-nadir Doppler radars: a multi-scale analysis in preparation for the SKIM mission

Abstract. Surface currents are poorly known over most of the world's oceans. Satellite-borne Doppler wave and current scatterometers (DWaCSs) are among the proposed techniques to fill this observation gap. The Sea surface KInematics Multiscale (SKIM) proposal is the first satellite concept built on a DWaCS design at near-nadir angles and was demonstrated to be technically feasible as part of the European Space Agency Earth Explorer program. This article describes preliminary results from a field experiment performed in November 2018 off the French Atlantic coast, with sea states representative of the open ocean and a well-known tide-dominated current regime, as part of the detailed design and feasibility studies for SKIM. This experiment comprised airborne measurements performed using Ku-band and Ka-band Doppler radars looking at the sea surface at near-nadir incidence in a real-aperture mode, i.e., in a geometry and mode similar to that of SKIM, as well as an extensive set of in situ instruments. The Ku-band Radar for Observation of Surfaces (KuROS) airborne radar provided simultaneous measurements of the radar backscatter and Doppler velocity in a side-looking configuration, with a horizontal resolution of about 5 to 10 m along the line of sight and integrated in the perpendicular direction over the real-aperture 3 dB footprint diameter (about 580 m). The Ka-band RADar for Ocean Current (KaRADOC) system, also operating in the side-looking configuration, had a much narrower beam, with a circular footprint only 45 m in diameter. Results are reported for two days with contrasting conditions, a strong breeze on 22 November 2018 (wind speed 11.5 m s−1, Hs 2.6 m) and gentle breeze on 24 November 2018 (wind speed 5.5 m s−1, Hs 1.7 m). The measured line-of-sight velocity signal is analyzed to separate a non-geophysical contribution linked to the aircraft velocity, a geophysical contribution due to the intrinsic motion of surface waves and the desired surface current contribution. The surface wave contribution is found to be well predicted by Kirchhoff scattering theory using as input parameters in situ measurements of the directional spectrum of long waves, complemented by the short wave spectrum of Elfouhaily et al. (1997). It is found to be closely aligned with the wind direction, with small corrections due to the presence of swell. Its norm is found to be weakly variable with wind speed and sea state, quite stable and close to C0=2.0ms-1 at the Ka band, and more variable and close to C0=2.4ms-1 at the Ku band. These values are 10 %–20 % smaller than previous theoretical estimates. The directional spread of the short gravity waves is found to have a marked influence on this surface wave contribution. Overall, the results of this study support the feasibility of near-nadir radar Doppler remote sensing of the ocean total surface current velocity (TSCV).

Diffracted wavefield decomposition and multidimensional site effects in the Argostoli valley, Greece

SUMMARYEffects of seismic ground motion induced by surface geology and geometry are known to be associated with the generation of a substantial proportion of surface waves. As a consequence, surface waves significantly contribute to ground-motion variability and site amplification. There is a growing body of literature recognizing that an understanding of physical patterns of the wavefield crossing a site is the key aspect to characterize and quantify them. However, this task remains technically challenging due to the complexity of such effects as well as the limitations of geophysical investigations, especially in case of small sedimentary valleys. The present study attempts to investigate the waves propagating across two 2-D dense seismic arrays from a number of earthquakes and explore the extent to which they are contributing to the multidimensional site effects. The arrays were deployed in the small-size, shallow alluvial valley of Koutavos-Argostoli, located in Cephalonia Island, Greece, and consisted of three-component velocimeters with interstation distances ranging from 5 to 160 m. A set of 46 earthquakes, with magnitudes between 2 and 5 and epicentral distances up to 200 km, was analysed by using an advanced seismic array processing technique, MUSIQUE. The phase velocity, backazimuth and energy of the dominant waves crossing the array were extracted, and their identification as Love or prograde/retrograde Rayleigh waves was obtained. The results clearly indicate a predominance of scattered surface waves (up to 60 per cent of total energy), mainly from the closest valley edges, above the fundamental frequency (∼1.5 Hz) of the valley. Love waves dominate the low-frequency wavefield (<3 Hz) while Rayleigh waves dominate some high-frequency bands. An excellent consistency is observed, in a given frequency range, among the dominance of the type of diffracted surface waves, group velocities estimated from the ground velocity structure and site amplification. The outcomes of this research provide a better understanding of the contribution of edge-diffracted surface waves and the 2-D/3-D site amplification at small and shallow alluvium valleys like Argostoli. The method applied here can be used to calibrate and validate 3-D models for simulating seismic ground motion.

Temporal Modulation of the Local Microseism in the North Sea

AbstractPrimary and secondary microseism originating in the world oceans and peaking at around 14 and 7 s, respectively, characterize the Earth's background noise in that frequency range. Microseism generated in marginal seas with partly shorter periods and higher spatial and temporal variability is less studied and requires stations in immediate proximity to the source to be observed. Such studies can help to elucidate the exact microseism generation areas and mechanisms in a constrained area. We analyze 15 years of broadband data recorded at the seismic station on Helgoland island in the marginal North Sea. In addition to remote primary (RPM) and secondary microseism (RSM) originating in the North Atlantic, we observe strong and dominant local secondary microseism (LSM) with on average higher frequencies above 0.2 Hz, in accordance with shorter wave periods of about 4–8 s in the shallow North Sea. During times with low RSM activity we observe local primary microseism (LPM) at frequencies in agreement with local ocean wave periods. The higher horizontal to vertical (H/V) ratio of LPM with respect to LSM indicates a major non‐Rayleigh wave contribution. LSM and LPM show a strong modulation with local semidiurnal ocean tides and microseism energy maxima preceding the water level maximum by 2.5 and 1.5 hr, respectively. This time shift might be influenced by stronger currents during rising than falling tides. Active sources of tide‐modulated microseism migrate along the North Sea coast in sync with the ocean tidal signal as evidenced by comparison of LSM maxima at stations distributed along the coast.

Thermocline Warming Induced Extreme Indian Ocean Dipole in 2019

AbstractThe 2019 positive Indian Ocean Dipole (IOD) was the strongest event since the 1960s which developed independently without coinciding El Niño. The dynamics is not fully understood. Here we show that in March–May, westward propagating oceanic Rossby waves, a remnant consequence of the weak 2018 Pacific warm condition, led to anomalous sea surface temperature warming in the southwest tropical Indian Ocean (TIO), inducing deep convection and anomalous easterly winds along the equator, which triggered the initial cooling in the east. In June–August, the easterly wind anomalies continued to evolve through ocean‐atmosphere coupling involving Bjerknes feedback and equatorial nonlinear ocean advection, until its maturity in September–November. This study clarifies the contribution of oceanic Rossby waves in the south TIO in different dynamic settings and reveals a new triggering mechanism for extreme IOD events that will help to understand IOD diversity.

Multitemporal analysis of ambient noise polarization to characterize site response in the town of Amatrice, shattered by the 2016 central Italy earthquake

SUMMARYThe Mw 6.0 earthquake that hit central Italy on 2016 August 24 caused an abnormally high level of destruction in the town of Amatrice. In order to clarify the role of site response in causing such a disaster, a series of ambient noise recordings acquired in the aftermath of the event are analysed here to identify site resonance properties from the ratios H/V between horizontal and vertical amplitudes of ground motion. Although the noise data acquisition was limited by the emergency management activities, the use of a new analysis technique, which isolates the contribution of Rayleigh waves to the noise wavefield and averages instantaneous estimates of H/V ratios, provided more stable results compared to the standard Nakamura's technique based on mean spectral ratios. The results demonstrated the occurrence of significant resonance phenomena, but without an obvious correlation with the spatial distribution of damage severity. It is apparent that the damage severity was also influenced by some additional local factors related to building vulnerability. Moreover, the time-series analysis revealed seasonal variations in the Rayleigh wave ellipticity curves likely related to the water content changes in the surface deposits and their influence on the Poisson coefficient. Finally, the new method proved capable of recognizing time-varying directions of Rayleigh wave propagation. This capability could be exploited to support other passive seismic methods (e.g. ReMi), whose results’ reliability is limited by the lack of control on wave direction origin.

Properties of nuclear matter in relativistic Brueckner–Hartree–Fock model with high-precision charge-dependent potentials

Properties of nuclear matter are investigated in the framework of relativistic Brueckner–Hartree–Fock model with the latest high-precision charge-dependent Bonn (pvCD-Bonn) potentials, where the coupling between pion and nucleon is adopted as a pseudovector form. These realistic pvCD-Bonn potentials are renormalized to effective nucleon–nucleon (NN) interactions, G matrices. They are obtained by solving the Blankenbecler-Sugar (BbS) equation in nuclear medium. Then, the saturation properties of symmetric nuclear matter are calculated with pvCD-Bonn A, B, C potentials. The energies per nucleon are around −10.72 MeV to −16.83 MeV at saturation densities, 0.139 fm−3 to 0.192 fm−3 with these three potentials, respectively. It clearly demonstrates that the pseudovector coupling between pion and nucleon can generate reasonable saturation properties comparing with pseudoscalar coupling. Furthermore, these saturation properties have strong correlations with the tensor components of NN potentials, i.e., the D-state probabilities of deuteron, PD to form a relativistic Coester band. The equations of state of pure neutron matter from pvCD-Bonn potentials are almost identical, since the prominent difference of pvCD Bonn potentials are the components of tensor force, which provides very weak contributions in the case of total isospin T = 1. In addition, the charge symmetry breaking (CSB) and charge independence breaking (CIB) effects are also discussed in nuclear matter from the partial wave contributions with these high-precision charge-dependent potentials. In general, the magnitudes of CSB from the differences between nn and pp potentials are about 0.05 MeV, while those of CIB are around 0.35 MeV from the differences between np and pp potentials.

Mid-IR photothermal measurement of substantial heat transport by surface waves of polar amorphous films supported on silicon

We present measurements of significant thermal diffusivity by surface electromagnetic waves of an ultra-thin polar and amorphous dielectric film deposited on silicon (Si). We used a photothermal-beam-deflection technique with a modulated mid-infrared heating source to excite and launch surface electromagnetic waves onto the surface of an amorphous silicon carbide (a-SiC) film deposited on Si and generate periodic temperature and refractive index gradients above the sample surface. These gradients are capable of periodically deflecting a probe beam, passing very close to the surface, at the modulation frequency of the heating beam. We have fitted the measured probe beam deflection to an analytical model for the mirage effect that takes into account the thermal anisotropy of the measured sample to infer the contribution of the surface electromagnetic waves of the a-SiC film to thermal diffusivity in the plane of the sample under study. We found that reducing the thickness of the a-SiC film promotes the interaction between the surface electromagnetic waves propagating on either side of the a-SiC film, which significantly enhances thermal diffusivity in the plane of the measured sample. We also found that in-plane thermal diffusivity by surface electromagnetic waves on an amorphous silicon carbide film a few nanometers thick is several orders of magnitude greater than thermal diffusivity by phonons in silicon. We believe that the results obtained provide a better understanding of the physics of electromagnetic waves confined to solid surfaces.

Dynamics of Energetic Electrons in the Slot Region During Geomagnetically Quiet Times: Losses Due to Wave‐Particle Interactions Versus a Source From Cosmic Ray Albedo Neutron Decay (CRAND)

AbstractEarth's slot region, lying between the outer and inner radiation belts, has been identified as due to a balance between inward radial diffusion and pitch angle (PA) scattering induced by waves. However, recent satellite observations and modeling studies indicate that cosmic ray albedo neutron decay (CRAND) may also play a significant role in energetic electron dynamics in the slot region. In this study, using a drift‐diffusion‐source model, we investigate the relative contribution of all significant waves and CRAND to the dynamics of energetic electrons in the slot region during July 2014, an extended period of quiet geomagnetic activity. The bounce‐averaged PA diffusion coefficients from three types of waves (hiss, lightning‐generated whistlers [LGW], and very low frequency [VLF] transmitters) are calculated based on quasi‐linear theory, while the CRAND source follows the results in Xiang et al. (2019, https://doi.org/10.1029/2018GL081730). The simulation results indicate that both LGW and VLF transmitter waves can enhance loss and weaken the top hat PA distribution induced by hiss waves. For 470 keV electrons at L = 2.5, simulation results without CRAND show a much quicker decrease than observations from the Van Allen Probes. After including CRAND, simulated electron flux variations reproduce satellite observations, suggesting that CRAND is an important source for hundreds of keV electrons in the slot region during quiet times. The balance between the CRAND source and loss due to wave‐particle interactions provides a lower limit to relativistic electron fluxes in the slot region, which can act as an important reference point for instrument calibration when a true background level is warranted.

Why Do Benguela Niños Lead Atlantic Niños?

AbstractWe investigate the lag between warm interannual Sea Surface Temperature (SST) events in the eastern‐equatorial Atlantic, the Atlantic Niños, and the occurrence of Benguela Niños along the southwestern Angolan coast. While it is commonly agreed that both events are associated with equatorial and subsequent coastal‐trapped wave propagations driven remotely by a relaxation of the trade‐winds, it is surprising that SST anomalies off Angola tend to precede the ones in the eastern‐equatorial sector by ~1 month. To explain this counterintuitive behavior, our methodology is based on the experimentation with a Tropical Atlantic Ocean model. Using idealized wind‐stress perturbations from a composite analysis, we trigger warm equatorial and coastal events over a stationary and then, seasonally varying ocean mean‐state. In agreement with the linear dynamics, our results show that when the interannual wind‐stress forcing is restricted to the western‐central equatorial Atlantic, the model yields equatorial events leading the coastal ones. This implies that neither the differences in the ocean stratification between the two regions (thermocline depths or modal wave contributions) nor the seasonal phasing of the events explains the observed temporal sequence. Only if wind‐stress anomalies are also prescribed in the coastal fringe, the coastal warming precedes the eastern‐equatorial SST anomaly peak, emphasizing the role of the local forcing in the phenology of Benguela Niños. A weaker South‐Atlantic Anticyclone initiates the coastal warming before the development of eastern‐equatorial SST anomalies. Then, equatorward coastal wind anomalies, driven by a convergent anomalous circulation located on the warm Atlantic Niño, stop the remotely forced coastal warming prematurely.

Airborne measurements and large-eddy simulations of small-scale gravity waves at the tropopause inversion layer over Scandinavia

Abstract. Coordinated airborne measurements were performed by two research aircraft – Deutsches Zentrum für Luft- und Raumfahrt (DLR) Falcon and High Altitude and Long Range Aircraft (HALO) – in Scandinavia during the GW-LCYCLE II (Investigation of the life cycle of gravity waves) campaign in 2016 to investigate gravity wave processes in the upper troposphere and lower stratosphere (UTLS) region. A mountain wave event was probed over southern Scandinavia on 28 January 2016. The collected dataset constitutes a valuable combination of in situ measurements and horizontal- and altitude-resolved Doppler wind lidar and water vapour measurements with the differential absorption lidar (DIAL). In situ data at different flight altitudes and downward-pointing wind lidar measurements show pronounced changes of the horizontal scales in the vertical velocity field and of the leg-averaged momentum fluxes (MFs) in the UTLS region. The vertical velocity field was dominated by small horizontal scales with a decrease from around 20 to < 10 km in the vicinity of the tropopause inversion layer (TIL). These small scales were also found in the water vapour data and backscatter data of the DIAL. The leg-averaged MF profile determined from the wind lidar data is characterized by a pronounced kink of positive fluxes in the TIL and negative fluxes below. The largest contributions to the MF are from waves with scales > 30 km. The combination of the observations and idealized large-eddy simulations revealed the occurrence of interfacial waves having scales < 10 km on the tropopause inversion during the mountain wave event. The contribution of the interfacial waves to the leg-averaged MF is basically zero due to the phase relationship of their horizontal and vertical velocity perturbations. Interfacial waves have already been observed on boundary-layer inversions but their concept has not been applied to tropopause inversions so far. Our idealized simulations reveal that the TIL affects the vertical trend of leg-averaged MF of mountain waves and that interfacial waves can occur also on tropopause inversions. Our analyses of the horizontal- and altitude-resolved airborne observations confirm that interfacial waves actually do occur in the TIL. As predicted by linear theory, the horizontal scale of those waves is determined by the wind and stability conditions above the inversion. They are found downstream of the main mountain peaks and their MF profile varies around zero and can clearly be distinguished from the MF profile of Kelvin–Helmholtz instability. Further, the idealized large-eddy simulations reveal that the presence of the TIL is crucial in producing this kind of trapped wave at tropopause altitude.

Spatial evolution of young wind waves: numerical modelling verified by experiments

Abstract

Impacts of surface gravity waves on a tidal front: A coupled model perspective

A set of realistic coastal coupled ocean-wave numerical simulations is used to study the impact of surface gravity waves on a tidal temperature front and surface currents. The processes at play are elucidated through analyses of the budgets of the horizontal momentum, the temperature, and the turbulence closure equations. The numerical system consists of a 3D coastal hydrodynamic circulation model (Model for Applications at Regional Scale, MARS3D) and the third generation wave model WAVEWATCH III (WW3) coupled with OASIS-MCT at horizontal resolutions of 500 and 1500 m, respectively. The models were run for a period of low to moderate southwesterly winds as observed during the Front de Marée Variable (FroMVar) field campaign in the Iroise Sea where a seasonal small-scale tidal sea surface temperature front is present. Over the 2 day period considered, long fetch waves grow gradually propagating north east and east.Contrasting a stand-alone ocean run with a coupled ocean-wave run shows that waves move the Ushant front offshore by up to 4 kilometres and cool the offshore stratified side of the front by up to 1.5°C. The analysis of the temperature budget shows that the change in advection is the dominant factor contributing to the frontal shift while the contribution of wave enhanced vertical temperature diffusion is secondary. Temperature, considered to be a tracer, is advected in the coupled run by the Lagrangian current resulting from the quasi-Eulerian and Stokes drift. Although the Stokes drift is directed shorewards, changes in the quasi-Eulerian current lead to a more offshore advection in the coupled than the stand-alone run. The quasi-Eulerian current is reduced (enhanced) during the ebb (flood) flow which correspond to periods of wave-following (-opposing) currents. This is due to wave breaking enhanced vertical mixing acting on the positive vertical gradient present in the quasi-Eulerian current during both ebb and flood tides. Partially coupled runs reveal that it is the surface flux of TKE associated to wave breaking that is key rather than the changes in the surface stress. They further elucidate the role of other modelled wave related processes. Although the contribution of the Stokes–Coriolis force and the wave breaking induced enhancement in vertical mixing to the quasi-Eulerian current are of similar magnitude and sign, it does not contribute significantly to the frontal shift. This is because it partially counters the Stokes drift advection which pushes the front shorewards. All Stokes drift related processes combined thus only lead to a very slight displacement of the front.

Contribution of Wave Setup to Projected Coastal Sea Level Changes

AbstractAlong open coasts, wind waves are a key driver of coastal changes and can be major contributors to coastal hazards. Wind wave characteristics are projected to change in response to climate change, notably due to changes in atmospheric circulation patterns and the associated surface winds. Here, a first‐order estimate of projected 20‐yr mean wave setup changes, excluding extreme events and subannual variability, is provided for sandy beaches along most of the world's coastline over the middle and end of the 21st century. Calculations are based on an ensemble of wave model projections under the representative concentration pathways (RCP) 8.5 and on empirical formulations for wave setup. Projected wave setup changes are compared to other contributors currently accounted for in regional sea level projections to extend existing projections of 21st century coastal sea level changes. Projected wave setup changes exhibit a clear spatial heterogeneity and mostly average out at global scale. However, at regional or local scale, wave setup changes are a small yet nonnegligible contributor to total coastal sea level 20‐yr mean changes (which include global mean sea level rise, GMSLR) over the middle and end of the 21st century. Wave setup can be a substantial contributor to local departures of coastal sea level changes from GMSLR. Wave setup changes should therefore be included in projections of regional patterns of coastal sea level changes. The reported long‐term changes in wave setup also advocate for the inclusion of nonstationary wave contributions to projected regional patterns of coastal sea level changes, including for studies on extreme events.

Variations in storm-induced bed level dynamics across intertidal flats

Hydrodynamic forces on intertidal flats vary over a range of temporal and spatial scales. These spatiotemporal inhomogeneities have implications for intertidal flat morphodynamics and ecology. We determine whether storm events are capable of altering the long-term morphological evolution of intertidal flats, and unravel the contributions of tidal flow, wind-driven flow, waves, and water depth on inhomogeneities in bed level dynamics (bed level changes over ~days) across these areas. We complement decades of bed level measurements on eight intertidal flats in two estuaries in the Netherlands with an extensive 1-month field campaign on one of those flats. Across this intertidal flat, the hydrodynamics and morphodynamics of a storm event were captured, including the post-storm recovery. We show that individual events can persistently alter the morphological evolution of intertidal flats; magnitudes of some bed level changes are even comparable to years of continuous evolution. The morphological impacts of events are largely controlled by the relative timing of the forcing processes, and not solely by their magnitudes. Spatiotemporal variations in bed level dynamics of intertidal flats are driven by a combination of: (1) the inhomogeneous distributions of the hydrodynamic forcing processes (including the under-explored role of the wind); and (2) the linear proportionality between bed level dynamics and the local bed slope.

Numerical determination of the cutoff frequency in solar models

Context. In stratified atmospheres, acoustic waves can only propagate if their frequency is higher than the cutoff value. The determination of the cutoff frequency is fundamental for several topics in solar physics, such as evaluating the contribution of the acoustic waves to the chromospheric heating or the application of seismic techniques. However, different theories provide different cutoff values. Aims. We developed an alternative method to derive the cutoff frequency in several standard solar models, including various quiet-Sun and umbral atmospheres. The effects of magnetic field and radiative losses on the cutoff are examined. Methods. We performed numerical simulations of wave propagation in the solar atmosphere using the code MANCHA. The cutoff frequency is determined from the inspection of phase-difference spectra computed between the velocity signal at two atmospheric heights. The process is performed by choosing pairs of heights across all the layers between the photosphere and the chromosphere to derive the vertical stratification of the cutoff in the solar models. Result. The cutoff frequency predicted by the theoretical calculations departs significantly from the measurements obtained from the numerical simulations. In quiet-Sun atmospheres, the cutoff shows a strong dependence on the magnetic field for adiabatic wave propagation. When radiative losses are taken into account, the cutoff frequency is greatly reduced and the variation of the cutoff with the strength of the magnetic field is lower. The effect of the radiative losses in the cutoff is necessary to understand recent quiet-Sun and sunspot observations. In the presence of inclined magnetic fields, our numerical calculations confirm that the cutoff frequency is reduced as a result of the reduced gravity experienced by waves that propagate along field lines. An additional reduction is also found in regions with significant changes in the temperature, which is due to the lower temperature gradient along the path of field-guided waves. Conclusions. Our results show solid evidence that the cutoff frequency in the solar atmosphere is stratified. The cutoff values are not correctly captured by theoretical estimates. In addition, most of the widely used analytical cutoff formulae neglect the effect of magnetic fields and radiative losses, whose role is critical for determining the evanescent or propagating nature of the waves.

Adaptation time to magnified flood hazards underestimated when derived from tide gauge records

Sea-level rise magnifies flood hazards, raising the question when adaptation measures need to be taken. Here, we quantify when the recurrence of extreme water level events will double due to projected sea-level rise. Reproducing the most common method based on extreme water levels observed with tide gauges, at least one third of the coastal locations are to expect a doubling of extremes within a decade. However, tide gauges are commonly placed in wave-sheltered harbours where the contribution of waves to water levels is much smaller than at nearby wave-exposed coastlines such as beaches and dikes. In this study, we quantify doubling times at a variety of idealised shorelines based on modelled tides, storm surges and waves. We apply an extreme value analysis that accounts for the joint probability of extreme storm surges and extreme waves. Our results indicate that doubling times at wave-exposed shorelines are longer than those in wave-sheltered harbours, allowing for more time to adapt to magnified flood hazards. The median doubling times of average water levels including parameterised wave set-up are 1.2 to 5 times longer than those of still water levels as observed with tide gauges. For instantaneous water levels including wave run-up, doubling times are an additional 30% to 100% longer. We conclude that tide gauge-based analyses underestimate adaptation times by underestimating the contribution of waves to extreme water levels, and provide a quantitative framework to guide adaptation policy at wave-exposed shorelines.

Evolution of the Mid-tropospheric Vortex during the Formation of Super Typhoon Megi (2010)

As a follow-up of a previously published article on the contribution of tropical waves, this study explores the evolution of the mid-tropospheric mesoscale cyclonic vortex (MV) during the formation of Typhoon Megi (2010) with a successful cloud-resolving simulation. It is found that the formation and intensification of the MV were related to the deep convection and subsequent stratiform precipitation, while the weakening of the MV was related to the shallow convection. Both the upward transport of vorticity related to the deep convection and the horizontal convergence associated with the stratiform precipitation contributed to the formation and intensification of the MV. Even though the latter was dominant, the former could not be ignored, especially in the early stage of the MV. The MV played dual roles in the formation of Megi. On the one hand, the formation and intensification of MV were primarily associated with the stratiform precipitation, which induced the low-level divergence inhibiting the spin-up of the near-surface cyclonic circulation. On the other hand, the coupled low-level cold core under the MV benefited the accumulation of the convective available potential energy (CAPE), which was favorable for the convective activity. A sensitivity experiment with the evaporative cooling turned off indicated that the development of the MV retarded the genesis process of Megi.

3D full-waveform inversion in time-frequency domain: Field data application

We present a 3D Gauss-Newton full-waveform inversion (3D GN-FWI) method in the time-frequency domain for detection of subsurface anomalies. The use of Gauss-Newton approach is particularly important for near-surface imaging, in which acquired seismic wavefields are dominated by Rayleigh waves. The inverse Hessian matrix utilized in this approach acts as a weighting function to reduce the dominancy of Rayleigh waves, increase the contribution of body waves, and thus help resolve deeper structures. However, the Gauss-Newton method requires a huge amount of computer memory to store derivative wavefields (Jacobian matrix). To address this issue, the presented 3D GN-FWI method exploits advantages of the combined time-frequency domain. Specifically, the forward wave simulation is done in the time-domain to generate wavefields at multiple frequencies simultaneously without the requirement of an inverse matrix solver, while inversion is conducted in the frequency-domain to significantly reduce the required computer memory. Synthetic and field experimental datasets are used to assess the capability of the presented waveform method. The synthetic result shows that a variable profile with a buried void is well recovered. For the field experiment, a large mobile shaker was used to induce wavefields at individual frequencies with consistent magnitudes required for the presented frequency-domain inversion. The wavefields were recorded with uniform 2D grids of sources and receivers on the ground surface, and analyzed to obtain 3D subsurface wave velocity profiles. The seismic result is consistent with an invasive standard penetration test (SPT), including the identified bedrock depth and buried void.