Source author record

Øyvind Breivik

Øyvind Breivik appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

17works
2topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2022arXiv

The Impact of a Reduced High-wind Charnock Parameter on Wave Growth With Application to the North Sea, the Norwegian Sea and the Arctic Ocean

As atmospheric models move to higher resolution and resolve smaller scales, the maximum modeled wind speed also tends to increase. Wave models tuned to coarser wind fields tend to overestimate the wave growth under strong winds. A recently developed semi-empirical parameterization of the Charnock parameter, which controls the roughness length over surface waves, substantially reduces the aerodynamic drag of waves in high winds (above a threshold of 30 m/s). Here we apply the formulation in a recent version of the wave model WAM (Cycle 4.7), which uses a modified version of the physics parameterizations by Ardhuin et al (2010) as well as subgrid obstructions for better performance around complex topography. The new Charnock formulation is tested with wind forcing from NORA3, a recently completed non-hydrostatic atmospheric downscaling of the global reanalysis ERA5 for the North Sea, the Norwegian Sea and the Barents Sea. Such high-resolution atmospheric model integrations tend to have stronger (and more realistic) upper-percentile winds than what is typically found in coarser atmospheric models. A two-year comparison (2011-2012) of a control run against the run with the modified Charnock parameter shows a dramatic reduction of the wave height bias in high-wind cases. The added computational cost of the new physics and the reduction of the Charnock parameter compared to the earlier WAM physics is modest (14%). A longer (1998-2020) hindcast integration with the new Charnock parameter is found to compare well against in situ and altimeter wave measurements both for intermediate and high sea states.

preprint2016arXiv

A Stokes drift approximation based on the Phillips spectrum

A new approximation to the Stokes drift velocity profile based on the exact solution for the Phillips spectrum is explored. The profile is compared with the monochromatic profile and the recently proposed exponential integral profile. ERA-Interim spectra and spectra from a wave buoy in the central North Sea are used to investigate the behaviour of the profile. It is found that the new profile has a much stronger gradient near the surface and lower normalized deviation from the profile computed from the spectra. Based on estimates from two open-ocean locations, an average value has been estimated for a key parameter of the profile. Given this parameter, the profile can be computed from the same two parameters as the monochromatic profile, namely the transport and the surface Stokes drift velocity.

preprint2015arXiv

Surface Wave Effects in the NEMO Ocean Model: Forced and Coupled Experiments

The NEMO general circulation ocean model is extended to incorporate three physical processes related to ocean surface waves, namely the surface stress (modified by growth and dissipation of the oceanic wave field), the turbulent kinetic energy flux from breaking waves, and the Stokes-Coriolis force. Experiments are done with NEMO in ocean-only (forced) mode and coupled to the ECMWF atmospheric and wave models. Ocean-only integrations are forced with fields from the ERA-Interim reanalysis. All three effects are noticeable in the extra-tropics, but the sea-state dependent turbulent kinetic energy flux yields by far the largest difference. This is partly because the control run has too vigorous deep mixing due to an empirical mixing term in NEMO. We investigate the relation between this ad hoc mixing and Langmuir turbulence and find that it is much more effective than the Langmuir parameterization used in NEMO. The biases in sea surface temperature as well as subsurface temperature are reduced, and the total ocean heat content exhibits a trend closer to that observed in a recent ocean reanalysis (ORAS4) when wave effects are included. Seasonal integrations of the coupled atmosphere-wave-ocean model consisting of NEMO, the wave model ECWAM and the atmospheric model of ECMWF similarly show that the sea surface temperature biases are greatly reduced when the mixing is controlled by the sea state and properly weighted by the thickness of the uppermost level of the ocean model. These wave-related physical processes were recently implemented in the operational coupled ensemble forecast system of ECMWF.

preprint2015arXiv

The International Workshop on Wave Hindcasting and Forecasting and the Coastal Hazards Symposium

Following the 13th International Workshop on Wave Hindcasting and Forecasting and 4th Coastal Hazards Symposium in October 2013 in Banff, Canada, a topical collection has appeared in recent issues of Ocean Dynamics. Here we give a brief overview of the history of the conference since its inception in 1986 and of the progress made in the fields of wind-generated ocean waves and the modelling of coastal hazards before we summarize the main results of the papers that have appeared in the topical collection.

preprint2014arXiv

Approximate Stokes Drift Profiles in Deep Water

