Introduction

The cardiovascular system undergoes a constellation of important hemodynamic changes in microgravity—such as fluid shift of about 2 l from the legs to upper body, total blood volume reduction of around 11%, decrease in cardiac muscle mass of over 10%, and reduced exercise capability of up to 20% in the first 30 days—which lead to cardiovascular deconditioning1,2,3,4,5,6,7,8. Understanding 0G configuration in detail is particularly crucial to prevent and mitigate disabling effects through targeted countermeasures, and thus ensure the full health and well-being of astronauts during spaceflight missions2. A better comprehension of deconditioning mechanisms can answer open questions on readaptation for spaceflight beyond 1 year. It is currently unknown whether very lengthy missions lead to an amplification of reversible changes already known or the emergence of unrecognized and irreversible alterations in cardiovascular function3. Moreover, even if cardiovascular deconditioning is not extreme during spaceflight, it can later become hazardous at the time of reentry on the Earth or partial gravity restoration (e.g., Moon or Mars landing), when immediate demands for physical capacity are required and cannot be met. In view of future Moon or Mars missions, long-term cardiovascular deconditioning will need to be addressed to: (i) evaluate the efficacy of countermeasures and to optimize them in order to guarantee the astronauts’ safety; and (ii) establish risk thresholds, especially for physical activity after landing. Additionally, as space science is an extraordinary investment multiplier (9:1 return on investment by way of spin-offs from space technology), understanding space deconditioning can have impacts on other clinical conditions and on aging research3.

Although the overall scenario of cardiovascular deconditioning is widely accepted, the extreme difficulty of performing clinical measurements, the unsystematic data collection on older missions, the limited number of human space missions, and the heterogeneity of the missions (in terms of reference baseline posture, crew, objectives, countermeasures, duration, etc.) make the details of cardiovascular responses controversial and restricted to a modest number of cardiac parameters9,10,11,12. To the best of our knowledge, hemodynamic parameters related to the right heart, oxygen consumption, and exercise tolerance are barely known in microgravity. Almost no information at the vascular level regarding organ blood supply and wave propagation mechanisms is available.

Ground-based analogs—such as bed rest studies and water immersion—are not fully representative of microgravity as they introduce artificial horizontal and hydrostatic pressure gradients respectively. Parabolic flight temporarily reproduces 0G effects with each session lasting only 20–30 s13.

Given the above difficulties, the computational approach is a recent and promising tool to study the human body’s response to microgravity conditions14,15,16,17,18, as it can be used to investigate variables and processes which are hard to measure, as well as predict different scenarios and optimal countermeasures. In particular, reduced order modeling seems to be the most appropriate strategy, as it is a reasonable balance between the suitable level of local detail and the overall response (see, among the most recent19,20,21,22,23). Multiscale modeling has already provided useful insights into aging24,25, coronary pathologies26,27, and cardiac arrhythmias28,29,30.

The present study aims to computationally shed light on the cardiovascular knowledge gaps previously described, following a comprehensive bibliographic investigation of more than 50 studies of cardiovascular deconditioning during microgravity. Among the available literature we only focused on the most recent long-term spaceflight data. This allowed us to extract information to characterize the steady-state cardiovascular condition for long-term (at least 5 months) spaceflight. In order to model the typical long-term response, we prioritized studies having only routine and customary countermeasures, and avoided spaceflights with specific/peculiar ad hoc countermeasures. The main cardiovascular features of the 0G adaptation point are: blood shift from lower to upper body, total blood volume reduction, reduced cardiac function and volume, increase of lower extremity venous compliance, arterial resistance variation, and baroreflex response alteration. Our goal is to quantitatively describe in detail the cardiovascular response through a comprehensive 0D–1D multiscale model, which has been successfully tested by heart pacing and open-loop response31. The 0D–1D heart-arterial hemodynamics of the complete model has been previously validated to subject-specific and pathological conditions25,28,32,33,34. The physical-based model combines a 1D description of the arterial tree together with a lumped parameterization of the remaining regions, i.e., venous return, heart chambers, pulmonary circulation, and baroreceptor regulation (see the Methods for more detail and for the limitations of the present approach). In particular, we compared 1G supine and 0G spaceflight conditions in order to: (i) understand the underlying mechanisms leading to cardiovascular deconditioning, and (ii) describe the hemodynamics of zones for which clinical data are not yet feasible and accurate. The adoption of the 1G supine condition as a reference baseline has a twofold justification: (i) supine is the reference medical position defining the benchmark for physiological cardiovascular behavior in clinical practice; (ii) over recent decades, cardiovascular modeling (see, among others,19,20,21,23,35) has been derived which simulates the posture of primary medical interest and has developed a robust, data-based, and well-accepted parameter setting, which serves as a solid reference point for parameter tuning. Changes occurring in 0G condition described hereafter are always reported with respect to the (chosen) 1G supine baseline. The present study provides computational-based insights into several cardiac parameters and hemodynamic variables (in terms of time-series, average distribution throughout the body, and waveform alteration), which are currently unexplored.

