Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational data assimilation”

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 433 records · Page 24

Long-term measurements in a mixed-grass prairie reveal a change in soil organic carbon recalcitrance and its environmental sensitivity under warming

Soil respiration, the major pathway for ecosystem carbon (C) loss, has the potential to enter a positive feedback loop with the atmospheric CO2 due to climate warming. For reliable projections of climate-carbon feedbacks, accurate quantification of soil respiration and identification of mechanisms that control its variability are essential. Process-based models simulate soil respiration as functions of belowground C input, organic matter quality, and sensitivity to environmental conditions. However, evaluation and calibration of process-based models against the long-term in situ measurements are rare. Here, we evaluate the performance of the Terrestrial ECOsystem (TECO) model in simulating total and heterotrophic soil respiration measured during a 16-year warming experiment in a mixed-grass prairie; calibrate model parameters against these and other measurements collected during the experiment; and explore whether the mechanisms of C dynamics have changed over the years. Calibrating model parameters against observations of individual years substantially improved model performance in comparison to pre-calibration simulations, explaining 79–86% of variability in observed soil respiration. Interannual variation of the calibrated model parameters indicated increasing recalcitrance of soil C and changing environmental sensitivity of microbes. Overall, we found that (1) soil organic C became more recalcitrant in intact soil compared to root-free soil; (2) warming offset the effects of increasing C recalcitrance in intact soil and changed microbial sensitivity to moisture conditions. Furthermore, these findings indicate that soil respiration may decrease in the future due to C quality, but this decrease may be offset by warming-induced changes in C cycling mechanisms and their responses to moisture conditions.

59 BASIC BIOLOGICAL SCIENCES↗

Specification of the near-Earth space environment with SHIELDS

Here, predicting variations in the near-Earth space environment that can lead to spacecraft damage and failure is one example of “space weather” and a big space physics challenge. A project recently funded through the Los Alamos National Laboratory (LANL) Directed Research and Development (LDRD) program aims at developing a new capability to understand, model, and predict Space Hazards Induced near Earth by Large Dynamic Storms, the SHIELDS framework. The project goals are to understand the dynamics of the surface charging environment (SCE), the hot (keV) electrons representing the source and seed populations for the radiation belts, on both macro- and micro-scale. Important physics questions related to particle injection and acceleration associated with magnetospheric storms and substorms, as well as plasma waves, are investigated. These challenging problems are addressed using a team of world-class experts in the fields of space science and computational plasma physics, and state-of-the-art models and computational facilities. A full two-way coupling of physics-based models across multiple scales, including a global MHD (BATS-R-US) embedding a particle-in-cell (iPIC3D) and an inner magnetosphere (RAM-SCB) codes, is achieved. New data assimilation techniques employing in situ satellite data are developed; these provide an order of magnitude improvement in the accuracy in the simulation of the SCE. SHIELDS also includes a post-processing tool designed to calculate the surface charging for specific spacecraft geometry using the Curvilinear Particle-In-Cell (CPIC) code that can be used for reanalysis of satellite failures or for satellite design.

54 ENVIRONMENTAL SCIENCES↗

Stratospheric Impact on the Northern Hemisphere Winter and Spring Ozone Interannual Variability in the Troposphere

In this study we use ozone and stratospheric ozone tracer simulations from the high-resolution (0.5°x0.5°) Goddard Earth Observing System, Version 5 (GEOS-5) in a replay mode to study the impact of stratospheric ozone on tropospheric ozone interannual variability (IAV). We use these simulations in conjunction with ozone sonde measurements from 1990 to 2016 during the winter and spring seasons. The simulations include a stratospheric ozone tracer (StratO3) to aid 15in the evaluation of the impact of stratospheric ozone IAV on the IAV of tropospheric ozone at different altitudes and location. The model is in good agreement with the observed interannual variation of tropospheric ozone, except for the post-Pinatubo period(1992-1994)over the region of North America. Ozone sonde data show a negative ozone anomaly in 1992-1994 following the Pinatubo eruption, with recovery thereafter. The simulated anomaly is only half the magnitude of that observed. Our analysis suggests that the simulated Stratosphere-troposphere exchange (STE)flux deduced from the analysis 20might be too strong over the North American (50°N-70°N) region after the Mt. Pinatubo eruption in the early 1990s,maskingthe impact of lower stratospheric ozone concentration on tropospheric ozone. European ozone sonde measurements show a similar but weaker ozone depletion after the Mt. Pinatubo eruption, which is fully reproduced by the model. Analysis based on the stratospheric ozone tracer (StratO3) identifies differences in strength and vertical extent of stratospheric ozone impact on the tropospheric ozone interannual variation (IAV)between North America and Europe. Over NorthAmerica25stations, the StratO3IAV has a significant impact on tropospheric ozone from the upper to lower troposphere and explains about 60% and 66% of the simulated ozone IAV at 400 hPa, and ~11% and 34% at 700 hPa in winter and spring respectively. Over European stations, the influence is limited to the middle to upper troposphere, and becomes much smaller at 700 hPa. The Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2)assimilated fields exhibit strong longitudinal variations over northern hemisphere (NH) mid-high latitudes, with lower tropopause height and 30lower geopotential height over North America than Europe. These variations associated with the relevant variations in the location of tropospheric jet flows are responsible for the longitudinal differences in the stratospheric ozone impact, with stronger effects over North America than over Europe

