\section{Introduction}
Machine learning techniques, including deep learning, often considered a subset of artificial intelligence (AI) technologies, are widely used in many research and business fields nowadays, especially in combinations with so-called ``big data,'' including satellite remote sensing data. Deep learning methods have rich expressive power comparable to or even superior to human cognitive abilities in many tasks\cite{lecun2015}.
Several applications of deep learning techniques using satellite data have been proposed, including astronomy, planetary sciences, and earth observations, and research efficiency and quality improvements are expected
\cite{ishikawa2021, cheng2020, zhang2016, kussul2017}.

In contrast, the inference process of deep learning models is often regarded as a black box, owing to the vast number of parameters and computational complexities.
This might obstruct the decision of production implementations of AI systems. Explainable artificial intelligence (XAI) is a field in AI research that explains and interprets machine-learning models to avoid this situation\cite{barredoarrieta2020, molnar2022}.
Discussions have also begun on application of XAI techniques to Earth observation \cite{gevaert2022, kakogeourgiou2021}.
When applying XAI methods to real-world tasks, we must be careful to clarify the purpose of the explanation.

Local Interpretable Model Explanations (LIME)\cite{tulioribeiro2016} and SHapley Additive exPlanations (SHAP)\cite{lundberg2017} are representative XAI methods for importance analyses of explanatory variables.
LIME explains models using a local linear approximation and outputs values indicating the contribution of each explanatory variable to a given inference.
SHAP also outputs the importance of each explanatory variable for the inference using Shapley values based on game theory.
For image data, Gradient-weighted Class Activation Mapping (Grad-CAM)\cite{selvaraju2016} is a widely used XAI method that shows areas of importance in an input image for each inference of a convolutional neural network (CNN) model by calculating the gradient of the last convolutional layer in the model. This method can also be categorized as an importance analysis XAI method.
These methods are valuable both for explaining the behavior of deep learning models to stakeholders and improving performance, while results from different methods are not always consistent.

In tasks where feedback from information obtained by the XAI is valuable, importance analysis methods for explanatory variables are useful, such as in cases of risk management\cite{wen2021} and quality control in industries\cite{senonar2021}.
Using the importance analysis methods, we can investigate which input data are significant for each inference and whether those inputs affect the results positively or negatively.
Following these explanations, we can control the input variables to avoid risks or improve the quality.
Additionally, risk or quality control of image data is possible by identifying and removing spots related to undesirable conditions using methods such as Grad-CAM.

Despite their advantages, we should note that these XAI methods only provide explanations or interpretations of inferences or models, not evidence or rationale for inferences.
Because the XAI explanation does not show the inference process, it is incorrect to conclude that the important explanatory variables are the reason why the model is inferred. Variable importance is calculated using approximations (e.g., LIME) or counterfactual assumptions (e.g., SHAP). Deep learning and other state-of-the-art models can reproduce complicated nonlinear relations, and they should not be undermined by the trivialization of the models. In this context, we point out several risks of utilizing XAI methods and suggest using interpretation-friendly models instead\cite{rudin2018}. We must understand what is actually calculated and what is not obtained from XAI methods.

The motivation for using XAI methods is not only for feedback to control the input data, but also to verify whether the inference results are reliable. To meet this demand, importance analyses of explanatory variables are insufficient because knowing which explanatory variables are important does not clarify how the model determines those variables to be important. In other words, we may have questions, such as \textit{why predict this when this variable is high.} In thfis view, one of the possible approaches to answer such questions is to show the data trained by the model. Compared with simplifying the model, it is more accurate to explain that the model is trained using this dataset and find an abstract pattern for inferences. This approach is related to a ``data-centric'' perspective of machine learning applications. The training dataset is an important factor in machine learning model performance along with model architectures and training algorithms, and the data-centric approach focuses on the dataset, unlike the conventional model-centric approach.
It would be valuable to show a summary of the training dataset, such as the number of instances and statistics for the explanatory and objective variables. When and how the data are collected could also be important reference information to verify whether we can assume that the training and input data to be inferred in the production release are from the same population.

Showing an example of instances from a training dataset would be a beneficial data-centric explanation of the model inferences \cite{molnar2022}. In particular, an example similar to the inferred target data would be a valuable reference to verify whether the model is trained on a dataset sufficiently close to the target data\cite{hanawa2021}. If the most similar example is too far away from the target, we can interpret that it is invalid to apply the model to the target data, and the inference result is suspicious.

In this study, we present a data-centric XAI method that shows an instance in the training dataset similar to the target data to be inferred by a deep learning model. We demonstrate our method using a remote-sensing image classification task and discuss effective cases.

\section{Methodology}
We present a method
named What I Know (WIK)
 for evaluating similarities in the data instances using the model outputs and calculated values of the last layer before the output layer in a deep-learning model. Although there are several metrics to quantify data similarities, such as norms and cosine similarity, our approach is not to measure similarities of input data themselves but to evaluate similarity after data abstraction by the model.
Therefore, WIK's results are model-dependent, and it is possible that the most similar data instances are not the same when we evaluate similarities with the different models. This corresponds to a similarity evaluation in the context of a task for which the model is trained. We will show the difference between the traditional similarity measurement and our context-based method by providing examples in the discussion section.

In this method, the similarity is measured by the $L^2$-norm between the vectors of the last layer before the output layer as
\begin{equation}
\sqrt{\sum _{i=1} ^{n} \left( f_i - g_i \right)^2}
\end{equation}
where $f_i$ and $g_i$ are the i-th components of the vectors and $n$ is the dimensionality of the vectors,
 by considering the part of the model before the last layer as a feature extractor.
 Figure~\ref{fig_concept} shows the concept of an idea that a deep learning model consists of a feature extractor and a simple linear classifier or regressor. In the layer before the output layer, features of the input instance for the task the model trained for are extracted, and we can compare input data instances using those features.
As shown in Fig.~\ref{fig_concept}, this concept is applicable to CNN models, in addition to fully connected networks.
 The idea of using the output of the layer before the last layer is used for some applications, such as the Fr\'echet Inception Distance \cite{heusel2017} for measuring the similarity of data distributions mainly for evaluating generative networks. We investigated several similarity metrics in a multidimensional space, such as cosine similarity, and concluded that the $L^2$-norm provides the most reasonable results. In addition, the $L^2$-norm evaluation achieved the best result in the feature similarity measurements compared by Hanawa et al.\cite{hanawa2021}.
We also provide results with other similarity metrics in the next section.
\begin{figure*}[!t]
\centering
\includegraphics[width=5.4in]{concepts_features.pdf}
\caption{Concept of the interpretation that the neural network consists of the feature extractor and a simple linear classifier/regressor for the extracted features,
for (a) a fully connected network and (b) a CNN.
}
\label{fig_concept}
\end{figure*}

