The parameterised description of subgrid-scale processes in the clear and cloudy boundary layer has a strong impact on the performance skill in any numerical weather prediction (NWP) or climate model and is still a prime source of uncertainty. Yet, improvement of this parameterised description is hard because operational models are highly optimised and contain numerous compensating errors. Therefore, improvement of a single parameterised aspect of the boundary layer often results in an overall deterioration of the model as a whole. In this paper, we will describe a comprehensive integral revision of three parameterisation schemes in the High Resolution Local Area Modelling - Aire Limitée Adaptation dynamique Développement InterNational (HIRLAM-ALADIN) Research on Mesoscale Operational NWP In Europe - Applications of Research to Operations at Mesoscale (HARMONIE-AROME) model that together parameterise the boundary layer processes: the cloud scheme, the turbulence scheme, and the shallow cumulus convection scheme. One of the major motivations for this revision is the poor representation of low clouds in the current model cycle. The newly revised parametric descriptions provide an improved prediction not only of low clouds but also of precipitation. Both improvements can be related to a stronger accumulation of moisture under the atmospheric inversion. The three improved parameterisation schemes are included in a recent update of the HARMONIE-AROME configuration, but its description and the insights in the underlying physical processes are of more general interest as the schemes are based on commonly applied frameworks. Moreover, this work offers an interesting look behind the scenes of how parameterisation development requires an integral approach and a delicate balance between physical realism and pragmatism.

@en

Four fluoride fuel salt samples (78LiF-22ThF₄) in graphite crucibles were irradiated in the HFR Petten for a duration of 508 Full Power Days under the name SALIENT-01 (SALt Irradiation ExperimeNT). Goal of the experiment was to gain experience with the design of liquid salt experiments and the handling of the salts before and after irradiation. Specific research goals for SALIENT-01 are (i) to confirm claims of good fission product retention in the salt, (ii) to obtain size distributions for noble metal particles using Transmission Electron Microscopy and (iii) to assess possible interactions between fuel salt and fine-grained nuclear graphite, as well as possible uptake of fission products by the graphite. Here the design and irradiation history of the experiment are discussed together with plans for post-irradiation examinations. Limitations in representativeness of this experiment and capsule irradiations in general are discussed as well as follow-up actions to improve the quality of future irradiations.

@en

NRG together with JRC Karlsruhe has set out to perform a series of molten fuel salt irradiations in the High Flux Reactor (HFR) Petten, to build up irradiation experience with molten salt samples, the handling of irradiated salt and the treatment of salt waste produced by these irradiations. As a first step, an irradiation of small 78LiF-22ThF₄ fuel samples in graphite crucibles is being conducted under the name SALIENT-01 (SALt Irradiation ExperimeNT). The specific goals for SALIENT-01 are (i) to confirm claims of good fission product retention in the salt by post-irradiation Knudsen-cell effusion by fission product release measurements, in particular with respect to Cs and I which dominate the radioactive release during thermal transient and fuel failure in conventional reactor fuel, as for example occurred in the Fukushima accident, (ii) to obtain size distributions for noble metal particles using Transmission Electron Microscopy and (iii) to assess possible interactions between fuel salt and fine-grained nuclear graphite, as well as possible uptake of fission products by the graphite. The SALIENT-01 irradiation has started in August 2017 and is projected to finish in August 2019. In this contribution the design, irradiation history and status of the experiment are treated, as well as plans for post-irradiation examinations. The limitations of the experiment will also be discussed, as well as the follow-up actions taken to assess its representativeness and to improve the quality of future molten salt irradiations.

@en

High-resolution large-eddy simulations of the Antarctic very stable boundary layer reveal a mechanism for systematic and periodic intermittent bursting. A nonbursting state with a boundary layer height of just 3 m is alternated by a bursting state with a height of ≈5 m. The bursts result from unstable wave growth triggered by a shear-generated Kelvin–Helmholtz instability, as confirmed by linear stability analysis. The shear at the top of the boundary layer is built up by two processes. The upper, quasi-laminar layer accelerates due to the combined effect of the pressure force and rotation by the Coriolis force, while the lower layer decelerates by turbulent friction. During the burst, this shear is eroded and the initial cause of the instability is removed. Subsequently, the interfacial shear builds up again, causing the entire sequence to repeat itself with a time scale of ≈10 min. Despite the clear intermittent bursting, the overall change of the mean wind profile is remarkably small during the cycle. This enables such a fast erosion and recovery of the shear. This mechanism for cyclic bursting is remarkably similar to the mechanism hypothesized by Businger in 1973, with one key difference. Whereas Businger proposes that the flow acceleration in the upper layer results from downward turbulent transfer of high-momentum flow, the current results indicate no turbulent activity in the upper layer, hence requiring another source of momentum. Finally, it would be interesting to construct a climatology of shear-generated intermittency in relation to large-scale conditions to assess the generality of this Businger mechanism.@en

