Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “NATURAL CONVECTION”

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 811 records · Page 45

The use of satellite data in understanding and predicting convective and large-scale dynamical processes

A two-layer truncated baroclinic spectral model was developed to study the long-term evolution of disturbances to a baroclinically unstable mean flow. Topography and crudely-parameterized radiative processes were accounted for. As a result of Robert Schlaak's discovery of the underlying barotropic nature of the index oscillation as well as reviewers suggestions about the original manuscript, the model has been revised to allow for barotropic as well as baroclinic wave-mean flow interactions. The form-drag exerted by the topography on the barotropic part of the mean flow is larger than on the baroclinic part and thus researchers anticipate significant changes from the original calculations on the index oscillation when it is strongly modulated by topography. Researchers believe that since the index oscillation accounts for a significant portion of atmospheric temporal variance, the long term predictability could be improved if reliable forecasts of the index oscillation were available. Two spectral models of the index oscillation, one barotropic and the other baroclinic, have been developed. The latter allows for moisture, radiation, land-sea temperature countrasts, and energy exchanges with the underlying surface.

Dutton, J. A.↗

On the contribution of the thermospheric neutral wind to high-latitude energetics

The significance of the neutral wind dynamo as a contributor to ionospheric energetics at high latitudes is addressed by determining separately the electrical energy flux from the magnetospheric dynamo and the maximum electrical energy flux available due to the neutral wind dynamo, neglecting the natural feedback processes between the two. The neutral wind dynamo is found to contribute significantly to the flux of energy exchanged between the magnetosphere and thermosphere, particularly in the central polar cap. In the region of the magnetospheric convection reversal, the amount of available energy flux from the neutral wind dynamo can exceed that provided by the magnetospheric dynamo making the neutrals a dominant contributor to local electrodynamics.

Thayer, J. P.↗

A Cartesian, cell-based approach for adaptively-refined solutions of the Euler and Navier-Stokes equations

A Cartesian, cell-based approach for adaptively-refined solutions of the Euler and Navier-Stokes equations in two dimensions is developed and tested. Grids about geometrically complicated bodies are generated automatically, by recursive subdivision of a single Cartesian cell encompassing the entire flow domain. Where the resulting cells intersect bodies, N-sided 'cut' cells are created using polygon-clipping algorithms. The grid is stored in a binary-tree structure which provides a natural means of obtaining cell-to-cell connectivity and of carrying out solution-adaptive mesh refinement. The Euler and Navier-Stokes equations are solved on the resulting grids using a finite-volume formulation. The convective terms are upwinded: a gradient-limited, linear reconstruction of the primitive variables is performed, providing input states to an approximate Riemann solver for computing the fluxes between neighboring cells. The more robust of a series of viscous flux functions is used to provide the viscous fluxes at the cell interfaces. Adaptively-refined solutions of the Navier-Stokes equations using the Cartesian, cell-based approach are obtained and compared to theory, experiment, and other accepted computational results for a series of low and moderate Reynolds number flows.

Coirier, William J.↗

A Cartesian, cell-based approach for adaptively-refined solutions of the Euler and Navier-Stokes equations

A Cartesian, cell-based approach for adaptively-refined solutions of the Euler and Navier-Stokes equations in two dimensions is developed and tested. Grids about geometrically complicated bodies are generated automatically, by recursive subdivision of a single Cartesian cell encompassing the entire flow domain. Where the resulting cells intersect bodies, N-sided 'cut' cells are created using polygon-clipping algorithms. The grid is stored in a binary-tree data structure which provides a natural means of obtaining cell-to-cell connectivity and of carrying out solution-adaptive mesh refinement. The Euler and Navier-Stokes equations are solved on the resulting grids using a finite-volume formulation. The convective terms are upwinded: A gradient-limited, linear reconstruction of the primitive variables is performed, providing input states to an approximate Riemann solver for computing the fluxes between neighboring cells. The more robust of a series of viscous flux functions is used to provide the viscous fluxes at the cell interfaces. Adaptively-refined solutions of the Navier-Stokes equations using the Cartesian, cell-based approach are obtained and compared to theory, experiment and other accepted computational results for a series of low and moderate Reynolds number flows.

Coirier, William J.↗