The similarity evaluation with the last layer before the output layer does not use the weights of the connections to the output layers, and we applied a two-step similarity evaluation
as shown in Fig.~\ref{fig_flow}.
These weights provide important information for deriving inference results and are necessary for understanding how the model interprets the extracted features. The first step of the similarity evaluation is verifying the inference results, which the model concludes. Here, we focus on classification tasks and set the criterion that similar instances should be inferred as an identical class. In our future work, we will present a method for regression tasks. In the second step, we measure the $L^2$-norm of the vectors from the last layer before the output layer. A method of similarity measurement with the $L^2$-norm of the layer before the output has already been shown in a previous study\cite{hanawa2021}. In that work, it was pointed out that data instances that are evaluated as similar are not always classified as the same class. While it is not evident that we must avoid the mismatch for similarity evaluation XAI methods, such a situation does not occur in our method, owing to step 1.
\begin{figure*}[!t]
\centering
\includegraphics[width=5.4in]{flow_method.pdf}
\caption{
Flow of selecting the data instance most similar to the target data among the training data.
}
\label{fig_flow}
\end{figure*}

A similarity value calculated by WIK can also be output to evaluate the distances between instances in the feature space.
It would be useful to compare similarity values for different instances in the test dataset to understand the difference between the training and test datasets.
The similarity values are also useful to investigate how a model update affect the inference of the target instances.
Note that the similarity values of models trained for different tasks should not be compared because the scales of the similarity values are model-dependent.

To demonstrate the WIK method shown above, we use an example of remote sensing image classification using the EuroSAT dataset with the images taken by the Sentinel-2 satellite\cite{helber2017}. EuroSAT is a dataset with multi-spectral images in a 64$\times$64 pixel format and is labeled as ten classes of scenes, such as forest, river, industrial, and residential. EuroSAT has two types of datasets with difference sets of wavelength bands. Here, we use the EuroSAT visible light observation dataset with red, green, and blue color channels. We constructed a standard LeNet-type CNN with three sets of two convolutional layers and a pooling layer in sequence, followed by fully connected layers. The output layer had ten nodes corresponding to 10 classes, and the last layer before the output layer had 32 nodes.

We split the visible light EuroSAT dataset into 80\% and 20\% for training and testing, respectively, and trained the model with the training dataset.
The data were randomly divided between the training and testing sets, and both were confirmed to be class-balanced.
We evaluated the model using the test dataset and obtained an accuracy of 95.7\%. While this performance is not as good as the state-of-the-art results shown in a recent study\cite{gomez2021}, the purpose of the demonstration in this article is to show an example of the usage of our XAI method for real-world applications. This is a realistic example of a model that can be constructed in finite time.
While it is true that models trained by fine-tuning or transfer learning from a model pre-trained by a large dataset are often used in practical applications, we selected a simpler model structure to highlight the XAI method presented here.

\section{Results}
We applied the example-based XAI method WIK presented in the previous section to the trained model. Examples of the input data instances to be inferred from the test dataset and those similar to the input data selected from the training dataset by the method are shown in Fig.~\ref{fig_result1}. The six images in the first row show the input data from the test dataset,
and the images from the second row to the fifth row are images from the training dataset that the model considered being similar with different similarity metrics, including cosine similarity
\begin{equation}
\frac{\vec{f} \cdot \vec{g}}{\|\vec{f}\| \| \vec{g}\|},
\end{equation}
$L^1$-norm
\begin{equation}
\sum _{i=1} ^{n} | f_i - g_i |,
\end{equation}
 $L^2$-norm, and maximum norm
 \begin{equation}
\max _{i} | f_i - g_i |.
\end{equation}
Apparently, the corresponding image pairs of the target data and the $L^2$-norm result (fourth row)
look similar to structures of similar spatial scales, colors, and textures, although the positional relationship is not exactly the same. The highway images show almost straight highways with similar widths surrounded by rural scenery, and the rivers in the river images are similar. The industrial images have several similar buildings with similar sizes, shapes, and colors, and the same relationship can be seen in the residential images. While the contrast of the forest images is different to some extent, the forest and permanent crop image pairs have similar spatial scales and colors. This result indicates that abstracted features, similar to text describing an image, are successfully extracted by the extractor part of the model, although it is difficult to evaluate such similarities using a predetermined mathematical formula. Owing to step 1 of our method, the predicted labels for the selected images are the same as those of the original images. Among the same-label instances in the training dataset, the selected images have the most similar features.
\begin{figure*}[!t]
\centering
\includegraphics[width=5.4in]{result1.pdf}
\caption{
Selected images by the presented XAI method with different similarity measurements (second to fifth rows) for the input images (first row) and images with the smallest $L^2$-norm of the pixel values (sixth row).
A label is shown for each image. All the predicted labels in the images shown here are correct. In the proposed method, the label of the selected image is necessarily the same as that of the original image.
}
\label{fig_result1}
\end{figure*}

A comparison of the results with the different similarity metrics shown in Fig.~\ref{fig_result1}, particularly for the river and industrial images, reveals that t$L^2$-norm provides the most reasonable result.
The widths and curvatures of the river images are similar to those of the target image for the $L^2$- and maximum norm results.
The sizes, shapes, and colors of the buildings in the industrial image with the $L^2$-norm are similar to those in the target data compared to the results with the other metrics.
While 3 out of 6 images are the same for the $L^2$-norm result and maximum norm, the spatial scales of the structures in the residential and industrial images are more similar to the target data by the $L^2$-norm.
Based on these results, we use the $L^2$-norm as the best similarity measurement method for WIK in the following discussion.

With the result from WIK, we can investigate what data were used for training the model and what data instance pairs the model captures as similar. The purpose of the XAI methods is to explain model inferences to humans, and we should consider how we receive the results. If we consider the presented pair to be reasonably similar, it shows that the model extracts features without being counterintuitive. In other words, we can verify whether abstraction by the model is reasonable and consistent with that of humans. In addition, we can verify whether the model was trained on the dataset sufficiently covering the feature space to estimate the positioning of the target. If we do not feel the pair is similar, we should suspect that the training dataset is insufficient in the feature space around the target or that the model is insufficiently trained. In general, the inference result is unreliable in such situations.

Figure~\ref{fig_result2} shows results similar to those in Fig.~\ref{fig_result1}, but for the cases in which the images do not look so similar.
As we can see, the images are not similar for other metrics as well as for the $L^2$-norm results.
These are rare cases with the model and the training dataset shown here.  If these results are dominant, we would suspect that the quality of the model is not good.
In the highway images, a junction is observed on the result in the fourth row.  The model might interpret that the highway target data shows a crossroad and it is similar to the highway junction.
The selected result for the annual crop image has a different color from the target, and we can see that colors are not important to the model for areas of this type.
For the target data showing a pasture, the inference result is incorrect and the data is classified as a river. This causes results to be selected from the incorrect class, thus leading to inconsistency between the target data and the results.
\begin{figure}[!t]
\centering
\includegraphics[width=3.1in]{result2.pdf}
\caption{
Selected images by the presented XAI method with different similarity measurements shown as in Fig.~\ref{fig_result1} for cases in which the images appear dissimilar.
A predicted label is shown for each image. The predictions of the images labelled as showing a river are incorrect; these images should have been correctly labeled as pasture.
}
\label{fig_result2}
\end{figure}

