ALBERT

Description: Droughts are widespread disasters worldwide and are concurrently influenced by multiple large-scale climate signals. This is particularly true over Japan, where drought has strong heterogeneity due to multiple factors such as monsoon, topography, and ocean circulations. Regional heterogeneity poses challenges for drought prediction and management. To overcome this difficulty, this study provides a comprehensive analysis of teleconnection between climate signals and homogeneous drought zones over Japan. First, droughts are characterized by simulated soil moisture from land surface model during 1958-2012. The Mclust toolkit, distinct empirical orthogonal function, and wavelet coherence analysis are used, respectively, to investigate the homogeneous drought zone, principal component of each homogeneous zone, and teleconnection between climate signals and drought. Results indicate that nine homogeneous drought zones with different characteristics are defined and quantified. Among these nine zones, zone-1 is dominated by extreme drought events. Zone-2 and zone-6 are typical representatives of spring droughts, while zone-7 is wet for most of the period. The Hokkaido region is divided into wetter zone-4 and drier zone-9. Zone-3, zone-5 and zone-8 are distinguished by the topography. The analyses also reveal almost nine zones have a high level of homogeneity, with more than 60% explained variance. Also, these nine zones are dominated by different large-scale climate signals: the Arctic Oscillation has the strongest impact on zone-1, zone-7, and zone-8; the influence of the North Atlantic Oscillation on zone-3, zone-4, and zone-6 is significant; zone-2 and zone-9 are both dominated by the Pacific Decadal Oscillation; El Niño-Southern Oscillation dominates zone-5. The results will be valuable for drought management and drought prevention.

Description: The all-sky assimilation of radiances from microwave instruments is developed in the 4D-EnVar analysis system at Environment and Climate Change Canada (ECCC). Assimilation of cloud-affected radiances from Advanced Microwave Sounding Unit A (AMSUA) temperature sounding channels 4 and 5 for non-precipitating scenes over the ocean surface is the focus of this study. Cloud-affected radiances are discarded in the ECCC operational data assimilation system due to the limitations of forecast model physics, radiative transfer models, and the strong non-linearity of the observation operator. In addition to using symmetric estimate of innovation standard deviation for quality control, a state-dependent observation error inflation is employed at the analysis stage. The background state clouds are scaled by a factor of 0.5 to compensate for a systematic overestimation by the forecast model, before being used in the observation operator. The changes in the fit of the background state to observations show mixed results. The number of AMSUA channels 4 and 5 assimilated observations in the all-sky experiment is 5-12% higher than in the operational system. The all-sky approach improves temperature analysis when verified against ECMWF operational analysis in the areas where the extra cloud-affected observations were assimilated. Statistically significant reductions in error standard deviation by 1-4% for the analysis and forecasts of temperature, specific humidity, and horizontal wind speed up to maximum 4 days were achieved in the all-sky experiment in the lower troposphere. These improvements result mainly from the use of cloud information for computing the observation-minus-background departures. The operational implementation of all-sky assimilation is planned for Fall 2021.

Description: The reported decreasing trend of the annual tropical cyclone (TC) landfalls in southern China and increasing trend in southeastern China in recent decades are confirmed to be an abrupt shift occurring at the end of the twentieth century, based on a statistical analysis. The opposite trends in the two adjacent regions are often considered to be a result of tropical cyclone landfalls in southern China being deflected northward. However, it is demonstrated in this study that they are phenomenally independent. In fact, the abrupt decrease of TC landfalls in southern China occurs as a result of an abrupt decrease of the westward events in the postpeak season (October–December), which in turn is a consequence of a significant decrease of the TC genesis frequency in the southeastern part of the western North Pacific (WNP) Ocean basin. On the other hand, the abrupt increase of TC landfalls in southeastern China occurs because of an abrupt increase of the northwest events in the peak season (July–September), as the consequence of a statistically westward shift of TC genesis. The relevant variations of TC genesis are shown to be mainly caused by decreased relative vorticity and increased vertical wind shear, which, however, are intrinsically related to the accelerated zonal atmospheric circulation driven by a La Niña–like sea surface warming pattern over the WNP that developed after the end of twentieth century.

Description: The reproducibility of precipitation in the early stages of forecasts, often called a spin-down or spin-up problem, has been a significant issue in numerical weather prediction. This problem is caused by moisture imbalance in the analysis data, and in the case of the Japan Meteorological Agency’s (JMA’s) mesoscale data assimilation system JNoVA, we found that the imbalance stems from the existence of unrealistic supersaturated states in the minimal solution of the cost function in JNoVA. Based on the theory of constrained optimization problems, we implemented an exterior penalty function method for the mixing ratio within JNoVA to suppress unrealistic supersaturated states. The advantage of this method is the simplicity of its theory and implementation. The results of twin data assimilation cycle experiments conducted for the Heavy Rain Event of July 2018 over Japan show that—with the new method—unrealistic supersaturated states are reduced successfully, negative temperature bias to the observations is alleviated, and a sharper distribution of the mixing ratio is obtained. These changes help to initiate the development of convection at the proper location and improve the fractions skill score (FSS) of precipitation in the early stages of the forecast. From these results, we conclude that the initial shock caused by moisture imbalance is mitigated by implementing the penalty function method, and the new moisture balance has a positive impact on the reproducibility of precipitation in the early stages of forecasts.

Description: As a key to modulate the negative feedback to tropical cyclone (TC) intensity, the TC-induced inner-core sea surface cooling (SSCIC) is poorly understood. Using a linear two-layer theory and OGCM experiments, this study illustrates that the pattern of the inner-core mixing can be well interpreted by the wind-driven currents in the mixed layer (ML). This interpretation is based on: 1) the mixing is triggered by the ML bulk shear instability; 2) the lag of upwelling makes the inner-core bulk shear equivalent to the inner-core wind-driven currents. Overall, the patterns of the inner-core bulk shear and mixing resemble the crescent body of a sickle. As an accumulative result of mixing, the SSCIC is clearly weaker than the maximum cold wake because of the weaker mixing ahead of the inner core and nearly zero mixing in a part of the inner core. The SSCIC induced by a rectilinear-track TC is mainly dominated by the inner-core mixing. Only for a slow-moving case, upwelling and horizontal advection can make minor contributions to the SSCIC by incorporating them with mixing. The SSCIC strength is inversely proportional to the moving speed, suggesting the mixing time rather than the mixing strength dominates the SSCIC. Despite inability in treating the mixing strength, this study elucidates the fundamental dynamical mechanisms of SSCIC, especially emphasizes the different roles of mixing, upwelling and horizontal advection for fast- and slow-moving TCs, and thus provides a good start point to understand SSCIC.

Description: Diurnal variation in surface latent heat flux (LHF) and the effects of diurnal variations in LHF-related variables on the climatological LHF are examined using observations from the Global Tropical Moored Buoy Array. The estimated amplitude of the climatological diurnal LHF over the Indo-Pacific warm pool and the equatorial Pacific and Atlantic cold tongues is remarkable, with maximum values exceeding 20.0 W m−2. Diurnal variability of sea surface skin temperature (SSTskin) is the primary contributor to the diurnal LHF amplitude. Because the diurnal SSTskin amplitude has an inverse relationship with surface wind speed over the tropical oceans, an inverse spatial pattern between the diurnal LHF amplitude and surface wind speed results. Resolving diurnal variations in the SSTskin and wind improves the estimate of the climatological LHF by properly capturing the daytime SSTskin and daily mean wind speed, respectively. The diurnal SSTskin-associated contribution is large over the warm pool and equatorial cold tongues where low wind speeds tend to cause strong diurnal SSTskin warming, while the magnitude associated with the diurnal winds is large over the highly dynamic environment of the Inter-Tropical Convergence Zone. The total diurnal contribution is about 9.0 W m−2 on average over the buoy sites. There appears to be a power function (linear) relationship between the diurnal SSTskin-associated (wind-associated) contribution and surface mean wind speed (wind speed enhancement from diurnal variability). The total contribution from diurnal variability can be estimated accurately from high-frequency surface wind measurements using these relationships.

Description: Accurate representation of stratospheric trace gas transport is important for ozone modeling and climate projection. Intermodel spread can arise from differences in the representation of transport by the diabatic (overturning) circulation vs. comparatively faster adiabatic mixing by breaking waves, or through numerical errors, primarily diffusion. This study investigates the impact of these processes on transport using an idealised tracer, the age-of-air. Transport is assessed in two state-of-the-art dynamical cores based on fundamentally different numerical formulations: finite volume and spectral element. Integrating the models in free-running and nudged tropical wind configurations reveals the crucial impact of tropical dynamics on stratospheric transport. Using age-budget theory, vertical and horizontal gradients of age allow comparison of the roles of the diabatic circulation, adiabatic mixing, and the numerical diffusive flux. Their respective contribution is quantified by connecting the full 3-d model to the tropical leaky pipe framework of Neu and Plumb (1999). Transport by the two cores varies significantly in the free-running integrations, with the age in the middle stratosphere differing by about 2 years primarily due to differences in adiabatic mixing. When winds in the tropics are constrained, the difference in age drops to about 0.5 years; in this configuration, more than half the difference is due to the representation of the diabatic circulation. Numerical diffusion is very sensitive to the resolution of the core, but does not play a significant role in differences between the cores when they are run at comparable resolution. It is concluded that fundamental differences rooted in dynamical core formulation can account for a substantial fraction of transport bias between climate models.

Description: This paper investigates the value of weather and climate information at different timescales for decision making in the Tanzanian disaster risk reduction sector using non-monetary approaches. Interviews and surveys were conducted with institutions responsible for disaster management at national, regional and district level. A range of values were identified including: 1) making informed decisions for disaster preparedness, response, recovery and restoration related activities; 2) tailoring of directives and actions based on sectoral impacts; 3) identification of hotspot areas for diseases outbreaks and surplus food production. However, while, a number of guidelines, policies, acts and regulations for disaster risk reduction exist it is not clear how well they promote the use of weather and climate information across climate sensitive sectors. Nonetheless, we find that well-structured disaster risk reduction coordination across sectors and institutions from the national to district level exists, although there is a need for further development of integrated Early Warning Systems, and a common platform to evaluate effectiveness and usefulness of weather warnings and advisories. Key challenges to address in increasing the uptake of weather warnings and advisories include language barriers, limited dissemination to rural areas, and limited awareness of forecasts. Based on the findings of this study, we recommend further quantitative evaluation of the skill of the severe weather warnings issued by the Tanzania Meteorological Authority, and an assessment of how decisions and actions are made by recipients of the warnings in the disaster risk reduction sector at different stages in the warning, response and recovery process.

Description: This paper reports the assimilation of cloud optical depth datasets into a variational data assimilation system to improve cloud ice, cloud water, rain, snow, and graupel analysis in extreme weather events for improving forecasts. A cloud optical depth forward operator was developed and implemented in the Space and Time Multiscale Analysis System (STMAS), a multiscale three-dimensional variational analysis system. Using this improved analysis system, the NOAA GOES-15 DCOMP (Daytime Cloud Optical and Microphysical Properties) cloud optical depth products were assimilated to improve the microphysical states. For an eight-day period of extreme weather events in September 2013 in Colorado, the United States, the impact of the cloud optical depth assimilation on the analysis results and forecasts was evaluated. The DCOMP products improved the cloud ice and cloud water predictions significantly in convective and lower levels. The DCOMP products also reduced errors in temperature and relative humidity data at the top (250–150 hPa) and bottom (850–700 hPa) layers. With the cloud ice improvement at higher layers, the DCOMP products provided better forecasts of cloud liquid at low layers (900–700 hPa), temperature and wind at all layers, and relative humidity at middle and bottom layers. Furthermore, for this extreme weather event, both equitable threat score (ETS) and bias were improved throughout the 12 h period, with the most significant improvement observed in the first 3 h. This study will raise the expectation of cloud optical depth product assimilation in operational applications.

Description: Complex-terrain locations often have repeatable near-surface wind patterns, such as synoptic gap flows and local thermally forced flows. An example is the Columbia River Valley in east-central Oregon-Washington, a significant wind-energy-generation region and the site of the Second Wind-Forecast Improvement Project (WFIP2). Data from three Doppler lidars deployed during WFIP2 define and characterize summertime wind regimes and their large-scale contexts, and provide insight into NWP model errors by examining differences in the ability of a model [NOAA’s High-Resolution Rapid-Refresh (HRRR-version1)] to forecast wind-speed profiles for different regimes. Seven regimes were identified based on daily time series of the lidar-measured rotor-layer winds, which then suggested two broad categories. First, in three regimes the primary dynamic forcing was the large-scale pressure gradient. Second, in two regimes the dominant forcing was the diurnal heating-cooling cycle (regional sea-breeze-type dynamics), including the marine intrusion previously described, which generates strong nocturnal winds over the region. The other two included a hybrid regime and a non-conforming regime. For the large-scale pressure-gradient regimes, HRRR had wind-speed biases of ~1 m s−1 and RMSEs of 2-3 m s−1. Errors were much larger for the thermally forced regimes, owing to the premature demise of the strong nocturnal flow in HRRR. Thus, the more dominant the role of surface heating in generating the flow, the larger the errors. Major errors could result from surface heating of the atmosphere, boundary-layer responses to that heating, and associated terrain interactions. Measurement/modeling research programs should be aimed at determining which modeled processes produce the largest errors, so those processes can be improved and errors reduced.

Description: Despite an increased understanding of environments favorable for tornadic supercells, it is still sometimes unknown why one favorable environment produces many long-tracked tornadic supercells and another seemingly equally-favorable environment produces only short-lived supercells. One relatively unexplored environmental parameter that may differ between such environments is the degree of backing or veering of the midlevel shear vector, especially considering that such variations may not be captured by traditional supercell or tornado forecast parameters. We investigate the impact of the 3-6 km shear vector orientation on simulated supercell evolution by systematically varying it across a suite of idealized simulations. We found that the orientation of the 3-6 km shear vector dictates where precipitation loading is maximized in the storms, and thus alters the storm-relative location of downdrafts and outflow surges. When the shear vector is backed, outflow surges generally occur northwest of an updraft, produce greater convergence beneath the updraft, and do not disrupt inflow, meaning that the storm is more likely to persist and produce more tornado-like vortices (TLVs). When the shear vector is veered, outflow surges generally occur north of an updraft, produce less convergence beneath the updraft, and sometimes undercut it with outflow, causing it to tilt at low levels, sometimes leading to storm dissipation. These storms are shorter lived and thus also produce fewer TLVs. Our simulations indicate that the relative orientation of the 3-6 km shear vector may impact supercell longevity and hence the time period over which tornadoes may form.

Description: This study investigates the stratosphere-troposphere coupling associated with the Scandinavian (SCA) pattern in boreal winter. The results indicate that the SCA impacts stratospheric circulation but that its positive and negative phases have different effects. The positive phase of the SCA (SCA+) pattern is restricted to the troposphere, but the negative phase (SCA−) extends to the upper stratosphere. The asymmetry between phases is also visible in the lead-lag evolution of the stratosphere and troposphere. Prominent stratospheric anomalies are found to be intensified following SCA+ events, but prior to SCA− events. Further analysis reveals that the responses are associated with upward propagation of planetary waves, especially wavenumber 1 which is asymmetric between SCA phases. The wave amplitudes in the stratosphere, originating from the troposphere, are enhanced after the SCA+ events and before the SCA− events. Furthermore, the anomalous planetary wave activity can be understood through its interference with climatological stationary waves. Constructive wave interference is accompanied by clear upward propagation in the SCA+ events, while destructive interference suppresses stratospheric waves in the SCA− events. Our results also reveal that the SCA+ events are more likely to be followed by sudden stratospheric warming (SSW) events, because of the deceleration of stratospheric westerlies following the SCA+ events.

Description: The Subantarctic Mode Water (SAMW) is a major water mass in the South Indian and Pacific oceans and plays an important role in the ocean uptake and anthropogenic heat and carbon. The characteristics, formation, and long-term evolution of the SAMW are investigated in the “historical” and “SSP245” scenario simulations of the sixth Coupled Models Intercomparison Project (CMIP6). Defined by the low potential vorticity, the simulated SAMW is consistently thinner, shallower, lighter, and warmer than in observations, due to biases in the winter mixed layer properties and spatial distribution. The biases are especially large in the South Pacific Ocean. The winter mixed layer bias can be attributed to unrealistic heat loss and stratification in the models. Nevertheless, the SAMW is presented better in the CMIP6 than CMIP5, regarding its volume, location, and physical characteristics. In warmer climate, the simulated SAMW in the South Indian Ocean consistently becomes lighter in density, with a reduced volume and a southward shift in the subduction region. The reduced heat loss, instead of the increased Ekman pumping induced by the poleward intensified westerly wind, dominates in the SAMW change. The winter mixed layer shoals in the northern outcrop region and the SAMW subduction shifts southward where the mixed layer remains deep. The projected reduction of the SAMW volume is likely to impact the heat and freshwater redistribution in the Southern Ocean.