Results

Results are organized showing the main cardiac parameters first. Whenever possible, to validate the modeling response we compared our results to available literature data of long-term spaceflights which employed routine countermeasures and which referred to the 1G supine condition. The cardiovascular system was divided into four areas: cerebral, cardio-thoracic, abdominal, and lower limbs. Special focus was given to the proximal-to-distal pathway of the arterial tree. Time-series, characteristic levels (i.e., mean, maximum, excursion) and waveform changes of pressures (P), flow rates (Q), and volumes (V), were then evaluated for representative sites of each cardiovascular region. In the following, supine 1G configuration on the Earth is depicted with blue curves, while long-term 0G spaceflight is depicted with red curves. For each configuration results are referred to a generic steady-state heartbeat, so that t [0, RR], where RR [s] is the cardiac beating period and HR = 60/RR [bpm] is the heart rate.

Cardiac parameters

The most important cardiac parameters are reported in Table 1 (see Supplementary Information for definitions). We first compared present outcomes with literature data. To this end, it should be noted that the detailed cardiovascular response in literature is quite discordant and debated. As comprehensively reported by Norsk12, for different hemodynamic parameters—such as heart rate, arterial pressure, cardiac output, stroke volume—there are conflicting results, which depend on the level of adopted countermeasures, the mission duration, and, most importantly, the baseline (upright/sitting/supine) reference 1G condition. In view of the above considerations very few assessed metrics were available for full comparison and, among these, we selected the most recent studies involving long-term spaceflights (around 5 months). Thus, collected measurements are the most reliable results which we adopted as reference values for the present study.

Table 1 Modeled cardiac parameters for the supine 1G configuration and during long-term 0G spaceflight (green color is adopted when spaceflight data are available, yellow otherwise).

A first validation comes from the left ventricular contractility parameters: they all reduced (−18.32% for SV, −9.30% for Vlved, −9.93% for EF) except for a slight increase of Vlves, +4.75%, and are in good agreement with spaceflight data4,5,36,37,38. A reduction of cardiac contractility is also in general accordance with cardiac atrophy observed during bed rest39,40,41. CO drop (−7.58%) is less evident than SV, due to the HR increase in the spaceflight configuration (+13%, see the Methods section). The CO decrease is in qualitative agreement with the variations observed with respect to 1G supine condition4,5,36,37,38. As previously mentioned, CO (together with SV) is one of the cardiac parameters whose variation in 0G greatly depends on the chosen 1G baseline configuration. As observed by ref. 12,42, in 0G CO increased with respect to the sitting 1G position. Compared to the supine 1G position, CO was usually observed to decrease4,5,36,37,38, although mild increases were also reported42. Mean central arterial pressure (MAP) was also tested for validation: MAP decreased by 9.92% in our model (as does brachial arterial pressure), falling within the range of measured data during long-term spaceflights4,6,43,44. Mild reduction of MAP with respect to seated preflight values was also observed in long-term spaceflight using a countermeasure protocol1,12. The green rows of Table 1 show that there is satisfactory agreement between currently available in vivo and our in silico outcomes. This confirms the correct 0G configuration setting and the model’s accuracy. On the other hand, the scarcity of previously measured parameters and the good predictive level of the model justifies its use to obtain cardiac estimates not previously measured directly. It should be recalled that these cardiac variables were not used to tune the model parameters, but to a posteriori verify the modeling response once the parameter setting was completed. The comparison with the reported literature represents the only validation (with data actually measured and related to the same simulated configuration) currently available in 0G conditions.

Other cardiac parameters were also modeled. A substantial decrease of the work done by the heart (−19.29% for SW/min) and a more limited reduction of oxygen consumption (−1.91% for RPP, −8.38% for TTI/min) were found. The slight decrease of myocardial oxygen consumption indexes represents an additional modeling validation, in agreement with the maximal oxygen uptake reduction observed in long-term spaceflight45. Central venous pressure (CVP), estimated as the mean value of the superior and inferior vena cava pressure signals averaged over the beat, decreased by 5.08%. This result is in qualitative agreement with the CVP reduction observed in early spaceflights12,46,47 and parabolic flights48,49, though it should be kept in mind that the time scales and physical mechanisms inducing CVP variations found in literature are not fully comparable to the present study. Central aortic pulse pressure, PPAA, decreased (−22.30%) due to a higher reduction of central aortic systolic pressure (−13.39%), with respect to the diastolic pressure decrement (−6.90%). In agreement with the pulse pressure behavior, central aortic augmentation index, AIAA, also declined (−41.18%).