We are unaware of any effective way to evaluate example-based XAI methods because the value of those methods depends on whether the users satisfy the outputs. While one of the approaches to verify the self-consistency of the methods is shown\cite{hanawa2021}, that metric inevitably becomes 100\% with our method because it is requested in step 1 of the process. For real-world applications, it is critical that stakeholders understand the method's positioning and are persuaded to use it.

\section{Discussion}
We successfully showed that the presented XAI method WIK could provide information to verify that the AI model is appropriately trained with the data-centric approach, and this result cannot be achieved by simply using mathematical distance metrics for input data. The images in the sixth row in Fig.~\ref{fig_result1} are also images from the training dataset but are selected by the criterion that the $L^2$-norm
of the pixel values of the images themselves,
not of the features extracted by the model used in the WIK method,
with respect to the target data is the smallest.
Although it is also a similarity metric, the selected images are different from those selected by the presented XAI method shown in the fourth row. The spatial structures and contrast are weaker than the original images and the image in the second row, and the images are relatively smooth and uniform. This is caused by the manner in which the $L^2$-norm is calculated.
For spatially structured features, the $L^2$-norm
of the pixel values rather than that of the feature vectors
is high if the pattern does not match perfectly.
It is unrealistic that there is an image with a perfectly matching structure pattern, and the $L^2$-norm cannot extract abstract features of images related to spatial structures. Therefore, the $L^2$-norm measurement of the original images is unsuitable for the XAI method. By comparing the images in the second and third rows, it is clearly shown that the XAI method described in the previous section can successfully select images with similar features in the context of the task to classify what is in the image.

From the results of WIK, as shown in Fig.~\ref{fig_result1}, we can check whether each inference is suspicious during the production phase of AI model applications. When the inference results are used for making decisions, we can use the XAI results as discussion material in addition to the inference results themselves. Even when a machine learning model is used as part of a system for automation, although there might not be a natural way to include the XAI method in the system itself, the XAI result is available for performance monitoring. This perspective is closely related to the idea of MLOps for discussions on how to maintain machine-learning models. How we can detect degradations of machine learning models is under discussion, and there is a possibility that example-based XAI methods, such as those presented here, can help investigate the relationship between the distribution of recent target data and the model.
As the performance of the model degrades, instances selected by the WIK method might become less intuitively similar owing to discrepancies between the latest input data and the training dataset acquired in the past. Similarity values may also worsen with degradation.

Because the purpose of WIK demonstrated here is to obtain additional information to verify the reliability of the model inferences, this method is effective in the production and proof-of-concept phase. We may have to suspect
that related instances that should not be separated, or even exactly the same instance, are included in both of the training and validation datasets
if the pair is \textit{too} similar.
Conversely, if the pair does not look similar, we could improve the generalization performance of the model by investigating its origin. There are several possible explanations for this finding. The variation of the training dataset might be insufficient, or hyperparameters might not be set to accurate values. In addition, there is a possibility that this was caused by only a simple mistake. The proposed XAI method can help avoid such inappropriate situations.

If a large number of instances, such as the entire test dataset, need to be evaluated with WIK, the result can be summarized by performing a visual inspection of the selected instances with the highest and lowest similarity values. Identical instances can be identified in the training and test datasets owing to their complete similarity. The lowest similarity result indicates an area with few training data instances in the feature space. Investigating the highest and lowest similarity instances for each class can help understand the distribution of the training data in the feature space.

As we have shown in the introduction, the WIK method presented in this article is not necessarily effective for all purposes of XAIs. In real-world applications, selecting an appropriate method for each case is important. If the analysis focuses on which explanatory variable is important, importance analysis methods for the input data, such as LIME or SHAP, would be an appropriate choice. In such a case, WIK is ineffective. Conversely, if the reliability of the AI model is a concern, the presented method provides a unique perspective for model validity investigation using a data-centric approach.
The importance analysis methods provide a component analytical approach, whereas WIK provides reference data the model trained on.
Therefore, the importance analysis methods and WIK are complementary in understanding machine learning models from multiple perspectives.

Decision-making based on AI predictions could be a major application of the XAI method. We must be careful about risks when making decisions, and it is dangerous to use only inference results in many cases. In addition to the fact that AI inferences are not always perfect, we should be careful if validations with performance metrics do not guarantee prediction accuracy through false implicit assumptions. This could be caused by the difference in the data distributions between the training dataset and the new data to be inferred. Few XAI methods effective for this situation have been suggested, and the method presented here is a candidate to overcome this difficulty.

\section{Conclusion}
We presented an XAI method WIK to investigate the reliability of deep learning models by showing a data instance that has the most similar features in the context of the task of the model.
Our approach differs from the XAI methods with the importance analysis of explanatory variables and helps to intuitively understand how the model infers.
This method can verify the training data, an important component of machine learning, along with the models, from a data-centric perspective.
Experimental results on a remote sensing image classification task with the EuroSAT dataset, we confirmed that the WIK method performs well by selecting reasonably similar instances compared to the input data.
With this method, applications such as AI reliability monitoring, validity checks for AI developments, and more careful decision-making are expected.

\section*{Acknowledgments}
We would like to thank Prof. Katsuki Fujisawa and graduate students in his laboratory to have a discussion and give us useful comments.

\bibliographystyle{elsarticle-num}


\section{Introduction}

For a smooth, closed and uniformly convex hypersurface $\Sigma$ in $\mathbb{R}^{n+1}$, the Gauss
map $n: \Sigma\rightarrow \mathbb{S}^n$ is a diffeomorphism. Then, the Weingarten matrix
\begin{eqnarray*}
\mathcal{W}=d n
\end{eqnarray*}
is positive definite. The $k$-th fundamental symmetric function
of the principal curvatures $\kappa_1, ..., \kappa_n$ (or the eigenvalues of $\mathcal{W}$)
\begin{eqnarray*}
H_k=\sum_{i_1<i_2< \cdot\cdot\cdot<
i_k}\kappa_{i_1}\cdot\cdot\cdot\kappa_{i_k}
\end{eqnarray*}
is called
the $k$-th mean curvature of $M$.
A fundamental question in classical
differential geometry concerns how much one can recover through the inverse
Gauss map when some information is prescribed on $\mathbb{S}^n$ \cite{Yau}. The most notable
example is probably the problem of finding
a closed and uniformly convex hypersurface in $\mathbb{R}^{n+1}$ whose
$k$-th mean curvature ($1\leq k\leq n$) is prescribed as a positive function $f$
defined on $\mathbb{S}^n$, i.e.
\begin{eqnarray}\label{GG}
H_{k}(n^{-1}(x))=f(x).
\end{eqnarray}
For $k=n$, this is the famous Minkowski problem (see \cite{Sch13} for a comprehensive introduction).
For $1\leq k<n$, this problem has been considered by Guan-Guan in \cite{Guan-g}.

