-
Emerging AI-based weather prediction models as downscaling tools
Authors:
Nikolay Koldunov,
Thomas Rackow,
Christian Lessig,
Sergey Danilov,
Suvarchal K. Cheedela,
Dmitry Sidorenko,
Irina Sandu,
Thomas Jung
Abstract:
The demand for high-resolution information on climate change is critical for accurate projections and decision-making. Presently, this need is addressed through high-resolution climate models or downscaling. High-resolution models are computationally demanding and creating ensemble simulations with them is typically prohibitively expensive. Downscaling methods are more affordable but are typically…
▽ More
The demand for high-resolution information on climate change is critical for accurate projections and decision-making. Presently, this need is addressed through high-resolution climate models or downscaling. High-resolution models are computationally demanding and creating ensemble simulations with them is typically prohibitively expensive. Downscaling methods are more affordable but are typically limited to small regions. This study proposes the use of existing AI-based numerical weather prediction systems (AI-NWP) to perform global downscaling of climate information from low-resolution climate models. Our results demonstrate that AI-NWP initalized from low-resolution initial conditions can develop detailed forecasts closely resembling the resolution of the training data using a one day lead time. We constructed year-long atmospheric fields using AI-NWP forecasts initialized from smoothed ERA5 and low-resolution CMIP6 models. Our analysis for 2-metre temperature indicates that AI-NWP can generate high-quality, long-term datasets and potentially perform bias correction, bringing climate model outputs closer to observed data. The study highlights the potential for off-the-shelf AI-NWP to enhance climate data downscaling, offering a simple and computationally efficient alternative to traditional downscaling techniques. The downscaled data can be used either directly for localized climate information or as boundary conditions for further dynamical downscaling.
△ Less
Submitted 25 June, 2024;
originally announced June 2024.
-
Implementation of implicit filter for spatial spectra extraction
Authors:
Kacper Nowak,
Sergey Danilov,
Vasco Müller,
Caili Liu
Abstract:
Scale analysis based on coarse-graining has been proposed recently as an alternative to Fourier analysis. It is now broadly used to analyze energy spectra and energy transfers in eddy-resolving ocean simulations. However, for data from unstructured-mesh models it requires interpolation to a regular grid. We present a high-performance Python implementation of an alternative coarse-graining method w…
▽ More
Scale analysis based on coarse-graining has been proposed recently as an alternative to Fourier analysis. It is now broadly used to analyze energy spectra and energy transfers in eddy-resolving ocean simulations. However, for data from unstructured-mesh models it requires interpolation to a regular grid. We present a high-performance Python implementation of an alternative coarse-graining method which relies on implicit filters using discrete Laplacians. This method can work on arbitrary (structured or unstructured) meshes and is applicable to the direct output of unstructured-mesh ocean circulation atmosphere models. The computation is split into two phases: preparation and solving. The first one is specific only to the mesh. This allows for auxiliary arrays that are then computed to be reused, significantly reducing the computation time. The second part consists of sparse matrix algebra and solving linear system. Our implementation is accelerated by GPUs to achieve unmatched performance and scalability. This results in processing data based on meshes with more than 10M surface vertices in a matter of seconds. As an illustration, the method is applied to compute spatial spectra of ocean currents from high-resolution FESOM2 simulations.
△ Less
Submitted 10 April, 2024;
originally announced April 2024.
-
Split-explicit external mode solver in finite volume sea ice ocean model FESOM2
Authors:
Tridib Banerjee,
Patrick Scholz,
Sergey Danilov,
Knut Klingbeil,
Dimitry Sidorenko
Abstract:
A novel split-explicit (SE) external mode solver for the Finite volumE Sea ice-Ocean Model (FESOM2) is presented. It is compared with the semi-implicit (SI) solver currently used in FESOM2. The SE solver utilises a dissipative asynchronous (forward-backward) time-stepping scheme. Its implementation with Arbitrary Lagrangian-Eulerian (ALE) vertical coordinates like Z-star and Z-tilde is explored. T…
▽ More
A novel split-explicit (SE) external mode solver for the Finite volumE Sea ice-Ocean Model (FESOM2) is presented. It is compared with the semi-implicit (SI) solver currently used in FESOM2. The SE solver utilises a dissipative asynchronous (forward-backward) time-stepping scheme. Its implementation with Arbitrary Lagrangian-Eulerian (ALE) vertical coordinates like Z-star and Z-tilde is explored. The comparisons are performed through multiple test cases involving idealised and realistic global simulations. The SE solver demonstrates lower phase errors and dissipation, but maintain a simulated mean ocean state very similar to the SI solver. The SE solver is also shown to possess better run-time performance and parallel scalability across all tested workloads.
△ Less
Submitted 14 November, 2023;
originally announced November 2023.
-
Discrete variance decay analysis of spurious mixing
Authors:
Tridib Banerjee,
Sergey Danilov,
Knut Klingbeil
Abstract:
Expressions for local discrete variance decay (DVD) rates are directly derived from discrete tracer equations without any assumptions on discrete fluxes of the second moment. Spurious mixing (SM) associated with numerical implementations of scalar advection and diffusion is thus estimated. The new framework is shown to avoid the need for second-moment flux definition when solved on finite-volume c…
▽ More
Expressions for local discrete variance decay (DVD) rates are directly derived from discrete tracer equations without any assumptions on discrete fluxes of the second moment. Spurious mixing (SM) associated with numerical implementations of scalar advection and diffusion is thus estimated. The new framework is shown to avoid the need for second-moment flux definition when solved on finite-volume cell edges but still invoke certain second-moment fluxes when the DVD rates are partitioned to participating cell nodes. These implied discrete fluxes are shown to differ from those proposed in earlier literature (but share the same dissipative part) and thus reveal the non-uniqueness of their nature. They are shown to be ambiguous for high-order advection schemes introducing uncertainty to the locality of any estimates produced by a DVD approach. Additional damping of flux divergence through temporal averaging or some coarse-graining is thus shown to be necessary. Through the application of this technique, SM is found to be correlated with the distribution of eddy kinetic energy. The contribution from vertical advection to SM is found to be relatively small and correlated with the distribution of buoyancy fluxes. The explored high-order schemes are found to demonstrate levels of spurious mixing which may locally exceed background physical mixing.
△ Less
Submitted 25 October, 2023;
originally announced October 2023.
-
Simulating linear kinematic features in viscous-plastic sea ice models on quadrilateral and triangular grids
Authors:
C. Mehlmann,
S. Danilov,
M. Losch,
J. F. Lemieux,
N. Hutter,
T. Richter,
P. Blain,
E. C. Hunke,
P Korn
Abstract:
Linear Kinematic Features (LKFs) are found everywhere in the Arctic sea ice cover. They are strongly localized deformations often associated with the formation of leads and pressure ridges. Viscous-plastic sea ice models start to generate LKFs at high spatial grid resolution, typically with a grid spacing below 5 km. Besides grid spacing, other aspects of a numerical implementation, such as discre…
▽ More
Linear Kinematic Features (LKFs) are found everywhere in the Arctic sea ice cover. They are strongly localized deformations often associated with the formation of leads and pressure ridges. Viscous-plastic sea ice models start to generate LKFs at high spatial grid resolution, typically with a grid spacing below 5 km. Besides grid spacing, other aspects of a numerical implementation, such as discretization details, may affect the number and definition of simulated LKFs. To explore these effects, the solutions of sea ice models with different grid spacings, mesh types, and numerical discretization techniques are compared in an idealized configuration, which could also serve as a benchmark problem in the future. The A, B, and C-grid discretizations of sea ice dynamics on quadrilateral meshes leads to a similar number of LKFs as the A-grid approximation on triangular meshes (with the same number of vertices). The discretization on an Arakawa CD-grid on both structured quadrilateral and triangular meshes resolves the same number of LKFs as conventional Arakawa A-grid, B-grid, and C-grid approaches, but on two times coarser meshes. This is due to the fact that the CD-grid approach has a higher number of degrees of freedom to discretize the velocity field. Due to its enhanced resolving properties, the CD-grid discretization is an attractive alternative to conventional discretizations.
△ Less
Submitted 7 March, 2021;
originally announced March 2021.
-
Observation of anomalously strong penetration of terahertz electric field through terahertz-opaque gold films into a GaAs/AlGaAs quantum well
Authors:
S. D. Ganichev,
S. N. Danilov,
M. Kronseder,
D. Schuh,
I. Gronwald,
D. Bougeard,
E. L. Ivchenko,
A. Ya. Shul'man
Abstract:
We observe an anomalously high electric field of terahertz (THz) radiation acting on a two-dimensional electron gas (2DEG) placed beneath a thin gold film, which, however, is supposed to be opaque at THz frequencies. We show that the anomalously strong penetration of the THz electric field through a very high conductive gold film emerges if two conditions are fulfilled simultaneously: (i) the film…
▽ More
We observe an anomalously high electric field of terahertz (THz) radiation acting on a two-dimensional electron gas (2DEG) placed beneath a thin gold film, which, however, is supposed to be opaque at THz frequencies. We show that the anomalously strong penetration of the THz electric field through a very high conductive gold film emerges if two conditions are fulfilled simultaneously: (i) the film's thickness is less than the skin depth and (ii) the THz electric field is measured beneath the film at distances substantially smaller than the radiation wavelength. We demonstrate that under these conditions the strength of the field acting on a 2DEG is almost the same as it would be in the absence of the gold film. The effect is detected for macroscopically homogeneous perforation-free gold films illuminated by THz-laser radiation with a spot smaller than the film area. This eliminates the near-field of the edge diffraction as a possible cause of the anomalous penetration. The microscopic origin of the effect remains unexplained in its details, yet. The observed effect can be used for the development of THz devices based on two-dimensional materials requiring robust highly conducting top gates placed at less than nanometer distance from the electron gas location.
△ Less
Submitted 3 June, 2020; v1 submitted 16 February, 2020;
originally announced February 2020.
-
Time resolution and dynamic range of field effect transistor based terahertz detectors
Authors:
Przemyslaw Zagrajek,
Sergey N. Danilov,
Jacek Marczewski,
Michal Zaborowski,
Cezary Kolacinski,
Dariusz Obrebski,
Pawel Kopyt,
Bartlomiej Salski,
Dmytro But,
Wojciech Knap,
Sergey D. Ganichev
Abstract:
We studied time resolution and response power dependence of three terahertz detectors based on significantly different types of field effect transistors. We analyzed the photoresponse of custom-made Si junctionless FETs, Si MOSFETs and GaAs-based high electron mobility transistors detectors. Applying monochromatic radiation of high power, pulsed, line-tunable molecular THz laser, which operated at…
▽ More
We studied time resolution and response power dependence of three terahertz detectors based on significantly different types of field effect transistors. We analyzed the photoresponse of custom-made Si junctionless FETs, Si MOSFETs and GaAs-based high electron mobility transistors detectors. Applying monochromatic radiation of high power, pulsed, line-tunable molecular THz laser, which operated at frequencies in the range from 0.6-3.3 THz, we demonstrated that all these detectors have at least nanosecond response time. We showed that detectors yield a linear response in a wide range of radiation power. At high powers the response saturates varying with radiation power P as $U = R_0 P/(1+P/P_s)$, where $R_0$ is the low power responsivity, $P_s$ is the saturation power. We demonstrated that the linear part response decreases with radiation frequency increase as $R_0 \propto f^{-3}$, whereas the power at which signal saturates increases as $P_s \propto f^3$. We discussed the observed dependences in the framework of the Dyakonov-Shur mechanism and detector-antenna impedance matching. Our study showed that FET transistors can be used as ultrafast room temperature detectors of THz radiation and that their dynamic range extends over many orders of magnitude of power of incoming THz radiation. Therefore, when embedded with current driven read out electronics they are very well adopted for operation with high power pulsed sources.
△ Less
Submitted 21 May, 2019;
originally announced May 2019.
-
Infrared/Terahertz Spectra of the Photogalvanic Effect in (Bi,Sb)Te based Three Dimensional Topological Insulators
Authors:
H. Plank,
J. Pernul,
S. Gebert,
S. N. Danilov,
J. König-Otto,
S. Winnerl,
M. Lanius,
J. Kampmeier,
G. Mussler,
I. Aguilera,
D. Grützmacher,
S. D. Ganichev
Abstract:
We report on the systematic study of infrared/terahertz spectra of photocurrents in (Bi,Sb)Te based three dimensional topological insulators. We demonstrate that in a wide range of frequencies, ranging from fractions up to tens of terahertz, the photocurrent is caused by the linear photogalvanic effect (LPGE) excited in the surface states. The photocurrent spectra reveal that at low frequencies th…
▽ More
We report on the systematic study of infrared/terahertz spectra of photocurrents in (Bi,Sb)Te based three dimensional topological insulators. We demonstrate that in a wide range of frequencies, ranging from fractions up to tens of terahertz, the photocurrent is caused by the linear photogalvanic effect (LPGE) excited in the surface states. The photocurrent spectra reveal that at low frequencies the LPGE emerges due to free carrier Drude-like absorption. The spectra allow to determine the room temperature carrier mobilities in the surface states despite the presents of thermally activate residual impurities in the material bulk. In a number of samples we observed an enhancement of the linear photogalvanic effect at frequencies between 30÷60 THz, which is attributed to the excitation of electrons from helical surface to bulk conduction band states. Under this condition and applying oblique incidence we also observed the circular photogalvanic effect driven by the radiation helicity.
△ Less
Submitted 30 November, 2017;
originally announced November 2017.