The observed variations of cardiac work, oxygen consumption, and contractility indexes, as well as central mean and pulse pressures signal an overall scenario of cardiac deconditioning. The cardiovascular system reaches a condition commonly seen in a sedentary lifestyle. This is achieved initially by direct fluid dynamics changes and is then maintained through continued presence in the less demanding 0G environment.

Pressures, flow rates and volumes

Pressure, P(t), flow rate, Q(t), and volume, V(t), time-series were analyzed together with the relative variations of the beat-averaged values (\(\overline{P}\), \(\overline{Q}\), and \(\overline{V}\)) during spaceflight with respect to the supine 1G configuration on the Earth. In addition to mean variation, relative variation of pulse pressure, PP (defined, in the general case, as the difference between maximum and minimum pressure values), and maximum pressure values, Pmax, were evaluated. In addition to being general indicators of health and being mechanical properties of the cardiovascular system, these latter quantities estimate the vertical stretching of the pressure signals. Figure 1 shows representative results for the four macro-regions considered, while P and Q variations at other significant sites are listed in Tables 2 and 3.

Fig. 1: Pressures and flow rates throughout the body in 1G supine and 0G conditions.
figure 1

Time-series of the pressure, P(t), and flow rate, Q(t), at different sites: left vertebral artery (cerebral), right ventricle and left brachial artery (cardio-thoracic region), left renal artery (abdominal region), leg veins (lower limbs). Blue: supine 1G configuration on the Earth; red: 0G spaceflight configuration. Mean pressure (\(\overline{P}\)), mean flow rate (\(\overline{Q}\)), together with maximum pressure (\({P}_{\max }\)), and pulse pressure (PP) relative variations between 1G and 0G configurations are reported in each panel.

Table 2 Beat-averaged levels of pressure \(\overline{P}={\int}_{RR}P(t)dt/RR\) [mmHg] and flow rate \(\overline{Q}={\int}_{RR}Q(t)dt/RR\) [ml/s] at different significant sites in supine 1G and 0G conditions, with relative variations [%].
Table 3 Maximum pressure, \({P}_{\max }\), and pulse pressure, PP, values [mmHg] at different significant sites in supine 1G and 0G conditions, with relative variations [%].

Both pressure and flow rate mean levels decreased in the spaceflight configuration. Pressure variations were quite homogeneous throughout the different cardiovascular regions, ranging between −8 and −12% (apart from the cardiac region and venous return, where variations are lower, between −2 and −6%). Flow rate changes were more heterogeneous. In particular, variations of the mean values decreased from the upper to the lower body, from −16 to −18% in the cerebral region to −2 to −6% in the lower limbs (apart from the anterior tibial artery, −11.94%). Pressure maximum variations are more heterogeneous (around −5 to −15%) than mean levels, and they are both smaller in the venous return than in the arterial circulation. PP variations remarkably showed high variability (ranging from +4.82 to −56.00%), which is not directly attributable to the location, type (venous/arterial), and size of the region. An example is represented by the superior and inferior venae cavae (see Table 3), both veins ending in the right atrium and similar in size, but showing completely different PP variations (+4.82 and −32.97%).

Volume time-series, V(t), and relative variations of the average values per beat (\(\overline{V}\)) during spaceflight with respect to the supine 1G condition are reported in Fig. 2 for the 0D compartments. Mean pressure and flow rate values are strictly related to blood volume variations induced by the prolonged exposure to microgravity, since V = V0 + CP, where C and V0 are the compliance and unstressed volume of each district, respectively. In the lower limbs and abdominal region below VIP (VIP, see the Methods for more details), mean volume variations were widely negative. In the upper body, volume reduction was smaller, resulting in volume increase in the upper body and pulmonary regions. The different upper-to-lower body volume repartition primarily resulted from the blood shift, which in the upper body partially counterbalances the effects of volume reduction. The further cardiac volume and contractility reduction mitigated the fluid shift effects in the cardiac region. As a result, heart chamber volumes were reduced more than the thoracic surroundings.

Fig. 2: Volumes throughout the body in 1G supine and 0G conditions.
figure 2

Time-series of the volume, V(t), at different 0D sites: venules (upper body), left ventricle and pulmonary artery (cardio-thoracic region), veins below VIP (abdominal region), capillaries (lower limbs). Blue: supine 1G configuration on the Earth; red: 0G spaceflight configuration. Mean volume (\(\overline{V}\)) relative variations between 1G and 0G configurations are reported for each region (colored values refer to the corresponding time-series shown).

The time-series of the hemodynamic variables shown in Figs. 1 and 2 demonstrate the whole signal differences between 1G supine and 0G conditions. Signals were not merely translated or scaled but, also due to the different HR, their waveforms changed. This aspect is discussed in more detail in the next section.