Description: Future projections of precipitation change over tropical land are often enhanced by vegetation responses to CO2 forcing in Earth system models. Projected decreases in rainfall over the Amazon basin and increases over the Maritime Continent are both stronger when plant physiological changes are modeled than if these changes are neglected, but the reasons for this amplification remain unclear. The responses of vegetation to increasing CO2 levels are complex and uncertain, including possible decreases in stomatal conductance and increases in leaf area index due to CO2 fertilization. Our results from an idealized atmospheric general circulation model show that the amplification of rainfall changes occurs even when we use a simplified vegetation parameterization based solely on CO2-driven decreases in stomatal conductance, indicating that this mechanism plays a key role in complex model projections. Based on simulations with rectangular continents we find that reducing terrestrial evaporation to zero with increasing CO2 notably leads to enhanced rainfall over a narrow island. Strong heating and ascent over the island trigger moisture advection from the surrounding ocean. In contrast, over larger continents rainfall depends on continental evaporation. Simulations with two rectangular continents representing South America and Africa reveal that the stronger decrease in rainfall over the Amazon basin seen in Earth system models is due to a combination of local and remote effects, which are fundamentally connected to South America’s size and its location with respect to Africa. The response of tropical rainfall to changes in evapotranspiration is thus connected to size and configuration of the continents.

Description: Modeling studies have shown that surface air temperature (SAT) increase in response to an increase in the atmospheric CO2 concentration is larger over land than over ocean. This so-called land–ocean warming contrast, φ, defined as the land–mean SAT change divided by the ocean-mean SAT change, is a striking feature of global warming. Small heat capacity over land is unlikely the sole cause because the land-ocean warming contrast is found in the equilibrium state of CO2 doubling experiments.Several different mechanisms have been proposed to explain the land–ocean warming contrast, but the comprehensive understanding has not yet been obtained. In Part I of this study, we propose a framework to diagnose φ based on energy budgets at the top of atmosphere and for the atmosphere, which enables the decomposition of contributions from effective radiative forcing (ERF), climate feedback, heat capacity, and atmospheric energy transport anomaly to φ. Using this framework, we analyzed the SAT response to an abrupt CO2 quadrupling using 15 Coupled Model Intercomparison Project Phase 6 (CMIP6) Earth system models. In the near-equilibrium state (years 121-150), φ is 1.49 ± 0.11, which is primarily induced by the land–ocean difference in ERF and heat capacity. We found that contributions from ERF, feedback, and energy transport anomaly tend to cancel each other, leading to a small inter-model spread of φ compared to the large spread of individual components. In the equilibrium state without heat capacity contribution, ERF and energy transport anomaly are the major contributors to φ, which shows a weak negative correlation with the equilibrium climate sensitivity.

Description: Based on observational data analyses and idealized modeling experiments, we investigated the distinctive impacts of central Pacific (CP-) El Niño and eastern Pacific (EP-) El Niño on the Antarctic sea ice concentration (SIC) in austral spring (September to November). The tropical heat sources associated with EP-El Niño and the co-occurred positive phase of Indian Ocean Dipole (IOD) excite two branches of Rossby wave trains that propagate southeastward, causing an anomalous anticyclone over the eastern Ross-Amundsen-Bellingshausen Seas. Anomalous northerly (southerly) wind west (east) of the anomalous anticyclone favor poleward (offshore) movements of sea ice, resulting in a sea ice loss (growth) in the eastern Ross-Amundsen Seas (the Bellingshausen-Weddell Seas). Meanwhile, the anomalous northerly (southerly) wind also advected warmer and wetter (colder and drier) air into the eastern Ross-Amundsen Seas (the Bellingshausen-Weddell Seas), causing surface warming (cooling) through the enhanced (reduced) surface heat fluxes and thus contributing to the sea ice melting (growth). CP-El Niño, however, forces a Rossby wave train that generates an anomalous anticyclone in the eastern Ross-Amundsen Seas, 20° west of that caused by EP-El Niño. Consequently, a positive SIC anomaly occurs in the Bellingshausen Sea. A dry version of the Princeton atmospheric general circulation model was applied to verify the roles of anomalous heating in the tropics. The result showed that EP-El Niño can remotely induce an anomalous anticyclone and associated dipole temperature pattern in the Antarctic region, whereas CP-El Niño generates a similar anticyclone pattern with its location shift westward by 20° in longitudes.

Description: The NOAA National Water Model (NWM), maintained and executed by the NOAA National Weather Service (NWS) Office of Water Prediction, provides operational hydrological guidance throughout the Contiguous United States. Based on the WRF-Hydro model architecture developed by the National Center for Atmospheric Research (NCAR), the NWM was recently modified for semi-arid domains, by permitting it to explicitly resolve infiltration from ephemeral channels into the underlying channel bed as an added model sink term. To analyze the added value of channel infiltration in semi-arid environments, we calibrated NWM v2.1 (with the channel infiltration function) to 56 independent basins in the western CONUS, following identical calibration methods as the pre-operational NWM v2.1 (not including channel infiltration). Calibration of the model consists of two parts, including 1) calibration of channel infiltration only with other parameters set to the calibrated parameters used for pre-operational NWM v2.1 and 2) calibration of all parameters including channel infiltration with settings otherwise equivalent to the calibration of NWM v2.1. The calibrated channel-infiltration enhanced NWM improves predictive skill compared to the control NWM in 85% of evaluated basins, for the calibration period. The current NWM settings for physical processes and the biases of the calibration scheme limit model performance in semi-arid environments. To explore whether channel infiltration paired with an alternative calibration scheme could address these limitations, NWM v2.1 was calibrated with a new objective function in selected basins. We found that this updated objective function could ameliorate model biases in some semi-arid environments.

Description: Forecasts of marine cold air outbreaks critically rely on the interplay of multiple parameterisation schemes to represent sub-grid scale processes, including shallow convection, turbulence, and microphysics. Even though such an interplay has been recognised to contribute to forecast uncertainty, a quantification of this interplay is still missing. Here, we investigate the tendencies of temperature and specific humidity contributed by individual parameterisation schemes in the operational weather prediction model AROME-Arctic. From a case study of an extensive marine cold air outbreak over the Nordic Seas, we find that the type of planetary boundary layer assigned by the model algorithm modulates the contribution of individual schemes and affects the interactions between different schemes. In addition, we demonstrate the sensitivity of these interactions to an increase or decrease in the strength of the parameterised shallow convection. The individual tendencies from several parameterisations can thereby compensate each other, sometimes resulting in a small residual. In some instances this residual remains nearly unchanged between the sensitivity experiments, even though some individual tendencies differ by up to an order of magnitude. Using the individual tendency output, we can characterise the subgrid-scale as well as grid-scale responses of the model and trace them back to their underlying causes. We thereby highlight the utility of individual tendency output for understanding process-related differences between model runs with varying physical configurations and for the continued development of numerical weather prediction models.

Description: This study explores the possibilities of employing machine learning algorithms to predict foehn occurrence in Switzerland at a north-Alpine (Altdorf) and south-Alpine (Lugano) station from its synoptic fingerprint in reanalysis data and climate simulations. This allows for an investigation on a potential future shift in monthly foehn frequencies. First, inputs from various atmospheric fields from the European Centre for Medium-Range Weather Forecasts (ECMWF) Re-Analysis-Interim (ERAI) were used to train an XGBoost model. Here, similar predictive performance to previous work was achieved, showing that foehn can accurately be diagnosed from the coarse synoptic situation. In the next step, the algorithm was generalized to predict foehn based on Community Earth System Model (CESM) ensemble simulations of a present-day and warming future climate. The best generalization between ERAI and CESM was obtained by including the present-day data in the training procedure and simultaneously optimizing two objective functions, namely the negative log loss and squared mean loss, on both datasets, respectively. It is demonstrated that the same synoptic fingerprint can be identified in CESM climate simulation data. Finally, predictions for present-day and future simulations were verified and compared for statistical significance. Our model is shown to produce valid output for most months, revealing that south foehn in Altdorf is expected to become more common during spring, while north foehn in Lugano is expected to become more common during summer.

Description: Tropical cyclones are associated with a variety of significant social hazards, including wind, rain, and storm surge. Despite this, most of the model validation effort has been directed toward track and intensity forecasts. In contrast, few studies have investigated the skill of state-of-the-art, high-resolution ensemble prediction systems in predicting associated TC hazards, which is crucial since TC position and intensity do not always correlate with the TC-related hazards, and can result in impacts far from the actual TC center. Furthermore, dynamic models can provide flow-dependent uncertainty estimates, which in turn can provide more specific guidance to forecasters than statistical uncertainty estimates based on past errors. This study validates probabilistic forecasts of wind speed and precipitation hazards derived from the HWRF ensemble prediction system and compares its skill to forecasts by the stochastically-based operational Monte Carlo Model (NHC), the IFS (ECMWF), and the GEFS (NOAA) in use 2017-2019. Wind and Precipitation forecasts are validated against NHC best track wind radii information and the National Stage IV QPE Product. The HWRF 34 kn wind forecasts have comparable skill to the global models up to 60 h lead time before HWRF skill decreases, possibly due to detrimental impacts of large track errors. In contrast, HWRF has comparable quality to its competitors for higher thresholds of 50 kn and 64 kn throughout 120 h lead time. In terms of precipitation hazards, HWRF performs similar or better than global models, but depicts higher, although not perfect, reliability, especially for events over 5 in120h−1. Post-processing, like Quantile Mapping, improves forecast skill for all models significantly and can alleviate reliability issues of the global models.

Description: This study examines historical simulations of ENSO in the E3SM-1-0, CESM2, and GFDL-CM4 climate models, provided by three leading U.S. modeling centers as part of the Coupled Model Intercomparison Project phase 6 (CMIP6). These new models have made substantial progress in simulating ENSO’s key features, including: amplitude; timescale; spatial patterns; phase-locking; spring persistence barrier; and recharge oscillator dynamics. However, some important features of ENSO are still a challenge to simulate. In the central and eastern equatorial Pacific, the models’ weaker-than-observed subsurface zonal current anomalies and zonal temperature gradient anomalies serve to weaken the nonlinear zonal advection of subsurface temperatures, leading to insufficient warm/cold asymmetry of ENSO’s sea surface temperature anomalies (SSTA). In the western equatorial Pacific, the models’ excessive simulated zonal SST gradients amplify their zonal temperature advection, causing their SSTA to extend farther west than observed. The models underestimate both ENSO’s positive dynamic feedbacks (due to insufficient zonal wind stress responses to SSTA) and its thermodynamic damping (due to insufficient convective cloud shading of eastern Pacific SSTA during warm events); compensation between these biases leads to realistic linear growth rates for ENSO, but for somewhat unrealistic reasons. The models also exhibit stronger-than-observed feedbacks onto eastern equatorial Pacific SSTAs from thermocline depth anomalies, which accelerates the transitions between events and shortens the simulated ENSO period relative to observations. Implications for diagnosing and simulating ENSO in climate models are discussed.

Description: In this study, we investigate the response of tropical cyclones (TCs) to climate change by using the Princeton environment-dependent probabilistic tropical cyclone (PepC) model and a statistical-deterministic method to downscale TCs using environmental conditions obtained from the Geophysical Fluid Dynamics Laboratory (GFDL) High-resolution Forecast-oriented Low Ocean Resolution (HiFLOR) model, under the Representative Concentration Pathway 4.5 (RCP4.5) emissions scenario for the North Atlantic basin. The downscaled TCs for the historical climate (1986-2005) are compared with those in the mid- (2016-35) and late-twenty-first century (2081-2100). The downscaled TCs are also compared with TCs explicitly simulated in HiFLOR. We show that while significantly more storms are detected in HiFLOR towards the end of the twenty-first century, the statistical-deterministic model projects a moderate increase in TC frequency, and PepC projects almost no increase in TC frequency. The changes in storm frequency in all three datasets are not significant in the mid-twenty-first century. All three project that storms will become more intense and the fraction of major hurricanes and Category 5 storms will significantly increase in the future climates. However, HiFLOR projects the largest increase in intensity while PepC projects the least. The results indicate that HiFLOR’s TC projection is more sensitive to climate change effects and statistical models are less sensitive. Nevertheless, in all three datasets, storm intensification and frequency increase lead to relatively small changes in TC threat as measured by the return level of landfall intensity.

Description: Tropical convection regimes range from deep organized to shallow convective systems. Mesoscale processes such as cold pools within tropical convective systems can play a significant role in the evolution of convection over land and open ocean. Although cold pools are widely observed, their diurnal properties are not well understood over tropical oceans and land. The oceanic cold pool identification metric applied herein uses the gradient feature (GF) technique and is compared with diurnally-resolved buoy-identified thermal cold pools. This study provides a first-ever diurnal climatology of GF number, area, and attributed TRMM 3B42 precipitation using a space-borne scatterometer (RapidScat). Buoy data over the Pacific, Atlantic, and Indian Ocean have been used to validate and examine the RapidScat-identified diurnal cycle of GF number and precipitation. Buoy-observed cold pool duration, precipitation, temperature, and wind speed is analyzed to understand the in situ cold pool properties over tropical oceans. GF- and buoy-observed cold pool number and precipitation exhibits a similar bimodal diurnal variability with a morning and afternoon maxima, thus establishing confidence in using GF as a proxy to observe cold pools over tropical oceans. The morning peak is attributed to cold pools associated with deep moist convection while the afternoon peak is related to shallower clouds in relatively drier environments resulting in smaller cold pools over global tropical oceans.

Description: A model diagnosis for the energy flux of off-equatorial Rossby waves in the atmosphere has previously been done using quasi-geostrophic equations and is singular at the equator. The energy flux of equatorial waves has been separately investigated in previous studies using a space-time spectral analysis or a ray theory. A recent analytical study has derived an exact universal expression for the energy flux which can indicate the direction of the group velocity for linear shallow water waves at all latitudes. This analytical result is extended in the present study to a height-dependent framework for three-dimensional waves in the atmosphere. This is achieved by investigating the classical analytical solution of both equatorial and off-equatorial waves in a Boussinesq fluid. For the horizontal component of the energy flux, the same expression has been obtained between equatorial waves and off-equatorial waves in the height-dependent framework, which is linked to a scalar quantity inverted from the isentropic perturbation of Ertel’s potential vorticity. The expression of the vertical component of the energy flux requires computation of another scalar quantity that may be obtained from the meridional integral of geopotential anomaly in a wavenumber-frequency space. The exact version of the universal expression is explored and illustrated for three-dimensional waves induced by an idealized Madden-Julian Oscillation forcing in a basic model experiment. The zonal and vertical fluxes manifest the energy transfer of both equatorial Kelvin waves and off-equatorial Rossby waves with a smooth transition at around 10°S and around 10°N. The meridional flux of wave energy represents connection between off-equatorial divergence regions and equatorial convergence regions.

Description: In the hydrological sciences, the outstanding challenge of regional modeling requires to capture common and event-specific hydrologic behaviors driven by rainfall spatial variability and catchment physiography during floods. The overall objective of this study is to develop robust understanding and predictive capability of how rainfall spatial variability influences flood peak discharge relative to basin physiography. A machine learning approach is used on a high-resolution dataset of rainfall and flooding events spanning 10 years, with rainfall events and basins of widely varying characteristics selected across the continental United States. It overcomes major limitations in prior studies that were based on limited observations or hydrological model simulations. This study explores first-order dependencies in the relationships between peak discharge, rainfall variability, and basin physiography, and it sheds light on these complex interactions using a multi-dimensional statistical modeling approach. Amongst different machine learning techniques, XGBoost is used to determine the significant physiographical and rainfall characteristics that influence peak discharge through variable importance analysis. A parsimonious model with low bias and variance is created which can be deployed in the future for flash flood forecasting. The results confirm that although the spatial organization of rainfall within a basin has a major influence on basin response, basin physiography is the primary driver of peak discharge. These findings have unprecedented spatial and temporal representativeness in terms of flood characterization across basins. An improved understanding of sub-basin scale rainfall spatial variability will aid in robust flash flood characterization as well as with identifying basins which could most benefit from distributed hydrologic modeling.

Description: Teleconnections from the Tropics energize variations of the North Pacific climate, but detailed diagnosis of this relationship has proven difficult. Simple univariate methods, such as regression on El Niño-Southern Oscillation (ENSO) indices, may be inadequate since the key dynamical processes involved -- including ENSO diversity in the Tropics, re-emergence of mixed layer thermal anomalies, and oceanic Rossby wave propagation in the North Pacific -- have a variety of overlapping spatial and temporal scales. Here we use a multivariate Linear Inverse Model to quantify tropical and extra-tropical multi-scale dynamical contributions to North Pacific variability, in both observations and CMIP6 models. In observations, we find that the Tropics are responsible for almost half of the seasonal variance, and almost three quarters of the decadal variance, along the North American coast and within the subtropical front region northwest of Hawaii. SST anomalies that are generated by local dynamics within the Northeast Pacific have much shorter time scales, consistent with transient weather forcing by Aleutian low anomalies. Variability within the Kuroshio-Oyashio Extension (KOE) region is considerably less impacted by the Tropics, on all time scales. Consequently, without tropical forcing the dominant pattern of North Pacific variability would be a KOE pattern, rather than the Pacific Decadal Oscillation (PDO). In contrast to observations, most CMIP6 historical simulations produce North Pacific variability that maximizes in the KOE region, with amplitude significantly higher than observed. Correspondingly, the simulated North Pacific in all CMIP6 models is shown to be relatively insensitive to the Tropics, with a dominant spatial pattern generally resembling the KOE pattern, not the PDO.