A deep-water approximation to the Stokes drift velocity profile is explored as an alternative to the monochromatic profile. The alternative profile investigated relies on the same two quantities required for the monochromatic profile, viz the Stokes transport and the surface Stokes drift velocity. Comparisons with parametric spectra and profiles under wave spectra from the ERA-Interim reanalysis and buoy observations reveal much better agreement than the monochromatic profile even for complex sea states. That the profile gives a closer match and a more correct shear has implications for ocean circulation models since the Coriolis-Stokes force depends on the magnitude and direction of the Stokes drift profile and Langmuir turbulence parameterizations depend sensitively on the shear of the profile. The alternative profile comes at no added numerical cost compared to the monochromatic profile.

preprint2014arXiv

Wind and Wave Extremes over the World Oceans from Very Large Ensembles

Global return values of marine wind speed and significant wave height are estimated from very large aggregates of archived ensemble forecasts at +240-h lead time. Long lead time ensures that the forecasts represent independent draws from the model climate. Compared with ERA-Interim, a reanalysis, the ensemble yields higher return estimates for both wind speed and significant wave height. Confidence intervals are much tighter due to the large size of the dataset. The period (9 yrs) is short enough to be considered stationary even with climate change. Furthermore, the ensemble is large enough for non-parametric 100-yr return estimates to be made from order statistics. These direct return estimates compare well with extreme value estimates outside areas with tropical cyclones. Like any method employing modeled fields, it is sensitive to tail biases in the numerical model, but we find that the biases are moderate outside areas with tropical cyclones.

preprint2013arXiv

Wave Extremes in the North East Atlantic from Ensemble Forecasts

A method for estimating return values from ensembles of forecasts at advanced lead times is presented. Return values of significant wave height in the North-East Atlantic, the Norwegian Sea and the North Sea are computed from archived +240-h forecasts of the ECMWF ensemble prediction system (EPS) from 1999 to 2009. We make three assumptions: First, each forecast is representative of a six-hour interval and collectively the data set is then comparable to a time period of 226 years. Second, the model climate matches the observed distribution, which we confirm by comparing with buoy data. Third, the ensemble members are sufficiently uncorrelated to be considered independent realizations of the model climate. We find anomaly correlations of 0.20, but peak events (>P97) are entirely uncorrelated. By comparing return values from individual members with return values of subsamples of the data set we also find that the estimates follow the same distribution and appear unaffected by correlations in the ensemble. The annual mean and variance over the 11-year archived period exhibit no significant departures from stationarity compared with a recent reforecast, i.e., there is no spurious trend due to model upgrades. EPS yields significantly higher return values than ERA-40 and ERA-Interim and is in good agreement with the high-resolution hindcast NORA10, except in the lee of unresolved islands where EPS overestimates and in enclosed seas where it is biased low. Confidence intervals are half the width of those found for ERA-Interim due to the magnitude of the data set.

preprint2013arXiv

Wind and Wave Extremes over the World Oceans From Very Large Forecast Ensembles

Global return value estimates of significant wave height and 10-m neutral wind speed are estimated from very large aggregations of archived ECMWF ensemble forecasts at +240-h lead time from the period 2003-2012. The upper percentiles are found to match ENVISAT wind speed better than ERA-Interim (ERA-I), which tends to be biased low. The return estimates are significantly higher for both wind speed and wave height in the extratropics and the subtropics than what is found from ERA-I, but lower than what is reported by Caires and Sterl (2005) and Vinoth and Young (2011). The highest discrepancies between ERA-I and ENS240 are found in the hurricane-prone areas, suggesting that the ensemble comes closer than ERA-I in capturing the intensity of tropical cyclones. The width of the confidence intervals are typically reduced by 70% due to the size of the data sets. Finally, non-parametric estimates of return values were computed from the tail of the distribution. These direct return estimates compare very well with Generalized Pareto estimates.

preprint2012arXiv

Advances in Search and Rescue at Sea

A topical collection on "Advances in Search and Rescue at Sea" has appeared in recent issues of Ocean Dynamics following the latest in a series of workshops on "Technologies for Search and Rescue and other Emergency Marine Operations" (2004, 2006, 2008 and 2011), hosted by IFREMER in Brest, France. Here we give a brief overview of the history of search and rescue at sea before we summarize the main results of the papers that have appeared in the topical collection. Keywords: Search and rescue (SAR), Trajectory modelling, Stochastic Lagrangian ocean models, Lagrangian measurement methods, ocean surface currents.

preprint2012arXiv

