Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “In situ model”

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

The Eruption of a Magnetic Flux Rope Observed By Solar Orbiter and Parker Solar Probe

Magnetic flux ropes are a key component of coronal mass ejections, forming the core of these eruptive phenomena. However, determining whether a flux rope is present prior to eruption onset and, if so, the rope’s handedness and the number of turns that any helical field lines make is difficult without magnetic field modelling or in-situ detection of the flux rope. We present two distinct observations of plasma flows along a filament channel on 4 and 5 September 2022 made using the Solar Orbiter spacecraft. Each plasma flow exhibited helical motions in a right-handed sense as the plasma moved from the source active region across the solar disk to the quiet Sun, suggesting that the magnetic configuration of the filament channel contains a flux rope with positive chirality and at least one turn. The length and velocity of the plasma flow increased from the first to the second observation, suggesting evolution of the flux rope, with the flux rope subsequently erupting within ∼5 hours of the second plasma flow. The erupting flux rope then passed over the Parker Solar Probe spacecraft during its Encounter 13, enabling in-situ diagnostics of the structure. Although complex and consistent with the flux rope erupting from underneath the heliospheric current sheet, the in-situ measurements support the inference of a right-handed flux rope from remote-sensing observations. These observations provide a unique insight into the eruption and evolution of a magnetic flux rope near the Sun.

Dynamics↗

Evolution of Martian Volcanism and Lithosphere: Geochemical Insights from Noachian to Amazonian Volcanic Terrains

The lithosphere of Mars has played a crucial role in shaping its geological features, influencing volcanic activity, climate shifts, and the planet's early habitability. Understanding the evolution of Mars' litho-sphere through different geological eons has been challenging due to weathering and resurfacing affecting the study of early Noachian volcanic chemistry. Therefore, in our study we have focused on recently discovered Noachian volcanic terranes [1-2], along with Hesperian and Amazonian volcanic terranes [3], to trace the evolution of the Martian lithosphere and thermal flux. The composition of igneous rock formed from the eruption of magma tends to preserve the record of thermal properties viz. pressure, temperature, or degree of partial melting at which it forms. Therefore, Martian volcanic provinces are of great geologic interest; they have been active throughout its history, from Noachian (>3.7 Ga) to the Late Amazonian (<500 Ma) [4]. Earlier studies focused on the evolved magmatism from Hesperian to Amazonian [4], but the type and style of Noachian-aged volcanism remain to be understood. We have investigated the geochemical compositions of Noachian volcanic provinces and compared them with relatively younger volcanic provinces on Mars using remote sensing (Mars Odyssey Gamma Ray and Neutron Spectrometer suite-GRS) and in-situ observations. Furthermore, petrologic modeling is performed to understand magmatic processes and thermal evolution. Our findings reveal distinct Pressure-Temperature (P-T) conditions across different Martian geological eras. In the Noachian terranes, we model P-T conditions ranging from 1.3-1.6 GPa and temperatures between 1350-1390°C. Conversely, the Hesperian volcanic terrane exhibits variations between 1.6-1.7 GPa and temperatures spanning 1370-1395°C. For the Amazonian volcanic terrane, P-T conditions ranged notably higher, from 1.9-2.8 GPa with temperatures between 1380-1415°C. Based on our modeling, we estimated the corresponding litho-spheric thicknesses and depths of melting for each eon. In the Noachian, the lithospheric thickness ranged from 110-135 km. Meanwhile, for the Hesperian terranes, it varied between 120-235 km, and for the Amazonian, it extended from 160-235 km. We also calculated partial melting percentages (F), observing a range of 8-12% for the Noachian, 10-11% for the Hesperian, and 10-12% for the Amazonian. Mantle potential temperatures (Tp) were calculated to illustrate variations across these geological eons. Our estimates of heat flux depicted temporal changes, showing ranges from 51-65 mW/m 2 for the Noachian, 38-45 mW/m 2 for the Hesperian, and 27-39 mW/m 2 for the Amazonian. Our study implies that the lithosphere remained consistently uniform until the Hesperian eon, indicating sustained weaker volcanic systems from mid-Noachian to Hesperian. This contrasts with fewer deep-seated plumes observed in the Amazonian era. These variations suggest potential implications for Martian climatic conditions which show a milder climate prevailed until the Hesperian, differing from Mars' current cold, arid conditions.

