Search NASA⌕ Search

SEARCH · Search NASA

Results for “Critical zone”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Hot corrosion behavior of 304 & P91 graded composite transition joint under molten sulfate salts

A novel graded composite transition joint (GCTJ) between AISI 304 stainless steel and ASTM A335 P91 steel, has been demonstrated remarkablely superior creep performance compared to the conventional dissimilar metal weldment (DMW) under equivalent conditions. However, this advance is challenged by hot corrosion under salt deposition at elevated temperatures. Here, this study investigates the hot corrosion behavior of 304&P91 GCTJ exposed to sulfate salts at 700°C. Compared to the 304 steel, corrosion attacks initiate in the P91 triangle, where the dual-phase microstructure of ferrite and tempered martensite significantly influences pitting initiation. Anodic dissolution mainly occurs within the tempered martensite due to more vulnerable sites within the tempered martensite. With prolonged exposure, corrosion propagates across both phases, and the corrosion depth within the P91 triangle is related to the exposed surface area ratio between 304 and P91. Electrochemical analysis reveals the occurrence of galvanic corrosion between the 304 and P91, with a positive linear relationship between anodic dissolution current density (I a ) and the exposed surface area ratio between 304 and P91 in molten salts, further emphasizing the critical impact of this ratio on the corrosion severity in the P91 triangle. These findings underscore the importance of transition zone design optimization in mitigating localized corrosion.

Galvanic corrosion↗
DOE OSTI2026

Enhancing microsegregation during rapid directional solidification through ternary microalloying

The nano-cellular dendritic microstructure formed during rapid directional solidification in powder bed fusion additive manufacturing creates unique properties such as simultaneous improvement in strength and ductility. However, process control of microsegregation features remains challenging due to low sensitivity of critical solidification mechanisms to process parameters. This study leverages microalloying to achieve large changes in dendrite composition, microstructure, and interdendritic zone width during laser powder bed fusion without modifying process parameters. CALPHAD simulations predict that the addition of Zr significantly steepens the solidus line of the dilute Cu-Cr alloy system, leading to enhanced Cr rejection into the melt and greater than 95% reduction in solubility of Cr in the solidified Cu matrix. Experimental validation using time-of-flight secondary ion mass spectrometry and Kelvin probe force microscopy reveals that the ternary alloy containing 0.01 wt% Zr exhibited wider interdendritic regions compared to the binary, a significantly higher number of Cr-rich particles within interdendritic regions, near-complete ejection of oxygen impurities from the matrix, and greater nanoscale work function contrast. These features indicate more aggressive Cr segregation in the presence of Zr and a purer Cu matrix and provide a potentially robust method for engineering the nano-cellular dendritic solidification microstructure.

CALPHAD↗
DOE OSTI2023

Stability Criteria for Self-Propagating Reaction Waves in Co/Al Multilayers

The propagation of self-sustained formation reactions in sputter-deposited Co/Al multilayers is known to exhibit a design-dependent instability. Multilayers having thin bilayers (<55 nm period) exhibit stable propagating waves, whereas those with a larger period react unstably. The specific two-dimensional (2D) instability observed involves the transverse propagation of a band in front of a stalled front commonly referred to as a “spin band.” Previous finite-element studies have shown that these instabilities are thermodynamically driven by the forward conduction of heat away from the flame front. However, the magnitude of that loss is inherently tied to the bilayer design in traditional bimetallic multilayers, which couples any proposed stability criteria to a varying critical diffusion distance. Here, this work utilizes a recently developed class of materials known as “inert-mediated reactive multilayers” to decouple the thermodynamic and kinetic contributions to propagating wave stability by reducing the stored chemical energy density in normally stable bilayer designs. By depositing an inert product phase (B2-CoAl) within the mid-plane of Co and Al reactant layers, spin instabilities arise as a function of both diluted volume and critical diffusion distance. From there, a stability criterion is determined for Co/Al multilayers based on enthalpy loss from the reaction zone, and its physical significance is explored.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Range extender mediates long-distance enhancer activity

