Search NASA⌕ Search

SEARCH · Search NASA

Results for “data parameterization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

Methods of testing parameterizations: Vertical ocean mixing

The ocean's velocity field is characterized by an exceptional variety of scales. While the small-scale oceanic turbulence responsible for the vertical mixing in the ocean is of scales a few centimeters and smaller, the oceanic general circulation is characterized by horizontal scales of thousands of kilometers. In oceanic general circulation models that are typically run today, the vertical structure of the ocean is represented by a few tens of discrete grid points. Such models cannot explicitly model the small-scale mixing processes, and must, therefore, find ways to parameterize them in terms of the larger-scale fields. Finding a parameterization that is both reliable and plausible to use in ocean models is not a simple task. Vertical mixing in the ocean is the combined result of many complex processes, and, in fact, mixing is one of the less known and less understood aspects of the oceanic circulation. In present models of the oceanic circulation, the many complex processes responsible for vertical mixing are often parameterized in an oversimplified manner. Yet, finding an adequate parameterization of vertical ocean mixing is crucial to the successful application of ocean models to climate studies. The results of general circulation models for quantities that are of particular interest to climate studies, such as the meridional heat flux carried by the ocean, are quite sensitive to the strength of the vertical mixing. We try to examine the difficulties in choosing an appropriate vertical mixing parameterization, and the methods that are available for validating different parameterizations by comparing model results to oceanographic data. First, some of the physical processes responsible for vertically mixing the ocean are briefly mentioned, and some possible approaches to the parameterization of these processes in oceanographic general circulation models are described in the following section. We then discuss the role of the vertical mixing in the physics of the large-scale ocean circulation, and examine methods of validating mixing parameterizations using large-scale ocean models.

Tziperman, Eli↗

Including auroral oval boundaries in the IRI model

At the last International Reference Ionosphere (IRI) Workshop it was decided that future editions of the model should include a representationof the aurroral oval boundaries. In this paper we review the different existing parameterizations of the auroral oval discussing their data base, boundary criteria, matematical formation, and overall usefulness for IRI. As a first candidate for incorporation into IRI we recomment the parameterization of the Feldstein ovals by Holzworth and Meng. Ways of implementing this model into IRI are discussed. We will also address adjustability with user-provided boundaries or boundary-related parameters, to better support strom-related studies.

Bilitza, A.↗

Understanding and Utilizing PBL Height Data from Multiple Observing Systems in the GEOS System