Mars↗

Stratospheric free chlorine measured by balloon-borne in situ resonance fluorescence

Eight balloon-borne in situ measurements of ClO in the stratosphere are analyzed and are compared with recent model calculations. While the use of in situ stratospheric studies of free radicals to test models by comparing observed and predicted concentration profiles is essential for a prognosis of changes in stratospheric ozone, resulting from future changes in stratospheric ozone, such studies provide only limited insight into the nature of stratospheric photochemistry, because natural variability and the large number of fast reactions which compete in the coupling among the key radicals frustrate a detailed comparison between a mean distribution provided by the models and an instantaneous distribution provided by a single observation.

Anderson, J. G.↗

Evaluation of EPIC Oxygen Bands Stability With Radiative Transfer Simulations Over the South Pole

The Earth Polychromatic Imaging Camera (EPIC) onboard the Deep Space Climate Observatory (DSCOVR) satellite orbiting the Sun at the Lagrange-1 point was launched without onboard calibration systems. Vicarious calibration is conducted for 8 of the 10 UV/VIS/NIR channels using other low earth orbiting satellite instruments, while its two O 2 bands (688 nm and 764nm) rely on indirect moon-view calibrations because the same narrow-band O 2 bands are not readily available from other in-flight instruments. This study compares EPIC measurements from the four O 2 bands aiming at examining sensor stability over a uniquely suited location, i.e., the permanently snow-covered South Pole. The study utilizes radiative transfer model simulations with in-situ atmospheric soundings taken at South Pole during months of December and January from 2015 to 2022. The absolute discrepancy between the model simulations and observations is less than 1.0% for the two reference bands, but 5.75% and 15.63% for the 688nm, and 764 nm absorption bands, respectively. The simulated A-band and B-band ratios are 16.09% and 4.74% higher than that from the observations. Various sensitivities are conducted to estimate possible contributions to the discrepancies from input atmospheric profiles, spectral surface albedos and surface BRDF. While none of the input uncertainties is likely to account for the large discrepancies in the oxygen absorption bands, a small shift in the instrument response function could be the main reason for these biases. On the other hand, the model simulations are able to capture systematic variations with observed angular measurements and explain the multi-year trends found in observed O 2 band ratios due to satellite orbit shifting. When model simulated contributions from the angle variations are deducted from the observed O 2 band ratios, the residual O 2 band ratios are found to be stable since 2015.

EPIC↗

Advances in Land Data Assimilation at the NASA Goddard Space Flight Center