Nearshore wave forecasting and hindcasting by dynamical and statistical downscaling

A high-resolution nested WAM/SWAN wave model suite aimed at rapidly establishing nearshore wave forecasts as well as a climatology and return values of the local wave conditions with Rapid Enviromental Assessment (REA) in mind is described. The system is targeted at regions where local wave growth and partial exposure to complex open-ocean wave conditions makes diagnostic wave modelling difficult. SWAN is set up on 500 m resolution and is nested in a 10 km version of WAM. A model integration of more than one year is carried out to map the spatial distribution of the wave field. The model correlates well with wave buoy observations (0.96) but overestimates the wave height somewhat (18%, bias 0.29 m). To estimate wave height return values a much longer time series is required and running SWAN for such a period is unrealistic in a REA setting. Instead we establish a direction-dependent transfer function between an already existing coarse open-ocean hindcast dataset and the high-resolution nested SWAN model. Return values are estimated using ensemble estimates of two different extreme-value distributions based on the full 52 years of statistically downscaled hindcast data. We find good agreement between downscaled wave height and wave buoy observations. The cost of generating the statistically downscaled hindcast time series is negligible and can be redone for arbitrary locations within the SWAN domain, although the sectors must be carefully chosen for each new location. The method is found to be well suited to rapidly providing detailed wave forecasts as well as hindcasts and return values estimates of partly sheltered coastal regions.

preprint2012arXiv

Real time assimilation of HF radar currents into a coastal ocean model

A real time assimilation and forecasting system for coastal currents is presented. The purpose of the system is to deliver current analyses and forecasts based on assimilation of high frequency radar surface current measurements. The local Vessel Traffic Service monitoring the ship traffic to two oil terminals on the coast of Norway received the analyses and forecasts in real time. A new assimilation method based on optimal interpolation is presented where spatial covariances derived from an ocean model are used instead of simplified mathematical formulations. An array of high frequency radar antennae provide the current measurements. A suite of nested ocean models comprise the model system. The observing system is found to yield good analyses and short range forecasts that are significantly improved compared to a model twin without assimilation. The system is fast; analysis and six hour forecasts are ready at the Vessel Traffic Service 45 minutes after acquisition of radar measurements.

preprint2012arXiv

The Leeway of Shipping Containers at Different Immersion Levels

The leeway of 20-foot containers in typical distress conditions is established through field experiments in a Norwegian fjord and in open-ocean conditions off the coast of France with wind speed ranging from calm to 14 m/s. The experimental setup is described in detail and certain recommendations given for experiments on objects of this size. The results are compared with the leeway of a scaled-down container before the full set of measured leeway characteristics are compared with a semi-analytical model of immersed containers. Our results are broadly consistent with the semi-analytical model, but the model is found to be sensitive to choice of drag coefficient and makes no estimate of the cross-wind leeway of containers. We extend the results from the semi-analytical immersion model by extrapolating the observed leeway divergence and estimates of the experimental uncertainty to various realistic immersion levels. The sensitivity of these leeway estimates at different immersion levels are tested using a stochastic trajectory model. Search areas are found to be sensitive to the exact immersion levels, the choice of drag coefficient and somewhat less sensitive to the inclusion of leeway divergence. We further compare the search areas thus found with a range of trajectories estimated using the semi-analytical model with only perturbations to the immersion level. We find that the search areas calculated without estimates of crosswind leeway and its uncertainty will grossly underestimate the rate of expansion of the search areas. We recommend that stochastic trajectory models of container drift should account for these uncertainties by generating search areas for different immersion levels and with the uncertainties in crosswind and downwind leeway reported from our field experiments.

preprint2011arXiv

A high-resolution hindcast of wind and waves for The North Sea, The Norwegian Sea and The Barents Sea

A combined high-resolution atmospheric downscaling and wave hindcast based on the ERA-40 reanalysis covering the Norwegian Sea, the North Sea and the Barents Sea is presented. The period covered is from September 1957 to August 2002. The dynamic atmospheric downscaling is performed as a series of short prognostic runs initialized from a blend of ERA-40 and the previous prognostic run to preserve the fine-scale surface features from the high-resolution model while maintaining the large-scale synoptic field from ERA-40. The nested WAM wave model hindcast consists of a coarse 50 km model covering the North Atlantic forced with ERA-40 winds and a nested 10-11 km resolution model forced with downscaled winds. A comparison against in situ and satellite observations of wind and sea state reveals significant improvement in mean values and upper percentiles of wind vectors and the significant wave height over ERA-40. Improvement is also found in the mean wave period. ERA-40 is biased low in wind speed and significant wave height, a bias which is not reproduced by the downscaling. The atmospheric downscaling also reproduces polar lows, which can not be resolved by ERA-40, but the lows are too weak and short-lived as the downscaling is not capable of capturing their full life cycle.