Description: The extratropical effect of the quasi-biennial oscillation (QBO), known as the Holton-Tan effect, is manifest as aweaker, warmer winter Arctic polar vortex during the east QBO phase. While previous studies have shown that the extratropical QBO signal is caused by the modified propagation of planetary waves in the stratosphere, the mechanism dominating the onset and seasonal development of the Holton-Tan effects remains unclear. Here, the governing wave-mean flow dynamics of the early winter extratropical QBO signal onset and its reversibility is investigated on a synoptic timescale with a finite-amplitude diagnostic using reanalysis and a chemistry-climate model. The extratropical QBO signal onset in October is found to primarily result from modulated stratospheric life-cycles of wave pulses entering the stratosphere from the troposphere, rather than from a modulation of their tropospheric wave source. A comprehensive analysis of the wave activity budget during fall, when the stratospheric winter polar vortex starts forming and waves start propagating up into the stratosphere, shows significant differences. During the east QBO phase, the deceleration of the mid-high latitude stratospheric zonal mean jet by the upward propagating wave pulses is less reversible, due to stronger dissipation processes, while during the west phase, a more reversible deceleration of the main polar vortex is found owing to the waves being dissipated at lower latitudes, accompanied by a weak but different response of the tropospheric subtropical jet. From this synoptic wave-event viewpoint, the early season onset of the Holton-Tan effect results from the cumulative effect of the QBO dependent wave-induced deceleration during the life cycle of individual upward wave pulses.

Description: The provision of climate services has the potential to generate adaptive capacity and help coffee farmers become or remain profitable by integrating climate information in a risk-management framework. Yet, in order to achieve this goal, it is necessary to identify the local demand for climate information, the relationships between coffee yield and climate variables, farmers’ perceptions, and to examine the potential actions that can be realistically put in place by farmers at the local level. In this study, we assessed the climate information demands from coffee farmers and their perception on the climate impacts to coffee yield in the Samalá watershed in Guatemala. After co-identifying the related candidate climate predictors, we propose an objective, flexible forecast system for coffee yield based on precipitation. The system, known as NextGen, analyzes multiple historical climate drivers to identify candidate predictors, and provides both deterministic and probabilistic forecasts for the target season. To illustrate the approach, a NextGen implementation is conducted in the Samalá watershed in southwestern Guatemala. The results suggest that accumulated June-July-August precipitation provides the highest predictive skill associated with coffee yield for this region. In addition to a formal cross-validated skill assessment, retrospective forecasts for the period 1989-2009 were compared to agriculturalists’ perception on the climate impacts to coffee yield at the farm level. We conclude with examples of how demand-based climate service provision in this location can inform adaptation strategies like optimum shade, pest control, and fertilization schemes months in advance. These potential adaptation strategies were validated by local agricultural technicians at the study site.

Description: Atmospheric circulation is a key driver of climate variability, and the representation of atmospheric circulation modes in regional climate models (RCMs) can enhance the credibility of regional climate projections. This study examines the representation of large‐scale atmospheric circulation modes in Coupled Model Inter‐comparison Project phase 5 RCMs once driven by ERA‐Interim, and by two general circulation models (GCMs). The study region is Western Europe and the circulation modes are classified using the Promax rotated T‐mode principal component analysis. The results indicate that the RCMs can replicate the classified atmospheric modes as obtained from ERA5 reanalysis, though with biases dependent on the data providing the lateral boundary condition and the choice of RCM. When the boundary condition is provided by ERA‐Interim that is more consistent with observations, the simulated map types and the associating time series match well with their counterparts from ERA5. Further, on average, the multi‐model ensemble mean of the analysed RCMs, driven by ERA‐Interim, indicated a slight improvement in the representation of the modes obtained from ERA5. Conversely, when the RCMs are driven by the GCMs that are models without assimilation of observational data, the representation of the atmospheric modes, as obtained from ERA5, is relatively less accurate compared to when the RCMs are driven by ERA‐Interim. This suggests that the biases stem from the GCMs. On average, the representation of the modes was not improved in the multi‐model ensemble mean of the five analysed RCMs driven by either of the GCMs. However, when the best‐performed RCMs were selected on average the ensemble mean indicated a slight improvement. Moreover, the presence of the North Atlantic Oscillation (NAO) in the simulated modes depends also on the lateral boundary conditions. The relationship between the modes and the NAO was replicated only when the RCMs were driven by reanalysis. The results indicate that the forcing model is the main factor in reproducing the atmospheric circulation.

Description: Regional and local wind systems are often complex, particularly near coastal areas with a highly variable orography. Thus, the realistic representation of regional wind systems in weather and climate models is of strong relevance. Here, we evaluate the ability of a 13‐year convection‐permitting climate simulation in reproducing the interaction of several regional summer wind systems over the complex orography in the eastern Mediterranean region. The COSMO‐CLM simulations are driven by hourly ERA‐5 reanalysis and have a spatial resolution of 2.8 and 7.0 km. The simulated near‐surface wind fields are compared with unique very high‐resolution wind observations collected within the “Dead Sea Research Venue” project (DESERVE) and data from the Israel Meteorological Service synop network. The high‐resolution COSMO‐CLM simulations largely reproduce the main characteristics of the regional wind systems (Mediterranean and Dead Sea breeze, slope winds in the Judean Mountains and winds along the Jordan Rift valley), whereas ERA‐5 is only able to represent the Mediterranean Sea breeze. The high‐resolution simulations substantially improve the representation of regional winds, particularly over complex orography. Indeed, the 2.8 km simulation outperforms the 7.0 km run, on 88% of the days. Two mid‐July 2015 case studies show that only the 2.8 simulation can realistically simulate the penetration of the Mediterranean Sea Breeze into the Jordan Rift valley and complex interactions with other wind systems like the Dead Sea breeze. Our results may have profound implications for regional weather and climate prediction since very high‐resolution information seems to be necessary to reproduce the main summertime climatic features in this region. We envisage that such simulations may also be required at other regions with complex orography.

Description: In this paper we show that COSMO‐CLM regional climate model simulations at 7.0 (CLM‐7.0) and 2.8km (CLM‐2.8) resolution can realistically reproduce near‐surface regional and local wind systems over the complex orography of the eastern Mediterranean as opposite to coarser resolutions (ERA‐5, 31 km). The Mediterranean and local Dead Sea breezes, slope winds over the Judean Mountains, and winds along the Jordan Rift valley are well represented both climatologically and on individual days. CLM‐2.8 captures the small‐scale variability of the wind field better than CLM‐7.0 particularly near the Dead Sea and on 88% of the days CLM‐2.8 represents wind speed even more realistically than CLM‐7.0. image

Description: German Helmholtz Association (“Changing Earth” program)

Description: AXA Research Fund http://dx.doi.org/10.13039/501100001961

Description: Ministry of Science, Research and Arts

Description: Helmholtz Association of German Research Centers

Description: Author Posting. © American Meteorological Society, 2022. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 52(12),(2022): 3199-3219, https://doi.org/10.1175/jpo-d-22-0009.1.

Description: The abyssal overturning circulation is thought to be primarily driven by small-scale turbulent mixing. Diagnosed water-mass transformations are dominated by rough topography “hotspots,” where the bottom enhancement of mixing causes the diffusive buoyancy flux to diverge, driving widespread downwelling in the interior—only to be overwhelmed by an even stronger upwelling in a thin bottom boundary layer (BBL). These water-mass transformations are significantly underestimated by one-dimensional (1D) sloping boundary layer solutions, suggesting the importance of three-dimensional physics. Here, we use a hierarchy of models to generalize this 1D boundary layer approach to three-dimensional eddying flows over realistically rough topography. When applied to the Mid-Atlantic Ridge in the Brazil Basin, the idealized simulation results are roughly consistent with available observations. Integral buoyancy budgets isolate the physical processes that contribute to realistically strong BBL upwelling. The downward diffusion of buoyancy is primarily balanced by upwelling along the sloping canyon sidewalls and the surrounding abyssal hills. These flows are strengthened by the restratifying effects of submesoscale baroclinic eddies and by the blocking of along-ridge thermal wind within the canyon. Major topographic sills block along-thalweg flows from restratifying the canyon trough, resulting in the continual erosion of the trough’s stratification. We propose simple modifications to the 1D boundary layer model that approximate each of these three-dimensional effects. These results provide local dynamical insights into mixing-driven abyssal overturning, but a complete theory will also require the nonlocal coupling to the basin-scale circulation.

Description: We acknowledge funding support from National Science Foundation Awards 1536515, 1736109, and 2149080. This material is based upon work supported by the National Science Foundation Graduate Research Fellowship Program under Grant 174530.

Description: 2023-05-18

Description: Author Posting. © American Meteorological Society, 2022. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 52(6), (2022): 1091–1110, https://doi.org/10.1175/JPO-D-21-0068.1.

Description: Hundreds of full-depth temperature and salinity profiles collected by Deepglider autonomous underwater vehicles (AUVs) in the North Atlantic reveal robust signals in eddy isopycnal vertical displacement and horizontal current throughout the entire water column. In separate glider missions southeast of Bermuda, subsurface-intensified cold, fresh coherent vortices were observed with velocities exceeding 20 cm s−1 at depths greater than 1000 m. With vertical resolution on the order of 20 m or less, these full-depth glider slant profiles newly permit estimation of scaled vertical wavenumber spectra from the barotropic through the 40th baroclinic mode. Geostrophic turbulence theory predictions of spectral slopes associated with the forward enstrophy cascade and proportional to inverse wavenumber cubed generally agree with glider-derived quasi-universal spectra of potential and kinetic energy found at a variety of locations distinguished by a wide range of mean surface eddy kinetic energy. Water-column average spectral estimates merge at high vertical mode number to established descriptions of internal wave spectra. Among glider mission sites, geographic and seasonal variability implicate bottom drag as a mechanism for dissipation, but also the need for more persistent sampling of the deep ocean.

Description: This work was funded by NSF Grant 1736217 and would not have been possible without the help of Kirk O’Donnell, James Bennett, Noel Pelland, and all contributors to Deepglider development. We additionally thank the captain crew of the R/V Atlantic Explorer and the BATS team at the Bermuda Institute of Ocean Sciences, particularly Rod Johnson, as well as Seakeepers International for their professionalism, capability, and generous assistance in deploying and recovering gliders.

Description: Climate model simulations typically exhibit a bias, which can be corrected using statistical approaches. In this study, a geostatistical approach for bias correction of daily precipitation at ungauged locations is presented. The method utilizes a double quantile mapping with dry day correction for future periods. The transfer function of the bias correction for the ungauged locations is established using distribution functions estimated by ordinary kriging with anisotropic variograms. The methodology was applied to the daily precipitation simulations of the entire CORDEX‐Africa ensemble for a study region located in the West African Sudanian Savanna. This ensemble consists of 23 regional climate models (RCM) that were run for three different future scenarios (RCP 2.6, RCP 4.5, and RCP 8.5). The evaluation of the approach for a historical 50‐year period (1950–2005) showed that the method can reduce the inherent strong precipitation bias of RCM simulations, thereby reproducing the main climatological features of the observed data. Moreover, the bias correction technique preserves the climate change signal of the uncorrected RCM simulations. However, the ensemble spread is increased due to an overestimation of the rainfall probability of uncorrected RCM simulations. The application of the bias correction method to the future period (2006–2100) revealed that annual precipitation increases for most models in the near (2020–2049) and far future (2070–2099) with a mean increase of up to 165mm⋅a−1 (18%). An analysis of the monthly and daily time series showed a slightly delayed onset and intensification of the rainy season.

Description: Adapting water management strategies to future precipitation projected by climate models is associated with high uncertainty in sparsely gauged catchments. Kriging was utilized to estimate distribution parameters for ungauged locations in a West African region to perform a bias correction of the CORDEX‐Africa ensemble. The application of the bias correction method revealed higher annual precipitation amounts and an intensifaction of the rainy season but only little change to the onset of the rainy season.

Description: German Federal Ministry of Education and Research, Bonn (BMBF), West African Science Service Centre on Climate Change and Adapted Land Use (WASCAL)

Description: Invasive alien species continue to spread and proliferate in waterways worldwide, but environmental drivers of invasion dynamics lack assessment. Knowledge gaps are pervasive in the Global South, where the frequent heavy human‐modification of rivers provides high opportunity for invasion. In southern Africa, the spatio‐temporal ecology of a widespread and high‐impact invasive alien snail, Tarebia granifera, and its management status is understudied. Here, an ecological assessment was conducted at seven sites around Nandoni Reservoir on the Luvuvhu River in South Africa. The distribution and densities of T. granifera were mapped and the potential drivers of population structure were explored. T. granifera was widespread at sites impacted to varying extents due to anthropogenic activity, with densities exceeding 500 individuals per square meter at the most impacted areas. T. granifera predominantly preferred shallow and sandy environments, being significantly associated with sediment (i.e., chlorophyll‐a, Mn, SOC, SOM) and water (i.e., pH, conductivity, TDS) variables. T. granifera seemed to exhibit two recruitment peaks in November and March, identified via size‐based stock assessment. Sediment parameters (i.e., sediment organic matter, sediment organic carbon, manganese) and water chemistry (i.e., pH, total dissolved solids, conductivity) were found to be important in structuring T. granifera populations, with overall snail densities highest during the summer season. We provide important autecological information and insights on the distribution and extent of the spread of T. granifera. This may help in the development of invasive alien snail management action plans within the region, as well as modelling efforts to predict invasion patterns elsewhere based on environmental characteristics.

Description: Alexander von Humboldt‐Stiftung http://dx.doi.org/10.13039/100005156

Description: National Research Foundation http://dx.doi.org/10.13039/501100001321

Description: University of Venda http://dx.doi.org/10.13039/501100008976

Description: Author Posting. © American Meteorological Society, 2022. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 52(8), (2022): 1705-1730, https://doi.org/10.1175/jpo-d-21-0243.1.

Description: Formation and evolution of barrier layers (BLs) and associated temperature inversions (TIs) were investigated using a 1-yr time series of oceanic and air–sea surface observations from three moorings deployed in the eastern Pacific fresh pool. BL thickness and TI amplitude showed a seasonality with maxima in boreal summer and autumn when BLs were persistently present. Mixed layer salinity (MLS) and mixed layer temperature (MLT) budgets were constructed to investigate the formation mechanism of BLs and TIs. The MLS budget showed that BLs were initially formed in response to horizontal advection of freshwater in boreal summer and then primarily maintained by precipitation. The MLT budget revealed that penetration of shortwave radiation through the mixed layer base is the dominant contributor to TI formation through subsurface warming. Geostrophic advection is a secondary contributor to TI formation through surface cooling. When the BL exists, the cooling effect from entrainment and the warming effect from detrainment are both significantly reduced. In addition, when the BL is associated with the presence of a TI, entrainment works to warm the mixed layer. The presence of BLs makes the shallower mixed layer more sensitive to surface heat and freshwater fluxes, acting to enhance the formation of TIs that increase the subsurface warming via shortwave penetration.

Description: SK is supported by JSPS Overseas Research Fellowships. JS and SK are supported by NASA Grant 80NSSC18K1500. JTF and the mooring deployment were funded by NASA Grants NNX15AG20G and 80NSSC18K1494. DZ is supported by NASA Grant 80NSSC18K1499. This publication is partially funded by the Cooperative Institute for Climate, Ocean, and Ecosystem Studies (CICOES) under NOAA Cooperative Agreement NA20OAR4320271, Contribution 2021-1152. This is PMEL Contribution 5268.

Description: 2023-01-27

