The recent European mega-heatwaves of 2003 and 2010 broke temperature records across Europe(1-5). Although events of this magnitude were unprecedented from a historical perspective, they are expected to become common by the end of the century(6,7). However, our understanding of extreme heatwave events is limited and their representation in climate models remains imperfect(8). Here we investigate the physical processes underlying recent mega-heatwaves using satellite and balloon measurements of land and atmospheric conditions from the summers of 2003 in France and 2010 in Russia, in combination with a soil-water-atmosphere model. We find that, in both events, persistent atmospheric pressure patterns induced land-atmosphere feedbacks that led to extreme temperatures. During daytime, heat was supplied by large-scale horizontal advection, warming of an increasingly desiccated land surface and enhanced entrainment of warm air into the atmospheric boundary layer. Overnight, the heat generated during the day was preserved in an anomalous kilometres-deep atmospheric layer located several hundred metres above the surface, available to re-enter the atmospheric boundary layer during the next diurnal cycle. This resulted in a progressive accumulation of heat over several days, which enhanced soil desiccation and led to further escalation in air temperatures. Our findings suggest that the extreme temperatures in mega-heatwaves can be explained by the combined multi-day memory of the land surface and the atmospheric boundary layer

Abstract. The Meso-NH Atmospheric Simulation System is a joint e ort of the Centre National de Recherches Me te orologiques and Laboratoire d'Ae rologie. It comprises several elements; a numerical model able to simulate the atmospheric motions, ranging from the large meso-alpha scale down to the micro-scale, with a comprehensive physical package, a¯exible ®le manager, an ensemble of facilities to prepare initial states, either idealized or interpolated from meteorological analyses or forecasts, a¯exible post-processing and graphical facility to visualize the results, and an ensemble of interactive procedures to control these functions. Some of the distinctive features of this ensemble are the following: the model is currently based on the Lipps and Hemler form of the anelastic system, but may evolve towards a more accurate form of the equations system. In the future, it will allow for simultaneous simulation of several scales of motion, by the so-called``interactive grid-nesting technique''. It allows for the in-line computation and accumulation of various terms of the budget of several quantities. It allows for the transport and di usion of passive scalars, to be coupled with a chemical module. It uses the relatively new Fortran 90 compiler. It is tailored to be easily implemented on any UNIX machine. Meso-NH is designed as a research tool for small and meso-scale atmospheric processes. It is freely accessible to the research community, and we have tried to make it as``user-friendly'' as possible, and as general as possible, although these two goals sometimes appear contradictory. The present paper presents a general description of the adiabatic formulation and some of the basic validation simulations. A list of the currently available physical parametrizations and initialization methods is also given. A more precise description of these aspects will be provided in a further paper.

Abstract. The current version of the Dutch AtmosphericLarge-Eddy Simulation (DALES) is presented. DALES is a large-eddy simulation code designed for studies of the physics of the atmospheric boundary layer, including convective and stable boundary layers as well as cloudy boundary layers. In addition, DALES can be used for studies of more specific cases, such as flow over sloping or heterogeneous terrain, and dispersion of inert and chemically active species. This paper contains an extensive description of the physical and numerical formulation of the code, and gives an overview of its applications and accomplishments in recent years.

Abstract. Due to the major role of the sun in heating the earth's surface, the atmospheric planetary boundary layer over land is inherently marked by a diurnal cycle. The afternoon transition, the period of the day that connects the daytime dry convective boundary layer to the night-time stable boundary layer, still has a number of unanswered scientific questions. This phase of the diurnal cycle is challenging from both modelling and observational perspectives: it is transitory, most of the forcings are small or null and the turbulence regime changes from fully convective, close to homogeneous and isotropic, toward a more heterogeneous and intermittent state.These issues motivated the BLLAST (Boundary-Layer Late Afternoon and Sunset Turbulence) field campaign that was conducted from 14 June to 8 July 2011 in southern France, in an area of complex and heterogeneous terrain. A wide range of instrumented platforms including full-size aircraft, remotely piloted aircraft systems, remote-sensing instruments, radiosoundings, tethered balloons, surface flux stations and various meteorological towers were deployed over different surface types. The boundary layer, from the earth's surface to the free troposphere, was probed during the entire day, with a focus and intense observation periods that were conducted from midday until sunset. The BLLAST field campaign also provided an opportunity to test innovative measurement systems, such as new miniaturized sensors, and a new technique for frequent radiosoundings of the low troposphere.Twelve fair weather days displaying various meteorological conditions were extensively documented during the field experiment. The boundary-layer growth varied from one day to another depending on many contributions including stability, advection, subsidence, the state of the previous day's residual layer, as well as local, meso-or synoptic scale conditions.Ground-based measurements combined with tetheredballoon and airborne observations captured the turbulence decay from the surface throughout the whole boundary layer and documented the evolution of the turbulence characteristic length scales during the transition period.Closely integrated with the field experiment, numerical studies are now underway with a complete hierarchy of models to support the data interpretation and improve the model representations.