Wind machines are used in the agricultural sector to prevent or mitigate the adverse effects of night frost in spring. In this study we aim to quantify the impact of wind machine operation on the local temperature field in an orchard. To this end, a field experiment is conducted and experimental analysis is combined with numerical simulation studies in order to assess the functional relations between wind machine performance and the dominating physical processes occurring during radiative frost events. Experimental observations showed that the temperature response strongly depends on the radial distance to the fan and the height above the surface. In agreement with previous studies, the wind machine was able to achieve rotation-averaged temperature increases of up to 50% of the inversion strength ( ≈ 3 K) in an area of 3–5 ha at 1 m height. Furthermore, it was observed that even weak ambient winds (<1 m/s) already may cause strong upwind-downwind asymmetries in the protected area, the downwind area being larger. The numerical model, inspired by the field experiment, showed similar spatial temperature responses as compared to observations. Interestingly, it was found that slower rotation times of the wind machine (3 to 6 min) lead to a significant increase of affected area, while the temperature enhancement itself stayed relatively constant. Variation of the horizontal tilt angle showed that, in our model, temperature enhancement was maximized between 8^∘ and 16^∘. This nearly horizontal flow already facilitates efficient vertical mixing of momentum and heat, presumably due to generation of shear instabilities at the lower edge of the jet. Finally, like in the observations also the numerical result showed strong upwind-downwind asymmetry in the affected area due to background wind.

@en

An overview is given of 50-year Cabauw observations and research on the structure and dynamics of the atmospheric boundary layer. It is shown that over time this research site with its 200-m meteorological tower has grown into an atmospheric observatory with a comprehensive observational program encompassing almost all aspects of the atmospheric column including its boundary conditions. This is accomplished by the Cabauw Experimental Site for Atmospheric Research (CESAR) a consortium of research institutes. CESAR plays an important role in the educational programs of the CESAR universities. The current boundary-layer observational program is described in detail, and other parts of the CESAR observational program discussed more briefly. Due to an open data policy the CESAR datasets are used by researchers all over the world. Examples are given of the use of the long time series for model evaluation, satellite validation, and process studies. The role of tall towers is discussed in relation to the development of more and better ground-based remote sensing techniques. CESAR is now incorporated into the Ruisdael observatory, the large-scale atmospheric research infrastructure in the Netherlands. With Ruisdael the embedding of the Dutch atmospheric community in national policy landscape, and in the European atmospheric research infrastructures is assured for the coming decade.@en

We present an extension of the dynamic global vegetation model, Lund-Potsdam-Jena Managed Land (LPJmL), to simulate planted forests intended for carbon (C) sequestration. We implemented three functional types to simulate plantation trees in temperate, tropical, and boreal climates. The parameters of these functional types were optimized to fit target growth curves (TGCs). These curves represent the evolution of stemwood C over time in typical productive plantations and were derived by combining field observations and LPJmL estimates for equivalent natural forests. While the calibrated model underestimates stemwood C growth rates compared to the TGCs, it represents substantial improvement over using natural forests to represent afforestation. Based on a simulation experiment in which we compared global natural forest versus global forest plantation, we found that forest plantations allow for much larger C uptake rates on the timescale of 100 years, with a maximum difference of a factor of 1.9, around 54 years. In subsequent simulations for an ambitious but realistic scenario in which 650Mha (14% of global managed land, 4.5% of global land surface) are converted to forest over 85 years, we found that natural forests take up 37PgC versus 48PgC for forest plantations. Comparing these results to estimations of C sequestration required to achieve the 2°C climate target, we conclude that afforestation can offer a substantial contribution to climate mitigation. Full evaluation of afforestation as a climate change mitigation strategy requires an integrated assessment which considers all relevant aspects, including costs, biodiversity, and trade-offs with other land-use types. Our extended version of LPJmL can contribute to such an assessment by providing improved estimates of C uptake rates by forest plantations.