Description: 〈jats:title〉Abstract〈/jats:title〉 〈jats:p〉Tipping points in the Earth system describe critical thresholds beyond which a single component, part of the system, or the system as a whole changes from one stable state to another. In the present-day Southern Ocean, the Weddell Sea constitutes an important dense-water formation site, associated with efficient deep-ocean carbon and oxygen transfer and low ice-shelf basal melt rates. Here, a regime shift will occur when continental shelves are continuously flushed with warm, oxygen-poor offshore waters from intermediate depth, leading to less efficient deep-ocean carbon and oxygen transfer and higher ice-shelf basal melt rates. We use a global ocean–biogeochemistry model including ice-shelf cavities and an eddy-permitting grid in the southern Weddell Sea to address the susceptibility of this region to such a system change for four 21〈jats:sup〉st〈/jats:sup〉-century emission scenarios. Assessing the projected changes in shelf–open ocean density gradients, bottom-water properties, and on-shelf heat transport, our results indicate that the Weddell Sea undergoes a regime shift by 2100 in the highest-emission scenario SSP5-8.5, but not yet in the lower-emission scenarios. The regime shift is imminent by 2100 in the scenarios SSP3-7.0 and SSP2-4.5, but avoidable under the lowest-emission scenario SSP1-2.6. While shelf-bottom waters freshen and acidify everywhere, bottom waters in the Filchner Trough undergo accelerated warming and deoxygenation following the system change, with implications for local ecosystems and ice-shelf basal melt. Additionally, deep-ocean carbon and oxygen transfer decline, implying that the local changes ultimately affect ocean circulation, climate, and ecosystems globally.〈/jats:p〉

Description: 〈title xmlns:mml="http://www.w3.org/1998/Math/MathML"〉Abstract〈/title〉〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉The Lagrangian representation of fluid flows offers a natural perspective to study many kinds of physical mechanisms. By contrast, the Eulerian representation is more convenient from a diagnostic point of view. This article attempts to combine elements of both worlds by proposing an Eulerian method that allows one to extract Lagrangian information about the atmospheric flow. The method is based on the offline advection of passive tracer fields and includes a relaxation term. The latter device allows one to run the integration in a continuous fashion without the need for reinitialization. As a result one obtains accumulated Lagrangian information, for example, about the recent parcel displacement or the recent parcel‐based diabatic heating, at each point of an Eulerian grid at any time step. The method is implemented with a pseudospectral algorithm suitable for gridded global atmospheric data and compared with the more traditional trajectory method. The method's utility is demonstrated on the basis of a few examples, which relate to cloud formation and the development of temperature anomalies. The examples highlight that the method provides a convenient diagnostic of parcel‐based changes, paving an intuitive way to explore the physical processes involved. Due to its gridpoint‐based nature, the proposed method can be applied to large data sets in a straightforward and computationally efficient manner, suggesting that the method is particularly useful for climatological analyses.〈/p〉

Description: 〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉The Lagrangian representation of fluid flows offers the most natural perspective to study many kinds of physical mechanisms; by contrast, the Eulerian representation is more convenient from a diagnostic point of view. This article attempts to combine elements of both worlds by proposing an Eulerian method that allows one to extract Lagrangian information about the atmospheric flow. The method enables one to easily produce a sequence of maps showing accumulated Lagrangian changes. 〈boxed-text position="anchor" id="qj4453-blkfxd-0001" content-type="graphic" xml:lang="en"〉〈graphic position="anchor" id="jats-graphic-1" xlink:href="urn:x-wiley:00359009:media:qj4453:qj4453-toc-0001"〉 〈/boxed-text〉〈/p〉

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: Numerical simulations allow us to gain a comprehensive understanding of the underlying mechanisms of past, present, and future climate changes. The mid-Holocene (MH) and the last interglacial (LIG) were the two most recent warm episodes of Earth’s climate history and are the focus of paleoclimate research. Here, we present results of MH and LIG simulations with two versions of the state-of-the-art Earth system model AWI-ESM. Most of the climate changes in MH and LIG compared to the preindustrial era are agreed upon by the two model versions, including 1) enhanced seasonality in surface temperature that is driven by the redistribution of seasonal insolation; 2) a northward shift of the intertropical convergence zone (ITCZ) and tropical rain belt; 3) a reduction in annual mean Arctic sea ice concentration; 4) weakening and northward displacement of the Northern Hemisphere Hadley circulation, which is related to the decrease and poleward shift of the temperature gradient from the subtropical to the equator in the Northern Hemisphere; 5) a westward shift of the Indo-Pacific Walker circulation due to anomalous warming over the Eurasia and North Africa during boreal summer; and 6) an expansion and intensification of Northern Hemisphere summer monsoon rainfall, with the latter being dominated by the dynamic component of moisture budget (i.e., the strengthening of wind circulation). However, the simulated responses of the Atlantic meridional overturning circulation (AMOC) in the two models yield different results for both the LIG and the MH. AMOC anomalies between the warm interglacial and preindustrial periods are associated with changes in North Atlantic westerly winds and stratification of the water column at the North Atlantic due to changes in ocean temperature, salinity, and density.

Description: 〈title xmlns:mml="http://www.w3.org/1998/Math/MathML"〉Abstract〈/title〉〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉Regional assessments of the wind erosion risk are rare and vary due to the methods used and the available data to be included. The adaptation of existing methods has the advantage that the results can be compared directly. We adopted an already successfully applied methodology (ILSWE—applied in East Africa), to investigate the spatiotemporal variability of the wind erosion risk between 2005 and 2019 in Southern Africa. The approach integrates climatic variables, a vegetation index, and soil properties to describe the potential impact of wind erosion at the landscape scale. The annual and seasonal variability is determined by the vegetation cover, whereas droughts and strong El Niño events had only regional effects. We estimated that 8.3% of the study area experiences a moderate to elevated wind erosion risk over the 15‐year period with annual and inter‐annual fluctuations showing a slight upward trend. In general, the desert and drylands in the west have the highest proportion of risk areas, the moist forests in the east are characterized by a very low risk of wind erosion, while the grasslands, shrublands, and croplands in the interior most likely react to changes of climatic conditions. The validation process is based on a comparison with the estimated frequency of dust storms derived from the aerosol optical depth and angstrom exponent and revealed an overall accuracy of 65%. The results of this study identify regions and yearly periods prone to wind erosion to prioritize for further analysis and conservation policies for mitigation and adaptation strategies.〈/p〉

Description: Bundesministerium für Bildung und Forschung http://dx.doi.org/10.13039/501100002347

Description: German Federal Ministry of Education and Research (BMBF)

Description: http://www.climatologylab.org/terraclimate

Description: https://www.isric.org/

Description: http://www.fao.org/soils-portal/data-hub/soil-maps-and-databases/harmonized-world-soil-database-v12

Description: https://land.copernicus.eu/global/products/fcover

Description: https://www.esa-landcover-cci.org/%20

Description: https://databasin.org/

Description: https://giovanni.gsfc.nasa.gov/giovanni/

Description: 〈title xmlns:mml="http://www.w3.org/1998/Math/MathML"〉Abstract〈/title〉〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉Many operational weather services use ensembles of forecasts to generate probabilistic predictions. Computational costs generally limit the size of the ensemble to fewer than 100 members, although the large number of degrees of freedom in the forecast model would suggest that a vastly larger ensemble would be required to represent the forecast probability distribution accurately. In this study, we use a computationally efficient idealised model that replicates key properties of the dynamics and statistics of cumulus convection to identify how the sampling uncertainty of statistical quantities converges with ensemble size. Convergence is quantified by computing the width of the 95% confidence interval of the sampling distribution of random variables, using bootstrapping on the ensemble distributions at individual time and grid points. Using ensemble sizes of up to 100,000 members, it was found that for all computed distribution properties, including mean, variance, skew, kurtosis, and several quantiles, the sampling uncertainty scaled as 〈mml:math id="jats-math-1" display="inline" overflow="scroll"〉〈mml:msup〉〈mml:mrow〉〈mml:mi〉n〈/mml:mi〉〈/mml:mrow〉〈mml:mrow〉〈mml:mo form="prefix"〉−〈/mml:mo〉〈mml:mn〉1〈/mml:mn〉〈mml:mo stretchy="false"〉/〈/mml:mo〉〈mml:mn〉2〈/mml:mn〉〈/mml:mrow〉〈/mml:msup〉〈/mml:math〉 for sufficiently large ensemble size 〈mml:math id="jats-math-2" display="inline" overflow="scroll"〉〈mml:mrow〉〈mml:mi〉n〈/mml:mi〉〈/mml:mrow〉〈/mml:math〉. This behaviour is expected from the Central Limit Theorem, which further predicts that the magnitude of the uncertainty depends on the distribution shape, with a large uncertainty for statistics that depend on rare events. This prediction was also confirmed, with the additional observation that such statistics also required larger ensemble sizes before entering the asymptotic regime. By considering two methods for evaluating asymptotic behaviour in small ensembles, we show that the large‐〈mml:math id="jats-math-3" display="inline" overflow="scroll"〉〈mml:mrow〉〈mml:mi〉n〈/mml:mi〉〈/mml:mrow〉〈/mml:math〉 theory can be applied usefully for some forecast quantities even for the ensemble sizes in operational use today.〈/p〉

Description: 〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉An idealised ensemble that replicates key properties of the dynamics and statistics of cumulus convection is used to identify how sampling uncertainty of statistical quantities converges with ensemble size. A universal asymptotic scaling for this convergence was found, which was dependent on the statistic and the distribution shape, with largest uncertainty for statistics that depend on rare events. This is demonstrated in the figure below for a Gaussian distributed model variable, where the sampling uncertainty (y‐axis) for 5 quantiles (red lines) indicates that after a certain ensemble size, it begins converging asymptotically (grey lines), and the more extreme the quantile, the more members it requires for this to be the case. 〈boxed-text position="anchor" id="qj4410-blkfxd-0001" content-type="graphic" xml:lang="en"〉〈graphic position="anchor" id="jats-graphic-1" xlink:href="urn:x-wiley:00359009:media:qj4410:qj4410-toc-0001"〉

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: Klaus Tschira Stiftung http://dx.doi.org/10.13039/501100007316

Description: 〈title xmlns:mml="http://www.w3.org/1998/Math/MathML"〉Abstract〈/title〉〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉Local ensemble transform Kalman filters (LETKFs) allow explicit calculation of the Kalman gain, and by this the contribution of individual observations to the analysis field. Though this is a known feature, the information on the analysis contribution of individual observations (partial analysis increment) has not been used as systematic diagnostic up to now despite providing valuable information. In this study, we demonstrate three potential applications based on partial analysis increments in the regional modelling system of Deutscher Wetterdienst and propose their use for optimising LETKF data assimilation systems, in particular with respect to satellite data assimilation and localisation. While exact calculation of partial analysis increments would require saving the large, five‐dimensional ensemble weight matrix in the analysis step, it is possible to compute an approximation from standard LETKF output. We calculate the Kalman gain based on ensemble analysis perturbations, which is an approximation in the case of localisation. However, this only introduces minor errors, as the localisation function changes very gradually among nearby grid points. On the other hand, the influence of observations always depends on the presence of other observations and settings for the observation error and for localisation. However, the influence of observations behaves approximately linearly, meaning that the assimilation of other observations primarily decreases the magnitude of the influence, but it does not change the overall structure of the partial analysis increments. This means that the calculation of partial analysis increments can be used as an efficient diagnostic to investigate the three‐dimensional influence of observations in the assimilation system. Furthermore, the diagnostic can be used to detect whether the influence of additional experimental observations is in accordance with other observations without conducting computationally expensive single‐observation experiments. Last but not least, the calculation can be used to approximate the influence an observation would have when applying different assimilation settings.〈/p〉

Description: Long believed to be insignificant, melt activity on the Northeast Greenland Ice Stream (NEGIS) has increased in recent years. Summertime Arctic clouds have the potential to strongly affect surface melt processes by regulating the amount of radiation received at the surface. However, the cloud effect over Greenland is spatially and temporally variable and high‐resolution information on the northeast is absent. This study aims at exploring the potential of a high‐resolution configuration of the polar‐optimized Weather Research & Forecasting Model (PWRF) in simulating cloud properties in the area of the Nioghalvfjerdsfjorden Glacier (79 N Glacier). Subsequently, the model simulations are employed to investigate the impact of Arctic clouds on the surface energy budget and on surface melting during the extensive melt event at the end of July 2019. Compared to automatic weather station (AWS) measurements and remote‐sensing data (Sentinel‐2A and the Moderate Resolution Imaging Spectroradiometer, MODIS), PWRF simulates cloud properties with sufficient accuracy. It appears that peak melt was caused by an increase in solar radiation and sensible heat flux (SHF) in response to a blocking anticyclone and foehn winds in the absence of clouds. Cloud warming over high‐albedo surfaces helped to precondition the surface and prolonged the melting as the anticyclone abated. The results are sensitive to the surface albedo and suggest spatiotemporal differences in the cloud effect as snow and ice properties change over the course of the melting season. This demonstrates the importance of including high‐resolution information on clouds in analyses of ice sheet dynamics.

Description: German Federal Ministry for Education and Research http://dx.doi.org/10.13039/501100002347

Description: https://doi.org/10.5065/EM0T-1D34

Description: https://cds.climate.copernicus.eu/cdsapp/#!/search?type=dataset

Description: https://ladsweb.modaps.eosdis.nasa.gov/search/

Description: 〈title xmlns:mml="http://www.w3.org/1998/Math/MathML"〉Abstract〈/title〉〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉For both the meso‐ and synoptic scales, reduced mathematical models give insight into their dynamical behaviour. For the mesoscale, the weak temperature gradient approximation is one of several approaches, while for the synoptic scale the quasigeostrophic theory is well established. However, the way these two scales interact with each other is usually not included in such reduced models, thereby limiting our current perception of flow‐dependent predictability and upscale error growth. Here, we address the scale interactions explicitly by developing a two‐scale asymptotic model for the meso‐ and synoptic scales with two coupled sets of equations for the meso‐ and synoptic scales respectively. The mesoscale equations follow a weak temperature gradient balance and the synoptic‐scale equations align with quasigeostrophic theory. Importantly, the equation sets are coupled via scale‐interaction terms: eddy correlations of mesoscale variables impact the synoptic potential vorticity tendency and synoptic variables force the mesoscale vorticity (for instance due to tilting of synoptic‐scale wind shear). Furthermore, different diabatic heating rates—representing the effect of precipitation—define different flow characteristics. With weak mesoscale heating relatable to precipitation rates of 〈mml:math id="jats-math-1" display="inline" overflow="scroll"〉〈mml:mrow〉〈mml:mi〉𝒪〈/mml:mi〉〈mml:mo stretchy="false"〉(〈/mml:mo〉〈mml:mn〉6〈/mml:mn〉〈mml:mspace width="0.3em"/〉〈mml:mtext〉mm〈/mml:mtext〉〈mml:mo〉·〈/mml:mo〉〈mml:msup〉〈mml:mrow〉〈mml:mi mathvariant="normal"〉h〈/mml:mi〉〈/mml:mrow〉〈mml:mrow〉〈mml:mo form="prefix"〉−〈/mml:mo〉〈mml:mn〉1〈/mml:mn〉〈/mml:mrow〉〈/mml:msup〉〈mml:mo stretchy="false"〉)〈/mml:mo〉〈/mml:mrow〉〈/mml:math〉, the mesoscale dynamics resembles two‐dimensional incompressible vorticity dynamics and the upscale impact of the mesoscale on the synoptic scale is only of a dynamical nature. With a strong mesosocale heating relatable to precipitation rates of 〈mml:math id="jats-math-2" display="inline" overflow="scroll"〉〈mml:mrow〉〈mml:mi〉𝒪〈/mml:mi〉〈mml:mo stretchy="false"〉(〈/mml:mo〉〈mml:mn〉60〈/mml:mn〉〈mml:mspace width="0.3em"/〉〈mml:mtext〉mm〈/mml:mtext〉〈mml:mo〉·〈/mml:mo〉〈mml:msup〉〈mml:mrow〉〈mml:mi mathvariant="normal"〉h〈/mml:mi〉〈/mml:mrow〉〈mml:mrow〉〈mml:mo form="prefix"〉−〈/mml:mo〉〈mml:mn〉1〈/mml:mn〉〈/mml:mrow〉〈/mml:msup〉〈mml:mo stretchy="false"〉)〈/mml:mo〉〈/mml:mrow〉〈/mml:math〉, divergent motions and three‐dimensional effects become relevant for the mesoscale dynamics and the upscale impact also includes thermodynamical effects.〈/p〉

Description: 〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉We develop a two‐scale asymptotic model for the meso‐ and synoptic scales following a weak temperature gradient balance and quasigeostrophic theory, but with explicit scale interactions and dependent on the mesoscale diabatic heating. With weak mesoscale heating, the mesoscale dynamics resembles 2D incompressible vorticity dynamics and the upscale impact on the synoptic scale is only of a dynamical nature. With strong mesoscale heating, divergent motions and 3D effects become relevant for the mesoscale and the upscale impact also includes thermodynamical effects. 〈boxed-text position="anchor" id="qj4456-blkfxd-0001" content-type="graphic" xml:lang="en"〉〈graphic position="anchor" id="jats-graphic-1" xlink:href="urn:x-wiley:00359009:media:qj4456:qj4456-toc-0001"〉

Description: German Research Foundation (DFG)

Description: Author Posting. © American Meteorological Society, 2022. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 52(12), (2022): 3221–3240, https://doi.org/10.1175/jpo-d-22-0010.1.