Research in land surface data assimilation has grown rapidly over the last decade. In this presentation we provide a brief overview of key research contributions by the NASA Goddard Space Flight Center (GSFC). The GSFC contributions to land assimilation primarily include the continued development and application of the Land Information System (US) and the ensemble Kalman filter (EnKF). In particular, we have developed a method to generate perturbation fields that are correlated in space, time, and across variables and that permit the flexible modeling of errors in land surface models and observations, along with an adaptive filtering approach that estimates observation and model error input parameters. A percentile-based scaling method that addresses soil moisture biases in model and observational estimates opened the path to the successful application of land data assimilation to satellite retrievals of surface soil moisture. Assimilation of AMSR-E surface soil moisture retrievals into the NASA Catchment model provided superior surface and root zone assimilation products (when validated against in situ measurements and compared to the model estimates or satellite observations alone). The multi-model capabilities of US were used to investigate the role of subsurface physics in the assimilation of surface soil moisture observations. Results indicate that the potential of surface soil moisture assimilation to improve root zone information is higher when the surface to root zone coupling is stronger. Building on this experience, GSFC leads the development of the Level 4 Surface and Root-Zone Soil Moisture (L4_SM) product for the planned NASA Soil-Moisture-Active-Passive (SMAP) mission. A key milestone was the design and execution of an Observing System Simulation Experiment that quantified the contribution of soil moisture retrievals to land data assimilation products as a function of retrieval and land model skill and yielded an estimate of the error budget for the SMAP L4_SM product. Terrestrial water storage observations from GRACE satellite system were also successfully assimilated into the NASA Catchment model and provided improved estimates of groundwater variability when compared to the model estimates alone. Moreover, satellite-based land surface temperature (LST) observations from the ISCCP archive were assimilated using a bias estimation module that was specifically designed for LST assimilation. As with soil moisture, LST assimilation provides modest yet statistically significant improvements when compared to the model or satellite observations alone. To achieve the improvement, however, the LST assimilation algorithm must be adapted to the specific formulation of LST in the land model. An improved method for the assimilation of snow cover observations was also developed. Finally, the coupling of LIS to the mesoscale Weather Research and Forecasting (WRF) model enabled investigations into how the sensitivity of land-atmosphere interactions to the specific choice of planetary boundary layer scheme and land surface model varies across surface moisture regimes, and how it can be quantified and evaluated against observations. The on-going development and integration of land assimilation modules into the Land Information System will enable the use of GSFC software with a variety of land models and make it accessible to the research community.

Reichle, Rolf↗

Mesospheric H2O and H2O2 densities inferred from in situ positive ion composition measurement

A model for production and loss of oxonium ions in the high-latitude D-region is developed, based on the observed excess of 34(+) which has been interpreted as H2O2(+). The loss mechanism suggested in the study is the attachment of N2 and/or CO2 in three-body reactions. Furthermore, mesospheric water vapor and H2O2 densities are inferred from measurements of four high-latitude ion compositions, based on the oxonium model. Mixing ratios of hydrogen peroxide of up to two orders of magnitude higher than previous values were obtained. A number of reactions, reaction constants, and a block diagram of the oxonium ion chemistry in the D-region are given.

Kopp, E.↗

Seawater Debye Model Function at L-Band and Its Impact on Salinity Retrieval From Aquarius Satellite Data

A model function of seawater, which specifies the dielectric constant of seawater as a function of salinity, temperature, and frequency, is important for the retrieval of sea surface salinity using satellite data. In 2017, a model function has been developed based on measurement data at 1.4134 GHz using a third-order polynomial expression in salinity ( S ) and temperature ( T ). Although the model showed improvements in salinity retrieval, it had an inconsistent behavior between partitioned salinities. To improve the stability of the model, new dielectric measurements of seawater have been made recently over a broad range of salinities and temperatures to expand the data set used for developing the model function. The structure of the model function has been changed from a polynomial expansion in S and T to a physics-based model consisting of a Debye molecular resonance term plus a conductivity term. Each unknown parameter is expressed in S and T based on the expanded measurement data set. Physical arguments have been used to limit the number of unknown coefficients in these expressions to improve the stability of the model function. The new model function has been employed in the retrieval algorithm of the Aquarius satellite mission to obtain a global salinity map. The retrieved salinity using a different model function is compared with in situ data collected by Argo floats to evaluate the impact and the performance of model functions. The results indicate that the new model function has significant improvements in salinity retrieval compared with other existing models.

Seawater↗

Assessment of Clear Sky Radiative Forcing in the Caribbean Region Using an Aerosol Dispersion Model and Ground Radiometry During Puerto Rico Dust Experiment

