Search

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.

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., Comment: 16 pages, 9 figures

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., Comment: Submitted to Geoscientific Model Development number: GMD-2023-208. Development code for FESOM2.5 with explicit subcycling available at: https://zenodo.org/doi/10.5281/zenodo.10040943. Name of funder: Deutsche Forschungsgemeinschaft. Grant agreement or award number: 274762653

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., Comment: Submitted to Ocean Modelling Manuscript number: OCEMOD-D-23-00145. Name of funder: Deutsche Forschungsgemeinschaft. Grant agreement or award number: 274762653

We report on the observation of complex nonlinear intensity dependence of the circular and linear photogalvanic currents induced by infrared radiation in compensated (Bi$_{0.3}$Sb$_{0.7}$)$_2$(Te$_{0.1}$Se$_{0.9}$)$_3$ 3D topological insulators. The photocurrents are induced by direct optical transitions between topological surface and bulk states. We show that an increase of the radiation intensity results first in a highly superlinear raise of the amplitude of both types of photocurrents, whereas at higher intensities the photocurrent saturates. Our analysis of the observed nonlinearities shows that the superlinear behavior of the photocurrents is caused by a heating of the electron gas, while the saturation is induced by a slow relaxation of the photoexcited carriers resulting in absorbance bleaching. The observed nonlinearities give access to the Fermi level position with respect to the Dirac point and the energy relaxation times of Dirac fermions providing an experimental room temperature probe for topological surface states., Comment: 12 pages, 7 figures

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.

We demonstrate that the terahertz/infrared radiation induced photogalvanic effect, which is sensitive to the surface symmetry and scattering details, can be applied to study the high frequency conductivity of the surface states in (Bi1-xSbx)2Te3 based three dimensional (3D) topological insulators (TI). In particular, measuring the polarization dependence of the photogalvanic current and scanning with a micrometre sized beam spot across the sample, provides access to (i) topographical inhomogeneity's in the electronic properties of the surface states and (ii) the local domain orientation. An important advantage of the proposed method is that it can be applied to study TIs at room temperature and even in materials with a high electron density of bulk carriers., Comment: 6 pages, 4 figures

We report on the observation of the reststrahl band assisted photocurrents in epitaxial graphene on SiC excited by infrared radiation. The peculiar spectral dependence for frequencies lying within the reststrahl band of the SiC substrate provides a direct and noninvasive way to probe the electric field magnitude at atomic distances from the material's surface. Furthermore our results reveal that nonlinear optical and optoelectronic phenomena in 2D crystals and other atomic scale structures can be giantly enhanced by a proper combination of the spectral range and substrate material.

Far-infrared cyclotron resonance photoconductivity (CRP) is investigated in HgTe quantum wells (QWs) of various widths grown on (013) oriented GaAs substrates. It is shown that CRP is caused by the heating of two-dimensional electron gas (2DEG). From the resonance magnetic field strength effective masses and their dependence on the carrier concentration is obtained. We found that the effective mass in each sample slightly increases from the value (0.0260 \pm 0.0005)m_0 at N_s = 2.2x10^11 cm^(-2) to (0.0335 \pm 0.0005)m_0 at N_s = 9.6x10^11 cm^(-2). Compared to determination of effective masses by the temperature dependence of magnitudes of the Shubnikov-de Haas (SdH) oscillations used so far in this material our measurements demonstrate that the CRP provides a more accurate (about few percents) tool. Combining optical methods with transport measurements we found that the transport time substantially exceeds the cyclotron resonance lifetime as well as the quantum lifetime which is the shortest., Comment: 3 pages, 2 figures

Spin-orbit coupling provides a versatile tool to generate and to manipulate the spin degree of freedom in low-dimensional semiconductor structures. The spin Hall effect, where an electrical current drives a transverse spin current and causes a nonequilibrium spin accumulation observed near the sample boundary, the spin-galvanic effect, where a nonequilibrium spin polarization drives an electric current, or the reverse process, in which an electrical current generates a nonequilibrium spin polarization, are all consequences of spin-orbit coupling. In order to observe a spin Hall effect a bias driven current is an essential prerequisite. The spin separation is caused via spin-orbit coupling either by Mott scattering (extrinsic spin Hall effect) or by Rashba or Dresselhaus spin splitting of the band structure (intrinsic spin Hall effect). Here we provide evidence for an elementary effect causing spin separation which is fundamentally different from that of the spin Hall effect. In contrast to the spin Hall effect it does not require an electric current to flow: It is spin separation achieved by spin-dependent scattering of electrons in media with suitable symmetry. We show that by free carrier (Drude) absorption of terahertz radiation spin separation is achieved in a wide range of temperatures from liquid helium up to room temperature. Moreover the experimental results give evidence that simple electron gas heating by any means is already sufficient to yield spin separation due to spin-dependent energy relaxation processes of nonequilibrium carriers., Comment: 19 pages, 4 figures, 1 table