Description: Small-scale mixing drives the diabatic upwelling that closes the abyssal ocean overturning circulation. Indirect microstructure measurements of in situ turbulence suggest that mixing is bottom enhanced over rough topography, implying downwelling in the interior and stronger upwelling in a sloping bottom boundary layer. Tracer release experiments (TREs), in which inert tracers are purposefully released and their dispersion is surveyed over time, have been used to independently infer turbulent diffusivities—but typically provide estimates in excess of microstructure ones. In an attempt to reconcile these differences, Ruan and Ferrari derived exact tracer-weighted buoyancy moment diagnostics, which we here apply to quasi-realistic simulations. A tracer’s diapycnal displacement rate is exactly twice the tracer-averaged buoyancy velocity, itself a convolution of an asymmetric upwelling/downwelling dipole. The tracer’s diapycnal spreading rate, however, involves both the expected positive contribution from the tracer-averaged in situ diffusion as well as an additional nonlinear diapycnal distortion term, which is caused by correlations between buoyancy and the buoyancy velocity, and can be of either sign. Distortion is generally positive (stretching) due to bottom-enhanced mixing in the stratified interior but negative (contraction) near the bottom. Our simulations suggest that these two effects coincidentally cancel for the Brazil Basin Tracer Release Experiment, resulting in negligible net distortion. By contrast, near-bottom tracers experience leading-order distortion that varies in time. Errors in tracer moments due to realistically sparse sampling are generally small (〈20%), especially compared to the O(1) structural errors due to the omission of distortion effects in inverse models. These results suggest that TREs, although indispensable, should not be treated as “unambiguous” constraints on diapycnal mixing.

Description: We acknowledge funding support from National Science Foundation Awards 1536515 and 1736109. This material is based upon work supported by the National Science Foundation Graduate Research Fellowship Program under Grant 174530. This research is also supported by the NOAA Climate and Global Change Postdoctoral Fellowship Program, administered by UCAR’s Cooperative Programs for the Advancement of Earth System Science (CPAESS) under Award NA18NWS4620043B.

Description: 2023-05-18

Description: Author Posting. © American Meteorological Society, 2022. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 52(11), (2022): 2841–2852, https://doi.org/10.1175/jpo-d-22-0025.1.

Description: Prediction of rapid intensification in tropical cyclones prior to landfall is a major societal issue. While air–sea interactions are clearly linked to storm intensity, the connections between the underlying thermal conditions over continental shelves and rapid intensification are limited. Here, an exceptional set of in situ and satellite data are used to identify spatial heterogeneity in sea surface temperatures across the inner core of Hurricane Sally (2020), a storm that rapidly intensified over the shelf. A leftward shift in the region of maximum cooling was observed as the hurricane transited from the open gulf to the shelf. This shift was generated, in part, by the surface heat flux in conjunction with the along- and across-shelf transport of heat from storm-generated coastal circulation. The spatial differences in the sea surface temperatures were large enough to potentially influence rapid intensification processes suggesting that coastal thermal features need to be accounted for to improve storm forecasting as well as to better understand how climate change will modify interactions between tropical cyclones and the coastal ocean.

Description: This research was made possible by the NOAA RESTORE Science Program (NA17NOS4510101 and NA19NOS4510194) and the NASA Physical Oceanography program (80NSSC21K0553 and WBS 281945.02.25.04.67) and NOAA IOOS program via GCOOS (NA16NOS0120018). The authors declare that they have no competing interests.

Description: Author Posting. © American Meteorological Society, 2022. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 52(8), (2022): 1797–1815, https://doi.org/10.1175/JPO-D-21-0288.1.

Description: Intruding slope water is a major source of nutrients to sustain the high biological productivity in the Gulf of Maine (GoM). Slope water intrusion into the GoM is affected by Gulf Stream warm-core rings (WCRs) impinging onto the nearby shelf edge. This study combines long-term mooring measurements, satellite remote sensing data, an idealized numerical ocean model, and a linear coastal-trapped wave (CTW) model to examine the impact of WCRs on slope water intrusion into the GoM through the Northeast Channel. Analysis of satellite sea surface height and temperature data shows that the slope sea region off the GoM is a hotspot of ring activities. A significant linear relationship is found between interannual variations of ring activities in the slope sea region off the GoM and bottom salinity at the Northeast Channel, suggesting the importance of WCRs in modulating variability of intruding slope water. Analysis of the mooring data reveals enhanced slope water intrusion through bottom-intensified along-channel flow following impingements of WCRs on the nearby shelf edge. Numerical simulations qualitatively reproduce the observed WCR impingement processes and associated episodic enhancement of slope water intrusion in the Northeast Channel. Diagnosis of the model result indicates that baroclinic CTWs excited by the ring–topography interaction are responsible for the episodically intensified subsurface along-channel inflow, which carries more slope water into the GoM. A WCR that impinges onto the shelf edge to the northeast of the Northeast Channel tends to generate stronger CTWs and cause stronger enhancement of the slope water intrusion into the GoM.

Description: This study is supported by the National Science Foundation through Grant OCE-1634965.

Description: Author Posting. © American Meteorological Society, 2022. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of the Atmospheric and Oceanic Technology 39(10), (2022): 1525–1539, https://doi.org/10.1175/jtech-d-21-0186.1.

Description: The static and dynamic performances of the RBRargo3 are investigated using a combination of laboratory-based and in situ datasets from floats deployed as part of an Argo pilot program. Temperature and pressure measurements compare well to co-located reference data acquired from shipboard CTDs. Static accuracy of salinity measurements is significantly improved using 1) a time lag for temperature, 2) a quadratic pressure dependence, and 3) a unit-based calibration for each RBRargo3 over its full pressure range. Long-term deployments show no significant drift in the RBRargo3 accuracy. The dynamic response of the RBRargo3 demonstrates the presence of two different adjustment time scales: a long-term adjustment O(120) s, driven by the temperature difference between the interior of the conductivity cell and the water, and a short-term adjustment O(5–10) s, associated to the initial exchange of heat between the water and the inner ceramic. Corrections for these effects, including dependence on profiling speed, are developed.

Description: During the last 20 years some very hot and dry summers affected Europe, resulting in regionally record‐breaking high temperature or low precipitation values. Long‐term changes of such extremely hot and dry summers are of great relevance for our society, as they are connected with manifold negative impacts on human society, natural ecosystems, and diverse economic sectors. Long‐term variations in drought and five record drought summer half years are studied based on 63 stations across Europe with high‐quality precipitation and temperature time series spanning the period 1901–2018. Eight drought indices are deployed to analyse drought intensity, frequency, and duration; four of them purely precipitation‐based and four integrating potential evapotranspiration in the computation. Additionally, three heavy precipitation indices and simultaneous increases in drought and heavy precipitation are studied. The five driest summer half years over Europe are identified (1947, 2018, 2003, 1921, and 1911). They are analysed by aggregating eight drought indices into the aggregated drought evaluation index (ADE) for five subregions. The ADE shows increasing summer drought conditions over most of Europe, except for some stations in northern Europe. The increase in drought conditions during the warm part of the year is particularly pronounced for indices integrating evapotranspiration in their definition. At the same time, the intensity of heavy precipitation events shows a positive trend, as well as an increased contribution to total precipitation. Several stations in central Europe show simultaneously increasing drought conditions and increasing heavy precipitation events. This increases the risks connected with precipitation extremes.

Description: Drought and heavy precipitation trends are studied for selected stations with long time series. Both extreme events are related to specific impacts on different economic sectors and thus society. Stations that simultaneously show increasing trends in drought and heavy precipitation mainly occur in central Europe. This indicates a probably higher exposure to these risks in central Europe as well as a demand for broader adaptation options in this region.

Description: Europe has been affected by record‐breaking heat waves in recent decades. Using station data and a gridded reanalysis as input, four commonly used heat wave indices, the heat wave magnitude index daily (HWMId), excess heat factor (EHF), wet‐bulb globe temperature (WBGT) and universal thermal climate index (UTCI), are computed. The extremeness of historical European heat waves between 1979 and 2019 using the four indices and different metrics is ranked. A normalisation to enable the comparison between the four indices is introduced. Additionally, a method to quantify the influence of the input parameters on heat wave magnitude is introduced. The spatio‐temporal behaviour of heat waves is assessed by spatial–temporal tracking. The areal extent, large‐scale intensity and duration are visualized using bubble plots. As expected, temperature explains the largest variance in all indices, but humidity is nearly as important in WBGT and wind speed plays a substantial role in UTCI. While the 2010 Russian heat wave is by far the most extreme event in duration and intensity in all normalized indices, the 2018 heat wave was comparable in size for EHF, WBGT and UTCI. Interestingly, the well‐known 2003 central European heat wave was only the fifth and tenth strongest in cumulative intensity in WBGT and UTCI, respectively. The June and July 2019 heat waves were very intense, but short‐lived, thus not belonging to the top heat waves in Europe when duration and areal extent are taken into account. Overall, the proposed normalized indices and the multi‐metric assessment of large‐scale heat waves allow for a more robust description of their extremeness and will be helpful to assess heat waves worldwide and in climate projections.

Description: Europe has been affected by record‐breaking heat waves in recent decades. Using station data and a gridded reanalysis, the extremeness of European heat waves between 1979 and 2019 is ranked using four indices: heat wave magnitude index daily (HWMId), excess heat factor (EHF), wet‐bulb globe temperature (WBGT) and universal thermal climate index (UTCI). In order to assess heatwaves worldwide and in climate projections, the spatial extent, large‐scale intensity and duration of heatwaves are visualized using bubble plots.

Description: AXA Research Fund http://dx.doi.org/10.13039/501100001961

Description: Bundesministerium für Bildung und Forschung http://dx.doi.org/10.13039/501100002347

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: Karlsruher Institut für Technologie http://dx.doi.org/10.13039/100009133

Description: Author Posting. © American Meteorological Society, 2022. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Physical Oceanography 52(8), (2022): 1927-1943, https://doi.org/10.1175/jpo-d-21-0124.1.

Description: The Galápagos Archipelago lies on the equator in the path of the eastward flowing Pacific Equatorial Undercurrent (EUC). When the EUC reaches the archipelago, it upwells and bifurcates into a north and south branch around the archipelago at a latitude determined by topography. Since the Coriolis parameter (f) equals zero at the equator, strong velocity gradients associated with the EUC can result in Ertel potential vorticity (Q) having sign opposite that of planetary vorticity near the equator. Observations collected by underwater gliders deployed just west of the Galápagos Archipelago during 2013–16 are used to estimate Q and to diagnose associated instabilities that may impact the Galápagos Cold Pool. Estimates of Q are qualitatively conserved along streamlines, consistent with the 2.5-layer, inertial model of the EUC by Pedlosky. The Q with sign opposite of f is advected south of the Galápagos Archipelago when the EUC core is located south of the bifurcation latitude. The horizontal gradient of Q suggests that the region between 2°S and 2°N above 100 m is barotropically unstable, while limited regions are baroclinically unstable. Conditions conducive to symmetric instability are observed between the EUC core and the equator and within the southern branch of the undercurrent. Using 2-month and 3-yr averages, e-folding time scales are 2–11 days, suggesting that symmetric instability can persist on those time scales.

Description: This work was supported by the National Science Foundation (Grants OCE-1232971 and OCE-1233282), the NASA Earth and Space Science Fellowship Program (Grant 80NSSC17K0443), and the Global Ocean Monitoring and Observing Program of the National Oceanographic and Atmospheric Administration (NA13OAR4830216). Color maps are from Thyng et al. (2016).

Description: 2023-02-01

Description: Author Posting. © American Meteorological Society, 2022. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of the Atmospheric and Oceanic Technology 39(8), (2022): 1183-1198, https://doi.org/10.1175/jtech-d-21-0068.1.

Description: Horizontal kinematic properties, such as vorticity, divergence, and lateral strain rate, are estimated from drifter clusters using three approaches. At submesoscale horizontal length scales O(1–10)km, kinematic properties become as large as planetary vorticity f, but challenging to observe because they evolve on short time scales O(hourstodays). By simulating surface drifters in a model flow field, we quantify the sources of uncertainty in the kinematic property calculations due to the deformation of cluster shape. Uncertainties arise primarily due to (i) violation of the linear estimation methods and (ii) aliasing of unresolved scales. Systematic uncertainties (iii) due to GPS errors, are secondary but can become as large as (i) and (ii) when aspect ratios are small. Ideal cluster parameters (number of drifters, length scale, and aspect ratio) are determined and error functions estimated empirically and theoretically. The most robust method—a two-dimensional, linear least squares fit—is applied to the first few days of a drifter dataset from the Bay of Bengal. Application of the length scale and aspect-ratio criteria minimizes errors (i) and (ii), and reduces the total number of clusters and so computational cost. The drifter-estimated kinematic properties map out a cyclonic mesoscale eddy with a surface, submesoscale fronts at its perimeter. Our analyses suggest methodological guidance for computing the two-dimensional kinematic properties in submesoscale flows, given the recently increasing quantity and quality of drifter observations, while also highlighting challenges and limitations.

Description: This research was supported by the Office of Naval Research (ONR) Departmental Research Initiative ASIRI under Grant N00014-13-1-0451 (SE and AM) and Grant N00014-13-1-0477 (VH and LC). The authors thank the captain and crew of the R/V Roger Revelle, and Andrew Lucas with the Multiscale Ocean Dynamics group at the Scripps Institution for Oceanography for providing the FastCTD data collected in 2015, which was supported by ONR Grant N00014-13-1-0489, as well as Eric D’Asaro for helpful discussions and Lance Braasch for assistance with the drifter dataset. AM and SE further thank NSF (Grant OCE-I434788) and ONR (Grant N00014-16-1-2470) for support. VH and LC were additionally supported by ONR Grants N00014-15-1-2286, N00014-14-1-0183, N00014-19-1-26-91 and NOAA Global Drifter Program (GDP) Grant NA15OAR4320071.

Description: 2023-02-01

Description: Author Posting. © American Meteorological Society, 2022. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 35(17), (2022): 5465-5482, https://doi.org/10.1175/jcli-d-21-0671.1.

Description: Understanding the contribution of ocean circulation to glacial–interglacial climate change is a major focus of paleoceanography. Specifically, many have tried to determine whether the volumes and depths of Antarctic- and North Atlantic–sourced waters in the deep ocean changed at the Last Glacial Maximum (LGM; ∼22–18 kyr BP) when atmospheric pCO2 concentrations were 100 ppm lower than the preindustrial. Measurements of sedimentary geochemical proxies are the primary way that these deep ocean structural changes have been reconstructed. However, the main proxies used to reconstruct LGM Atlantic water mass geometry provide conflicting results as to whether North Atlantic–sourced waters shoaled during the LGM. Despite this, a number of idealized modeling studies have been advanced to describe the physical processes resulting in shoaled North Atlantic waters. This paper aims to critically assess the approaches used to determine LGM Atlantic circulation geometry and lay out best practices for future work. We first compile existing proxy data and paleoclimate model output to deduce the processes responsible for setting the ocean distributions of geochemical proxies in the LGM Atlantic Ocean. We highlight how small-scale mixing processes in the ocean interior can decouple tracer distributions from the large-scale circulation, complicating the straightforward interpretation of geochemical tracers as proxies for water mass structure. Finally, we outline promising paths toward ascertaining the LGM circulation structure more clearly and deeply.

Description: S.K.H. was supported by the Investment in Science Fund at WHOI and the John E. and Anne W. Sawyer Endowed Fund in Support of Scientific Staff. F.J.P. was supported by a Stanback Postdoctoral Fellowship at Caltech.

Description: 2023-02-01

Description: Mechanisms behind the phenomenon of Arctic amplification are widely discussed. To contribute to this debate, the (AC)3 project was established in 2016 (www.ac3-tr.de/). It comprises modeling and data analysis efforts as well as observational elements. The project has assembled a wealth of ground-based, airborne, shipborne, and satellite data of physical, chemical, and meteorological properties of the Arctic atmosphere, cryosphere, and upper ocean that are available for the Arctic climate research community. Short-term changes and indications of long-term trends in Arctic climate parameters have been detected using existing and new data. For example, a distinct atmospheric moistening, an increase of regional storm activities, an amplified winter warming in the Svalbard and North Pole regions, and a decrease of sea ice thickness in the Fram Strait and of snow depth on sea ice have been identified. A positive trend of tropospheric bromine monoxide (BrO) column densities during polar spring was verified. Local marine/biogenic sources for cloud condensation nuclei and ice nucleating particles were found. Atmospheric–ocean and radiative transfer models were advanced by applying new parameterizations of surface albedo, cloud droplet activation, convective plumes and related processes over leads, and turbulent transfer coefficients for stable surface layers. Four modes of the surface radiative energy budget were explored and reproduced by simulations. To advance the future synthesis of the results, cross-cutting activities are being developed aiming to answer key questions in four focus areas: lapse rate feedback, surface processes, Arctic mixed-phase clouds, and airmass transport and transformation.