@en

Conventional in situ observations of meteorological variables are restricted to a limited number of levels near the surface, with the lowest observation often made around 1-m height. This can result in missed observations of both shallow fog, and the initial growth stage of thicker fog layers. At the same time, numerical experiments have demonstrated the need for high vertical grid resolution in the near-surface layer to accurately simulate the onset of fog; this requires correspondingly high-resolution observational data for validation. A two-week field campaign was conducted in November 2017 at the Cabauw Experimental Site for Atmospheric Research (CESAR) in the Netherlands. The aim was to observe the growth of shallow fog layers and assess the possibility of obtaining very high-resolution observations near the surface during fog events. Temperature and relative humidity were measured at centimetre resolution in the lowest 7 m using distributed temperature sensing. Further, a novel approach was employed to estimate visibility in the lowest 2.5 m using a camera and an extended light source. These observations were supplemented by the existing conventional sensors at the site, including those along a 200-m tall tower. Comparison between the increased-resolution observations and their conventional counterparts show the errors to be small, giving confidence in the reliability of the techniques. The increased resolution of the observations subsequently allows for detailed investigations of fog growth and evolution. This includes the observation of large temperature inversions in the lowest metre (up to 5 K) and corresponding regions of (super)saturation where the fog formed. Throughout the two-week observation period, fog was observed twice at the conventional sensor height of 2.0 m. Two additional low-visibility events were observed in the lowest 0–0.5 m using the camera-based observations, but were missed by the conventional sensors. The camera observations also showed the growth of shallow radiation fog, forming in the lowest 0.5 m as early as two hours before it was observed at the conventional height of 2 m.@en

Conventional in situ observations of visibility and other meteorological variables are restricted to a limited number of heights near the surface, with the lowest observation often made above 1 m. This can result in missed observations of shallow fog as well as the initial growth stage of thicker fog layers. At the same time, numerical experiments have demonstrated the need for high vertical grid resolution in the near-surface layer to accurately simulate the onset of fog; this requires correspondingly high-resolution observational data for validation. In November 2017, a field experiment was conducted at the Cabauw Experimental Site for Atmospheric Research (CESAR) in the Netherlands with the aim of observing the growth of shallow fog from the ground up, assessing the applicability of emerging high-resolution methods for observing shallow fog. Two innovative, high-resolution techniques were employed: distributed temperature sensing (DTS), providing temperature and relative humidity observations at vertical resolutions as fine as 1 cm, and a novel camera-LED method to observe near-surface visibility below the conventional sensor height of 2.0 m. These observations were supplemented by the existing observations at the site, including those along a 200-m tall tower. Comparison between the high-resolution observations and their conventional counterparts shows the errors to be small, giving confidence to the reliability of the techniques. The high resolution of the observations subse- quently allows for detailed investigations of near surface processes. The growth of fog layers from the ground up was observed with very strong temperature inversions in the lowest metre (up to 5 K), and corresponding region of (super)saturation where the fog formed and grew. Throughout the two-week observation period, fog was observed twice at the conventional sensor height of 2.0 m, but up to four times in the lowest 0-0.5 m using the camera estimates, with the shallow fog also forming up to two hours before it was observed by the conventional sensor. The observations are supplemented by high-resolution numerical simulations of the experimental period, highlighting the sensitivity of the fog layer to surface properties and ambient conditions, providing greater insight into what drives the growth of a very shallow fog layer (i.e. < 1 m) into a deeper, and therefore more dangerous, layer.@en

In this work we study the dynamics of the surface-based temperature inversion over the Antarctic Plateau during the polar winter. Using 6 years of observations from the French–Italian Antarctic station Concordia at Dome C, we investigate sudden regime transitions in the strength of the near-surface temperature inversion. Here we define “near-surface” as being within the domain of the 45-m measuring tower. In particular, we consider the strongly nonlinear relation between the 10-m inversion strength (T _10m – T _s ) and the 10-m wind speed. To this end, all individual events for which the 10-m inversion strength increases or decreases continuously by more than 15 K in time are considered. Composite time series and vertical profiles of wind and temperature reveal specific characteristics of the transition from weak to very strong inversions and vice versa. In contrast to midlatitudes, the largest variations in temperature are not found at the surface but at a height of 10 m. A similar analysis was performed on results from an atmospheric single-column model (SCM). Overall, the SCM results reproduce the observed characteristics of the transitions in the near-surface inversion remarkably well. Using model output, the underlying mechanisms of the regime transitions are identified. The nonlinear relation between inversion strength and wind speed at a given level is explained by variations in the geostrophic wind speed, changes in the depth of the turbulent layer and the vertical divergence of turbulent fluxes. Moreover, the transitions between different boundary layer regimes cannot be explained without considering the contribution of subsidence heating.