Hemodynamic signals: waveform alteration

We focused on the waveform alteration of pressure and flow rate signals in terms of the normalized signal difference, NSD, which is a dimensionless measure of signaling similarity (ranging in the interval [0.2]) between the 1G supine and 0G spaceflight configurations (see the Methods section for more details). Representative results of each cardiovascular region are shown in Fig. 3, while NSD values for both \(P^{\prime}\) and \(Q^{\prime}\) at different sites can be found in Table 4.

Fig. 3: Pressure and flow rate waveform alterations.
figure 3

Normalized time-series \(P^{\prime} (t^{\prime} )\) and \(Q^{\prime} (t^{\prime} )\) at different sites: left internal carotid artery (cerebral), left ventricle and inferior vena cava (cardio-thoracic region), celiac artery (abdominal region), posterior tibial artery (lower limbs). Blue: supine 1G configuration on the Earth; red: 0G spaceflight configuration. NSD (normalized signal difference) values for \(Q^{\prime}\) and \(P^{\prime}\) are reported in each panel.

Table 4 NSD values at different significant sites for pressure and flow rate signals.

It is immediately apparent that the normalized signals \(P^{\prime}\) e \(Q^{\prime}\) in 0G (red curves) in all cardiovascular sites were delayed with respect to the 1G supine condition (blue curves), as a result of the horizontal stretching due to the faster RR beating rate in 0G. However, the delay was not merely due to HR variations, as it was not uniformly found in all the regions, but varied without showing a significant trend from upper to lower regions. This means that waveform alteration acts differently throughout the body and as well as at the same site. For instance, by comparing 0G and 1G, note the different shape and phase of the \(Q^{\prime}\) signals in the inferior vena cava, and the different oscillation amplitude for the carotid \(P^{\prime}\) and \(Q^{\prime}\) signals. NSD, which was introduced to quantify signal stretching besides pure horizontal translation, was generally not negligible and often beyond 0.5. Flow rates were more affected by shape variations than pressures, and proximal regions were less prone to waveform variations than distal regions (see Table 4).

Arterial proximal-to-distal tree

Special attention was given to the arterial proximal-to-distal tree, focusing on the arterial 1D system from the ascending aorta to the anterior tibial artery. The scenario emerging from Fig. 4 (and detailed in Tables 2 to 4) is multi-faceted and not always ascribable to the macro-regions crossed by the pathway. Mean flow rate variations were close to −9% in the ascending aorta, increasing to approximately −6% in the last part of the abdominal aorta, and then decreasing again towards the anterior tibial artery to around −12%. Mean pressure variations were instead quite constant, in the range of −10%, along the pathway, similarly to what happens for the Pmax (−11 to −13%). The overall proximal-to-distal trend of PP variations was a slight decrease, while NSD variations increased remarkably. In particular, NSD at least doubled for \(P^{\prime}\) and \(Q^{\prime}\), and \(Q^{\prime}\) variations were also larger than for \(P^{\prime}\). What was previously reported for the \(P^{\prime}\), \(Q^{\prime}\), and NSD behavior from proximal-to distal areas is also true here, where only arterial sites were considered.

Fig. 4: Proximal-to-distal arterial tree: pressures and flow rates in 1G supine and 0G conditions.
figure 4

Time-series P(t), Q(t), \(P^{\prime} (t^{\prime} )\), and \(Q^{\prime} (t^{\prime} )\) at the extremes of the 1D arterial proximal-to-distal tree: ascending aorta and anterior tibial artery. Blue: supine 1G configuration on the Earth; red: 0G spaceflight configuration. NSD values for \(Q^{\prime}\) and \(P^{\prime}\), together with mean pressure (\(\overline{P}\)), mean flow rate (\(\overline{Q}\)), max pressure (\({P}_{\max }\)), and pulse pressure (PP) relative variations between 1G and 0G configurations are reported in the panels.

Discussion