Description: The formation of low stratus cloud over idealized hills is investigated using numerical model simulations. The main driver for the cloud formation is radiative cooling due to outgoing longwave radiation. Despite a purely horizontal flow, the advection terms in the prognostic equations for heat and moisture produce vertical mixing across the upper cloud edge, leading to a loss of cloud water content. This behavior is depicted via a budget analysis. More precisely, this spurious mixing is caused by the diffusive error of the advection scheme in regions where the sloping surfaces of the terrain‐following vertical coordinate intersect the cloud top. This study shows that the intensity of the (spurious) numerical diffusion depends strongly on the horizontal resolution, the order of the advection schemes, and the choice of scalar advection scheme. A large‐eddy simulation with 4‐m horizontal resolution serves as a reference. For horizontal resolutions of a few hundred meters and simulations carried out with a model setup as used in numerical weather prediction, a strong reduction of the simulated liquid‐water path is observed. In order to keep the (spurious) numerical diffusion at coarser resolutions small, at least a fifth‐order advection scheme should be used. In the present case, a weighted essentially nonoscillatory scalar advection scheme turns out to increase the numerical diffusion along a sharp cloud edge compared with an upwind scheme. Furthermore, the choice of vertical coordinate has a strong impact on the simulated liquid‐water path over orography. With a modified definition of the sigma coordinate, it is possible to produce cloud water where the classical sigma coordinate does not allow any cloud formation.

Description: Diffusive errors of the advection scheme reduce the cloud water content of low stratus over idealized hills. This is due to the terrain‐following vertical coordinate and depends strongly on the horizontal resolution. Orographic features should be represented by at least 𝒪(10) grid points and a fifth‐order advection scheme (or higher) should be used. A weighted essentially nonoscillatory scalar advection scheme increases numerical diffusion along a sharp cloud edge compared with an upwind scheme. Modifying the definition of the sigma coordinate leads to a strong gain in the simulated liquid‐water path.

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: Hans Ertel Centre for Weather Research

Description: 〈title xmlns:mml="http://www.w3.org/1998/Math/MathML"〉Abstract〈/title〉〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉In addition to their ecological importance, rivers and streams have always been used in diverse ways by humans, resulting in the development of settlements and their connected built environments along many of the world's watercourses. During heavy rainfall, buildings, traffic infrastructure and water‐related infrastructure are exposed to potential hazards in the form of (flash) floods. In contrast to near‐natural watercourses, anthropogenically modified channels in urban areas are particularly susceptible to damage by flooding. Previous damage assessments have highlighted the need to forecast such damage to watercourses in order to identify critical areas and justify the selection and expansion of adaptation measures. Within the scope of the current study, we have developed a method based on the hydro‐morphological properties of watercourses to make transferable estimates of the economic damage potential based on ecologically‐relevant parameters. Using a scale‐specific cause‐effect analysis, we have identified characteristics of the watercourse type and adjacent structures as well as construction‐related properties of reinforcements that can increase the damage potential during flooding. In this way, we are able to show that several influencing factors determine the vulnerability of watercourses: in addition to the specific longitudinal gradient and size (macroscale) of various watercourse types, damage‐relevant boundary conditions in watercourse sections (mesoscale) and the resistance of typical bed and bank constructions are also important, reflecting the specific structural conditions. Taking rivers in Germany and the Czech Republic as case studies, in the following, we review the local identification of critical areas and describe the necessary data management. The presented “Hydro‐morphological based Vulnerability Assessment‐Concept (HyVAC)” can contribute to the flood damage prevention at watercourses by utilizing existing basic data to the greatest possible extent and thus is suitable for preliminary investigations according to the EC Flood Risk Management Directive.〈/p〉

Description: STRIMA II

Description: EU‐funded research project

Description: 〈title xmlns:mml="http://www.w3.org/1998/Math/MathML"〉Abstract〈/title〉〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉〈italic toggle="no"〉Aeolus〈/italic〉 is the first satellite mission to acquire vertical profiles of horizontal line‐of‐sight winds globally and thus fills an important gap in the Global Observing System, most notably in the Tropics. This study explores the impact of this dataset on analyses and forecasts from the European Centre for Medium‐Range Weather Forecasts (ECMWF) and Deutscher Wetterdienst (DWD), focusing specifically on the West African Monsoon (WAM) circulation during the boreal summers of 2019 and 2020. The WAM is notoriously challenging to forecast and is characterized by prominent and robust large‐scale circulation features such as the African Easterly Jet North (AEJ‐North) and Tropical Easterly Jet (TEJ). Assimilating 〈italic toggle="no"〉Aeolus〈/italic〉 generally improves the prediction of zonal winds in both forecasting systems, especially for lead times above 24 h. These improvements are related to systematic differences in the representation of the two jets, with the AEJ‐North weakened at its southern flank in the western Sahel in the ECMWF analysis, while no obvious systematic differences are seen in the DWD analysis. In addition, the TEJ core is weakened in the ECMWF analysis and strengthened on its southern edge in the DWD analysis. The regions where the influence of 〈italic toggle="no"〉Aeolus〈/italic〉 on the analysis is greatest correspond to the Intertropical Convergence Zone (ITCZ) region for ECMWF and generally the upper troposphere for DWD. In addition, we show the presence of an altitude‐ and orbit‐dependent bias in the Rayleigh‐clear channel, which causes the zonal winds to speed up and slow down diurnally. Applying a temperature‐dependent bias correction to this channel contributes to a more accurate representation of the diurnal cycle and improved prediction of the WAM winds. These improvements are encouraging for future investigations of the influence of 〈italic toggle="no"〉Aeolus〈/italic〉 data on African Easterly Waves and associated Mesoscale Convective Systems.〈/p〉

Description: 〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉Forecasting in tropical Africa is hampered by large model errors and low availability of conventional observations. The assimilation of 〈italic〉Aeolus〈/italic〉 wind data into the operational ECMWF system leads to a consistent root‐mean‐square error (RMSE) reduction of the order of 2% in +48 h zonal wind forecasts over the region during boreal summer 2019, including the African and Tropical Easterly Jets (AEJ, TEJ) and subtropical jets (STJ). 〈boxed-text position="anchor" id="qj4442-blkfxd-0001" content-type="graphic" xml:lang="en"〉〈graphic position="anchor" id="jats-graphic-1" xlink:href="urn:x-wiley:00359009:media:qj4442:qj4442-toc-0001"〉 〈/graphic〉 〈/boxed-text〉〈/p〉

Description: https://aeolus-ds.eo.esa.int/oads/access/collection

Description: 〈title xmlns:mml="http://www.w3.org/1998/Math/MathML"〉Abstract〈/title〉〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉Owing to the complicated spatial–temporal characteristics of East Asian precipitation (EAP), climate models have limited skills in simulating the modern Asian climate. This consequently leads to large uncertainties in simulations of the past EAP variation and future projections. Here, we explore the performance of the newly developed Alfred Wegener Institute Climate Model, version 3 (AWI‐CM3) in simulating the climatological summer EAP. To test whether the model's skill depends on its atmosphere resolution, we design two AWI‐CM3 simulations with different horizontal resolutions. The result shows that both simulations have acceptable performance in simulating the summer mean EAP, generally better than the majority of individual models participating in the Coupled Model Intercomparison Project (CMIP6). However, for the monthly EAP from June to August, AWI‐CM3 exhibits a decayed skill, which is due to the subseasonal movement of the western Pacific subtropical high bias. The higher‐resolution AWI‐CM3 simulation shows an overall improvement relative to the one performed at a relatively lower resolution in all aspects taken into account regarding the EAP. We conclude that AWI‐CM3 is a suitable tool for exploring the EAP for the observational period. Having verified the model's skill for modern climate, we suggest employing the AWI‐CM3, especially with high atmosphere resolution, both for applications in paleoclimate studies and future projections.〈/p〉

Description: 〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉This figure shows the skill scores of AWI‐CM3 and CMIP6 models in simulating the climatological summer East Asian precipitation (EAP), which indicates that AWI‐CM3 simulations perform better than most CMIP6 individual models for the summer mean EAP, while AWI‐CM3's skills decay from June to August.〈boxed-text position="anchor" content-type="graphic" id="joc8075-blkfxd-0001" xml:lang="en"〉 〈graphic position="anchor" id="jats-graphic-1" xlink:href="urn:x-wiley:08998418:media:joc8075:joc8075-toc-0001"〉 〈alt-text〉image〈/alt-text〉 〈/graphic〉 〈/boxed-text〉〈/p〉

Description: National Natural Science Foundation of China http://dx.doi.org/10.13039/501100001809

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: Helmholtz Climate Initiative REKLIM

Description: Helmholtz Program

Description: Bundesministerium für Bildung und Forschung http://dx.doi.org/10.13039/501100002347

Description: China Scholarship Council http://dx.doi.org/10.13039/501100004543

Description: https://opendata.dwd.de/climate_environment/GPCC/html/fulldata-monthly_v2022_doi_download.html

Description: https://crudata.uea.ac.uk/cru/data/hrg/cru_ts_4.05

Description: http://aphrodite.st.hirosaki-u.ac.jp/products.html

Description: https://jra.kishou.go.jp/JRA-55/index_en.html

Description: https://esgf-node.llnl.gov/search/cmip6

Description: 〈title xmlns:mml="http://www.w3.org/1998/Math/MathML"〉Abstract〈/title〉〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉In this study, a new multilayer urban canopy parameterization for high‐resolution (∼1 km) atmospheric models using the nudging approach to represent the impacts of urban canopies on airflow is presented. In our parameterization, a nudging term is added to the momentum equations and a source term to the turbulent kinetic energy equation to account for building effects. The challenge of this parameterization lies in defining appropriate values for the nudging coefficient and the weighting function used to reflect canopy effects. Values of both are derived and the parameterization developed is implemented and tested for idealized cases in the Mesoscale Transport and Stream model (METRAS). Comparison data are taken from obstacle‐resolving microscale model results. Results show that the parameterization using the nudging approach can simulate aerodynamic effects induced within the canopy by obstacles well, in terms of reduction of wind speeds and production of additional turbulent kinetic energy. Thus, models with existing nudging can use this approach as an efficient and effective method to parameterize dynamic urban canopy effects.〈/p〉

Description: 〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉In this study, a new multilayer urban canopy parameterization for high‐resolution (∼1 km) atmospheric models using the nudging approach to represent the impacts of urban canopies on airflow is presented. Results show that the parameterization developed can simulate aerodynamic effects induced within the canopy by obstacles well, in terms of reduction of wind speeds and production of additional turbulent kinetic energy. Models with existing nudging can use this approach as an efficient and effective method to parameterize dynamic urban canopy effects. 〈boxed-text position="anchor" id="qj4524-blkfxd-0001" content-type="graphic" xml:lang="en"〉〈graphic position="anchor" id="jats-graphic-1" xlink:href="urn:x-wiley:00359009:media:qj4524:qj4524-toc-0001"〉 〈/graphic〉 〈/boxed-text〉〈/p〉

Description: Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany's Excellence Strategy‐EXC 2037 'CLICCS‐Climate, Climatic Change, and Society'

Description: 〈title xmlns:mml="http://www.w3.org/1998/Math/MathML"〉Abstract〈/title〉〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉Reliable prediction of heavy precipitation events causing floods in a world of changing climate is crucial for the development of appropriate adaption strategies. Many attempts to provide such predictions have already been conducted but there is still much potential for improvement left. This is particularly true for statistical downscaling of heavy precipitation due to changes present in the corresponding atmospheric drivers. In this study, a circulation pattern (CP) conditional downscaling to the station level is proposed which considers occurring frequency changes of CPs. Following a strict circulation‐to‐environment approach we use atmospheric predictors to derive CPs. Subsequently, precipitation observations are used to derive CP conditional cumulative distribution functions (CDFs) of daily precipitation. Raw precipitation time series are sampled from these CDFs. Bias correction is applied to the sampled time series with quantile mapping (QM) and parametric transfer functions (PTFs) as methods being tested. The added value of this CP conditional downscaling approach is evaluated against the corresponding common non‐CP conditional approach. The performance evaluation is conducted by using Kling–Gupta Efficiency (KGE), root mean squared error (RMSE), and mean absolute error (MAE) metrics. In both cases the applied bias correction is identical. Potential added value can therefore only be attributed to the CP conditioning. It can be shown that the proposed CP conditional downscaling approach is capable of yielding more reliable and accurate downscaled daily precipitation time series in comparison to a non‐CP conditional approach. This can be seen in particular for the extreme parts of the distribution. Above the 95th percentile, an average performance gain of +0.24 and a maximum gain of +0.6 in terms of KGE is observed. These findings support the assumption of conserving and utilizing atmospheric information through CPs can be beneficial for more reliable statistical precipitation downscaling. Due to the availability of these atmospheric predictors in climate model output, the presented method is potentially suitable for downscaling precipitation projections.〈/p〉

Description: Bundesministerium für Bildung und Forschung http://dx.doi.org/10.13039/501100002347

Description: https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-pressure-levels?tab=overview

Description: https://cdc.dwd.de/portal/

Description: 〈title xmlns:mml="http://www.w3.org/1998/Math/MathML"〉Abstract〈/title〉〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉Riverbed clogging is key to assessing vertical connectivity in the hyporheic zone and is often quantified using single‐parameter or qualitative approaches. However, clogging is driven by multiple, interacting physical and bio‐geochemical parameters, which do not allow for a conclusive assessment of hyporheic connectivity with single‐parameter approaches. In addition, existing qualitative assessments lack transparency and repeatability. This study introduces a Multi‐Parameter Approach to quantify Clogging and vertical hyporheic connectivity (MultiPAC), which builds on standardized measurements of physical (grain size characteristics, porosity, hydraulic conductivity) and bio‐geochemical (interstitial dissolved oxygen) parameters. We apply MultiPAC at three gravel‐bed rivers and show how the set of parameters provides a representative appreciation of physical riverbed clogging, thus quantifying vertical hyporheic connectivity. However, more parameters are required to fully characterize biological clogging. In addition, MultiPAC locates clogged layers in the hyporheic zone through multi‐parameter vertical profiles over the riverbed depth. The discussion outlines the relevance of MultiPAC to guide field surveys.〈/p〉

Description: https://github.com/Ecohydraulics/kf-converter-w-flopy

Description: 〈title xmlns:mml="http://www.w3.org/1998/Math/MathML"〉Abstract〈/title〉〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉Projected changes in summer precipitation deficits partly depend on alterations in synoptic circulations. Here, the automated Jenkinson–Collison classification is used to assess the ability of 21 global climate models (GCMs) to capture the frequency of recurring circulation types (CTs) and their implications for European daily precipitation amounts in summer (JJA). The ability of the GCMs to reproduce the observed present‐day climate features is evaluated first. Most GCMs capture the observed links between the mean CTs directional flow characteristics and the occurrence of dry days and related dry months. The most robust relationships are found for anticyclonic and easterly CTs which are generally associated with higher‐than‐average occurrences of dry conditions. Future changes in summer CTs' frequencies are estimated in the high‐emission SSP5‐8.5 scenario for the sake of a high signal‐to‐noise ratio. Our results reveal consistent changes, mainly in the zonal CTs. A robust decrease in frequency of the westerlies and an increase in the frequency of easterly CTs favour more continental, dry and warm air masses over central Europe. These dynamical changes are shown to enhance the projected summer drying over central and southern Europe.〈/p〉

Description: 〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉Summer large‐scale circulations are derived over Europe using an automated classification. Spatial characteristics of the patterns and their influence on dry days are investigated. Future changes are explored based on global climate models. The predicted drier summers in Europe are found to be influenced by consistent changes in west‐easterly circulations.〈boxed-text position="anchor" content-type="graphic" id="joc8033-blkfxd-0001" xml:lang="en"〉 〈graphic position="anchor" id="jats-graphic-1" xlink:href="urn:x-wiley:08998418:media:joc8033:joc8033-toc-0001"〉 〈/graphic〉 〈/boxed-text〉〈/p〉

Description: EU International Training Network (ITN) Climate Advanced Forecasting of sub‐seasonal Extremes (CAFE)

Description: H2020 Marie Skłodowska‐Curie Actions

Description: https://github.com/PedroLormendez/jcclass

Description: Wind energy is a key option in global dialogues about climate change mitigation. Here, we combined observations from surface wind stations, reanalysis datasets, and state‐of‐the‐art regional climate models from the Coordinated Regional Climate Downscaling Experiment (CORDEX Africa) to study the current and future wind energy potential in Zambia. We found that winds are dominated by southeasterlies and are rarely strong with an average speed of 2.8 m·s〈sup〉−1〈/sup〉. When we converted the observed surface wind speed to a turbine hub height of 100 m, we found a ~38% increase in mean wind speed for the period 1981–2000. Further, both simulated and observed wind speed data show statistically significant increments across much of the country. The only areas that divert from this upward trend of wind speeds are the low land terrains of the Eastern Province bordering Malawi. Examining projections of wind power density (WPD), we found that although wind speed is increasing, it is still generally too weak to support large‐scale wind power generation. We found a meagre projected annual average WPD of 46.6 W·m〈sup〉−2〈/sup〉. The highest WPDs of ~80 W·m−2 are projected in the northern and central parts of the country while the lowest are to be expected along the Luangwa valley in agreement with wind speed simulations. On average, Zambia is expected to experience minor WPD increments of 0.004 W·m〈sup〉−2〈/sup〉 per year from 2031 to 2050. We conclude that small‐scale wind turbines that accommodate cut‐in wind speeds of 3.8 m·s〈sup〉−1〈/sup〉 are the most suitable for power generation in Zambia. Further, given the limitations of small wind turbines, they are best suited for rural and suburban areas of the country where obstructions are few, thus making them ideal for complementing the government of the Republic of Zambia's rural electrification efforts.