@en

Observations of two typical contrasting weakly stable and very stable boundary layers from the winter at Dome C station, Antarctica, are used as a benchmark for two centimetre-scale-resolution large-eddy simulations. By taking the Antarctic winter, the effects of the diurnal cycle are eliminated, enabling the study of the long-lived steady stable boundary layer. With its homogeneous, flat snow surface, and extreme stabilities, the location is a natural laboratory for studies on the long-lived stable boundary layer. The two simulations differ only in the imposed geostrophic wind speed, which is identified as the main deciding factor for the resulting regime. In general, a good correspondence is found between the observed and simulated profiles of mean wind speed and temperature. Discrepancies in the temperature profiles are likely due to the exclusion of radiative transfer in the current simulations. The extreme stabilities result in a considerable contrast between the stable boundary layer at the Dome C site and that found at typical mid-latitudes. The boundary-layer height is found to range from approximately 50m to just 5m in the most extreme case. Remarkably, heating of the boundary layer by subsidence may result in thermal equilibrium of the boundary layer in which the associated heating is balanced by the turbulent cooling towards the surface. Using centimetre-scale resolutions, accurate large-eddy simulations of the extreme stabilities encountered in Antarctica appear to be possible. However, future simulations should aim to include radiative transfer and sub-surface heat transport to increase the degree of realism of these types of simulations.@en

This poster will be presented at EMS 2019 in Copenhagen. Fog - in particular, the associated reduction is visibility - presents a hazard to airport operations. Although tech- nology has improved to allow greater safety during fog events, protocol still requires more time between aircraft movements, often resulting in significant delays and cancellations. Yet, observations at Amsterdam’s Schiphol International Airport in the Netherlands (one of Europe’s busiest airports) suggest that the airport buildings and aircraft operations themselves may help to alleviate some of the fog hazard. Meteorological data from from a network of weather stations at and around Schiphol airport are used to assess the occurrence and severity of fog events. Runways located closer to airport terminals are shown to experience both fewer and shorter fog events than those at greater distances from the main terminal complex. Further, large aircraft - in particular Boeing 747s - are observed to raise local runway temperature and wind speed by as much as a few degrees under nocturnal stable boundary layer conditions, which can have a significant impact on the formation of local fog. These find- ings present an interesting look at the local influence of airport construction and aircraft operations, suggesting that they may ultimately lead to "built-in" fog mitigation beyond the natural, undisturbed state. The benefit likely grows with the size of the airport. The larger and busier the airport (such as Heathrow, or Paris-Charles de Gaulle), for example, the more damaging a fog event can be, but the greater potential for disruption of fog formation due to greater number of aircraft movements.@en

We present a conceptual model for the diurnal cycle of the dry atmospheric boundary layer (ABL). It may serve as a framework for future numerical studies on the transitional dynamics that characterize the ABL over land. The conceptual model enables us to define expressions for relevant physical scales as a function of the most prominent forcing parameters and the low degree of complexity facilitates a dimensionless description. This is useful to help generalize boundary layer dynamics that occur on a diurnal time scale. Further, the model’s application for numerical studies is illustrated herein with two examples: a single-column-model study that assesses the effect of wind forcing on the main characteristics of the diurnal cycle, and a large-eddy-simulation study on the daily evolution of turbulence under weak-wind-forcing conditions. The results from these studies sketch the general evolution of the present set of diurnal-cycle systems in more detail. We discuss how the setups are able to reproduce well-known dynamical features of the ABL and also highlight limitations, where the simple conceptual system is unable to describe realistic ABL behavior. We conclude that the present conceptual model has an interesting balance between model-system complexity and physical realism, such that it is useful for future idealized studies on the diurnal cycle of the ABL.

@en