This study investigates the surface and top of the atmosphere solar radiative forcing by long-range transport of Saharan dust. The calculations of radiative forcing are based on measurements collected in the Puerto Rico Dust Experiment (PRIDE) carried out during July, 2000. The purpose of the experiment was the characterization of the Saharan dust plume, which frequently reaches the Caribbean region during the summer. The experiment involved the use of three approaches to study the plume: space and ground based remote sensing, airborne and ground based in-situ measurements and aerosol dispersion modeling. The diversity of measuring platforms provides an excellent opportunity for determination of the direct effect of dust on the clear sky radiative forcing. Specifically, comparisons of heating rates, surface and TOA fluxes derived from the Navy global aerosol dispersion model NAAPS (NRL Aerosol Analysis and Prediction System) and actual measurements of fluxes from ground and space based platforms are shown. In addition, the direct effect of dust on the clear sky radiative forcing is modeled. The extent and time of evolution of the radiative properties of the plume are computed with the aerosol concentrations modeled by NAAPS. Standard aerosol parameterizations, as well as in-situ composition and size distributions measured during PRIDE, are utilized to compute the aerosol optical depth, single scattering albedo and asymmetry factor. Radiative transfer computations are done with an in-house modified spectral radiative transfer code (Fu-Liou). The code includes gas absorption and cloud particles (ice and liquid phase) and it allows the input of meteorological data. The code was modified to include modules for the aerosols contribution to the calculated fluxes. This comparison study helps to narrow the current uncertainty in the dust direct radiative forcing, as recently reported in the 2001 IPCC assessment.

Gasso, Santiago↗

Watershed and Hydrodynamic Modeling for Evaluating the Impact of Land Use Change on Submerged Aquatic Vegetation and Seagrasses in Mobile Bay

There is a continued need to understand how human activities along the northern Gulf of Mexico coast are impacting the natural ecosystems. The gulf coast is experiencing rapid population growth and associated land cover/land use change. Mobile Bay, AL is a designated pilot region of the Gulf of Mexico Alliance (GOMA) and is the focus area of many current NASA and NOAA studies, for example. This is a critical region, both ecologically and economically to the entire United States because it has the fourth largest freshwater inflow in the continental USA, is a vital nursery habitat for commercially and recreational important fisheries, and houses a working waterfront and port that is expanding. Watershed and hydrodynamic modeling has been performed for Mobile Bay to evaluate the impact of land use change in Mobile and Baldwin counties on the aquatic ecosystem. Watershed modeling using the Loading Simulation Package in C++ (LSPC) was performed for all watersheds contiguous to Mobile Bay for land use Scenarios in 1948, 1992, 2001, and 2030. The Prescott Spatial Growth Model was used to project the 2030 land use scenario based on observed trends. All land use scenarios were developed to a common land classification system developed by merging the 1992 and 2001 National Land Cover Data (NLCD). The LSPC model output provides changes in flow, temperature, sediments and general water quality for 22 discharge points into the Bay. These results were inputted in the Environmental Fluid Dynamics Computer Code (EFDC) hydrodynamic model to generate data on changes in temperature, salinity, and sediment concentrations on a grid with four vertical profiles throughout the Bay s aquatic ecosystems. The models were calibrated using in-situ data collected at sampling stations in and around Mobile bay. This phase of the project has focused on sediment modeling because of its significant influence on light attenuation which is a critical factor in the health of submerged aquatic vegetation. The impact of land use change on sediment concentrations was evaluated by analyzing the LSPC and EFDC sediment simulations for the four land use scenarios. Such analysis was also performed for storm and non-storm periods. In- situ data of total suspended sediments (TSS) and light attenuation were used to develop a regression model to estimate light attenuation from TSS. This regression model was used to derive marine light attenuation estimates throughout Mobile bay using the EFDC TSS outputs. The changes in sediment concentrations and associated impact on light attenuation in the aquatic ecosystem were used to perform an ecological analysis to evaluate the impact on seagreasses and Submerged Aquatic Vegetation (SAV) habitat. This is the key product benefiting the Mobile Bay coastal environmental managers that integrates the influences of sediments due to land use driven flow changes with the restoration potential of SAVs.

Estes, Maurice G.↗

The Backscattering Contribution of Soybean Pods at L-band