Climatology and natural variability of the global hydrologic cycle in the GLA atmospheric general circulation model

Time average climatology and low-frequency variabilities of the global hydrologic cycle (GHC) in the Goddard Laboratory for Atmospheres (GLA) general circulation model (GCM) were investigated in the present work. A 730-day experiment was conducted with the GLA GCM forced by insolation, sea surface temperature, and ice-snow undergoing climatological annual cycles. Ifluences of interactive soil moisture on time average climatology and natural variability of the GHC were also investigated by conducting 365-day experiments with and without interactive soil moisture. Insolation, sea surface temperature, and ice-snow were fixed at their July levels in the latter two experiments. Results show that the model's time average hydrologic cycle variables for July in all three experiments agree reasonably well with observations. Except in the case of precipitable water, the zonal average climates of the annual cycle experiment and the two perpetual July experiments are alike, i.e., their differences are within limits of the natural variability of the model's climate. Statistics of various components of the GHC, i.e., water vapor, evaporation, and precipitation, are significantly affected by the presence of interactive soil moisture. A long-term trend is found in the principal empirical modes of variability of ground wetness, evaporation, and sensible heat. Dominant modes of variability of these quantities over land are physically consistent with one another and with land surface energy balance requirements. The dominant mode of precipitation variability is found to be closely related to organized convection over the tropical western Pacific Ocean. The precipitation variability has timescales in the range of 2 to 3 months and can be identified with the stationary component of the Madden-Julian Oscillation. The precipitation mode is not sensitive to the presence of interactive soil moisture but is closely linked to both the rotational and divergent components of atmospheric moisture transport. The present results indicate that globally coherent natural variability of the GHC in the GLA GCM has two basic timescales in the absence of annual cycles of external forcings: a long-term trend associated with atmosphere-soil moisture interaction which affects the model atmosphere mostly over midlatitude continental regions and a large-scale 2- to 3-month variability associated with atmospheric moist processes over the western Pacific Ocean.

Lau, K.-M.↗

Thermal and Reliability Characterization of an Epoxy Resin-Based Double-Side Cooled Power Module

Wide-Band Gap (WBG) power devices have become a promising option for high-power applications due to the superior material properties over traditional Silicon. To not limit WBG devices’ mother nature, a rugged and high-performance power device packaging solution is necessary. This study proposes a Double-Side Cooled (DSC) 1.2 kV half-bridge power module having dual epoxy resin insulated metal substrate (eIMS) for solving convectional power module challenges and providing a cost-effective solution. The thermal performance outperforms traditional Alumina (Al2O3) Direct Bonded Copper (DBC) DSC power module due to moderate thermal conductivity (10 W/mK) and thin (120 mm) epoxy resin composite dielectric working as the IMS insulation layer. This novel organic dielectric can withstand high voltage (5 kVAC @ 120 μm) and has a Glass Transition Temperature (Tg) of 300°C, which is suitable for high-power applications. In the thermal-mechanical modeling, the organic DSC power module can pass the thermal cycling test over 1,000 cycles by optimizing the mechanical properties of the encapsulant material. In conclusion, this article not only proposes a competitive organic-based power module but also a methodology of evaluation for thermal and mechanical performance.

42 ENGINEERING↗

The chemical composition and climatology of the earth's early atmosphere

The earth's climate as it relates to the evolution of life is discussed.. Seven fundamental characteristics of the early evolutionary environment are examined, including a carbon dioxide and water vapor atmosphere, atmospheric mass between 500 and 1000 mb, a global hydrosphere, lowered solar luminosity, hospitable average global temperatures, a convectively active atmosphere, and trace gases. The influence of the early earth's extensive hydrosphere on the origin of life is considered. The warming of that hydrosphere due to radiative fluxes and the greenhouse effect is examined, and the nature of the feedback between clouds and climate is addressed.

Henderson-Sellers, A.↗

Nighttime observations of thunderstorm electrical activity from a high altitude airplane