The Netherlands is characterized by highly variable land use within a small area, and a strong influence of the North Sea on national climate. Devoid of significant topography, it is an excellent location for assessing the relative influence of various factors on fog occurrence in the absence of terrain effects. Using observations from a dense network of weather stations throughout the country, the climatology of fog in the Netherlands is assessed over a period of 45 years. On a national scale, interannual variability is linked to changes in synoptic pressure-gradient forcing. Within the country, a comprehensive in-depth analysis of regional differences between fog occurrence is made, together with an assessment of local physical factors which could bias fog formation in one location over another. Regional variability is shown to be strongly related to the mesoscale influences of urbanization and the North Sea. In fact, some locations experience over twice as much fog as others. From this finding, a simple index is presented, which combines the water and urban fraction surrounding a station. This “Regionally Weighted Index” (RWI) is able to accurately sort the stations according to their relative fogginess. Its practical use is encouraged for assessing a given site's climatological favourability, even when in situ meteorological observations are unavailable.

@en

The reduction in visibility that accompanies fog events presents a hazard to human safety and navigation. However, accurate fog prediction remains elusive, with numerical methods often unable to capture the conditions of fog formation, and observational methods having high false-alarm rates in order to obtain high hit rates of prediction. In this work, 5 years of observations from the Cabauw Experimental Site for Atmospheric Research are used to further investigate how false alarms may be reduced using the statistical method for diagnosing radiation-fog events from observations developed by Menut et al. (Boundary-Layer Meteorol 150:277–297, 2014). The method is assessed for forecast lead times of 1–6 h and implementing four optimization schemes to tune the prediction for different needs, compromising between confidence and risk. Prediction scores improve significantly with decreased lead time, with the possibility of achieving a hit rate of over 90% and a false-alarm rate of just 13%. In total, a further 31 combinations of predictive variables beyond the original combination are explored (including mostly, e.g., variables related to moisture and static stability of the boundary layer). Little change to the prediction scores indicates any appropriate combination of variables that measure saturation, turbulence, and near-surface cooling can be used. The remaining false-alarm periods are manually assessed, identifying the lack of spatio–temporal information (such as the temporal evolution of the local conditions and the advective history of the airmass) as the ultimate limiting factor in the methodology’s predictive capabilities. Future observational studies are recommended that investigate the near-surface evolution of fog and the role of non-local heterogeneity on fog formation.

@en

Investigation of meteorological measurements along a 45 m tower at Dome C on the high East Antarctic Plateau revealed two distinct stable boundary layer (SBL) regimes at this location. The first regime is characterized by strong winds and continuous turbulence. It results in full vertical coupling of temperature, wind magnitude and wind direction in the SBL. The second regime is characterized by weak winds, associated with weak turbulent activity and very strong temperature inversions reaching up to 25 K in the lowest 10 m. Vertical temperature profiles are generally exponentially shaped (convex) in the first regime and ‘convex–concave–convex’ in the second. The transition between the two regimes is particularly abrupt when looking at the near-surface temperature inversion and it can be identified by a 10 m wind-speed threshold. With winds under this threshold, the turbulent heat supply toward the surface becomes significantly lower than the net surface radiative cooling. The threshold value (including its range of uncertainty) appears to agree with recent theoretical predictions from the so-called ‘minimum wind speed for sustainable turbulence’ (MWST) theory. For the quasi-steady, clear-sky winter cases, the relation between the near-surface inversion amplitude and the wind speed takes a characteristic ‘S’ shape. Closer analysis suggests that this relation corresponds to a ‘critical transition’ between a steady turbulent and a steady ‘radiative’ regime, with a dynamically unstable branch in the transition zone. These fascinating characteristics of the Antarctic boundary layer challenge present and future numerical models to represent this region in a physically correct manner.

@en