The accuracy of PBL height simulation is a key issue in many applications including forecasting near surface meteorology and air quality, however, it is a very challenging problem due to the lack of not only comprehensive, global Planetary Boundary Layer (PBL) observations but also a strategy and infrastructure to utilize PBL height data from a variety of sensors. Following the designation of PBL as an incubation class observable in the 2017 Decadal Survey, the PBL Incubation Study Team Report [14] made clear that “a future global PBL observing system requires modeling and data assimilation as essential components.” There is an urgent need for global modeling development in order to utilize Program of Record (POR) observations, assess their impacts, and identify gaps to be filled by future PBL missions. Our overall objective is to develop PBL data assimilation capabilities in the NASA Global Earth Observing System (GEOS), focusing on PBL height from multiple observing systems, to support the assessment and use of future PBL observations. The NASA GEOS system is composed of the GEOS global atmospheric general circulation model (AGCM) and the atmospheric data assimilation system (ADAS). The PBL parameterizations include the “Lock” K-profile scheme driven by surface and cloud-top buoyancy fluxes ([4]), and the “Louis” local scheme for stable conditions based on the Richardson number ([5]). Above the mixed layer defined by the Lock surface plume, shallow cumulus convection is represented by the mass flux scheme of [9]. Additional parameterizations are summarized in [1]. The ADAS employs the hybrid 4D Ensemble- Variational (EnVar) configuration ([15]), with the ensemble providing flow-dependent background error covariance information. The resultant analysis increments are fed back to the forecast model through the 4D incremental analysis update (IAU) approach ([11]). In this study, PBL height data are being or have been generated from radiosondes, GNSS RO, satellite (CATS, CALIPSO and ICESat-2) and ground-based (MPLNET) lidars, and wind profiler. Investigations have been conducted to specify quality marks for PBL height retrievals for the data assimilation purpose. These PBL height data have different strengths and weaknesses ([2], [3], [6], [7], [8], [10]), and the satellite PBL height data provide better global coverage and complement in-situ PBL height data. Radiosondes offer high accuracy and in situ measurement of temperature and humidity profiles, but with poor spatio-temporal sampling. The in-situ observing systems like MPLNET and wind profiler provide long history of PBL height records at each station. The GNSS RO based PBL height is retrieved based on the sharp gradients in refractivity profile that represent the fine vertical structure of temperature and moisture changes above the PBL. However, not all RO refractivity profiles reach the surface depending on location and regime, and RO refractivity retrievals can be negatively biased below 2km. The PBL height data from satellite lidars provide high resolution along track PBL height retrievals, but over land they are affected by previous day convective PBL aerosol and strongly associated with mixing layer and retrievals cannot be made below thick, attenuating clouds. A successful assimilation of PBL height data requires a thorough understanding of the observing method and the retrieval algorithm for each observing system in order to use the PBL height data from multiple observing systems properly. Due to the sensitivity of PBL height data to the observing method and choice of algorithm, it is important to use a model definition appropriate for each observation type to compute differences between PBL height data and model PBL height (OmFs). The GEOS model currently includes two PBL height definitions suitable for direct comparison with observed PBL height, and additional definitions are being added in this study. Evaluation of different model PBL height definitions is underway. Meanwhile, efforts have been made in the GEOS data assimilation system to develop PBL height data assimilation capability. PBL height data can be assimilated using two different approaches. The traditional approach is to construct an observation operator and its tangent linear and adjoint, which link control variables to PBL height data from each observing system. This observation operator can be very complicated, e.g., the lidar-based PBL height observation operator includes the backscatter lidar forward observation operator, the algorithm to derive PBL height from attenuated total backscatter, interpolation, and calculations handling the mismatch between observed and model scales. The other approach is to augment PBL height to the control variable vector, and it is adopted in this study. The latter approach was also used in previous studies, e.g., the assimilation of PBL height data from radiosonde and aircraft in the Real Time Mesoscale Analysis (RTMA) system for a dispersion modelling study ([13]); the PBL height assimilation study using lidar PBL height data at Greensburg, Kansas for a field campaign ([12]). The PBL height assimilation from multiple observing systems in this study allows us to take advantage of the diverse PBL height data that provide much better global coverage collectively under different meteorological conditions and with different temporal and spatial scales. As all the PBL heights are tightly coupled with the PBL thermodynamic variables, the strong correlations, which are provided by the 4D ensemble forecast, enable PBL height data from various sources to interact and combine coherently and provide additional information for PBL temperature and moisture fields. The results of comparisons among PBL height data from different sources and the evaluation of the model PBL height definitions with the PBL height data will be presented, and the PBL height data synergy strategies and preliminary results will also be discussed at the conference.

Y. Zhu↗

Standard ozone profiles from balloon and satellite data sets

Standard profiles based on upper level averaged profiles from BUV and lower level averaged profiles from balloon measurements are presented in a parametric representation as a function of time of year and latitude. The representation is a simple 4-parameter function representing the ozone amount (m-atm-cm) in each of 12 atmospheric layers defined following the standard Umkehr convention. The same parameterization is applied to the Nimbus-7 SBUV data and is compared to the BUV/balloon parameterization. The ozone variance unaccounted for by the representation is presented and discussed. The season-latitude representation reduces considerably the ozone variance at all levels and explains much of the correlation between layers. This simple representation and corresponding covariance matrix have been used as a priori information in the ozone vertical profile inversion of the Nimbus-7 SBUV experimental measurements.

Klenk, K. F.↗

Parameterized Cross Sections for Pion Production in Proton-Proton Collisions

An accurate knowledge of cross sections for pion production in proton-proton collisions finds wide application in particle physics, astrophysics, cosmic ray physics, and space radiation problems, especially in situations where an incident proton is transported through some medium and knowledge of the output particle spectrum is required when given the input spectrum. In these cases, accurate parameterizations of the cross sections are desired. In this paper much of the experimental data are reviewed and compared with a wide variety of different cross section parameterizations. Therefore, parameterizations of neutral and charged pion cross sections are provided that give a very accurate description of the experimental data. Lorentz invariant differential cross sections, spectral distributions, and total cross section parameterizations are presented.