In this paper, we want to consider the similar problem for uniformly $h$-convex hypersurfaces in hyperbolic space.
In order to formulate this problem, at first, we will briefly describe $h$-convex geometry in hyperbolic space
followed by Section 5 in \cite{And20} and Section 2 in \cite{Li-Xu}. Their works deeply declared
interesting formal similarities between the geometry of
$h$-convex domains in hyperbolic space and that of
convex Euclidean bodies.

We shall work
in the hyperboloid model of $\mathbb{H}^{n+1}$. For that, consider the
Minkowski space $\mathbb{R}^{n+1,1}$ with canonical coordinates
$X=(X^1, . . . , X^{n+1}, X^0)$ and the Lorentzian metric
\begin{eqnarray*}
\langle X, X\rangle=\sum_{i=1}^{n+1}(X^i)^2-(X^0)^2.
\end{eqnarray*}
$\mathbb{H}^{n+1}$ is the future time-like hyperboloid in Minkowski
space $\mathbb{R}^{n+1,1}$, i.e.
\begin{eqnarray*}
\mathbb{H}^{n+1}=\Big\{X=(X^1, \cdot\cdot\cdot, X^{n+1}, X^0) \in
\mathbb{R}^{n+1,1}: \langle X, X\rangle=-1, X^0>0\Big\}.
\end{eqnarray*}
The horospheres are hypersurfaces in $\mathbb{H}^{n+1}$ whose
principal curvatures equal to $1$ everywhere. In the hyperboloid model of $\mathbb{H}^{n+1}$,
they can be parameterized by $\mathbb{S}^n \times \mathbb{R}$
\begin{eqnarray*}
H_{x}(r)=\{X \in \mathbb{H}^{n+1}: \langle X, (x, 1)\rangle=-e^r\},
\end{eqnarray*}
where $x \in \mathbb{S}^n$ and $r \in \mathbb{R}$ represents the signed geodesic distance from the ``north pole" $N =(0, 1)\in
\mathbb{H}^{n+1}$.
The interior of the horosphere is called the horo-ball and we denote by
\begin{eqnarray*}
B_{x}(r)=\{X \in \mathbb{H}^{n+1}: 0>\langle X, (x,
1)\rangle>-e^r\}.
\end{eqnarray*}
If we use the Poincar\'e ball model $\mathbb{B}^{n+1}$ of $\mathbb{H}^{n+1}$, then $B_x(r)$ corresponds to
an $(n+1)$-dimensional ball which tangents to $\partial \mathbb{B}^{n+1}$ at $x$. Furthermore, $B_x(r)$ contains the origin for $r>0$.
\begin{definition}
A compact domain $\Omega\subset \mathbb{H}^{n+1}$ (or its boundary $\partial \Omega$) is horospherically convex
(or $h$-convex for short) if every boundary point $p$ of
$\partial \Omega$ has a supporting horo-ball, i.e. a horo-ball $B$
such that $\Omega\subset \overline{B}$ and $p \in \partial B$. When $\Omega$ is smooth, it is $h$-convex if
and only if the principal curvatures of $\partial \Omega$ are greater
than or equal to $1$.

For a smooth compact domain $\Omega$, we say
$\Omega$ (or $\partial \Omega$) is uniformly $h$-convex if its principal curvatures are greater than $1$.
\end{definition}

\begin{definition}\label{HG}
Let $\Omega\subset \mathbb{H}^{n+1}$ be a $h$-convex compact domain. For each $X \in \partial \Omega$,
$\partial \Omega$ has a supporting horo-ball $B_{x}(r)$
for some $r \in \mathbb{R}$ and $x \in \mathbb{S}^n$. Then the horospherical
Gauss map $G: \partial \Omega\rightarrow S_{\infty}=(\mathbb{S}^n, g_{\infty})$ of $\Omega$ (or $\partial \Omega$)
is defined by
\begin{eqnarray*}
G(X)=x, \quad g_{\infty}(X)=e^{2r}\sigma,
\end{eqnarray*}
where $\sigma$ is the canonical $\mathbb{S}^n$ metric.
\end{definition}

Note that the canonical $\mathbb{S}^n$ metric $\sigma$ is used in order to measure geometric quantities associated to the
Euclidean Gauss map of a hypersurface in $\mathbb{R}^{n+1}$. However, Espinar-G\'alvez-Mirain
explained in detail why we use the horospherical metric $g_{\infty}$ on $\mathbb{S}^n$
for measuring geometrical quantities with respect to the horospherical Gauss map (see section 2 in \cite{Esp09}).
Let $M$ be a $h$-convex hypersurface in $\mathbb{H}^{n+1}$.
For each $X \in M$, $M$  has a supporting horo-ball $B_{x}(r)$.
Then, we have (see (5.3) in \cite{And20} or (2.2) in \cite{Li-Xu})
\begin{eqnarray*}
X-\nu=e^{-r}\langle x, 1\rangle,
\end{eqnarray*}
where $\nu$ is the unit outward vector of $M$. Differentiating the above equation gives
\begin{eqnarray}\label{dG}
\langle d G, d G\rangle_{g_{\infty}}=(\mathcal{W}-I)^2,
\end{eqnarray}
where $\mathcal{W}$ is Weingarten matrix of $M$ and $I$ is the identity matrix.

The relation \eqref{dG} declares that the matrix $\mathcal{W}-I$ plays the role in $h$-convex hyperbolic geometry
of the matrix $\mathcal{W}$ in convex Euclidean geometry. This fact
motivates Andrews-Chen-Wei \cite{And20} to define
the shifted Weingarten matrix by
$\widetilde{\mathcal{W}}:=\mathcal{W}-I$.
Clearly, a smooth hypersurface $M\subset \mathbb{H}^{n+1}$ is uniformly $h$-convex
if and only if its shifted Weingarten matrix $\widetilde{\mathcal{W}}$ is positive definite. Thus, $G$ is a diffeomorphism from $M$ to $\mathbb{S}^n$ if $M$ is uniformly $h$-convex.

Let $\kappa_1, . . ., \kappa_n$ be principal curvatures of $M$, the shifted principal curvatures is defined by \cite{And20}
\begin{eqnarray*}
(\tilde{\kappa}_1, . . . , \tilde{\kappa}_n):=(\kappa_1 -1, . . ., \kappa_n -1),
\end{eqnarray*}
which are eigenvalues of the shifted Weingarten matrix
$\widetilde{\mathcal{W}}$.  Thus, the hyperbolic curvature radii take the form (see Definition 8 in \cite{Esp09})
\begin{eqnarray*}
\mathcal{R}_i:=\frac{1}{\kappa_i-1}.
\end{eqnarray*}
Espinar-G\'alvez-Mirain showed the hyperbolic curvature radii plays the role in hyperbolic space of the Euclidean
curvature radii from several different perspectives and used it to extend the
Christoffel problem \cite{Chr65, Fi} to hyperbolic space.