L-band (1.25 GHz) radar measurements of a soybean canopy indicate that the emergence of seed pods is a significant contributor to the backscatter during the late stages of the growing season. In order to validate the measured data, a realistic scattering model of the soybean canopy is developed. The parameters of the soybean canopy and underlying soil used in the model vary over the growing season based on in situ measurements. Scattering amplitudes for soybean leaves are modeled analytically by using a thin disk approximation; stem and pods are jointly modeled using a numerical electromagnetic field solver. These scattering amplitudes are together incorporated into a coherent scattering model to obtain the backscattering coefficient for VV- and HH-polarizations. The modeling results show good agreement with the radar field measurements, having RMSEs of 0.51 dB for VV-pol and 1.1 dB for HH-pol. Both measured data and modeled results show that the change of soil moisture can be accurately monitored by L-band backscatter. It is also found that the difference between HH- and VV-polarized backscatter increases as the size of the soybean pods becomes larger. A method is developed here to estimate the number of pods in a soybean canopy based on polarimetric radar backscatter at L-band.

Radar backscatter↗

Land Surface Verification Toolkit (LVT) - A Generalized Framework for Land Surface Model Evaluation

Model evaluation and verification are key in improving the usage and applicability of simulation models for real-world applications. In this article, the development and capabilities of a formal system for land surface model evaluation called the Land surface Verification Toolkit (LVT) is described. LVT is designed to provide an integrated environment for systematic land model evaluation and facilitates a range of verification approaches and analysis capabilities. LVT operates across multiple temporal and spatial scales and employs a large suite of in-situ, remotely sensed and other model and reanalysis datasets in their native formats. In addition to the traditional accuracy-based measures, LVT also includes uncertainty and ensemble diagnostics, information theory measures, spatial similarity metrics and scale decomposition techniques that provide novel ways for performing diagnostic model evaluations. Though LVT was originally designed to support the land surface modeling and data assimilation framework known as the Land Information System (LIS), it also supports hydrological data products from other, non-LIS environments. In addition, the analysis of diagnostics from various computational subsystems of LIS including data assimilation, optimization and uncertainty estimation are supported within LVT. Together, LIS and LVT provide a robust end-to-end environment for enabling the concepts of model data fusion for hydrological applications. The evolving capabilities of LVT framework are expected to facilitate rapid model evaluation efforts and aid the definition and refinement of formal evaluation procedures for the land surface modeling community.

Kumar, Sujay V.↗

Earth System Modeling and Field Experiments in the Arctic-Boreal Zone - Report from a NASA Workshop

Early climate modeling studies predicted that the Arctic Ocean and surrounding circumpolar land masses would heat up earlier and faster than other parts of the planet as a result of greenhouse gas-induced climate change, augmented by the sea-ice albedo feedback effect. These predictions have been largely borne out by observations over the last thirty years. However, despite constant improvement, global climate models have greater difficulty in reproducing the current climate in the Arctic than elsewhere and the scatter between projections from different climate models is much larger in the Arctic than for other regions. Biogeochemical cycle (BGC) models indicate that the warming in the Arctic-Boreal Zone (ABZ) could lead to widespread thawing of the permafrost, along with massive releases of CO2 and CH4, and large-scale changes in the vegetation cover in the ABZ. However, the uncertainties associated with these BGC model predictions are even larger than those associated with the physical climate system models used to describe climate change. These deficiencies in climate and BGC models reflect, at least in part, an incomplete understanding of the Arctic climate system and can be related to inadequate observational data or analyses of existing data. A workshop was held at NASA/GSFC, May 22-24 2012, to assess the predictive capability of the models, prioritize the critical science questions; and make recommendations regarding new field experiments needed to improve model subcomponents. This presentation will summarize the findings and recommendations of the workshop, including the need for aircraft and flux tower measurements and extension of existing in-situ measurements to improve process modeling of both the physical climate and biogeochemical cycle systems. Studies should be directly linked to remote sensing investigations with a view to scaling up the improved process models to the Earth System Model scale. Data assimilation and observing system simulation studies should be used to guide the deployment pattern and schedule for inversion studies as well. Synthesis and integration of previously funded Arctic-Boreal projects (e.g., ABLE, BOREAS, ICESCAPE, ICEBRIDGE, ARCTAS) should also be undertaken. Such an effort would include the integration of multiple remotely sensed products from the EOS satellites and other resources.