Blattnig, Steve R.↗

An Improved Cryosat-2 Sea Ice Freeboard Retrieval Algorithm Through the Use of Waveform Fitting

We develop an empirical model capable of simulating the mean echo power cross product of CryoSat-2 SAR and SAR In mode waveforms over sea ice covered regions. The model simulations are used to show the importance of variations in the radar backscatter coefficient with incidence angle and surface roughness for the retrieval of surfaceelevation of both sea ice floes and leads. The numerical model is used to fit CryoSat-2 waveforms to enable retrieval of surface elevation through the use of look-up tables and a bounded trust region Newton least squares fitting approach. The use of a model to fit returns from sea ice regions offers advantages over currently used threshold retrackingmethods which are here shown to be sensitive to the combined effect of bandwidth limited range resolution and surface roughness variations. Laxon et al. (2013) have compared ice thickness results from CryoSat-2 and IceBridge, and found good agreement, however consistent assumptions about the snow depth and density of sea ice werenot used in the comparisons. To address this issue, we directly compare ice freeboard and thickness retrievals from the waveform fitting and threshold tracker methods of CryoSat-2 to Operation IceBridge data using a consistent set of parameterizations. For three IceBridge campaign periods from March 20112013, mean differences (CryoSat-2 IceBridge) of 0.144m and 1.351m are respectively found between the freeboard and thickness retrievals using a 50 sea ice floe threshold retracker, while mean differences of 0.019m and 0.182m are found when using the waveform fitting method. This suggests the waveform fitting technique is capable of better reconciling the seaice thickness data record from laser and radar altimetry data sets through the usage of consistent physical assumptions.

CryoSat-2↗

Parameterized desert/clear atmosphere limb-darkening model derived from earth radiation budget satellite along-track measurements

A parameterized desert/clear atmosphere limb-darkening model was derived using longwave measurements from the Earth Radiation Budget Satellite operating in a unique along-track mode that allows all points along the ground track to be viewed over an entire range of viewing zenith angles at essentially constant solar zenith angle. Application of the model to radiances for scenes defined as clear desert by the Earth Radiation Budget Experiment data analysis algorithms shows that this parameterized model reflects the geographical and diurnal behavior expected for the limb-darkening phenomenon.

Brooks, David R.↗

A Data Library of Liquid Clouds Modelled With a Large Eddy Simulation Framework

We describe a library of atmospheric large eddy simulations (LES) of liquid-phase boundary layer clouds constructed to enable aerosol–cloud–turbulence interaction studies, support parameterization evaluation and development, and provide training data for machine learning applications. The simulations use a modern LES framework designed for high numerical accuracy, coupled to a detailed spectral bin microphysical scheme. Case studies are configured to represent observed conditions in four key global cloud regions—the Northeastern Atlantic, Northeastern Pacific, Continental United States and Southern Ocean—following a semi-idealised approach. The library also includes aerosol concentration halving and doubling experiments to expose the sensitivities of the case studies to aerosol perturbations. Simulation results are compared to observations on a case-by-case basis, then the library's coverage is evaluated in terms of spreads in meteorological factors and atmospheric boundary layer attributes.

aerosol↗

Retrieval of Polar Stratospheric Cloud Microphysical Properties from Lidar Measurements: Dependence on Particle Shape Assumptions

A retrieval algorithm has been developed for the microphysical analysis of polar stratospheric cloud (PSC) optical data obtained using lidar instrumentation. The parameterization scheme of the PSC microphysical properties allows for coexistence of up to three different particle types with size-dependent shapes. The finite difference time domain (FDTD) method has been used to calculate optical properties of particles with maximum dimensions equal to or less than 2 mu m and with shapes that can be considered more representative of PSCs on the scale of individual crystals than the commonly assumed spheroids. Specifically. these are irregular and hexagonal crystals. Selection of the optical parameters that are input to the inversion algorithm is based on a potential data set such as that gathered by two of the lidars on board the NASA DC-8 during the Stratospheric Aerosol and Gas Experiment 0 p (SAGE) Ozone Loss Validation experiment (SOLVE) campaign in winter 1999/2000: the Airborne Raman Ozone and Temperature Lidar (AROTEL) and the NASA Langley Differential Absorption Lidar (DIAL). The 0 microphysical retrieval algorithm has been applied to study how particle shape assumptions affect the inversion of lidar data measured in leewave PSCs. The model simulations show that under the assumption of spheroidal particle shapes, PSC surface and volume density are systematically smaller than the FDTD-based values by, respectively, approximately 10-30% and approximately 5-23%.