These results depict a complex picture where changes imposed by the 0G configuration created remarkably heterogenous stresses on the hemodynamic variables of the cardiovascular system. Volume reduction and baroreflex responses acted on global levels of pressure, volume, and flow rate, while blood shift, resistance and compliance variations altered the local mean levels of the hemodynamic variables throughout the body. The most important alteration in volume was the combination of blood shift and reduction: in the upper body the incoming blood volume partially compensated for the volume lost, resulting in either small reductions or positive variations (in particular, at the pulmonary and cerebral levels). In the lower body two contributions added up, leading to a more evident volume decrease. However, the upper body blood volume overload did not translate into pressure and flow rate increases due to the complex interplay between resistance and compliance variations within each region. Moreover, it should be kept in mind that (negative) variations were referred to the 1G supine condition, while in the cerebral circulation positive variations are expected if referred to the 1G standing position, as suggested by ref. 50. The increased venous compliance observed in the legs confined lower body mean pressure variations to a range which is similar to that of the upper body, that is −8 to −12% (see Fig. 1 and Table 2). The smallest variations (−5%) occurred in the capillary–venous districts of the legs. Cardiac compliance increased accounting for the reduced contractility, and this in turn limited pressure variations at the central level (−1 to −8%). If arterial resistance had not changed, following homogeneous pressure variations we would have expected homogeneous flow rate variations throughout the body, around −8 to −12%. Instead, all lower body resistances reduced, while carotid/vertebral resistances increased. Therefore, flow rates towards upper body and cerebral areas were reduced (−17 to −19%), while lower body flow rate was on average boosted (−2 to −5%, with the anterior tibial artery as only exception). The cardio-thoracic region experienced flow rate reductions in between lower and upper body values, namely −8 to −10%.

Mean flow rate variations decreased from ascending aorta to the inferior abdominal aorta, consistent with what happened from upper towards lower body, i.e., due to resistance changes flow rate was more hampered in the upper thoracic region with respect to the lower abdomen. Note also that ascending aorta variations were damped by the cardiac volume and contractility reduction. Mean flow rate variations then increased again towards the inner iliac-femoral-posterior tibial path, for two possible reasons: (i) conduction 1D arteries (e.g., femoral and external iliac arteries) are less influenced by resistance variations than their bifurcating branches (e.g., deep femoral and inner iliac arteries); (ii) on the Earth, resistance of the posterior tibial artery is higher than the anterior tibial. Equal relative reduction of resistance (−10%) for both arteries led to an absolute resistance reduction which is more important for the posterior tibial artery. In the latter, flow rate was thus more enhanced (i.e., lower mean Q variations) than in the anterior tibial artery (higher mean Q variations).

The hemodynamic and cardiac parameter variability reveals a pattern which, overall, resembles the spaceflight configuration of a relaxed, sedentary lifestyle. This is directly due to both fluid dynamics effects as well as reduced physical activity and lower cardiac work to counter the force of gravity3. Cardiac work, oxygen consumption, and contractility indexes were reduced, as were central mean and pulse pressures. Almost all the hemodynamic and cardiac parameters resulted from variables pertaining to the thoracic-cardiac region, for which the above considerations about P and Q mean values still hold. The hemodynamic mechanisms mainly responsible for the observed cardiac deconditioning are the global volume reduction together with the decrease of cardiac volume and contractility. The slight increase of Vlves can be interpreted as a possible sign of a reduced ability of the ventricle to empty and ventricular overload due to fluid shift (pulmonary artery and vein volumes increased on average by 14.61% and 21.39%, respectively).

Maximum and mean pressure value variations were smaller in the venous than in the arterial compartments (see Tables 2 and 3). This behavior is related to the higher compliance of the venous regions, which was further increased in the lower limbs during spaceflight. Venous vessels were thus more able to absorb and reduce pressure changes than their arterial counterparts. In general, PP variations were higher than mean pressure and Pmax changes: apart from one exception (+4.82%), pressure signals tended to reduce their ranges to a greater extent than their means. In general signals were less pulsatile (between −10 and −30% from 1G supine to 0G) and, in each region, the greatest PP variations (e.g., >20%) usually occurred in the distal areas, where PP on the Earth is already limited and among the lowest in that region: arterioles, venules and veins were thus the most affected areas. However, the contrary was not in general true: the smallest PP values on the Earth (e.g., right atrium and superior vena cava) did not necessarily imply the largest PP variations during spaceflight.

PP and NSD variations showed great variability throughout the body, meaning that signals were stretched differently in vertical and horizontal directions depending on the site. PP and NSD alterations were inverted along the 1D arterial tree (see Fig. 4 and Tables 3, 4). At one extreme (ascending aorta) PP variation was high and signal was more squeezed than in distal regions, but NSD variations were limited. At the other extreme (anterior tibial artery), pressure signal was stretched less but experienced a remarkable waveform change in terms of NSD.

A possible interpretation for the NSD increase towards the distal regions is that 0G-induced waveform alterations mainly happen peripherally, where the traveling wave is already modified in 1G by vessel tapering and non-linearities51, and the great majority of reflections occur52. In addition, since distal regions are most closely influenced by arterial resistance and venous compliance variations, more pronounced alterations can be expected there. Thus, in the proximal regions the interaction of many backward waves arriving at different times makes them dampen and partially elide each other. As a result, the total wave was shaped mainly by the forward contribution, and waveform variation induced by the 0G configuration was less evident (NSD < 0.2). When instead, in the peripheral zones, forward and backward waves are comparable due to the proximity to reflection sources (such as bifurcations, discontinuities, or arteriolar interface), NSD was up to or beyond 0.5. The backward wave being most influenced during the 0G configuration explains the greater NSD variations for Q than P. In fact, the flow rate backward wave is upside down with respect to the corresponding forward wave, while the pressure backward wave is reflected upright. Thus, the same backward wave variation (in absolute amplitude) has a stronger impact on the total flow rate wave rather than on pressure.