Sellers, Piers↗

Studies of the intermediate and deep circulation in the western equatorial Atlantic

This proposal concerns the preparation and design of an experiment, the objective of which is to improve our knowledge of the intermediate and deep circulation in the western equatorial Atlantic Ocean. We shall focus on the description of the western boundary currents, of their crossing with the equator, on the estimation of their mass and heat fluxes, and their seasonal and interannual variations. We will use satellite altimetric data, tomographic measurements, and in situ observations (current measurements, hydrology, and floaters). We propose a feasibility study and the definition of a strategy based on a high-resolution Geophysical Fluid Dynamics Laboratory (GFDL) numerical model to define which in situ measurements are necessary to optimally complete the altimetric observations.

Desaubies, Yves↗

Assessment of climate variability of the Greenland ice sheet: Integration of in situ and satellite data

The proposed research involves the application of multispectral satellite data in combination with ground truth measurements to monitor surface properties of the Greenland ice sheet which are essential for describing the energy and mass of the ice sheet. Several key components of the energy balance are parameterized using satellite data and in situ measurements. The analysis will be done for a ten year time period in order to get statistics on the seasonal and interannual variations of the surface processes and the climatology. Our goal is to investigate to what accuracy and over what geographic areas large scale snow properties and radiative fluxes can be derived based upon a combination of available remote sensing and meterological data sets. Data analysis showed the following results: (1)cloud classification based on longwave sky radiation revealed that overcast sky occurred for 25% of the time in winter, and for 15% in spring and summer respectively (winter and summer both show the same occurrence of clear sky of approximately 26%); (2) comparison of aerodynamic profile method with eddy correlation method to derive sensible and latent heat flux showed good agreement in the diurnal cycle and the turbulent fluxes were underestimated with the aerodynamic method by 10 - 30% as compared to the in situ eddy flux method; (3) the katabatic wind shows a distinct diurnal cycle with a maximum in the morning (7-9 h solar time) and a minimum in the later afternoon (18 h solar time); (4) snow grain size was modeled with a surface energy balance model (SNTHERM) and compared with in situ measurements. Sharp decreases in the modeled snow grain size, caused by accumulation events such as precipitation and deposition, could be verified with observational data; (4) radiative transfer modeling of firn supports our beliefs that the observed trends in 18 and 19 GHz passive microwave brightness temperatures are attributable to accumulation rate changes (modeling also indicates the above relationship is detectable because of the presence of depth hoar; (5) snow melt can be detected by a distinct signal in the passive microwave cross-polarized gradient ratio (19h-37v)/(19h+37v) and has been used for wet/dry snow classification; (6) top of the atmosphere (TOA) broadband albedos were derived from AVHRR visible and near infrared reflectances for the entire ice sheet from May 1990 - June 1991, and the highest albedo values are found along the southeast coast of the ice sheet which is consistent with the summer peak of precipitation due to onshore flow loaded with high water vapor content (TOA albedo values dropped to around 40% along the south-western coast during July and August due to bare ice surfaces); and (7) the net all-wave radiation balance at the top of the atmosphere is negative over the entire ice sheet except for the summer months - June-July-August (in June, the net radiation balance is slightly positive over the dry snow areas--15 W/m2).

Steffen, K.↗

Clotho Tessera, Venus: A fragment of Fortuna Tessera

Clotho Tessera, adjacent to southeast Lakshmi Planum, may provide additional evidence for lateral crustal motions, and a model for the origin of small tessera fragments. Clotho Tessera and Lakshmi Planum are so noticeably different, and in such close proximity, it is difficult to derive a reasonable model of their formation in situ. Squeezing of material out from beneath Lakshmi has been suggested as an origin for Moira Tessera, which is also adjacent to Lakshmi and 1400 km west of Clotho. However, a logical model of juxtaposition of the two different terrains, originally from points once distant, can be made for Clotho and Lakshmi (and perhaps other small tesserae as well). It is suggested that Clotho Tessera was once part of Fortuna Tessera, but was cut off by a transcurrent fault zone (the DLZ) striking perpendicular to the Sigrun rift and carried westward where it collided with Lakshmi Planum (forming Danu Montes). A gravity anomaly along the southern border of Lakshmi, in the area of Danu Montes, was interpreted as indicating subduction there, providing additional supporting evidence for the collision hypothesis. Diffusion of the DLZ with proximity to Sigrun Fossae may be due to either higher ductility near the postulated Sigrun rift, or to burial by flows away from the rift nearer to Valkyrie Fossae. Other possible examples of migrating tesserae occur elsewhere: small pieces of Ananke Tessera can be fit back together as though they had rifted apart, and the spreading apart of Ananke and Virilis Tesserae has been suggested because of their symmetric locations about the axis of an inferred spreading zone. Other tessera fragments appear to have been isolated by rifting, with little, if any, significant lateral motion (e.g., Meni and Tellus Tesserae, and Thethus and Fortuna Tesserae). The migrating terrain model for Clotho Tessera supports Sukhanov's interpretation of tesseral fragments as rafts of lighter crustal material.

Kozak, Richard C.↗

Technical Report Series on Global Modeling and Data Assimilation: Global Surface Ocean Carbon Estimates in a Model Forced by MERRA - Volume 31

MERRA products were used to force an established ocean biogeochemical model to estimate surface carbon inventories and fluxes in the global oceans. The results were compared to public archives of in situ carbon data and estimates. The model exhibited skill for ocean dissolved inorganic carbon (DIC), partial pressure of ocean CO2 (pCO2) and air-sea fluxes (FCO2). The MERRA-forced model produced global mean differences of 0.02% (approximately 0.3 microns) for DIC, -0.3% (about -1.2 (micro) atm; model lower) for pCO2, and -2.3% (-0.003 mol C/sq m/y) for FCO2 compared to in situ estimates. Basin-scale distributions were significantly correlated with observations for all three variables (r=0.97, 0.76, and 0.73, P<0.05, respectively for DIC, pCO2, and FCO2). All major oceanographic basins were represented as sources to the atmosphere or sinks in agreement with in situ estimates. However, there were substantial basin-scale and local departures.

Gregg, Watson W.↗

Backscatter model for the unusual radar properties of the Greenland Ice Sheet

A number of planetary objects exhibit large radar reflectivity and polarization ratios, and more recently, a similar behavior has been observed over a vast portion of the Earth's surface: the percolation facies of the Greenland Ice Sheet. Surface-based ranging radar data and snow stratigraphy studies demonstrated that the radar properties of that portion of Greenland are caused by enhanced scattering from massive, large, solid-ice bodies buried in the top few meters of the dry, cold, clean snowy surface of the ice sheet and created by seasonal melting and refreezing events. Here, we model the icy inclusions as randomly oriented, discrete, noninteracting, dielectric cylinders embedded in a transparent snow medium. An exact analytical solution is used to compute the scattered field from the cylinders. Using this model, we correctly predict the polarimetric radar observations gathered by an airborne imaging system at three wavelengths (5.6, 24, and 68 cm), between 19 deg and 65 deg incidence angle. The diameter and number density of the cylinders that are inferred from the radar data using the model are consistent with in situ observations of the icy inclusions. The large radar reflectivity and polarization ratios are interpreted as arising internal reflections of the radar signals in the icy inclusions that first-order external reflection models fail to predict. The results compare favorably with predictions from the coherent backscatter or weak localization theory and may provide a complementary framework for interpreting exotic radar echoes from other planetary objects.

Rignot, Eric↗