Reichardt, J.↗

Guiding Principles for Geochemical/Thermodynamic Model Development and Validation in Nuclear Waste Disposal: A Close Examination of Recent Thermodynamic Models for H + —Nd 3+ —NO 3 - (—Oxalate) Systems

Development of a defensible source-term model (STM), usually a thermodynamical model for radionuclide solubility calculations, is critical to a performance assessment (PA) of a geologic repository for nuclear waste disposal. Such a model is generally subjected to rigorous regulatory scrutiny. In this article, we highlight key guiding principles for STM model development and validation in nuclear waste management. We illustrate these principles by closely examining three recently developed thermodynamic models with the Pitzer formulism for aqueous H + —Nd 3+ —NO 3 - (—oxalate) systems in a reverse alphabetical order of the authors: the XW model developed by Xiong and Wang, the OWC model developed by Oakes et al., and the GLC model developed by Guignot et al., among which the XW model deals with trace activity coefficients for Nd(III), while the OWC and GLC models are for concentrated Nd(NO 3 ) 3 electrolyte solutions. The principles highlighted include the following: (1) Principle 1. Validation against independent experimental data: A model should be validated against experimental data or field observations that have not been used in the original model parameterization. We tested the XW model against multiple independent experimental data sets including electromotive force (EMF), solubility, water vapor, and water activity measurements. The results show that the XW model is accurate and valid for its intended use for predicting trace activity coefficients and therefore Nd solubility in repository environments. (2) Principle 2. Testing for relevant and sensitive variables: Solution pH is such a variable for an STM and easily acquirable. All three models are checked for their ability to predict pH conditions in Nd(NO 3 ) 3 electrolyte solutions. The OWC model fails to provide a reasonable estimate for solution pH conditions, thus casting serious doubt on its validity for a source-term calculation. In contrast, both the XW and GLC models predict close-to-neutral pH values, in agreement with experimental measurements. (3) Principle 3. Honoring physical constraints: Upon close examination, it is found that the Nd(III)-NO 3 association schema in the OWC model suffers from two shortcomings. Firstly, its second stepwise stability constant for Nd(NO 3 ) 2+ (log K 2 ) is much higher than the first stepwise stability constant for NdNO 3 2+ (log K 1 ), thus violating the general rule of (log K 2 –log K 1 ) < 0, or $\frac{K1}{K2}$>1. Secondly, the OWC model predicts abnormally high activity coefficients for Nd(NO 3 ) 2 + (up to ~900) as the concentration increases. (4) Principle 4. Minimizing degrees of freedom for model fitting: The OWC model with nine fitted parameters is compared with the GLC model with five fitted parameters, as both models apply to the concentrated region for Nd(NO 3 ) 3 electrolyte solutions. The latter appears superior to the former because the latter can fit osmotic coefficient data equally well with fewer model parameters. The work presented here thus illustrates the salient points of geochemical model development, selection, and validation in nuclear waste management.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Estimation of new production in the ocean by compound remote sensing

Compound remote sensing, a technique for estimating oceanic new production using remotely sensed data on ocean color and temperature, is discussed. This approach, which depends on parameterizations developed from ship observations as well as on satellite data, yields more representative estimates of the large-scale average new production than those calculated from ship data alone. The approach is illustated with data for the Georges Bank.

Sathyendranath, Shubha↗

Observations of Middle Atmosphere CO from the UARS ISAMS during the Early Northern Winter 1991/92