Although most mammalian transcriptional enhancers regulate their cognate promoters over distances of tens of kilobases, some enhancers act over distances in the megabase range1. The sequence features that enable such long-distance enhancer–promoter interactions remain unclear. Here we used in vivo enhancer-replacement experiments at the mouse Shh locus to show that short- and medium-range limb enhancers cannot initiate gene expression at long-distance range. We identify a cis-acting element, range extender (REX), that confers long-distance regulatory activity and is located next to a long-range limb enhancer of Sall1. The REX element has no endogenous enhancer activity. However, addition of the REX to other short- and mid-range limb enhancers substantially increases their genomic interaction range. In the most extreme example observed, addition of REX increased the range of an enhancer by an order of magnitude from its native 73 kb to 848 kb. The REX element contains highly conserved [C/T]AATTA homeodomain motifs that are critical for its activity. These motifs are enriched in long-range limb enhancers genome-wide, including the ZRS (zone of polarizing activity (ZPA) regulatory sequence), a benchmark long-range limb enhancer of Shh2. The ZRS enhancer with mutated [C/T]AATTA motifs maintains limb activity at short range, but loses its long-range activity, resulting in severe limb reduction in knock-in mice. In summary, we identify a sequence signature associated with long-range enhancer–promoter interactions and describe a prototypical REX element that is necessary and sufficient to confer long-distance activation by remote enhancers.

Bower, Grace↗
NASA NTRS1986

Lithologic discrimination of volcanic and sedimentary rocks by spectral examination of Landsat TM data from the Puma, Central Andes Mountains

The Central Andes are widely used as a modern example of noncollisional mountain-building processes. The Puna is a high plateau in the Chilean and Argentine Central Andes extending southward from the altiplano of Bolivia and Peru. Young tectonic and volcanic features are well exposed on the surface of the arid Puna, making them prime targets for the application of high-resolution space imagery such as Shuttle Imaging Radar B and Landsat Thematic Mapper (TM). Two TM scene quadrants from this area are analyzed using interactive color image processing, examination, and automated classification algorithms. The large volumes of these high-resolution datasets require significantly different techniques than have been used previously for the interpretation of Landsat MSS data. Preliminary results include the determination of the radiance spectra of several volcanic and sedimentary rock units and the use of the spectra for automated classification. Structural interpretations have revealed several previously unknown folds in late Tertiary strata, and key zones have been targeted to be investigated in the field. The synoptic view of space imagery is already filling a critical gap between low-resolution geophysical data and traditional geologic field mapping in the reconnaissance study of poorly mapped mountain frontiers such as the Puna.

Fielding, E. J.↗
NASA NTRS2000

Detached and Floating-Zone Growth of Semiconductor Crystals on the ISS

Understanding the mechanism of detached Bridgman growth and establishing the growth of large scale germanium-silicon crystals by the float-zone technique are the key points of the project "RDGS - Reduction of Defects in Germanium-Silicon". The contact angle of the melt and the growth angle of the crystal are essential parameters which allow a controlled use of detached growth. The contact angle was determined for a variety of different substrates and melt compositions; pBN showed the highest value for pure germanium as well as for germanium-rich GeSi melts. The growth angle of Ge(sub 0.95) Si(sub 0.05) was measured to be 8.5-10.5 degrees which concurs with the values of pure germanium and silicon, respectively. The temperature dependence and the concentration dependence of the surface tension were determined for concentrations up to 10at% silicon (partial derivative (gamma)/partial derivative T=-0.08 (raised dot) 10(exp -3)N/m (raised dot) K, partial derivative (gamma)/partial derivative (C)=2.2 (raised dot) 10(exp -3)N/m (raised dot) at%). Using these values, the critical Marangoni number indicating the transition to time-dependent thermocapillary flow will be exceeded for the growth of large scale float-zone crystals onboard the ISS. Therefore, suitable tools for flow control are required.

Dold, P.↗
DOE OSTI2026

Correlation between microstructure and residual stress formation in friction stir welded armor steels characterized by neutron diffraction