Nocturnal thunderstorms were observed from above and features of cloud structure and lightning which are not generally visible from the ground are discussed. Most, lightning activity seems to be associated with clouds with strong convective cauliflower tops. In both of the storms lightning channels were visible in the clear air above the cloud. It is shown that substances produced by thunderstorm electrical discharges can be introduced directly into the stratosphere. The cause and nature of the discharges above the cloud are not clear. They may be produced by accumulations of space charge in the clear air above the cloud. The discharges may arise solely because of the intense electric fields produced by charges within the cloud. In the latter case the ions introduced by these discharges will increase the electrical conductivity of the air above the cloud and increase the conduction current that flows from the cloud to the electrosphere. More quantitative data at higher resolution may show significant spectral differences between cloud to ground and intracloud strokes. It is shown that electric field change data taken with an electric field change meter mounted in an airplane provide data on lightning discharges from above that are quite similar to those obtained from the ground in the past. The optical signals from dart leaders, from return strokes, and from continuing currents are recognizable, can be used to provide information on the fine structure of lightning, and can be used to distinguish between cloud to ground and intracloud flashes.

Brook, M.↗

Retarding potential analyzer research

Studies were continued in the areas of ionospheric electrodynamics, plasma irregularities, and ion-neutral coupling. The spectral characteristics of structure in the ion concentration and the electric field at high latitudes were examined. It was found that in the region of velocity shears the relationship between the slopes of the two spectra are dependent on the existence of large field-aligned currents and apparently large ionospheric conductivity gradients. The relationship between the Joule dissipation and the locally measured value of the ion temperature were studied by examining the horizontal ion drift velocity. The importance of ion-neutral coupling was investigated in additional studies involving the nature of the coupling at high latitude boundaries and during times of magnetic disturbances. Much of this behavior can be modeled using globular circulation models. In order to study the storm time response of the thermosphere the global ion convection pattern was derived from DE-2 measurements and used as input in a circulation model. The calculated variables of neutral velocity and composition and temperature were compared with measurements to validate the technique.

Hanson, W. B.↗

Electrochemical growth of linear conducting crystals in microgravity

Much attention has been given to the synthesis of linear conducting materials. These inorganic, organic, and polymeric materials have some very interesting electrical and optical properties, including low temperature superconductivity. Because of the anisotropic nature of these compounds, impurities and defects strongly influences the unique physical properties of such crystals. Investigations have demonstrated that electrochemical growth has provided the most reproducible and purest crystals. Space, specifically microgravity, eliminates phenomena such as buoyancy driven convection, and could permit formation of crystals many times purer than the ones grown to date. Several different linear conductors were flown on Get Away Special G-007 on board the Space Shuttle Columbia, STS 61-C, the first of a series of Project Explorer payloads. These compounds were grown by electrochemical methods, and the growth was monitored by photographs taken throughout the mission. Due to some thermal problems, no crystals of appreciable size were grown. The experimental results will be incorporated into improvements for the next 2 missions of Project Explorer. The results and conclusions of the first mission are discussed.

Cronise, Raymond J., IV↗

On the drag of model dendrite fragments at low Reynolds number

An experimental study of low Reynolds number drag on laboratory models of dendrite fragments has been conducted. The terminal velocities of the dendrites undergoing free fall along their axis of symmetry were measured in a large Stokes flow facility. Corrections for wall interference give nearly linear drag vs Reynolds number curves. Corrections for both wall interference and inertia effects show that the dendrite Stokes settling velocities are always less than that of a sphere of equal mass and volume. In the Stokes limit, the settling speed ratio is found to correlate well with the primary dendrite arm aspect ratio and a second dimensionless shape parameter which serves as a measure of the fractal-like nature of the dendrite models. These results can be used to estimate equiaxed grain velocities and distance of travel in metal castings. The drag measurements may be used in numerical codes to calculate the movement of grains in a convecting melt in an effort to determine macrosegregation patterns caused by the sink/float mechanism.

Zakhem, R.↗

On the drag of model dendrite fragments at low Reynolds number

An experimental study of low Reynolds number drag on laboratory models of dendrite fragments has been conducted. The terminal velocities of the dendrites undergoing free fall along their axis of symmetry were measured in a large Stokes flow facility. Corrections for wall interference give nearly linear drag vs Reynolds number curves. Corrections for both wall interference and inertia effects show that the dendrite Stokes settling velocities are always less than that of a sphere of equal mass and volume. In the Stokes limit, the settling speed ratio is found to correlate well with primary dendrite arm aspect ratio and a second dimensionless shape paremeter which serves as a measure of the fractal-like nature of the dendrite models. These results can be used to estimate equiaxed grain velocities and distance of travel in metal castings. The drag measurements may be used in numerical codes to calculate the movement of grains in a convecting melt in an effort to determine macrosegregation patterns caused by the sink/float mechanism.