tropospheric ozone IAV, Model, ozonesonde, stratos↗

Comparison of Global Precipitation Estimates across a Range of Temporal and Spatial Scales

Characteristics of precipitation estimates for rate and amount from three global high-resolution precipitation products (HRPPs), four global climate data records (CDRs), and four reanalyses are compared. All datasets considered have at least daily temporal resolution. Estimates of global precipitation differ widely from one product to the next, with some differences likely due to differing goals in producing the estimates. HRPPs are intended to produce the best snapshot of the precipitation estimate locally. CDRs of precipitation emphasize homogeneity over instantaneous accuracy. Precipitation estimates from global reanalyses are dynamically consistent with the large-scale circulation but tend to compare poorly to rain gauge estimates since they are forecast by the reanalysis system and precipitation is not assimilated. Regional differences among the estimates in the means and variances are as large as the means and variances, respectively. Even with similar monthly totals, precipitation rates vary significantly among the estimates. Temporal correlations among datasets are large at annual and daily time scales, suggesting that compensating bias errors at annual and random errors at daily time scales dominate the differences. However, the signal-to-noise ratio at intermediate (monthly) time scales can be large enough to result in high correlations overall. It is shown that differences on annual time scales and continental regions are around 0.8 mm day -1 , which corresponds to 23 W m -2 . These wide variations in the estimates, even for global averages, highlight the need for better constrained precipitation products in the future.

54 ENVIRONMENTAL SCIENCES↗

Systematic Differences between Satellite-Based Presipitation Climatologies over the Tropical Oceans

Since the beginning of the World Climate Research Program's Global Precipitation Climatology Project (GPCP) satellite remote sensing of precipitation has made dramatic improvements, particularly for tropical regions. Data from microwave and infrared sensors now form the most critical input to precipitation data sets and can be calibrated with surface gauges to so that the strengths of each data source can be maximized in some statistically optimal sense. It is clear however that there still remain significant uncertainties with satellite precipitation retrievals which limit their usefulness for many purposes. Systematic differences i'A tropical precipitation estimates have been brought to light in comparison activities such as the GPCP Algorithm Intercomparison Project and more recent Wetnet Precipitation Intercomparison Project 3. These uncertainties are assuming more importance because of the demands for validation associated with global climate modeling and data assimilation methodologies. The objective of the present study is to determine the physical basis for systematic differences in spatial structure of tropical precipitation as portrayed by several different satellite-based data sets. The study is limited to oceanic regions only and deals primarily with aspects of spatial variability. We are specifically interested in why MSU channel 1 and GPI precipitation differences are so striking over the Eastern Pacific ITCZ and why they both differ from other microwave emission-based precipitation estimates from SSM/I and a scattering-based deep convective ice index from MSU channel 2. Our results to date have shown that MSU channel I precipitation estimates are biased high over the Eastern Pacific ITCZ because of two factors: (1) the hypersensitivity of this frequency to cloud water in contrast to falling rain drops, and (2) unaccounted for scattering effects by precipitation-size ice which depresses the signal of the liquid water emission. Likewise, cold cloud top climatologies such as the GPI show an excess (a deficit) in estimated rainfall over the E. Pacific ITCZ (Warm Pool region). We show that these algorithms need to account for regionally varying heights (or temperatures) at which tropical convection detrains to form cirrus shields. A second objective we pursue is to identify variations in the macroscale cloud physical and thermodynamic properties of precipitation regimes" and relate these differences to tropical dynamical mechanisms of tropical heat and moisture balance. Finally, we interpret the algorithm differences and their associations with tropical dynamics in terms of WCRP GPCP goals for constructing precipitation climatologies.