Using the shifted principal curvatures, the $k$-th shifted mean curvature for $1\leq k\leq n$ can be defined by
\begin{eqnarray*}
\widetilde{H}_k=\sum_{i_1<i_2< \cdot\cdot\cdot<
i_k}\tilde{\kappa}_{i_1}\cdot\cdot\cdot\tilde{\kappa}_{i_k},
\end{eqnarray*}
which is used by Andrews-Chen-Wei \cite{And20} to define
the modified quermassintegrals and the corresponding Alexandrov-Fenchel type inequalities
have proved in \cite{And20, Hu20} by using shifted curvature flows in hyperbolic space.
Later, Wang-Wei-Zhou \cite{Wang20} studied
inverse shifted curvature flow in hyperbolic space.

In this paper, we consider prescribed shifted Weingarten curvatures problem in hyperbolic space
which is motivated by the work of Guan-Guan \cite{Guan-g} on
the similar problem \eqref{GG} in Euclidean space .

\begin{problem}\label{LX}
Let $1\leq k\leq n$ be a fixed integer. For a given smooth positive function $\tilde{f}(x)$ defined on $\mathbb{S}^n$,
does there exist a smooth, closed and uniformly $h$-convex hypersurface $M\subset \mathbb{H}^{n+1}$ satisfying
\begin{eqnarray}\label{PC}
\widetilde{H}_{k}(G^{-1}(x))=\tilde{f}(x).
\end{eqnarray}
\end{problem}

This problem is also a special case of the generalized Christoffel problem (see (5.17) in \cite{Esp09}).
For $k=n$, this problem is prescribed shifted Gauss problem which has been studied in \cite{Li-Xu, Chen23}.
In this paper, we solve Problem \ref{LX} when $\tilde{f}$ is an even function, i.e.
$\tilde{f}(x)=\tilde{f}(-x)$ for all $x \in \mathbb{S}^n$.
\begin{theorem}\label{Main}
Assume $1\leq k<n$, there exists a smooth, closed, origin-symmetric and uniformly $h$-convex hypersurface in
$\mathbb{H}^{n+1}$
satisfying the equation \eqref{PC} for any smooth positive even function $\tilde{f}$ defined on $\mathbb{S}^n$.
\end{theorem}

In Sect. 2, we will show that Problem \ref{LX} is reduced to solve a Hessian quotient
equation on $\mathbb{S}^n$.  After establishing the a priori estimates for solutions to the Hessian quotient
equation in Sect. 3, we will use the degree theory to prove Theorem \ref{Main} in Sect. 4.

\section{The Hessian quotient
equation associated to Problem \ref{LX}}

In this section, using the the horospherical support function of a
$h$-convex hypersurface in $\mathbb{H}^{n+1}$, we can reduce Problem \ref{LX} to solve a Hessian quotient
equation on $\mathbb{S}^n$.

We will continue to review $h$-convex geometry in hyperbolic space in Section 1
followed by Section 5 in \cite{And20} and Section 2 in \cite{Li-Xu}.
We also work
in the hyperboloid model of $\mathbb{H}^{n+1}$ as Section 1.
Let $\Omega$ be a $h$-convex compact domain in $\mathbb{H}^{n+1}$. Then
for each $x\in \mathbb{S}^n$ we define the horospherical support function
of $\Omega$ (or $\partial \Omega$) in direction $x$ by
\begin{eqnarray*}
u(x):=\inf\{s \in \mathbb{R}: \Omega\subset \overline{B}_{x}(s)\}.
\end{eqnarray*}
We also have the alternative characterisation
\begin{eqnarray}\label{SD}
u(x)=\sup\{\log(-\langle X, (x, 1)\rangle): X \in \Omega\}.
\end{eqnarray}
The support function completely determines a
$h$-convex compact domain $\Omega$, as an intersection of horo-balls:
\begin{eqnarray*}
\Omega=\bigcap_{x \in \mathbb{S}^n}\overline{B}_{x}(u(x)).
\end{eqnarray*}

Since $G$ is a diffeomorphism from $\partial \Omega$ to $\mathbb{S}^n$ for a compact uniformly $h$-convex domain $\Omega$. Then,
$\overline{X}=G^{-1}$ is a smooth embedding from $\mathbb{S}^n$ to $\partial \Omega$ and
$\overline{X}$ can be written in terms of the support function $u$, as follows:
\begin{eqnarray}\label{X}
\overline{X}(x)=\frac{1}{2}\varphi(-x, 1)+\frac{1}{2}
\Big(\frac{|D \varphi|^2}{\varphi}+\frac{1}{\varphi}\Big)(x, 1)-(D\varphi, 0),
\end{eqnarray}
where $\varphi=e^u$ and
$D$ is the Levi-Civita connection of the standard metric $\sigma$ of $\mathbb{S}^{n}$. Then,
after choosing normal coordinates around $x$ on $\mathbb{S}^{n+1}$, we
express the shifted Weingarten matrix in the horospherical support function
(see (1.16) in \cite{And20}, Lemma 2.2 in \cite{Li-Xu})
\begin{eqnarray*}
\widetilde{\mathcal{W}}=\Big(\varphi A[\varphi]\Big)^{-1},
\end{eqnarray*}
where
\begin{eqnarray*}
A[\varphi]=D^2\varphi-\frac{1}{2}\frac{|D\varphi|^2}{\varphi}I+\frac{1}{2}\Big(\varphi-\frac{1}{\varphi}\Big)I
\end{eqnarray*}
Thus, $\Omega \subset \mathbb{H}^{n+1}$ is  uniformly $h$-convex if and only if
the matrix $A[\varphi]$ is positive definite.

So, Problem \ref{LX} is equivalent to find a smooth positive solution $\varphi(x)$ with $A[\varphi(x)]>0$ for all $x \in \mathbb{S}^n$
to the equation
\begin{eqnarray}\label{MA}
\frac{\sigma_n(A[\varphi])}{\sigma_{n-k}(A[\varphi])}=\varphi^{-k}f(x),
\end{eqnarray}
where $f=\tilde{f}^{-1}$ and $\sigma_k(A)$ is the $k$-th elementary symmetric function of a symmetric matrix $A$. Next, we will give the definition of the elementary symmetric functions and review their basic properties which could be found in
\cite{L96}.

