Suspended particulate matter (SPM) concentrations in estuaries have been observed to vary strongly over the spring–neap cycle through complex interactions between trapping and re-suspension. However, a systematic framework for analysing the processes causing this spring–neap SPM variability in general is missing. In this study we set up such a framework, consisting of three tiers. First, by studying the sediment transport capacity, how the locations of sediment trapping change over the spring–neap cycle is identified. Second, how the transport capacity affects the sediment stock and bottom pool of sediment is studied. This bottom pool only adapts gradually to the changing transport conditions, incorporating a lag or memory effect. Using a two-timescale analysis it is shown that this slow movement of the bottom pool is the leading source of such lag effects. Third, the SPM concentration is explained from an almost instantaneously balanced exchange between the bottom pool and the water column through re-suspension and deposition.

We demonstrate the use of this framework on two model cases implemented in the idealised width-averaged iFlow model: an idealised test case where the sediment dynamics do not affect the water motion and a case representative of the Loire estuary, with strong feedback between sediment and the water motion through sediment-induced damping of turbulence. The first is illustrative as it allows a full understanding in terms of cause and effect between water motion, transport, and SPM concentration. In the more realistic Loire case, the SPM dynamics cannot be explained in terms of cause and effect but can explain the trapping locations and timing of maximum concentrations in a systematic way in terms of the governing physical mechanisms.@en

This study examines the local, intratidal effects of suspended sediment concentrations (SSCs) on the hydrodynamics and vertical mixing in the Ems Estuary, located on the border between Germany and The Netherlands, during summer and winter seasons when the estuary turbidity maximum (ETM) is located upstream and adjacent to the study site, respectively. Measurements of density, SSCs, turbulent kinetic energy dissipation, and current velocity were collected and analyzed over a semi-diurnal tidal cycle in August of 2018 and January of 2019 as part of the collaborative Ems-Dollard Measurement (EDoM) campaign. During August, the estuary turbidity maximum was located 25 km upstream from the measurement site and local SSCs were low. Results revealed that under these conditions, suspended sediment minimally impacted vertical mixing by stabilizing density near-bottom during flood tide, while typical salinity-induced tidal straining patterns dominated. During January, the ETM was located only 5 km upstream of the measurement site leading to higher local sediment concentrations. Salinity-induced straining of the density occurred on early flood tide, creating stratification that suppressed vertical mixing. The suppression was enhanced by the contribution of vertical gradients in SSC to density, as signified by the gradient Richardson number. Suppression of vertical mixing by sediment-enhanced stratification was most significant within the hour following maximum flood currents when elevated velocity shear occurred. The variability observed between the local dynamics during August and January were attributed to greater sediment concentrations due to the ETM proximity in January. The intratidal asymmetry of vertical mixing observed under higher SSCs likely has implications for sediment transport.

@en

The water motion computed using 3D and 2DH models in tidally dominated shallow waters can, in some cases, differ significantly. In 2DH models, bed friction is typically parametrised in terms of the depth-averaged velocity, whereas in 3D models, typically the near-bed velocity is used. This difference causes the bed shear stress in 2DH models to point towards the depth-averaged velocity, whereas in 3D models, it points towards the near-bed velocity, which are not necessarily the same. Focussing on linearised barotropic models, we derive an exact friction parametrisation for 2DH models such that the same depth-averaged dynamics are described as in the corresponding 3D model. The result is a convolutional friction formulation where the instantaneous friction depends on the present and past velocities, thus modifying the traditional 2DH friction formulation that only depends on the present depth-averaged velocity. In the case of harmonic (tidal) waves, this parametrisation has a clear physical interpretation and shows that the near-bed velocity should be parametrised as a rotated, deformed and phase shifted variant of the depth-averaged velocity. We demonstrate that in certain regions of the parameter space, it may be impossible to calibrate a 2DH model that uses a traditional friction law to reproduce the water levels from a 3D model, showing that the 3D friction formulation can be crucial to capture the 3D dynamics within a depth-averaged model. This phenomenon is explored in detail in a narrow well-mixed estuary.

@en