The NSD increase observed in the peripheral regions is likely due to reflection mechanisms within the traveling wave combined with the HR increase which, during spaceflight, exacerbated the waveform alteration. Moreover, due to the HR increase in 0G the systolic interval was longer (with respect to the RR period) than in 1G. Even if the imposed HR increase was not worrisome (+13%), diastole was shortened and diastolic function (e.g., ventricular filling) could be partly affected. As previously noted, waveform variations induced by reflection mechanisms were not absorbed downstream in the arterial tree, leading to evident signal variations (e.g., NSD values around 0.4–0.5) in arteriolar–capillary districts, leg veins, and venae cavae. It is reasonable to expect that these regions downstream of the arterial system do nothing but undergo and maintain the changes introduced upstream. The normal wave pattern was therefore modified and all the physiological phenomena at the capillary–venous level—such as regular perfusion, mean pressure per beat, and average nutrient supply at the cellular level—can experience important alterations.

Despite the fact that hemodynamic changes are important and diversified the observed variations in the steady-state 0G spaceflight configuration are not dramatic per se. For instance, although EF reduction was not negligible (−10%), myocardial contractility reduction is considered clinically relevant when EF drops 20% or more, which has not been observed to date during any spaceflight53. After long-term spaceflight a 0G adaptation point imposed by a less demanding environment is indeed reached by the cardiovascular system. The scenario becomes potentially hazardous at the time of reentry on the Earth or partial gravity restoration (e.g., Moon/Mars landing). Leaving the 0G environment imparts a backward blood shift from upper-to-lower body. Recalling that mean pressure at central, cerebral, and organ levels was decreased by 10%, backward blood shift combined with an overall hypovolemic condition (−11.5% for the total blood volume) can impair normal perfusion above VIP. Moreover, compliance increase and resistance decrease in the lower limbs amplify the venous pooling already present during reentry and landing. All these mechanisms converge to produce the well-known condition of orthostatic intolerance, which is further aggravated by the reduced and slowed baroreceptor response. In fact, as baroreceptors are not able to promptly control the vasculature and regulate chronotropic and inotropic effects, a normotensive state is hardly guaranteed3.

Exercise capability was also markedly reduced. A hypothetical Moon or Mars landing (with partial gravity restored) requires a high work capacity and exercise tolerance, which were instead damped with respect to preflight conditions. The heart got used to a reduced pumping requirement and reduced work, with subsequent need for less oxygen. The less demanding spaceflight environment caused a central PP reduction of 22.30%. In 1G conditions PP is lower for sedentary than resistance-trained people54,55,56,57, thus PP can be taken as a marker of training level and physical activity. Gravity is itself a natural source of loading for the human body: the lower spine, lower extremities, and postural muscles are loaded up to 4–5 times the body weight simply walking on the Earth2. Therefore, in 1G, strength training is naturally provided by the gravity vector on the body. During spaceflight this natural load is absent and, without any particular countermeasure program, the exercise tolerance of a spaceflight traveler becomes comparable to that of an untrained person with a sedentary lifestyle.

All these aspects, along with other spaceflight-induced effects such as muscular atrophy and bone demineralization, need to be properly addressed especially in long-term human space missions where prompt physical capability after partial or complete gravity restoration is required. The outcomes of the present study highlighted a wealth of cardiovascular information otherwise not yet available in vivo, showing how computational hemodynamics is particularly valuable in environmental conditions such as microgravity and spaceflight where even basic clinical measures are rarely available. It should be also noted that given identification (sex and age) and anthropometric parameters (such as weight and height), together with basic cardiovascular information (such as brachial cuff pressure), the model can be made patient-specific. That would provide good estimates of cardiovascular fitness related to the heterogeneity of the crew, similar to what was already done in 1G physiological and pathological conditions32,34,58. Cardiovascular deconditioning without ad hoc countermeasures, as modeled here, is of fundamental importance to have a clear and valid baseline for testing against the cardiovascular response with specific countermeasures. By comparing and contrasting the two conditions the efficacy of a single or a combination of countermeasures can be extensively evaluated. A future step might include the effect of one or more countermeasures in the model to observe the response to and help design, implement, and improve such countermeasures. Future modeling developments can also estimate the grade of orthostatic intolerance and exercise tolerance by simulating the return to a gravitational environment, starting from the cardiovascular deconditioned configuration here described. Eventually, understanding the mechanisms of cardiovascular deconditioning during spaceflight will have a significant impact on the knowledge of aging physiology on the Earth3. Since deconditioning in spaceflight by gravity deprivation is strongly analogous to deconditioning on the Earth by gravity withdrawal, as in sedentary aging, long-term spaceflight physiology can offer precious hints for delaying or preventing modern lifestyle medical disorders related to increased longevity.