Robertson, Frankin R.↗

Five-Year (2004-2009)Observations of Upper Tropospheric Water Vapor and Cloud Ice from MLS and Comparisons with GEOS-5 Analyses

This paper gives an overview of August 2004 through July 2009 upper tropospheric (UT) water vapor (H2O) and ice water content (IWC) from the Aura Microwave Limb Sounder (MLS) and comparisons with outputs from the NASA Goddard Earth Observing System Version 5 (GEOS-5) data assimilation system. Both MLS and GEOS-5 show that high values of H2O and IWC at 215 to 147 hPa are associated with areas of deep convection. They exhibit good (within approximately 15%) agreement in IWC at these altitudes, but GEOS-5 H2O is approximately 50% (215 hPa) to approximately 30% (147 hPa) larger than MLS, possibility due to its higher temperatures at these altitudes. GOES-5 produces a weaker intertropical convergence zone than MLS, while a seasonally-migrating band of tropical deep convection is clearly evident in both the MLS and GEOS-5 UT H2O and IWC. MLS and GEOS-5 both show spatial anti-correlation between IWC and H2O at 100 hPa, where less H2O is associated with low temperatures in regions of tropical convection. At 100 hPa, GEOS-5 produces 50% less IWC and 15% less H2O in the tropics, and approximately 20% more H2O in the extra-tropics, than does MLS. Behavior of the 100 hPa H2O, which exhibits a quasi-biennial oscillation, appears consistent with it being controlled by temperature. The seasonal cycle in the vertical transport of tropical mean H2O from approximately 147 hPa to approximately 10 hPa appears much stronger in MLS than in GEOS-5. The UT IWC and H2O interannual variations, from both MLS and GEOS-5, show clear imprints of the El Nino-Southern Oscillation.

Jiang, Jonathan H.↗

A Comparison of Satellite Based, Modeled Derived Daily Solar Radiation Data with Observed Data for the Continental US

Many applications of simulation models and related decision support tools for agriculture and natural resource management require daily meteorological data as inputs. Availability and quality of such data, however, often constrain research and decision support activities that require use of these tools. Daily solar radiation (SRAD) data are especially problematic because the instruments require electronic integrators, accurate sensors are expensive, and calibration standards are seldom available. The Prediction Of Worldwide Energy Resources (NASA/POWER; power.larc.nasa.gov) project at the NASA Langley Research Center estimates daily solar radiation based on data that are derived from satellite observations of outgoing visible radiances and atmospheric parameters based upon satellite observations and assimilation models. The solar data are available for a global 1 degree x 1 degree coordinate grid. SRAD can also be estimated based on attenuation of extraterrestrial radiation (Q0) using daily temperature and rainfall data to estimate the optical thickness of the atmosphere. This study compares daily solar radiation data from NASA/POWER (SRADNP) with instrument readings from 295 stations (SRADOB), as well as with values that were estimated with the WGENR solar generator. WGENR was used both with daily temperature and precipitation records from the stations reporting solar data and records from the NOAA Cooperative Observer Program (COOP), thus providing two additional sources of solar data, SRADWG and SRADCO. Values of SRADNP for different grid cells consistently showed higher correlations (typically 0.85 to 0.95) with SRADOB data than did SRADWG or SRADCO for sites within the corresponding cells. Mean values of SRADOB, SRADWG and SRADNP for sites within a grid cell usually were within 1 MJm-2d-1 of each other, but NASA/POWER values averaged 1.1 MJm-2d-1 lower than SRADOB. The magnitude of this bias was greater at lower latitudes and during summer months and may be at least partially explained by assumptions in ambient aerosol properties. Overall, the NASA/POWER solar radiation data are a promising resource for regional modeling studies where realistic accounting of historic variation is required.

White, Jeffrey W.↗

Satellite-Assisted Particulate Matter (SAPM) for the Models, In situ, and Remote Sensing of Aerosols (MIRA) Working Group