Structure and kinematics of carbon monoxide in the upper stratosphere and lower mesosphere (10-0.03 hPa) are studied for the early northern winter 1991/92 using the Upper Atmosphere Research Satellite Improved Stratospheric and Mesospheric Sounder (ISAMS) measurements. The study is aided by data from a 6-week parameterized-chemistry run of the Goddard Space Flight Center 3D Chemistry and Transport Model (CTM), initialized on 8 December 1991. Generally, CO mixing ratios increase with height due to the increasing source contribution from CO, photolysis. In the tropical upper stratosphere. however, a local maximum in CO mixing ratio occurs. A simple photochemical model is used to show that this feature results largely from methane oxidation. In the extratropics the photochemical lifetime of CO is long, and therefore its evolution is dictated by large-scale motion of air. evidenced by strong correlation with Ertel potential vorticity. This makes CO one of the few useful observable tracers at the stratopause level and above. Thus CO maps are used to study the synoptic evolution of the polar vortex in early January 1992. Modified Lagrangian mean mixing diagnostics are applied to ISAMS and CTM data to examine the strength of the mixing barrier at the polar vortex edge. It is demonstrated that planetary wave activity weakens the barrier. promoting vortex erosion. The vortex erosion first appears in the lower mesosphere and subsequently descends through the upper stratosphere. and is attributed to effects of planetary wave dissipation. Agreement between ISAMS and CTM is good in the horizontal distribution of CO throughout the examined period, but vertical CO gradients in the CTM weaken with time relative to the ISAMS observations.

Allen, D. R.↗

Tail Lobe Revisited: Magnetic Field Modeling Based on Plasma Data

Plasma data from the ISEE-1 and -2 spacecraft during 1977-1980 have been used to determine the distribution of data points in the magnetotail in the range of distances -20 < XGSM < --15, i.e. which of the records that were located in the current sheet, in the tail lobe, in the magnetosheath and in the boundary layers respectively. The ISEE-1 and -2 magnetic field data for the records in the tail lobe were then used to model the tail lobe magnetic field dependence on the solar wind dynamic pressure, on the Interplanetary Magnetic Field (IMF) and on the Dst index. The tail lobe magnetic field was assumed to be dependent on the square root of the dynamic pressure based on the balance between the total magnetic pressure in the tail lobes and the dynamic pressure of the solar wind. The IMF dependent terms, added to the pressure term, were sought in many different forms while the Dst dependence of the tail lobe magnetic field was assumed to be linear. The field shows a strong dependence on the square root of the dynamic pressure and the different IMF dependent terms all constitute a significant contribution to the total field. However, the dependence on the Dst index turned out to be very weak at those down-tail distances. The results of this study are intended to be used for parameterizing future versions of the data-based models of the global magnetospheric magnetic field.

Karlsson, S. B. P.↗

Application of Spectral Analysis Techniques in the Intercomparison of Aerosol Data: Part III. Using Combined PCA to Compare Spatiotemporal Variability of MODIS, MISR and OMI Aerosol Optical Depth

Satellite measurements of global aerosol properties are very useful in constraining aerosol parameterization in climate models. The reliability of different data sets in representing global and regional aerosol variability becomes an essential question. In this study, we present the results of a comparison using combined principal component analysis (CPCA), applied to monthly mean, mapped (Level 3) aerosol optical depth (AOD) product from Moderate Resolution Imaging Spectroradiometer (MODIS), Multiangle Imaging Spectroradiometer (MISR), and Ozone Monitoring Instrument (OMI). This technique effectively finds the common space-time variability in the multiple data sets by decomposing the combined AOD field. The results suggest that all of the sensors capture the globally important aerosol regimes, including dust, biomass burning, pollution, and mixed aerosol types. Nonetheless, differences are also noted. Specifically, compared with MISR and OMI, MODIS variability is significantly higher over South America, India, and the Sahel. MODIS deep blue AOD has a lower seasonal variability in North Africa, accompanied by a decreasing trend that is not found in either MISR or OMI AOD data. The narrow swath of MISR results in an underestimation of dust variability over the Taklamakan Desert. The MISR AOD data also exhibit overall lower variability in South America and the Sahel. OMI does not capture the Russian wild fire in 2010 nor the phase shift in biomass burning over East South America compared to Central South America, likely due to cloud contamination and the OMI row anomaly. OMI also indicates a much stronger (boreal) winter peak in South Africa compared with MODIS and MISR.

satellite observation↗

Atmospheric Form Drag Coefficients Over Arctic Sea Ice Using Remotely Sensed Ice Topography Data, Spring 2009-2015

