Restructuring and serving web-accessible streamflow data from the NOAA National Water Model historic simulations.

In 2016, the National Oceanic and Atmospheric Administration deployed the first iteration of an operational National Water Model (NWM) to forecast the water cycle in the continental United States. With many versions, an hourly, multi-decadal historic simulation is made available to the public. In all released to date, the files containing simulated streamflow contain a snapshot of model conditions across the entire domain for a single timestep which makes accessing  time series a technical and resource-intensive challenge. In the most recent release, extracting a complete streamflow time series for a single location requires managing 367,920 files (~16.2 TB). In this work we describe a reproducable process for restructuring a sequential set of NWM steamflow files for efficient time series access and provide restructured datasets for versions 1.2 (1993-2018), 2.0 (1993-2020), and 2.1 (1979-2022). These datasets have been made accessible via an OPeNDAP enabled THREDDS data server for public use and a brief analysis highlights the latest version of the model should not be assumed best for all locations. Laslty, we describe an R package that expedites data retrieval with examples for multiple use-cases.

Social vulnerability of the people exposed to wildfires in U.S. West Coast states.

Understanding of the vulnerability of populations exposed to wildfires is limited. We used an index from the U.S. Centers for Disease Control and Prevention to assess the social vulnerability of populations exposed to wildfire from 2000-2021 in California, Oregon, and Washington, which accounted for 90% of exposures in the western United States. The number of people exposed to fire from 2000-2010 to 2011-2021 increased substantially, with the largest increase, nearly 250%, for people with high social vulnerability. In Oregon and Washington, a higher percentage of exposed people were highly vulnerable (>40%) than in California (~8%). Increased social vulnerability of populations in burned areas was the primary contributor to increased exposure of the highly vulnerable in California, whereas encroachment of wildfires on vulnerable populations was the primary contributor in Oregon and Washington. Our results emphasize the importance of integrating the vulnerability of at-risk populations in wildfire mitigation and adaptation plans.

A novel approach to 3D bone creation in minutes : 3D ultrasound.