Models, In situ, and Remote sensing of Aerosols (MIRA) is an international working group that provides a forum for collaborations amongst these three atmospheric aerosol communities. MIRA consists of a collection of interdisciplinary and independently funded Topic/Project groups with clear goals and generally characterized by requests for additional scientific data. The Satellite-Assisted Particulate Matter (SAPM) Topic group focuses on the study of fine particulate matter (PM2.5), as it is a major contributor to air pollution and negatively impacts human health. Our group aims to provide intercomparisons of various methods and techniques for estimating surface PM2.5 assisted by satellite remote sensors (passive and active), global aerosol models, and in situ aerosol measurements. The overall motivation of our work is to enhance PM2.5 coverage beyond in situ ground stations, which can be limited. In this presentation, we provide an overview of various PM2.5 estimation approaches by current SAPM team members, as each approach has its own strengths and limitations. One approach uses near-surface aerosol extinction retrievals from the spaceborne Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) and an assumed value of mass extinction efficiency to derive PM2.5 concentrations (Toth et al 2019). A study using this method found a promising agreement (R=0.60; slope=0.89) with U.S. Environmental Protection Agency (EPA) in situ PM2.5 measurements over the contiguous United States (CONUS) for a 12-year (2007-2018) period (Toth et al 2022). In a different approach combining spaceborne lidar and a model, the Cloud Aerosol Transport System (CATS) lidar and Goddard Earth Observing System (GEOS) model are blended in an ensemble variational assimilation scheme to retrieve aerosol extinction and convert speciated mass concentrations to total PM2.5. A recent study (Matus et al in prep.) found agreement within 2 μg/m3 over the CONUS on the 2016 annual mean when compared to EPA ground-based measurements. Another study explored trends in city aerosol optical depth (AOD) from a spaceborne passive remote sensor, Moderate Resolution Imaging Spectroradiometer (MODIS), and found that these trends agree well with trends in city surface PM2.5, which provides confidence in the use of city AOD for city PM2.5 trend studies (Vohra et al 2021). In a follow-on study, steep and significant trends (~2.5 to 7.8% a-1) in AOD are found for several South Asian cities (e.g., Bangalore and Hyderabad), suggesting rapid growth in PM2.5 (Vohra et al 2022). Also in this presentation, we report results on our current collective SAPM study of India, a country that exhibits high levels of PM2.5 pollution. We show the temporal (i.e., seasonal) and spatial variability of PM2.5 over India using the CALIOP and CATS+model approaches, comparisons between the estimated PM2.5 from the retrieval schemes, and validation with in situ PM2.5 ground-based observations from the Central Pollution Control Board in India. This region provides us with an excellent case study to test our PM2.5 retrievals in heavily polluted scenes and over a complex and varying topography. The SAPM Topic group actively seeks international participants/collaborators in our group, including those working with in situ aerosol measurements (e.g., ACTRIS). We are interested in aerosol datasets in order to either improve our PM2.5 estimates from various approaches (e.g., using in situ mass scattering/absorption coefficient and aerosol hygroscopic properties for various aerosol species) and/or validate the PM2.5 estimates (e.g., using in situ ground-based PM2.5 concentrations).

Travis D Toth↗

Highlights of TOMS Version 9 Total Ozone Algorithm