A nonlinear Harmonic Navier-Stokes (HNS) framework is introduced for simulating instabilities in laminar spanwise-invariant shear layers, featuring sharp and smooth wall surface protuberances. While such cases play a critical role in the process of laminar-to-turbulent transition, classical stability theory analyses such as parabolized or local stability methods fail to provide (accurate) results, due to their underlying assumptions. The generalized incompressible Navier-Stokes (NS) equations are expanded in perturbed form, using a spanwise and temporal Fourier ansatz for flow perturbations. The resulting equations are discretized using spectral collocation in the wall-normal direction and finite-difference methods in the streamwise direction. The equations are then solved using a direct sparse-matrix solver. The nonlinear mode interaction terms are converged iteratively. The solution implementation makes use of a generalized domain transformation to account for geometrical smooth surface features, such as humps. No-slip conditions can be embedded in the interior domain to account for the presence of sharp surface features such as forward- or backward-facing steps. Common difficulties with Navier-Stokes solvers, such as the treatment of the outflow boundary and convergence of nonlinear terms, are considered in detail. The performance of the developed solver is evaluated against several cases of representative boundary layer instability growth, including linear and nonlinear growth of Tollmien-Schlichting waves in a Blasius boundary layer and stationary crossflow instabilities in a swept flat-plate boundary layer. The latter problem is also treated in the presence of a geometrical smooth hump and a sharp forward-facing step at the wall. HNS simulation results, such as perturbation amplitudes, growth rates, and shape functions, are compared to benchmark flow stability analysis methods such as Parabolized Stability Equations (PSE), Adaptive Harmonic Linearized Navier-Stokes (AHLNS), or Direct Numerical Simulations (DNS). Good agreement is observed in all cases. The HNS solver is subjected to a grid convergence study and a simple performance benchmark, namely memory usage and computational cost. The computational cost is found to be considerably lower than high-fidelity DNS at comparable grid resolutions.@en

Tidal basins are connected to the outer sea by one or more tidal inlets. In this presentation, we focus on so-called double inlet systems, i.e., tidal inlet systems with two connecting channels. A typical example is the Marsdiep-Vlie system, located in the Dutch Wadden Sea. Using an idealized morphodynamic model, we aim at directly finding morphodynamic equilibria of double inlet systems and assess their stability. For unstable configurations, the linear stability mechanism will be presented.@en