A conceptual model is used in combination with observational analysis to understand regime transitions of near-surface temperature inversions at night as well as in Arctic conditions. The model combines a surface energy budget with a bulk parameterization for turbulent heat transport. Energy fluxes or feedbacks due to soil and radiative heat transfer are accounted for by a "lumped parameter closure," which represents the "coupling strength" of the system. Observations from Cabauw, Netherlands, and Dome C, Antarctica, are analyzed. As expected, inversions are weak for strong winds, whereas large inversions are found under weak-wind conditions. However, a sharp transition is found between those regimes, as it occurs within a narrow wind range. This results in a typical S-shaped dependency. The conceptual model explains why this characteristic must be a robust feature. Differences between the Cabauw and Dome C cases are explained from differences in coupling strength (being weaker in the Antarctic). For comparison, a realistic column model is run. As findings are similar to the simple model and the observational analysis, it suggests generality of the results. Theoretical analysis reveals that, in the transition zone near the critical wind speed, the response time of the system to perturbations becomes large. As resilience to perturbations becomes weaker, it may explain why, within this wind regime, an increase of scatter is found. Finally, the so-called heat flux duality paradox is analyzed. It is explained why numerical simulations with prescribed surface fluxes show a dynamical response different from more realistic surface-coupled systems.

@en

Geostrophic wind speed data, derived from pressure observations, are used in combination with tower measurements to investigate the nocturnal stable boundary layer at Cabauw, the Netherlands. Since the geostrophic wind speed is not directly influenced by local nocturnal stability, it may be regarded as an external forcing parameter of the nocturnal stable boundary layer. This is in contrast to local parameters such as in situ wind speed, the Monin-Obukhov stability parameter (z/L), or the local Richardson number. To characterize the stable boundary layer, ensemble averages of clear-sky nights with similar geostrophic wind speeds are formed. In this manner, the mean dynamical behavior of near-surface turbulent characteristics and composite profiles of wind and temperature are systematically investigated. The classification is found to result in a gradual ordering of the diagnosed variables in terms of the geostrophic wind speed. In an ensemble sense the transition from the weakly stable to very stable boundary layer is more gradual than expected. Interestingly, for very weak geostrophic winds, turbulent activity is found to be negligibly small while the resulting boundary cooling stays finite. Realistic numerical simulations for those cases should therefore have a comprehensive description of other thermodynamic processes such as soil heat conduction and radiative transfer.

@en

The performance of an atmospheric single-column model (SCM) is studied systematically for stably-stratified conditions. To this end, 11 years (2005–2015) of daily SCM simulations were compared to observations from the Cabauw observatory, The Netherlands. Each individual clear-sky night was classified in terms of the ambient geostrophic wind speed with a (Formula presented.) bin-width. Nights with overcast conditions were filtered out by selecting only those nights with an average net radiation of less than (Formula presented.). A similar procedure was applied to the observational dataset. A comparison of observed and modelled ensemble-averaged profiles of wind speed and potential temperature and time series of turbulent fluxes showed that the model represents the dynamics of the nocturnal boundary layer (NBL) at Cabauw very well for a broad range of mechanical forcing conditions. No obvious difference in model performance was found between near-neutral and strongly-stratified conditions. Furthermore, observed NBL regime transitions are represented in a natural way. The reference model version performs much better than a model version that applies excessive vertical mixing as is done in several (global) operational models. Model sensitivity runs showed that for weak-wind conditions the inversion strength depends much more on details of the land-atmosphere coupling than on the turbulent mixing. The presented results indicate that in principle the physical parametrizations of large-scale atmospheric models are sufficiently equipped for modelling stably-stratified conditions for a wide range of forcing conditions.

@en

This study aims to find the typical growth rate of the temperature inversion during the onset of the stable boundary layer around sunset.The sunset transition is a very challenging period for numerical weather prediction, since neither accepted theories for the convective boundary layer nor those for the stable boundary layer appear to be applicable. To gainmore insight in this period, a systematic investigation of the temperature inversion growth rate is conducted. A statistical procedure is used to analyze almost 16 years of observations from the Cabauw observational tower, supported by observations from two additional sites (DomeCandKarlsruhe). The results show that, on average, the growth rate of the temperature inversion (normalized by the maximum inversion during the night) weakly declines with increasing wind speed. The observed growth rate is quantitatively consistent among the sites, and it appears insensitive to various other parameters. The results were also insensitive to the afternoon decay rate of the net radiation except when this decay rate was very weak. These observations are compared to numerical solutions of three models with increasing complexity: a bulk model, an idealized single-column model (SCM), and an operational-level SCM. It appears only the latter could reproduce qualitative features of the observations using a first-order closure. Moreover, replacing this closure with a prognostic TKE scheme substantially improved the quantitative performance. This suggests that idealized models assuming instantaneous equilibrium flux-profile relations may not aid in understanding this period, since history effects may qualitatively affect the dynamics.

@en