The fundamental basis of TOMS total ozone algorithm was developed some 45 years ago by Dave and Mateer. It was designed to estimate total ozone from satellite measurements of the backscattered UV radiances at few discrete wavelengths in the Huggins ozone absorption band (310-340 nm). Over the years, as the need for higher accuracy in measuring total ozone from space has increased, several improvements to the basic algorithms have been made. They include: better correction for the effects of aerosols and clouds, an improved method to account for the variation in shape of ozone profiles with season, latitude, and total ozone, and a multi-wavelength correction for remaining profile shape errors. These improvements have made it possible to retrieve total ozone with just 3 spectral channels of moderate spectral resolution (approx. 1 nm) with accuracy comparable to state-of-the-art spectral fitting algorithms like DOAS that require high spectral resolution measurements at large number of wavelengths. One of the deficiencies of the TOMS algorithm has been that it doesn't provide an error estimate. This is a particular problem in high latitudes when the profile shape errors become significant and vary with latitude, season, total ozone, and instrument viewing geometry. The primary objective of the TOMS V9 algorithm is to account for these effects in estimating the error bars. This is done by a straightforward implementation of the Rodgers optimum estimation method using a priori ozone profiles and their error covariances matrices constructed using Aura MLS and ozonesonde data. The algorithm produces a vertical ozone profile that contains 1-2.5 pieces of information (degrees of freedom of signal) depending upon solar zenith angle (SZA). The profile is integrated to obtain the total column. We provide information that shows the altitude range in which the profile is best determined by the measurements. One can use this information in data assimilation and analysis. A side benefit of this algorithm is that it is considerably simpler than the present algorithm that uses a database of 1512 profiles to retrieve total ozone. These profiles are tedious to construct and modify. Though conceptually similar to the SBUV V8 algorithm that was developed about a decade ago, the SBUV and TOMS V9 algorithms differ in detail. The TOMS algorithm uses 3 wavelengths to retrieve the profile while the SBUV algorithm uses 6-9 wavelengths, so TOMS provides less profile information. However both algorithms have comparable total ozone information and TOMS V9 can be easily adapted to use additional wavelengths from instruments like GOME, OMI and OMPS to provide better profile information at smaller SZAs. The other significant difference between the two algorithms is that while the SBUV algorithm has been optimized for deriving monthly zonal means by making an appropriate choice of the a priori error covariance matrix, the TOMS algorithm has been optimized for tracking short-term variability using month and latitude dependent covariance matrices.

Bhartia, Pawan↗

Improving E3SM Land Model Photosynthesis Parameterization via Satellite SIF, Machine Learning, and Surrogate Modeling

The parameterization of key photosynthesis parameters is one of the key uncertain sources in modeling ecosystem gross primary productivity (GPP). Solar-induced chlorophyll fluorescence (SIF) offers a good proxy for GPP since it marks the actual process of photosynthesis; while machine learning (ML) provides a robust approach to model the GPP-SIF relationship. Here, we trained the boosted regressing tree (BRT) and the Random Forest ML models with Greenhouse Gases Observing Satellite SIF data and in situ GPP observations from 49 eddy covariance towers. These trained ML GPP-SIF models were fed into the Energy Exascale Earth System Model (E3SM) Land Model (ELM) to generate ELM-simulated global SIF estimates, which were then benchmarked against satellite SIF observations with a surrogate modeling approach. Our results indicated good modeling performance of the ML-based GPP-SIF relationship. The ELM model when fed with the ML GPP-SIF models also can well predict the spatial-temporal variations in SIF. We also found high model accuracy for the surrogate modeling. Model parameter sensitivity analysis suggested that the fraction of leaf nitrogen in RuBisCO (flnr) is the most sensitive parameter to the SIF; other sensitive parameters include the Ball-Berry stomatal conductance slope (mbbopt) and the vcmax entropy (vcmaxse). The posterior uncertainty in simulated GPP was greatly reduced after benchmarking, and the model produced improved spatial patterns of mean GPP relative to FLUXCOM GPP. Our integrated approach provides a new avenue for improving land models and using remote-sensing SIF, which can be further improved in the future with more ground- and satellite-based observations.

54 ENVIRONMENTAL SCIENCES↗

Do dynamic global vegetation models capture the seasonality of carbon fluxes in the Amazon basin? A data-model intercomparison

To predict forest response to long-term climate change with high confidence requires that dynamic global vegetation models (DGVMs) be successfully tested against ecosystem response to short-term variations in environmental drivers, including regular seasonal patterns. Here, we used an integrated dataset from four forests in the Brasil flux network, spanning a range of dry-season intensities and lengths, to determine how well four state-of-the-art models (IBIS, ED2, JULES, and CLM3.5) simulated the seasonality of carbon exchanges in Amazonian tropical forests. We found that most DGVMs poorly represented the annual cycle of gross primary productivity (GPP), of photosynthetic capacity (Pc), and of other fluxes and pools. Models simulated consistent dry-season declines in GPP in the equatorial Amazon (Manaus K34, Santarem K67, and Caxiuanã CAX); a contrast to observed GPP increases. Model simulated dry-season GPP reductions were driven by an external environmental factor, ‘soil water stress’ and consequently by a constant or decreasing photosynthetic infrastructure (Pc), while observed dry-season GPP resulted from a combination of internal biological (leaf-flush and abscission and increased Pc) and environmental (incoming radiation) causes. Moreover, we found models generally overestimated observed seasonal net ecosystem exchange (NEE) and respiration (Re) at equatorial locations. In contrast, a southern Amazon forest (Jarú RJA) exhibited dry-season declines in GPP and Re consistent with most DGVMs simulations. While water limitation was represented in models and the primary driver of seasonal photosynthesis in southern Amazonia, changes in internal biophysical processes, light-harvesting adaptations (e.g., variations in leaf area index (LAI) and increasing leaf-level assimilation rate related to leaf demography), and allocation lags between leaf and wood, dominated equatorial Amazon carbon flux dynamics and were deficient or absent from current model formulations. In conclusion, correctly simulating flux seasonality at tropical forests requires a greater understanding and the incorporation of internal biophysical mechanisms in future model developments.