Abstract. Using a Large-Eddy Simulation model, we have systematically studied the inability of boundary layer turbulence to efficiently mix reactive species. This creates regions where the species are accumulated in a correlated or anti-correlated way, thereby modifying the mean reactivity. We quantify this modification by the intensity of segregation, I S , and analyse the driving mechanisms: heterogeneity of the surface moisture and heat fluxes, various background wind patterns and non-uniform isoprene emissions. The heterogeneous surface conditions are characterized by cool and wet forested patches with high isoprene emissions, alternated with warm and dry patches that represents pasture with relatively low isoprene emissions. For typical conditions in the Amazon rain forest, applying homogeneous surface forcings and in the absence of free tropospheric NO x , the isoprene-OH reaction rate is altered by less than 10 %. This is substantially smaller than the previously assumed I S of 50 % in recent large-scale model analyses of tropical rain forest chemistry. Spatial heterogeneous surface emissions enhance the segregation of species, leading to alterations of the chemical reaction rates up to 20 %. The intensities of segregation are enhanced when the background wind direction is parallel to the borders between the patches and reduced in the case of a perpendicular wind direction. The effects of segregation on trace gas concentrations vary per species. For the highly reactive OH, the differences in concentration averaged over the boundary layer are less than 2 % compared to homogeneous surface conditions, while the isoprene concentration is increased by as much as 12 % due to the reduced chemical reaction rates. These processes take place at the sub-grid scale of chemistry transport models and therefore need to be parameterized.

Ambient total OH reactivity was measured at the Finnish boreal forest station SMEAR II in Hyytiälä (Latitude 61&deg;51' N; Longitude 24&deg;17' E) in July and August 2010 using the Comparative Reactivity Method (CRM). The CRM – total OH reactivity method – is a direct, in-situ determination of the total loss rate of hydroxyl radicals (OH) caused by all reactive species in air. During the intensive field campaign HUMPPA-COPEC 2010 (<b>H</b>yytiälä <b>U</b>nited <b>M</b>easurements of <b>P</b>hotochemistry and <b>P</b>articles in <b>A</b>ir – <b>C</b>omprehensive <b>O</b>rganic <b>P</b>recursor <b>E</b>mission and <b>C</b>oncentration study) the total OH reactivity was monitored both inside (18 m) and directly above the forest canopy (24 m) for the first time. The comparison between these two total OH reactivity measurements, absolute values and the temporal variation have been analyzed here. Stable boundary layer conditions during night and turbulent mixing in the daytime induced low and high short-term variability, respectively. The impact on total OH reactivity from biogenic emissions and associated photochemical products was measured under "normal" and "stressed" (i.e. prolonged high temperature) conditions. The advection of biomass burning emissions to the site caused a marked change in the total OH reactivity vertical profile. By comparing the OH reactivity contribution from individually measured compounds and the directly measured total OH reactivity, the size of any unaccounted for or "missing" sink can be deduced for various atmospheric influences. For "normal" boreal conditions a missing OH reactivity of 58%, whereas for "stressed" boreal conditions a missing OH reactivity of 89% was determined. Various sources of not quantified OH reactive species are proposed as possible explanation for the high missing OH reactivity

No abstract

The Meso-NH Atmospheric Simulation System is a joint eort of the Centre National de Recherches Me te orologiques and Laboratoire d'Ae rologie. It comprises several elements; a numerical model able to simulate the atmospheric motions, ranging from the large meso-alpha scale down to the micro-scale, with a comprehensive physical package, a¯exible ®le manager, an ensemble of facilities to prepare initial states, either idealized or interpolated from meteorological analyses or forecasts, a¯exible post-processing and graphical facility to visualize the results, and an ensemble of interactive procedures to control these functions. Some of the distinctive features of this ensemble are the following: the model is currently based on the Lipps and Hemler form of the anelastic system, but may evolve towards a more accurate form of the equations system. In the future, it will allow for simultaneous simulation of several scales of motion, by the so-called``interactive grid-nesting technique''. It allows for the in-line computation and accumulation of various terms of the budget of several quantities. It allows for the transport and diusion of passive scalars, to be coupled with a chemical module. It uses the relatively new Fortran 90 compiler. It is tailored to be easily implemented on any UNIX machine. Meso-NH is designed as a research tool for small and meso-scale atmospheric processes. It is freely accessible to the research community, and we have tried to make it as``user-friendly'' as possible, and as general as possible, although these two goals sometimes appear contradictory. The present paper presents a general description of the adiabatic formulation and some of the basic validation simulations. A list of the currently available physical parametrizations and initialization methods is also given. A more precise description of these aspects will be provided in a further paper.