Zakhem, R.↗

Solar Rossby Wave 'Hills' Identified As Supergranules

We explore the nature of 'hills' observed on the solar surface which had previously been attributed to Rossby waves. We investigate the sol ar hills phenomenon by analyzing the output from a synthetic model ba sed solely on the observed solar photospheric convection spectrum. We show that the characteristics of these hills can be explained by the corrugation of the surface produced by the radial flows of the conve ction. The hills in our simulations are dominated by supergranules, a well-known component of solar convection. Rossby waves have been predicted to exist within the Sun and may play an important role in the d ynamics of the solar interior, including the Sun's differential rotat ion and magnetic dynamo. Our study suggests, however, that the hills observed at the solar limb do not confirm the existence of solar Ross by waves.

Williams, P. E.↗

Surface Biology and Geology (SBG) Thermal Infrared (TIR) Free -Flyer Concept

Scientists and engineers from the United States and Italy are developing a thermal infrared (TIR) free-flyer concept that forms an essential element of the NASA Surface Biology & Geology (SBG) designated observable and the greater NASA Earth System Observatory (ESO), which also includes observations of aerosols, clouds, convection, and precipitation; mass change; and surface deformation & change.In response to the recommendations from the 2017 National Academies decadal survey for Earth science, NASA initiated the SBG designated observable with five key research and applications focus areas: ecosystems and natural resources, hydrology, weather, climate, and solid Earth. SBG includes spaceborne measurements of hyperspectral imagery in the visible to shortwave infrared (0.4-2.5 um) and multispectral imagery in the mid and thermal infrared (3-12 um) that provide the remote sensing data needed to inform each research and applications area. High-level TIR data products include Earth surface temperature & emissivity, evapotranspiration, substrate composition, volcanic plumes, and high-temperature features. While the science, applications, and technology build on prior mission studies conducted for more than a decade, it has only been recently that a joint team of scientists and engineers from the NASA Jet Propulsion Laboratory (JPL), Agenzia Spaziale Italiana (ASI), Istituto Nazionale Geofisica e Volcanologia (INGV), and the Istituto Nazionale Astrofisica (INAF) have developed a compelling concept for a TIR free-flyer. In this concept, the TIR instrument is an eight-band radiometer. Seven of the bands are between 3 and 12 um with an additional band at 1.65 um. The instrument’s 68.8 deg total viewing angle and the free-flyer’s 665 km operational altitude results in a ground sampling distance (GSD) of <60 m at nadir, produces a swath width of 935 km, and allows for a 2-3-day revisit time. A two-band visible and near infrared (VNIR) camera, with <30 m GSD at nadir and 935 km swath width, will complement the TIR instrument. VNIR high-level data products include top & bottom of atmosphere reflectance and normalized difference vegetation indices (NDVI). The TIR instrument and VNIR camera will be mounted on a PLATiNO+ spacecraft, a multi-application, reconfigurable, and scalable platform. The launcher selected for the TIR free-flyer is the VEGA-C. To maximize the science and applications benefits the SBG-TIR free-flyer team are collaborating with personnel from the ESA Land Surface Temperature Monitoring mission (LSTM) as well individuals from the Thermal Infrared Imaging Satellite for High-resolution Natural resource Assessment (TRISHNA), a joint mission by CNES and ISRO. Pre-Phase A concept studies for the SBG-TIR free-flyer commenced in the spring 2021. This allows for the concept to be refined; for technical risks to be identified and mitigated, as necessary; and strategies for reducing the development schedule to be defined and assessed before conducting a NASA-ASI Mission Concept Review (MCR) in the spring of 2022 and a NASA Key Decision Point-A (KDP-A) as early as the summer of 2022. These are key milestones in being able to launch and begin operating the SBG-TIR free-flyer in the second half of the decade, so that many of the questions posed by the National Academies in the 2017 decadal survey can be answered.

Fabrizia-Buongiorno, Maria↗

Radiative Heat Transport and Optical Characterization of High Temperature Molten Salts