Description: Wind speed is very slow in Zambia. It is increasing but remains unlikely to support large commercial wind farms especially not at the 100‐m hub height. Any efforts for wind power generation in Zambia should be towards ultra‐tall wind turbines fitted with larger rotors.

Description: Alexander von Humboldt‐Stiftung

Description: 〈title xmlns:mml="http://www.w3.org/1998/Math/MathML"〉Abstract〈/title〉〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉The Arctic near‐surface air temperature increases most strongly during the cold season, and ocean heat storage has often been cited as a crucial component in linking the ice‐albedo radiative feedback, which is active in summer, and near‐surface air temperature increase in winter, when the lapse rate feedback contributes to Arctic warming. Here, we first estimate how much local heat storage and ocean heat transport contribute to net surface energy fluxes on a seasonal scale in CMIP6 models. We then compare contributions in a base state under weak anthropogenic forcing to a near‐present‐day state in which significant Arctic amplification is simulated. Our analysis indicates that, in a few regions, ocean heat transport plays a larger role for cold‐season net surface energy fluxes compared with local heat storage. Analyzing differences between past and near‐present‐day conditions suggests that the lapse rate feedback, which mainly acts during the cold season in warm water inflow regions, may be more strongly influenced than previously thought by increased ocean heat transport from lower latitudes.〈/p〉

Description: 〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉Arctic Ocean net upward surface energy fluxes in the cold season were decomposed into contributions from local heat storage (yellow, see schematic) and ocean heat transport (red). Our analysis of CMIP6 model output suggests that, in a few inflow regions, ocean heat transport contributes more to cold‐season net surface energy fluxes compared with local heat storage. In parts of these inflow regions, the relative contribution of ocean heat transport increased with time. 〈boxed-text position="anchor" id="qj4496-blkfxd-0001" content-type="graphic" xml:lang="en"〉〈graphic position="anchor" id="jats-graphic-1" xlink:href="urn:x-wiley:00359009:media:qj4496:qj4496-toc-0001"〉 〈/graphic〉 〈/boxed-text〉〈/p〉

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: https://esgf-node.llnl.gov/projects/cmip6/

Description: https://nsidc.org/data/g10010

Description: The novel Aeolus satellite, which carries the first Doppler wind lidar providing profiles of horizontal line‐of‐sight (HLOS) winds, addresses a significant gap in direct wind observations in the global observing system. The gap is particularly critical in the tropical upper troposphere and lower stratosphere (UTLS). This article validates the Aeolus Rayleigh–clear wind product and short‐range forecasts of the European Centre for Medium‐Range Weather Forecasts (ECMWF) with highly accurate winds from the Loon super pressure balloon network at altitudes between 16 and 20 km. Data from 229 individual balloon flights are analysed, applying a collocation criterion of 2 hr and 200 km. The comparison of Aeolus and Loon data shows systematic and random errors of -0.31 and 6.37 m·s〈sup〉-1〈/sup〉, respectively, for the Aeolus Rayleigh–clear winds. The horizontal representativeness error of Aeolus HLOS winds (nearly the zonal wind component) in the UTLS ranges from 0.6–1.1 m·s〈sup〉-1〈/sup〉 depending on the altitude. The comparison of Aeolus and Loon datasets against ECMWF model forecasts suggests that the model systematically underestimates the HLOS winds in the tropical UTLS by about 1 m·s〈sup〉-1〈/sup〉. While Aeolus winds are currently considered as point winds by the ECMWF data assimilation system, the results of the present study demonstrate the need for a more realistic HLOS wind observation operator for assimilating Aeolus winds.

Description: El Niño–Southern Oscillation (ENSO) is one of the most important modes of climate variability on interannual timescales. We aim to find out whether a change in ENSO frequency can be predicted for the nearer future. We analyse the unforced pre‐industrial control run and the forced 1%/year CO〈sub〉2〈/sub〉 increase run for an ensemble of 43 general circulation models that participated in the Coupled Model Intercomparison Project Phase 6 (CMIP6). We assume that the uncertainty of ENSO frequency trend estimates from an ensemble is caused by apparent trends as well as model differences. The part of the uncertainty caused by apparent trends is estimated from the pre‐industrial control simulations. As a measure for ENSO frequency, we use the number of El Niño‐ and La Niña‐like months in a moving 30‐year time window. Its linear decadal trend is calculated for every member. The multimember mean of the trend for both experiments is less than 0.7 events per decade. Given that the standard error is of the same order of magnitude, we consider this a negligible trend. The uncertainties are large in both experiments and we can attribute most of the intermember variability to apparent trends due to natural variability rather than different model reactions to CO〈sub〉2〈/sub〉 forcing. This means that the impact of intermodel differences might have been overstated in previous studies. Apparent trends make it very difficult to make reliable predictions of changes in ENSO frequency based on 120‐year time series.

Description: The 1pctCO2 and piControl ensembles from CMIP6 are analysed for 43 models with a focus on changes in ENSO frequency. We find that most of the intermember variability can be attributed to natural variability instead of model differences. Therefore, the uncertainty can only marginally be reduced and it is very difficult to reliably predict changes in ENSO frequency on a timescale of 150 years.

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: https://doi.org/10.5281/zenodo.6841964

Description: 〈title xmlns:mml="http://www.w3.org/1998/Math/MathML"〉Abstract〈/title〉〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉The transboundary region of the Iishana system in the western Cuvelai Basin, between southern Angola and northern Namibia, is frequently affected by floods at irregular intervals. As a result, the predominantly rural, subsistence farming population has experienced crop failures, human, and economic losses. To date, very little is known about the generation of floods, flood concentration, and stormwater drainage dynamics in this region. In this study, 2D‐hydrodynamic modeling was applied to reconstruct one of the latest major flood events during the rainy season from November 2008 to March 2009 in order to study the runoff behavior and interconnectivity of the Iishana system. The model focused on the eastern part of the Iishana system, which was most affected by floods and flood damage due to the high population density in and around Oshakati, the regional capital. Two main streams were identified noteworthy because they merge and subsequently affect Oshakati. Regarding the simulated flood event water depths vary from 0.1 m to 14 m, with an average of 0.2 m, while water depths above 5 m were attributed to borrow pits. The inundation area ranged up to 1860 km〈sup〉2〈/sup〉 and the amount of water left after the rainy season on March 25th, 2009, was determined between 0.116 and 0.547 km〈sup〉3〈/sup〉, depending on the amount of evapotranspiration considered in the model. Thus, in the Angolan part of the Iishana system, significantly larger quantities of water are available for longer periods of time during the subsequent dry season, whereas the system in Namibia stores less water, resulting in a shorter water retention period.〈/p〉

Description: Deutsche Hydrologische Gesellschaft (DHG)

Description: Freie Universität Berlin http://dx.doi.org/10.13039/501100007537

Description: https://dx.doi.org/10.17169/refubium-35737

Description: Heat stress is projected to intensify with global warming, causing significant socioeconomic impacts and threatening human health. Wet-bulb temperature (WBT), which combines temperature and humidity effects, is a useful indicator for assessing regional and global heat stress variability and trends. However, the variations of European WBT and their underlying mechanisms remain unclear. Using observations and reanalysis datasets, we demonstrate a remarkable warming of summer WBT during the period 1958–2021 over Europe. Specifically, the European summer WBT has increased by over 1.08C in the past 64 years. We find that the increase in European summer WBT is driven by both near-surface warming temperatures and increasing atmospheric moisture content. We identify four dominant modes of European summer WBT variability and investigate their linkage with the large-scale atmospheric circulation and sea surface temperature anomalies. The first two leading modes of the European WBT variability exhibit prominent interdecadal to long-term variations, mainly driven by a circumglobal wave train and concurrent sea surface temperature variations. The last two leading modes of European WBT variability mainly show interannual variations, indicating a direct and rapid response to large-scale atmospheric dynamics and nearby sea surface temperature variations. Further analysis shows the role of global warming and changes in midlatitude circulations in the variations of summer WBT. Our findings can enhance the understanding of plausible drivers of heat stress in Europe and provide valuable insights for regional decision-makers and climate adaptation planning.

Description: 〈title xmlns:mml="http://www.w3.org/1998/Math/MathML"〉Abstract〈/title〉〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉The prediction skill of sub‐seasonal forecast models is evaluated for seven year‐round weather regimes in the Atlantic–European region. Reforecasts based on models from three prediction centers are considered and verified against weather regimes obtained from ERA‐Interim reanalysis. Results show that predicting weather regimes as a proxy for the large‐scale circulation outperforms the prediction of raw geopotential height. Greenland blocking tends to have the longest year‐round skill horizon for all three models, especially in winter. On the other hand, the skill is lowest for the European blocking regime for all three models, followed by the Scandinavian blocking regime. Furthermore, all models struggle to forecast flow situations that cannot be assigned to a weather regime (so‐called no regime), in comparison with weather regimes. Related to this, variability in the occurrence of no regime, which is most frequent in the transition seasons, partly explains the predictability gap between transition seasons and winter and summer. We also show that models have difficulties in discriminating between related regimes. This can lead to misassignments in the predicted regime during flow situations in which related regimes manifest. Finally, we document the changes in skill between model versions, showing important improvements for the ECMWF and NCEP models. This study is the first multi‐model assessment of year‐round weather regimes in the Atlantic–European domain. It advances our understanding of the predictive skill for weather regimes, reveals strengths and weaknesses of each model, and thus increases our confidence in the forecasts and their usefulness for decision‐making.〈/p〉

Description: 〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉This study is the first sub‐seasonal multi‐model assessment of seven year‐round weather regimes in the Atlantic–European domain. Greenland blocking tends to have the longest year‐round skill horizon for all models, especially in winter. The skill is lowest for the European blocking regime for all models, followed by Scandinavian blocking. Variability in the occurrence of no regime partly explains the predictability gap between the transition seasons and winter and summer. 〈boxed-text position="anchor" id="qj4512-blkfxd-0001" content-type="graphic" xml:lang="en"〉〈graphic position="anchor" id="jats-graphic-1" xlink:href="urn:x-wiley:00359009:media:qj4512:qj4512-toc-0001"〉 〈/graphic〉 〈/boxed-text〉〈/p〉

Description: Helmholtz Association http://dx.doi.org/10.13039/501100001656

Description: AXPO Solutions AGN/A

Description: 〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉Flow‐dependent errors in tropical analyses and short‐range forecasts are analysed using global observing‐system simulation experiments assimilating only temperature, only winds, and both data types using the ensemble Kalman filter (EnKF) Data Assimilation Research Testbed (DART) and a perfect model framework. The idealised, homogeneous observation network provides profiles of wind and temperature data from the nature run for January 2018 using the National Center for Atmospheric Research (NCAR) Community Earth System Model (CESM) forced by the observed sea‐surface temperature. The results show that the assimilation of abundant wind observations in a perfect model makes the temperature data in the Tropics largely uninformative. Furthermore, the assimilation of wind data reduces the background errors in specific humidity twice as much as the assimilation of temperature observations. In all experiments, the largest analysis uncertainties and the largest short‐term forecast errors are found in regions of strong vertical and longitudinal gradients in the background wind, especially in the upper troposphere and lower stratosphere over the Indian Ocean and Maritime Continent. The horizontal error correlation scales are on average short throughout the troposphere, just several hundred km. The correlation scales of the wind variables in precipitating regions are half of those in nonprecipitating regions. In precipitating regions, the correlations are elongated vertically, especially for the wind variables. Strong positive cross‐correlations between temperature and specific humidity in the precipitating regions are explained using the Clausius–Clapeyron equation.〈/p〉

Description: China Scholarship Council http://dx.doi.org/10.13039/501100004543

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: This work focuses on the potential of a network of Doppler lidars for the improvement of short‐term forecasts of low‐level wind. For the impact assessment, we developed a new methodology that is based on ensemble sensitivity analysis (ESA). In contrast to preceding network design studies using ESA, we calculate the explicit sensitivity including the inverse of the background covariance B matrix to account directly for the localization scale of the assimilation system. The new method is applied to a pre‐existing convective‐scale 1,000‐member ensemble simulation to mitigate effects of spurious correlations. We evaluate relative changes in the variance of a forecast metric, that is, the low‐level wind components averaged over the Rhein–Ruhr metropolitan area in Germany. This setup allows us to compare the relative variance change associated with the assimilation of hypothetical observations from a Doppler wind lidar with respect to the assimilation of surface‐wind observations only. Furthermore, we assess sensitivities of derived variance changes to a number of settings, namely observation errors, localization length scale, regularization factor, number of instruments in the network, and their location, as well as data availability of the lidar measurements. Our results demonstrate that a network of 20–30 Doppler lidars leads to a considerable variance reduction of the forecast metric chosen. On average, an additional network of 25 Doppler lidars can reduce the 1–3 hr forecast error by a factor of 1.6–3.3 with respect to 10‐m wind observations only. The results provide the basis for designing an operational network of Doppler lidars for the improvement of short‐term low‐level wind forecasts that could be especially valuable for the renewable energy sector.

Description: This study presents the potential of a Doppler lidar network to improve short‐term low‐level wind forecasts. The approach used in this study does not require real observations and can provide valuable information for designing an operational network. The study is based on a convective‐scale 1,000‐member ensemble simulation over Germany. The results show that Doppler lidars lead to considerable variance reduction and should be considered for future observational networks.

Description: Hans‐Ertel‐Centre for Weather Research funded by the German Federal Ministry for Transportation and Digital Infrastructure

Description: https://doi.org/10.5281/zenodo.6331758

Description: Land surface heterogeneity in conjunction with ambient winds influences the convective atmospheric boundary layer by affecting the distribution of incoming solar radiation and forming secondary circulations. This study performed coupled large‐eddy simulation (ICON‐LEM) with a land surface model (TERRA‐ML) over a flat river corridor mimicked by soil moisture heterogeneity to investigate the impact of ambient winds on secondary circulations. The coupled model employed double‐periodic boundary conditions with a spatial scale of 4.8 km. All simulations used the same idealized initial atmospheric conditions with constant incident radiation of 700 W⋅m〈sup〉−2〈/sup〉 and various ambient winds with different speeds (0 to 16 m⋅s〈sup〉−1〈/sup〉) and directions (e.g., cross‐river, parallel‐river, and mixed). The atmospheric states are decomposed into ensemble‐averaged, mesoscale, and turbulence. The results show that the secondary circulation structure persists under the parallel‐river wind conditions independently of the wind speed but is destroyed when the cross‐river wind is stronger than 2 m⋅s〈sup〉−1〈/sup〉. The soil moisture and wind speed determine the influence on the surface energy distribution independent of the wind direction. However, secondary circulations increase advection and dispersive heat flux while decreasing turbulent energy flux. The vertical profiles of the wind variance reflect the secondary circulation, and the maximum value of the mesoscale vertical wind variance indicates the secondary circulation strength. The secondary circulation strength positively scales with the Bowen ratio, stability parameter (−Z〈sub〉i〈/sub〉/L), and thermal heterogeneity parameter under cross‐river wind and mixed wind conditions. The proposed similarity analyses and scaling approach provide a new quantitative perspective on the impact of the ambient wind under heteronomous soil moisture conditions on secondary circulation.

Description: Deutsche Forschungsgemeinschaft http://dx.doi.org/10.13039/501100001659

Description: 〈title xmlns:mml="http://www.w3.org/1998/Math/MathML"〉Abstract〈/title〉〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉The usually short lifetime of convective storms and their rapid development during unstable weather conditions makes forecasting these storms challenging. It is necessary, therefore, to improve the procedures for estimating the storms' expected life cycles, including the storms' lifetime, size, and intensity development. We present an analysis of the life cycles of convective cells in Germany, focusing on the relevance of the prevailing atmospheric conditions. Using data from the radar‐based cell detection and tracking algorithm KONRAD of the German Weather Service, the life cycles of isolated convective storms are analysed for the summer half‐years from 2011 to 2016. In addition, numerous convection‐relevant atmospheric ambient variables (e.g., deep‐layer shear, convective available potential energy, lifted index), which were calculated using high‐resolution COSMO‐EU assimilation analyses (0.0625°), are combined with the life cycles. The statistical analyses of the life cycles reveal that rapid initial area growth supports wider horizontal expansion of a cell in the subsequent development and, indirectly, a longer lifetime. Specifically, the information about the initial horizontal cell area is the most important predictor for the lifetime and expected maximum cell area during the life cycle. However, its predictive skill turns out to be moderate at most, but still considerably higher than the skill of any ambient variable is. Of the latter, measures of midtropospheric mean wind and vertical wind shear are most suitable for distinguishing between convective cells with short lifetime and those with long lifetime. Higher thermal instability is associated with faster initial growth, thus favouring larger and longer living cells. A detailed objective correlation analysis between ambient variables, coupled with analyses discriminating groups of different lifetime and maximum cell area, makes it possible to gain new insights into their statistical connections. The results of this study provide guidance for predictor selection and advancements of nowcasting applications.〈/p〉