An extensive field campaign, the Ems-Dollard Measurements (EDoM), was executed in the Ems Estuary, bordering the Netherlands and Germany, aimed at better understanding the mechanisms that drive the exchange of water and sediments between a relatively exposed outer estuary and a hyper-turbid tidal river. More specifically, the reasons for the large up-estuary sediment accumulation rates and the role of the tidal river on the turbidity in the outer estuary were insufficiently understood. The campaign was designed to unravel the hydrodynamic and sedimentary exchange mechanisms, comprising two hydrographic surveys during contrasting environmental conditions using eight concurrently operating ships and 10 moorings measuring for at least one spring–neap tidal cycle. All survey locations were equipped with sensors measuring flow velocity, salinity, and turbidity (and with stationary ship surveys taking water samples), while some of the survey ships also measured turbulence and sediment settling properties. These observations have provided important new insights into horizontal sediment fluxes and density-driven exchange flows, both laterally and longitudinally. An integral analysis of these observations suggests that large-scale residual transport is surprisingly similar during periods of high and low discharge, with higher river discharge resulting in both higher seaward-directed fluxes near the surface and landward-directed fluxes near the bed. Sediment exchange seems to be strongly influenced by a previously undocumented lateral circulation cell driving residual transport. Vertical density-driven flows in the outer estuary are influenced by variations in river discharge, with a near-bed landward flow being most pronounced in the days following a period with elevated river discharge. The study site is more turbid during winter conditions, when the estuarine turbidity maximum (ETM) is pushed seaward by river flow, resulting in a more pronounced impact of suspended sediments on hydrodynamics. All data collected during the EDoM campaign, but also standard monitoring data (waves, water levels, discharge, turbidity, and salinity) collected by Dutch and German authorities are made publicly available at 4TU Centre for Research Data (https://doi.org/10.4121/c.6056564.v3; van Maren et al., 2022).@en

Structure and intensity of estuarine exchange flow depend significantly on the eddy viscosity A_v profile which is dynamically linked to various forces (e.g., gravitational, tidal, wind-driven). The impact of winds on the exchange flow is complex due to its direct (local and remote changes in shear and density stratification) and indirect (modifications to A_v profiles) contributions. This study aims (i) to include wind entrainment effects in the tidally averaged A_v parameterization; (ii) to develop an analytical one-dimensional model for the wind driven exchange flow by using this novel parameterization and assess the tidally averaged dynamics over a relevant physical parameter-space, subdomains of which have not yet been explored numerically. This one-dimensional model is based on a balance between frictional forces and pressure gradient, calibrated with a tidally-resolving one-dimensional water-column model with second-moment closure. Structure and intensity of the resulting exchange flow profiles are analyzed with respect to three dimensionless parameters (the unsteadiness of boundary layer mixing U_n, scaled-directional wind stress W, and horizontal stratification S_i). While down-estuarine winds enhance the gravitational circulation, up-estuarine winds result in either a two-layer inverted circulation opposing the gravitational circulation, or a three-layer flow (favored by relatively strong S_i, weak W, and moderate U_n) that is up-estuarine at the surface with classical two-layer circulation underneath. Relative thickness of surface and bottom boundary layers affect both the intensity and the inflection depth of the exchange flow layers. Up-estuarine winds with W≳0.5 yield unstable stratification and reduce the exchange flow intensity with increasing W.

@en

An idealized width-averaged model is employed to study the influence of wind stress on subtidal salt intrusion and stratification in well-mixed and partially stratified estuaries. We show that even in mild conditions, wind forcing can influence the estuarine salinity structure in a substantial way. By studying the role of wind forcing on dominant salt transport balances and associated salt transport regimes, we unify and clarify ambiguous observations from previous authors regarding the influence of wind stress: the response of the estuarine salinity structure to wind forcing is different depending on the underlying dominant salt transport balance, which in turn was found to determine whether wind-induced salinity shear, wind-induced modulation of the longitudinal salt distribution, or wind-induced mixing dominates. SIGNIFICANCE STATEMENT: The purpose of this idealized study is to better understand how wind influences the salinity distribution in estuaries on large time scales. This is important because a change in winds can move saline water further inland, threatening freshwater availability and the natural balance of delicate ecosystems. We clarify the sometimes ambiguous observations regarding the influence of wind on the salt distribution and highlight the importance of including average wind forcing in analyses of estuarine dynamics on large time scales.

@en

Tidally averaged transport of salt in estuaries is controlled by various subtidal and tidal processes. In this study, we show the relative importance of various subtidal and tidal transport processes in a width-averaged sense. This is done for a large range of forcing and geometric parameters, which describe well-mixed to salt wedge estuaries. To this end, we develop a width-averaged process-based model aimed at conducting and analyzing a large number of experiments (∼40,000). We find that the salt transport is dominated by one of seven salt transport balances, or regimes. Four of these regimes are dominated by subtidal processes, while the other three are dominated by tidal processes. Which regime occurs in a part of an estuary depends on four dimensionless parameters, representing local geometry, and forcing conditions. One of the regimes features salt import by correlations between the depth-averaged tidal velocity and salinity. While this mechanism was previously only associated with along-channel geometric variations, we find it can also be a dominant mechanism in a significant part of the parameter space due to river-induced tidal asymmetry, independent of river geometry. We apply our classification to a case study of part of the Dutch Rhine delta and compare to decomposition results of a fully realistic three-dimensional model. We find the estuary features two regimes, with import dominated by subtidal shear transport in the seaward part of the estuary and by depth-averaged tidal correlations in the landward part of the estuary.

@en

Channel–shoal patterns are often observed in the back–barrier basins of inlet systems and are important from both an economical and ecological point of view. Focussing on double–inlet systems, the initial formation of these patterns is investigated using an idealized model. The model is governed by the depth–averaged shallow water equations, a depth–integrated concentration equation and a tidally–averaged bottom evolution equation. Focussing on rectangular basins and neglecting the effects of earth rotation, it is found that laterally uniform morphodynamic equilibria can become linearly unstable, resulting in initial patterns that resemble channels and shoals. When the water motion is only forced by an M₂ tidal constituent, the existence of (laterally uniform) morphodynamic equilibria for which both inlets are connected strongly depends on the relative phase and amplitudes of the tidal forcing. If such equilibria exist, they can be either stable against small perturbations or linearly unstable. If these equilibria are linearly unstable, two instability mechanisms can be identified, the first related to the convergences and divergences of diffusive transports, the second mechanism related to a combination of advective and diffusive transports. In the former case, all eigenvalues are real and the bedforms grow exponentially in time. In the latter case, the eigenvalues are complex, resulting in bedforms that both migrate and grow in time. In case external overtides and a time–independent discharge are included, no diffusive instabilities are found anymore for the parameters considered in this paper. This implies that all instabilities are migrating in time. In all cases considered, the bed perturbations have only an appreciable amplitude at locations where the underlying laterally uniform equilibrium has a local minimum in water depth. This is consistent with observations from numerical models and laboratory experiments.

@en

The salinity structure in estuaries is classically described in terms of the salinity structure as well mixed, partially mixed, or salt wedge. The existing knowledge about the processes that result in such salinity structures comes from highly idealized models that are restricted to either well-mixed and partially mixed cases or subtidal salt wedge estuaries. Hence, there is still little knowledge about the processes driving transitions between these different salinity structures and the estuarine parameters at which such a transition is found. As an important step toward a unified description of the dominant processes driving well-mixed, partially mixed, and salt wedge estuaries, a subtidal width-averaged model applicable to all these salinity structures is developed and systematically analyzed. Using our model, we identify four salinity regimes, resulting from different balances of dominant processes. It is shown that each regime is uniquely determined by two dimensionless parameters: an estuarine Froude and Rayleigh number, representing freshwater discharge and tidal mixing, respectively, resulting in a classification of the regimes in terms of these two parameters. Furthermore, analytical expressions to approximate the salt intrusion length in each regime are developed. These expressions are used to illustrate that the salt intrusion length in different regimes responds in a highly different manner to changes in depth and freshwater discharge. As one of the key results, we show that there are only very weak relations between the process-based regime of an estuary and the salt intrusion length and top-bottom stratification. This implies that the salinity structure of an estuary cannot be uniquely matched to a regime.

@en

The classification diagram developed by Hansen and Rattray is one of the classic papers on classification of estuarine salinity dynamics. However, we found several inconsistencies in both their stratification–circulation and estuarine classification diagrams. These findings considerably change the interpretation of their work. Furthermore, while their classification includes salt wedge estuaries, the model used to derive this is only applicable to well-mixed and partially mixed estuaries. Here, we identify and solve these inconsistencies, and we propose new adjusted and extended stratification–circulation and classification diagrams. To this end, we summarize the model of Hansen and Rattray and extend their work to find analytical model solutions and an adjusted stratification–circulation diagram. Using this new diagram, it is shown that Hansen and Rattray incorrectly discussed the behavior of dispersion-dominated estuaries and that several parts of the diagram correspond to physically unrealistic model solutions. This is then used to demonstrate that several estuarine classes identified by Hansen and Rattray correspond to physically unrealistic model solutions and can therefore not be interpreted. A new and extended classification is proposed by using a recently developed model that extends the work of Hansen and Rattray to salt wedge estuaries. This results in an extended estuarine classification including examples of the location of 12 estuaries in this new diagram.@en

A new parameterization for the estuarine turbulent eddy viscosity coefficient is developed considering the influence of wind forcing and feedback between stratification and shear. The emerging tidally averaged eddy viscosity profile A_v^T is parameterized as parabolic under well-mixed conditions, and is composed of a skewed-Gaussian-like form for the upper layer, and a parabolic form for the bottom layer under stratified conditions. The precise shape of the profiles depends parametrically on the bottom boundary layer thickness, the bulk Richardson number, and the Wedderburn number. The parameterized A_v^T profiles show excellent agreement with profiles obtained from numerical models. To explore the importance of vertically varying A_v^T with regard to exchange processes, an analytical model is designed. This one-dimensional model is based on a balance between frictional forces and pressure gradient. The resulting exchange flow is analyzed over the relevant parameter space that is associated with horizontal and vertical stratification through the bulk Richardson number, and the bi-directional wind stress via the Wedderburn number. Down-estuary wind enhances the up-estuary flow near the bottom and down-estuary flow near the surface driving an exchange flow pattern typically associated with gravitational circulation. Up-estuary wind results in either a two-layer inverted circulation opposing the gravitational circulation, or a three-layer flow that is up-estuarine at the surface with classical two-layer circulation underneath. Three-layer flow emerges with a weak wind. With increasing runoff velocity, three-layer flow transitions to a single layer flow under weak stratification conditions.

@en

Asymmetric tidal turbulence (ATT) strongly influences estuarine health and functioning. However, its impact on the three-dimensional estuarine dynamics and the feedback of water motion and salinity distribution on ATT remain poorly understood, especially for short estuaries (estuarine length ≪ tidal wavelength). This study systematically investigates the abovementioned interactions in a short estuary for the first time, considering periodically weakly stratified conditions. This is done by developing a three-dimensional semi-analytical model (combining perturbation method with finite element method) that allows a dissection of the contributions of different processes to ATT, estuarine circulation, and salt transport. The generation of ATT is dominated by (i) strain-induced periodic stratification and (ii) asymmetric bottomshear-generated turbulence, and their contributions to ATT are different both in amplitude and phase. The magnitude of the residual circulation related to ATT and the eddy viscosity–shear covariance (ESCO) is about half of that of the gravitational circulation (GC) and shows a ‘‘reversed’’ pattern as compared to GC. ATT generated by strain-induced periodic stratification contributes to an ESCO circulation with a spatial structure similar to GC. This circulation reduces the longitudinal salinity gradients and thus weakens GC. Contrastingly, the ESCO circulation due to asymmetric bottom-sheargenerated turbulence shows patterns opposite to GC and acts to enhance GC. Concerning the salinity dynamics at steady state, GC and tidal pumping are equally important to salt import, whereas ESCO circulation yields a significant seaward salt transport. These findings highlight the importance of identifying the sources of ATT to understand its impact on estuarine circulation and salt distribution.

@en

Many estuaries exhibit seasonality in the estuary-scale distribution of suspended particulate matter (SPM). This SPM distribution depends on various factors, including freshwater discharge, salinity intrusion, erodibility, and the ability of cohesive SPM to flocculate into larger aggregates. Various authors indicate that biotic factors, such as the presence of algae and their excretion of sticky transparent exopolymer particles (TEP), affect the flocculation and erosion processes. Consequently, seasonality in these biotic factors may play a role in the observed seasonality in SPM. Whereas the impact of abiotic factors on seasonality in SPM is well studied, the relative contribution of biotically induced seasonality is largely unknown. In this study, we employ two approaches to assess the aggregated importance of biotically induced seasonality in flocculation and erosion on seasonality in SPM in the Scheldt estuary. In the first approach, we focus on seasonality of in situ observations in the Scheldt estuary of turbidity, floc size, Chlorophyll-a, and TEP, showing that the abiotic parameters show seasonality, while seasonality in TEP is ambiguous. The second approach concerns a reverse engineering method to calibrate biotically affected parameters of a coupled sediment transport-flocculation model to turbidity observations, allowing us to compare the modeled SPM concentrations to the observations. Driven by seasonality in freshwater discharge, the model captures the observed seasonality in SPM without requiring biotically induced seasonality in flocculation and erosion, which is supported by the absence of seasonality in TEP.

@en

In finite element methods, the accuracy of the solution cannot increase indefinitely since the round-off error related to limited computer precision increases when the number of degrees of freedom (DoFs) is large enough. Because a priori information of the highest attainable accuracy is of great interest, we construct an innovative method to obtain the highest attainable accuracy given the order of the elements. In this method, the truncation error is extrapolated when it converges at the asymptotic rate, and the bound of the round-off error follows from a generically valid error estimate, obtained and validated through extensive numerical experiments. The highest attainable accuracy is obtained by minimizing the sum of these two types of errors. We validate this method using a one-dimensional Helmholtz equation in space. It shows that the highest attainable accuracy can be accurately predicted, and the CPU time required is much smaller compared with that using successive grid refinement.

@en

Climate and human pressures can influence the evolution of estuarine sediment dynamics concurrently, but the understanding and quantification of their cause–effect relationships are still challenging due to the occurrence of complex hydro-morpho-sedimentary feedbacks. The Garonne Tidal River (GTR, upper Gironde Estuary, France) is a clear example of a system stressed by both anthropogenic and climate change, as it has been subject to decreasing river discharges, natural morphological changes, and gravel extraction. To understand the relative effect of each hydrological and geomorphological pressure on the turbidity maximum zone (TMZ), the sediment dynamics in the GTR over the last six decades was evaluated using the width-averaged idealized iFlow model. Model results show a gradual increase in tidal amplitude and currents over the decades that has led to the upstream shift of the landward sediment-transport capacity components (external M₄ tide, spatial settling lag, and tidal return flow). The upstream displacement of the TMZ between the 1950s and the 2010s was estimated to be at least 19 km, of which about three fourth was induced by geomorphological changes and one fourth by hydrological changes. Concerning the geomorphological changes, the natural evolution of the lower Gironde morphology was the main pressure inducing the displacement of the TMZ in the GTR. Anthropogenic and natural changes in morphology and bed roughness in the GTR itself also contributed to this evolution. The natural geomorphological changes were, in turn, probably promoted by the evolution of sediment dynamics, so this study reveals the closed circle that governs the intensification of the TMZ.

@en

The existence of morphodynamic equilibria of double-inlet systems is investigated using a cross-sectionally averaged morphodynamic model. The number of possible equilibria and their stability strongly depend on the forcing conditions and geometry considered. This is illustrated by considering a rectangular double-inlet system forced by M₂ tidal constituents only. Depending on the M₂ amplitudes and phases at both entrances, no equilibrium, one equilibrium or multiple morphodynamic equilibria may exist. In case no equilibrium is found, the minimum water depth becomes zero somewhere in the system, reducing the double-inlet system to two single-inlet systems. In the other cases, the location of the minimum water depth and the direction of the tidally-averaged sediment transport, as well as their actual values, depend strongly on the M₂ tidal characteristics. Such parameter sensitivity is also observed when including the residual and M₄ forcing contributions to the water motion, and when allowing for width variations. This suggests that, when considering a specific system, the number and stability of morphodynamic equilibria, as well as the characteristics of these quantities, can only be assessed by investigating that specific system in detail. As an example, the Marsdiep-Vlie inlet system in the Dutch Wadden Sea is considered. It is found that, by using parameter values and a geometry characteristic for this system, the water motion and bathymetry in morphodynamic equilibrium are qualitatively reproduced. Also the direction and order of magnitude of the tidally-averaged suspended sediment transport compare well with those obtained from a high-complexity numerical model.

@en

A depth-averaged (2DH) exploratory model is developed to identify morphodynamic equilibria in short mesotidal inlet systems with arbitrary planform geometries. The water motion is forced by an M ₂ tidal constituent at the seaward entrance and is described by the depth-averaged shallow water equations, whereas the depth-integrated suspended sediment concentration follows from a diffusion equation, taking into account local inertia, horizontal eddy diffusion and topographically induced diffusive effects, erosion, and deposition. Based on a scaling analysis, it follows that the fine sandy bed evolution is dominated by the depth-integrated diffusive sediment transport. The depth-integrated advective contributions are one order smaller and therefore neglected. This observation also allows for the neglect of the advective terms in the governing equations. The associated morphodynamic equilibria are directly identified based on a continuation approach. By means of the exploratory model, the morphodynamic equilibria are studied in basins with a planform geometry characterized by width variations as a function of the distance to the seaward boundary. The model results show that in the case of a sufficient degree of widening in the landward direction, the equilibrium bed level exhibits significant lateral structures, characterized by shallow zones and deeper channels. The first channel bifurcation, as observed in many short tidal inlet systems, is forced by the planform geometry of the basin, and the associated physical mechanisms are explained. Furthermore, two mechanisms inducing asymmetric morphodynamic equilibria are investigated, of which the effect of an asymmetric basin planform seems to be dominant over that of the Coriolis force.

@en

The initial stage of the laminar–turbulent transition of semi-infinite flows can be characterized as either an absolute or convective instability, naturally associated with localized wave packets. A convective instability is directly linked to an absolute instability in a different reference frame. Therefore, our aim is to determine the absolute stability of a flow in a given but arbitrary reference frame, which can only be directly inferred from the absolute eigenvalue spectrum. If advective processes are present, the associated absolute eigenfunctions grow exponentially in space in the advective direction. The eigenvalue spectrum is usually computed numerically, which requires truncating the domain and prescribing artificial boundary conditions at these truncation boundaries. For separated boundary conditions, the resulting spectrum approaches the absolute spectrum as the domain length tends to infinity. Since advective processes result in spatially exponentially growing eigenfunctions, it becomes increasingly difficult to represent these functions numerically as the domain length increases. Hence, a naive numerical implementation of the eigenvalue problem may result in a computed spectrum that strongly deviates from the (mathematically correct) absolute spectrum due to numerical errors. To overcome these numerical inaccuracies, we employ a weighted method ensuring the convergence to the absolute spectrum. From a physical point of view, this method removes the advection-induced spatial exponential growth from the eigenfunctions. The resulting (absolute) spectrum allows for a direct interpretation of the character of the pertinent perturbations and the eigensolutions can be used to construct and analyse the evolution of localized wave packets in an efficient way.@en