The goal of this project was to advance the field of molten fluoride radiative heat transfer (RHT), through optical property measurements, heat transfer loop construction and operation, and computational fluid dynamics simulations. The following outcomes have been accomplished. A comprehensive review has been conducted to compile all available RHT-relevant optical absorption data for molten fluorides. The various measurement methods have also been compiled and critiqued. Thermal emission measurement methods have been refined and documented. A unique thermal emissivity measurement facility has been constructed and used to measure the emissivity of flowing FLiBe-exposed SS316 samples. A slight increase has been observed in the emissivity compared to unexposed SS316. A dual reflectance-emission setup was partially demonstrated for the purpose of molten fluoride absorption measurement. Preliminary infrared edge data for FLiNaK was obtained. A visible near-infrared (NIR) measurement setup has been constructed and will be used to measure absorption of corrosion products in FLiBe. A new multi-band thermal radiation method using blackbody-weighted spectral bands for the absorption-coefficient was deployed to investigate the thermal performance of the UW Natural Circulation FLiBe Loop under various power levels and operating conditions. The results showed an increasing and decreasing trend in the overall heat transfer with optical thickness with a maximum at optical thickness τ D ≈Ο(1). A forced convection FLiBe loop was constructed for the purpose of studying corrosion and thermal hydraulic phenomena. A vertical cantilever pump is used to achieve up to 5 gpm and 5 m/s flow velocity in hot and cold test sections. A methodology has been established for modeling radiative heat transfer (RHT) in molten salt coolants for the purpose of quantifying the contribution of thermal radiation in conditions relevant to FHR or MSR systems. The comparison between the results obtained for forced and mixed convection showed that in the studied ranges of optical thickness, the maximum effect of RHT on overall heat transfer was 60% improvement observed at z/D = 250 for forced convection; whereas for mixed convection, that value was higher at 170% observed at z/D = 30. The highest effect of RHT on the overall energy distribution was observed at τ D ≈3.4. Negligible effect (<20%) was observed above τ D =18 and below τ D ≈0.002. A CFD code-to-code comparison was conducted using COMSOL Multiphysics at UW-Madison and STAR-CCM at MIT, for a forced and mixed convection laminar flow for a grey medium in a vertical heated using both COMSOL Multiphysics at UW-Madison and STAR-CCM at MIT. The relative error was the largest at τ D =40 for the radiative heat flux, bulk temperature, and Nusselt number. On the other hand, the relative error was the least at τ D =0.002 for the conductive heat flux, bulk temperature, and Nusselt number.

Anderson, Mark↗

Leveraging ARM Data to Improve Models for Predictive Understanding of Energy and Security Challenges