Description: 〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉Based on a combination of data of the cell tracking algorithm KONRAD of the German Weather Service and COSMO‐EU model analyses for the summer half‐years from 2011 to 2016, statistical relationships between storm attributes (lifetime and maximum horizontal area), and ambient variables as well as the storms' history are quantified. The initial growth of the cell area is a better indicator of the lifetime and maximum area than ambient variables are. Of the latter, measures of the midtropospheric wind and vertical wind shear, in particular, are most suitable for distinguishing between convective cells with short and long lifetimes, whereas higher convective instability favours larger cells. 〈boxed-text position="anchor" id="qj4505-blkfxd-0001" content-type="graphic" xml:lang="en"〉〈graphic position="anchor" id="jats-graphic-1" xlink:href="urn:x-wiley:00359009:media:qj4505:qj4505-toc-0001"〉 〈/graphic〉 〈/boxed-text〉〈/p〉

Description: Bundesministerium für Digitales und Verkehr http://dx.doi.org/10.13039/100008383

Description: A fundamental statistic of climate variability is its spatiotemporal correlation function. Its complex structure can be concisely summarized by a frequency-dependent measure of the effective spatial degrees of freedom (ESDOF). Here we present, for the first time, frequency-dependent ESDOF estimates of global natural surface temperature variability from purely instrumental measurements, using the HadCRUT4 dataset (1850-2014). The approach is based on a newly developed method for estimating the frequency-dependent spatial correlation function from gappy data fields. Results reveal a multicomponent structure of the spatial correlation function, including a large-amplitude short-distance component (with weak time scale dependence) and a small-amplitude long-distance component (with increasing relative amplitude toward the longer time scales). Two frequency-dependent ESDOF measures are applied, each responding mainly to either of the two components. Both measures exhibit a significant ESDOF reduction from monthly to multidecadal time scales, implying an increase of the effective spatial scale of natural surface temperature fluctuations. Moreover, it is found that a good approximation to the global number of equally spaced samples needed to estimate the variance of global mean temperature is given, at any frequency, by the greater one of the two ESDOF measures, decreasing from ;130 at monthly to ;30 at multidecadal time scales. Finally, the multicomponent structure of the correlation function together with the detected ESDOF scaling properties indicate that the ESDOF reduction toward the longer time scales cannot be explained simply by diffusion acting on stochastically driven anomalies, as it might be suggested f rom simple stochastic-diffusive energy balance models.

Description: Climate variability occurs over wide ranges of spatial and temporal scales. It exhibits a complex spatial covariance structure, which depends on geographic location (e.g., tropics vs extratropics) and also consists of a superposition of (i) components with gradually decaying positive correlation functions and (ii) teleconnections that often involve anticorrelations. In addition, there are indications that the spatial covariance structure depends on frequency. Thus, a comprehensive assessment of the spatiotemporal covariance structure of climate variability would require an extensive set of statistical diagnostics. Therefore, it is often desirable to characterize the covariance structure by a simple summarizing metric that is easy to compute from datasets. Such summarizing metrics are useful, for example, in the context of comparisons between climate models or between models and observations. Here we introduce a frequency-dependent version of a simple measure of the effective spatial degrees of freedom. The measure is based on the temporal variance of the global average of some climate variable, and its novel aspect consists in its frequency dependence. We also provide a clear geometric interpretation of the measure. Its easy applicability is demonstrated using near-surface temperature and precipitation fields obtained from a paleoclimate model simulation. This application reveals a distinct scaling behavior of the spatial degrees of freedom as a function of frequency, ranging from monthly to millennial scales.

Description: Coastal upwelling, driven by alongshore winds and characterized by cold sea surface temperatures and high upper-ocean nutrient content, is an important physical process sustaining some of the oceans’ most productive ecosystems. To fully understand the ocean properties in eastern boundary upwelling systems, it is important to consider the depth of the source waters being upwelled, as it affects both the SST and the transport of nutrients toward the surface. Here, we construct an upwelling source depth distribution for parcels at the surface in the upwelling zone. We do so using passive tracers forced at the domain boundary for every model depth level to quantify their contributions to the upwelled waters. We test the dependence of this distribution on the strength of the wind stress and stratification using high-resolution regional ocean simulations of an idealized coastal upwelling system. We also present an efficient method for estimating the mean upwelling source depth. Furthermore, we show that the standard deviation of the upwelling source depth distribution increases with increasing wind stress and decreases with increasing stratification. These results can be applied to better understand and predict how coastal upwelling sites and their surface properties have and will change in past and future climates.

Description: 〈title xmlns:mml="http://www.w3.org/1998/Math/MathML"〉Abstract〈/title〉〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉As projected by multiple climate models, short‐duration heavy precipitation events (SDHPEs) are expected to intensify particularly quickly under the changing climate posing substantial risk to natural and human systems. Yet over the years, SDHPEs have received less scientific attention than long‐duration heavy precipitation events (LDHPEs), mainly due to the limitations of measurement systems. Our aim is to provide insight into spatial and temporal variability of SDHPEs detected by the radar network of the 〈italic toggle="no"〉Deutscher Wetterdienst〈/italic〉 (DWD) in Germany from 2001 to 2020 as well as to explore their links to circulation patterns (CPs). The study is based on the Catalogue of Radar‐based heavy Rainfall Events (CatRaRE) generated using reprocessed gauge‐adjusted data of the DWD radar network as well as a new numerical method for classifying CPs over Central Europe called “〈italic toggle="no"〉Großwetterlagen〈/italic〉 for Reanalyses” (GWL‐REA). The results have demonstrated that SDHPEs, which are defined based on either locally valid precipitation values with a return period of 5 years (CatRaRE T5) or absolute precipitation values equal to DWD Warning Level 3 (CatRaRE W3), are common phenomena occurring most frequently in the afternoon hours of the summer season. They constitute up to 90% of all heavy precipitation events included in the catalogues covering relatively small areas—the median area of SDHPEs ranges from 22 km〈sup〉2〈/sup〉 (CatRaRE T5) to 24 km〈sup〉2〈/sup〉 (CatRaRE W3), while the median area of LDHPEs ranges from 175 km〈sup〉2〈/sup〉 (CatRaRE W3) to 184 km〈sup〉2〈/sup〉 (CatRaRE T5). As compared to LDHPEs, SDHPEs are generated by a wider spectrum of circulation conditions, including not only cyclonic but also anticyclonic CPs. In the warm season, the anticyclonic CPs, often accompanied by air mass advection from the south, can induce high thermal instability leading to the development of relatively small, isolated convective cells, which often cannot be captured by rain gauge stations.〈/p〉

Description: Federal Ministry for Digital and Transport (BMDV)

Description: https://www.dwd.de/DE/leistungen/catrare/catrare.html

Description: 〈title xmlns:mml="http://www.w3.org/1998/Math/MathML"〉Abstract〈/title〉〈p xmlns:mml="http://www.w3.org/1998/Math/MathML" xml:lang="en"〉The Intertropical Convergence Zone (ITCZ) is a central component of the atmospheric general circulation, but remarkably little is known about the dynamical and thermodynamical structure of the convergence zone itself. This is true even for the structure of the low‐level convergence that gives the ITCZ its name. Following on from the major international field campaigns in the 1960s and 1970s, we performed extensive atmospheric profiling of the Atlantic ITCZ during a ship‐based measurement campaign aboard the research vessel 〈italic toggle="no"〉SONNE〈/italic〉 in summer 2021. Combining data collected during our north–south crossing of the ITCZ with reanalysis data shows the ITCZ to be a meridionally extended region of intense precipitation, with enhanced surface convergence at its edges rather than in the center. Based on the location of these edges, we construct a composite view of the structure of the Atlantic ITCZ. The ITCZ, far from being simply a region of enhanced deep convection, has a rich inner life, that is, a rich dynamical and thermodynamic structure that changes throughout the course of the year, and has a northern edge that differs systematically from the southern edge.〈/p〉

Description: Bundesministerium für Bildung und Forschung http://dx.doi.org/10.13039/501100002347

Description: Horizon 2020 Framework Programme CONSTRAIN http://dx.doi.org/10.13039/100010661

Description: https://doi.org/10.5281/ZENODO.7051674

Description: https://doi.org/10.24381/cds.adbb2d47

Description: NORP-SORP Workshop on Polar Fresh Water: Sources, Pathways and Impacts of Freshwater in Northern and Southern Polar Oceans and Seas (SPICE-UP) What: Up to 60 participants at a time and more than twice as many registrants in total from 20 nations and across experience levels met to discuss the current status of research on freshwater in both polar regions, future directions, and synergies between the Arctic and Southern Ocean research communities When: 19–21 September 2022 Where: Online

Description: 〈jats:title〉Abstract〈/jats:title〉 〈jats:p〉—J. BLUNDEN, T. BOYER, AND E. BARTOW-GILLIES〈/jats:p〉 〈jats:p〉Earth’s global climate system is vast, complex, and intricately interrelated. Many areas are influenced by global-scale phenomena, including the “triple dip” La Niña conditions that prevailed in the eastern Pacific Ocean nearly continuously from mid-2020 through all of 2022; by regional phenomena such as the positive winter and summer North Atlantic Oscillation that impacted weather in parts the Northern Hemisphere and the negative Indian Ocean dipole that impacted weather in parts of the Southern Hemisphere; and by more localized systems such as high-pressure heat domes that caused extreme heat in different areas of the world. Underlying all these natural short-term variabilities are long-term climate trends due to continuous increases since the beginning of the Industrial Revolution in the atmospheric concentrations of Earth’s major greenhouse gases.〈/jats:p〉 〈jats:p〉In 2022, the annual global average carbon dioxide concentration in the atmosphere rose to 417.1±0.1 ppm, which is 50% greater than the pre-industrial level. Global mean tropospheric methane abundance was 165% higher than its pre-industrial level, and nitrous oxide was 24% higher. All three gases set new record-high atmospheric concentration levels in 2022.〈/jats:p〉 〈jats:p〉Sea-surface temperature patterns in the tropical Pacific characteristic of La Niña and attendant atmospheric patterns tend to mitigate atmospheric heat gain at the global scale, but the annual global surface temperature across land and oceans was still among the six highest in records dating as far back as the mid-1800s. It was the warmest La Niña year on record. Many areas observed record or near-record heat. Europe as a whole observed its second-warmest year on record, with sixteen individual countries observing record warmth at the national scale. Records were shattered across the continent during the summer months as heatwaves plagued the region. On 18 July, 104 stations in France broke their all-time records. One day later, England recorded a temperature of 40°C for the first time ever. China experienced its second-warmest year and warmest summer on record. In the Southern Hemisphere, the average temperature across New Zealand reached a record high for the second year in a row. While Australia’s annual temperature was slightly below the 1991–2020 average, Onslow Airport in Western Australia reached 50.7°C on 13 January, equaling Australia's highest temperature on record.〈/jats:p〉 〈jats:p〉While fewer in number and locations than record-high temperatures, record cold was also observed during the year. Southern Africa had its coldest August on record, with minimum temperatures as much as 5°C below normal over Angola, western Zambia, and northern Namibia. Cold outbreaks in the first half of December led to many record-low daily minimum temperature records in eastern Australia.〈/jats:p〉 〈jats:p〉The effects of rising temperatures and extreme heat were apparent across the Northern Hemisphere, where snow-cover extent by June 2022 was the third smallest in the 56-year record, and the seasonal duration of lake ice cover was the fourth shortest since 1980. More frequent and intense heatwaves contributed to the second-greatest average mass balance loss for Alpine glaciers around the world since the start of the record in 1970. Glaciers in the Swiss Alps lost a record 6% of their volume. In South America, the combination of drought and heat left many central Andean glaciers snow free by mid-summer in early 2022; glacial ice has a much lower albedo than snow, leading to accelerated heating of the glacier. Across the global cryosphere, permafrost temperatures continued to reach record highs at many high-latitude and mountain locations.〈/jats:p〉 〈jats:p〉In the high northern latitudes, the annual surface-air temperature across the Arctic was the fifth highest in the 123-year record. The seasonal Arctic minimum sea-ice extent, typically reached in September, was the 11th-smallest in the 43-year record; however, the amount of multiyear ice—ice that survives at least one summer melt season—remaining in the Arctic continued to decline. Since 2012, the Arctic has been nearly devoid of ice more than four years old.〈/jats:p〉 〈jats:p〉In Antarctica, an unusually large amount of snow and ice fell over the continent in 2022 due to several landfalling atmospheric rivers, which contributed to the highest annual surface mass balance, 15% to 16% above the 1991–2020 normal, since the start of two reanalyses records dating to 1980. It was the second-warmest year on record for all five of the long-term staffed weather stations on the Antarctic Peninsula. In East Antarctica, a heatwave event led to a new all-time record-high temperature of −9.4°C—44°C above the March average—on 18 March at Dome C. This was followed by the collapse of the critically unstable Conger Ice Shelf. More than 100 daily low sea-ice extent and sea-ice area records were set in 2022, including two new all-time annual record lows in net sea-ice extent and area in February.〈/jats:p〉 〈jats:p〉Across the world’s oceans, global mean sea level was record high for the 11th consecutive year, reaching 101.2 mm above the 1993 average when satellite altimetry measurements began, an increase of 3.3±0.7 over 2021. Globally-averaged ocean heat content was also record high in 2022, while the global sea-surface temperature was the sixth highest on record, equal with 2018. Approximately 58% of the ocean surface experienced at least one marine heatwave in 2022. In the Bay of Plenty, New Zealand’s longest continuous marine heatwave was recorded.〈/jats:p〉 〈jats:p〉A total of 85 named tropical storms were observed during the Northern and Southern Hemisphere storm seasons, close to the 1991–2020 average of 87. There were three Category 5 tropical cyclones across the globe—two in the western North Pacific and one in the North Atlantic. This was the fewest Category 5 storms globally since 2017. Globally, the accumulated cyclone energy was the lowest since reliable records began in 1981. Regardless, some storms caused massive damage. In the North Atlantic, Hurricane Fiona became the most intense and most destructive tropical or post-tropical cyclone in Atlantic Canada’s history, while major Hurricane Ian killed more than 100 people and became the third costliest disaster in the United States, causing damage estimated at $113 billion U.S. dollars. In the South Indian Ocean, Tropical Cyclone Batsirai dropped 2044 mm of rain at Commerson Crater in Réunion. The storm also impacted Madagascar, where 121 fatalities were reported.〈/jats:p〉 〈jats:p〉As is typical, some areas around the world were notably dry in 2022 and some were notably wet. In August, record high areas of land across the globe (6.2%) were experiencing extreme drought. Overall, 29% of land experienced moderate or worse categories of drought during the year. The largest drought footprint in the contiguous United States since 2012 (63%) was observed in late October. The record-breaking megadrought of central Chile continued in its 13th consecutive year, and 80-year record-low river levels in northern Argentina and Paraguay disrupted fluvial transport. In China, the Yangtze River reached record-low values. Much of equatorial eastern Africa had five consecutive below-normal rainy seasons by the end of 2022, with some areas receiving record-low precipitation totals for the year. This ongoing 2.5-year drought is the most extensive and persistent drought event in decades, and led to crop failure, millions of livestock deaths, water scarcity, and inflated prices for staple food items.〈/jats:p〉 〈jats:p〉In South Asia, Pakistan received around three times its normal volume of monsoon precipitation in August, with some regions receiving up to eight times their expected monthly totals. Resulting floods affected over 30 million people, caused over 1700 fatalities, led to major crop and property losses, and was recorded as one of the world’s costliest natural disasters of all time. Near Rio de Janeiro, Brazil, Petrópolis received 530 mm in 24 hours on 15 February, about 2.5 times the monthly February average, leading to the worst disaster in the city since 1931 with over 230 fatalities.〈/jats:p〉 〈jats:p〉On 14–15 January, the Hunga Tonga-Hunga Ha'apai submarine volcano in the South Pacific erupted multiple times. The injection of water into the atmosphere was unprecedented in both magnitude—far exceeding any previous values in the 17-year satellite record—and altitude as it penetrated into the mesosphere. The amount of water injected into the stratosphere is estimated to be 146±5 Terragrams, or ∼10% of the total amount in the stratosphere. It may take several years for the water plume to dissipate, and it is currently unknown whether this eruption will have any long-term climate effect.〈/jats:p〉