\begin{definition}
For any $k=1,2,\cdots,n$, we set
$$\sigma_{k}(\lambda)=\sum\limits_{1\le i_1<i_2<\cdots < i_k\le n}
\lambda_{{i}_{1}}\lambda_{{i}_{2}}\cdots\lambda_{{i}_{k}},\leqno(1.1.1)$$
for any $\lambda=(\lambda_{1},\cdots,\lambda_{n})\in\mathbb{R}^{n}$ and set $\sigma_0(\lambda)=1$.
Let $\lambda_1(A)$, ..., $\lambda_n(A)$ be the eigenvalues of $A$ and
denote by $\lambda(A)=(\lambda_1(A), ..., \lambda_n(A)$. We define by $\sigma_{k}(A)=\sigma_{k}(\lambda(A))$.
\end{definition}

We recall that the Garding's cone is defined as
$$\Gamma_{k}=\{\lambda\in \mathbb{R}^{n}:\sigma_{i}(\lambda)>0,\forall 1\le i\le k\}.$$

The following two propositions will be used later.
\begin{proposition}(Generalized Newtown-MacLaurin
inequality) For $\lambda\in\Gamma_{k}$ and $k>l\geq0$, $r>s\geq0$,
$k\geq{r}$, $l\geq{s}$, we have
\begin{eqnarray}\label{GNM}
\bigg[\frac{\sigma_{k}(\lambda)/C_{n}^{k}}{\sigma_{l}(\lambda)/C_{n}^{l}}\bigg]^{\frac{1}{k-l}}
\leq\bigg[\frac{\sigma_{r}(\lambda)/C_{n}^{r}}{\sigma_{s}(\lambda)/
C_{n}^{s}}\bigg]^{\frac{1}{r-s}}
\end{eqnarray}
and the equality holds if $\lambda_{1}=...=\lambda_{n}>0.$
\end{proposition}

\begin{proposition}\label{con}
For $\lambda\in\Gamma_{k}$, we have for $n\geq k>l\geq0$
\begin{eqnarray}\label{el}
\frac{\partial}{\partial \lambda_i}\bigg[\frac{\sigma_{k}(\lambda)}{\sigma_{l}(\lambda)}\bigg]>0
\end{eqnarray}
for $1\leq i\leq n$
and
\begin{eqnarray*}
\bigg[\frac{\sigma_{k}(\lambda)}{\sigma_{l}(\lambda)}\bigg]^{\frac{1}{k-l}}
\end{eqnarray*}
is a concave function.
Moreover, we have
\begin{eqnarray}\label{sum}
\sum_{i=1}^{n}\frac{\partial[\frac{\sigma_{k}}{\sigma_{l}}]
^{\frac{1}{k-l}}(\lambda)}{\partial\lambda_{i}}
\geq[\frac{C_{n}^{k}}{C_{n}^{l}}]^{\frac{1}{k-l}}.
\end{eqnarray}
\end{proposition}

\section{The a priori estimates}

For convenience, in the following of this paper, we always assume that $f$ is
a smooth positive, even function on $\mathbb{S}^n$ and $\varphi$ is a smooth even solution
to the equation \eqref{MA} with $A[\varphi]>0$.
Moreover, let $M$ be the smooth, closed and uniformly $h$-convex hypersurface in $\mathbb{H}^{n+1}$ with the horospherical support function $u=\log \varphi$.
Clearly, $M$ is symmetry with the origin and $\varphi(x)> 1$ for $x \in \mathbb{S}^n$.

The following easy and important equality is key for
the $C^0$ estimate.
\begin{lemma}
We have
\begin{eqnarray}\label{c0-12}
\frac{1}{2}\Big(\max_{\mathbb{S}^n}\varphi+\frac{1}{\max_{\mathbb{S}^n}\varphi}\Big)\leq \min_{\mathbb{S}^n}\varphi.
\end{eqnarray}
\end{lemma}

\begin{proof}
The inequality can be found in the proof of Lemma 7.2 in \cite{Li-Xu}. For completeness, we give a proof here.
Assume that $\varphi(x_1)=\max_{\mathbb{S}^n}\varphi$ and denote $\overline{X}(x_1)=G^{-1}(x_1)$ as before. Then, we have for any $x \in \mathbb{S}^n$ by
the definition of the horospherical support function \eqref{SD}
\begin{eqnarray*}
-\langle \overline{X}(x_1), (x, 1)\rangle\leq \varphi(x), \quad \forall x \in \mathbb{S}^n.
\end{eqnarray*}
Substituting the expression \eqref{X} for $\overline{X}$ into the above equality yields
\begin{eqnarray}\label{D0-3}
\frac{1}{2}\varphi(x_1)(1+\langle x_1, x\rangle)+\frac{1}{2}\frac{1}{\varphi(x_1)}(1-\langle x_1, x\rangle)\leq \varphi(x),
\end{eqnarray}
where we used the fact $D\varphi(x_1)=0$. Note that $\varphi(x_1)\geq 1$, we find from \eqref{D0-3}
\begin{eqnarray}\label{D0-4}
\frac{1}{2}\Big(\varphi(x_1)+\frac{1}{\varphi(x_1)}\Big)\leq \varphi(x) \quad \mbox{for} \quad \langle x, x_1\rangle\geq 0.
\end{eqnarray}
Since $\varphi$ is even, we can assume
that the minimum point $x_0$ of $\varphi(x)$ satisfies $\langle x_0, x_1\rangle\geq 0$. Thus,
the equality \eqref{c0-12} follows that from \eqref{D0-4}.
\end{proof}

Now, we use the maximum principle to get the $C^0$-estimate.
\begin{lemma}\label{C-C0}
We have
\begin{eqnarray}\label{C0}
0<\frac{1}{C}\leq u(x)\leq C, \quad \forall \ x \in \mathbb{S}^n,
\end{eqnarray}
where $C$ is a positive constant depending on $k$, $n$ and $f$.
\end{lemma}

\begin{proof}
Applying the maximum principle, we have from the equation \eqref{MA}
\begin{eqnarray*}
(\max_{\mathbb{S}^n}\varphi)^{2k}\frac{1}{2^k}\bigg[1-\frac{1}{(\max_{\mathbb{S}^n}\varphi)^2}\bigg]^k\geq C>0
\end{eqnarray*}
and
\begin{eqnarray*}
(\min_{\mathbb{S}^n}\varphi)^{2k}\frac{1}{2^k}\bigg[1-\frac{1}{(\min_{\mathbb{S}^n}\varphi)^2}\bigg]^k\leq C.
\end{eqnarray*}
Since the function
\begin{eqnarray*}
g(x)=x^{2k}\frac{1}{2^k}\bigg(1-\frac{1}{x^2}\bigg)^k
\end{eqnarray*}
is increasing in $[1, +\infty)$, $g(1)=0$ and $g(+\infty)=+\infty$, we obtain
\begin{eqnarray}\label{c0-11}
\min_{\mathbb{S}^n}\varphi\leq C, \quad \mbox{and} \quad \max_{\mathbb{S}^n}\varphi\geq C>1.
\end{eqnarray}
Combining \eqref{c0-11} and \eqref{c0-12}, we find
\begin{eqnarray*}
1<C\leq\min_{\mathbb{S}^n}\varphi\leq \max_{\mathbb{S}^n}\varphi\leq C^{\prime},
\end{eqnarray*}
which implies that
\begin{eqnarray*}
0<\frac{1}{C}\leq\min_{\mathbb{S}^n}u\leq \max_{\mathbb{S}^n}u\leq C.
\end{eqnarray*}
So, we complete the proof.
\end{proof}

As a corollary, we have the gradient estimate from Lemma 7.3 in \cite{Li-Xu}.

\begin{corollary}\label{C-C1}
We have
\begin{eqnarray}\label{C1}
|D\varphi(x)|\leq C, \quad \forall \ x \in \mathbb{S}^n,
\end{eqnarray}
where $C$ is a positive constant depending only on the constant in Lemma \ref{C-C0}.
\end{corollary}

We give some notations before considering the $C^2$ estimate.
Denote by
\begin{eqnarray*}
U_{ij}=\varphi_{ij}-\frac{1}{2}\frac{|D\varphi|^2}{\varphi}\delta_{ij}+\frac{1}{2}(\varphi-\frac{1}{\varphi})\delta_{ij}
\end{eqnarray*}
and
\begin{eqnarray*}
F(U)=\Big[\frac{\sigma_n(U)}{\sigma_{n-k}(U)}\Big]^{\frac{1}{k}}, \quad F^{ij}=\frac{\partial F}{\partial U_{ij}},
\quad F^{ij, st}=\frac{\partial^2 F}{\partial U_{ij}\partial U_{st}}.
\end{eqnarray*}

\begin{lemma}
We have for $1\leq i\leq n$
\begin{eqnarray}\label{C2}
\lambda_i(U(x))\leq C, \quad \forall \ x \in \mathbb{S}^n,
\end{eqnarray}
where $\lambda_1(U), ..., \lambda_n(U)$ are eigenvalues of the matrix $U$ and $C$ is a positive constant depending
only on the constant in Lemma \ref{C-C0} and Corollary \ref{C-C1}.
\end{lemma}

\begin{proof}
Since
\begin{eqnarray*}
\lambda_i(U)\leq \mathrm{tr} U=\Delta \varphi-\frac{n}{2\varphi}|D\varphi|^2+\frac{n}{2}(\varphi-\frac{1}{\varphi}), \quad \forall 1\leq i\leq n,
\end{eqnarray*}
it is sufficient to prove $\Delta\varphi\leq C$ in view of $C^0$ estimate \eqref{C0} and $C^1$ estimate \eqref{C1}.
Moreover, these two estimates \eqref{C0} and \eqref{C1} together with the positivity of the matrix $U$ imply $\lambda_{i}(D^2\varphi)\geq -C$ for $1\leq i\leq n$.
Thus, we find
\begin{eqnarray}\label{C2-3}
|\lambda_i(D^2\varphi)|\leq C \max \{ \max_{\mathbb{S}^n}\Delta \varphi, 1\}, \quad \forall 1\leq i\leq n.
\end{eqnarray}
We take the auxiliary function
$$W(x)=\Delta \varphi.$$
Assume $x_0$ is the maximum point of $W$. After an appropriate
choice of the normal frame at $x_0$, we further assume $U_{ij}$, hence $\varphi_{ij}$ and $F^{ij}$ is diagonal at the point $x_0$. Then,
\begin{equation}\label{W1}
W_i(x_0)=\sum_{q}\varphi_{qqi}=0,
\end{equation}
and
\begin{equation}\label{W2}
W_{ii}(x_0)=\sum_{q}\varphi_{qqii}\leq0.
\end{equation}
Using \eqref{W2} and the positivity of $F^{ij}$ given by \eqref{el}, we arrive at $x_0$ if $\Delta \varphi$ is large enough
\begin{eqnarray*}
0&\ge&\sum_iF^{ii}W_{ii}=\sum_iF^{ii}\sum_{q}\varphi_{qqii}\geq\sum_iF^{ii}\sum_{q}\big(\varphi_{iiqq}-C\Delta \varphi\big),
\end{eqnarray*}
where we use Ricci identity and the equality \eqref{C2-3} to get the last inequality. Thus it
follows from the definition of $U$, \eqref{C0}, \eqref{C1} and \eqref{W1},
\begin{eqnarray}\label{C2-1}
0&\ge&\sum_iF^{ii}\sum_{q}\bigg[U_{iiqq}+\Big(\frac{1}{2\varphi}|D\varphi|^2\Big)_{qq}-
\frac{1}{2}\Big(\varphi-\frac{1}{\varphi}\Big)_{qq}\bigg]
-C\sum_{i}F^{ii}\Delta \varphi\nonumber\\
&\ge&\sum_iF^{ii}\sum_{q}\bigg[U_{iiqq}+\frac{1}{\varphi}(\varphi_{qq})^2-C\Delta\varphi-C\bigg]
-C\sum_iF^{ii}\Delta \varphi.
\end{eqnarray}
Differentiating the equation \eqref{MA} twice gives
\begin{eqnarray*}
F^{ii}U_{iiqq}+F^{ij, st}U_{ijq}U_{stq}=(\varphi^kf)_{qq}.
\end{eqnarray*}
Since $F$ is concave (see Proposition \ref{con}), it yields
\begin{eqnarray}\label{C2-2}
F^{ii}U_{iiqq}\geq -C\Delta \varphi-C,
\end{eqnarray}
where we used \eqref{C0} and \eqref{C1}.
Substituting \eqref{C2-2} into \eqref{C2-1} and using $(\Delta \varphi)^2 \leq n \sum_{q}(\varphi_{qq})^2$, we have
\begin{eqnarray}\label{Q}
0\ge\bigg(C(\Delta \varphi)^2-C\Delta \varphi-C\bigg)\sum_{i}F^{ii}-C\Delta \varphi-C.
\end{eqnarray}
Applying the inequality \eqref{sum}, we find
\begin{eqnarray}\label{C2-8}
\sum_{i}F^{ii}\geq C.
\end{eqnarray}
Then we conclude at $x_0$ by combining the inequalities \eqref{C2-8} and \eqref{Q}
\begin{eqnarray*}
C\geq |\Delta \varphi|^2
\end{eqnarray*}
if $\Delta \varphi$ is chosen large enough.
So, we complete the proof.
\end{proof}

\section{The proof of the main theorem}

In this section, we use the degree theory for nonlinear elliptic
equations developed in \cite{Li89} to prove Theorem \ref{Main}. Such approach was also used
in prescribed curvature problem for star-shaped hypersurfaces \cite{An, Li02, Jin, Li-Sh}, prescribed curvature problem for
convex hypersurfaces
\cite{Guan02, Guan21} and the Gaussian
Minkowski type problem \cite{Huang-Xi21, Liu22, Feng1, Feng2}.

For the use of the degree theory, the uniqueness of constant solutions to the equation
\eqref{MA} is important for us.

\begin{lemma}\label{U-C}
The $h$-convex solutions to the equation
\begin{eqnarray}\label{MA-c}
\varphi^k\frac{\sigma_n}{\sigma_{n-k}}(A[\varphi(x)])=\gamma
\end{eqnarray}
with $\varphi>1$ are given by
\begin{eqnarray*}
\varphi(x)=\Big(1+2\gamma^{\frac{1}{k}}\Big)^{\frac{1}{2}}\Big(\sqrt{|x_0|^2+1}-\langle x_0, x\rangle\Big),
\end{eqnarray*}
where $x_0 \in \mathbb{R}^{n+1}$. In particular, $\varphi(x)=\Big(1+2\gamma^{\frac{1}{k}}\Big)^{\frac{1}{2}}$ is the unique even solution.
\end{lemma}

\begin{proof}
This lemma is a corollary of Proposition 8.1 in \cite{Li-Xu}, its proof is similar to that of Theorem 8.1 (7).
\end{proof}

Now, we begin to use the degree theory to prove Theorem \ref{Main}.
After establishing the  a priori estimates \eqref{C0}, \eqref{C1} and \eqref{C2} and noting
that $\sigma_n(U)\geq C>0$ which is given by using generalized Newtown-MacLaurin
inequality \eqref{GNM} and the equation \eqref{MA}, we know that the
equation \eqref{MA} is uniformly elliptic, i.e.
\begin{eqnarray}\label{C2+++}
\lambda_i(U(x))\geq C>0, \quad \forall \ x \in \mathbb{S}^n,
\end{eqnarray}
where $\lambda_1(U), ..., \lambda_n(U)$ are eigenvalues of the matrix $U$.
From Evans-Krylov estimates \cite{Eva82, Kry83} and Schauder estimates \cite{GT}, we have
\begin{eqnarray}\label{C2+}
|\varphi|_{C^{4,\alpha}(\mathbb{S}^n)}\leq C
\end{eqnarray}
for any smooth, even and uniformly $h$-convex solution $\varphi$ to the equation \eqref{MA}.
We define
\begin{eqnarray*}
\mathcal{B}^{2,\alpha}(\mathbb{S}^n)=\{\varphi \in
C^{2,\alpha}(\mathbb{S}^n): \varphi \ \mbox{is even}\}
\end{eqnarray*}
and
\begin{eqnarray*}
\mathcal{B}_{0}^{4,\alpha}(\mathbb{S}^n)=\{\varphi \in
C^{4,\alpha}(\mathbb{S}^n): A[\varphi]>0 \ \mbox{and} \ \varphi \ \mbox{is even}\}.
\end{eqnarray*}
Let us consider $$\mathcal{L}(\cdot, t): \mathcal{B}_{0}^{4,\alpha}(\mathbb{S}^n)\rightarrow
\mathcal{B}^{2,\alpha}(\mathbb{S}^n),$$ which is defined by
\begin{eqnarray*}
\mathcal{L}(\varphi, t)=\frac{\sigma_n}{\sigma_{n-k}}(U)-\varphi^{-k} [(1-t)\gamma+t f],
\end{eqnarray*}
where the constant $\gamma$ will be chosen later and $U$ is denoted as before
\begin{eqnarray*}
U=D^2\varphi-\frac{1}{2}\frac{|D\varphi|^2}{\varphi}I+\frac{1}{2}\Big(\varphi-\frac{1}{\varphi}\Big)I.
\end{eqnarray*}

Let $$\mathcal{O}_R=\{\varphi \in \mathcal{B}_{0}^{4,\alpha}(\mathbb{S}^n):
1+\frac{1}{R}< \varphi, \ \frac{1}{R} I< U, \ |\varphi|_{C^{4,\alpha}(\mathbb{S}^n)}<R\},$$ which clearly is an open
set of $\mathcal{B}_{0}^{4,\alpha}(\mathbb{S}^n)$. Moreover, if $R$ is
sufficiently large, $\mathcal{L}(\varphi, t)=0$ has no solution on $\partial
\mathcal{O}_R$ by the a priori estimates established in \eqref{C0}, \eqref{C2+++} and \eqref{C2+}.
Therefore the degree $\deg(\mathcal{L}(\cdot, t), \mathcal{O}_R, 0)$ is
well-defined for $0\leq t\leq 1$. Using the homotopic invariance of
the degree (Proposition 2.2 in \cite{Li89}), we have
\begin{eqnarray}\label{hot}
\deg(\mathcal{L}(\cdot, 1), \mathcal{O}_R, 0)=\deg(\mathcal{L}(\cdot, 0), \mathcal{O}_R, 0).
\end{eqnarray}

Lemma \ref{U-C} tells us that
$\varphi=c$ is the unique even solution for $\mathcal{L}(\varphi, 0)=0$ in $\mathcal{O}_R$.
Direct calculation show that the linearized operator of $\mathcal{L}$ at
$\varphi=c$ is
\begin{eqnarray*}
L_{c}(\psi)=a(n, k, c)(\Delta_{\mathbb{S}^n}+n)\psi,
\end{eqnarray*}
where
\begin{eqnarray*}
a(n, k, c)=\frac{k}{n}\frac{2}{c-c^{-1}}\frac{\sigma_n(\frac{1}{2}(c-c^{-1})I)}{\sigma_{n-k}(\frac{1}{2}(c-c^{-1})I)}.
\end{eqnarray*}
Since
$\Delta_{\mathbb{S}^n}\psi+n\psi=0$
has the unique even solution $\psi=0$,
$L_{c}$ is an invertible operator. So, we have by Proposition
2.3 in \cite{Li89}
\begin{eqnarray*}
\deg(\mathcal{L}(\cdot, 0), \mathcal{O}_R, 0)=\deg(L_{c_0}, \mathcal{O}_R, 0).
\end{eqnarray*}
Because
the eigenvalues of the Beltrami-Laplace operator $\Delta$ on $\mathbb{S}^n$ are strictly less than
$-n$ except for the first two eigenvalues $0$ and $-n$,
there is only one positive eigenvalue $na(n, k, c)$ of $L_{c}$
with multiplicity $1$.
Then we have by Proposition
2.4 in \cite{Li89}
\begin{eqnarray*}
\deg(\mathcal{L}(\cdot, 0), \mathcal{O}_R, 0)=\deg(L_{c_0}, \mathcal{O}_R, 0)=-
1.
\end{eqnarray*}
Therefore, it follows from \eqref{hot}
\begin{eqnarray*}
\deg(\mathcal{L}(\cdot, 1), \mathcal{O}_R; 0)=\deg(\mathcal{L}(\cdot, 0), \mathcal{O}_R, 0)=-
1.
\end{eqnarray*}
So, we obtain a solution at $t=1$. This completes the proof of
Theorem \ref{Main}.

\bigskip

\bigskip