54 ENVIRONMENTAL SCIENCES↗

Quality and Control of Water Vapor Winds

Water vapor imagery from the geostationary satellites such as GOES, Meteosat, and GMS provides synoptic views of dynamical events on a continual basis. Because the imagery represents a non-linear combination of mid- and upper-tropospheric thermodynamic parameters (three-dimensional variations in temperature and humidity), video loops of these image products provide enlightening views of regional flow fields, the movement of tropical and extratropical storm systems, the transfer of moisture between hemispheres and from the tropics to the mid- latitudes, and the dominance of high pressure systems over particular regions of the Earth. Despite the obvious larger scale features, the water vapor imagery contains significant image variability down to the single 8 km GOES pixel. These features can be quantitatively identified and tracked from one time to the next using various image processing techniques. Merrill et al. (1991), Hayden and Schmidt (1992), and Laurent (1993) have documented the operational procedures and capabilities of NOAA and ESOC to produce cloud and water vapor winds. These techniques employ standard correlation and template matching approaches to wind tracking and use qualitative and quantitative procedures to eliminate bad wind vectors from the wind data set. Techniques have also been developed to improve the quality of the operational winds though robust editing procedures (Hayden and Veldon 1991). These quality and control approaches have limitations, are often subjective, and constrain wind variability to be consistent with model derived wind fields. This paper describes research focused on the refinement of objective quality and control parameters for water vapor wind vector data sets. New quality and control measures are developed and employed to provide a more robust wind data set for climate analysis, data assimilation studies, as well as operational weather forecasting. The parameters are applicable to cloud-tracked winds as well with minor modifications. The improvement in winds through use of these new quality and control parameters is measured without the use of rawinsonde or modeled wind field data and compared with other approaches.

Jedlovec, Gary J.↗

20th century changes in carbon isotopes and water-use efficiency: tree-ring-based evaluation of the CLM4.5 and LPX-Bern models

Measurements of the stable carbon isotope ratio ( δ 13 C) on annual tree rings offer new opportunities to evaluate mechanisms of variations in photosynthesis and stomatal conductance under changing CO 2 and climate conditions, especially in conjunction with process-based biogeochemical model simulations. The isotopic discrimination is indicative of the ratio between the CO 2 partial pressure in the intercellular cavities and the atmosphere ( c i / c a ) and of the ratio of assimilation to stomatal conductance, termed intrinsic water-use efficiency (iWUE). We performed isotope-enabled simulations over the industrial period with the land biosphere module (CLM4.5) of the Community Earth System Model and the Land Surface Processes and Exchanges (LPX-Bern) dynamic global vegetation model. Results for C3 tree species show good agreement with a global compilation of δ 13 C measurements on leaves, though modeled 13 C discrimination by C3 trees is smaller in arid regions than measured. A compilation of 76 tree-ring records, mainly from Europe, boreal Asia, and western North America, suggests on average small 20th century changes in isotopic discrimination and in c i / c a and an increase in iWUE of about 27% since 1900. LPX-Bern results match these century-scale reconstructions, supporting the idea that the physiology of stomata has evolved to optimize trade-offs between carbon gain by assimilation and water loss by transpiration. In contrast, CLM4.5 simulates an increase in discrimination and in turn a change in iWUE that is almost twice as large as that revealed by the tree-ring data. Factorial simulations show that these changes are mainly in response to rising atmospheric CO 2 . The results suggest that the downregulation of c i / c a and of photosynthesis by nitrogen limitation is possibly too strong in the standard setup of CLM4.5 or that there may be problems associated with the implementation of conductance, assimilation, and related adjustment processes on long-term environmental changes.

54 ENVIRONMENTAL SCIENCES↗

A Comparison of the Age-Spectra from Data Assimilation Models