Extreme weather and natural hazards can disrupt the energy sector, affecting demand, generation, transmission, distribution, consumption and operational planning at regional and national scales. These disruptions stem from a broad range of atmospheric phenomena, including winter storms, freezing rain, wet snow loading, severe convection, flooding and landslides, wildfires, prolonged heat, and drought. Many of these same phenomena can also affect national security through impacts to transportation and infrastructure. To support the U.S. Department of Energy (DOE) focus on energy resilience and national security, the Atmospheric Radiation Measurement (ARM) User Facility is uniquely positioned to contribute measurement data, analyses, and modeling frameworks that can significantly improve predictive understanding of these hazards to mitigate their effects. To explore this opportunity, ARM convened a two-part virtual workshop in November 2025. The workshop engaged interdisciplinary experts in atmospheric science, energy systems, modeling, and operations. The goal of the meeting was to engage with these interdisciplinary experts to address three questions: • What are examples of atmospheric processes that represent significant risks to energy security or national security and where are those risks greatest? • What measurements or measurement strategies would improve ARM’s capacity to address these issues? • How can ARM and users of the ARM facility better work with the Energy Exascale Earth System Model (E3SM) and multi-sector modeling communities to apply ARM data to improving E3SM simulations of these phenomena? Participants were asked to submit white papers ahead of the meeting to initiate thinking on these themes and to help organize discussions. Workshop sessions were then organized around themes identified in the white papers. First from the white papers and then through subsequent discussions, workshop participants identified many examples that address the three questions listed above. Participants called out energy system vulnerabilities to weather phenomena such as the impact of freezing rain, strong winds, and excessive heat on power grids. They also noted the effects that weather phenomena could have on energy demand or supply (e.g., through effects of extreme temperatures). They called out security vulnerabilities such as impacts to crops from aerosol-borne pathogens and risks to industry due to melting permafrost in the Arctic. In all, over a dozen meteorological phenomena were linked to energy or security vulnerabilities. For many of the identified phenomena, participants pointed out where ARM was well poised to address issues (e.g., through measurements of cloud microphysics to inform studies of freezing rain) but also noted needs for additional measurements or modified measurement strategies. For example, adaptive scanning of severe weather would be valuable for probing winter storms or severe convection. Participants pointed out the value in integrating external observations with ARM measurements and with applying artificial intelligence (AI) to ARM observation analysis and they advocated for using model simulations to help optimize measurement strategies through Observing System Simulation Experiments (OSSEs). It was clear from the workshop that there are many ways that ARM observations can be used to mitigate energy and security concerns, but meeting participants were also asked to identify what they considered to be the greatest opportunities by ranking issues pertaining to the three workshop questions. This was accomplished through a survey administered to participants between the two virtual sessions. The highest-priority phenomena identified were winter storms, severe convection, and arctic processes. Discussion in the second session, therefore, focused primarily on these three areas, which were most fully developed in exploring ARM opportunities. Nevertheless, it was also clear that ARM has opportunities to contribute to all the identified topics. This report describes the workshop, including input from discussion and white papers (Sections 2 and 3) and a list of priority recommendations (section 4). Many other ideas for ARM contributions are discussed in individual white papers (Appendix D).

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Bursty bulk flows in the inner central plasma sheet

High-speed flows in the inner central plasma sheet (first reported by Baumjohann et al. (1990) are studied, together with the concurrent behavior of the plasma and magnetic field, by using AMPTE/IRM data from about 9 to 19 R(E) in the earth magnetotail. The conclusions drawn from the detailed analysis of a representative event are reinforced by a superposed epoch analysis applied on two years of data. The high-speed flows organize themselves in 10-min time scale flow enhancements called here bursty-bulk flow (BBF) events. Both temporal and spatial effects are responsible for their bursty nature. The flow velocity exhibits peaks of very large amplitude with a characteristic time scale of the order of a minute, which are usually associated with magnetic field dipolarizations and ion temeperature increases. The BBFs represent intervals of enhanced earthward convection and energy transport per unit area in the y-z GSM direction of the order of 5 x 10 exp 19 ergs/R(E-squared).

Angelopoulos, V.↗

Investigation of the Short-Time Variability of Tropical Tropospheric Ozone

Since 1998, a ground based tropospheric ozone lidar has been running at Reunion Island and has been involved with a daily measurement campaign that was performed in the latter part of the biomass burning season, during November-December 1999. The averaged ozone profile obtained during November-December 1 999 agrees well with averaged ozone profile obtained from ozonesondes launch at Reunion during November-December (1992- 2001). Comparing weekly sonde launches (part of the Southern Hemisphere Additional Ozonesondes: SHADOZ program) with the daily ground-based lidar observations shows that some striking features of the day to day variability profiles are not observed in the sonde measurements. Ozone profiles respond to the nature of disturbances which vary from the one day to the next. The vertical ozone distribution at Reunion is examined as a function of prevailing atmospheric circulation. Backtrajectories show that most of the enhanced ozone crossed over biomass burning and convectively active regions in Madagascar and the southern African continent. The analyses of the meteorological data show that ozone stratification profiles are in agreement with the movement of the synoptical situations in November-December 1999. Three different sequences of transport are explained using wind fields. The first sequence from 23 to 25 November is characterized by Northerly transport, the second sequence from 26 to 30 November, the air masses are influenced by meridional transport. The third sequence from 2 to 6 December is characterized by westerly transport associated with the subtropical jet stream. The large standard deviations of lidar profiles in the middle and upper troposphere are in agreement with the upper wind variabilities which evidence passing ridge and trough disturbances. During the transition period between the dry season and the wet season, multiple ozone sources including stratosphere-troposphere exchanges, convection and biomass burning contribute to tropospheric ozone at Reunion Island through sporadic events characterized by a large spatial and temporal variability.

Randriambelo, Tantely↗