Friction stir welding (FSW) is a solid-state joining process that minimizes the heat-affected zone (HAZ) compared with fusion-based arc welding, making it well suited for joining martensitic armor steels where hardness and ballistic resistance are critical. This study investigates residual stress formation in three defect-free FSW butt-joint configurations relevant to armored-vehicle fabrication: similar rolled homogeneous armor (RHA–RHA, Case 1), similar high-hardness armor (HHA–HHA, Case 2), and dissimilar HHA–RHA (Case 3) joints produced under temperature-controlled conditions (770 °C). Neutron diffraction was employed to quantify the magnitude and spatial distribution of residual stresses in the longitudinal, transverse, and normal directions and to correlate them with weld microstructure and hardness. Tensile residual stresses were concentrated in the softened HAZ, reaching approximately 300 MPa for Case 2 and 400 MPa for Case 1 (≈50–70 % of the base-metal yield strength; ∼581 MPa for RHA and ∼566 MPa for HHA), while compressive residual stresses dominated the stir zone. The spatial extent of tensile stresses scaled with the width of the softened HAZ, which was largest in the dissimilar HHA–RHA joint and smallest in the HHA–HHA joint. Full-width-at-half-maximum (FWHM) analysis revealed low microstrain in overtempered HAZ regions and high microstrain in the stir zone associated with severe plastic deformation and fresh martensite formation. This work demonstrates that residual stress evolution in FSW of martensitic armor steels is governed not primarily by peak temperature or thermal contraction, as inferred from fusion-welding analogies, but by the competition between transformation-induced volumetric expansion and tempering-induced stress relaxation. The relative dominance of these mechanisms is shown to depend on alloy hardenability and local thermal history, leading to more extensive HAZ softening and broader tensile stress regions in the lower-hardenability RHA steel. These findings establish a transferable mechanistic framework for optimizing solid-state joining strategies in high-strength steels and other transformation-hardening alloys beyond armor applications.

Friction stir welding↗
DOE OSTI2025

Increased flood exposure in the Pacific Northwest following earthquake-driven subsidence and sea-level rise

Climate-driven sea-level rise is increasing the frequency of coastal flooding worldwide, exacerbated locally by factors like land subsidence from groundwater and resource extraction. However, a process rarely considered in future sea-level rise scenarios is sudden (over minutes) land subsidence associated with great (>M8) earthquakes, which can exceed 1 m. Along the Washington, Oregon, and northern California coasts, the next great Cascadia subduction zone earthquake could cause up to 2 m of sudden coastal subsidence, dramatically raising sea level, expanding floodplains, and increasing the flood risk to local communities. Here, we quantify the potential expansion of the 1% floodplain (i.e., the area with an annual flood risk of 1%) under low (~0.5 m), medium (~1 m), and high (~2 m) earthquake-driven subsidence scenarios at 24 Cascadia estuaries. If a great earthquake occurred today, floodplains could expand by 90 km 2 (low), 160 km 2 (medium), or 300 km 2 (high subsidence), more than doubling the flooding exposure of residents, structures, and roads under the high subsidence scenario. By 2100, when climate-driven sea-level rise will compound the hazard, a great earthquake could expand floodplains by 170 km 2 (low), 240 km 2 (medium), or 370 km 2 (high subsidence), more than tripling the flooding exposure of residents, structures, and roads under the high subsidence scenario compared to the 2023 floodplain. Our findings can support decision-makers and coastal communities along the Cascadia subduction zone as they prepare for compound hazards from the earthquake cycle and climate-driven sea-level rise and provide critical insights for tectonically active coastlines globally.

Science & Technology - Other Topics↗
NASA NTRS1976

Excitation of pulsations in the CNO ionization zone of luminous stars

Homogeneous models of luminous hydrogen-and-helium burning stars were constructed on the basis of Carson's new radiative opacities. These opacities exhibit a large 'bump' at moderate temperatures and low densities as a result of the ultimate ionization of the CNO group of elements and they induce in the envelopes of the more massive stars a strong local convection zone, a high central condensation and pulsational instability by means of the kappa-mechanism above a certain stellar mass. This critical mass for Population I hydrogen-burning stars is probably greater than 45 solar masses for the fundamental radial mode, with the overtones being more stable, while for the homogeneous helium-burning stars, the critical mass lies in the range 2-6 solar masses for all modes up to at least the third overtone. Convection is alternately treated by means of the mixing length theory and the assumption of strict adiabaticity.

Stothers, R.↗
NASA NTRS1994

Variations of a global constraint factor in cracked bodies under tension and bending loads