AIMS: The objective of this study is to assess the use of ultrasound (US) as a radiation-free imaging modality to reconstruct 3D anatomy of the knee for use in preoperative templating in knee arthroplasty. METHODS: Using an US system, which is fitted with an electromagnetic (EM) tracker that is integrated into the US probe, allows 3D tracking of the probe, femur, and tibia. The raw US radiofrequency (RF) signals are acquired and, using real-time signal processing, bone boundaries are extracted. Bone boundaries and the tracking information are fused in a 3D point cloud for the femur and tibia. Using a statistical shaping model, the patient-specific surface is reconstructed by optimizing bone geometry to match the point clouds. An accuracy analysis was conducted for 17 cadavers by comparing the 3D US models with those created using CT. US scans from 15 users were compared in order to examine the effect of operator variability on the output. RESULTS: The results revealed that the US bone models were accurate compared with the CT models (root mean squared error (RM)S: femur, 1.07 mm (SD 0.15); tibia, 1.02 mm (SD 0.13). Additionally, femoral landmarking proved to be accurate (transepicondylar axis: 1.07° (SD 0.65°); posterior condylar axis: 0.73° (SD 0.41°); distal condylar axis: 0.96° (SD 0.89°); medial anteroposterior (AP): 1.22 mm (SD 0.69); lateral AP: 1.21 mm (SD 1.02)). Tibial landmarking errors were slightly higher (posterior slope axis: 1.92° (SD 1.31°); and tubercle axis: 1.91° (SD 1.24°)). For implant sizing, 90% of the femora and 60% of the tibiae were sized correctly, while the remainder were only one size different from the required implant size. No difference was observed between moderate and skilled users. CONCLUSION: The 3D US bone models were proven to be closely matched compared with CT and suitable for preoperative planning. The 3D US is radiation-free and offers numerous clinical opportunities for bone visualization rapidly during clinic visits, to enable preoperative planning with implant sizing. There is potential to extend its application to 3D dynamic ligament balancing, and intraoperative registration for use with robots and navigation systems. Cite this article: Bone Joint J 2021;103-B(6 Supple A):81-86.

Three-dimensional analysis of the tibial resection plane relative to the arthritic tibial plateau in total knee arthroplasty.

BACKGROUND: Kinematically aligned total knee arthroplasty strives to correct the arthritic deformity by restoring the native tibial joint line. However, the precision of such surgical correction needs to be quantified in order to reduce recuts of the resection and to design assisting instrumentation. This study describes a method for novel three-dimensional analysis of tibial resection parameters in total knee arthroplasty. Pre-operative versus post-operative differences in the slopes of the varus-valgus and flexion-extension planes and the proximal-distal level between the tibia resection and the arthritic tibial joint line can reliably be measured using the three-dimensional models of the tibia and fibula. This work uses the proposed comparison method to determine the parameters for resecting the tibia in kinematically aligned total knee arthroplasty. METHODS: Three-dimensional shape registration was performed between arthritic surface models segmented from pre-operative magnetic resonance imaging scans and resected surface models segmented from post-operative computed tomography scans. Mean, standard deviation and 95% confidence intervals were determined for all measurements.  RESULTS: Results indicate that kinematically aligned total knee arthroplasty consistently corrects the varus deformity and restores the slope of the flexion-extension plane and the proximal-distal level of the arthritic tibial joint line. The slope of the varus-valgus plane is most precisely associated with the overall arthritic slope after approximately 3° of correction and the posterior slope is biased towards the overall arthritic plateau, though less precisely than the varus correlation. CONCLUSIONS: Use of this analysis on a larger population can quantify the effectiveness of the tibial resection for correcting pathologies, potentially reduce imprecisions in the surgical technique, and enable development of instrumentation that reduces the risk of resection recuts. The kinematic alignment technique consistently corrects varus deformities.

Clinical and statistical correlation of various lumbar pathological conditions.

Current clinical evaluations often rely on static anatomic imaging modalities for diagnosis of mechanical low back pain, which provide anatomic snapshots and a surrogate analysis of a functional disease. Three dimensional in vivo motion is available with the use of digital fluoroscopy, which was used to capture kinematic data of the lumbar spine in order to identify coefficients of motion that may assist the physician in differentiating patient pathology. Forty patients distributed among 4 classes of lumbar degeneration, from healthy to degenerative, underwent CT, MRI, and digital x-ray fluoroscopy. Each patient underwent diagnosis by a neurosurgeon. Fluoroscopy was taken as the patient performed lateral bending (LB), axial rotation (AR) and flexion-extension (FE). Patient specific models were registered with the fluoroscopy images to obtain in vivo kinematic data. Motion coefficients, C(LB), C(AR), C(FE), were calculated as the ratio of in-plane motion to total out-of-plane motion. Range of motion (ROM) was calculated about the axis of motion for each exercise. Inter- and Intra- group statistics were examined for each coefficient and a flexible Bayesian classifier was used to differentiate patients with degeneration. The motion coefficients C(LB) and C(FE) were significantly different (p<0.05) in 4 of 6 group comparisons. In plane motion, ROM(LB), was significantly different in only 1 of 6 group comparisons. The classifier achieved 95% sensitivity and specificity using (C(FE), C(LB), ROM(LB)) as input features, and 40% specificity and 80% sensitivity using ROM variables. The new coefficients were better correlated with patient pathology than ROM measures. The coefficients suggest a relationship between pathology and measured motion which has not been reported previously.

Salt marsh ecosystem biogeochemical responses to nutrient enrichment: a paired 15N tracer study.

We compared processing and fate of dissolved NO3- in two New England salt marsh ecosystems, one receiving natural flood tide concentrations of approximately 1-4 micromol NO3-/ L and the other receiving experimentally fertilized flood tides containing approximately 70-100 micromol NO3-/ L. We conducted simultaneous 15NO3- (isotope) tracer additions from 23 to 28 July 2005 in the reference (8.4 ha) and fertilized (12.4 ha) systems to compare N dynamics and fate. Two full tidal cycles were intensively studied during the paired tracer additions. Resulting mass balances showed that essentially 100% (0.48-0.61 mol NO3-N.ha(-1).h(-1)) of incoming NO3- was assimilated, dissimilated, sorbed, or sedimented (processed) within a few hours in the reference system when NO3- concentrations were 1.3-1.8 micromol/L. In contrast, only 50-60% of incoming NO3- was processed in the fertilized system when NO3- concentrations were 84-96 micromol/L; the remainder was exported in ebb tidewater. Gross NO3- processing was approximately 40 times higher in the fertilized system at 19.34-24.67 mol NO3-N.ha(-1).h(-1). Dissimilatory nitrate reduction to ammonium was evident in both systems during the first 48 h of the tracer additions but <1% of incoming 15NO3- was exported as 15NH4+. Nitrification rates calculated by 15NO3- dilution were 6.05 and 4.46 mol.ha(-1).h(-1) in the fertilized system but could not be accurately calculated in the reference system due to rapid (<4 h) NO3- turnover. Over the five-day paired tracer addition, sediments sequestered a small fraction of incoming NO3-, although the efficiency of sequestration was 3.8% in the reference system and 0.7% in the fertilized system. Gross sediment N sequestration rates were similar at 13.5 and 12.6 mol.ha(-1).d(-1), respectively. Macrophyte NO3- uptake efficiency, based on tracer incorporation in aboveground tissues, was considerably higher in the reference system (16.8%) than the fertilized system (2.6%), although bulk uptake of NO3- by plants was lower in the reference system (1.75 mol NO3-.ha(-1).d(-1)) than the fertilized system (approximately 10 mol NO3-.ha(-1).d(-1)). Nitrogen processing efficiency decreased with NO3- load in all pools, suggesting that the nutrient processing capacity of the marsh ecosystem was exceeded in the fertilized marsh.