Sea ice topography significantly impacts turbulent energy/momentum exchange, e.g., atmospheric (wind) drag, over Arctic sea ice. Unfortunately, observational estimates of this contribution to atmospheric drag variability are spatially and temporally limited. Here we present new estimates of the neutral atmospheric form drag coefficient over Arctic sea ice in early spring, using high-resolution Airborne Topographic Mapper elevation data from NASA's Operation IceBridge mission. We utilize a new three-dimensional ice topography data set and combine this with an existing parameterization scheme linking surface feature height and spacing to form drag. To be consistent with previous studies investigating form drag, we compare these results with those produced using a new linear profiling topography data set. The form drag coefficient from surface feature variability shows lower values [less than 0.5-1 × 10(exp. −3)] in the Beaufort/Chukchi Seas, compared with higher values [greater than 0.5-1 ×10(exp. −3)] in the more deformed ice regimes of the Central Arctic (north of Greenland and the Canadian Archipelago), which increase with coastline proximity. The results show moderate interannual variability, including a strong increase in the form drag coefficient from 2013 to 2014/2015 north of the Canadian Archipelago. The form drag coefficient estimates are extrapolated across the Arctic with Advanced Scatterometer satellite radar backscatter data, further highlighting the regional/interannual drag coefficient variability. Finally, we combine the results with existing parameterizations of form drag from floe edges (a function of ice concentration) and skin drag to produce, to our knowledge, the first pan-Arctic estimates of the total neutral atmospheric drag coefficient (in early spring) from 2009 to 2015.

sea ice↗

Evaluation of a Reduced-Order Model for IBR Fault Response Representation via OEM Blackbox Models: Preprint

Driven by the need to capture the electromagnetic transients of transmission lines, inverter switching behavior, and detailed control systems, electromagnetic transient (EMT) studies have become increasingly important in industry, such as IBR interconnection study and fault study. However, original equipment manufacturer (OEM) inverter models typically include extensive parameters and proprietary settings that are unavailable to protection engineers. This paper introduces a data-driven, reduced-order model (ROM) developed as a PSCAD library component for use in EMT-based fault studies. The ROM replicates key OEM model behaviors without requiring detailed knowledge of control design or parameterization. The accompanying Python automation scripts streamline data generation, parameter fitting, and validation. The ROM's performance is demonstrated through comparison with both IEEE 2800-compliant and non-compliant OEM models in a real-world power system. Relay responses show nearly identical results, while simulation runtime is reduced by an average of 32.8\%, highlighting the ROM's practicality for protection engineers.

14 SOLAR ENERGY↗

CERES Fast Longwave And SHortwave Radiative Flux (FLASHFlux): Research to Operation

The Clouds and Earth’s Radiant Energy System (CERES) project provides satellite-based observations of how Earth’s energy flows are varying in time and space and how clouds and aerosols are affecting the Earth’s radiation budget. Nominally, CERES data products require months of validation and calibration before releasing a climate quality data. The Fast Longwave and Shortwave Radiative Flux (FLASHFlux) data product was developed as to provide key radiative flux data products within a week of observations for scientific, applied science research and educational usage. FLASHFlux achieves this by using simplified calibration process and optimized CERES data fusion production system, an operation meteorology product from Global Modeling and Assimilation Office(GMAO), and surface parameterized radiative flux models. The CERES FLASHFlux provides two data products: 1) an hourly Level 2 Single Scanner Footprint (SSF) data separately for Aqua and Terra observations, and 2) a daily Level 3 Time Interpolated and Spatially Averaged (TISA) 1ox1o gridded data that combines Aqua and Terra observations.The current operational Version 3C is being transitioned to Version4A. We compare the FLASHFlux Version4A fluxes using the specialized calibration to Version3C and to the climate quality CERES Edition4A TOA and surface fluxes over time for Terra and Aqua to determine the uncertainties of the calibration process. In addition, we also compare the impact of the fluxes by using the near real-time GMAO GEOS FP-IT instead of GMAO G5-CERES meteorology data.We also compare both SSF overpass products (for Terra and Aqua) and TISA radiative flux products to a set of surface measurement sites that are globally distributed and assess agreement.Lastly, we introduce the first CERES FF SSF data products from the CERES instrument on board NOAA-20.

P Sawaengphokhai↗