Elastic-plastic finite-element analyses were used to calculate stresses and displacements around a crack in finite-thickness plates for an elastic-perfectly plastic material. Middle- and edge-crack specimens were analyzed under tension and bending loads. Specimens were 1.25 to 20 mm thick with various widths and crack lengths. A global constraint factor alpha(sub g), an averaged normal-stress to flow-stress ratio over the plastic region, was defined to simulate three-dimensional (3D) effects in two-dimensional (2D) models. For crack lengths and uncracked ligament lengths greater than four times the thickness, the global constraint factor was found to be nearly a unique function of a normalized stress-intensity factor (related to plastic-zone size to thickness ratio) from small- to large-scale yielding conditions for various specimen types and thickness. For crack length-to-thickness ratios less than four, the global constraint factor was specimen type, crack length and thickness dependent. Using a 2D strip-yield model and the global constraint factors, plastic-zone sizes and crack-tip displacements agreed reasonably well with the 3D analyses. For a thin sheet aluminum alloy, the critical crack-tip-opening angle during stable tearing was found to be independent of specimen type and crack length for crack length-to-thickness ratios greater than 4.

Newman, J. C., Jr.↗
NASA NTRS2010

Genesis of Pre-Hurricane Felix (2007): The Role of the Easterly Wave Critical Layer - Part 1

The formation of pre Hurricane Felix (2007) in a tropical easterly wave is examined in a two-part study using the Weather Research and Forecasting (WRF) model with a high-resolution nested grid configuration that permits the representation of cloud system processes. The simulation commences during the wave stage of the precursor African easterly-wave disturbance. Here the simulated and observed developments are compared, while in Part II of the study various large-scale analyses, physical parameterizations, and initialization times are explored to document model sensitivities. In this first part the authors focus on the wave/vortex morphology, its interaction with the adjacent intertropical convergence zone complex, and the vorticity balance in the neighborhood of the developing storm. Analysis of the model simulation points to a bottom-up development process within the wave critical layer and supports the three new hypotheses of tropical cyclone formation proposed recently by Dunkerton, Montgomery, and Wang. It is shown also that low-level convergence associated with the ITCZ helps to enhance the wave signal and extend the "wave pouch" from the jet level to the top of the atmospheric boundary layer. The region of a quasi-closed Lagrangian circulation within the wave pouch provides a focal point for diabatic merger of convective vortices and their vortical remnants. The wave pouch serves also to protect the moist air inside from dry air intrusion, providing a favorable environment for sustained deep convection. Consistent with the authors' earlier findings, the tropical storm forms near the center of the wave pouch via system-scale convergence in the lower troposphere and vorticity aggregation. Components of the vorticity balance are shown to be scale dependent, with the immediate effects of cloud processes confined more closely to the storm center than the overturning Eliassen circulation induced by diabatic heating, the influence of which extends to larger radii.

Wang, Zhuo↗
NASA NTRS2010

Defining an Abrasion Index for Lunar Surface Systems as a Function of Dust Interaction Modes and Variable Concentration Zones

Unexpected issues were encountered during the Apollo era of lunar exploration due to detrimental abrasion of materials upon exposure to the fine-grained, irregular shaped dust on the surface of the Moon. For critical design features involving contact with the lunar surface and for astronaut safety concerns, operational concepts and dust tolerance must be considered in the early phases of mission planning. To systematically define material selection criteria, dust interaction can be characterized by two-body or three-body abrasion testing, and subcategorically by physical interactions of compression, rolling, sliding and bending representing specific applications within the system. Two-body abrasion occurs when a single particle or asperity slides across a given surface removing or displacing material. Three-body abrasion occurs when multiple particles interact with a solid surface, or in between two surfaces, allowing the abrasives to freely rotate and interact with the material(s), leading to removal or displacement of mass. Different modes of interaction are described in this paper along with corresponding types of tests that can be utilized to evaluate each configuration. In addition to differential modes of abrasion, variable concentrations of dust in different zones can also be considered for a given system design and operational protocol. These zones include: (1) outside the habitat where extensive dust exposure occurs, (2) in a transitional zone such as an airlock or suitport, and (3) inside the habitat or spacesuit with a low particle count. These zones can be used to help define dust interaction frequencies, and corresponding risks to the systems and/or crew can be addressed by appropriate mitigation strategies. An abrasion index is introduced that includes the level of risk, R, the hardness of the mineralogy, H, the severity of the abrasion mode, S, and the frequency of particle interactions, F.

Kobrick, Ryan L.↗
DOE OSTI2023

Critical Minerals and Rare Earth Elements in Powder River Basin Coal and Associated Sediments, Wyoming, USA