We use kinematic and diabatic back trajectory calculations, driven by winds from a general circulation model (GCM) and two different data assimilation systems (DAS), to compute the age spectrum at three latitudes in the lower stratosphere. The age-spectra are compared to chemical transport model (CTM) calculations, and the mean ages from all of these studies are compared to observations. The age spectra computed using the GCM winds show a reasonably well-isolated tropics in good agreement with observations; however, the age spectra determined from the DAS differ from the GCM spectra. For the diabatic trajectory calculations, the age spectrum is too broad as a result of too much exchange between the tropics and mid-latitudes. The age spectrum determined using the kinematic trajectory calculation is less broad and lacks an age offset; both of these features are due to excessive vertical dispersion of parcels. The tropical and mid-latitude mean age difference between the diabatically and kinematically determined age-spectra is about one year, the former being older. The CTM calculation of the age spectrum using the DAS winds shows the same dispersive characteristics of the kinematic trajectory calculation. These results suggest that the current DAS products will not give realistic trace gas distributions for long integrations; they also help explain why the mean ages determined in a number of previous DAS driven CTM's are too young compared with observations. Finally, we note trajectory-generated age spectra show significant age anomalies correlated with the seasonal cycles, and these anomalies can be linked to year-to-year variations in the tropical heating rate. These anomalies are suppressed in the CTM spectra suggesting that the CTM transport is too diffusive.

Schoeberl, Mark R.↗

Evaluation of High Mountain Asia-Land Data Assimilation System (version 1) from 2003 to 2016: 2. The impact of assimilating satellite-based snow cover and freeze/thaw observations into a land surface model

This second paper of the two-part series focuses on demonstrating the impact of assimilating satellite-based snow cover and freeze/thaw observations into the hyper-resolution, offline terrestrial modeling system used for the High Mountain Asia (HMA) region from 2003 to 2016. To this end, this study systematically evaluates a total of six sets of 0.01° (∼1 km) model simulations forced by different precipitation forcings, with and without the dual assimilation scheme enabled, at point-scale, basin-scale, and domain-scale. The key variables of interest include surface net shortwave radiation, surface net longwave radiation, skin temperature, near-surface soil temperature, snow depth, snow water equivalent (SWE), and total runoff. First, the point-scale assessment is mainly conducted via evaluating against ground-based measurements. In general, the assimilation enabled estimates are better than no-assimilation counterparts. Second, the basin-scale runoff assessment demonstrates that across three snow-dominated basins, the assimilation enabled experiment yields systematic improvements in all goodness-of-fit statistics through mitigating the negative effects brought by the fixed long-term precipitation correction factors. For example, when forced by the bias-corrected precipitation, the assimilation-enabled experiment improves the bias by 69%, the root-mean-squared error by 30%, and the unbiased root-mean-squared error by 18% (relative to the no-assimilation counterpart). Finally, the domainscale assessment is conducted via evaluating against satellite-based SWE and skin temperature products. Both sets of domain-scale analysis further corroborate the findings in the point-scale evaluations. Overall, this study suggests the benefits of the proposed multi-variate assimilation system in improving the cryospherichydrological process within a land surface model for use in HMA.

Yuan Xue↗

Leaf functional traits highlight phenotypic variation of two tree species in the urban environment

Urbanization is transforming landscapes globally, altering environmental conditions that affect ecosystem functioning, particularly in urban areas where trees are crucial for regulating microclimates, improving air quality, and sustaining biodiversity. This study investigates the environmental differences and tree leaf structure and morphology in urban and suburban sites in the Chicago Metropolitan Region. The leaf functional traits of Norway Maple and Little − leaved Linden were studied in three locations in the summer of 2023: an urban park (University of Illinois Chicago, Chicago, IL), a suburban park (Morton Arboretum, Lisle, IL), and a suburban residential site (Lombard, IL). The urban site had higher daytime and nighttime air, and land surface temperatures compared to the suburban sites with significant fluctuations observed across the sites. Cumulative growing degree days, a measure of potential photosynthetically active days, were also higher in the urban park than in the suburban sites between March and August. Norway Maple trees growing in the urban site displayed higher specific leaf area (SLA) and lower leaf dry matter content (LDMC) than in the suburban sites, resulting in thinner leaves. Similarly, Little−leaved Linden trees in the suburban residential site displayed higher SLA and lower LDMC than those in the suburban park. The values of gas exchange traits − namely photosynthetic assimilation, transpiration rates, and stomatal conductance − of Norway Maple were higher at the urban site compared to suburban sites as temperatures increased during the summer. Norway Maple gas exchange values decreased as the growing season progressed, as expected by ontogeny. In contrast, Little−leaved Linden maintained similar leaf gas exchange values throughout the growing season. Both species in the urban site exhibited lower instantaneous water use efficiency and reduced LDMC, suggesting greater water loss in response to elevated temperatures compared to suburban park and residential sites. Comparisons with existing global trait databases emphasize the need for localized data to accurately capture site−specific responses. Although some traits aligned with database values, others deviated significantly, underscoring the importance of comprehensive, site−specific datasets for robust ecosystem modeling and management strategies.