preprint2011arXiv

An Operational Search and Rescue Model for the Norwegian Sea and the North Sea

A new operational, ensemble-based search and rescue model for the Norwegian Sea and the North Sea is presented. The stochastic trajectory model computes the net motion of a range of search and rescue objects. A new, robust formulation for the relation between the wind and the motion of the drifting object (termed the leeway of the object) is employed. Empirically derived coefficients for 63 categories of search objects compiled by the US Coast Guard are ingested to estimate the leeway of the drifting objects. A Monte Carlo technique is employed to generate an ensemble that accounts for the uncertainties in forcing fields (wind and current), leeway drift properties, and the initial position of the search object. The ensemble yields an estimate of the time-evolving probability density function of the location of the search object, and its envelope defines the search area. Forcing fields from the operational oceanic and atmospheric forecast system of The Norwegian Meteorological Institute are used as input to the trajectory model. This allows for the first time high-resolution wind and current fields to be used to forecast search areas up to 60 hours into the future. A limited set of field exercises show good agreement between model trajectories, search areas, and observed trajectories for liferafts and other search objects. Comparison with older methods shows that search areas expand much more slowly using the new ensemble method with high resolution forcing fields and the new leeway formulation. It is found that going to higher-order stochastic trajectory models will not significantly improve the forecast skill and the rate of expansion of search areas.

preprint2011arXiv

BAKTRAK: Backtracking drifting objects using an iterative algorithm with a forward trajectory model

The task of determining the origin of a drifting object after it has been located is highly complex due to the uncertainties in drift properties and environmental forcing (wind, waves and surface currents). Usually the origin is inferred by running a trajectory model (stochastic or deterministic) in reverse. However, this approach has some severe drawbacks, most notably the fact that many drifting objects go through nonlinear state changes underway (e.g., evaporating oil or a capsizing lifeboat). This makes it difficult to naively construct a reverse-time trajectory model which realistically predicts the earliest possible time the object may have started drifting. We propose instead a different approach where the original (forward) trajectory model is kept unaltered while an iterative seeding and selection process allows us to retain only those particles that end up within a certain time-space radius of the observation. An iterative refinement process named BAKTRAK is employed where those trajectories that do not make it to the goal are rejected and new trajectories are spawned from successful trajectories. This allows the model to be run in the forward direction to determine the point of origin of a drifting object. The method is demonstrated using the Leeway stochastic trajectory model for drifting objects due to its relative simplicity and the practical importance of being able to identify the origin of drifting objects. However, the methodology is general and even more applicable to oil drift trajectories, drifting ships and hazardous material that exhibit non-linear state changes such as evaporation, chemical weathering, capsizing or swamping. The backtracking method is tested against the drift trajectory of a life raft and is shown to predict closely the initial release position of the raft and its subsequent trajectory.

preprint2011arXiv

Wind-induced drift of objects at sea: the leeway field method

A method for conducting leeway field experiments to establish the drift properties of small objects (0.1-25 m) is described. The objective is to define a standardized and unambiguous procedure for condensing the drift properties down to a set of coefficients that may be incorporated into existing stochastic trajectory forecast models for drifting objects of concern to search and rescue operations and other activities involving vessels lost at sea such as containers with hazardous material. An operational definition of the slip or wind and wave-induced motion of a drifting object relative to the ambient current is proposed. This definition taken together with a strict adherence to 10 m wind speed allows us to refer unambiguously to the leeway of a drifting object. We recommend that all objects if possible be studied using what we term the direct method, where the object's leeway is studied directly using an attached current meter. We divide drifting objects into four categories, depending on their size. For the smaller objects (less than 0.5 m), an indirect method of measuring the object's motion relative to the ambient current must be used. For larger objects, direct measurement of the motion through the near-surface water masses is strongly recommended. Larger objects are categorized according to the ability to attach current meters and wind monitoring systems to them. The leeway field method proposed here is illustrated with results from field work where three objects were studied in their distress configuration; a 1:3.3 sized model of a 40-ft Shipping container, a World War II mine and a 220 l (55-gallon) oil drum.