The demand for Rare Earth Elements (REE) and other critical minerals (CM) required for consumer products, defense-related applications, and low-carbon energy technology is increasing rapidly. Unconventional sources of REE/CM, such as coal and associated sediments, could prove very important in building resilient and ethical energy supply chains. Utilization of existing infrastructure and the highly trained energy workforce in traditionally coal-producing regions in emerging REE/CM industries could provide an economic boost to coal communities. The Powder River Basin (PRB) of Wyoming and Montana, USA, is a prime candidate to investigate the feasibility of extracting REE/CM from coal and associated sediments. The PRB hosts thick (>50 ft) coal seams that are mined at the surface and more than 40% of the coal produced in the US comes from the PRB. Geochemical data from multiple locations in PRB coal systems show REE enrichments at the top and bottom margins and at internal partings in coal seams. This trend is exhibited in two cores drilled at Peabody’s North Antelope Rochelle mine in the east central PRB, the largest coal mine in the world. The two cores include the Wyodak Anderson coal zone as well as the over- and underlying shale units. The Wyodak Anderson coal zone is part of the Paleocene Tongue River Member of the Fort Union Formation in which most of the coal resources in the PRB reside. To understand REE/CM enrichment, a total of 188 samples were analyzed for their major and trace element chemistry from both cores, sampled at one foot or smaller intervals. The two cores show distinct REE enrichments in the uppermost and lowermost 3 to 12 feet of the Wyodak Anderson coal zone, as well as enrichments in bounding carbonaceous shale units. Total REE+Y (REY) concentrations as high as 2510 ppm (all concentrations reported on an ash basis), or 15 times average upper continental crust values (Taylor and McLennan, 1995), were identified in the coal. Partings within the coal zone also show relative REE enrichment, with concentrations up to 485 ppm REY. The remaining interior portions of the coal zone contain lower concentrations of REE, resulting in an average REY for all samples from both cores of 288 ppm. The proportion of high-value critical REE (Nd, Eu, Tb, Dy, Er, Y) compared to REY averages 36%, which is higher than the critical REE proportion in average upper continental crust of 33%. This finding highlights the relative enrichment of middle REE compared to light REE, in contrast to many conventional REE deposits. Other critical trace elements that show enrichments above average upper continental crust include Ga, Nb, and V. Within this sample suite, Ga concentrations range from 2.8 to 213 ppm, with an average of 29 ppm; Nb from 2.9 to 166 ppm with an average of 26 ppm; and V from 27 to 1695 ppm with an average of 188 ppm. None of these trace elements exhibit the same distinct pattern of enrichment as REE in the upper and lower bounding layers and internal partings of the coal zone. Major element chemistry indicates CaO concentrations averaging ~25% in the coal samples from both cores, which has important implications for the extractability of REE from coal. Calcium-rich PRB coal has been shown to be more amenable to REE extraction than lower calcium coals by some methods (Stuckman et al., 2019; Taggart et al., 2016), highlighting the importance of REE extractability in addition to REE concentration in assessing the value of potential feedstocks. The original coal resource in the PRB is estimated at 1.16 trillion short tons (Luppens et al., 2015). Thus, the PRB represents an important potential unconventional source for REE and other CM. Ongoing research to fully characterize the REE/CM resource, identify enrichment mechanisms, and further develop and scale extraction technologies is necessary to understand the full potential of this resource.

Phillips, Erin↗
NASA NTRS2002

Plant Production Systems for Microgravity: Critical Issues in Water, Air, and Solute Transport Through Unsaturated Porous Media

This NASA Technical Memorandum is a compilation of presentations and discussions in the form of minutes from a workshop entitled 'Plant Production Systems for Microgravity: Critical Issues in Water, Air, and Solute Transport Through Unsaturated Porous Media' held at NASA's Johnson Space Center, July 24-25, 2000. This workshop arose from the growing belief within NASA's Advanced Life Support Program that further advances and improvements in plant production systems for microgravity would benefit from additional knowledge of fundamental processes occurring in the root zone. The objective of the workshop was to bring together individuals who had expertise in various areas of fluid physics, soil physics, plant physiology, hardware development, and flight tests to identify, discuss, and prioritize critical issues of water and air flow through porous media in microgravity. Participants of the workshop included representatives from private companies involved in flight hardware development and scientists from universities and NASA Centers with expertise in plant flight tests, plant physiology, fluid physics, and soil physics.

Steinberg, Susan L.↗
DOE OSTI2025

U.S. Solar Siting Regulation and Zoning Ordinances (2025)