Cho, Ahram↗

Evaluation of High Mountain Asia-Land Data Assimilation System (version 1) from 2003 to 2016: 2. The impact of assimilating satellite-based snow cover and freeze/thaw observations into a land surface model

This second paper of the two-part series focuses on demonstrating the impact of assimilating satellite-based snow cover and freeze/thaw observations into the hyper-resolution, offline terrestrial modeling system used for the High Mountain Asia (HMA) region from 2003 to 2016. To this end, this study systematically evaluates a total of six sets of 0.01° (∼1 km) model simulations forced by different precipitation forcings, with and without the dual assimilation scheme enabled, at point-scale, basin-scale, and domain-scale. The key variables of interest include surface net shortwave radiation, surface net longwave radiation, skin temperature, near-surface soil temperature, snow depth, snow water equivalent (SWE), and total runoff. First, the point-scale assessment is mainly conducted via evaluating against ground-based measurements. In general, the assimilation enabled estimates are better than no-assimilation counterparts. Second, the basin-scale runoff assessment demonstrates that across three snow-dominated basins, the assimilation enabled experiment yields systematic improvements in all goodness-of-fit statistics through mitigating the negative effects brought by the fixed long-term precipitation correction factors. For example, when forced by the bias-corrected precipitation, the assimilation-enabled experiment improves the bias by 69%, the root-mean-squared error by 30%, and the unbiased root-mean-squared error by 18% (relative to the no-assimilation counterpart). Finally, the domain-scale assessment is conducted via evaluating against satellite-based SWE and skin temperature products. Both sets of domain-scale analysis further corroborate the findings in the point-scale evaluations. Overall, this study suggests the benefits of the proposed multi-variate assimilation system in improving the cryospheric-hydrological process within a land surface model for use in HMA.

Yuan Xue↗

Leaf traits and canopy structure together explain canopy functional diversity: an airborne remote sensing approach

Plant functional diversity is strongly connected to photosynthetic carbon assimilation in terrestrial ecosystems. However, many of the plant functional traits that regulate photosynthetic capacity, including foliar nitrogen concentration and leaf mass per area, vary significantly between and within plant functional types and vertically through forest canopies, resulting in considerable landscape-scale heterogeneity in three dimensions. Hyperspectral imagery has been used extensively to quantify functional traits across a range of ecosystems but is generally limited to providing information for top-of-canopy leaves only. On the other hand, lidar data can be used to retrieve the vertical structure of forest canopies. Because these data are rarely collected at the same time, there are unanswered questions about the effect of forest structure on the 3-D spatial patterns of functional traits across ecosystems. In the United States, the National Ecological Observatory Network's Airborne Observation Platform (NEON AOP) provides an opportunity to address this structure-function relationship by collecting lidar and hyperspectral data together across a variety of ecoregions. With a fusion of hyperspectral and lidar data from the NEON AOP and field-collected foliar trait data, we assessed the impacts of forest structure on spatial patterns of N. In addition, we examine the influence of abiotic gradients and management regimes on top of canopy %N and total canopy N (i.e. the total amount of N (g/m 2 ) within a forest canopy) at a NEON site consisting of a mosaic of open longleaf pine and dense broadleaf deciduous forests. Here, our resulting maps suggest that in contrast with top-of-canopy values, total canopy N variation is dampened across this landscape resulting in relatively homogeneous spatial patterns. At the same time, we found that leaf functional diversity and canopy structural diversity showed distinct dendritic patterns related to the spatial distribution of plant functional types.

54 ENVIRONMENTAL SCIENCES↗