Methods

Mathematical model

The present multiscale modeling of the cardiovascular system combines a 1D description of the arterial tree together with a lumped parameterization of the remaining regions, that is, venous return, heart chambers, pulmonary circulation and baroreceptor regulation (see Fig. 5). The model has been tested using heart pacing and open-loop response, resulting in good agreement with measured data31. The hemodynamics of large-to-medium sized arteries relies on the model first proposed and validated by Guala et al.32, which has been employed to study aging24,25 and atrial fibrillation effects28,29.

Fig. 5: Scheme of the multiscale model.
figure 5

The sketch includes the arterial tree (right), the structure of the 0D compartments (left), the arteriolar 0D–1D interface (dashed box), and the baroreceptor mechanisms (top). Details of the 0D compartments are depicted in the bottom panels. According to VIP, the model is divided into upper body (RH right heart, PC pulmonary circulation, LH left heart, U ABD and HA upper body arteries and arterioles; upper abdomen, head and arms venous return, SVC superior vena cava,) and lower body (L ABD and LEGS lower body arteries and arterioles; lower abdomen and legs venous return, IVC inferior vena cava). Modeling parameters and variables are defined in the Supplementary Information, see Supplementary Table 1 for the legend.

Mass and momentum balance equations were expressed in one-dimensional form, considering dissipation and inertial terms, and were combined with a constitutive equation accounting for the nonlinear viscoelastic behavior of the tapered arterial walls. Arterioles, capillaries, venous return (venules, veins, and venae cavae), cardiac, and pulmonary (arterial and venous) circulations were modeled through RLC 0D compartments, each characterizing local hemodynamics. The RLC network accounts for viscous (resistance, R), inertial (inductance, L), and distensibility/contractility (compliance, C, or elastance, E) effects. All four cardiac chambers were contractile, as modeled through time-varying elastance functions, including the atrial kick feature. Valvular dynamics were well captured, dominated by the main forces acting on the valves, such as pressure difference across the valve, frictional effects from neighboring tissue resistance, the dynamic motion effect of the blood acting on the valve leaflet and the action of the vortex downstream of the valve. The 1D → 0D interface was at the arteriolar level, with a characteristic impedance per artery put before the corresponding arteriole. The 0D → 1D interface was at the aortic valve level. The position of the volume indifference point (VIP, see the Supplementary Information for more details) split the model into upper and lower body regions (see Fig. 5). The short-term baroregulation mechanisms were also included in the model. The chronotropic and inotropic effects of both ventricles, as well as the control of the systemic vasculature (peripheral arterial resistances, unstressed volumes of the venous system, and venous compliances), were taken into account.

Simulations were run until the steady-state solution was achieved: for both 1G supine and 0G conditions the transient dynamics were completely extinguished after 100 cardiac periods. Thus, results are referred to the generic steady-state RR beat. Details of the governing equations and model parameters are offered in the Supplementary Information.

Definition and nomenclature for the common cardiac parameters are provided in the Supplementary Information. We defined a new metric in the present study, NSD. For the waveform analysis of signals, we exploited the dimensionless time \(t^{\prime} =t/RR\) (\(t^{\prime} \in\) [0, 1]), where RR [s] is the cardiac beating period. We normalized the time-series as \(y^{\prime} (t^{\prime} )=(y(t^{\prime} )-{\mu }_{y})/{\sigma }_{y}\), where μ and σ are the mean and standard deviation values of y(t) over the cardiac beat, respectively (y(t) = P(t) or Q(t)). We then defined the NSD averaged per heartbeat as

$${\rm{NSD}}=\mathop{\int}\nolimits_{0}^{1}| y^{\prime} {(t^{\prime} )}_{1{\rm{G}}}-y^{\prime} {(t^{\prime} )}_{0{\rm{G}}}| dt^{\prime} ,\quad {\rm{NSD}}\in [0,2],$$