A machine readable collection of documented solar siting ordinances at the state and local (e.g., county, township) level throughout the United States. The data were compiled using the Infrastructure Continuous Ordinance Mapping for Planning and Siting Systems (INFRA-COMPASS) tool, which leverages Large Language Models (LLMs) to automate the collection of local codes and ordinances applicable to energy infrastructure. URLs for the ordinance source documents are included in the Solar Ordinances spreadsheet. The GeoPackage file included below contains the jurisdiction shapes for each ordinance. Note that the GeoPackage file is formatted for ingestion by NLR's reVX setbacks tool and therefore does not contain any of the state-level regulations. NOTE: This data was collected with the help of generative AI. The Large Language Models used for this effort make mistakes. Always validate the data for critical use cases. This data is an update to a previously developed database of wind ordinances found in OEDI Submission 5734: see the "U.S. Solar Siting Regulation and Zoning Ordinances 2022" link below. INFRA-COMPASS version used for collection: v0.11.3 LLMs used for collection: GPT-4.1, GPT-4.1 mini, GPT-4.1 nano

14 SOLAR ENERGY↗
DOE OSTI2025

U.S. Wind Siting Regulation and Zoning Ordinances (2025)

A machine readable collection of documented wind siting ordinances at the state and local (e.g., county, township) level throughout the United States. The data were compiled using the Infrastructure Continuous Ordinance Mapping for Planning and Siting Systems (INFRA-COMPASS) tool, which leverages Large Language Models (LLMs) to automate the collection of local codes and ordinances applicable to energy infrastructure. URLs for the ordinance source documents are included in the Wind Ordinances spreadsheet. The GeoPackage file included below contains the jurisdiction shapes for each ordinance. Note that the GeoPackage file is formatted for ingestion by NREL's reVX setbacks tool and therefore does not contain any of the state-level regulations. NOTE: This data was collected with the help of generative AI. The Large Language Models used for this effort make mistakes. Always validate the data for critical use cases. This data is an update to a previously developed database of wind ordinances found in OEDI Submission 5733: see the "U.S. Wind Siting Regulation and Zoning Ordinances 2022" link below. INFRA-COMPASS version used for collection: v0.8.2 LLMs used for collection: GPT-4.1, GPT-4.1 mini, GPT-4.1 nano, GPT-4o mini

17 WIND ENERGY↗
NASA NTRS1992

Large-scale flows and solar luminosity variations

The influence of a variation of the large-scale flows present in the solar convection zone on solar luminosity is examined. The convection zone is modeled as a plane-parallel fluid layer through which heat is transported by small-scale convection as described by the mixing-length theory, and the large-scale flow is treated as a small perturbation superposed on the layer. The excess energy transported by the flow is calculated and is found to depend critically on the mixing length. The excess luminosity is found to be carried as small-scale convective heat flow rather than bulk enthalpy transport throughout all but the top few percent of the convection zone. The results are applied to the specific examples of supergranules and giant cells.

Arendt, S.↗
DOE OSTI2021

An Improved Zhang's Dynamic Water Balance Model Using Budyko-Based Snow Representation for Better Streamflow Predictions

Understanding the water balance of a catchment in relation to its regional climate forcings and catchment characteristics is critical for predicting current and future water resources amid changing climate and land cover. This study intends to improve Zhang's monthly water balance model (a physics-based conceptual hydrologic model) that reflects the physical partitioning process of the hydrological cycle at the basin level based on regional climate and catchment characteristics. The existing model does not include snow process and has confronted evident limitations in snow-affected areas, which is a critical aspect since snowmelt water has been a significant source of water resources for many regions, especially in the temperate and frigid zones. We introduce a snow module based on surface energy balance and Budyko-limits on melting and combine it with the existing water balance equations. Moreover, monthly parameterization is applied to the model to better explain the time-varying hydrological characteristics of a catchment. The proposed model involves five different monthly parameters, which determine the physical partitioning process of the hydrological cycle, and they are regionally calibrated and validated under Budyko-type constraints. The model is applied to 1,210 basins across the continental United States (CONUS), and the simulated streamflow is compared to the observed data. The proposed model significantly outperformed the original model, improving the median NSE by 31% (from 0.51 to 0.67) and increasing the number of catchments with an acceptable NSE by 58%. The spatial variability of the basin characteristics across the CONUS is also investigated based on the calibrated parameters.

54 ENVIRONMENTAL SCIENCES↗