where \(y^{\prime} (t^{\prime} )\) can be either \(Q^{\prime} (t^{\prime} )\) or \(P^{\prime} (t^{\prime} )\) (having zero mean and unitary standard deviation values), and the subscripts 1G and 0G refer to the supine 1G and 0G spaceflight configurations, respectively. As hemodynamic time-series \(P^{\prime} (t^{\prime} )\) and \(Q^{\prime} (t^{\prime} )\) were normalized in vertical (amplitude) and horizontal (time) directions, signals in 1G and 0G can be fully compared: NSD quantitatively accounts for the shape and net waveform variation between them. The lower bound value (NSD = 0) is reached when the two signals overlap. Thus, if NSD is close to 0, no appreciable waveform variation is detected. On the contrary, the upper bound value (NSD = 2) can be reached only if the two signals are step functions (centered in \(t^{\prime}\) = 0.5 and \(y^{\prime}\) = 0 and with half-width equal to 1), one complementary to the other. The closer NSD is to 2, the more significant is the alteration of the waveform.

Spaceflight setting: criteria, cardiovascular mechanisms, and adopted changes

The setting reproducing simulated long-term (at least 5 months) spaceflight conditions without ad hoc countermeasures relies on an extensive bibliographic investigation of more than 50 studies (listed in the Supplementary Information) of human cardiovascular changes induced by microgravity. A general overview of the adopted hemodynamic 0G setting is sketched, while a detailed description is offered in the Supplementary Information. Among hemodynamic data available in literature we focused on those measured during spaceflights rather than ground-based experiments. When no useful information during actual spaceflight was available—such as for blood shift—we considered parabolic flight data, which is more informative since blood shift is extinguished in few hours. Within the subset of spaceflight data priority has been given to long-term and recent space missions. By analyzing the selected literature we based the model setting and validation primarily on long-term data having routine and involuntary countermeasures, as well as countermeasures not expressly declared. We did not directly use (for setting or validation) studies adopting particular ad hoc countermeasures, but accounted for them only in order to better contextualize and qualitatively compare our results. Despite the great caution in collecting and managing spaceflight information, it should be kept in mind that: (i) spaceflight literature data often reports variations with respect to the 1G upright condition, while the present model has been derived in 1G supine condition for clinical interest. For this reason, when possible, further estimates were done to convert 1G upright condition measures to 1G supine conditions; (ii) older data are usually without any countermeasure program while the more recent ones adopt a wide variety of different countermeasures. Although more recent studies were preferred, we took advantage of older measurements (usually without countermeasures) to properly weight the possible effects of countermeasures; (iii) there are overall few long-term missions and data are not always in agreement (usually because of a different type of experimental setting and reference position). This last aspect made the spaceflight setting even more challenging, since a significant part of the analyzed data has been considered not sufficiently reliable for present purposes.

The main cardiovascular changes due to weightlessness introduced into the model were related to blood shift, blood volume reduction, leg venous compliance, cardiac function, arterial resistance, and baroreflex response. Table 5 summarizes the mechanisms considered, the resulting cardiovascular changes and details, where necessary, about their modeling implementation. For blood shift, once VIP was located, blood volume was shifted from lower to upper body using parabolic flight data and based on both the vessel distensibility and the distance from VIP. Total blood volume was then reduced by 11.5%, in accordance with long-term spaceflight data. To account for the contractility and volume reductions ventricular elastances were modified and cardiac volume was decreased. Leg venous compliance was increased by 27%, while arterial resistances were changed according to the body region (+10% for vertebral and carotid arteries; −10% for the lower body), as observed after long exposure to weightlessness. For baroreflex mechanisms only baseline values were altered as we were interested in long-term spaceflight in steady-state conditions: HR was increased by 13% and average aortic-carotid sinus pressure was decreased by 10%. An extended description of the hemodynamic considerations and the resulting estimated variations for each cardiovascular mechanism is detailed in the Supplementary Information, together with the corresponding bibliographic references supporting the adopted changes.

Table 5 Summary of the macroscopic mechanisms considered, the adopted cardiovascular changes with respect to the 1G supine condition and details of the corresponding modeling implementation.

Limitations

A limiting aspect of the present model is that it does not directly control interstitial fluid shift, since we can only act directly on blood volume and blood shift. Nevertheless, interstitial fluid shift and other long-term mechanisms—such as muscle atrophy and blood volume reduction—were intrinsically accounted for in the setting of the long-term spaceflight configuration by modifying the total and unstressed volumes of all cardiovascular compartments. In this regard, effects—such as net diuresis and fluid re-equilibration, fluid shift from intravascular to interstitial space, and the possibility of reduced thirst—were not individually accounted for, but overall were included in the blood volume reduction mechanism.

The model currently does not explicitly include a gravity force term since its effects were modeled based on the available microgravity literature. A gravity force term can be added to the model equations and adjusted to appropriately reproduce lunar, Martian, or other gravity conditions. At present this use would be purely predictive but very helpful in planning future interplanetary missions and extraterrestrial human occupancy.

Reporting summary

Further information on experimental design is available in the Nature Research Reporting Summary linked to this article.